Biomedical subjects

Rui Liu

Publications and source records attributed to Rui Liu.

8 recordsLinked to original sources

PubMed2026

Toxoplasma gondii infection disrupts secondary bile acid transformation in feline gut microbiota.

UNLABELLED: Bile acid (BA) transformation relies on gut microbiota and is vulnerable to Toxoplasma gondii infection, yet feline microbial BA-transforming capacity upon toxoplasmosis remains unclear. Here, we constructed a catalog of 2,474 nonredundant feline gut microbial genomes and integrated serum metabolomic data to verify BA transformation alterations. The results revealed that the feline gut microbiome harbored widespread genetic potential for BA transformation but lacked a complete 7α-dehydroxylation pathway due to the absence of the key gene baiE. The BA transformation-related genomes (2,045 in total) were predominantly from the phyla Bacillota_A and Actinomycetota, among which only 37 encoded baiB, all belonging to Bacillota_A. The distribution of BA transformation-related genes varied across intestinal regions: genes encoding 7α-HSDH were primarily enriched in the small intestine, whereas genes encoding 3α-HSDH, baiCD, and baiH were more abundant in the large intestine. Additionally, the abundance of genes encoding BSH and 3α-HSDH increased significantly in the small intestine on day 3 post-infection, accompanied by increases in the phylum Bacillota_C and genera such as Blautia_A, Enterococcus_E, and Ligilactobacillus. Serum metabolomics revealed a significant increase in cholesterol levels post-infection, supporting the impact of T. gondii infection on intestinal BA transformation. These findings illustrated that the feline gut microbiota played an important role in BA transformation and that T. gondii infection disrupted the microbial potential for secondary BA transformation. This study provided new insights into gut microbiota-associated metabolic perturbations during feline toxoplasmosis. IMPORTANCE: Bile acid (BA) transformation plays a critical role in host metabolism and immune regulation. Although studies on BA transformation are increasing, the capacity for BA transformation within the feline gut microbiota and the impact of Toxoplasma gondii infection on this capacity remain unclear. To bridge this gap, we constructed a catalog of 2,474 nonredundant feline gut microbial genomes and integrated serum metabolomic data to verify BA transformation alterations. Our findings revealed that the feline gut microbiome lacked a complete 7α-dehydroxylation pathway, and the specific functions involved in BA transformation may differ between the small and large intestines. Furthermore, integrated metagenomic and serum metabolomic analyses suggested that T. gondii infection disrupted BA transformation capacity in the small intestine. This study provided new insights into gut microbiota-associated metabolic perturbations during feline toxoplasmosis.

Toxoplasma gondii↗

PubMed2026

Interkingdom remodeling of the intestinal bacteriome and virome during Toxoplasma gondii infection in rats.

Toxoplasma gondii infection is associated with intestinal microbiome disruption, but its effects on genome-resolved bacterial populations, the gut virome, and bacteriome-virome relationships remain poorly understood. Using previously generated shotgun metagenomic datasets from 36 intestinal samples collected from 18 Sprague-Dawley rats across control, acute, and chronic infection groups, we reconstructed 294 quality-filtered, non-redundant bacterial metagenome-assembled genomes (MAGs) and identified 899 medium-to-high-quality viral operational taxonomic units (vOTUs) from assembled metagenomic contigs. Infection was associated with reduced bacterial richness in the small intestine during both acute and chronic stages and lower Shannon diversity during chronic infection. In contrast, large-intestinal α-diversity remained stable despite significant compositional reorganization. Taxonomic changes included increased Lactobacillus intestinalis, Limosilactobacillus reuteri, and Prevotella sp900547005, together with decreased Rothia sp002492045 and Akkermansia muciniphila. Functional profiling revealed region- and stage-specific changes in predicted bacterial metabolic potential, including reduced energy-related pathways and carbohydrate-active enzyme abundance. The virome also showed significant compositional changes in both intestinal regions. Quimbyviridae and Podoviridae_crAss-like viruses decreased in the small intestine during chronic infection, while Quimbyviridae, Flandersviridae, and Podoviridae_crAss-like viruses showed stage-specific decreases in the large intestine. Predicted bacterial hosts were assigned to 48.39% of vOTUs, with Lachnospiraceae and Ruminococcaceae being the most frequently linked families. Trans-kingdom networks further revealed region-specific positive and negative abundance correlations between bacterial and viral taxa. These findings extend previous microbiota-metabolome observations by integrating genome-resolved bacteriome analysis with contig-based virome profiling, providing a foundation for future mechanistic studies of toxoplasmosis-associated microbiome remodeling.

Gut virome↗

PubMed2026

Genome-resolved analysis reveals disruption of gut microbial vitamin B and K2 biosynthesis during Toxoplasma gondii infection in mice.

UNLABELLED: Toxoplasma gondii infection remodels the gut microbiome, yet its impact on microbial vitamin biosynthetic potential and host redox metabolism remains unclear. Here, we integrated mouse gut metagenomes with publicly available metagenome-assembled genomes (MAGs) to construct a genome-resolved atlas of B-vitamin and vitamin K2 biosynthesis. From 45,697 MAGs, we curated 4,771 representative genomes, of which 2,682 met high-quality criteria (completeness ≥90%, contamination <5%). Functional annotation identified 229,717 vitamin-related genes corresponding to 177 Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs across de novo pathways for eight B vitamins, thiamine (B1), riboflavin (B2), niacin (B3), pantothenate (B5), pyridoxine (B6), biotin (B7), folate (B9), cobalamin (B12), and vitamin K2. Among the high-quality genomes, 1,665 encoded complete de novo pathways for at least one vitamin, highlighting functional specialization and community-level complementarity. Transcripts per million-normalized metagenomic read counts revealed significant differences in KEGG ortholog abundances across six of the nine vitamin pathways. Reanalysis of metagenomic data from infected mice (acute, chronic, and control; n = 10 per group) revealed a stage-dependent reduction in α-diversity of vitamin biosynthesis pathways during acute infection, and a clear β-diversity separation from chronic and control groups. Core niacin biosynthesis genes (nadB, nadA, nadC) displayed phylum-specific redistribution, indicating selective remodeling of microbial NAD+ precursor production under infection-induced metabolic stress. These results suggest that T. gondii infection disrupts cooperative vitamin biosynthetic networks while specifically modulating niacin pathways linked to host NAD+ metabolism. IMPORTANCE: Gut microbes can synthesize essential vitamins, but how infection alters this function is poorly understood. By integrating mouse gut metagenomes with genome-resolved microbial data, we show that Toxoplasma gondii infection reshapes the vitamin biosynthetic potential of the gut microbiome in a stage-dependent manner. Acute infection reduces the diversity of vitamin biosynthesis pathways and shifts the taxonomic distribution of key niacin biosynthesis genes involved in microbial NAD+ precursor production. These findings identify vitamin metabolism, especially niacin-related pathways, as a sensitive functional axis of microbiome remodeling during infection. Our work links microbial taxonomic changes to functional metabolic consequences and suggests that microbiome-mediated regulation of NAD+-related metabolism may contribute to host redox adaptation during T. gondii infection.

B vitamins↗

PubMed2026

Chromosome-Scale Genome Analysis Reveals Locus-Specific Disruption of the Citrinin-Associated Region in a Furu-Derived Monascus ruber Strain BC20.

Monascus species are widely used in traditional fermented foods for pigment and flavor formation, but citrinin contamination remains a major safety concern that limits broader food applications. Therefore, this study aimed to evaluate the citrinin risk of a furu-derived Monascus ruber strain, BC20, by integrating phenotypic screening across food-relevant matrices with genome-resolved analysis. After 14 days of cultivation across eight matrices, including fungal media as well as dairy-, cereal-, and bran-based substrates, citrinin was not detected by immunoaffinity cleanup combined with HPLC-FLD (LOD, 4 μg/kg; LOQ, 12 μg/kg). To investigate the genetic basis of this phenotype, we generated a chromosome-scale genome assembly for BC20 and conducted comparative analyses across a total of 19 Monascus genomes. ANI analysis and phylogenomic inference consistently placed BC20 within the ruber-pilosus clade. Comparative synteny analysis showed that the citrinin-associated locus in BC20 no longer retained an intact cluster configuration but instead exhibited a remnant-locus architecture, and similar patterns were also observed in several related genomes from the same clade. By contrast, the monacolin K (mk) locus remained syntenically conserved in BC20, supporting locus-specific structural disturbance rather than assembly-derived pseudo-absence. Additionally, its antifungal susceptibility was determined. Overall, BC20 represents a M. ruber candidate strain with undetectable citrinin, and this study provides a practical analytical framework for citrinin risk screening in food-related Monascus isolates.

biosynthetic gene cluster↗

PubMed2026

Polygenic risk scores in major depressive disorder: A systematic review across diagnostic, treatment, course/severity, and subtype domains.

BACKGROUND: Major depressive disorder (MDD) is heterogeneous across diagnostic, treatment-related, course/severity, and subtype domains. Polygenic risk score (PRS) studies have examined these domains, but differences in PRS sources, samples, methods, and endpoint definitions have fragmented the evidence. We synthesised findings and examined potential contributors to heterogeneity. METHODS: PubMed/MEDLINE, Embase, PsycINFO, and Web of Science were searched for studies published from January 2016 through 25 November 2025. Result records were synthesised using SWiM, and certainty was assessed with an adapted GRADE framework. RESULTS: Sixty studies contributed 493 retained records; 450 were descriptively classified as positive, null, or reverse, although records were not independent. Positive findings accounted for 44/56 diagnostic, 61/273 treatment-related, 64/100 course/severity, and 14/21 subtype records. For MDD/depression-derived PRSs and case-control MDD status, all 10 contributing studies showed higher liability in cases (exploratory exact sign test p = 0.002; FDR q = 0.004). The same PRS group showed positive findings for overall depressive symptom severity (14/18), although the study-level test was imprecise (5/5 studies; p = 0.063). Pharmacological response/remission findings for these PRSs were mostly null or directionally mixed (10 positive, 18 null, and 9 reverse). Treatment-resistant depression (TRD) findings differed by operational definition. Atypical and psychotic subtype signals arose mainly from single-study PRS and endpoint contrasts. CONCLUSIONS: PRS evidence was clearest for MDD diagnostic status and showed a tentative pattern for overall symptom burden. Treatment and subtype findings were less consistent or less replicated. Larger, ancestrally diverse studies with standardised endpoints and transparent PRS methods are needed.

Humans↗

PubMed2026

CTRP9 ameliorates heart failure with preserved ejection fraction by regulating lipid metabolism.

BACKGROUND: Heart failure with preserved ejection fraction (HFpEF) is a major clinical challenge, with cardiac lipotoxicity emerging as a key driver of disease progression. Despite CTRP9’s role in lipid metabolism and cardioprotective properties, its therapeutic potential in HFpEF remains unexplored. This study aimed to investigate whether CTRP9 ameliorates HFpEF by regulating cardiac lipid metabolism and to identify the underlying molecular mechanisms. METHODS: In the established two-hit HFpEF mouse model (induced by a high-fat diet and L-NAME), the mice were treated with either CTRP9 or saline. Cardiac function was evaluated by echocardiography, while hypertrophy, fibrosis, and lipid accumulation were assessed using histology and molecular assays. Proteomic sequencing was further employed to identify downstream targets of CTRP9. RESULTS: CTRP9 treatment significantly improved diastolic function and attenuated cardiac hypertrophy and fibrosis in HFpEF mice. Myocardial lipid accumulation was substantially reduced, accompanied by enhanced fatty acid oxidation. Proteomic analysis identified GPD1 as a key downstream target upregulated by CTRP9. Cardiac-specific knockdown of GPD1 partly abolished the therapeutic benefits of CTRP9. CONCLUSION: Our data suggest that CTRP9 ameliorates HFpEF through GPD1-mediated regulation of cardiac lipid metabolism, identifying the CTRP9-GPD1 axis as a promising therapeutic target for HFpEF.

Animals↗

PubMed2025

Significance in scale space for Hi-C data.

MOTIVATION: Hi-C technology has been developed to profile genome-wide chromosome conformation. So far Hi-C data have been generated from a large compendium of different cell types and different tissue types. Among different chromatin conformation units, chromatin loops were found to play a key role in gene regulation across different cell types. While many different loop calling algorithms have been developed, most loop callers identified shared loops as opposed to cell-type-specific loops. RESULTS: We propose SSSHiC, a new loop calling algorithm based on significance in scale space, which can be used to understand data at different levels of resolution. By applying SSSHiC to neuronal and glial Hi-C data, we detected more loops that are potentially engaged in cell-type-specific gene regulation. Compared with other loop callers, such as Mustache, these loops were more frequently anchored to gene promoters of cellular marker genes and had better APA scores. Therefore, our results suggest that SSSHiC can effectively capture loops that contain more gene regulatory information. AVAILABILITY AND IMPLEMENTATION: The Hi-C data used in this study can be accessed through the PsychENCODE Knowledge Portal at https://www.synapse.org/#! Synapse: syn21760712. The code utilized for Curvature SSS cited in this study is available at https://github.com/jsmarron/MarronMatlabSoftware/blob/master/Matlab9/Matlab9Combined.zip. All custom code used in this research can be found in the GitHub repository: https://github.com/jerryliu01998/HiC. The code has also been submitted to Code Ocean with the doi: 10.24433/CO.1912913.v1.

Algorithms↗

PubMed2025

ESMpHLA: Evolutionary Scale Model-Based Deep Learning Prediction of HLA Class I Binding Peptides.

The recognition of endogenous peptides by HLA class I plays a crucial role in CD8+ T cell immune responses and human adaptive cell immune. Thus, the prediction of HLA class I-peptide binding affinities is always the core issue for the research of immune recognition and vaccine development. In this study, an evolutionary scale model (ESM) combined with parallel CNN blocks and a cross attention mechanism was used to construct a novel ESMpHLA model for predicting HLA class I binding peptides. Based on the 91,560 binding peptides of 41 HLA-A alleles, 56,731 of 50 HLA-B alleles and 2444 of 10 HLA-C alleles, the ESMpHLA model was successfully established and achieved satisfying prediction performances with the overall accuracy and AUC values of 0.874 and 0.938 for the test dataset. The results indicate that the ESMpHLA model performs well in dealing with different HLA class I 2-field alleles as well as the peptides with different lengths. Then, the generalisation ability of the ESMpHLA model was validated by an independent test dataset compiled from recent IEDB weekly benchmark datasets. The results showed that the ESMpHLA model achieved the highest ROC-AUC and PR-AUC values when compared with the latest BVMHC, CapsNet-MHC, STMHCpan and BVLSTM models. In addition, two ensemble models were also established by integrating the above 5 deep learning models using soft-voting and hard-voting strategies.

Humans↗