Biomedical subjects

Yue Li

Publications and source records attributed to Yue Li.

7 recordsLinked to original sources

PubMed2026

Integrated analysis uncovers exogenous induction and molecular regulation of erinacine A accumulation in Hericium erinaceus.

Erinacine A, a cyathane-type diterpenoid mainly from Hericium erinaceus mycelia, exhibits prominent neurotrophic and neuroprotective activities, making it a promising candidate for managing neurodegenerative diseases. However, its low abundance and unclear genetic regulatory mechanisms hinder its application as a nutraceutical. This study aimed to decipher its regulatory mechanisms and enhance production. Four exogenous inducers were screened, with salicylic acid (SA) and ergosterol (ERG) significantly increasing erinacine A content by 62.21% and 146.70% at 20 days, respectively. Transcriptome and WGCNA of inducer-treated sample identified darkorange and magenta modules associated with erinacine A biosynthesis, with the eri gene cluster enriched in the darkorange module and eriG and eriF as hub genes. Forward genetic analysis via QTL mapping of the HeD127 dikaryon population revealed significant phenotypic variation in erinacine A content (0.341-13.085 mg/g) and identified two loci (erA-1 and erA-2) explaining 18.63% of phenotypic variation. Integrating these forward and reverse genetic analyses revealed that salicylic acid and ergosterol synergistically regulate core carbon metabolic pathways to augment acetyl-CoA supply for the mevalonate pathway, suppressed competitive metabolism, enhanced diterpene skeleton construction and structural modification. These results deepen our understanding of the genetic and molecular basis governing accumulation of erinacine A, and facilitate its application in neuroprotective pharmaceuticals.

Diterpenes↗

PubMed2026

Haplotype-resolved 3D genome maps reveal RNAPII-mediated allelic regulation in hybrid rice.

To understand how the two parental genomes coordinate transcription in hybrids, chromatin architecture must be resolved at the haplotype level. Here, using phased Bridge-Linker Hi-C, we reconstructed a haplotype-resolved three-dimensional (3D) genome of the elite hybrid rice (Oryza sativa) line Shanyou 63 (SY63). We identified extensive allele-specific chromatin conformations. Furthermore, we generated allele-resolved RNAPII ChIA-PET maps and phased transcriptomes to explore how chromatin interactions contribute to allelic regulation. Although maternal and paternal homologs share broadly similar chromatin features, we detected widespread haplotype-biased RNAPII binding and chromatin looping at high resolution. These allele-specific RNAPII-mediated contacts were significantly associated with biased expression. Stronger RNAPII binding on one haplotype promoted the formation of long-range regulatory loops with distal genes, thereby contributing to allele-biased transcription at a subset of loci, even when promoter-proximal RNAPII occupancy was comparable between alleles. These results demonstrate that subtle differences in RNAPII engagement and 3D regulatory wiring between parental haplotypes can reshape transcriptional output in hybrids, providing new insights into the mechanisms underlying the allelic regulation of gene expression.

Allele-specific chromatin interactions↗

PubMed2026

DNA-aware evaluation and debiasing of sequence-to-function models.

MOTIVATION: Genome sequence-to-function (S2F) models are widely used to interpret base-resolution functional genomics assays. Most S2F models are trained and evaluated against observed counts and profile-shapes using statistical objectives and fidelity metrics. These choices are well motivated, but they are DNA-independent. At the same time, experimental measurements arise from DNA-dependent assays with distinct characteristics. This mismatch motivates a complementary DNA-aware evaluation of S2F-predicted and experimental functional genomic tracks. RESULTS: We study DNA-dependency of experimental and S2F-predicted tracks using track-conditional genome language models (cgLMs). cgLMs predict masked nucleotides from a conditioning track under controlled DNA visibility. Across ATAC-seq and TF ChIP-seq peaks from GM12878 and K562, cgLM-probing reveals a consistent masked DNA-decodability gap between many experimental and S2F-predicted tracks. In particular, single-task (e.g. BPNet) and multi-task (e.g. AlphaGenome) S2F-predicted tracks enabled cgLMs to recover masked nucleotides with significantly higher accuracy and confidence than matched experimental tracks. Analyses of nonpeak and dinucleotide-shuffled sequences show that this gap is not confined to peaks and is not captured by standard DNA-agnostic profile-shape fidelity metrics alone. ChromBPNet Tn5-denoised predictions were an exception and behaved closer to the experimental regime, suggesting that staged training may reduce the gap. We then convert this diagnostic into a critic-derived objective, DNA-dependency matching (DDM), using a frozen multi-headed cgLM critic. We introduce Critic-Guided Profile-Shape Editing (CGPSE), a preliminary post hoc debiasing framework for frozen S2F models. In GM12878 ATAC-seq, CGPSE partially reduces the masked DNA-decodability gap for AlphaGenome and BPNet predictions, while exposing a tradeoff with profile-shape fidelity. AVAILABILITY AND IMPLEMENTATION: https://github.com/li-lab-mcgill/dna-aware-s2f-eval.

DNA↗

PubMed2025

Sparse polygenic risk score inference with the spike-and-slab LASSO.

MOTIVATION: Large-scale biobanks, with rich phenotypic and genomic data across hundreds of thousands of samples, provide ample opportunities to elucidate the genetics of complex traits and diseases. Consequently, there is growing demand for robust and scalable methods for disease risk prediction from genotype data. Inference in this setting is challenging due to the high-dimensionality of genomic data, especially when coupled with smaller sample sizes. Popular Polygenic Risk Score (PRS) inference methods address this challenge by adopting sparse Bayesian priors or penalized regression techniques, such as the Least Absolute Shrinkage and Selection Operator (LASSO). However, the former class of methods are not as scalable and do not produce exact sparsity, while the latter tends to over-shrink large coefficients. RESULTS: In this study, we present SSLPRS, a novel PRS method based on the Spike-and-Slab LASSO (SSL) prior, which offers a theoretical bridge between the two frameworks. We extend previous work to derive a coordinate-ascent inference algorithm that operates on GWAS summary statistics, which is orders-of-magnitude more efficient than corresponding individual-level-based implementations. To illustrate the statistical properties of the proposed model, we conducted experiments involving nine simulation configurations and nine quantitative phenotypes from the UK Biobank. Our results demonstrate that SSLPRS is competitive with state-of-the-art methods in terms of prediction accuracy and exhibits superior variable selection performance, especially in sparse genetic architectures. In simulations, this translates to upwards of 50% improvement in positive predictive value. In analysis of real phenotypes, we show that selected variants are highly enriched for meaningful genomic annotations and have better replication rates in larger meta-analyses. AVAILABILITY AND IMPLEMENTATION: SSLPRS is available in the open-source package https://github.com/li-lab-mcgill/penprs.

Multifactorial Inheritance↗

PubMed2025

Spatial-Temporal Diversity of Extrachromosomal DNA Shapes Urothelial Carcinoma Evolution and Tumor-Immune Microenvironment.

Extrachromosomal DNA (ecDNA) presents a promising target for cancer therapy; however, its spatial-temporal diversity and influence on tumor evolution and the immune microenvironment remain largely unclear. We apply computational methods to analyze ecDNA from whole-genome sequencing data of 595 urothelial carcinoma (UC) patients. We demonstrate that ecDNA drives clonal evolution through structural rearrangements during malignant transformation and recurrence of UC. This supports a model wherein tumors evolve via the selective expansion of ecDNA-bearing cells. Through multi-regional sampling of tumors, we demonstrate that ecDNA contributes to the evolution of multifocality and increased intratumoral heterogeneity. EcDNA is present in 36% of UC tumors and correlates with an immunosuppressive phenotype and poor prognosis. Single-cell RNA sequencing analyses reveal that ecDNA+ malignant cells exhibit diminished expression of major histocompatibility complex class I molecules, enabling them to evade T-cell immunity. Finally, we show that sequencing of urinary sediment-derived DNA has excellent specificity in detecting ecDNA.

Journal Article↗

PubMed2024

SLC40A1-related hemochromatosis associated with a p.Y333H mutation in mainland China: a pedigree report and literature review.

BACKGROUND: Haemochromatosis is a genetic disease characterized by the excessive deposition of iron in various tissues and organs, eventually results in organ damage including cirrhosis, diabetes, cardiomyopathy, etc. SLC40A1-related haemochromatosis is associated with gain-of-function mutations in the SLC40A1 gene, which encodes ferroportin. While sporadic reports of this condition exist in mainland China, the understanding of the phenotype and genetic pattern associated with the SLC40A1 p.Y333H mutation remains incomplete. CASE PRESENTATION: We report a pedigree with heterozygous p.Y333H mutation in Chinese Han population. The proband is a 64-year-old man complaining of persistent abnormality of liver enzyme levels for 1 year, with a history of knee joint pain, diabetes and skin pigmentation. He displayed markedly elevated serum ferritin level and transferrin saturation. Magnetic resonance imaging showed iron deposition in the liver, spleen, and pancreas, along with cirrhosis and splenomegaly. Whole exome sequencing identified a heterozygous allelic variant c.997T > C (p.Y333H). Genetic screening of family members identified four first-degree relatives and three second-degree relatives having the same mutation. Additional cases with this mutation from two published studies were included. Among the probands and screened relatives, all eight males aged over 30 y had ferritin level > 1000 µg/L, transferrin saturation > 90%. Four patients with organ damage in the present study received therapeutic phlebotomy, alleviating clinical symptoms and improving in transferrin saturation and serum ferritin. CONCLUSIONS: This study reports the largest pedigree with heterozygous SLC40A1 p.Y333H mutation in the Chinese population to date. In Chinese families, males over 30 years old with hemochromatosis due to SLC40A1 p.Y333H mutation exhibit severe iron overload phenotypes.

Humans↗

PubMed2024

Protocol to perform integrative analysis of high-dimensional single-cell multimodal data using an interpretable deep learning technique.

The advent of single-cell multi-omics sequencing technology makes it possible for researchers to leverage multiple modalities for individual cells. Here, we present a protocol to perform integrative analysis of high-dimensional single-cell multimodal data using an interpretable deep learning technique called moETM. We describe steps for data preprocessing, multi-omics integration, inclusion of prior pathway knowledge, and cross-omics imputation. As a demonstration, we used the single-cell multi-omics data collected from bone marrow mononuclear cells (GSE194122) as in our original study. For complete details on the use and execution of this protocol, please refer to Zhou et al.1.

Deep Learning↗