Publications

2026

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.

2025

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.

Cao Y, Guo Y, Yu Z, Nie H, Yang J, Qiu D, Li Q, Xin X, Cheng C, Li Y, et al. Polyadenylation Complex CFII Recognizes Downstream Cis-element for Pre-mRNA Polyadenylation Through Interaction with an RNA-Binding Protein in Arabidopsis.. Advanced science (Weinheim, Baden-Wurttemberg, Germany). 2025:e04562. doi:10.1002/advs.202504562

Alternative polyadenylation of pre-mRNA generates mRNAs with alternative 3' ends, thereby increasing the complexity of the transcriptome and proteome. This process is carried out by pre-mRNA 3' end processing complexes, with the specificity and site selection of polyadenylation primarily defined by the interaction between these complexes and conserved pre-mRNA cis-elements. However, the cis-element recognized by the polyadenylation CFII complex and the downstream cis-element of the poly(A) site remain unknown. Here, FPA, a conserved Spen family RNA-binding protein, is identified to physically interacts with the CFII complex in Arabidopsis. FPA specifically binds to a GA-rich cis-element downstream of the proximal poly(A) site, which is required for the CFII complex to promote proximal poly(A) site usage and prevent 3' extension. Thus, the CFII complex recognizes the downstream cis-element of the poly(A) site through its interaction with FPA, highlighting the role of RNA-binding proteins as accessory factors in alternative polyadenylation.

Jin X, Li J, Lu W, Deng X, Wei Y, Shu Y, Liu B, Liu Z, Long Y, Zhu X, et al. Pre-mRNA processing factors differentially impact coordination between co-transcriptional cleavage and transcription termination.. Nature communications. 2025;16(1):7086. doi:10.1038/s41467-025-62555-7

Co-transcriptional cleavage and transcription termination are closely related processes during mRNA maturation, yet their coordination remains poorly understood due to difficulties in detecting these transient events. Here, we applied single-molecule nascent RNA sequencing to simultaneously capture the cleavage status and readthrough distance on the same nascent RNA molecules and characterize 14 mutants of various pre-mRNA processing factors in Arabidopsis. Our results reveal diverse roles for these processing factors in coordinating cleavage and termination: core components of CPSF and CstF complex stimulate both cleavage and termination, facilitating access to exoribonuclease AtXRN3; mutations in nuclear poly(A) polymerase PAPS1 and AtXRN3 caused delayed termination with minimal effects on cleavage, suggesting their roles are further downstream; BORDER proteins facilitate termination while simultaneously inhibiting cleavage, suggesting a complex interplay between these two actions; the phosphatase SSU72 specifically promotes efficient termination without affecting cleavage. Our method also enables us to distinguish cleaved readthrough transcripts from full-length readthrough, and we found termination factor FPA specifically promotes termination of cleaved readthrough, suggesting FPA facilitates access of AtXRN3 to the 3' cleavage product. Our comprehensive datasets reveal cleavage and termination are highly coordinated during pre-mRNA processing.