Biomedical subjects

Huan Liu

Publications and source records attributed to Huan Liu.

11 recordsLinked to original sources

PubMed2026

Adolescent depression as a systemic multimorbidity catalyst: integrated genetic and metabolic pathway analysis.

BACKGROUND: Although adolescent depression has been linked to individual chronic conditions, its broader role in shaping multimorbidity risk remains understudied. METHODS: A total of 87,562 UK Biobank participants were included, of whom 18,851 had documented adolescent depression. Cox proportional hazards models were applied to evaluate associations between adolescent depression and 24 chronic diseases, followed by stratified analyses by sex and age. Two-sample Mendelian randomization (MR) was then conducted to infer causality for diseases showing significant associations. Genomic colocalization analyses were performed using relevant GWAS data to identify shared causal variants. Mediation analyses were performed to detect possible mediating factors, including the frailty index, KDM biological age acceleration, allostatic load and 30 circulating biomarkers. RESULTS: Adolescent depression was associated with elevated risk for 12 chronic diseases, with strongest associations for hypothyroidism (HR = 1.29 [1.18-1.42]), diabetes (HR = 1.25 [1.13-1.38]) and chronic obstructive pulmonary disease (COPD) (HR = 1.74 [1.50-2.01]). Risks were notably higher among females and younger adults. MR confirmed likely causal relationships for hypothyroidism (OR = 1.45 [1.03-2.05]), diabetes (OR = 1.01 [1.01-1.02]) and COPD (OR = 1.04 [1.02-1.06]). Genomic colocalization revealed a shared genetic signal at the CDSN/PSORS1C1 locus between adolescent depression and hypothyroidism. Mediation analyses revealed disease-specific pathways: creatinine for hypothyroidism, testosterone for diabetes, KDM biological ageing for COPD and frailty index across all three conditions. CONCLUSIONS: Adolescent depression confers systemic vulnerability through genetic and metabolic mechanisms, with amplified risks in females and individuals aged ≤55 years. These findings support early, integrated interventions to mitigate long-term multimorbidity.

Humans↗

PubMed2026

Genetic regulation of CPEB3-mediated alternative polyadenylation associated with survival of patients with hepatocellular carcinoma.

BACKGROUND: Alternative polyadenylation (APA) is a key post-transcriptional mechanism that regulates gene expression by modulating 3'UTR length, its dysregulation has been implicated in carcinogenesis. How genetic variants influence APA to affect hepatocellular carcinoma (HCC) prognosis remains unclear. METHODS: Prognosis-APA quantitative trait loci (apaQTL) were performed using genotype and APA profiling from TCGA data. A two-stage survival analysis in 848 Chinese and 369 TCGA LIHC patients and functional validation were used to identify prognostic apaQTL in HCC progression. RESULTS: A total of 2,025 and 817 significant APA events were identified in Chinese and TCGA cohort, respectively. Besides, 859 events were associated with poor prognosis in HCC and enriched in RNA splicing / metabolism pathways. We detected 32,034 significant apaQTLs, predominantly enriched in 3'UTRs and RBP-binding regions. CPEB3 was prioritized as a key APA regulator RBP; its low expression correlated with poor patient survival and promoted proliferation, migration, and invasion in HCC cells. Notably, a functional apaQTL variant rs2037547, located in GSK3B and mediated by CPEB3, demonstrated a poor survival of HCC patients in both cohort (pooled HR=1.29, p=0.016). Mechanistically, rs2037547 promoted aberrant APA at proximal poly(A) sites of GSK3B through CPEB3, leading to increased expression of short 3'UTR isoform. This regulatory alteration enhanced HCC cell proliferation, invasion, and migration, and contributed to HCC progression. CONCLUSION: These findings elucidated the distinct role of apaQTL-mediated APA dysregulation in HCC prognosis, providing insights for prognostic stratification and potential targets for personalized therapy in HCC.

RNA-binding proteins↗

PubMed2026

Pan-cancer characterization of HMGA1 reveals its oncogenic role in tumor microenvironment and stemness: functional validation in pancreatic cancer migration and invasion.

BACKGROUND: HMGA1 is a chromatin-associated oncogenic factor implicated in tumor progression, epithelial-mesenchymal transition (EMT), stemness, and metastasis. However, its pan-cancer expression and prognostic patterns, epigenetic activation, and relationship with malignant-cell stemness/plasticity and tumor microenvironment (TME) remodeling in pancreatic adenocarcinoma (PAAD) remain incompletely defined. This study aimed to systematically characterize HMGA1 across cancers and clarify its clinical and biological relevance in PAAD. METHODS: Pan-cancer transcriptomic, clinical, genetic, methylation, immune, and stemness data were integrated from multiple public databases. PAAD single-cell RNA sequencing data were analyzed to localize HMGA1 expression, infer malignant-cell pseudotime, calculate a stemness module score, and assess ligand-receptor communication using CellChat. Public HMGA1-knockdown RNA sequencing data were reanalyzed to evaluate transcriptional remodeling. The Cancer Genome Atlas (TCGA)-PAAD expression and methylation data were used to assess TME-remodeling, immune-suppression, stemness/plasticity, cytokine/chemokine, checkpoint, and promoter-methylation features. HMGA1 expression and function were further examined using immunohistochemistry (IHC), quantitative real-time polymerase chain reaction, Western blotting, wound-healing assays, and Transwell migration and invasion assays. RESULTS: HMGA1 was upregulated in most tumor types, and high expression was associated with unfavorable survival in multiple cancers, including PAAD. In PAAD, HMGA1 was enriched in malignant epithelial cells and positively correlated with pseudotime (Spearman's rho =0.594), while the stemness module score increased along pseudotime (rho =0.748). HMGA1-high malignant cells showed markedly stronger CellChat-inferred outgoing communication, predominantly involving extracellular matrix (ECM)-receptor, adhesion-related, and selected immunomodulatory ligand-receptor axes. HMGA1 knockdown was associated with broad remodeling of EMT, TGF-β, Hedgehog, IL6/JAK/STAT3, KRAS, and cancer stem cell/stemness-related programs rather than uniform suppression of these programs. HMGA1 promoter methylation was inversely correlated with HMGA1 expression (rho =-0.633) and the TME-remodeling score (rho =-0.347). HMGA1 was associated with selected mediators, including PPIA, PLAU, ANXA1, LGALS9, TGFB1, CD276, and CD47, but not with a generalized checkpoint-high phenotype. Functionally, HMGA1 knockdown significantly reduced pancreatic cancer (PC) cell migration and invasion. CONCLUSIONS: These findings support an association-based model in which promoter hypomethylation-associated HMGA1 activation is linked to malignant epithelial stemness/plasticity, ECM/adhesion-dominant TME remodeling, selected immunomodulatory programs, and aggressive PAAD phenotypes. Further mechanistic and clinical validation is required before HMGA1 can be used for therapeutic stratification or immunotherapy-response prediction.

HMGA1↗

PubMed2026

Genetic heterogeneity affects the risk of incident depression, comorbidity, and response to environment: A prospective trajectory study.

BACKGROUND: Depression exhibits significant heterogeneity in its genetic underpinnings. The role of genetic components in the development of depression and its comorbidities remains insufficiently explored. METHODS: First, depression risk loci from a large-scale genome-wide meta-analysis were annotated to Gene Ontology (GO) terms by functional enrichment. GO-based polygenic risk scores (GO-PRS) were then calculated for individuals in the UK Biobank. Principal component analysis (PCA) was applied for dimensionality reduction, followed by cluster analysis to identify genetic subtypes of depression. Multistate models were applied to assess the impact of genetic patterns on the trajectory from healthy status to incident depression, and depression to 26 subsequent diseases, as well as the associations between environmental factors and disease trajectories across genetic subtypes. RESULTS: Participants were categorized into three genetic subtypes: immune-dominant, neuro-dominant, and comprehensive-risk. Significant differences in risk of depression and subsequent diseases, and susceptibility to environmental factors were observed across subtypes. Comprehensive-risk subtype showed higher risks of depression compared to immune-dominant (HR: 1.10, 95% CI: 1.05-1.15) and neuro-dominant subtype (HR: 1.12, 95% CI: 1.08-1.16). Comprehensive-risk subtype exhibited higher risks of transition from depression to subsequent diseases, such as anemia compared to immune-dominant subtype, and diseases of the digestive system compared to neuro-dominant subtype. Environmental factors were more strongly associated with the transition from depression to subsequent diseases in immune-dominant and comprehensive-risk subtypes, including cardiovascular, respiratory, and metabolic diseases. CONCLUSIONS: Our findings highlight the genetic heterogeneity of depression and comorbidities, and shed light on how genetic components modulate responses to environmental factors.

Humans↗

PubMed2026

Three-dimensional source apportionment and quantitative characterization of horizontal and vertical transport fluxes of O3 and its precursors in the Beijing-Tianjin-Hebei region, China.

Persistent surface ozone (O3) pollution in the Beijing-Tianjin-Hebei (BTH) region is driven by coupled precursor emissions and multi-scale transport, yet its altitude-dependent transport and source contributions remain insufficiently quantified. Here we integrated the Weather Research and Forecasting and the Comprehensive Air Quality Model with Extensions with the Ozone Source Apportionment Technology and a quantitative transport-flux framework to characterize three-dimensional source apportionment and horizontal/vertical fluxes of O3, Volatile Organic Compounds‌ (VOCs), and Nitrogen Oxides (NOx) across dynamic meteorological scenarios. Simulations showed that VOCs and NOx were dominated by local emissions near the surface (73.61 %-82.18 %), whereas surface O3 was primarily controlled by regional transport, with local contributions of only 11.01 %-13.75 %. Notably, the transport dominance further strengthened with altitude, exceeding 93 % at 1.8 km. Industrial and transportation emissions together contributed more than 75 % of precursor emissions and account for approximately 80 % of O3 formation, while favorable/unfavorable meteorological years modulated long-range transport efficiency and the vertical distribution of contributions. Horizontal flux analysis highlighted three major pathways (Northwest-Southeast, Southeast-Northwest, and Southwest-Northeast), with Shijiazhuang serving as a critical pollutant "sink" across altitude layers. Vertical fluxes revealed an altitude transition near 600 m: net downward transport dominated below 600 m, whereas enhanced summer convection promoted upward transport above 600 m. These results support altitude-dependent, scenario-specific strategies for coordinated regional O3 mitigation in the BTH region.

Ozone↗

PubMed2026

A single-nucleus and spatial transcriptomic atlas of poplar leaves reveals the regulation of leaf polarity and cuticle deposition.

Leaf adaxial-abaxial polarity is fundamental for plant morphogenesis and environmental adaptation through asymmetric cell differentiation. Emerging evidence reveals dorsoventral metabolic gradients act downstream of transcriptional networks to fine-tune cellular specialization. While conserved transcription factors (e.g., HD-ZIP III and KANADI) establish initial polarity, the molecular networks driving position-specific cellular differentiation and their integration with metabolic adaptation remain unclear. Leveraging single-nucleus and spatial transcriptomics, we resolve major cell classes (mesophyll, epidermal, and vascular-associated) and their adaxial-abaxial subtypes, revealing dorsoventral polarity in transcriptional profiles and metabolic pathways. Adaxial cells are enriched in phenylpropanoid/flavonoid biosynthesis, while abaxial cells show preferential activation of stress and hormone signaling. Notably, we identify MYC2 as a key regulator of adaxial cuticle biosynthesis, binding to promoters of lipid biosynthetic and transport genes (e.g., CER10 and LTPG1) and promoting cuticle thickening. Our study uncovers how positional identity shapes transcriptional and metabolic polarity in leaves, with MYC2 emerging as a central regulator coordinating organ-specific adaptations. These findings provide insights into the spatial regulation of plant development and stress resilience, offering potential strategies for engineering stress-tolerant woody crops.

Plant Leaves↗

PubMed2025

Polygenic enrichment analysis in multi-omics levels identifies cell/tissue specific associations with schizophrenia based on single-cell RNA sequencing data.

OBJECTIVE: Understanding the specific cellular origin and tissue heterogeneity in schizophrenia is critically important for exploring the disease etiology. This study aims to investigate these aspects by performing multiple analyses based on omics data. METHOD: We performed single-cell disease relevance score (scDRS) algorithm to link brain single-cell RNA sequencing (scRNA-seq) with schizophrenia risk across multi-omics scales at single-cell resolution. This approach identified cell types with overexpression of schizophrenia-related genes implicated by multi-omics panels (ATAC-seq, RNA-seq, TWAS, and GWAS). Schizophrenia-related genes from these multi-omics panels were extracted and combined with scRNA-seq data to calculate scDRS. Subsequently, the cell-type vs. disease association and tissue heterogeneity were assessed using scDRS for each omics panel. RESULTS: We identified two novel cell subpopulations in the brain that differentially express SCUBE3 (59 cells, 7.0 %) and FN1 (21 cells, 2.5 %). At the individual cell level, schizophrenia-associated cell subpopulations included microglial cell associated with ATAC-seq panel (Passociation = 0.002, Pheterogeneity = 0.009) and deep layer neuron suggestively associated with GWAS panel (Passociation = 0.033, Pheterogeneity = 0.017). At the brain tissue level, microglial cell was significantly associated with cortical plate in ATAC-seq panel (Passociation = 0.002, Pheterogeneity = 0.011). Gene level analysis identified several genes associated with schizophrenia across multi-omics panels. CONCLUSIONS: Our study outlines the signature of cell subpopulations, brain regions, and disease risk genes in schizophrenia at single-cell resolution across multi-omics scales. These findings provide a reference for future precision medicine approaches targeting specific cell types and brain regions in schizophrenia.

Schizophrenia↗

PubMed2025

Glutamate molecular structure and protein affect the inhibition of breast cancer cell metastasis: Cell-derived exosomes inhibitory effects through the MAPK signaling pathway.

The aim of this study was to investigate the inhibitory effect of glutamate molecular structure and protein on breast cancer cell metastasis and the potential inhibitory mechanism of cell-derived exosomes via MAPK signaling pathway. Breast cancer cell lines with high metastatic potential were selected by in vitro cell culture technique. The effects of specific inhibitors of glutamic acid on the proliferation and metastasis of breast cancer cells were studied. Changes in protein expression profiles were analyzed by proteomics techniques to identify key proteins associated with breast cancer metastasis. Breast cancer cells were treated with inhibitors of the MAPK signaling pathway to evaluate their effect on cell metastasis and compare with exosome treatment. The results showed that the specific inhibitors of glutamate molecular structure could significantly inhibit the proliferation and metastasis of breast cancer cells. Proteomic analysis revealed several down-regulated proteins that are closely related to breast cancer metastasis.

Humans↗

PubMed2025

Perinatal dysfunction of innate immunity in cystic fibrosis.

In patients with cystic fibrosis (CF), repeated cycles of infection and inflammation eventually lead to fatal lung damage. Although diminished mucus clearance can be restored by highly effective CFTR modulator therapy, inflammation and infection often persist. To elucidate the role of the innate immune system in CF etiology, we investigated a CF pig model and compared these results with those for preschool children with CF. In newborn CF pigs, we observed changes in lung immune cell composition before the onset of infection that were dominated by increased monocyte infiltration, whereas neutrophil numbers remained constant. Flow cytometric and transcriptomic profiling revealed that the infiltrating myeloid cells displayed a more immature status. Cells with comparably immature transcriptomic profiles were enriched in the blood of CF pigs at birth as well as in preschool children with CF. This pattern coincided with decreased CD16 expression in the myeloid cells of both pigs and humans, which translated into lower phagocytic activity and reduced production of reactive oxygen species in both species. These results were indicative of a congenital, translationally conserved, and functionally relevant aberration of the immune system in CF. In newborn wild-type pigs, CFTR transcription in immune cells, including lung-derived and circulating monocytes, isolated from the bone marrow, thymus, spleen, and blood was below the detection limits of highly sensitive assays, suggesting an indirect etiology of the observed effects. Our findings highlight the need for additional immunological treatments to target innate immune deficits in patients with CF.

Cystic Fibrosis↗

PubMed2025

A SlEIN2-centered epigenetic network equilibrates fruit ripening and innate immunity in tomato.

Ethylene and DNA/RNA methylation serve as essential factors in controlling fruit ripening. In tomato, the mRNA N6-methyladenosine (m6A) demethylase SlALKBH2 regulates mRNA stability of the DNA 5-methylcytosine demethylase gene SlDML2 via modulating m6A modifications. However, the interplay between ethylene and these epigenetic marks remains unclear. Here, we show that SlDML2 expression is significantly inhibited in slein2 fruits, but remains unchanged in the high-order sleil mutant (sleil1 sleil2 sleil3/SlEIL3 sleil4 and sleil1 sleil2/SlEIL2 sleil3 sleil4) fruits, indicative of post-transcriptional regulation of SlDML2 expression by SlEIN2, a core ethylene signaling component acting upstream of the master transcription factors SlEILs. Interestingly, SlEIN2 preferentially regulates the asymmetric CHH methylation in promoters of several key ripening regulator genes. Mechanistically, SlEIN2 physically interacts with SlALKBH2, which promotes SlDML2 expression in a SlEIN2-dependent manner. Furthermore, SlAGO4A and SlAGO4B, components of the RNA-directed DNA methylation pathway, were upregulated in slein2 fruits. Silencing SlAGO4A/B in wild-type fruit caused precocious ripening with necrosis, indicative of hyperimmunity. Conversely, SlAGO4A/B silencing in slein2 markedly delayed this hyperimmunity. Taken together, our study reveals that ethylene, beyond transcriptional regulation, employs an elaborate epigenetic machinery mediated by the SlAGO4A/B-SlEIN2-SlALKBH2 module to balance fruit ripening and innate immunity.

Solanum lycopersicum↗

PubMed2016

Matrix stiffness promotes cartilage endplate chondrocyte calcification in disc degeneration via miR-20a targeting ANKH expression.

The mechanical environment is crucial for intervertebral disc degeneration (IDD). However, the mechanisms underlying the regulation of cartilage endplate (CEP) calcification by altered matrix stiffness remain unclear. In this study, we found that matrix stiffness of CEP was positively correlated with the degree of IDD, and stiff matrix, which mimicked the severe degeneration of CEP, promoted inorganic phosphate-induced calcification in CEP chondrocytes. Co-expression analysis of the miRNA and mRNA profiles showed that increasing stiffness resulted in up-regulation of miR-20a and down-regulation of decreased ankylosis protein homolog (ANKH) during inorganic phosphate-induced calcification in CEP chondrocytes. Through a dual luciferase reporter assay, we confirmed that miR-20a directly targets 3'-untranslated regions of ANKH. The inhibition of miR-20a attenuated the calcium deposition and calcification-related gene expression, whereas the overexpression of miR-20a enhanced calcification in CEP chondrocytes on stiff matrix. The rescue of ANKH expression restored the decreased pyrophosphate efflux and inhibited calcification. In clinical samples, the levels of ANKH expression were inversely associated with the degeneration degree of CEP. Thus, our findings demonstrate that the miR-20a/ANKH axis mediates the stiff matrix- promoted CEP calcification, suggesting that miR-20a and ANKH are potential targets in restraining the progression of IDD.

3' Untranslated Regions↗