Research

microscope

1. Mechanisms of mRNA polyadenylation

mRNA polyadenylation is a crucial step in eukaryotic gene expression, influencing mRNA stability, translatability, and cytoplasmic translocation. This process also impacts pre-mRNA splicing and transcription termination. Recent research has linked mRNA polyadenylation to various cellular processes, including responses to environmental stimuli, pathogen invasion, and cancer development.

Our lab aims to uncover the molecular mechanisms of mRNA polyadenylation at genetic, transcriptomic, and genomic levels, examining both tissue-specific and single-cell contexts. We use functional genomic tools to study these mechanisms in model plants such as Arabidopsis and rice, focusing on the biochemical, genetic, and molecular functions of polyadenylation factors and their roles in plant gene expression during development and stress responses.

Leveraging advanced RNA-sequencing technologies, we conduct large-scale analyses of alternative polyadenylation (APA) to map the comprehensive landscape of APA across single cells and entire organisms. With over 70% of genes utilizing APA for expression regulation, the role of APA in cellular and organismal functions is just beginning to be understood. Our multi-omics approach reveals that APA is linked to various cellular processes, including circadian regulation, epigenetic control, and quantitative traits, highlighting its role as a central hub for fine-tuning gene expression at the transcriptome level.

Test tubes

2. Gene expression in neuronal cell communication

Our research collaboration extends to gene expression regulation in human induced pluripotent stem cells (iPSCs), particularly in neuronal disease models like Alzheimer’s disease (AD). AD and related dementias present a global challenge, characterized by extensive neuronal loss, yet the mechanisms behind this loss remain unclear.

Microglia, the central nervous system’s resident immune cells, are suspected to play a role in AD pathogenesis. We explore how extracellular vesicles (EVs) — lipid bilayer membrane vesicles 30 to 200 nm in diameter, generated by the endosomal system — influence cell-to-cell communication. EVs contain proteins, mRNA, and miRNA that can modulate gene expression in recipient cells. Notably, microglia-derived microvesicles have been shown to regulate neuronal excitability, representing a promising area of research in neurodegenerative diseases.

We employ multi-omics tools to analyze the components of microglia-derived EVs in AD-like models derived from patient and control samples. Identifying these potential markers and their roles in regulating gene expression, especially concerning neuronal degeneration, could provide new insights into neuronal gene expression regulation and identify potential targets for future drug development.

Lab bench

3. Gene expression regulation in development and responses

In collaboration with various research groups, we study gene function and expression related to plant development and responses to environmental stimuli.

Mangroves, for example, are woody plants that thrive in extreme coastal conditions. Their unique adaptations, such as viviparous seed reproduction, pneumatophores for submergence and hypoxia, and various salinity responses, are of significant interest due to their potential applications in crop breeding programs. We investigate these adaptations through gene expression analysis to understand the molecular mechanisms underlying these features.

Additionally, we explore transgenerational transmission of parental phenotypes and epigenetic modifications, which may offer new insights into how organisms adapt to acute stresses and how such “memories” are passed to future generations. Our past research also includes efforts in plant disease resistance through genetic engineering of antimicrobial peptides to improve crop resilience.

Publications

  • Yang D, Li X, He H, He N, Li QQ, Fu H. Alternative Polyadenylation Dynamics During the Rice Blast Immune Response.. Molecular plant pathology. 2026;27(7):e70301. doi:10.1111/mpp.70301

    Alternative polyadenylation (APA) is a widespread post-transcriptional regulation that generates transcripts with variable 3' untranslated region (UTR) lengths. In this study, poly(A)-tag sequencing (PAT-seq) was performed during the process of Magnaporthe oryzae infection of rice. Genome-wide dynamic changes of APA profiles were identified during the response, with widespread shortening of gene transcripts. Shortened genes were found to function in biological processes associated with stress responses. Importantly, the changes of 3' UTR length were negatively related to the expression levels of the corresponding genes. In addition, we found that the expression levels of miRNAs that bind to the 3' UTRs of the APA genes were negatively correlated with the expression levels of the corresponding APA genes. This correlation suggests a potential regulatory strategy where immune-related genes might evade miRNA-mediated repression via 3' UTR shortening. To assess the biological impact of these APA dynamics, we functionally characterised a representative candidate, Os05g0509500, that exhibits complex and dynamic APA site switching upon infection. Knockout of Os05g0509500 via CRISPR/Cas9 revealed its positive regulatory role in rice blast, offering initial insights into the function of APA in this disease. Several core polyadenylation protein factors were significantly differentially expressed during the rice response. Our results revealed that the precise transcriptome-wide poly(A) site selection of genes during rice blast challenge, and the relationships between miRNAs and targeted APA genes, are tightly regulated. Such a mechanism provides a new perspective for designing novel strategies to control rice blast disease.

  • Li Y, Zhao Z, Lin A, Ye C, Qu H, Li QQ. PAPS1-associated alternative polyadenylation changes correlate with pollen development and flowering time in Arabidopsis.. Plant physiology. 2026. doi:10.1093/plphys/kiag393

    Alternative polyadenylation (APA) serves as a critical co-transcriptional regulatory mechanism that shapes mRNA fate and protein function. POLY(A) POLYMERASE 1 (PAPS1) contributes to poly(A) tail synthesis; however, its relationship to poly(A) site choice is less well characterized. Here, we profiled PAPS1-associated poly(A) site usage across Arabidopsis (Arabidopsis thaliana) tissues using UMI-quantified poly(A) tag sequencing (qPAT-seq). Across tissues, paps1-4 showed tissue-dependent shifts in poly(A) site usage, with pronounced effects in pollen. In the paps1-4 mutant, polyadenylated transcripts were biased toward longer isoforms with significantly longer 3' UTRs, and poly(A) signal usage was altered in distinct ways: a U-to-A shift immediately downstream of poly(A) sites occurred in leaves and flower buds, whereas pollen showed more complex motif combinations with reduced usage of canonical cis-elements. In pollen, DE-APAGs included multiple genes previously linked to a pollen developmental regulatory module involving AT-RICH INTERACTING DOMAIN-CONTAINING PROTEIN 1 (ARID1), RETINOBLASTOMA RELATED 1 (RBR1), and DUO POLLEN1 (DUO1), showing coordinated changes in APA and expression. Flowering-time-related terms were also enriched among genes with significantly altered APA in paps1-4, including regulators connected to FLOWERING LOCUS C (FLC) and photoperiod responses, and these changes were more consistent with poly(A) site selection differences than with widespread poly(A) tail length changes. Notably, multiple poly(A) factors in pollen switched between single- and multi-poly(A) site usage patterns in paps1-4, representing a major component of PAPS1-associated regulation. Together, these findings reveal tissue-specific APA alterations associated with PAPS1 and highlight pollen and flowering pathways as key contexts of its regulatory influence.

  • Weng Y, Chen T, Zhou L, Zhang Y-Y, Shen Y, Li QQ. Divergent Evolutionary Routes of Hypoxia Tolerance in Mangroves.. Molecular ecology. 2026;35(9):e70372. doi:10.1111/mec.70372

    Extreme environments impose conditions that are lethal to most organisms, yet a limited number of evolutionary lineages have adapted to chronic stress. However, whether distinct lineages exposed to the similar stressors rely on shared adaptive mechanisms remains poorly understood. Mangroves offer a compelling system to address this question, as they thrive in shared harsh environments. Here, we examine two mangrove species, Kandelia obovata and Avicennia marina, both inhabiting intertidal zones but deploying contrasting strategies to withstand hypoxia. Phylogenomic analyses revealed that K. obovata experienced an ancestral whole-genome duplication (WGD) prior to intertidal colonization, whereas A. marina underwent two rounds of WGD that expanded its repertoire of hypoxia-responsive genes. Comparative genomic analyses showed that K. obovata underwent contraction of hypoxia-related gene families and harbours fewer gene copies, but each copy is enriched with cis-regulatory elements. In contrast, A. marina retained significantly more hypoxia-related genes derived from recent WGD. Transcriptomic profiling under controlled hypoxia gradients further showed that K. obovata relies on rapid but stable gene expression responses, while A. marina exhibits delayed but hierarchical gene regulatory networks. Together, these results demonstrate that even under identical extreme selective pressures, plant lineages can follow distinct evolutionary routes, shaped by differences in genome evolution and regulatory architecture.

  • Zhao Z, Zhou J, Zhao D, Lu X, Li QQ. Alternative polyadenylation dynamics shape pollen development at single-cell resolution.. The Plant cell. 2026. doi:10.1093/plcell/koag128

    Alternative polyadenylation (APA) is a widespread co-transcriptional mechanism that regulates gene expression in growth, development and environmental responses. Pollen development is essential for the reproductive success of flowering plants, yet the contribution of APA to this process remains poorly understood. Here, we combine bulk RNA-seq in multiple tissues with single-nucleus transcriptomics across pollen developmental stages to systematically characterize APA dynamics during Arabidopsis thaliana pollen development. We show that mature pollen exhibits the most tissue-specific APA profile among the examined tissues, characterized by widespread 3' untranslated region (3' UTR) shortening. At single-nucleus resolution, APA patterns display pronounced temporal and cell-type specificity, particularly during the transition from bicellular to tricellular pollen and during vegetative nucleus maturation. Sperm nuclei exhibit the most distinct poly(A) site usage patterns. Moreover, genetic analyses of representative genes showed that altered poly(A) site usage is associated with changes in transcript abundance and pollen development phenotypes. Consistent with these observations, in vivo reporter assays showed that 3' UTR configurations are sufficient to modulate gene expression at the transcript level. Together, our study establishes APA as a structured co-transcriptional regulatory layer during pollen development and provides a framework for understanding 3' end-mediated gene regulation in male gametophytes.

  • Wang W, Wang M, Wang R, Lin S, Huang F, Jin Y, He N, Cheng Z, Li QQ, Yang D. Screening and Validation of Interacting Proteins of Receptor-like Cytoplasmic Kinase OsRLCK118 Involved in Rice Blast Resistance.. Journal of fungi (Basel, Switzerland). 2026;12(2). doi:10.3390/jof12020148

    Rice blast, caused by Magnaporthe oryzae, severely threatens global rice production. Although the receptor-like cytoplasmic kinase OsRLCK118 positively regulates rice immunity, its downstream signaling mechanism remains unknown. To systematically identify OsRLCK118-interacting proteins, we performed immunoprecipitation-mass spectrometry (IP-MS) and a yeast library screen, yielding 781 and 287 candidates, respectively, with 35 overlapping hits. Among these, OsSAMS1, a known positive regulator of blast resistance, was selected for validation. Membrane yeast two-hybrid, split-luciferase complementation, and co-immunoprecipitation assays confirmed the physical interaction between OsRLCK118 and OsSAMS1. Furthermore, in vitro kinase assays showed that OsRLCK118 specifically phosphorylates OsSAMS1. These results uncover a novel signaling axis connecting pathogen recognition to ethylene biosynthesis via OsRLCK118-dependent phosphorylation of OsSAMS1, providing both mechanistic insight into rice immunity and potential genetic targets for resistance breeding.

  • Short AW, Liu J, Li QQ, Ye C, Cao K, Wee AKS. Transcriptome wide evidence of interspecies differences in the regulation of photoprotection between chilling-tolerant and -sensitive mangrove species.. The Plant journal : for cell and molecular biology. 2026;125(1):e70611. doi:10.1111/tpj.70611

    Comparing closely related species with differing abiotic stress tolerances can reveal physiological and genetic mechanisms that underlie environmental adaptations. Mangroves, found primarily in tropical and subtropical coasts, are expanding poleward due to rising minimum temperatures, for which species-specific cold tolerance likely influences expansion limits. In the Indo-West Pacific, Kandelia obovata tolerates chilling temperatures below 10°C, while its sister species, Kandelia candel, is highly sensitive to chilling and restricted to tropical regions. The discrepancy in chilling tolerance offers unique materials for understanding chilling tolerance mechanisms. This study compares chilling-induced changes in photochemical efficiency, enzymatic activity and gene expression between the two species to reveal their differences in photoprotective strategies. We specifically focus on the impact of chilling on the light reactions of photosynthesis as one of the primary impacts of chilling on plant species is chilling-induced photoinhibition. Despite maintaining higher photochemical efficiency, K. obovata exhibited greater reactive oxygen species (ROS) levels and lower antioxidant enzyme activity than K. candel during chilling at 10°C. Furthermore, interspecies comparison of the expression of genes involved in photoprotection during chilling stress revealed differential expression of PGR5 and CHL, which regulate cyclic electron flow and non-photochemical quenching. The upregulation of these genes in K. obovata and their downregulation in K. candel suggest that differences in photoprotective responses may contribute to the contrasting chilling tolerances of these species. These findings highlight a potential mechanism that contributes to the cold tolerance and biogeographical distributions of subtropical and tropical plant species, emphasizing the need for further comparative studies across a broader range of taxa.

  • Yang D, Jin Y, He N, Lin S, Cheng Z, Huang F, Zhang H, Li QQ, Yu W. OsMAPKKK69 Negatively Regulates Resistance to Blast and Bacterial Blight Diseases in Rice (Oryza sativa L.).. Plants (Basel, Switzerland). 2025;14(16). doi:10.3390/plants14162566

    Rice blast is one of the main diseases of rice, causing severe economic losses to agricultural production; thus, the search for blast resistance is a top priority for rice breeding. When challenged by the blast causal fungus Magnaporthe oryzae the expression level of OsMAPKKK69 gene in rice cultivar Nipponbar was found to increase significantly. Such an induction was also found in a different genetic material, cultivar Shufanggaonuo, indicating that OsMAPKKK69 plays an important role in blast disease response. However, the function of OsMAPKKK69 remains unclear. In this study, wild type ZH11 was selected as the background material to investigate the expression and functions of OsMAPKKK69 in rice disease resistance by constructing knockout mutants. The results showed that OsMAPKKK69 is mainly expressed in four-week-old shoots and localized in cell membrane, cytoplasm, and nucleus. The two allelic knockout mutants, osmapkkk69-1 and osmapkkk69-2, were more resistant to M. oryzae and bacterial blight Xanthomonas oryzae pv. Oryzae (Xoo). Further agronomic trait analysis revealed that the osmapkkk69-1 and osmapkkk69-2 mutants had reduced plant height, smaller grain size, a significant increase in tillering number, but also a significant increase in yield per plant. Our results show that OsMAPKKK69 is involved in the immune response of rice by negatively regulating the resistance to rice blast and blight diseases, and in regulating important agronomic traits. This study lays a foundation for revealing the molecular mechanism of OsMAPKKK69 in the immune response to rice diseases and provides novel genetic resources for rice breeding.

  • Ge S, Li J, Ma H, Sánchez-Bermúdez M, Wu S, Lv Q, Guo F, Dong J, Ma G, Li QQ, et al. A cis-natural antisense RNA regulates alternative polyadenylation of SlSPX5 under Pi starvation in tomato.. Nature communications. 2025;16(1):7981. doi:10.1038/s41467-025-63406-1

    Alternative polyadenylation (APA) generates transcript diversity by producing mRNA isoforms with distinct 3' ends. Despite the critical roles that APA plays in various biological processes, the mechanisms regulating APA in response to stresses have remained poorly understood in plants. Here, we perform comprehensive analysis of APA in tomato, and focus on a phosphate (Pi)- regulated APA gene SlSPX5, encoding a putative Pi sensor protein. SlSPX5 interacts with and sequesters the transcription factor SlPHL1 in the cytosol, thereby inhibiting the expression of Pi starvation inducible genes. We discover that a cis-natural antisense RNA (cis-NAT) is activated from SlSPX5 to promote its proximal polyadenylation under Pi-depleted conditions. The transcription of this cis-NAT induces RNA Polymerase II pausing, generating Ser2 phosphorylation signals that recruit polyadenylation machinery to the 5' end of SlSPX5. Our findings demonstrate that a cis-NAT regulates APA of its cognate gene in response to Pi starvation.