Causality Investigation between Gut Microbiome and Sleep-Related Traits: A Bidirectional Two-Sample Mendelian Randomization Study.

Recent research has highlighted associations between sleep and microbial taxa and pathways. However, the causal effect of these associations remains unknown. To investigate this, we performed a bidirectional two-sample Mendelian randomization (MR) analysis using summary statistics of genome-wide association studies (GWAS) from 412 gut microbiome traits (N = 7738) and GWAS studies from seven sleep-associated traits (N = 345,552 to 386,577). We employed multiple MR methods to assess causality, with Inverse Variance Weighted (IVW) as the primary method, alongside a Bonferroni correction ((p < 2.4 × 10-4) to determine significant causal associations. We further applied Cochran's Q statistical analysis, MR-Egger intercept, and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) for heterogeneity and pleiotropy assessment. IVW estimates revealed 79 potential causal effects of microbial taxa and pathways on sleep-related traits and 45 inverse causal relationships, with over half related to pathways, emphasizing their significance. The results revealed two significant causal associations: genetically determined relative abundance of pentose phosphate decreased sleep duration (p = 9.00 × 10-5), and genetically determined increase in fatty acid level increased the ease of getting up in the morning (p = 8.06 × 10-5). Sensitivity analyses, including heterogeneity and pleiotropy tests, as well as a leave-one-out analysis of single nucleotide polymorphisms, confirmed the robustness of these relationships. This study explores the potential causal relationships between sleep and microbial taxa and pathways, offering novel insights into their complex interplay.

Human lineage mutations regulate RNA-protein binding of conserved genes NTRK2 and ITPR1 involved in human evolution.

Background: The role of human lineage mutations (HLMs) in human evolution through post-transcriptional modification is unclear. Aims: To investigate the contribution of HLMs to human evolution through post-transcriptional modification. Methods: We applied a deep learning model Seqweaver to predict how HLMs impact RNA-binding protein affinity. Results: We found that only 0.27% of HLMs had significant impacts on RNA-binding proteins at the threshold of the top 1% of human common variations. These HLMs enriched in a set of conserved genes highly expressed in adult excitatory neurons and prenatal Purkinje neurons, and were involved in synapse organisation and the GTPase pathway. These genes also carried excess damaging coding mutations that caused neurodevelopmental disorders, ataxia and schizophrenia. Among these genes, NTRK2 and ITPR1 had the most aggregated evidence of functional importance, suggesting their essential roles in cognition and bipedalism. Conclusions: Our findings suggest that a small subset of human-specific mutations have contributed to human speciation through impacts on post-transcriptional modification of critical brain-related genes.

Haplotype function score improves biological interpretation and cross-ancestry polygenic prediction of human complex traits.

We propose a new framework for human genetic association studies: at each locus, a deep learning model (in this study, Sei) is used to calculate the functional genomic activity score for two haplotypes per individual. This score, defined as the Haplotype Function Score (HFS), replaces the original genotype in association studies. Applying the HFS framework to 14 complex traits in the UK Biobank, we identified 3619 independent HFS-trait associations with a significance of p < 5 × 10-8. Fine-mapping revealed 2699 causal associations, corresponding to a median increase of 63 causal findings per trait compared with single-nucleotide polymorphism (SNP)-based analysis. HFS-based enrichment analysis uncovered 727 pathway-trait associations and 153 tissue-trait associations with strong biological interpretability, including 'circadian pathway-chronotype' and 'arachidonic acid-intelligence'. Lastly, we applied least absolute shrinkage and selection operator (LASSO) regression to integrate HFS prediction score with SNP-based polygenic risk scores, which showed an improvement of 16.1-39.8% in cross-ancestry polygenic prediction. We concluded that HFS is a promising strategy for understanding the genetic basis of human complex traits.

Scattered throughout the human genome are variations in the genetic code that make individuals more or less likely to develop certain traits. To identify these variants, scientists carry out Genome-wide association studies (GWAS) which compare the DNA variants of large groups of people with and without the trait of interest. This method has been able to find the underlying genes for many human diseases, but it has limitations. For instance, some variations are linked together due to where they are positioned within DNA, which can result in GWAS falsely reporting associations between genetic variants and traits. This phenomenon, known as linkage equilibrium, can be avoided by analyzing functional genomics which looks at the multiple ways a gene's activity can be influenced by a variation. For instance, how the gene is copied and decoded in to proteins and RNA molecules, and the rate at which these products are generated. Researchers can now use an artificial intelligence technique called deep learning to generate functional genomic data from a particular DNA sequence. Here, Song et al. used one of these deep learning models to calculate the functional genomics of haplotypes, groups of genetic variants inherited from one parent. The approach was applied to DNA samples from over 350 thousand individuals included in the UK BioBank. An activity score, defined as the haplotype function score (or HFS for short), was calculated for at least two haplotypes per individual, and then compared to various complex traits like height or bone density. Song et al. found that the HFS framework was better at finding links between genes and specific traits than existing methods. It also provided more information on the biology that may be underpinning these outcomes. Although more work is needed to reduce the computer processing times required to calculate the HFS, Song et al. believe that their new method has the potential to improve the way researchers identify links between genes and human traits.

Rhesus monkeys exhibiting spontaneous ritualistic behaviors resembling obsessive-compulsive disorder.

Obsessive-compulsive disorder (OCD) is a chronic and debilitating psychiatric disorder that affects ∼2%-3% of the population globally. Studying spontaneous OCD-like behaviors in non-human primates may improve our understanding of the disorder. In large rhesus monkey colonies, we found 10 monkeys spontaneously exhibiting persistent sequential motor behaviors (SMBs) in individual-specific sequences that were repetitive, time-consuming and stable over prolonged periods. Genetic analysis revealed severely damaging mutations in genes associated with OCD risk in humans. Brain imaging showed that monkeys with SMBs had larger gray matter (GM) volumes in the left caudate nucleus and lower fractional anisotropy of the corpus callosum. The GM volume of the left caudate nucleus correlated positively with the daily duration of SMBs. Notably, exposure to a stressor (human presence) significantly increased SMBs. In addition, fluoxetine, a serotonergic medication commonly used for OCD, decreased SMBs in these monkeys. These findings provide a novel foundation for developing better understanding and treatment of OCD.

Emden: A novel method integrating graph and transformer representations for predicting the effect of mutations on clinical drug response.

Precision medicine based on personalized genomics provides promising strategies to enhance the efficacy of molecular-targeted therapies. However, the clinical effectiveness of drugs has been severely limited due to genetic variations that lead to drug resistance. Predicting the impact of missense mutations on clinical drug response is an essential way to reduce the cost of clinical trials and understand genetic diseases. Here, we present Emden, a novel method integrating graph and transformer representations that predicts the effect of missense mutations on drug response through binary classification with interpretability. Emden utilized protein sequences-based features and drug structures as inputs for rapid prediction, employing competitive representation learning and demonstrating strong generalization capabilities and robustness. Our study showed promising potential for clinical drug guidance and deep insight into computer-assisted precision medicine. Emden is freely available as a web server at https://www.psymukb.net/Emden.

Inferring the Effects of Protein Variants on Protein-Protein Interactions with Interpretable Transformer Representations.

Identifying pathogenetic variants and inferring their impact on protein-protein interactions sheds light on their functional consequences on diseases. Limited by the availability of experimental data on the consequences of protein interaction, most existing methods focus on building models to predict changes in protein binding affinity. Here, we introduced MIPPI, an end-to-end, interpretable transformer-based deep learning model that learns features directly from sequences by leveraging the interaction data from IMEx. MIPPI was specifically trained to determine the types of variant impact (increasing, decreasing, disrupting, and no effect) on protein-protein interactions. We demonstrate the accuracy of MIPPI and provide interpretation through the analysis of learned attention weights, which exhibit correlations with the amino acids interacting with the variant. Moreover, we showed the practicality of MIPPI in prioritizing de novo mutations associated with complex neurodevelopmental disorders and the potential to determine the pathogenic and driving mutations. Finally, we experimentally validated the functional impact of several variants identified in patients with such disorders. Overall, MIPPI emerges as a versatile, robust, and interpretable model, capable of effectively predicting mutation impacts on protein-protein interactions and facilitating the discovery of clinically actionable variants.

A Comprehensive Study of De Novo Mutations on the Protein-Protein Interaction Interfaces Provides New Insights into Developmental Delay.

Mutations, especially those at the protein-protein interaction (PPI) interface, have been associated with various diseases. Meanwhile, though de novo mutations (DNMs) have been proven important in neuropsychiatric disorders, such as developmental delay (DD), the relationship between PPI interface DNMs and DD has not been well studied. Here we curated developmental delay DNM datasets from the PsyMuKB database and showed that DD patients showed a higher rate and deleteriousness in DNM missense on the PPI interface than sibling control. Next, we identified 302 DD-related PsychiPPIs, defined as PPIs harboring a statistically significant number of DNM missenses at their interface, and 42 DD candidate genes from PsychiPPI. We observed that PsychiPPIs preferentially affected the human protein interactome network hub proteins. When analyzing DD candidate genes using gene ontology and gene spatio-expression, we found that PsychiPPI genes carrying PPI interface mutations, such as FGFR3 and ALOX5, were enriched in development-related pathways and the development of the neocortex, and cerebellar cortex, suggesting their potential involvement in the etiology of DD. Our results demonstrated that DD patients carried an excess burden of PPI-truncating DNM, which could be used to efficiently search for disease-related genes and mutations in large-scale sequencing studies. In conclusion, our comprehensive study indicated the significant role of PPI interface DNMs in developmental delay pathogenicity.

An Integrative Analysis of Identified Schizophrenia-Associated Brain Cell Types and Gene Expression Changes.

Schizophrenia (SCZ) is a severe mental disorder that may result in hallucinations, delusions, and extremely disordered thinking. How each cell type in the brain contributes to SCZ occurrence is still unclear. Here, we leveraged the human dorsolateral prefrontal cortex bulk RNA-seq data, then used the RNA-seq deconvolution algorithm CIBERSORTx to generate SCZ brain single-cell RNA-seq data for a comprehensive analysis to understand SCZ-associated brain cell types and gene expression changes. Firstly, we observed that the proportions of brain cell types in SCZ differed from normal samples. Among these cell types, astrocyte, pericyte, and PAX6 cells were found to have a higher proportion in SCZ patients (astrocyte: SCZ = 0.163, control = 0.145, P.adj = 4.9 × 10-4, effect size = 0.478; pericyte: SCZ = 0.057, control = 0.066, P.adj = 1.1 × 10-4, effect size = 0.519; PAX6: SCZ = 0.014, control = 0.011, P.adj = 0.014, effect size = 0.377), while the L5/6_IT_CAR3 cells and LAMP5 cells are the exact opposite (L5/6_IT_Car3: SCZ = 0.102, control = 0.108, P.adj = 0.016, effect size = 0.369; LAMP5: SCZ = 0.057, control = 0.066, P.adj = 2.2 × 10-6, effect size = 0.617). Next, we investigated gene expression in cell types and functional pathways in SCZ. We observed chemical synaptic transmission dysregulation in two types of GABAergic neurons (PVALB and LAMP5), and immune reaction involvement in GABAergic neurons (SST) and non-neuronal cell types (endothelial and oligodendrocyte). Furthermore, we observed that some differential expression genes from bulk RNA-seq displayed cell-type-specific abnormalities in the expression of molecules in SCZ. Finally, the cell types with the SCZ-related transcriptomic changes could be considered to belong to the same module since we observed two major similar coordinated transcriptomic changes across these cell types. Together, our results offer novel insights into cellular heterogeneity and the molecular mechanisms underlying SCZ.

Evaluation of the Effectiveness of Derived Features of AlphaFold2 on Single-Sequence Protein Binding Site Prediction.

Though AlphaFold2 has attained considerably high precision on protein structure prediction, it is reported that directly inputting coordinates into deep learning networks cannot achieve desirable results on downstream tasks. Thus, how to process and encode the predicted results into effective forms that deep learning models can understand to improve the performance of downstream tasks is worth exploring. In this study, we tested the effects of five processing strategies of coordinates on two single-sequence protein binding site prediction tasks. These five strategies are spatial filtering, the singular value decomposition of a distance map, calculating the secondary structure feature, and the relative accessible surface area feature of proteins. The computational experiment results showed that all strategies were suitable and effective methods to encode structural information for deep learning models. In addition, by performing a case study of a mutated protein, we showed that the spatial filtering strategy could introduce structural changes into HHblits profiles and deep learning networks when protein mutation happens. In sum, this work provides new insight into the downstream tasks of protein-molecule interaction prediction, such as predicting the binding residues of proteins and estimating the effects of mutations.

Mendelian Randomization and GWAS Meta Analysis Revealed the Risk-Increasing Effect of Schizophrenia on Cancers.

The causal relationship between cancer and Schizophrenia (SCZ) remains controversial. Some researchers have found that SCZ is a cancer-preventive factor in cohort studies or meta-analyses, whereas others have found the opposite. To understand more about this issue, we used two-sample Mendelian randomization (2SMR) on available GWAS summary results to evaluate potential genetic connections between SCZ and 13 cancers. We discovered that the genetic susceptibility to schizophrenia lead to an increasing risk of breast cancer (odds ratio [OR] per log-odds increase in schizophrenia risk: 1.049, 95% confidence interval [CI]:1.023-1.075; p = 0.00012; FDR = 0.0017), ovarian cancer (OR, 1.326; 95% CI, 1.267-1.387; p = 0.0007; FDR = 0.0045), and thyroid cancer (OR, 1.575; 95% CI, 1.048-2.365; p = 0.0285; FDR = 0.123). Secondly, we performed a meta-analysis based on the GWAS summary statistics of SCZ and the three significant cancers. Next, we associated genetic variants to genes using two gene mapping strategies: (a) positional mapping based on genomic proximity and (b) expression quantitative trait loci (eQTL) mapping based on gene expression linkage across multiple tissues. As a result, we identified 114 shared loci and 437 shared genes in three groups, respectively. Functional enrichment analysis shows that the most enriched biological pathways are related to epigenetic modification. In addition, we noticed that SCZ would affect the level of thyroid-stimulating hormone (OR, 1.095; 95% CI, 1.006-1.191; p = 0.0354; FDR = 0.177), which may further affect the level of estrogen and the risk of the above three cancers. In conclusion, our findings under the 2SMR assumption provide crucial insights into the risk-increasing effect of SCZ on three cancers' risk. Furthermore, these results may provide insights into understanding the genetic predisposition and underlying biological pathways of comorbid SCZ and cancers.

Locus-level antagonistic selection shaped the polygenic architecture of human complex diseases.

BACKGROUND: We aimed to evaluate the potential role of antagonistic selection in polygenic diseases: if one variant increases the risk of one disease and decreases the risk of another disease, the signals of genetic risk elimination by natural selection will be distorted, which leads to a higher frequency of risk alleles. METHODS: We applied local genetic correlations and transcriptome-wide association studies to identify genomic loci and genes adversely associated with at least two diseases. Then, we used different population genetic metrics to measure the signals of natural selection for these loci and genes. RESULTS: First, we identified 2120 cases of antagonistic pleiotropy (negative local genetic correlation) among 87 diseases in 716 genomic loci (antagonistic loci). Next, by comparing with non-antagonistic loci, we observed that antagonistic loci explained an excess proportion of disease heritability (median 6%), showed enhanced signals of balancing selection, and reduced signals of directional polygenic adaptation. Then, at the gene expression level, we identified 31,991 cases of antagonistic pleiotropy among 98 diseases at 4368 genes. However, evidence of altered signals of selection pressure and heritability distribution at the gene expression level is limited. CONCLUSION: We conclude that antagonistic pleiotropy is widespread among human polygenic diseases, and it has distorted the evolutionary signal and genetic architecture of diseases at the locus level.

Genome-wide identification of the shared genetic basis of cannabis and cigarette smoking and schizophrenia implicates NCAM1 and neuronal abnormality.

OBJECTIVES: Confirming the existence and composition of the shared genetic basis of Schizophrenia and cannabis and cigarette smoking has critical values for the clinical prevention and intervention of psychosis. METHODS: To achieve this goal, we leveraged Genome-Wide summary statistics of Schizophrenia (n = 99,934), cigarette smoking (n = 518,633) and cannabis usage (n = 162,082). We applied Causal Analysis Using Summary Effect Estimates (CAUSE) and genomic structural equation modeling (GenomicSEM) to quantify the contribution of a common genetic factor of cannabis and cigarette smoking and schizophrenia (referred to as SCZ_SMO), then identified genome-wide loci that made up SCZ_SMO. RESULTS: We estimated that SCZ_SMO explained 8.6% of Schizophrenia heritability (Z score <-2.5 in CAUSE, p<10-20 in Genomic SEM). There were 20 independent loci showing association with SCZ_SMO at the genome-wide threshold of p<5 × 10-8. At the top locus on chromosome 11, fine-mapping identified rs7945073 (posterior inclusion probability =0.12, p = 2.24 × 10-32) as the top risk variants. Gene-level association and fine-mapping highlighted NCAM1, PHC2, and SEMA6D as risk genes of SCZ_SMO. Other risk genes were enriched in cortex, neuron, and dendritic spines (adjusted p<0.05). SCZ_SMO showed significant positive correlation (p<10-6) with the genetic risk of attention deficit hyperactivity disorder (r = 0.50), lifestyle problems (r = 0.83), social deprivation (r = 0.58) and all-cause pregnant loss (r = 0.60). CONCLUSION: Our result provided new evidence on the shared genetic basis model for the association between Schizophrenia and smoking and provided genetic and biological insights into their shared mechanism.

An Integrative Segmentation Framework for Cell Nucleus of Fluorescence Microscopy.

Nucleus segmentation of fluorescence microscopy is a critical step in quantifying measurements in cell biology. Automatic and accurate nucleus segmentation has powerful applications in analyzing intrinsic characterization in nucleus morphology. However, existing methods have limited capacity to perform accurate segmentation in challenging samples, such as noisy images and clumped nuclei. In this paper, inspired by the idea of cascaded U-Net (or W-Net) and its remarkable performance improvement in medical image segmentation, we proposed a novel framework called Attention-enhanced Simplified W-Net (ASW-Net), in which a cascade-like structure with between-net connections was used. Results showed that this lightweight model could reach remarkable segmentation performance in the BBBC039 testing set (aggregated Jaccard index, 0.90). In addition, our proposed framework performed better than the state-of-the-art methods in terms of segmentation performance. Moreover, we further explored the effectiveness of our designed network by visualizing the deep features from the network. Notably, our proposed framework is open source.

De novo mutations identified by whole-genome sequencing implicate chromatin modifications in obsessive-compulsive disorder.

Obsessive-compulsive disorder (OCD) is a chronic anxiety disorder with a substantial genetic basis and a broadly undiscovered etiology. Recent studies of de novo mutation (DNM) exome-sequencing studies for OCD have reinforced the hypothesis that rare variation contributes to the risk. We performed, to our knowledge, the first whole-genome sequencing on 53 parent-offspring families with offspring affected with OCD to investigate all rare de novo variants and insertions/deletions. We observed higher mutation rates in promoter-anchored chromatin loops (empirical P = 0.0015) and regions with high frequencies of histone marks (empirical P = 0.0001). Mutations affecting coding regions were significantly enriched within coexpression modules of genes involved in chromatin modification during human brain development. Four genes­SETD5, KDM3B, ASXL3, and FBL­had strong aggregated evidence and functionally converged on transcription's epigenetic regulation, suggesting an important OCD risk mechanism. Our data characterized different genome-wide DNMs and highlighted the contribution of chromatin modification in the etiology of OCD.

Psychological factors as the risk factor of mouth ulcers: A two-sample Mendelian randomization study.

To examine whether psychological traits (PT) had causal effects on Mouth Ulcers (MU), we applied two-sample Mendelian randomization (MR) to genetics association summary statistics of eleven PT and MU. After the adjustment of outlier variants, genetic correlations and multiple testing, well-being (WB) spectrum PT like life satisfactory (odds ratio [OR] = 0.638 per one standard deviation increment of PT score) had protective effects on MU. Reverse WB traits like neuroticism (OR = 1.60) increased the risk of MU. The lack of well-being characteristics may increase the risk of MU, which highlighted the value of preventive oral care for people who have a reverse mental condition.

A Comprehensive Evaluation of Cross-Omics Blood-Based Biomarkers for Neuropsychiatric Disorders.

The identification of peripheral multi-omics biomarkers of brain disorders has long been hindered by insufficient sample size and confounder influence. This study aimed to compare biomarker potential for different molecules and diseases. We leveraged summary statistics of five blood quantitative trait loci studies (N = 1980 to 22,609) and genome-wide association studies (N = 9725 to 500,199) from 14 different brain disorders, such as Schizophrenia (SCZ) and Alzheimer's Disease (AD). We applied summary-based and two-sample Mendelian Randomization to estimate the associations between blood molecules and brain disorders. We identified 524 RNA, 807 methylation sites, 29 proteins, seven cytokines, and 22 metabolites having a significant association with at least one of 14 brain disorders. Simulation analyses indicated that a cross-omics combination of biomarkers had better performance for most disorders, and different disorders could associate with different omics. We identified an 11-methylation-site model for SCZ diagnosis (Area Under Curve, AUC = 0.74) by analyzing selected candidate markers in published datasets (total N = 6098). Moreover, we constructed an 18-methylation-sites model that could predict the prognosis of elders with mild cognitive impairment (hazard ratio = 2.32). We provided an association landscape between blood cross-omic biomarkers and 14 brain disorders as well as a suggestion guide for future clinical discovery and application.

A selection pressure landscape for 870 human polygenic traits.

Characterizing the natural selection of complex traits is important for understanding human evolution and both biological and pathological mechanisms. We leveraged genome-wide summary statistics for 870 polygenic traits and attempted to quantify signals of selection on traits of different forms in European ancestry across four periods in human history and evolution. We found that 88% of these traits underwent polygenic change in the past 2,000-3,000 years. Recent selection was associated with ancient selection signals in the same trait. Traits related to pigmentation, body measurement and nutritional intake exhibited strong selection signals across different time scales. Our findings are limited by our use of exclusively European data and the use of genome-wide association study data, which identify associations between genetic variants and phenotypes that may not be causal. In sum, we provide an overview of signals of selection on human polygenic traits and their characteristics across human evolution, based on a European subset of human genetic diversity. These findings could serve as a foundation for further populational and medical genetic studies.

Label propagation-based semi-supervised feature selection on decoding clinical phenotypes with RNA-seq data.

BACKGROUND: Clinically, behavior, cognitive, and mental functions are affected during the neurodegenerative disease progression. To date, the molecular pathogenesis of these complex disease is still unclear. With the rapid development of sequencing technologies, it is possible to delicately decode the molecular mechanisms corresponding to different clinical phenotypes at the genome-wide transcriptomic level using computational methods. Our previous studies have shown that it is difficult to distinguish disease genes from non-disease genes. Therefore, to precisely explore the molecular pathogenesis under complex clinical phenotypes, it is better to identify biomarkers corresponding to different disease stages or clinical phenotypes. So, in this study, we designed a label propagation-based semi-supervised feature selection approach (LPFS) to prioritize disease-associated genes corresponding to different disease stages or clinical phenotypes. METHODS: In this study, we pioneering put label propagation clustering and feature selection into one framework and proposed label propagation-based semi-supervised feature selection approach. LPFS prioritizes disease genes related to different disease stages or phenotypes through the alternative iteration of label propagation clustering based on sample network and feature selection with gene expression profiles. Then the GO and KEGG pathway enrichment analysis were carried as well as the gene functional analysis to explore molecular mechanisms of specific disease phenotypes, thus to decode the changes in individual behavioral and mental characteristics during neurodegenerative disease progression. RESULTS: Large amounts of experiments were conducted to verify the performance of LPFS with Huntington's gene expression data. Experimental results shown that LPFS performs better in comparison with the-state-of-art methods. GO and KEGG enrichment analysis of key gene sets shown that TGF-beta signaling pathway, cytokine-cytokine receptor interaction, immune response, and inflammatory response were gradually affected during the Huntington's disease progression. In addition, we found that the expression of SLC4A11, ZFP474, AMBP, TOP2A, PBK, CCDC33, APSL, DLGAP5, and Al662270 changed seriously by the development of the disease. CONCLUSIONS: In this study, we designed a label propagation-based semi-supervised feature selection model to precisely selected key genes of different disease phenotypes. We conducted experiments using the model with Huntington's disease mice gene expression data to decode the mechanisms of it. We found many cell types, including astrocyte, microglia, and GABAergic neuron, could be involved in the pathological process.