Genomes of the Banyan Tree and Pollinator Wasp Provide Insights into Fig-Wasp Coevolution.

Banyan trees are distinguished by their extraordinary aerial roots. The Ficus genus includes species that have evolved a species-specific mutualism system with wasp pollinators. We sequenced genomes of the Chinese banyan tree, F. microcarpa, and a species lacking aerial roots, F. hispida, and one wasp genome coevolving with F. microcarpa, Eupristina verticillata. Comparative analysis of the two Ficus genomes revealed dynamic karyotype variation associated with adaptive evolution. Copy number expansion of auxin-related genes from duplications and elevated auxin production are associated with aerial root development in F. microcarpa. A male-specific AGAMOUS paralog, FhAG2, was identified as a candidate gene for sex determination in F. hispida. Population genomic analyses of Ficus species revealed genomic signatures of morphological and physiological coadaptation with their pollinators involving terpenoid- and benzenoid-derived compounds. These three genomes offer insights into and genomic resources for investigating the geneses of aerial roots, monoecy and dioecy, and codiversification in a symbiotic system.

Benchmarking spatial clustering methods with spatially resolved transcriptomics data.

Spatial clustering, which shares an analogy with single-cell clustering, has expanded the scope of tissue physiology studies from cell-centroid to structure-centroid with spatially resolved transcriptomics (SRT) data. Computational methods have undergone remarkable development in recent years, but a comprehensive benchmark study is still lacking. Here we present a benchmark study of 13 computational methods on 34 SRT data (7 datasets). The performance was evaluated on the basis of accuracy, spatial continuity, marker genes detection, scalability, and robustness. We found existing methods were complementary in terms of their performance and functionality, and we provide guidance for selecting appropriate methods for given scenarios. On testing additional 22 challenging datasets, we identified challenges in identifying noncontinuous spatial domains and limitations of existing methods, highlighting their inadequacies in handling recent large-scale tasks. Furthermore, with 145 simulated data, we examined the robustness of these methods against four different factors, and assessed the impact of pre- and postprocessing approaches. Our study offers a comprehensive evaluation of existing spatial clustering methods with SRT data, paving the way for future advancements in this rapidly evolving field.

The roles of nuclear orphan receptor NR2F6 in anti-viral innate immunity.

Proper transcription regulation by key transcription factors, such as IRF3, is critical for anti-viral defense. Dynamics of enhancer activity play important roles in many biological processes, and epigenomic analysis is used to determine the involved enhancers and transcription factors. To determine new transcription factors in anti-DNA-virus response, we have performed H3K27ac ChIP-Seq and identified three transcription factors, NR2F6, MEF2D and MAFF, in promoting HSV-1 replication. NR2F6 promotes HSV-1 replication and gene expression in vitro and in vivo, but not dependent on cGAS/STING pathway. NR2F6 binds to the promoter of MAP3K5 and activates AP-1/c-Jun pathway, which is critical for DNA virus replication. On the other hand, NR2F6 is transcriptionally repressed by c-Jun and forms a negative feedback loop. Meanwhile, cGAS/STING innate immunity signaling represses NR2F6 through STAT3. Taken together, we have identified new transcription factors and revealed the underlying mechanisms involved in the network between DNA viruses and host cells.

LAFEM: A Scoring Model to Evaluate Functional Landscape of Lysine Acetylome.

Protein lysine acetylation is a critical post-translational modification involved in a wide range of biological processes. To date, about 20,000 acetylation sites of Homo sapiens were identified through mass spectrometry-based proteomic technology, but more than 95% of them have unclear functional annotations because of the lack of existing prioritization strategy to assess the functional importance of the acetylation sites on large scale. Hence, we established a lysine acetylation functional evaluating model (LAFEM) by considering eight critical features surrounding lysine acetylation site to high-throughput estimate the functional importance of given acetylation sites. This was achieved by selecting one of the random forest models with the best performance in 10-fold cross-validation on undersampled training dataset. The global analysis demonstrated that the molecular environment of acetylation sites with high acetylation functional scores (AFSs) mainly had the features of larger solvent-accessible surface area, stronger hydrogen bonding-donating abilities, near motif and domain, higher homology, and disordered degree. Importantly, LAFEM performed well in validation dataset and acetylome, showing good accuracy to screen out fitness directly relevant acetylation sites and assisting to explain the core reason for the difference between biological models from the perspective of acetylome. We further used cellular experiments to confirm that, in nuclear casein kinase and cyclin-dependent kinase substrate 1, acetyl-K35 with higher AFS was more important than acetyl-K9 with lower AFS in the proliferation of A549 cells. LAFEM provides a prioritization strategy to large scale discover the fitness directly relevant acetylation sites, which constitutes an unprecedented resource for better understanding of functional acetylome.

ARL6IP1 mediates small-molecule-induced alleviation of Alzheimer pathology through FXR1-dependent BACE1 translation initiation.

Exploring the potential lead compounds for Alzheimer's disease (AD) remains one of the challenging tasks. Here, we report that the plant extract conophylline (CNP) impeded amyloidogenesis by preferentially inhibiting BACE1 translation via the 5' untranslated region (5'UTR) and rescued cognitive decline in an animal model of APP/PS1 mice. ADP-ribosylation factor-like protein 6-interacting protein 1 (ARL6IP1) was then found to mediate the effect of CNP on BACE1 translation, amyloidogenesis, glial activation, and cognitive function. Through analysis of the 5'UTR-targetd RNA-binding proteins by RNA pulldown combined with LC-MS/MS, we found that FMR1 autosomal homolog 1 (FXR1) interacted with ARL6IP1 and mediated CNP-induced reduction of BACE1 by regulating the 5'UTR activity. Without altering the protein levels of ARL6IP1 and FXR1, CNP treatment promoted ARL6IP1 interaction with FXR1 and inhibited FXR1 binding to the 5'UTR both in vitro and in vivo. Collectively, CNP exhibited a therapeutic potential for AD via ARL6IP1. Through pharmacological manipulation, we uncovered a dynamic interaction between FXR1 and the 5'UTR in translational control of BACE1, adding to the understanding of the pathophysiology of AD.

AP2S1 regulates APP degradation through late endosome-lysosome fusion in cells and APP/PS1 mice.

AP2S1 is the sigma 2 subunit of adaptor protein 2 (AP2) that is essential for endocytosis. In this study, we investigated the potential role of AP2S1 in intracellular processing of amyloid precursor protein (APP), which contributes to the pathogenesis of Alzheimer disease (AD) by generating the toxic ß-amyloid peptide (Aß). We found that knockdown or overexpression of AP2S1 decreased or increased the protein levels of APP and Aß in cells stably expressing human full-length APP695, respectively. This effect was unrelated to endocytosis but involved lysosomal degradation. Morphological studies revealed that silencing of AP2S1 promoted the translocalization of APP from RAB9-positive late endosomes (LE) to LAMP1-positive lysosomes, which was paralleled by the enhanced LE-lysosome fusion. In support, silencing of vacuolar protein sorting-associated protein 41 (VPS41) that is implicated in LE-lyso fusion prevented AP2S1-mediated regulation of APP degradation and translocalization. In APP/PS1 mice, an animal model of AD, AAV-mediated delivery of AP2S1 shRNA in the hippocampus significantly reduced the protein levels of APP and Aß, with the concomitant APP translocalization, LE-lyso fusion and the improved cognitive functions. Taken together, these data uncover a LE-lyso fusion mechanism in APP degradation and suggest a novel role for AP2S1 in the pathophysiology of AD.

Treponema pallidum promoted microglia apoptosis and prevented itself from clearing by human microglia via blocking autophagic flux.

Treponema pallidum (Tp) has a well-known ability to evade the immune system and can cause neurosyphilis by invading the central nervous system (CNS). Microglia are resident macrophages of the CNS that are essential for host defense against pathogens, this study aims to investigate the interaction between Tp and microglia and the potential mechanism. Here, we found that Tp can exert significant toxic effects on microglia in vivo in Tg (mpeg1: EGFP) transgenic zebrafish embryos. Single-cell RNA sequencing results showed that Tp downregulated autophagy-related genes in human HMC3 microglial cells, which is negatively associated with apoptotic gene expression. Biochemical and cell biology assays further established that Tp inhibits microglial autophagy by interfering with the autophagosome-lysosome fusion process. Transcription factor EB (TFEB) is a master regulator of lysosome biogenesis, Tp activates the mechanistic target of rapamycin complex 1 (mTORC1) signaling to inhibit the nuclear translocation of TFEB, leading to decreased lysosomal biogenesis and accumulated autophagosome. Importantly, the inhibition of autophagosome formation reversed Tp-induced apoptosis and promoted microglial clearance of Tp. Taken together, these findings show that Tp blocks autophagic flux by inhibiting TFEB-mediated lysosomal biosynthesis in human microglia. Autophagosome accumulation was demonstrated to be a key mechanism underlying the effects of Tp in promoting apoptosis and preventing itself from clearing by human microglia. This study offers novel perspectives on the potential mechanism of immune evasion employed by Tp within CNS. The results not only establish the pivotal role of autophagy dysregulation in the detrimental effects of Tp on microglial cells but also bear considerable implications for the development of therapeutic strategies against Tp, specifically involving mTORC1 inhibitors and autophagosome formation inhibitors, in the context of neurosyphilis patients.

Dense Bicoid hubs accentuate binding along the morphogen gradient.

Morphogen gradients direct the spatial patterning of developing embryos; however, the mechanisms by which these gradients are interpreted remain elusive. Here we used lattice light-sheet microscopy to perform in vivo single-molecule imaging in early Drosophila melanogaster embryos of the transcription factor Bicoid that forms a gradient and initiates patterning along the anteroposterior axis. In contrast to canonical models, we observed that Bicoid binds to DNA with a rapid off rate throughout the embryo such that its average occupancy at target loci is on-rate-dependent. We further observed Bicoid forming transient "hubs" of locally high density that facilitate binding as factor levels drop, including in the posterior, where we observed Bicoid binding despite vanishingly low protein levels. We propose that localized modulation of transcription factor on rates via clustering provides a general mechanism to facilitate binding to low-affinity targets and that this may be a prevalent feature of other developmental transcription factors.

Identifying the kind behind SMILES-anatomical therapeutic chemical classification using structure-only representations.

Anatomical Therapeutic Chemical (ATC) classification for compounds/drugs plays an important role in drug development and basic research. However, previous methods depend on interactions extracted from STITCH dataset which may make it depend on lab experiments. We present a pilot study to explore the possibility of conducting the ATC prediction solely based on the molecular structures. The motivation is to eliminate the reliance on the costly lab experiments so that the characteristics of a drug can be pre-assessed for better decision-making and effort-saving before the actual development. To this end, we construct a new benchmark consisting of 4545 compounds which is with larger scale than the one used in previous study. A light-weight prediction model is proposed. The model is with better explainability in the sense that it is consists of a straightforward tokenization that extracts and embeds statistically and physicochemically meaningful tokens, and a deep network backed by a set of pyramid kernels to capture multi-resolution chemical structural characteristics. Its efficacy has been validated in the experiments where it outperforms the state-of-the-art methods by 15.53% in accuracy and by 69.66% in terms of efficiency. We make the benchmark dataset, source code and web server open to ease the reproduction of this study.

Establishment of an in vitro model of monocyte-like THP-1 cells for trained immunity induced by bacillus Calmette-Guérin.

BACKGROUND: Mycobacteria bloodstream infections are common in immunocompromised people and usually have disastrous consequences. As the primary phagocytes in the bloodstream, monocytes and neutrophils play critical roles in the fight against bloodstream mycobacteria infections. In contrast to macrophages, the responses of monocytes infected with the mycobacteria have been less investigated. RESULTS: In this study, we first established a protocol for infection of non-adherent monocyte-like THP-1 cells (i.e. without the differentiation induced by phorbol 12-myristate 13-acetate (PMA) by bacillus Calmette-Guérin (BCG). Via the protocol, we were then capable of exploring the global transcriptomic profiles of non-adherent THP-1 cells infected with BCG, and found that NF-κB, MAPK and PI3K-Akt signaling pathways were enhanced, as well as some inflammatory chemokine/cytokine genes (e.g. CCL4, CXCL10, TNF and IL-1ß) were up-regulated. Surprisingly, the Akt-HIF-mTOR signaling pathway was also activated, which induces trained immunity. In this in vitro infection model, increased cytokine responses to lipopolysaccharides (LPS) restimulation, higher cell viability, and decreased Candida albicans loads were observed. CONCLUSIONS: We have first characterized the transcriptomic profiles of BCG-infected non-adherent THP-1 cells, and first developed a trained immunity in vitro model of the cells.

Ligand binding properties of three odorant-binding proteins in striped flea beetle Phyllotreta striolata towards two phthalate esters.

Odorant-binding proteins (OBPs) initiate insect olfactory perception and mediate specific binding and selection of odorants via uncertain binding mechanisms. We characterized the binding characteristics of four OBPs from the striped flea beetle Phyllotreta striolata (SFB), a major cruciferous crop pest. Tissue expression analysis revealed that the two ABPII OBPs (PstrOBP12 and PstrOBP19) were highly expressed mainly in the antenna, whereas the two minus-C OBPs (PstrOBP13 and PstrOBP16) showed a broad expression pattern. Competitive binding assays of cruciferous plant volatiles showed that PstrOBP12, PstrOBP16 and PstrOBP19 had very strong binding capacities for only two phthalate esters (Ki < 20 µM), and PstrOBP13 specifically bound to four aromatic volatiles (Ki < 11 µM). Fluorescence quenching assays displayed that two phthalate esters bound to three PstrOBPs via different quenching mechanisms. PstrOBP12/PstrOBP16-diisobutyl phthalate and PstrOBP19-bis(6-methylheptyl) phthalate followed static quenching, while PstrOBP12/PstrOBP16-bis(6-methylheptyl) phthalate and PstrOBP19-diisobutyl phthalate followed dynamic quenching. Homology modelling and molecular docking displayed that PstrOBP12-diisobutyl phthalate was driven by H-bonding and van der Waals interactions, while PstrOBP16-diisobutyl phthalate and PstrOBP19-bis(6-methylheptyl) phthalate followed hydrophobic interactions. Finally, behavioural activity analysis demonstrated that phthalate esters exhibited different behavioural activities of SFB at different doses, with low doses attracting and high doses repelling. Overall, we thus revealed the different binding properties of the three PstrOBPs to two phthalate esters, which was beneficial in shedding light on the ligand-binding mechanisms of OBPs.

Interference-induced generation of a chirp-free short isolated attosecond pulse in the water window region with multicolor laser fields.

Compensating for the intrinsic attosecond chirp (atto-chirp) of wideband high-order harmonics in the water window region is a significant challenge, in order to obtain isolated attosecond pulses (IAPs) with a width of tens of attoseconds (as). Here, we propose to realize the generation of IAP with duration as short as 20 as, central energy of 365 eV, and bandwidth exceeding 150 eV from chirp-free high harmonics generated by a four-color driving laser, without the necessity for atto-chirp compensation with natural materials. Unlike any other gating methods that an IAP arises from only one electron ionization event, we take advantage of the interference between harmonic radiation produced by multiple ionizing events. We further demonstrate that such chirp-free short IAP survives after taking account of macroscopic propagation effects. Given that the synthesized multicolor laser field can also effectively increase the harmonic flux, this work provides a practical way for experiments to generate the broad bandwidth chirp-free IAPs in the water window region.

The m6A methyltransferase METTL3 drives thyroid cancer progression and lymph node metastasis by targeting LINC00894.

BACKGROUND: Long noncoding RNAs (lncRNAs) are significant contributors to various human malignancies. The aberrant expression of lncRNA LINC00894 has been reported in various human malignancies. We aimed to illustrate the role of LINC00894 and its underlying mechanism in the development of papillary thyroid carcinoma (PTC). METHODS: We performed bioinformatics analysis of differentially expressed RNAs from TCGA and GEO datasets and selected the target lncRNA LINC00894. SRAMP analysis revealed abundant M6A modification sites in LINC00894. Further analysis of StarBase, GEPIA, and TCGA datasets was performed to identify the related differentially expressed genes METTL3. Colony formation and CCK-8 assays confirmed the relationship between LINC00894, METTL3, and the proliferative capacity of PTC cells. The analysis of AnnoLnc2, Starbase datasets, and meRIP-PCR and qRT‒PCR experiments confirmed the influence of METTL3-mediated m6A modification on LINC00894. The study employed KEGG enrichment analysis as well as Western blotting to investigate the impact of LINC00894 on the expression of proteins related to the Hippo signalling pathway. RESULTS: LINC00894 downregulation was detected in PTC tissues and cells and was even further downregulated in PTC with lymphatic metastasis. LINC00894 inhibits the lymphangiogenesis of vascular endothelial cells and the proliferation of cancer cells. METTL3 enhances PTC progression by upregulating LINC00894 by enhancing LINC00894 mRNA stability through the m6A-YTHDC2-dependent pathway. LINC00894 may inhibit PTC malignant phenotypes through the Hippo signalling pathway. CONCLUSION: The METTL3-YTHDC2 axis stabilizes LINC00894 mRNA in an m6A-dependent manner and subsequently inhibits tumour malignancy through the Hippo signalling pathway.

Extended Black Hole Thermodynamics from Extended Iyer-Wald Formalism.

In recent years, there has been significant interest in the field of extended black hole thermodynamics, where the cosmological constant and/or other coupling parameters are treated as thermodynamic variables. Drawing inspiration from the Iyer-Wald formalism, which reveals the intrinsic and universal structure of conventional black hole thermodynamics, we illustrate that a proper extension of this formalism also unveils the underlying theoretical structure of extended black hole thermodynamics. As a remarkable consequence, for any gravitational theory described by a diffeomorphism invariant action, it is always possible to construct a consistent extended thermodynamics using this extended formalism.

Preliminary study of finger temperature recovery in patients with diabetes mellitus following cold stimulation.

OBJECTIVE: To explore the difference in temperature recovery following cold stimulation between participants with and without diabetes mellitus (DM). MATERIALS AND METHODS: The participants without (control group; n = 25) and with (DM group; n = 26) DM were subjected to local cold stimulation (10º C for 90 s). The thermal images of their hands were continuously captured using a thermal camera within 7 min following cold stimulation, and the highest temperature of each fingertip was calculated. According to the temperature values at different timepoints, the temperature recovery curves were drawn, and the baseline temperature (T-base), initial temperature after cooling (T0), temperature decline amplitude (T-range), and area under the temperature recovery curve > T0 (S) were calculated. Finally, symmetry differences between the two groups were analysed. RESULTS: No statistical differences in the T-base, T0, and T-range were observed between the DM and control groups. After drawing the rewarming curve according to the temperature of the fingertips of the patients following cold stimulation, the S in the DM group was significantly lower than that in the control group (p < 0.05). Furthermore, the asymmetry of the base temperature of the hand was observed in the DM group. CONCLUSIONS: Following cold stimulation, the patients with DM exhibited a different rewarming pattern than those without DM. Thus, cold stimulation tests under infrared thermography may contribute to the early screening of diabetic peripheral neuropathy in future.

Probing the Electron Accepting Ability of Phosphaphenalenes.

Phosphaphenalenes, extended π conjugates with the incorporation of phosphorus, are attractive avenues towards molecular materials for the applications in organic electronics, but their electron accepting ability have not been investigated. In this study, we present systematic studies on the reductive behavior of a representative phosphaphenalene and its oxide by chemical and electrochemical methods. The chemical reduction of the phosphaphenalene by alkali metals reveals the facile P‒C bond cleavage to form phosphaphenalenide anion, which functions as a transfer block for structure modification on the phosphorus atom. In contrast, the pentavalent P-oxide reacts with one or two equivalents of elemental sodium to form stable radical anion and dianion salts, respectively.

pH Mapping of Gliomas Using Quantitative Chemical Exchange Saturation Transfer MRI: Quasi-Steady-State, Spillover-, and MT-Corrected Omega Plot Analysis.

BACKGROUND: Quantitative in-situ pH mapping of gliomas is important for therapeutic interventions, given its significant association with tumor progression, invasion, and metastasis. Although chemical exchange saturation transfer (CEST) offers a noninvasive way for pH imaging based on the pH-dependent exchange rate (ksw ), the reliable quantification of ksw in glioma remains constrained due to technical challenges. PURPOSE: To quantify the pH of gliomas by measuring the proton exchange rate through optimized omega plot analysis. STUDY TYPE: Prospective. PHANTOMS/ANIMAL MODEL/SUBJECTS: Creatine and murine brain lysates phantoms, six rats with glioma xenograft model, and three patients with World Health Organization grade 2-4 gliomas. FIELD STRENGTH/SEQUENCE: 11.7 T, 7.0 T, CEST imaging, T2 -weighted (T2 W) imaging, and T1 -mapping. ASSESSMENT: Omega plot analysis, quasi-steady-state (QUASS) analysis, multi-pool Lorentzian fitting, amine and amide concentration-independent detection, pH enhanced method with the combination of amide and guanidyl (pHenh ), and magnetization transfer ratio (MTR) were utilized for pH metric quantification. The clinical outcomes were determined through radiologic follow-up and histopathological analysis. STATISTICAL TESTS: Mann-Whitney U test was performed to compare glioma with normal tissue, and Pearson's correlation analysis was used to assess the relationship between ksw and other parameters. RESULTS: In vitro experiments reveal that the determined ksw at 2 ppm increases exponentially with pH (creatine phantoms: ksw = 106 + 0.147 × 10(pH-4.198) ; lysates: ksw = 185.1 + 0.101 × 10(pH-3.914) ). Omega plot analysis exhibits a linear correlation between 1/MTRRex and 1/ω1 2 in the glioma xenografts (R2 > 0.98) and glioma patients (R2 > 0.99). The exchange rate in the rat glioma decreases compared to the contralateral normal tissue (349.46 ± 30.40 s-1 vs. 403.54 ± 51.01 s-1 , P = 0.025), while keeping independence from changes in concentration (r = 0.5037, P = 0.095). Similar pattern was observed in human data. DATA CONCLUSION: Utilizing QUASS-based, spillover-, and MT-corrected omega plot analysis for the measurement of exchange rates, offers a feasible method for quantifying pH within glioma. LEVEL OF EVIDENCE: NA TECHNICAL EFFICACY: Stage 1.

Serological characteristics and clinical implications of IgG subclasses in visceral leishmaniasis.

OBJECTIVES: Visceral leishmaniasis (VL) represents the most severe form of Leishmaniasis infection, often resulting in fatality without timely treatment. Previous studies have found that immunosuppression increases the risk of VL disease progression and mortality, and the total immunoglobulin G (IgG) levels in peripheral blood vary before and after treatment. However, the distinct levels and roles of IgG subclasses in VL have not been documented yet. This study aims to elucidate the characteristics and clinical significance of IgG subclasses in VL. METHODS: A total of 43 cases newly-diagnosed with VL were enrolled in the cohort. We measured the levels of IgG subclasses before and after standard treatment and conducted assessments of bone marrow features. In addition, we analysed other haematological indices and examined the variations in IgG subclasses, as well as their correlation with clinical and laboratory factors. RESULTS: The levels of total IgG, IgG1, and the ratios of both IgG1/IgG and IgG1/IgG2 decreased significantly after treatment, whereas the ratios of IgG2/ IgG showed an obvious increase. The VL patients without hyperglobulinemia displayed significant lower IgG1/IgG2 ratios, but higher IgG2/IgG ratios compared with those with hyperglobulinemia. In addition, VL patients with positive bone marrow amastigotes had significant higher IgG1/IgG and IgG1/IgG2 ratios, but lower IgG2/IgG ratios. IgG subclasses were correlated with abnormal blood test results, particularly immunological elements including IgM and Complement 4 (C4). CONCLUSIONS: IgG1 and IgG2 exhibited contrasting changes after treatment in VL patients. The features of bone marrow and laboratory tests indicated that IgG1 and IgG2 serve different roles in the progression of VL. The ratios of IgG subclasses may be more precise indicators to evaluate immune reaction in VL than traditional total IgG.

Sodium thiosulfate: A donor or carrier signaling molecule for hydrogen sulfide?

Sodium thiosulfate has been used for decades in the treatment of calciphylaxis and cyanide detoxification, and has recently shown initial therapeutic promise in critical diseases such as neuronal ischemia, diabetes mellitus, heart failure and acute lung injury. However, the precise mechanism of sodium thiosulfate remains incompletely defined and sometimes contradictory. Although sodium thiosulfate has been widely accepted as a donor of hydrogen sulfide (H2S), emerging findings suggest that it is the executive signaling molecule for H2S and that its effects may not be dependent on H2S. This article presents an overview of the current understanding of sodium thiosulfate, including its synthesis, biological characteristics, and clinical applications of sodium thiosulfate, as well as the underlying mechanisms in vivo. We also discussed the interplay of sodium thiosulfate and H2S. Our review highlights sodium thiosulfate as a key player in sulfide signaling with the broad clinical potential for the future.

ALKBH5 facilitates the progression of infantile hemangioma by increasing FOXF1 expression in a m6A-YTHDF2 dependent manner to activate HK-2 signaling.

Alkylation repair homolog protein 5 (ALKBH5) is reported to participate in infantile hemangioma (IH) progression. However, the underlying mechanism of ALKBH5 in IH remains unclear. Using qRT-PCR and Western blotting, ALKBH5, forkhead box F1 (FOXF1) and hexokinase 2 (HK-2) expressions in IH tissues and IH-derived endothelial cells XPTS-1 were assessed. The Me-RIP assay was used to analyze FOXF1 m6A level. CCK8, colony formation, flow cytometry and transwell assays were employed to determine IH cell viability, proliferation, apoptosis, migration and invasion. The interactions between YTH (YT521-B homology) domain 2 (YTHDF2), FOXF1 and HK-2 were analyzed by RIP, dual luciferase reporter gene assay and/or ChIP assay. The in vivo IH growth was evaluated in immunocompromised mice. FOXF1 was overexpressed in IH tissues, and its silencing inhibited IH cell proliferation, migration and invasion whereas promoting cell apoptosis in vitro. ALKBH5 upregulation facilitated FOXF1 mRNA stability and expression in IH cells in a m6A-YTHDF2-dependent manner. FOXF1 downregulation reversed the impact of ALKBH5 upregulation on IH cellular phenotypes. It also turned out that FOXF1 positively regulated HK-2 expression in IH cells through interacting with the HK-2 promoter. HK-2 upregulation abolished FOXF1 knockdown's inhibition on IH cell aggressive behaviors. ALKBH5 or FOXF1 silencing suppressed IH tumor development via HK-2 signaling in immunocompromised mice. ALKBH5 promoted FOXF1 expression m6A-YTHDF2 dependently, which in turn elevated HK-2 expression, thereby accelerating IH development.