Caloric Restriction Reprograms the Single-Cell Transcriptional Landscape of Rattus Norvegicus Aging.

Aging causes a functional decline in tissues throughout the body that may be delayed by caloric restriction (CR). However, the cellular profiles and signatures of aging, as well as those ameliorated by CR, remain unclear. Here, we built comprehensive single-cell and single-nucleus transcriptomic atlases across various rat tissues undergoing aging and CR. CR attenuated aging-related changes in cell type composition, gene expression, and core transcriptional regulatory networks. Immune cells were increased during aging, and CR favorably reversed the aging-disturbed immune ecosystem. Computational prediction revealed that the abnormal cell-cell communication patterns observed during aging, including the excessive proinflammatory ligand-receptor interplay, were reversed by CR. Our work provides multi-tissue single-cell transcriptional landscapes associated with aging and CR in a mammal, enhances our understanding of the robustness of CR as a geroprotective intervention, and uncovers how metabolic intervention can act upon the immune system to modify the process of aging.

Very-low-density lipoprotein receptor-enhanced lipid metabolism in pancreatic stellate cells promotes pancreatic fibrosis.

Lipoprotein disorder is a common feature of chronic pancreatitis (CP); however, the relationship between lipoprotein disorder and pancreatic fibrotic environment is unclear. Here, we investigated the occurrence and mechanism of pancreatic stellate cell (PSC) activation by lipoprotein metabolites and the subsequent regulation of type 2 immune responses, as well as the driving force of fibrotic aggressiveness in CP. Single-cell RNA sequencing revealed the heterogeneity of PSCs and identified very-low-density lipoprotein receptor (VLDLR)+ PSCs that were characterized by a higher lipid metabolism. VLDLR promoted intracellular lipid accumulation, followed by interleukin-33 (IL-33) expression and release in PSCs. PSC-derived IL-33 strongly induced pancreatic group 2 innate lymphoid cells (ILC2s) to trigger a type 2 immune response accompanied by the activation of PSCs, eventually leading to fibrosis during pancreatitis. Our findings indicate that VLDLR-enhanced lipoprotein metabolism in PSCs promotes pancreatic fibrosis and highlight a dominant role of IL-33 in this pro-fibrotic cascade.

A PARylation-phosphorylation cascade promotes TOPBP1 loading and RPA-RAD51 exchange in homologous recombination.

The efficiency of homologous recombination (HR) in the repair of DNA double-strand breaks (DSBs) is closely associated with genome stability and tumor response to chemotherapy. While many factors have been functionally characterized in HR, such as TOPBP1, their precise regulation remains unclear. Here, we report that TOPBP1 interacts with the RNA-binding protein HTATSF1 in a cell-cycle- and phosphorylation-dependent manner. Mechanistically, CK2 phosphorylates HTATSF1 to facilitate binding to TOPBP1, which promotes S-phase-specific TOPBP1 recruitment to damaged chromatin and subsequent RPA/RAD51-dependent HR, genome integrity, and cancer-cell viability. The localization of HTATSF1-TOPBP1 to DSBs is potentially independent of the transcription-coupled RNA-binding and processing capacity of HTATSF1 but rather relies on the recognition of poly(ADP-ribosyl)ated RPA by HTATSF1, which can be blunted with PARP inhibitors. Together, our study provides a mechanistic insight into TOPBP1 loading at HR-prone DSB sites via HTATSF1 and reveals how RPA-RAD51 exchange is tuned by a PARylation-phosphorylation cascade.

CHIT1-positive microglia drive motor neuron ageing in the primate spinal cord.

Ageing is a critical factor in spinal-cord-associated disorders1, yet the ageing-specific mechanisms underlying this relationship remain poorly understood. Here, to address this knowledge gap, we combined single-nucleus RNA-sequencing analysis with behavioural and neurophysiological analysis in non-human primates (NHPs). We identified motor neuron senescence and neuroinflammation with microglial hyperactivation as intertwined hallmarks of spinal cord ageing. As an underlying mechanism, we identified a neurotoxic microglial state demarcated by elevated expression of CHIT1 (a secreted mammalian chitinase) specific to the aged spinal cords in NHP and human biopsies. In the aged spinal cord, CHIT1-positive microglia preferentially localize around motor neurons, and they have the ability to trigger senescence, partly by activating SMAD signalling. We further validated the driving role of secreted CHIT1 on MN senescence using multimodal experiments both in vivo, using the NHP spinal cord as a model, and in vitro, using a sophisticated system modelling the human motor-neuron-microenvironment interplay. Moreover, we demonstrated that ascorbic acid, a geroprotective compound, counteracted the pro-senescent effect of CHIT1 and mitigated motor neuron senescence in aged monkeys. Our findings provide the single-cell resolution cellular and molecular landscape of the aged primate spinal cord and identify a new biomarker and intervention target for spinal cord degeneration.

HALL: a comprehensive database for human aging and longevity studies.

Diverse individuals age at different rates and display variable susceptibilities to tissue aging, functional decline and aging-related diseases. Centenarians, exemplifying extreme longevity, serve as models for healthy aging. The field of human aging and longevity research is rapidly advancing, garnering significant attention and accumulating substantial data in recent years. Omics technologies, encompassing phenomics, genomics, transcriptomics, proteomics, metabolomics and microbiomics, have provided multidimensional insights and revolutionized cohort-based investigations into human aging and longevity. Accumulated data, covering diverse cells, tissues and cohorts across the lifespan necessitates the establishment of an open and integrated database. Addressing this, we established the Human Aging and Longevity Landscape (HALL), a comprehensive multi-omics repository encompassing a diverse spectrum of human cohorts, spanning from young adults to centenarians. The core objective of HALL is to foster healthy aging by offering an extensive repository of information on biomarkers that gauge the trajectory of human aging. Moreover, the database facilitates the development of diagnostic tools for aging-related conditions and empowers targeted interventions to enhance longevity. HALL is publicly available at https://ngdc.cncb.ac.cn/hall/index.

CK2-HTATSF1-TOPBP1 signaling axis modulates tumor chemotherapy response.

Homologous recombination (HR) plays a key role in maintaining genomic stability, and the efficiency of the HR system is closely associated with tumor response to chemotherapy. Our previous work reported that CK2 kinase phosphorylates HIV Tat-specific factor 1 (HTATSF1) Ser748 to facilitate HTATSF1 interaction with TOPBP1, which in turn, promotes RAD51 recruitment and HR repair. However, the clinical implication of the CK2-HTATSF1-TOPBP1 pathway in tumorigenesis and chemotherapeutic response remains to be elucidated. Here, we report that the CK2-HTATSF1-TOPBP1 axis is generally hyperactivated in multiple malignancies and renders breast tumors less responsive to chemotherapy. In contrast, deletion mutations of each gene in this axis, which also occur in breast and lung tumor samples, predict higher HR deficiency scores, and tumor cells bearing a loss-of-function mutation of HTATSF1 are vulnerable to poly(ADP-ribose) polymerase inhibitors or platinum drugs. Taken together, our study suggests that the integrity of the CK2-HTATSF1-TOPBP1 axis is closely linked to tumorigenesis and serves as an indicator of tumor HR status and modulates chemotherapy response.

Targeting RORγ inhibits the growth and metastasis of hepatocellular carcinoma.

Approximately 80%-90% of hepatocellular carcinomas (HCC) occur in a premalignant environment of fibrosis and abnormal extracellular matrix (ECM), highlighting an essential role of ECM in the tumorigenesis and progress of HCC. However, the determinants of ECM in HCC are poorly defined. Here, we show that nuclear receptor RORγ is highly expressed and amplified in HCC tumors. RORγ functions as an essential activator of the matrisome program via directly driving the expression of major ECM genes in HCC cells. Elevated RORγ increases fibronectin-1 deposition, cell-matrix adhesion, and collagen production, creating a favorable microenvironment to boost liver cancer metastasis. Moreover, RORγ antagonists effectively inhibit tumor growth and metastasis in multiple HCC xenografts and immune-intact models, and they effectively sensitize HCC tumors to sorafenib therapy in mice. Notably, elevated RORγ expression is associated with ECM remodeling and metastasis in patients with HCC. Taken together, we identify RORγ as a key player of ECM remodeling in HCC and as an attractive therapeutic target for advanced HCC.

Lineage Landscape: a comprehensive database that records lineage commitment across species.

Commitment to specific cell lineages is critical for mammalian embryonic development. Lineage determination, differentiation, maintenance, and organogenesis result in diverse life forms composed of multiple cell types. To understand the formation and maintenance of living individuals, including human beings, a comprehensive database that integrates multi-omic information underlying lineage differentiation across multiple species is urgently needed. Here, we construct Lineage Landscape, a database that compiles, analyzes and visualizes transcriptomic and epigenomic information related to lineage development in a collection of species. This landscape draws together datasets that capture the ongoing changes in cell lineages from classic model organisms to human beings throughout embryonic, fetal, adult, and aged stages, providing comprehensive, open-access information that is useful to researchers of a broad spectrum of life science disciplines. Lineage Landscape contains single-cell gene expression and bulk transcriptomic, DNA methylation, histone modifications, and chromatin accessibility profiles. Using this database, users can explore genes of interest that exhibit dynamic expression patterns at the transcriptional or epigenetic levels at different stages of lineage development. Lineage Landscape currently includes over 6.6 million cells, 15 million differentially expressed genes and 36 million data entries across 10 species and 34 organs. Lineage Landscape is free to access, browse, search, and download at http://data.iscr.ac.cn/lineage/#/home.

ADAR1 links R-loop homeostasis to ATR activation in replication stress response.

Unscheduled R-loops are a major source of replication stress and DNA damage. R-loop-induced replication defects are sensed and suppressed by ATR kinase, whereas it is not known whether R-loop itself is actively involved in ATR activation and, if so, how this is achieved. Here, we report that the nuclear form of RNA-editing enzyme ADAR1 promotes ATR activation and resolves genome-wide R-loops, a process that requires its double-stranded RNA-binding domains. Mechanistically, ADAR1 interacts with TOPBP1 and facilitates its loading on perturbed replication forks by enhancing the association of TOPBP1 with RAD9 of the 9-1-1 complex. When replication is inhibited, DNA-RNA hybrid competes with TOPBP1 for ADAR1 binding to promote the translocation of ADAR1 from damaged fork to accumulate at R-loop region. There, ADAR1 recruits RNA helicases DHX9 and DDX21 to unwind R-loops, simultaneously allowing TOPBP1 to stimulate ATR more efficiently. Collectively, we propose that the tempo-spatially regulated assembly of ADAR1-nucleated protein complexes link R-loop clearance and ATR activation, while R-loops crosstalk with blocked replication forks by transposing ADAR1 to finetune ATR activity and safeguard the genome.

A Universal Biomacromolecule-Enabled Assembly Strategy for Constructing Multifunctional Aerogels with 90% Inorganic Mass Loading from Inert Nano-Building Blocks.

Inert inorganic nano-building blocks, such as carbon nanotubes (CNTs) and boron nitride (BN) nanosheets, possess excellent physicochemical properties. However, it remains challenging to build aerogels with these inert nanomaterials unless they are chemically modified or compounded with petrochemical polymers, which affects their intrinsic properties and is usually not environmentally friendly. Here, a universal biomacromolecule-enabled assembly strategy is proposed to construct aerogels with 90 wt% ultrahigh inorganic loading. The super-high inorganic content is beneficial for exploiting the inherent properties of inert nanomaterials in multifunctional applications. Taking chitosan-CNTs aerogel as a proof-of-concept demonstration, it delivers sensitive pressure response as a pressure sensor, an ultrahigh sunlight absorption (94.5%) raising temperature under light (from 25 to 71 °C within 1 min) for clean-up of crude oil spills, and superior electromagnetic interference shielding performance of up to 68.9 dB. This strategy paves the way for the multifunctional application of inert nanomaterials by constructing aerogels with ultrahigh inorganic loading.

Unveiling the landscape predictors of resilient vegetation in coastal wetlands to inform conservation in the face of climate extremes.

Unveiling spatial variation in vegetation resilience to climate extremes can inform effective conservation planning under climate change. Although many conservation efforts are implemented on landscape scales, they often remain blind to landscape variation in vegetation resilience. We explored the distribution of drought-resilient vegetation (i.e., vegetation that could withstand and quickly recover from drought) and its predictors across a heterogeneous coastal landscape under long-term wetland conversion, through a series of high-resolution satellite image interpretations, spatial analyses, and nonlinear modelling. We found that vegetation varied greatly in drought resilience across the coastal wetland landscape and that drought-resilient vegetation could be predicted with distances to coastline and tidal channel. Specifically, drought-resilient vegetation exhibited a nearly bimodal distribution and had a seaward optimum at ~2 km from coastline (corresponding to an inundation frequency of ~30%), a pattern particularly pronounced in areas further away from tidal channels. Furthermore, we found that areas with drought-resilient vegetation were more likely to be eliminated by wetland conversion. Even in protected areas where wetland conversion was slowed, drought-resilient vegetation was increasingly lost to wetland conversion at its landward optimum in combination with rapid plant invasions at its seaward optimum. Our study highlights that the distribution of drought-resilient vegetation can be predicted using landscape features but without incorporating this predictive understanding, conservation efforts may risk failing in the face of climate extremes.

Polymer-based cascaded waveguide Bragg grating for blood glucose sensing.

Waveguide Bragg grating (WBG) blood glucose sensing, as a biological sensing technology with broad application prospects, plays an important role in the fields of health management and medical treatment. In this work, a polymer-based cascaded WBG is applied to glucose detection. We investigated photonic devices with two different grating structures cascaded-a crossed grating and a bilateral grating-and analyzed the effects of the crossed grating period, bilateral grating period, and number of grating periods on the sensing performance of the glucose sensor. Finally, the spectral reflectance characteristics, response time, and sensing specificity of the cascaded WBG were evaluated. The experimental results showed that the glucose sensor has a sensitivity of 175 nm/RIU in a glucose concentration range of 0-2 mg/ml and has the advantages of high integration, a narrow bandwidth, and low cost.

Repurposing endogenous type II CRISPR-Cas9 system for genome editing in Streptococcus thermophilus.

Streptococcus thermophilus has been extensively used in industrial milk fermentation. However, lack of efficient genetic manipulation approaches greatly hampered the industrial application of this species. Here, we repurposed the endogenous CRISPR1 and CRISPR3 systems, both belong to type II-A CRISPR-Cas9, by delivering a self-targeting CRISPR array with DNA repair template into S. thermophilus LMD-9. We achieved 785-bp deletion in lacZ gene by repurposing CRISPR1 and CRISPR3 systems with efficiencies of 35% and 59%, respectively, when 1-kb DNA repair template was provided. While providing with 1.5-kb repair template, the editing efficiency for deletion in lacZ gene reached 90% using CRISPR3 systems. Diverse editing outcomes encompassing a stop code insertion and single nucleotide variation within lacZ, as well as a 234-bp DNA fragment insertion upstream of ster_0903, were generated with high efficiencies of 75%-100% using the CRISPR3 system. Harnessing the customized endogenous CRISPR3 system to target six genes of eps gene cluster, we obtained six single-gene knockout mutants with efficiencies of 29%-80%, and proved that the epsA, epsE, and epsG were the key genes affecting exopolysaccharides biosynthesis in S. thermophilus LMD-9. Altogether, repurposing the native type II-A CRISPR-Cas9 can be served as a toolkit for precise genome engineering in S. thermophilus for biotechnological applications.

Enhancing Phosphorus Release from Sewage Sludge in Anaerobic Digestion via Thermal Hydrolysis Pretreatment: Insights from Phosphorus Speciation and Molecular Biological Pathways.

This study explores the mechanisms enhancing phosphorus (P) release from sludge in anaerobic digestion (AD) with thermal hydrolysis pretreatment (THP) using sequential chemical extraction, X-ray absorption near-edge structure spectroscopy (XANES), 31P NMR, and multiomics. THP-treated sludge notably increased liquid-phase P by 53.8% over 3 days compared to sewage sludge (SS), identifying solid-phase Fe-P as the primary P source. The THP+AD also provided a higher abundance of bacteria that contributed to P release through multiple pathways (MPRPB), whereas SS+AD enriched some microbial species with single P release pathway. Moreover, species co-occurrence network analysis underlined the pivotal role of P-releasing bacteria in THP+AD, with 8 out of 16 keystones being P-releasers. Among the 63 screened genes that were related to P transformations and release, the poly beta-hydroxybutyrate (PHB) synthesis genes associated with polyphosphate bacteria-mediated P release were more abundant in THP+AD than in SS+AD. Furthermore, the upregulation of genes involved in methyl phosphonate metabolism in the THP-treated sludge enhanced the methane production potential of the AD process. These findings suggested that MPRPB were indeed the main contributors to P release, and enrichment in the THP+AD process enhanced their capability for P liberation.

Deciphering the Impact of ZnO Nanoparticles and a Sunscreen Product Containing ZnO on Phosphorus Dynamics and Release in Chlorella pyrenoidosa in Aquatic Systems.

Zinc oxide nanoparticles (ZnO NPs) expedite the conversion of organic phosphorus (OP) into PO4-P (Pi), facilitating phosphorus (P) absorption by algae. Our study explored the mechanisms of converting OP (2-aminoethylphosphonic acid (AEP) and ß-glycerol phosphate (ß-GP)) into Pi in Chlorella pyrenoidosa under P deficiency with sunscreen and ZnO NPs. Cell density followed the order of K2HPO4 > ß-GP+ZnO > ß-GP > AEP+ZnO > AEP > P-free. ZnO NPs promoted the conversion of ß-GP, containing C-O-P bonds (0.028-0.041 mg/L), into Pi more efficiently than AEP, which possesses C-P bonds (0.022-0.037 mg/L). Transcriptomics revealed Pi transport/metabolism (phoB (3.99-12.01 fold), phoR (2.20-5.50 fold), ppa (4.49-10.40 fold), and ppk (2.50-5.40 fold)) and phospholipid metabolism (SQD1 (1.85-2.79 fold), SQD2 (2.60-6.53 fold), MGD (2.13-3.21 fold), and DGD (4.08-7.56 fold)) were up-regulated compared to K2HPO4. 31P nuclear magnetic resonance spectroscopy identified intracellular P as polyphosphate, orthophosphate, and pyrophosphate. Synchrotron radiation-based X-ray near-edge structure spectroscopy indicated that K2HPO4 and Zn3(PO4)2 in ß-GP+ZnO were increased by 8.09% and 7.28% compared to AEP+ZnO, suggesting superior P storage in ß-GP+ZnO. Overall, ZnO NPs improved photoinduced electron-hole pair separation and charge separation efficiency and amplified the ·OH and ·O2- levels, promoting OP photoconversion into Pi and algae growth.

Association between exposure to per- and perfluoroalkyl substances (PFAS) and reproductive hormones in human: A systematic review and meta-analysis.

BACKGROUND: Per- and polyfluoroalkyl substances (PFAS) is persistent endocrine disrupting chemicals. Previous evidence suggests that exposure to PFAS is associated with reproductive hormone levels, but the results of relevant studies are inconsistent. The objective of our study is to determine the association between exposure to PFAS and reproductive hormone levels in gender-specific general population. METHOD: Based on scientific search strategies, we systematically searched PubMed, Web of Science, Embase, Medline, and Scopus to obtain the eligible studies published before January 21, 2023. The quality of the included articles was assessed using the Office of Health Assessment and Translation (OHAT) Risk of Bias tool. We combined the ß coefficient and 95% confidence intervals (CI) using Stata.17 with random-effect model or fixed-effect model. We also performed subgroup analysis, sensitivity analysis, and Begger's and Egger's tests. RESULTS: Eleven studies involving 7714 participants were included. Meta-analysis showed that PFHxS exposure was positively associated with estradiol (E2) levels in female [ß = 0.030, 95% CI: (0.013, 0.046), P = 0.000]. A negative association was found between PFOA [ß = -0.012, 95% CI: (-0.023, -0.002), P = 0.017] and PFOS [ß = -0.011; 95% CI: (-0.021, -0.000), P = 0.042] exposure with male testosterone (TT) levels. In the subgroup analysis, there were stronger associations in children than in adults. And the high heterogeneity was mainly due to the cross-sectional studies. Publication bias was not found in most of the analyses. CONCLUSION: Our study showed that PFAS exposure was significantly associated with reproductive hormone levels. Further related studies are needed to identify the association and potential mechanism in the future.

Regeneration Roadmap: database resources for regenerative biology.

Regeneration plays an instrumental role in biological development and damage repair by constructing and replacing cells, tissues, and organs. Since regenerative capacity declines with age, promoting regeneration is heralded as a potential strategy for delaying aging. On this premise, mechanisms that regulate regeneration have been extensively studied across species and in different tissues. However, an open and comprehensive database collecting and standardizing the abundant data generated in regeneration research, such as high-throughput sequencing data, remains to be developed. In this work, we constructed Regeneration Roadmap to systematically and comprehensively collect such information over 2.38 million data entries across 11 species and 36 tissues, including regeneration-related genes, bulk and single-cell transcriptomics, epigenomics, and pharmacogenomics data. In this database, users can explore regulatory and expression changes of regeneration-associated genes in different species and tissues. Regeneration Roadmap provides the research community with a long-awaited and valuable data resource featuring convenient computing and visualizing tools, which is publicly available at https://ngdc.cncb.ac.cn/regeneration/index.

Investigation of radiomics models for predicting biochemical recurrence of advanced prostate cancer on pretreatment MR ADC maps based on automatic image segmentation.

OBJECTIVES: To develop radiomics models based on automatic segmentation of the pretreatment apparent diffusion coefficient (ADC) maps for predicting the biochemical recurrence (BCR) of advanced prostate cancer (PCa). METHODS: A total of 100 cases with pathologically confirmed PCa were retrospectively included in this study. These cases were randomly divided into training (n = 70) and test (n = 30) datasets. Two predictive models were constructed based on the combination of age, prostate specific antigen (PSA) level, Gleason score, and clinical staging before therapy and the prostate area (Model_1) or PCa area (Model_2). Another two predictive models were constructed based on only prostate area (Model_3) or PCa area (Model_4). The area under the receiver operating characteristic curve (ROC AUC) and precision-recall (PR) curve analysis were used to analyze the models' performance. RESULTS: Sixty-five patients without BCR (BCR-) and 35 patients with BCR (BCR+) were confirmed. The age, PSA, volume, diameter and ADC value of the prostate and PCa were not significantly different between the BCR- and BCR+ groups or between the training and test datasets (all p > 0.05). The AUCs were 0.637 (95% CI: 0.434-0.838), 0.841 (95% CI: 0.695-0.940), 0.840 (95% CI: 0.698-0.983), and 0.808 (95% CI: 0.627-0.988) for Model_1 to Model_4 in the test dataset without significant difference. The 95% bootstrap confidence intervals for the areas under the PR curve of the four models were not statistically different. CONCLUSION: The radiomics models based on automatically segmented prostate and PCa areas on the pretreatment ADC maps developed in our study can be promising in predicting BCR of advanced PCa.

Remote Sensing Monitoring of Grassland Locust Density Based on Machine Learning.

The main aim of this study was to utilize remote sensing data to establish regression models through machine learning to predict locust density in the upcoming year. First, a dataset for monitoring grassland locust density was constructed based on meteorological data and multi-source remote sensing data in the study area. Subsequently, an SVR (support vector regression) model, BP neural network regression model, random forest regression model, BP neural network regression model with the PCA (principal component analysis), and deep belief network regression model were built on the dataset. The experimental results show that the random forest regression model had the best prediction performance among the five models. Specifically, the model achieved a coefficient of determination (R2) of 0.9685 and a root mean square error (RMSE) of 1.0144 on the test set, which were the optimal values achieved among all the models tested. Finally, the locust density in the study area for 2023 was predicted and, by comparing the predicted results with actual measured data, it was found that the prediction accuracy was high. This is of great significance for local grassland ecological management, disaster warning, scientific decision-making support, scientific research progress, and sustainable agricultural development.

Telecoupling between urban expansion and forest ecosystem service loss through cultivated land displacement: A case study of Zhejiang Province, China.

Urbanization can either directly occupy forests or indirectly lead to forest loss elsewhere through cultivated land displacement, resulting in further forest fragmentation and ecosystem service (ES) loss. However, the effects of urban expansion on forest area and ESs are unknown, and this is especially true for indirect effects. Taking Zhejiang Province, China, a typical deforested province, as an example, this study quantified the direct and indirect effects of urban expansion on forest area and five ESs (timber yield, water yield, carbon sequestration, soil conservation, and biodiversity) from 2000 to 2020, explored the relationship between forest structure (forest proportion, mean patch area, edge density, and mean euclidean nearest neighbor distance) change and ESs, and revealed the telecoupling of urban expansion and forest loss and cascade effects among urbanization, deforestation, forest structure, and ESs. The results indicated that the indirect forest loss (4.30%-6.15%) caused by cultivated land displacement due to urban expansion was larger than the direct forest loss (2.42%). Urban expansion has a greater negative impact on carbon sequestration (6.40%-8.20%), water yield (6.08%-7.78%), and biodiversity (5.79%-7.44%) than on timber yield (4.77%-6.17%) and soil conservation (4.43%-5.77%). The indirect forest ES loss was approximately 2.83-4.34 times greater than the direct forest ES loss. Most forest ESs showed a nonlinear significant positive correlation with changes in forest proportion and mean patch area and a significant nonlinear negative correlation with changes in edge density and mean Euclidean nearest neighbor distance (p < 0.05). There is telecoupling between urban expansion in one region and forest ES loss in other distant regions. This study contributes to guiding sustainable forest conservation and management globally.