RETROSPECTIVE VARYING COEFFICIENT ASSOCIATION ANALYSIS OF LONGITUDINAL BINARY TRAITS: APPLICATION TO THE IDENTIFICATION OF GENETIC LOCI ASSOCIATED WITH HYPERTENSION.

Many genetic studies contain rich information on longitudinal phenotypes that require powerful analytical tools for optimal analysis. Genetic analysis of longitudinal data that incorporates temporal variation is important for understanding the genetic architecture and biological variation of complex diseases. Most of the existing methods assume that the contribution of genetic variants is constant over time and fail to capture the dynamic pattern of disease progression. However, the relative influence of genetic variants on complex traits fluctuates over time. In this study, we propose a retrospective varying coefficient mixed model association test, RVMMAT, to detect time-varying genetic effect on longitudinal binary traits. We model dynamic genetic effect using smoothing splines, estimate model parameters by maximizing a double penalized quasi-likelihood function, design a joint test using a Cauchy combination method, and evaluate statistical significance via a retrospective approach to achieve robustness to model misspecification. Through simulations we illustrated that the retrospective varying-coefficient test was robust to model misspecification under different ascertainment schemes and gained power over the association methods assuming constant genetic effect. We applied RVMMAT to a genome-wide association analysis of longitudinal measure of hypertension in the Multi-Ethnic Study of Atherosclerosis. Pathway analysis identified two important pathways related to G-protein signaling and DNA damage. Our results demonstrated that RVMMAT could detect biologically relevant loci and pathways in a genome scan and provided insight into the genetic architecture of hypertension.

Computationally inferred cell-type specific epigenome-wide DNA methylation analysis unveils distinct methylation patterns among immune cells for HIV infection in three cohorts.

BACKGROUND: Epigenome-wide association studies (EWAS) have identified CpG sites associated with HIV infection in blood cells in bulk, which offer limited knowledge of cell-type specific methylation patterns associated with HIV infection. In this study, we aim to identify differentially methylated CpG sites for HIV infection in immune cell types: CD4+ T-cells, CD8+ T-cells, B cells, Natural Killer (NK) cells, and monocytes. METHODS: Applying a computational deconvolution method, we performed a cell-type based EWAS for HIV infection in three independent cohorts (Ntotal = 1,382). DNA methylation in blood or in peripheral blood mononuclear cells (PBMCs) was profiled by an array-based method and then deconvoluted by Tensor Composition Analysis (TCA). The TCA-computed CpG methylation in each cell type was first benchmarked by bisulfite DNA methylation capture sequencing in a subset of the samples. Cell-type EWAS of HIV infection was performed in each cohort separately and a meta-EWAS was conducted followed by gene set enrichment analysis. RESULTS: The meta-analysis unveiled a total of 2,021 cell-type unique significant CpG sites for five inferred cell types. Among these inferred cell-type unique CpG sites, the concordance rate in the three cohorts ranged from 96% to 100% in each cell type. Cell-type level meta-EWAS unveiled distinct patterns of HIV-associated differential CpG methylation, where 74% of CpG sites were unique to individual cell types (false discovery rate, FDR <0.05). CD4+ T-cells had the largest number of unique HIV-associated CpG sites (N = 1,624) compared to any other cell type. Genes harboring significant CpG sites are involved in immunity and HIV pathogenesis (e.g. CD4+ T-cells: NLRC5, CX3CR1, B cells: IFI44L, NK cells: IL12R, monocytes: IRF7), and in oncogenesis (e.g. CD4+ T-cells: BCL family, PRDM16, monocytes: PRDM16, PDCD1LG2). HIV-associated CpG sites were enriched among genes involved in HIV pathogenesis and oncogenesis that were enriched among interferon-α and -γ, TNF-α, inflammatory response, and apoptotic pathways. CONCLUSION: Our findings uncovered computationally inferred cell-type specific modifications in the host epigenome for people with HIV that contribute to the growing body of evidence regarding HIV pathogenesis.

tRFtarget 2.0: expanding the targetome landscape of transfer RNA-derived fragments.

tRFtarget 1.0 (http://trftarget.net/) is a platform consolidating both computationally predicted and experimentally validated binding sites between transfer RNA-derived fragments (tRFs) and target genes (or transcripts) across multiple organisms. Here, we introduce a newly released version of tRFtarget 2.0, in which we integrated 6 additional tRF sources, resulting in a comprehensive collection of 2614 high-quality tRF sequences spanning across 9 species, including 1944 Homo sapiens tRFs and one newly incorporated species Rattus norvegicus. We also expanded target genes by including ribosomal RNAs, long non-coding RNAs, and coding genes >50 kb in length. The predicted binding sites have surged up to approximately 6 billion, a 20.5-fold increase than that in tRFtarget 1.0. The manually curated publications relevant to tRF targets have increased to 400 and the gene-level experimental evidence has risen to 232. tRFtarget 2.0 introduces several new features, including a web-based tool that identifies potential binding sites of tRFs in user's own datasets, integration of standardized tRF IDs, and inclusion of external links to contents within the database. Additionally, we enhanced website framework and user interface. With these improvements, tRFtarget 2.0 is more user-friendly, providing researchers a streamlined and comprehensive platform to accelerate their research progress.

Transformer with convolution and graph-node co-embedding: An accurate and interpretable vision backbone for predicting gene expressions from local histopathological image.

Inferring gene expressions from histopathological images has long been a fascinating yet challenging task, primarily due to the substantial disparities between the two modality. Existing strategies using local or global features of histological images are suffering model complexity, GPU consumption, low interpretability, insufficient encoding of local features, and over-smooth prediction of gene expressions among neighboring sites. In this paper, we develop TCGN (Transformer with Convolution and Graph-Node co-embedding method) for gene expression estimation from H&E-stained pathological slide images. TCGN comprises a combination of convolutional layers, transformer encoders, and graph neural networks, and is the first to integrate these blocks in a general and interpretable computer vision backbone. Notably, TCGN uniquely operates with just a single spot image as input for histopathological image analysis, simplifying the process while maintaining interpretability. We validate TCGN on three publicly available spatial transcriptomic datasets. TCGN consistently exhibited the best performance (with median PCC 0.232). TCGN offers superior accuracy while keeping parameters to a minimum (just 86.241 million), and it consumes minimal memory, allowing it to run smoothly even on personal computers. Moreover, TCGN can be extended to handle bulk RNA-seq data while providing the interpretability. Enhancing the accuracy of omics information prediction from pathological images not only establishes a connection between genotype and phenotype, enabling the prediction of costly-to-measure biomarkers from affordable histopathological images, but also lays the groundwork for future multi-modal data modeling. Our results confirm that TCGN is a powerful tool for inferring gene expressions from histopathological images in precision health applications.

Brain Registration and Evaluation for Zebrafish (BREEZE)-mapping: A pipeline for whole-brain structural and activity analyses.

Here, we present Brain Registration and Evaluation for Zebrafish (BREEZE)-mapping, a user-friendly pipeline for the registration and analysis of whole-brain images in larval zebrafish. We describe steps for pre-processing, registration, quantification, and visualization of whole-brain phenotypes in zebrafish mutants of genes associated with neurodevelopmental and neuropsychiatric disorders. By utilizing BioImage Suite Web, an open-source software package originally developed for processing human brain imaging data, we provide a highly accessible whole-brain mapping protocol developed for users with general computational proficiency. For complete details on the use and execution of this protocol, please refer to Weinschutz Mendes et al. (2023).1.

Detecting time-varying genetic effects in Alzheimer's disease using a longitudinal GWAS model.

Background: The development and progression of Alzheimer's disease (AD) is a complex process that can change over time, during which genetic influences on phenotypes may also fluctuate. Incorporating longitudinal phenotypes in genome wide association studies (GWAS) could help unmask genetic loci with time-varying effects. In this study, we incorporated a varying coefficient test in a longitudinal GWAS model to identify single nucleotide polymorphisms (SNPs) that may have time- or age-dependent effects in AD. Methods: Genotype data from 1,877 participants in the Alzheimer's Neuroimaging Data Initiative (ADNI) were imputed using the Haplotype Reference Consortium (HRC) panel, resulting in 9,573,130 SNPs. Subjects' longitudinal impairment status at each visit was considered as a binary and clinical phenotype. Participants' composite standardized uptake value ratio (SUVR) derived from each longitudinal amyloid PET scan was considered as a continuous and biological phenotype. The retrospective varying coefficient mixed model association test (RVMMAT) was used in longitudinal GWAS to detect time-varying genetic effects on the impairment status and SUVR measures. Post-hoc analyses were performed on genome-wide significant SNPs, including 1) pathway analyses; 2) age-stratified genotypic comparisons and regression analyses; and 3) replication analyses using data from the National Alzheimer's Coordinating Center (NACC). Results: Our model identified 244 genome-wide significant SNPs that revealed time-varying genetic effects on the clinical impairment status in AD; among which, 12 SNPs on chromosome 19 were successfully replicated using data from NACC. Post-hoc age-stratified analyses indicated that for most of these 244 SNPs, the maximum genotypic effect on impairment status occurred between 70 to 80 years old, and then declined with age. Our model further identified 73 genome-wide significant SNPs associated with the temporal variation of amyloid accumulation. For these SNPs, an increasing genotypic effect on PET-SUVR was observed as participants' age increased. Functional pathway analyses on significant SNPs for both phenotypes highlighted the involvement and disruption of immune responses- and neuroinflammation-related pathways in AD. Conclusion: We demonstrate that longitudinal GWAS models with time-varying coefficients can boost the statistical power in AD-GWAS. In addition, our analyses uncovered potential time-varying genetic variants on repeated measurements of clinical and biological phenotypes in AD.

A hybrid machine learning and regression method for cell type deconvolution of spatial barcoding-based transcriptomic data.

Spatial barcoding-based transcriptomic (ST) data require cell type deconvolution for cellular-level downstream analysis. Here we present SDePER, a hybrid machine learning and regression method, to deconvolve ST data using reference single-cell RNA sequencing (scRNA-seq) data. SDePER uses a machine learning approach to remove the systematic difference between ST and scRNA-seq data (platform effects) explicitly and efficiently to ensure the linear relationship between ST data and cell type-specific expression profile. It also considers sparsity of cell types per capture spot and across-spots spatial correlation in cell type compositions. Based on the estimated cell type proportions, SDePER imputes cell type compositions and gene expression at unmeasured locations in a tissue map with enhanced resolution. Applications to coarse-grained simulated data and four real datasets showed that SDePER achieved more accurate and robust results than existing methods, suggesting the importance of considering platform effects, sparsity and spatial correlation in cell type deconvolution.

Cis-meQTL for cocaine use-associated DNA methylation in an HIV-positive cohort show pleiotropic effects on multiple traits.

BACKGROUND: Cocaine use (CU) is associated with psychiatric and medical diseases. Little is known about the mechanisms of CU-related comorbidities. Findings from preclinical and clinical studies have suggested that CU is associated with aberrant DNA methylation (DNAm) that may be influenced by genetic variants [i.e., methylation quantitative trait loci (meQTLs)]. In this study, we mapped cis-meQTLs for CU-associated DNAm sites (CpGs) in an HIV-positive cohort (Ntotal = 811) and extended the meQTLs to multiple traits. RESULTS: We conducted cis-meQTL analysis for 224 candidate CpGs selected for their association with CU in blood. We identified 7,101 significant meQTLs [false discovery rate (FDR) < 0.05], which mostly mapped to genes involved in immunological functions and were enriched in immune pathways. We followed up the meQTLs using phenome-wide association study and trait enrichment analyses, which revealed 9 significant traits. We tested for causal effects of CU on these 9 traits using Mendelian Randomization and found evidence that CU plays a causal role in increasing hypertension (p-value = 2.35E-08) and decreasing heel bone mineral density (p-value = 1.92E-19). CONCLUSIONS: These findings suggest that genetic variants for CU-associated DNAm have pleiotropic effects on other relevant traits and provide new insights into the causal relationships between cocaine use and these complex traits.

iDESC: identifying differential expression in single-cell RNA sequencing data with multiple subjects.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) technology has enabled assessment of transcriptome-wide changes at single-cell resolution. Due to the heterogeneity in environmental exposure and genetic background across subjects, subject effect contributes to the major source of variation in scRNA-seq data with multiple subjects, which severely confounds cell type specific differential expression (DE) analysis. Moreover, dropout events are prevalent in scRNA-seq data, leading to excessive number of zeroes in the data, which further aggravates the challenge in DE analysis. RESULTS: We developed iDESC to detect cell type specific DE genes between two groups of subjects in scRNA-seq data. iDESC uses a zero-inflated negative binomial mixed model to consider both subject effect and dropouts. The prevalence of dropout events (dropout rate) was demonstrated to be dependent on gene expression level, which is modeled by pooling information across genes. Subject effect is modeled as a random effect in the log-mean of the negative binomial component. We evaluated and compared the performance of iDESC with eleven existing DE analysis methods. Using simulated data, we demonstrated that iDESC had well-controlled type I error and higher power compared to the existing methods. Applications of those methods with well-controlled type I error to three real scRNA-seq datasets from the same tissue and disease showed that the results of iDESC achieved the best consistency between datasets and the best disease relevance. CONCLUSIONS: iDESC was able to achieve more accurate and robust DE analysis results by separating subject effect from disease effect with consideration of dropouts to identify DE genes, suggesting the importance of considering subject effect and dropouts in the DE analysis of scRNA-seq data with multiple subjects.

Detection of interactions between genetic marker sets and environment in a genome-wide study of hypertension.

As human complex diseases are influenced by the interplay of genes and environment, detecting gene-environment interactions (G×E) can shed light on biological mechanisms of diseases and play an important role in disease risk prediction. Development of powerful quantitative tools to incorporate G×E in complex diseases has potential to facilitate the accurate curation and analysis of large genetic epidemiological studies. However, most of existing methods that interrogate G×E focus on the interaction effects of an environmental factor and genetic variants, exclusively for common or rare variants. In this study, we proposed two tests, MAGEIT_RAN and MAGEIT_FIX, to detect interaction effects of an environmental factor and a set of genetic markers containing both rare and common variants, based on the MinQue for Summary statistics. The genetic main effects in MAGEIT_RAN and MAGEIT_FIX are modeled as random or fixed, respectively. Through simulation studies, we illustrated that both tests had type I error under control and MAGEIT_RAN was overall the most powerful test. We applied MAGEIT to a genome-wide analysis of gene-alcohol interactions on hypertension in the Multi-Ethnic Study of Atherosclerosis. We detected two genes, CCNDBP1 and EPB42, that interact with alcohol usage to influence blood pressure. Pathway analysis identified sixteen significant pathways related to signal transduction and development that were associated with hypertension, and several of them were reported to have an interactive effect with alcohol intake. Our results demonstrated that MAGEIT can detect biologically relevant genes that interact with environmental factors to influence complex traits.

Differences in Mortality Among Patients With Asthma and COPD Hospitalized With COVID-19.

BACKGROUND: It remains unclear whether patients with asthma and/or chronic obstructive pulmonary disease (COPD) are at increased risk for severe coronavirus disease 2019 (COVID-19). OBJECTIVE: Compare in-hospital COVID-19 outcomes among patients with asthma, COPD, and no airway disease. METHODS: A retrospective cohort study was conducted on 8,395 patients admitted with COVID-19 between March 2020 and April 2021. Airway disease diagnoses were defined using International Classification of Diseases, 10th Revision codes. Mortality and sequential organ failure assessment (SOFA) scores were compared among groups. Logistic regression analysis was used to identify and adjust for confounding clinical features associated with mortality. RESULTS: The median SOFA score in patients without airway disease was 0.32 and mortality was 11%. In comparison, asthma patients had lower SOFA scores (median 0.15; P < .01) and decreased mortality, even after adjusting for age, diabetes, and other confounders (odds ratio 0.65; P = .01). Patients with COPD had higher SOFA scores (median 0.86; P < .01) and increased adjusted odds of mortality (odds ratio 1.40; P < .01). Blood eosinophil count of 200 cells/µL or greater, a marker of type 2 inflammation, was associated with lower mortality across all groups. Importantly, patients with asthma showed improved outcomes even after adjusting for eosinophilia, indicating that noneosinophilic asthma was associated with protection as well. CONCLUSIONS: COVID-19 severity was increased in patients with COPD and decreased in those with asthma, eosinophilia, and noneosinophilic asthma, independent of clinical confounders. These findings suggest that COVID-19 severity may be influenced by intrinsic immunological factors in patients with airway diseases, such as type 2 inflammation.

Systematically identifying genetic signatures including novel SNP-clusters, nonsense variants, frame-shift INDELs, and long STR expansions that potentially link to unknown phenotypes existing in dog breeds.

BACKGROUND: In light of previous studies that profiled breed-specific traits or used genome-wide association studies to refine loci associated with characteristic morphological features in dogs, the field has gained tremendous genetic insights for known dog traits observed among breeds. Here we aim to address the question from a reserve perspective: whether there are breed-specific genotypes that may underlie currently unknown phenotypes. This study provides a complete set of breed-specific genetic signatures (BSGS). Several novel BSGS with significant protein-altering effects were highlighted and validated. RESULTS: Using the next generation whole-genome sequencing technology coupled with unsupervised machine learning for pattern recognitions, we constructed and analyzed a high-resolution sequence map for 76 breeds of 412 dogs. Genomic structures including novel single nucleotide polymorphisms (SNPs), SNP clusters, insertions, deletions (INDELs) and short tandem repeats (STRs) were uncovered mutually exclusively among breeds. We also partially validated some novel nonsense variants by Sanger sequencing with additional dogs. Four novel nonsense BSGS were found in the Bernese Mountain Dog, Samoyed, Bull Terrier, and Basset Hound, respectively. Four INDELs resulting in either frame-shift or codon disruptions were found in the Norwich Terrier, Airedale Terrier, Chow Chow and Bernese Mountain Dog, respectively. A total of 15 genomic regions containing three types of BSGS (SNP-clusters, INDELs and STRs) were identified in the Akita, Alaskan Malamute, Chow Chow, Field Spaniel, Keeshond, Shetland Sheepdog and Sussex Spaniel, in which Keeshond and Sussex Spaniel each carried one amino-acid changing BSGS in such regions. CONCLUSION: Given the strong relationship between human and dog breed-specific traits, this study might be of considerable interest to researchers and all. Novel genetic signatures that can differentiate dog breeds were uncovered. Several functional genetic signatures might indicate potentially breed-specific unknown phenotypic traits or disease predispositions. These results open the door for further investigations. Importantly, the computational tools we developed can be applied to any dog breeds as well as other species. This study will stimulate new thinking, as the results of breed-specific genetic signatures may offer an overarching relevance of the animal models to human health and disease.

High-throughput functional analysis of autism genes in zebrafish identifies convergence in dopaminergic and neuroimmune pathways.

Advancing from gene discovery in autism spectrum disorders (ASDs) to the identification of biologically relevant mechanisms remains a central challenge. Here, we perform parallel in vivo functional analysis of 10 ASD genes at the behavioral, structural, and circuit levels in zebrafish mutants, revealing both unique and overlapping effects of gene loss of function. Whole-brain mapping identifies the forebrain and cerebellum as the most significant contributors to brain size differences, while regions involved in sensory-motor control, particularly dopaminergic regions, are associated with altered baseline brain activity. Finally, we show a global increase in microglia resulting from ASD gene loss of function in select mutants, implicating neuroimmune dysfunction as a key pathway relevant to ASD biology.

Deep learning-based morphological feature analysis and the prognostic association study in colon adenocarcinoma histopathological images.

Colorectal cancer (CRC) is now the third most common malignancy to cause mortality worldwide, and its prognosis is of great importance. Recent CRC prognostic prediction studies mainly focused on biomarkers, radiometric images, and end-to-end deep learning methods, while only a few works paid attention to exploring the relationship between the quantitative morphological features of patients' tissue slides and their prognosis. However, existing few works in this area suffered from the drawback of choosing the cells randomly from the whole slides, which contain the non-tumor region that lakes information about prognosis. In addition, the existing works, which tried to demonstrate their biological interpretability using patients' transcriptome data, failed to show the biological meaning closely related to cancer. In this study, we proposed and evaluated a prognostic model using morphological features of cells in the tumor region. The features were first extracted by the software CellProfiler from the tumor region selected by Eff-Unet deep learning model. Features from different regions were then averaged for each patient as their representative, and the Lasso-Cox model was used to select the prognosis-related features. The prognostic prediction model was at last constructed using the selected prognosis-related features and was evaluated through KM estimate and cross-validation. In terms of biological meaning, Gene Ontology (GO) enrichment analysis of the expressed genes that correlated with the prognostically significant features was performed to show the biological interpretability of our model.With the help of tumor segmentation, our model achieved better statistical significance and better biological interpretability compared to the results without tumor segmentation. Statistically, the Kaplan Meier (KM) estimate of our model showed that the model using features in the tumor region has a higher C-index, a lower p-value, and a better performance on cross-validation than the model without tumor segmentation. In addition, revealing the pathway of the immune escape and the spread of the tumor, the model with tumor segmentation demonstrated a biological meaning much more related to cancer immunobiology than the model without tumor segmentation. Our prognostic prediction model using quantitive morphological features from tumor regions was almost as good as the TNM tumor staging system as they had a close C-index, and our model can be combined with the TNM tumor stage system to make a better prognostic prediction. And to the best of our knowledge, the biological mechanisms in our study were the most relevant to the immune mechanism of cancer compared to the previous studies.

Single-cell multi-omics reveals dyssynchrony of the innate and adaptive immune system in progressive COVID-19.

Dysregulated immune responses against the SARS-CoV-2 virus are instrumental in severe COVID-19. However, the immune signatures associated with immunopathology are poorly understood. Here we use multi-omics single-cell analysis to probe the dynamic immune responses in hospitalized patients with stable or progressive course of COVID-19, explore V(D)J repertoires, and assess the cellular effects of tocilizumab. Coordinated profiling of gene expression and cell lineage protein markers shows that S100Ahi/HLA-DRlo classical monocytes and activated LAG-3hi T cells are hallmarks of progressive disease and highlights the abnormal MHC-II/LAG-3 interaction on myeloid and T cells, respectively. We also find skewed T cell receptor repertories in expanded effector CD8+ clones, unmutated IGHG+ B cell clones, and mutated B cell clones with stable somatic hypermutation frequency over time. In conclusion, our in-depth immune profiling reveals dyssynchrony of the innate and adaptive immune interaction in progressive COVID-19.

Testing gene-environment interactions in the presence of confounders and mismeasured environmental exposures.

Interest in investigating gene-environment (GxE) interactions has rapidly increased over the last decade. Although GxE interactions have been extremely investigated in large studies, few such effects have been identified and replicated, highlighting the need to develop statistical GxE tests with greater statistical power. The reverse test has been proposed for testing the interaction effect between continuous exposure and genetic variants in relation to a binary disease outcome, which leverages the idea of linear discriminant analysis, significantly increasing statistical power comparing to the standard logistic regression approach. However, this reverse approach did not take into consideration adjustment for confounders. Since GxE interaction studies are inherently nonexperimental, adjusting for potential confounding effects is critical for valid evaluation of GxE interactions. In this study, we extend the reverse test to allow for confounders. The proposed reverse test also allows for exposure measurement errors as typically occurs. Extensive simulation experiments demonstrated that the proposed method not only provides greater statistical power under most simulation scenarios but also provides substantive computational efficiency, which achieves a computation time that is more than sevenfold less than that of the standard logistic regression test. In an illustrative example, we applied the proposed approach to the Veterans Aging Cohort Study (VACS) to search for genetic susceptibility loci modifying the smoking-HIV status association.

Impact of Neoantigen Expression and T-Cell Activation on Breast Cancer Survival.

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.

A novel transcriptional risk score for risk prediction of complex human diseases.

Recently polygenetic risk score (PRS) has been successfully used in the risk prediction of complex human diseases. Many studies incorporated internal information, such as effect size distribution, or external information, such as linkage disequilibrium, functional annotation, and pleiotropy among multiple diseases, to optimize the performance of PRS. To leverage on multiomics datasets, we developed a novel flexible transcriptional risk score (TRS), in which messenger RNA expression levels were imputed and weighted for risk prediction. In simulation studies, we demonstrated that single-tissue TRS has greater prediction power than LDpred, especially when there is a large effect of gene expression on the phenotype. Multitissue TRS improves prediction accuracy when there are multiple tissues with independent contributions to disease risk. We applied our method to complex traits, including Crohn's disease, type 2 diabetes, and so on. The single-tissue TRS method outperformed LDpred and AnnoPred across the tested traits. The performance of multitissue TRS is trait-dependent. Moreover, our method can easily incorporate information from epigenomic and proteomic data upon the availability of reference datasets.

G2S3: A gene graph-based imputation method for single-cell RNA sequencing data.

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.