Biomedical subjects

Chao Li

Publications and source records attributed to Chao Li.

11 recordsLinked to original sources

PubMed2026

Genome-wide barriers to gene flow reveal the genetic basis of viviparity evolution in a lizard.

Viviparity (live-bearing) is a major evolutionary transition repeatedly linked with ecological and evolutionary diversification throughout vertebrates. Live-bearing reproduction entails a novel suite of phenotypes and life history traits, but the genetic processes by which such a reproductive innovation evolves are unknown. Remarkable among amniotes, the common lizard (Zootoca vivipara) has extant oviparous (egg-laying) and viviparous lineages and a to-date unresolved history of parity mode emergence. This species represents an ideal model to reconstruct the evolutionary and genetic mechanisms of how viviparity arises. By analyzing whole genomes of individuals from across the species' distribution, we robustly show that viviparity evolved once. However, gene flow from oviparous to viviparous populations is found to be long-term and extensive, causing pronounced gene tree discordance. We inferred signals of selection for viviparity in many independent regions across the genome, and these were recruited over considerable time. Genomic barriers to gene flow between oviparity and viviparity were found genome wide. These are enriched for regions under selection for parity mode and for genes known to be involved in pregnancy and parturition in squamates and mammals. Further implicating their functional role in viviparity, we show that genes in genomic regions under selection and resisting gene flow are more highly expressed in the uterus of viviparous lizards during pregnancy. Our study demonstrates that viviparity in an amniote evolved by selection in the face of gene flow and primarily by the genome-wide accumulation of functional regulatory variants. These results reveal how complex adaptive innovations can arise and be maintained.

Animals↗

PubMed2026

TNF-NF-κB signaling mediates immune-biomineralization crosstalk during shell repair under ocean acidification in Mytilus edulis.

Ocean acidification (OA) impairs biomineralization in bivalves, but its effects on immune-biomineralization crosstalk during shell repair remain unknown. Here, we exposed adult Mytilus edulis bearing standardized shell perforations to three pH levels (8.1, 7.9, and 7.7) for up to 40 days. OA slowed early repair and caused microstructural disorganization and an approximately 87% reduction of compressive strength at pH 7.7, yet the damaged area appeared largely closed by day 15, suggesting a decoupling between morphological closure and functional recovery. In addition, transcriptomic profiling of hemocytes and mantle tissue, based on an average of 6.5 Gb of clean reads per sample mapped to the M. edulis reference genome (NCBI Assembly GCF_000511035.1), revealed that these shell-level defects were accompanied by coordinated immune and metabolic reprogramming. Hemocytes, the primary immune effector cells of bivalves, exhibited pH- and time-dependent shifts with moderate acidification (pH 7.9) promoting inflammatory transcripts, whereas severe acidification (pH 7.7) suppressed these signals while upregulating stress-associated pathways; both treatments consistently downregulated lysosomal proteases and NF-κB negative regulators. The mantle, a primarily mineralizing organ, paradoxically upregulated immune-related genes while suppressing oxidative phosphorylation and extracellular matrix pathways. This tissue-level imbalance, with hemocytes recruited but functionally constrained and mantle metabolically suppressed yet immunologically activated, points to TNF-NF-κB pathway modulation as a key mediator of shell repair under acidification. Our findings demonstrate that visible shell closure masks underlying structural and mechanical failure, and that immune regulation, rather than simple suppression or activation, critically shapes the repair outcome. These results advocate for multifunctional indicators beyond closure area to assess shell integrity in acidified marine environments.

Animals↗

PubMed2026

Molecular characterization and biological characteristics of a highly pathogenic recombinant ALV-J strain (HUE2023) with cross-clade gp85 recombination.

Avian leukosis virus subgroup J (ALV-J) has undergone extensive diversification into phylogenetically distinct clades, yet whether recombination between these clades within the gp85 envelope glycoprotein generates variants with altered pathogenicity has received little direct investigation. A field strain (HUE2023) was recovered from breeding roosters displaying vascular tumors. The viral genome was sequenced and subjected to phylogenetic and recombination analyses. The three-dimensional structure of gp85 was predicted with AlphaFold3; electrostatic surface potentials and surface hydrophobicity were computed using the Adaptive Poisson-Boltzmann Solver and the Eisenberg hydrophobicity scale, respectively. Pathogenicity and immunosuppressive effects were assessed in Hy-Line Brown chickens. Recombination analysis revealed that HUE2023 is an inter-clade recombinant derived from Clade 1.1 (major parent: JS14NT01) and Clade 1.2 (minor parent: JS09GY3). A single-residue deletion at position 61 within receptor-binding domain 1 (RBD-1), unique to the recombinant, induced a localized conformational rearrangement that generated a concentrated electronegative surface patch and a contiguous hydrophobic pocket not observed in either parental gp85. Animal challenge showed that HUE2023 is highly pathogenic: female chickens in the high-dose group reached only 61% survival and displayed significant growth retardation (P < 0.05) together with marked immunosuppression. The recombination in the RBD-1 led to local conformational rearrangement, resulting in a concentrated and negatively charged surface area as well as a continuous hydrophobic pocket, which were never present in any of the parental gp85 sequences. These results indicate that gp85 recombination across clades can yield variants with fundamentally altered receptor-binding surfaces and argue for integrating structural surveillance into ALV-J monitoring programmes.

Animals↗

PubMed2026

High-frequency contralesional dorsal premotor cortex and low-frequency contralesional primary motor cortex rTMS in subacute stroke with severe upper limb impairment: comparable motor outcomes and differential regional degree centrality changes.

BACKGROUND: The contralesional dorsal premotor cortex has been proposed as a potential neuromodulatory target for patients with severe upper limb impairment due to subacute ischemic stroke. This proof-of-concept study aimed to compare behavioral outcomes and resting-state neuroimaging findings between high-frequency repetitive transcranial magnetic stimulation (rTMS) over the contralesional dorsal premotor cortex and guideline-supported low-frequency stimulation over the contralesional primary motor cortex. METHODS: In this randomized trial, 46 patients with severe upper limb impairment in the subacute stage after ischemic stroke were randomly assigned to receive either high-frequency rTMS over the contralesional dorsal premotor cortex or low-frequency rTMS over the contralesional primary motor cortex. Low-frequency stimulation over the contralesional primary motor cortex served as an evidence-supported active comparator for poststroke upper limb motor recovery. Stimulation was administered five times per week for two weeks using magnetic resonance imaging-guided neuronavigation. All participants received concurrent standard rehabilitation therapy. The primary outcome was the Fugl-Meyer Assessment for Upper Extremity. Secondary outcomes included the Arm Subscore of the Motricity Index, the Hong Kong version of the Functional Test for the Hemiplegic Upper Extremity, the Modified Barthel Index, and resting-state functional magnetic resonance imaging-derived degree centrality. RESULTS: Both groups showed significant improvements in the primary and secondary behavioral measures (p < 0.01), with no significant between-group differences in the magnitude of change (p > 0.05). In neuroimaging analyses, patients receiving high-frequency rTMS over the contralesional dorsal premotor cortex showed significantly greater degree centrality changes in the ipsilesional middle occipital gyrus, contralesional medial superior frontal gyrus, and contralesional middle frontal gyrus than those receiving low-frequency rTMS over the contralesional primary motor cortex (p < 0.05). Within the high-frequency stimulation group, degree centrality changes in the ipsilesional middle occipital gyrus were positively correlated with improvements in the Fugl-Meyer Assessment for Upper Extremity (r = 0.619, false discovery rate-corrected p = 0.018). CONCLUSIONS: High-frequency rTMS over the contralesional dorsal premotor cortex produced behavioral improvements comparable to guideline-supported low-frequency rTMS over the contralesional primary motor cortex, without establishing superiority or formal non-inferiority. Exploratory neuroimaging analyses showed greater degree centrality changes in the ipsilesional middle occipital gyrus after high-frequency premotor stimulation, and these changes correlated with upper-limb motor improvement. These findings support further investigation of contralesional dorsal premotor cortex-targeted high-frequency rTMS for severe subacute post-stroke upper limb impairment. REGISTRATION: URL: http://www.chictr.org.cn; Unique identifier: ChiCTR2000038049.

Humans↗

PubMed2026

Full-length single-cell spatial transcriptomics reveals spatial and cell-type-specific transcript isoforms in the primate brain.

The primate brain exhibits complex RNA alternative splicing heterogeneity crucial for functional complexity, yet systematic spatial isoform characterization has been lacking. We developed Fullscope-seq, a full-length single-molecule large field-of-view spatial transcriptomics sequencing method at single-cell resolution, based on programmed concatenation cDNA for multiple long-read sequencing platforms. Applying Fullscope-seq to the macaque brain, we uncovered thousands of genes exhibiting differential transcript usage (DTU) across cortical layers, cell types and brain regions. Fullscope-seq resolved hundreds of major isoform switches across distinct brain regions and identified DTUs between superficial and deep cortical layers. Cortical layer-specific DTUs showed cell-composition dependence, whereas regional DTUs were regulated according to both cellular composition and spatial contexts. These isoform variations showed substantial enrichment for neuropsychiatric disorder-associated genes and were conserved across platforms and species. Our study establishes a scalable framework for spatial isoform analysis and provides a resource for understanding transcriptomic diversity in complex tissues.

Animals↗

PubMed2026

Construction of circRNA-miRNA-mRNA regulatory networks in the intestine of turbot (Scophthalmus maximus) following Vibrio anguillarum infection.

Circular RNAs (circRNAs) play pivotal roles in post-transcriptional regulation by acting as molecular sponges for microRNAs (miRNAs) within the competitive endogenous RNA (ceRNA) network. However, the regulatory mechanisms in teleost immune responses remain poorly understood. In this study, circRNA-miRNA-mRNA networks were investigated in turbot (Scophthalmus maximus) following Vibrio anguillarum infection to elucidate host-pathogen interactions. Through high-throughput sequencing of intestinal tissues, a total of 50 differentially expressed circRNAs (DE-circRNAs) (18 at 2 hpi, 16 at 12 hpi, 16 at 48 hpi), 212 DE-miRNAs (11 at 2 hpi, 70 at 12 hpi, 15 at 48 hpi), and 1774 DE-mRNAs were identified. Functional enrichment analyses (GO/KEGG) revealed significant associations with immune pathways, including the MAPK signaling pathway and gap junction. An integrated circRNA-miRNA-mRNA regulatory network was constructed, highlighting key interactions including novel_circ_0002573/DE-miR-27a-3p/FGB and novel_circ_0002423/novel_347/GNE, which may regulate inflammatory and antibacterial responses. The expression patterns of selected circRNAs, miRNAs and mRNAs were validated using qRT-PCR, confirming the reliability of the sequencing results. Importantly, fibrinogen beta chain (FGB) and CXCR4/CXCL12 signaling were identified as critical immune modulators. These findings provide insights of the ceRNA regulatory networks involved in teleost intestinal immunity and provide potential molecular targets for selective breeding of disease resistance in this species.

Animals↗

PubMed2026

Efficient CRISPR/Cas-SF01 genome editing tools with high editing efficiency in allotetraploid oilseed rape.

CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-Cas9 has been widely utilized for plant genome editing, but the protospacer adjacent motif (PAM) requirement limits its editing scope. CRISPR/Cas12i3 belongs to the type-VI Cas system that has gained extensive attention due to its smaller size and less restricted canonical TTN PAM sequence. In this study, we explored the newly developed Cas-SF01 system (Cas12i3 variant) for genome editing in oilseed rape. We established an efficient protoplast transformation system in oilseed rape to compare editing efficiency between Cas-SF01 and Cas9. Cas-SF01 shows cleavage activities at the tested 5'-TTN-3' PAM sites with editing outcomes sharing considerable similarities with the CRISPR-Cas9 system in protoplast. Cas-SF01 also induces high efficiency mutagenesis for multiple target sites in stable transformed oilseed rape lines, generating mutants with multilocular silique and male sterile phenotypes. Furthermore, Cas-SF01-derived cytosine base editors (CBEs) were developed to produce targeted C-to-T base edits. Compared to SpCas9, Cas-SF01 has an expanded PAM range and effectively recognizes TTN PAMs, which has substantially broadened the scope of editable sites within the rapeseed genome. No mutations were identified at the putative off-target sites among the edited plants. This study developed a robust, first-of-its-kind Cas12 system in the allotetraploid Brassica napus, expanding the scope of editing and enriching genome-editing toolkits for biological research and genetic improvement.

Brassica napus↗

PubMed2026

Phosphoproteomics identification of ERK-dependent activation of Rps6kb1 in cardiac hypertrophy.

Cardiomyocyte growth is tightly controlled by multiple signaling pathways. Identification of master kinases in this process is essential in exploring potential targets for the treatment of pathological cardiac hypertrophy and heart failure. Here we identified the mTOR-independent activation of ribosomal protein S6 kinase b1 (Rps6kb1) during cardiomyocyte growth. By utilizing phosphoproteomics in primary neonatal rat ventricular myocytes, we revealed Rps6kb1 as one of most activated kinases under growth stimulation. We further demonstrated the role of Rps6kb1 phosphorylation in pathological cardiac hypertrophy and heart failure. We showed that the phosphorylation of multiple sites in Rps6kb1, including T367 in the kinase domain and S418/T421/S424 in the C-terminal domain, is not directly regulated by the activity of mTOR but coupled with the activation of the MEK1/ERK axis. In mice, cardiomyocyte-specific deletion of Rps6kb1 significantly inhibited both constitutively active ERK- and pressure overload-induced cardiac hypertrophy. In contrast, cardiomyocyte-specific overexpression of wild-type Rps6kb1, rather than the phosphorylation-defective mutant, elevated cardiac hypertrophy and augmented pressure overload-induced heart failure. In conclusion, our findings reveal that the MEK/ERK axis primes Rps6kb1 activation through phosphorylation of 2 separate domains of Rps6kb1, which may play an essential role in cardiac hypertrophy and heart failure under hemodynamic stress.

Animals↗

PubMed2025

Engineer the eukaryotic OMEGA-Fanzor systems for genome editing in plants.

The activity of the eukaryotic OMEGA-Fanzor genome editing system remains limited in plants. We engineered the Fanzor nucleases SpuFz1, GtFz1, NlovFz2, and MmeFz2 in plants, with NlovFz2 being the most efficient, achieving up to 50.0% editing in regenerated rice plants, making it a promising tool for plant genome editing.

Oryza↗

PubMed2025

Multi-ancestry genome-wide and transcriptome-wide association analyses identified new risk loci and genes for inflammatory bowel disease.

To advance genetic understanding of inflammatory bowel disease (IBD), we conducted genome-wide association meta-analyses of 63,415 IBD cases of European and East Asian descendants and identified 90 previously unknown risk loci. Integrating multi-ancestry transcriptome-wide association studies (TWAS), cell type-specific TWAS, alternative splicing (AS-WAS), and alternative polyadenylation (APA-WAS) analyses using RNA-seq data from normal colon tissues of 707 European and 364 East Asian individuals, we uncovered 506 high-confidence IBD risk genes, including 384 not previously reported. These genes converge on immune regulation, microbial interaction, and other pathways central to IBD pathogenesis, with over half showing transcriptional dysregulation supported by single-cell and spatial omics analyses. Notably, 46 risk genes are targeted by 225 drugs that have been approved or in Phase II/III trials, including sulfasalazine already used in IBD therapy. Our study findings deepen the understanding of IBD genetics and support the development of precision medicine for its prevention and treatment.

GWAS↗

PubMed2024

Is there a causal relationship between resistin levels and bone mineral density, fracture occurrence? A mendelian randomization study.

BACKGROUND: In a great many of observational studies, whether there is a relevance of resistin levels on bone mineral density (BMD) and fracture occurrence has been inconsistently reported, and the causality is unclear. METHODS: We aim to assess the resistin levels on BMD and fracture occurrence within a Mendelian randomization (MR) analysis. Exposure and outcome data were derived from the Integrative Epidemiology Unit (IEU) Open genome wide association studies (GWAS) database. Screening of instrumental variables (IVs) was performed subject to conditions of relevance, exclusivity, and independence. Inverse variance weighting (IVW) was our primary method for MR analysis based on harmonized data. Weighted median and MR-Egger were chosen to evaluate the robustness of the results of IVW. Simultaneously, heterogeneity and horizontal pleiotropy were also assessed and the direction of potential causality was detected by MR Steiger. Multivariable MR (MVMR) analysis was used to identify whether confounding factors affected the reliability of the results. RESULTS: After Bonferroni correction, the results showed a suggestively positive causality between resistin levels and total body BMD (TB-BMD) in European populations over the age of 60 [β(95%CI): 0.093(0.021, 0.165), P = 0.011]. The weighted median [β(95%CI): 0.111(0.067, 0.213), P = 0.035] and MR-Egger [β(95%CI): 0.162(0.025, 0.2983), P = 0.040] results demonstrate the robustness of the IVW results. No presence of pleiotropy or heterogeneity was detected between them. MR Steiger supports the causal inference result and MVMR suggests its direct effect. CONCLUSIONS: In European population older than 60 years, genetically predicted higher levels of resistin were associated with higher TB-BMD. A significant causality between resistin levels on BMD at different sites, fracture in certain parts of the body, and BMD in four different age groups between 0-60 years of age was not found in our study.

Bone Density↗