Discovery of antimicrobial peptides in the global microbiome with machine learning.

Novel antibiotics are urgently needed to combat the antibiotic-resistance crisis. We present a machine-learning-based approach to predict antimicrobial peptides (AMPs) within the global microbiome and leverage a vast dataset of 63,410 metagenomes and 87,920 prokaryotic genomes from environmental and host-associated habitats to create the AMPSphere, a comprehensive catalog comprising 863,498 non-redundant peptides, few of which match existing databases. AMPSphere provides insights into the evolutionary origins of peptides, including by duplication or gene truncation of longer sequences, and we observed that AMP production varies by habitat. To validate our predictions, we synthesized and tested 100 AMPs against clinically relevant drug-resistant pathogens and human gut commensals both in vitro and in vivo. A total of 79 peptides were active, with 63 targeting pathogens. These active AMPs exhibited antibacterial activity by disrupting bacterial membranes. In conclusion, our approach identified nearly one million prokaryotic AMP sequences, an open-access resource for antibiotic discovery.

Towards the biogeography of prokaryotic genes.

Microbial genes encode the majority of the functional repertoire of life on earth. However, despite increasing efforts in metagenomic sequencing of various habitats1-3, little is known about the distribution of genes across the global biosphere, with implications for human and planetary health. Here we constructed a non-redundant gene catalogue of 303 million species-level genes (clustered at 95% nucleotide identity) from 13,174 publicly available metagenomes across 14 major habitats and use it to show that most genes are specific to a single habitat. The small fraction of genes found in multiple habitats is enriched in antibiotic-resistance genes and markers for mobile genetic elements. By further clustering these species-level genes into 32 million protein families, we observed that a small fraction of these families contain the majority of the genes (0.6% of families account for 50% of the genes). The majority of species-level genes and protein families are rare. Furthermore, species-level genes, and in particular the rare ones, show low rates of positive (adaptive) selection, supporting a model in which most genetic variability observed within each protein family is neutral or nearly neutral.

STAB2: an updated spatio-temporal cell atlas of the human and mouse brain.

The brain is constituted of heterogeneous types of neuronal and non-neuronal cells, which are organized into distinct anatomical regions, and show precise regulation of gene expression during development, aging and function. In the current database release, STAB2 provides a systematic cellular map of the human and mouse brain by integrating recently published large-scale single-cell and single-nucleus RNA-sequencing datasets from diverse regions and across lifespan. We applied a hierarchical strategy of unsupervised clustering on the integrated single-cell transcriptomic datasets to precisely annotate the cell types and subtypes in the human and mouse brain. Currently, STAB2 includes 71 and 61 different cell subtypes defined in the human and mouse brain, respectively. It covers 63 subregions and 15 developmental stages of human brain, and 38 subregions and 30 developmental stages of mouse brain, generating a comprehensive atlas for exploring spatiotemporal transcriptomic dynamics in the mammalian brain. We also augmented web interfaces for querying and visualizing the gene expression in specific cell types. STAB2 is freely available at https://mai.fudan.edu.cn/stab2.

An augmented Mendelian randomization approach provides causality of brain imaging features on complex traits in a single biobank-scale dataset.

Mendelian randomization (MR) is an effective approach for revealing causal risk factors that underpin complex traits and diseases. While MR has been more widely applied under two-sample settings, it is more promising to be used in one single large cohort given the rise of biobank-scale datasets that simultaneously contain genotype data, brain imaging data, and matched complex traits from the same individual. However, most existing multivariable MR methods have been developed for two-sample setting or a small number of exposures. In this study, we introduce a one-sample multivariable MR method based on partial least squares and Lasso regression (MR-PL). MR-PL is capable of considering the correlation among exposures (e.g., brain imaging features) when the number of exposures is extremely upscaled, while also correcting for winner's curse bias. We performed extensive and systematic simulations, and demonstrated the robustness and reliability of our method. Comprehensive simulations confirmed that MR-PL can generate more precise causal estimates with lower false positive rates than alternative approaches. Finally, we applied MR-PL to the datasets from UK Biobank to reveal the causal effects of 36 white matter tracts on 180 complex traits, and showed putative white matter tracts that are implicated in smoking, blood vascular function-related traits, and eating behaviors.

Deciphering the genetic architecture of human brain structure and function: a brief survey on recent advances of neuroimaging genomics.

Brain imaging genomics is an emerging interdisciplinary field, where integrated analysis of multimodal medical image-derived phenotypes (IDPs) and multi-omics data, bridging the gap between macroscopic brain phenotypes and their cellular and molecular characteristics. This approach aims to better interpret the genetic architecture and molecular mechanisms associated with brain structure, function and clinical outcomes. More recently, the availability of large-scale imaging and multi-omics datasets from the human brain has afforded the opportunity to the discovering of common genetic variants contributing to the structural and functional IDPs of the human brain. By integrative analyses with functional multi-omics data from the human brain, a set of critical genes, functional genomic regions and neuronal cell types have been identified as significantly associated with brain IDPs. Here, we review the recent advances in the methods and applications of multi-omics integration in brain imaging analysis. We highlight the importance of functional genomic datasets in understanding the biological functions of the identified genes and cell types that are associated with brain IDPs. Moreover, we summarize well-known neuroimaging genetics datasets and discuss challenges and future directions in this field.

A survey on computational strategies for genome-resolved gut metagenomics.

Recovering high-quality metagenome-assembled genomes (HQ-MAGs) is critical for exploring microbial compositions and microbe-phenotype associations. However, multiple sequencing platforms and computational tools for this purpose may confuse researchers and thus call for extensive evaluation. Here, we systematically evaluated a total of 40 combinations of popular computational tools and sequencing platforms (i.e. strategies), involving eight assemblers, eight metagenomic binners and four sequencing technologies, including short-, long-read and metaHiC sequencing. We identified the best tools for the individual tasks (e.g. the assembly and binning) and combinations (e.g. generating more HQ-MAGs) depending on the availability of the sequencing data. We found that the combination of the hybrid assemblies and metaHiC-based binning performed best, followed by the hybrid and long-read assemblies. More importantly, both long-read and metaHiC sequencings link more mobile elements and antibiotic resistance genes to bacterial hosts and improve the quality of public human gut reference genomes with 32% (34/105) HQ-MAGs that were either of better quality than those in the Unified Human Gastrointestinal Genome catalog version 2 or novel.

Engineering of cytosine base editors with DNA damage minimization and editing scope diversification.

Cytosine base editors (CBEs), which enable precise C-to-T substitutions, have been restricted by potential safety risks, including DNA off-target edits, RNA off-target edits and additional genotoxicity such as DNA damages induced by double-strand breaks (DSBs). Though DNA and RNA off-target edits have been ameliorated via various strategies, evaluation and minimization of DSB-associated DNA damage risks for most CBEs remain to be resolved. Here we demonstrate that YE1, an engineered CBE variant with minimized DNA and RNA off-target edits, could induce prominent DSB-associated DNA damage risks, manifested as γH2AX accumulation in human cells. We then perform deaminase engineering for two deaminases lamprey LjCDA1 and human APOBEC3A, and generate divergent CBE variants with eliminated DSB-associated DNA damage risks, in addition to minimized DNA/RNA off-target edits. Furthermore, the editing scopes and sequence preferences of APOBEC3A-derived CBEs could be further diversified by internal fusion strategy. Taken together, this study provides updated evaluation platform for DSB-associated DNA damage risks of CBEs and further generates a series of safer toolkits with diversified editing signatures to expand their applications.

MorbidGCN: prediction of multimorbidity with a graph convolutional network based on integration of population phenotypes and disease network.

Exploring multimorbidity relationships among diseases is of great importance for understanding their shared mechanisms, precise diagnosis and treatment. However, the landscape of multimorbidities is still far from complete due to the complex nature of multimorbidity. Although various types of biological data, such as biomolecules and clinical symptoms, have been used to identify multimorbidities, the population phenotype information (e.g. physical activity and diet) remains less explored for multimorbidity. Here, we present a graph convolutional network (GCN) model, named MorbidGCN, for multimorbidity prediction by integrating population phenotypes and disease network. Specifically, MorbidGCN treats the multimorbidity prediction as a missing link prediction problem in the disease network, where a novel feature selection method is embedded to select important phenotypes. Benchmarking results on two large-scale multimorbidity data sets, i.e. the UK Biobank (UKB) and Human Disease Network (HuDiNe) data sets, demonstrate that MorbidGCN outperforms other competitive methods. With MorbidGCN, 9742 and 14 010 novel multimorbidities are identified in the UKB and HuDiNe data sets, respectively. Moreover, we notice that the selected phenotypes that are generally differentially distributed between multimorbidity patients and single-disease patients can help interpret multimorbidities and show potential for prognosis of multimorbidities.

An integrated brain-specific network identifies genes associated with neuropathologic and clinical traits of Alzheimer's disease.

Alzheimer's disease (AD) has a strong genetic predisposition. However, its risk genes remain incompletely identified. We developed an Alzheimer's brain gene network-based approach to predict AD-associated genes by leveraging the functional pattern of known AD-associated genes. Our constructed network outperformed existing networks in predicting AD genes. We then systematically validated the predictions using independent genetic, transcriptomic, proteomic data, neuropathological and clinical data. First, top-ranked genes were enriched in AD-associated pathways. Second, using external gene expression data from the Mount Sinai Brain Bank study, we found that the top-ranked genes were significantly associated with neuropathological and clinical traits, including the Consortium to Establish a Registry for Alzheimer's Disease score, Braak stage score and clinical dementia rating. The analysis of Alzheimer's brain single-cell RNA-seq data revealed cell-type-specific association of predicted genes with early pathology of AD. Third, by interrogating proteomic data in the Religious Orders Study and Memory and Aging Project and Baltimore Longitudinal Study of Aging studies, we observed a significant association of protein expression level with cognitive function and AD clinical severity. The network, method and predictions could become a valuable resource to advance the identification of risk genes for AD.

SemiBin2: self-supervised contrastive learning leads to better MAGs for short- and long-read sequencing.

MOTIVATION: Metagenomic binning methods to reconstruct metagenome-assembled genomes (MAGs) from environmental samples have been widely used in large-scale metagenomic studies. The recently proposed semi-supervised binning method, SemiBin, achieved state-of-the-art binning results in several environments. However, this required annotating contigs, a computationally costly and potentially biased process. RESULTS: We propose SemiBin2, which uses self-supervised learning to learn feature embeddings from the contigs. In simulated and real datasets, we show that self-supervised learning achieves better results than the semi-supervised learning used in SemiBin1 and that SemiBin2 outperforms other state-of-the-art binners. Compared to SemiBin1, SemiBin2 can reconstruct 8.3-21.5% more high-quality bins and requires only 25% of the running time and 11% of peak memory usage in real short-read sequencing samples. To extend SemiBin2 to long-read data, we also propose ensemble-based DBSCAN clustering algorithm, resulting in 13.1-26.3% more high-quality genomes than the second best binner for long-read data. AVAILABILITY AND IMPLEMENTATION: SemiBin2 is available as open source software at https://github.com/BigDataBiology/SemiBin/ and the analysis scripts used in the study can be found at https://github.com/BigDataBiology/SemiBin2_benchmark.

mBodyMap: a curated database for microbes across human body and their associations with health and diseases.

mBodyMap is a curated database for microbes across the human body and their associations with health and diseases. Its primary aim is to promote the reusability of human-associated metagenomic data and assist with the identification of disease-associated microbes by consistently annotating the microbial contents of collected samples using state-of-the-art toolsets and manually curating the meta-data of corresponding human hosts. mBodyMap organizes collected samples based on their association with human diseases and body sites to enable cross-dataset integration and comparison. To help users find microbes of interest and visualize and compare their distributions and abundances/prevalence within different body sites and various diseases, the mBodyMap database is equipped with an intuitive interface and extensive graphical representations of the collected data. So far, it contains a total of 63 148 runs, including 14 401 metagenomes and 48 747 amplicons related to health and 56 human diseases, from within 22 human body sites across 136 projects. Also available in the database are pre-computed abundances and prevalence of 6247 species (belonging to 1645 genera) stratified by body sites and diseases. mBodyMap can be accessed at: https://mbodymap.microbiome.cloud.

GMrepo v2: a curated human gut microbiome database with special focus on disease markers and cross-dataset comparison.

GMrepo (data repository for Gut Microbiota) is a database of curated and consistently annotated human gut metagenomes. Its main purposes are to increase the reusability and accessibility of human gut metagenomic data, and enable cross-project and phenotype comparisons. To achieve these goals, we performed manual curation on the meta-data and organized the datasets in a phenotype-centric manner. GMrepo v2 contains 353 projects and 71,642 runs/samples, which are significantly increased from the previous version. Among these runs/samples, 45,111 and 26,531 were obtained by 16S rRNA amplicon and whole-genome metagenomics sequencing, respectively. We also increased the number of phenotypes from 92 to 133. In addition, we introduced disease-marker identification and cross-project/phenotype comparison. We first identified disease markers between two phenotypes (e.g. health versus diseases) on a per-project basis for selected projects. We then compared the identified markers for each phenotype pair across datasets to facilitate the identification of consistent microbial markers across datasets. Finally, we provided a marker-centric view to allow users to check if a marker has different trends in different diseases. So far, GMrepo includes 592 marker taxa (350 species and 242 genera) for 47 phenotype pairs, identified from 83 selected projects. GMrepo v2 is freely available at: https://gmrepo.humangut.info.

Unnatural Direct Interspecies Electron Transfer Enabled by Living Cell-Cell Click Chemistry.

Direct interspecies electron transfer (DIET) is essential for maintaining the function and stability of anaerobic microbial consortia. However, only limited natural DIET modes have been identified and DIET engineering remains highly challenging. In this study, an unnatural DIET between Shewanella oneidensis MR-1 (SO, electron donating partner) and Rhodopseudomonas palustris (RP, electron accepting partner) was artificially established by a facile living cell-cell click chemistry strategy. By introducing alkyne- or azide-modified monosaccharides onto the cell outer surface of the target species, precise covalent connections between different species in high proximity were realized through a fast click chemistry reaction. Remarkably, upon covalent connection, outer cell surface C-type cytochromes mediated DIET between SO and RP was achieved and identified, although this was never realized naturally. Moreover, this connection directly shifted the natural H2 mediated interspecies electron transfer (MIET) to DIET between SO and RP, which delivered superior interspecies electron exchange efficiency. Therefore, this work demonstrated a naturally unachievable DIET and an unprecedented MIET shift to DIET accomplished by cell-cell distance engineering, offering an efficient and versatile solution for DIET engineering, which extends our understanding of DIET and opens up new avenues for DIET exploration and applications.

Compared to histamine-2 receptor antagonist, proton pump inhibitor induces stronger oral-to-gut microbial transmission and gut microbiome alterations: a randomised controlled trial.

OBJECTIVE: We aim to compare the effects of proton pump inhibitors (PPIs) and histamine-2 receptor antagonists (H2RAs) on the gut microbiota through longitudinal analysis. DESIGN: Healthy volunteers were randomly assigned to receive either PPI (n=23) or H2RA (n=26) daily for seven consecutive days. We collected oral (saliva) and faecal samples before and after the intervention for metagenomic next-generation sequencing. We analysed intervention-induced alterations in the oral and gut microbiome including microbial abundance and growth rates, oral-to-gut transmissions, and compared differences between the PPI and H2RA groups. RESULTS: Both interventions disrupted the gut microbiota, with PPIs demonstrating more pronounced effects. PPI usage led to a significantly higher extent of oral-to-gut transmission and promoted the growth of specific oral microbes in the gut. This led to a significant increase in both the number and total abundance of oral species present in the gut, including the identification of known disease-associated species like Fusobacterium nucleatum and Streptococcus anginosus. Overall, gut microbiome-based machine learning classifiers could accurately distinguish PPI from non-PPI users, achieving an area under the receiver operating characteristic curve (AUROC) of 0.924, in contrast to an AUROC of 0.509 for H2RA versus non-H2RA users. CONCLUSION: Our study provides evidence that PPIs have a greater impact on the gut microbiome and oral-to-gut transmission than H2RAs, shedding light on the mechanism underlying the higher risk of certain diseases associated with prolonged PPI use. TRIAL REGISTRATION NUMBER: ChiCTR2300072310.

nMAGMA: a network-enhanced method for inferring risk genes from GWAS summary statistics and its application to schizophrenia.

MOTIVATION: Annotating genetic variants from summary statistics of genome-wide association studies (GWAS) is crucial for predicting risk genes of various disorders. The multimarker analysis of genomic annotation (MAGMA) is one of the most popular tools for this purpose, where MAGMA aggregates signals of single nucleotide polymorphisms (SNPs) to their nearby genes. In biology, SNPs may also affect genes that are far away in the genome, thus missed by MAGMA. Although different upgrades of MAGMA have been proposed to extend gene-wise variant annotations with more information (e.g. Hi-C or eQTL), the regulatory relationships among genes and the tissue specificity of signals have not been taken into account. RESULTS: We propose a new approach, namely network-enhanced MAGMA (nMAGMA), for gene-wise annotation of variants from GWAS summary statistics. Compared with MAGMA and H-MAGMA, nMAGMA significantly extends the lists of genes that can be annotated to SNPs by integrating local signals, long-range regulation signals (i.e. interactions between distal DNA elements), and tissue-specific gene networks. When applied to schizophrenia (SCZ), nMAGMA is able to detect more risk genes (217% more than MAGMA and 57% more than H-MAGMA) that are involved in SCZ compared with MAGMA and H-MAGMA, and more of nMAGMA results can be validated with known SCZ risk genes. Some disease-related functions (e.g. the ATPase pathway in Cortex) are also uncovered in nMAGMA but not in MAGMA or H-MAGMA. Moreover, nMAGMA provides tissue-specific risk signals, which are useful for understanding disorders with multitissue origins.

Identifying age-specific gene signatures of the human cerebral cortex with joint analysis of transcriptomes and functional connectomes.

The human cerebral cortex undergoes profound structural and functional dynamic variations across the lifespan, whereas the underlying molecular mechanisms remain unclear. Here, with a novel method transcriptome-connectome correlation analysis (TCA), which integrates the brain functional magnetic resonance images and region-specific transcriptomes, we identify age-specific cortex (ASC) gene signatures for adolescence, early adulthood and late adulthood. The ASC gene signatures are significantly correlated with the cortical thickness (P-value <2.00e-3) and myelination (P-value <1.00e-3), two key brain structural features that vary in accordance with brain development. In addition to the molecular underpinning of age-related brain functions, the ASC gene signatures allow delineation of the molecular mechanisms of neuropsychiatric disorders, such as the regulation between ARNT2 and its target gene ETF1 involved in Schizophrenia. We further validate the ASC gene signatures with published gene sets associated with the adult cortex, and confirm the robustness of TCA on other brain image datasets. Availability: All scripts are written in R. Scripts for the TCA method and related statistics result can be freely accessed at https://github.com/Soulnature/TCA. Additional data related to this paper may be requested from the authors.

Disrupted long-range gene regulations elucidate shared tissue-specific mechanisms of neuropsychiatric disorders.

Neurological and psychiatric disorders have overlapped phenotypic profiles, but the underlying tissue-specific functional processes remain largely unknown. In this study, we explore the shared tissue-specificity among 14 neuropsychiatric disorders through the disrupted long-range gene regulations by GWAS-identified regulatory SNPs. Through Hi-C interactions, averagely 38.0% and 17.2% of the intergenic regulatory SNPs can be linked to target protein-coding genes in brain and non-brain tissues, respectively. Interestingly, while the regulatory target genes in the brain tend to enrich in nervous system development related processes, those in the non-brain tissues are inclined to interfere with synapse and neuroinflammation related processes. Compared to psychiatric disorders, neurological disorders present more prominently the neuroinflammatory processes in both brain and non-brain tissues, indicating an intrinsic difference in mechanisms. Through tissue-specific gene regulatory networks, we then constructed disorder similarity networks in two brain and three non-brain tissues, highlighting both known disorder clusters (e.g. the neurodevelopmental disorders) and unexpected disorder clusters (e.g. Parkinson's disease is consistently grouped with psychiatric disorders). We showcase the potential pharmaceutical applications of the small bowel and its disorder clusters, illustrated by the known drug targets NR1I3 and NFACT1, and their small bowel-specific regulatory modules. In conclusion, disrupted long-range gene regulations in both brain and non-brain tissues contribute to the similarity among distinct clusters of neuropsychiatric disorders, and the tissue-specifically shared functions and regulators for disease clusters may provide insights for future therapeutic investigations.

Common and disorder-specific cortical thickness alterations in internalizing, externalizing and thought disorders during early adolescence: an Adolescent Brain and Cognitive Development study.

BACKGROUND: A growing body of neuroimaging studies has reported common neural abnormalities among mental disorders in adults. However, it is unclear whether the distinct disorder-specific mechanisms operate during adolescence despite the overlap among disorders. METHODS: We studied a large cohort of more than 11 000 preadolescent (age 9-10 yr) children from the Adolescent Brain and Cognitive Development cohort. We adopted a regrouping approach to compare cortical thickness (CT) alterations and longitudinal changes between healthy controls (n = 4041) and externalizing (n = 1182), internalizing (n = 1959) and thought disorder (n = 347) groups. Genome-wide association study (GWAS) was performed on regional CT across 4468 unrelated European youth. RESULTS: Youth with externalizing or internalizing disorders exhibited increased regional CT compared with controls. Externalizing (p = 8 × 10-4, Cohen d = 0.10) and internalizing disorders (p = 2 × 10-3, Cohen d = 0.08) shared thicker CT in the left pars opercularis. The somatosensory and the primary auditory cortex were uniquely affected in externalizing disorders, whereas the primary motor cortex and higher-order visual association areas were uniquely affected in internalizing disorders. Only youth with externalizing disorders showed decelerated cortical thinning from age 10-12 years. The GWAS found 59 genome-wide significant associated genetic variants across these regions. Cortical thickness in common regions was associated with glutamatergic neurons, while internalizing-specific regional CT was associated with astrocytes, oligodendrocyte progenitor cells and GABAergic neurons. LIMITATIONS: The sample size of the GWAS was relatively small. CONCLUSION: Our study provides strong evidence for the presence of specificity in CT, developmental trajectories and underlying genetic underpinnings among externalizing and internalizing disorders during early adolescence. Our results support the neurobiological validity of the regrouping approach that could supplement the use of a dimensional approach in future clinical practice.

STAB: a spatio-temporal cell atlas of the human brain.

The human brain is the most complex organ consisting of billions of neuronal and non-neuronal cells that are organized into distinct anatomical and functional regions. Elucidating the cellular and transcriptome architecture underlying the brain is crucial for understanding brain functions and brain disorders. Thanks to the single-cell RNA sequencing technologies, it is becoming possible to dissect the cellular compositions of the brain. Although great effort has been made to explore the transcriptome architecture of the human brain, a comprehensive database with dynamic cellular compositions and molecular characteristics of the human brain during the lifespan is still not available. Here, we present STAB (a Spatio-Temporal cell Atlas of the human Brain), a database consists of single-cell transcriptomes across multiple brain regions and developmental periods. Right now, STAB contains single-cell gene expression profiling of 42 cell subtypes across 20 brain regions and 11 developmental periods. With STAB, the landscape of cell types and their regional heterogeneity and temporal dynamics across the human brain can be clearly seen, which can help to understand both the development of the normal human brain and the etiology of neuropsychiatric disorders. STAB is available at http://stab.comp-sysbio.org.

mMGE: a database for human metagenomic extrachromosomal mobile genetic elements.

Extrachromosomal mobile genetic elements (eMGEs), including phages and plasmids, that can move across different microbes, play important roles in genome evolution and shaping the structure of microbial communities. However, we still know very little about eMGEs, especially their abundances, distributions and putative functions in microbiomes. Thus, a comprehensive description of eMGEs is of great utility. Here we present mMGE, a comprehensive catalog of 517 251 non-redundant eMGEs, including 92 492 plasmids and 424 759 phages, derived from diverse body sites of 66 425 human metagenomic samples. About half the eMGEs could be further grouped into 70 074 clusters using relaxed criteria (referred as to eMGE clusters below). We provide extensive annotations of the identified eMGEs including sequence characteristics, taxonomy affiliation, gene contents and their prokaryotic hosts. We also calculate the prevalence, both within and across samples for each eMGE and eMGE cluster, enabling users to see putative associations of eMGEs with human phenotypes or their distribution preferences. All eMGE records can be browsed or queried in multiple ways, such as eMGE clusters, metagenomic samples and associated hosts. The mMGE is equipped with a user-friendly interface and a BLAST server, facilitating easy access/queries to all its contents easily. mMGE is freely available for academic use at: https://mgedb.comp-sysbio.org.