Calcium transients on the ER surface trigger liquid-liquid phase separation of FIP200 to specify autophagosome initiation sites.

The mechanism that initiates autophagosome formation on the ER in multicellular organisms is elusive. Here, we showed that autophagy stimuli trigger Ca2+ transients on the outer surface of the ER membrane, whose amplitude, frequency, and duration are controlled by the metazoan-specific ER transmembrane autophagy protein EPG-4/EI24. Persistent Ca2+ transients/oscillations on the cytosolic ER surface in EI24-depleted cells cause accumulation of FIP200 autophagosome initiation complexes on the ER. This defect is suppressed by attenuating ER Ca2+ transients. Multi-modal SIM analysis revealed that Ca2+ transients on the ER trigger the formation of dynamic and fusion-prone liquid-like FIP200 puncta. Starvation-induced Ca2+ transients on lysosomes also induce FIP200 puncta that further move to the ER. Multiple FIP200 puncta on the ER, whose association depends on the ER proteins VAPA/B and ATL2/3, assemble into autophagosome formation sites. Thus, Ca2+ transients are crucial for triggering phase separation of FIP200 to specify autophagosome initiation sites in metazoans.

Neoantigen-driven B cell and CD4 T follicular helper cell collaboration promotes anti-tumor CD8 T cell responses.

CD4 T follicular helper (TFH) cells support B cells, which are critical for germinal center (GC) formation, but the importance of TFH-B cell interactions in cancer is unclear. We found enrichment of TFH cell transcriptional signature correlates with GC B cell signature and with prolonged survival in individuals with lung adenocarcinoma (LUAD). We further developed a murine LUAD model in which tumor cells express B cell- and T cell-recognized neoantigens. Interactions between tumor-specific TFH and GC B cells, as well as interleukin (IL)-21 primarily produced by TFH cells, are necessary for tumor control and effector CD8 T cell function. Development of TFH cells requires B cells and B cell-recognized neoantigens. Thus, tumor neoantigens can regulate the fate of tumor-specific CD4 T cells by facilitating their interactions with tumor-specific B cells, which in turn promote anti-tumor immunity by enhancing CD8 T cell effector functions.

Pan-cancer ion transport signature reveals functional regulators of glioblastoma aggression.

Ion channels, transporters, and other ion-flux controlling proteins, collectively comprising the "ion permeome", are common drug targets, however, their roles in cancer remain understudied. Our integrative pan-cancer transcriptome analysis shows that genes encoding the ion permeome are significantly more often highly expressed in specific subsets of cancer samples, compared to pan-transcriptome expectations. To enable target selection, we identified 410 survival-associated IP genes in 33 cancer types using a machine-learning approach. Notably, GJB2 and SCN9A show prominent expression in neoplastic cells and are associated with poor prognosis in glioblastoma, the most common and aggressive brain cancer. GJB2 or SCN9A knockdown in patient-derived glioblastoma cells induces transcriptome-wide changes involving neuron projection and proliferation pathways, impairs cell viability and tumor sphere formation in vitro, perturbs tunneling nanotube dynamics, and extends the survival of glioblastoma-bearing mice. Thus, aberrant activation of genes encoding ion transport proteins appears as a pan-cancer feature defining tumor heterogeneity, which can be exploited for mechanistic insights and therapy development.

Integration of expression QTLs with fine mapping via SuSiE.

Genome-wide association studies (GWASs) have achieved remarkable success in associating thousands of genetic variants with complex traits. However, the presence of linkage disequilibrium (LD) makes it challenging to identify the causal variants. To address this critical gap from association to causation, many fine-mapping methods have been proposed to assign well-calibrated probabilities of causality to candidate variants, taking into account the underlying LD pattern. In this manuscript, we introduce a statistical framework that incorporates expression quantitative trait locus (eQTL) information to fine-mapping, built on the sum of single-effects (SuSiE) regression model. Our new method, SuSiE2, connects two SuSiE models, one for eQTL analysis and one for genetic fine-mapping. This is achieved by first computing the posterior inclusion probabilities (PIPs) from an eQTL-based SuSiE model with the expression level of the candidate gene as the phenotype. These calculated PIPs are then utilized as prior inclusion probabilities for risk variants in another SuSiE model for the trait of interest. By prioritizing functional variants within the candidate region using eQTL information, SuSiE2 improves SuSiE by increasing the detection rate of causal SNPs and reducing the average size of credible sets. We compared the performance of SuSiE2 with other multi-trait fine-mapping methods with respect to power, coverage, and precision through simulations and applications to the GWAS results of Alzheimer's disease (AD) and body mass index (BMI). Our results demonstrate the better performance of SuSiE2, both when the in-sample linkage disequilibrium (LD) matrix and an external reference panel is used in inference.

Decoding transcriptomic signatures of cysteine string protein alpha-mediated synapse maintenance.

Synapse maintenance is essential for generating functional circuitry, and decrement in this process is a hallmark of neurodegenerative disease. Yet, little is known about synapse maintenance in vivo. Cysteine string protein α (CSPα), encoded by the Dnajc5 gene, is a synaptic vesicle chaperone that is necessary for synapse maintenance and linked to neurodegeneration. To investigate the transcriptional changes associated with synapse maintenance, we performed single-nucleus transcriptomics on the cortex of young CSPα knockout (KO) mice and littermate controls. Through differential expression and gene ontology analysis, we observed that both neurons and glial cells exhibit unique signatures in the CSPα KO brain. Significantly, all neuronal classes in CSPα KO brains show strong signatures of repression in synaptic pathways, while up-regulating autophagy-related genes. Through visualization of synapses and autophagosomes by electron microscopy, we confirmed these alterations especially in inhibitory synapses. Glial responses varied by cell type, with microglia exhibiting activation. By imputing cell-cell interactions, we found that neuron-glia interactions were specifically increased in CSPα KO mice. This was mediated by synaptogenic adhesion molecules, with the classical Neurexin1-Neuroligin 1 pair being the most prominent, suggesting that communication of glial cells with neurons is strengthened in CSPα KO mice to preserve synapse maintenance. Together, this study provides a rich dataset of transcriptional changes in the CSPα KO cortex and reveals insights into synapse maintenance and neurodegeneration.

SDPRX: A statistical method for cross-population prediction of complex traits.

Polygenic risk score (PRS) has demonstrated its great utility in biomedical research through identifying high-risk individuals for different diseases from their genotypes. However, the broader application of PRS to the general population is hindered by the limited transferability of PRS developed in Europeans to non-European populations. To improve PRS prediction accuracy in non-European populations, we develop a statistical method called SDPRX that can effectively integrate genome wide association study summary statistics from different populations. SDPRX automatically adjusts for linkage disequilibrium differences between populations and characterizes the joint distribution of the effect sizes of a variant in two populations to be both null, population specific, or shared with correlation. Through simulations and applications to real traits, we show that SDPRX improves the prediction performance over existing methods in non-European populations.

mbDecoda: a debiased approach to compositional data analysis for microbiome surveys.

Potentially pathogenic or probiotic microbes can be identified by comparing their abundance levels between healthy and diseased populations, or more broadly, by linking microbiome composition with clinical phenotypes or environmental factors. However, in microbiome studies, feature tables provide relative rather than absolute abundance of each feature in each sample, as the microbial loads of the samples and the ratios of sequencing depth to microbial load are both unknown and subject to considerable variation. Moreover, microbiome abundance data are count-valued, often over-dispersed and contain a substantial proportion of zeros. To carry out differential abundance analysis while addressing these challenges, we introduce mbDecoda, a model-based approach for debiased analysis of sparse compositions of microbiomes. mbDecoda employs a zero-inflated negative binomial model, linking mean abundance to the variable of interest through a log link function, and it accommodates the adjustment for confounding factors. To efficiently obtain maximum likelihood estimates of model parameters, an Expectation Maximization algorithm is developed. A minimum coverage interval approach is then proposed to rectify compositional bias, enabling accurate and reliable absolute abundance analysis. Through extensive simulation studies and analysis of real-world microbiome datasets, we demonstrate that mbDecoda compares favorably with state-of-the-art methods in terms of effectiveness, robustness and reproducibility.

LDER-GE estimates phenotypic variance component of gene-environment interactions in human complex traits accurately with GE interaction summary statistics and full LD information.

Gene-environment (GE) interactions are essential in understanding human complex traits. Identifying these interactions is necessary for deciphering the biological basis of such traits. In this study, we review state-of-art methods for estimating the proportion of phenotypic variance explained by genome-wide GE interactions and introduce a novel statistical method Linkage-Disequilibrium Eigenvalue Regression for Gene-Environment interactions (LDER-GE). LDER-GE improves the accuracy of estimating the phenotypic variance component explained by genome-wide GE interactions using large-scale biobank association summary statistics. LDER-GE leverages the complete Linkage Disequilibrium (LD) matrix, as opposed to only the diagonal squared LD matrix utilized by LDSC (Linkage Disequilibrium Score)-based methods. Our extensive simulation studies demonstrate that LDER-GE performs better than LDSC-based approaches by enhancing statistical efficiency by ~23%. This improvement is equivalent to a sample size increase of around 51%. Additionally, LDER-GE effectively controls type-I error rate and produces unbiased results. We conducted an analysis using UK Biobank data, comprising 307 259 unrelated European-Ancestry subjects and 966 766 variants, across 217 environmental covariate-phenotype (E-Y) pairs. LDER-GE identified 34 significant E-Y pairs while LDSC-based method only identified 23 significant E-Y pairs with 22 overlapped with LDER-GE. Furthermore, we employed LDER-GE to estimate the aggregated variance component attributed to multiple GE interactions, leading to an increase in the explained phenotypic variance with GE interactions compared to considering main genetic effects only. Our results suggest the importance of impacts of GE interactions on human complex traits.

Eltrombopag plus diacerein vs. eltrombopag in patients with ITP: a multicenter, randomized, open-label, phase 2 trial.

This study aimed to compare the efficacy and safety of eltrombopag plus diacerein vs. eltrombopag alone in patients with primary immune thrombocytopenia (ITP) who were previously unresponsive to 14 days of eltrombopag treatment at the full dose. Recruited patients were randomly assigned 1:1 to receive either eltrombopag plus diacerein (n=50) or eltrombopag monotherapy (n=52). Overall response rate, defined as a platelet count at or above 30×109/L, at least doubling of the baseline platelet count, and no bleeding, was reached in 44% of patients in the eltrombopag plus diacerein group compared with 13% in the eltrombopag group at day 15 (P = .0009), and reached in 42% of patients in the combination group compared with 12% in the monotherapy group at day 28 (P = .0006). The addition of diacerein to eltrombopag also led to a longer duration of response (P = .0004). The two most common treatment-emergent adverse events were respiratory infection and gastrointestinal reactions in the combination group, and fatigue and respiratory infection in the eltrombopag group. In conclusion, eltrombopag plus diacerein was well tolerated, and induced higher overall response rates and longer duration of response than eltrombopag alone, offering a rejuvenating salvage therapy for ITP patients unresponsive to 14 days of full dosage eltrombopag. Our work has the potential to enhance the care of patients treated with thrombopoietin receptor agonists, reducing the need for rapid transitions to less-preferable therapies. This study is registered at ClinicalTrials.gov as NCT04917679.

A statistical framework to identify cell types whose genetically regulated proportions are associated with complex diseases.

Finding disease-relevant tissues and cell types can facilitate the identification and investigation of functional genes and variants. In particular, cell type proportions can serve as potential disease predictive biomarkers. In this manuscript, we introduce a novel statistical framework, cell-type Wide Association Study (cWAS), that integrates genetic data with transcriptomics data to identify cell types whose genetically regulated proportions (GRPs) are disease/trait-associated. On simulated and real GWAS data, cWAS showed good statistical power with newly identified significant GRP associations in disease-associated tissues. More specifically, GRPs of endothelial and myofibroblasts in lung tissue were associated with Idiopathic Pulmonary Fibrosis and Chronic Obstructive Pulmonary Disease, respectively. For breast cancer, the GRP of blood CD8+ T cells was negatively associated with breast cancer (BC) risk as well as survival. Overall, cWAS is a powerful tool to reveal cell types associated with complex diseases mediated by GRPs.

Leveraging LD eigenvalue regression to improve the estimation of SNP heritability and confounding inflation.

Heritability is a fundamental concept in genetic studies, measuring the genetic contribution to complex traits and bringing insights about disease mechanisms. The advance of high-throughput technologies has provided many resources for heritability estimation. Linkage disequilibrium (LD) score regression (LDSC) estimates both heritability and confounding biases, such as cryptic relatedness and population stratification, among single-nucleotide polymorphisms (SNPs) by using only summary statistics released from genome-wide association studies. However, only partial information in the LD matrix is utilized in LDSC, leading to loss in precision. In this study, we propose LD eigenvalue regression (LDER), an extension of LDSC, by making full use of the LD information. Compared to state-of-the-art heritability estimating methods, LDER provides more accurate estimates of SNP heritability and better distinguishes the inflation caused by polygenicity and confounding effects. We demonstrate the advantages of LDER both theoretically and with extensive simulations. We applied LDER to 814 complex traits from UK Biobank, and LDER identified 363 significantly heritable phenotypes, among which 97 were not identified by LDSC.

Benchmarking of local genetic correlation estimation methods using summary statistics from genome-wide association studies.

Local genetic correlation evaluates the correlation of additive genetic effects between different traits across the same genetic variants at a genomic locus. It has been proven informative for understanding the genetic similarities of complex traits beyond that captured by global genetic correlation calculated across the whole genome. Several summary-statistics-based approaches have been developed for estimating local genetic correlation, including $\rho$-hess, SUPERGNOVA and LAVA. However, there has not been a comprehensive evaluation of these methods to offer practical guidelines on the choices of these methods. In this study, we conduct benchmark comparisons of the performance of these three methods through extensive simulation and real data analyses. We focus on two technical difficulties in estimating local genetic correlation: sample overlaps across traits and local linkage disequilibrium (LD) estimates when only the external reference panels are available. Our simulations suggest the likelihood of incorrectly identifying correlated regions and local correlation estimation accuracy are highly dependent on the estimation of the local LD matrix. These observations are corroborated by real data analyses of 31 complex traits. Overall, our findings illuminate the distinct results yielded by different methods applied in post-genome-wide association studies (post-GWAS) local correlation studies. We underscore the sensitivity of local genetic correlation estimates and inferences to the precision of local LD estimation. These observations accentuate the vital need for ongoing refinement in methodologies.

HBV-infected hepatocellular carcinoma can be robustly classified into three clinically relevant subgroups by a novel analytical protocol.

Liver cancer is the third leading cause of cancer-related death worldwide, and hepatocellular carcinoma (HCC) accounts for a relatively large proportion of all primary liver malignancies. Among the several known risk factors, hepatitis B virus (HBV) infection is one of the important causes of HCC. In this study, we demonstrated that the HBV-infected HCC patients could be robustly classified into three clinically relevant subgroups, i.e. Cluster1, Cluster2 and Cluster3, based on consistent differentially expressed mRNAs and proteins, which showed better generalization. The proposed three subgroups showed different molecular characteristics, immune microenvironment and prognostic survival characteristics. The Cluster1 subgroup had near-normal levels of metabolism-related proteins, low proliferation activity and good immune infiltration, which were associated with its good liver function, smaller tumor size, good prognosis, low alpha-fetoprotein (AFP) levels and lower clinical stage. In contrast, the Cluster3 subgroup had the lowest levels of metabolism-related proteins, which corresponded with its severe liver dysfunction. Also, high proliferation activity and poor immune microenvironment in Cluster3 subgroup were associated with its poor prognosis, larger tumor size, high AFP levels, high incidence of tumor thrombus and higher clinical stage. The characteristics of the Cluster2 subgroup were between the Cluster1 and Cluster3 groups. In addition, MCM2-7, RFC2-5, MSH2, MSH6, SMC2, SMC4, NCPAG and TOP2A proteins were significantly upregulated in the Cluster3 subgroup. Meanwhile, abnormally high phosphorylation levels of these proteins were associated with high levels of DNA repair, telomere maintenance and proliferative features. Therefore, these proteins could be identified as potential diagnostic and prognostic markers. In general, our research provided a novel analytical protocol and insights for the robust classification, treatment and prevention of HBV-infected HCC.

A novel Bayesian framework for harmonizing information across tissues and studies to increase cell type deconvolution accuracy.

Computational cell type deconvolution on bulk transcriptomics data can reveal cell type proportion heterogeneity across samples. One critical factor for accurate deconvolution is the reference signature matrix for different cell types. Compared with inferring reference signature matrices from cell lines, rapidly accumulating single-cell RNA-sequencing (scRNA-seq) data provide a richer and less biased resource. However, deriving cell type signature from scRNA-seq data is challenging due to high biological and technical noises. In this article, we introduce a novel Bayesian framework, tranSig, to improve signature matrix inference from scRNA-seq by leveraging shared cell type-specific expression patterns across different tissues and studies. Our simulations show that tranSig is robust to the number of signature genes and tissues specified in the model. Applications of tranSig to bulk RNA sequencing data from peripheral blood, bronchoalveolar lavage and aorta demonstrate its accuracy and power to characterize biological heterogeneity across groups. In summary, tranSig offers an accurate and robust approach to defining gene expression signatures of different cell types, facilitating improved in silico cell type deconvolutions.

Statistical methods for assessing the effects of de novo variants on birth defects.

With the development of next-generation sequencing technology, de novo variants (DNVs) with deleterious effects can be identified and investigated for their effects on birth defects such as congenital heart disease (CHD). However, statistical power is still limited for such studies because of the small sample size due to the high cost of recruiting and sequencing samples and the low occurrence of DNVs. DNV analysis is further complicated by genetic heterogeneity across diseased individuals. Therefore, it is critical to jointly analyze DNVs with other types of genomic/biological information to improve statistical power to identify genes associated with birth defects. In this review, we discuss the general workflow, recent developments in statistical methods, and future directions for DNV analysis.

DNA methylation profiles of cancer-related fatigue associated with markers of inflammation and immunometabolism.

Cancer patients are commonly affected by fatigue. Herein, we sought to examine epigenetic modifications (i.e., DNA methylation) related to fatigue in peripheral blood among patients during and after treatment for head and neck cancer (HNC). Further, we determined whether these modifications were associated with gene expression and inflammatory protein markers, which we have previously linked to fatigue in HNC. This prospective, longitudinal study enrolled eligible patients with data collected at pre-radiotherapy, end of radiotherapy, and six months and one-year post-radiotherapy. Fatigue data were reported by patients using the Multidimensional Fatigue Inventory (MFI)-20. DNA methylation (Illumina MethylationEPIC) and gene expression (Applied Biosystems Clariom S) arrays and assays for seven inflammatory markers (R&D Systems multiplex) were performed. Mixed models and enrichment analyses were applied to establish the associations. A total of 386 methylation loci were associated with fatigue among 145 patients (False Discovery Rate [FDR] < 0.05). Enrichment analyses showed the involvement of genes related to immune and inflammatory responses, insulin and lipid metabolism, neuropsychological disorders, and tumors. We further identified 16 methylation-gene expression pairs (FDR < 0.05), which were linked to immune and inflammatory responses and lipid metabolism. Ninety-one percent (351) of the 386 methylation loci were also significantly associated with inflammatory markers (e.g., interleukin 6, c-reactive protein; FDR < 0.05), which further mediated the association between methylation and fatigue (FDR < 0.05). These data suggest that epigenetic modifications associated with inflammation and immunometabolism, in conjunction with relevant gene expression and protein markers, are potential targets for treating fatigue in HNC patients. The findings also merit future prospective studies in other cancer populations as well as interventional investigations.

Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder.

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed health problems, prior pharmacological treatments, and polygenic scores (PGS) has potential to inform risk stratification. We examined self-reported SB and ideation using the Columbia Suicide Severity Rating Scale (C-SSRS) among 3,942 SCZ and 5,414 BPI patients receiving care within the Veterans Health Administration (VHA). These cross-sectional data were integrated with electronic health records (EHRs), and compared across lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. PGS were constructed using available genomic data for related traits. Genome-wide association studies were performed to identify and prioritize specific loci. Only 20% of the veterans who reported SB had a corroborating ICD-9/10 EHR code. Among those without prior SB, more than 20% reported new-onset SB at follow-up. SB were associated with a range of additional clinical diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking initiation, suicide attempt, and major depressive disorder were associated with SB. The GWAS for SB yielded no significant loci. Among individuals with a diagnosed mental illness, self-reported SB were strongly associated with clinical variables across several EHR domains. Analyses point to sequelae of substance-related and psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in health records, underscoring the value of regular screening with direct, in-person assessments, especially among high-risk individuals.

The ER-Localized Transmembrane Protein EPG-3/VMP1 Regulates SERCA Activity to Control ER-Isolation Membrane Contacts for Autophagosome Formation.

During autophagosome formation in mammalian cells, isolation membranes (IMs; autophagosome precursors) dynamically contact the ER. Here, we demonstrated that the ER-localized metazoan-specific autophagy protein EPG-3/VMP1 controls ER-IM contacts. Loss of VMP1 causes stable association of IMs with the ER, thus blocking autophagosome formation. Interaction of WIPI2 with the ULK1/FIP200 complex and PI(3)P contributes to the formation of ER-IM contacts, and these interactions are enhanced by VMP1 depletion. VMP1 controls contact formation by promoting SERCA (sarco[endo]plasmic reticulum calcium ATPase) activity. VMP1 interacts with SERCA and prevents formation of the SERCA/PLN/SLN inhibitory complex. VMP1 also modulates ER contacts with lipid droplets, mitochondria, and endosomes. These ER contacts are greatly elevated by the SERCA inhibitor thapsigargin. Calmodulin acts as a sensor/effector to modulate the ER contacts mediated by VMP1/SERCA. Our study provides mechanistic insights into the establishment and disassociation of ER-IM contacts and reveals that VMP1 modulates SERCA activity to control ER contacts.

Prioritized candidate causal haplotype blocks in plant genome-wide association studies.

Genome wide association studies (GWAS) can play an essential role in understanding genetic basis of complex traits in plants and animals. Conventional SNP-based linear mixed models (LMM) that marginally test single nucleotide polymorphisms (SNPs) have successfully identified many loci with major and minor effects in many GWAS. In plant, the relatively small population size in GWAS and the high genetic diversity found in many plant species can impede mapping efforts on complex traits. Here we present a novel haplotype-based trait fine-mapping framework, HapFM, to supplement current GWAS methods. HapFM uses genotype data to partition the genome into haplotype blocks, identifies haplotype clusters within each block, and then performs genome-wide haplotype fine-mapping to prioritize the candidate causal haplotype blocks of trait. We benchmarked HapFM, GEMMA, BSLMM, GMMAT, and BLINK in both simulated and real plant GWAS datasets. HapFM consistently resulted in higher mapping power than the other GWAS methods in high polygenicity simulation setting. Moreover, it resulted in smaller mapping intervals, especially in regions of high LD, achieved by prioritizing small candidate causal blocks in the larger haplotype blocks. In the Arabidopsis flowering time (FT10) datasets, HapFM identified four novel loci compared to GEMMA's results, and the average mapping interval of HapFM was 9.6 times smaller than that of GEMMA. In conclusion, HapFM is tailored for plant GWAS to result in high mapping power on complex traits and improved on mapping resolution to facilitate crop improvement.

Network assisted analysis of de novo variants using protein-protein interaction information identified 46 candidate genes for congenital heart disease.

De novo variants (DNVs) with deleterious effects have proved informative in identifying risk genes for early-onset diseases such as congenital heart disease (CHD). A number of statistical methods have been proposed for family-based studies or case/control studies to identify risk genes by screening genes with more DNVs than expected by chance in Whole Exome Sequencing (WES) studies. However, the statistical power is still limited for cohorts with thousands of subjects. Under the hypothesis that connected genes in protein-protein interaction (PPI) networks are more likely to share similar disease association status, we developed a Markov Random Field model that can leverage information from publicly available PPI databases to increase power in identifying risk genes. We identified 46 candidate genes with at least 1 DNV in the CHD study cohort, including 18 known human CHD genes and 35 highly expressed genes in mouse developing heart. Our results may shed new insight on the shared protein functionality among risk genes for CHD.