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.

circ_WASF2 regulates ferroptosis by miR-634/ GPX4 signaling in pancreatic cancer.

PURPOSE: Pancreatic cancer (PC) is one of the most lethal malignant gastrointestinal tumors (GI) characterized by a poor prognosis. Ferroptosis is an emerging programmed cell death that plays an essential role in the progression of various cancers. Ferroptosis is driven by iron-dependent phospholipid peroxidation and is regulated by mitochondrial activity, lipid peroxidation, and reactive oxygen species (ROS). The function and mechanism of ferroptosis in PC need more research. METHODS: The levels of circRNAs, miRNAs, and mRNAs were detected by quantitative real-time polymerase chain reaction (qRT-PCR). Western blot was used for protein detection. CCK8 assays were used to detect cell proliferation. Cell death, lipid peroxidation, ROS, and Fe2+ were detected by indicted kits. Dual-luciferase reporter and RNA pull-down assays were conducted to confirm the interaction between circRNAs, miRNAs, and mRNAs. RESULTS: In this research, we found that circular RNA hsa_circ_0000003(circ_WASF2) was upregulated in pancreatic cancer cells. The silence of circ_WASF2 inhibited cancer proliferation and increased cell death by increasing ferroptosis accompanied by up-regulation of lipid peroxidation, ROS, and Fe2+. Further studies showed that circ_WASF2 could attenuate ferroptosis by targeting miR-634 and the downstream glutathione peroxidase 4 (GPX4). GPX4 has been well-reported as a central factor in ferroptosis. Our research revealed a new pathway for regulating ferroptosis in PC. CONCLUSION: In summary, we have determined that circ_WASF2/miR-634/GPX4 contributed to ferroptosis-induced cell death, and provided a possible therapeutic target in PC.

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. Here, 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 via 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 would extend our understanding of DIET and open up new avenue for DIET exploration and applications.

Evaluation of models to predict prognosis in patients with advanced hepatocellular carcinoma treated with TACE combined with apatinib.

BACKGROUND: The HAP, Six-and-Twelve, Up to Seven, and ALBI scores have been substantiated as reliable prognostic markers in patients presenting with intermediate and advanced hepatocellular carcinoma (HCC) undergoing transarterial chemoembolization (TACE) treatment. Given this premise, our research aims to assess the predictive efficacy of these models in patients with intermediate and advanced HCC receiving a combination of TACE and Apatinib. Additionally, we have conducted a meticulous comparative analysis of these four scoring systems to discern their respective predictive capacities and efficacies in combined therapy. METHODS: Performing a retrospective analysis on the clinical data from 200 patients with intermediate and advanced HCC, we studied those who received TACE combined with Apatinib at the First Affiliated Hospital of the University of Science and Technology of China between June 2018 and December 2022. To identify the factors affecting survival, the study performed univariate and multivariate Cox regression analyses, with calculations of four different scores: HAP, Six-and-Twelve, Up to Seven, and ALBI. Lastly, Harrell's C-index was employed to compare the prognostic abilities of these scores. RESULTS: Cox proportional hazards model results revealed that the ALBI score, presence of portal vein tumor thrombus (PVTT, )and tumor size are independent determinants of prognostic survival. The Kaplan-Meier analyses showed significant differences in survival rates among patients classified by the HAP, Six-and-Twelve, Up to Seven, and ALBI scoring methods. Of the evaluated systems, the HAP scoring demonstrated greater prognostic precision, with a Harrell's C-index of 0.742, surpassing the alternative models (P < 0.05). In addition, an analysis of the area under the AU-ROC curve confirms the remarkable superiority of the HAP score in predicting short-term survival outcomes. CONCLUSION: Our study confirms the predictive value of HAP, Six-and-Twelve, Up to Seven, and ALBI scores in intermediate to advanced Hepatocellular Carcinoma (HCC) patients receiving combined Transarterial Chemoembolization (TACE) and Apatinib therapy. Notably, the HAP model excels in predicting outcomes for this specific HCC subgroup.

MiR-221-5p regulates blood-brain barrier dysfunction through the angiopoietin-1/-2/Tie-2 signaling axis after subarachnoid hemorrhage.

AIM: To explore the potential role of microRNA miR-221-5p on the angiopoietin-1 (Ang-1)/Ang-2/Tie-2 signaling axis after subarachnoid hemorrhage (SAH) in a rat model. METHODS: Aspects of the rat's behavior were measured using the Kaoutzanis scoring system to test neurological responses. This included feeding behavior, body contraction, motor, and eye-opening responses. Brain sections were studied using transmission electron microscopy and Evans blue extravasation. Levels of Ang-1, Ang-2, and Tie-2 were determined by Western blot, while miR-221-5p was quantified using stem-loop real-time quantitative PCR (RT-qPCR). RESULTS: The SAH group responded worse to the neurological response test than the sham-operated group. The intercellular space was widened in the SAH group, but not in the sham-operated group. Evans blue dye leaked significantly more into brain tissue cells of the SAH group. Stem-loop qRT-PCR showed elevated miR-221-5p levels. Additionally, Ang-1 and Tie-2 were reduced but Ang-2 expression was increased after SAH. This led to a significant reduction of the Ang-1/Ang-2 ratio in the brain tissue, which was associated with the destruction of the blood-brain barrier. CONCLUSION: The data indicate that miR-221-5p might regulate blood-brain barrier dysfunction through the Ang-1/Ang-2/Tie-2 signaling axis, suggesting that it should be further investigated as a potential novel biomarker.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing.

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Associations of fatty acids with the risk of biliary tract calculus and inflammation: a Mendelian randomization study.

BACKGROUND: The presence of gallstones in both the gallbladder and bile ducts is referred to as cholelithiasis. The prevalence of cholecystolithiasis and bile duct stones differs. Observational and Mendelian randomization (MR) studies have elucidated the significant contributing role of numerous fatty acids (FAs) in the development of cholelithiasis. Despite numerous studies about cholelithiasis, evidence on the relationship between serum FA levels and cholecystolithiasis, as well as bile duct stones with or without inflammation, remains insufficient. METHODS: A two-sample MR study was designed to clarify the impact of serum FA levels on various bile duct inflammatory diseases. The summary statistics of single nucleotide polymorphisms (SNPs) associated with fatty acids were obtained from the UK Biobank (UKB) and included data from 114,999 participants. The researchers obtained GWAS summary statistics for cholecystolithiasis and bile duct stones in 463,010 and 361,194 European participants, including cases with and without inflammation. No sample overlap between the exposure and outcome was verified through the "mr-lap" package. The SNPs were screened to identify instrumental variables (IVs). Cochran's Q test was applied for heterogeneity assessment. Inverse variance weighting (IVW) (fixed effects or random effects), MR-Egger regression and weighted median methods were used for MR. Multivariable MR was applied to determine the direct effect of each exposure on the outcome. A false discovery rate (FDR) was applied to adjust for multiple testing correction based on the Benjamini-Hochberg method. Finally, the FinnGen Consortium was used to validate some results. RESULTS: The overall concentration of polyunsaturated fatty acids (PUFAs) in the serum was negatively associated with the risk of calculus of the gallbladder with acute cholecystitis (IVW, OR = 0.996, P = 0.038, CI 0.992-0.999; weighted median, OR = 0.995, P = 0.025, CI 0.991-0.999). The percentage of PUFAs to total monounsaturated fatty acids(MUFAs) (IVW, OR = 0.998, P = 0.045, CI 0.997-0.999) and the percentage of PUFAs to total FAs (IVW, OR = 0.997, P = 0.025, CI 0.995-0.999) had a protective role against cholecystitis. The percentage of PUFAs to total FAs had a protective role against calculus of the gallbladder without cholecystitis (IVW, OR = 0.995, P = 0.026, CI 0.990-0.999; MR Egger, OR = 0.99, P = 0.03, CI 0.982-0.998; weighted median, OR = 0.991, P = 5.41e-06, CI 0.988-0.995). Conversely, the percentage of MUFAs to total FAs increased the risk for cholecystitis (IVW, OR = 1.001, P = 0.034, CI 1.0001-1.002). However, there were no causal effects of the above exposures on the outcomes through multivariable MR and multiple testing correction. Finally, the causal effects of the above exposures on cholecystitis were validated in the FinnGen Consortium, which suggested that the percentage of PUFAs to total FAs (IVW, OR = 0.744, P = 0.021, CI 0.579-0.957) had a protective role against cholecystitis. CONCLUSION: These Mendelian randomization findings suggested that more attention should be focused on people who have low serum PUFA levels, which may have a potential role in the occurrence of calculus of the gallbladder or cholecystitis rather than calculus of the bile duct without cholangitis or cholecystitis.

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.

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.

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.

Characterization of the complete mitochondrial genome and phylogenetic analysis of bean thrips Megalurothrips usitatus (Bagnall, 1913) (Thysanoptera: thripidae).

Megalurothrips usitatus (Bagnall, 1913) (Thysanoptera: Thripidae) is a widely distributed pest in Asia that primarily affects the production of snap beans and cowpea. The complete mitochondrial genome of Megalurothrips usitatus has been sequenced and annotated in this study, which is 17,209 bp long and contains 13 protein-coding genes (PCGs), two rRNAs, and 22 tRNA genes. Most of the protein-coding genes (PCGs) start with ATG except ND4 using TTG. Meanwhile, eight PCGs stop with TAA, four PCGs have an incomplete stop codon, and the gene Cytb ends with TAG. Phylogenetic analysis showed that M. usitatus is closely related to Frankliniella intonsa and F. occidentalis, providing a basis for the study of the mitochondrial evolution of Thripinae.

Gut microbiota impacts bone via Bacteroides vulgatus-valeric acid-related pathways.

Although the gut microbiota has been reported to influence osteoporosis risk, the individual species involved, and underlying mechanisms, remain largely unknown. We performed integrative analyses in a Chinese cohort of peri-/post-menopausal women with metagenomics/targeted metabolomics/whole-genome sequencing to identify novel microbiome-related biomarkers for bone health. Bacteroides vulgatus was found to be negatively associated with bone mineral density (BMD), which was validated in US white people. Serum valeric acid (VA), a microbiota derived metabolite, was positively associated with BMD and causally downregulated by B. vulgatus. Ovariectomized mice fed B. vulgatus demonstrated increased bone resorption and poorer bone micro-structure, while those fed VA demonstrated reduced bone resorption and better bone micro-structure. VA suppressed RELA protein production (pro-inflammatory), and enhanced IL10 mRNA expression (anti-inflammatory), leading to suppressed maturation of osteoclast-like cells and enhanced maturation of osteoblasts in vitro. The findings suggest that B. vulgatus and VA may represent promising targets for osteoporosis prevention/treatment.

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.

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.

Computational exploration of the global microbiome for antibiotic discovery.

Novel antibiotics are urgently needed to combat the antibiotic-resistance crisis. We present a machine learning-based approach to predict prokaryotic antimicrobial peptides (AMPs) by leveraging a vast dataset of 63,410 metagenomes and 87,920 microbial genomes. This led to the creation of AMPSphere, a comprehensive catalog comprising 863,498 non-redundant peptides, the majority of which were previously unknown. We observed that AMP production varies by habitat, with animal-associated samples displaying the highest proportion of AMPs compared to other habitats. Furthermore, within different human-associated microbiota, strain-level differences were evident. To validate our predictions, we synthesized and experimentally tested 50 AMPs, demonstrating their efficacy against clinically relevant drug-resistant pathogens both in vitro and in vivo. These AMPs exhibited antibacterial activity by targeting the bacterial membrane. Additionally, AMPSphere provides valuable insights into the evolutionary origins of peptides. In conclusion, our approach identified AMP sequences within prokaryotic microbiomes, opening up new avenues for the discovery of antibiotics.

Prospective Exploratory Study of the Effects of Postoperative Endocrine Medication on the Endometrium in Breast Cancer Patients.

Objective: This study aimed to investigate the pathological effects of long-term postoperative endocrine medication on the endometrium in breast cancer patients. Methods: Data of 99 patients with primary breast cancer who underwent hysteroscopy and obtained endometrial biopsy from 1 June 2018 to 31 December 2021 at the Beijing Obstetrics and Gynecology Hospital, Capital Medical University. Beijing Maternal and Child Health Care Hospital were prospectively collected. Results: Hysteroscopy was performed followed by endometrial histopathological examination in 99 breast cancer patients, including 44 taking tamoxifen (TAM), 26 taking other endocrine drugs, and 29 not taking endocrine drugs. The endometrial thickness in the TAM group was significantly higher than that in the other endocrine drug groups and the group not taking endocrine drugs (p=0.017). The receiver operating characteristic curves for the abnormal premenopausal endometrial thickening were plotted in this study; an endometrial thickness of 15.5 mm seen on ultrasound could be used as the most accurate ultrasound diagnostic threshold for the diagnosis of abnormal premenopausal endometrial hyperplasia, with an area under the curve of 0.888 (95% CI: 0.716, 1.000), a sensitivity of 100%, and a specificity of 75%, which was consistent with the results of our previous retrospective study. An endometrial thickness of ≥5 mm in postmenopausal women with breast cancer taking TAM was still used as the cut-off value for routine ultrasound diagnosis of abnormal postmenopausal endometrial hyperplasia. Conclusion: An ultrasound endometrial thickness (proliferative phase) of >15 mm in premenopausal patients can be used as the most accurate ultrasound diagnostic threshold for the diagnosis of abnormal endometrial hyperplasia. After menopause, an ultrasound endometrial thickness of ≥5 mm is still used as the criterion for diagnosing abnormal endometrial hyperplasia. Older patients should be monitored for signs of vaginal bleeding and fluid discharge, and hysteroscopy should be performed if necessary to ascertain the endometrial condition.

Enterotypes of the human gut mycobiome.

BACKGROUND: The fungal component of the human gut microbiome, also known as the mycobiome, plays a vital role in intestinal ecology and human health. However, the overall structure of the gut mycobiome as well as the inter-individual variations in fungal composition remains largely unknown. In this study, we collected a total of 3363 fungal sequencing samples from 16 cohorts across three continents, including 572 newly profiled samples from China. RESULTS: We identify and characterize four mycobiome enterotypes using ITS profiling of 3363 samples from 16 cohorts. These enterotypes exhibit stability across populations and geographical locations and significant correlation with bacterial enterotypes. Particularly, we notice that fungal enterotypes have a strong age preference, where the enterotype dominated by Candida (i.e., Can_type enterotype) is enriched in the elderly population and confers an increased risk of multiple diseases associated with a compromised intestinal barrier. In addition, bidirectional mediation analysis reveals that the fungi-contributed aerobic respiration pathway associated with the Can_type enterotype might mediate the association between the compromised intestinal barrier and aging. CONCLUSIONS: We show that the human gut mycobiome has stable compositional patterns across individuals and significantly correlates with multiple host factors, such as diseases and host age. Video Abstract.

Prioritizing genes associated with brain disorders by leveraging enhancer-promoter interactions in diverse neural cells and tissues.

BACKGROUND: Prioritizing genes that underlie complex brain disorders poses a considerable challenge. Despite previous studies have found that they shared symptoms and heterogeneity, it remained difficult to systematically identify the risk genes associated with them. METHODS: By using the CAGE (Cap Analysis of Gene Expression) read alignment files for 439 human cell and tissue types (including primary cells, tissues and cell lines) from FANTOM5 project, we predicted enhancer-promoter interactions (EPIs) of 439 cell and tissue types in human, and examined their reliability. Then we evaluated the genetic heritability of 17 diverse brain disorders and behavioral-cognitive phenotypes in each neural cell type, brain region, and developmental stage. Furthermore, we prioritized genes associated with brain disorders and phenotypes by leveraging the EPIs in each neural cell and tissue type, and analyzed their pleiotropy and functionality for different categories of disorders and phenotypes. Finally, we characterized the spatiotemporal expression dynamics of these associated genes in cells and tissues. RESULTS: We found that identified EPIs showed activity specificity and network aggregation in cell and tissue types, and enriched TF binding in neural cells played key roles in synaptic plasticity and nerve cell development, i.e., EGR1 and SOX family. We also discovered that most neurological disorders exhibit heritability enrichment in neural stem cells and astrocytes, while psychiatric disorders and behavioral-cognitive phenotypes exhibit enrichment in neurons. Furthermore, our identified genes recapitulated well-known risk genes, which exhibited widespread pleiotropy between psychiatric disorders and behavioral-cognitive phenotypes (i.e., FOXP2), and indicated expression specificity in neural cell types, brain regions, and developmental stages associated with disorders and phenotypes. Importantly, we showed the potential associations of brain disorders with brain regions and developmental stages that have not been well studied. CONCLUSIONS: Overall, our study characterized the gene-enhancer regulatory networks and genetic mechanisms in the human neural cells and tissues, and illustrated the value of reanalysis of publicly available genomic datasets.

Long-Read Sequencing Reveals Extensive DNA Methylations in Human Gut Phagenome Contributed by Prevalently Phage-Encoded Methyltransferases.

DNA methylation plays a crucial role in the survival of bacteriophages (phages), yet the understanding of their genome methylation remains limited. In this study, DNA methylation patterns are analyzed in 8848 metagenome-assembled high-quality phages from 104 fecal samples using single-molecule real-time sequencing. The results demonstrate that 97.60% of gut phages exhibit methylation, with certain factors correlating with methylation densities. Phages with higher methylation densities appear to have potential viability advantages. Strikingly, more than one-third of the phages possess their own DNA methyltransferases (MTases). Increased MTase copies are associated with higher genome methylation densities, specific methylation motifs, and elevated prevalence of certain phage groups. Notably, the majority of these MTases share close homology with those encoded by gut bacteria, suggesting their exchange during phage-bacterium interactions. Furthermore, these MTases can be employed to accurately predict phage-host relationships. Overall, the findings indicate the widespread utilization of DNA methylation by gut DNA phages as an evasion mechanism against host defense systems, with a substantial contribution from phage-encoded MTases.

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.