Insights Gained from Single-Cell Analysis of Immune Cells in the Tumor Microenvironment.

Understanding tumor immune microenvironments is critical for identifying immune modifiers of cancer progression and developing cancer immunotherapies. Recent applications of single-cell RNA sequencing (scRNA-seq) in dissecting tumor microenvironments have brought important insights into the biology of tumor-infiltrating immune cells, including their heterogeneity, dynamics, and potential roles in both disease progression and response to immune checkpoint inhibitors and other immunotherapies. This review focuses on the advances in knowledge of tumor immune microenvironments acquired from scRNA-seq studies across multiple types of human tumors, with a particular emphasis on the study of phenotypic plasticity and lineage dynamics of immune cells in the tumor environment. We also discuss several imminent questions emerging from scRNA-seq observations and their potential solutions on the horizon.

Metabolic regulation of homologous recombination repair by MRE11 lactylation.

Lactylation is a lactate-induced post-translational modification best known for its roles in epigenetic regulation. Herein, we demonstrate that MRE11, a crucial homologous recombination (HR) protein, is lactylated at K673 by the CBP acetyltransferase in response to DNA damage and dependent on ATM phosphorylation of the latter. MRE11 lactylation promotes its binding to DNA, facilitating DNA end resection and HR. Inhibition of CBP or LDH downregulated MRE11 lactylation, impaired HR, and enhanced chemosensitivity of tumor cells in patient-derived xenograft and organoid models. A cell-penetrating peptide that specifically blocks MRE11 lactylation inhibited HR and sensitized cancer cells to cisplatin and PARPi. These findings unveil lactylation as a key regulator of HR, providing fresh insights into the ways in which cellular metabolism is linked to DSB repair. They also imply that the Warburg effect can confer chemoresistance through enhancing HR and suggest a potential therapeutic strategy of targeting MRE11 lactylation to mitigate the effects.

TFPI is a colonic crypt receptor for TcdB from hypervirulent clade 2 C. difficile.

The emergence of hypervirulent clade 2 Clostridioides difficile is associated with severe symptoms and accounts for >20% of global infections. TcdB is a dominant virulence factor of C. difficile, and clade 2 strains exclusively express two TcdB variants (TcdB2 and TcdB4) that use unknown receptors distinct from the classic TcdB. Here, we performed CRISPR/Cas9 screens for TcdB4 and identified tissue factor pathway inhibitor (TFPI) as its receptor. Using cryo-EM, we determined a complex structure of the full-length TcdB4 with TFPI, defining a common receptor-binding region for TcdB. Residue variations within this region divide major TcdB variants into 2 classes: one recognizes Frizzled (FZD), and the other recognizes TFPI. TFPI is highly expressed in the intestinal glands, and recombinant TFPI protects the colonic epithelium from TcdB2/4. These findings establish TFPI as a colonic crypt receptor for TcdB from clade 2 C. difficile and reveal new mechanisms for CDI pathogenesis.

Landscape of Infiltrating T Cells in Liver Cancer Revealed by Single-Cell Sequencing.

Systematic interrogation of tumor-infiltrating lymphocytes is key to the development of immunotherapies and the prediction of their clinical responses in cancers. Here, we perform deep single-cell RNA sequencing on 5,063 single T cells isolated from peripheral blood, tumor, and adjacent normal tissues from six hepatocellular carcinoma patients. The transcriptional profiles of these individual cells, coupled with assembled T cell receptor (TCR) sequences, enable us to identify 11 T cell subsets based on their molecular and functional properties and delineate their developmental trajectory. Specific subsets such as exhausted CD8+ T cells and Tregs are preferentially enriched and potentially clonally expanded in hepatocellular carcinoma (HCC), and we identified signature genes for each subset. One of the genes, layilin, is upregulated on activated CD8+ T cells and Tregs and represses the CD8+ T cell functions in vitro. This compendium of transcriptome data provides valuable insights and a rich resource for understanding the immune landscape in cancers.

COLD6-OSM1 module senses chilling for cold tolerance via 2',3'-cAMP signaling in rice.

While it is known that temperature sensors trigger calcium (Ca2+) signaling to confer cold tolerance in cells, less is known about sensors that couple with other secondary messengers. Here, we identify a cold sensor complex of CHILLING-TOLERANCE DIVERGENCE 6 (COLD6) and osmotin-like 1 (OSM1), which triggers 2',3'-cyclic adenosine monophosphate (2',3'-cAMP) production to enhance cold tolerance in rice. COLD6, which is encoded by a major quantitative trait locus (QTL) gene, interacts with the rice G protein α subunit (RGA1) at the plasma membrane under normal conditions. Upon exposure to chilling, cold-induced OSM1 binds to COLD6, kicking out RGA1 from interaction. This triggers an elevation of 2',3'-cAMP levels for enhancing chilling tolerance. Genetic data show that COLD6 negatively regulates cold tolerance and functionally depends on OSM1 in chilling stress. COLD6 alleles were selected during rice domestication. Knockout and natural variation of COLD6 in hybrid rice enhanced chilling tolerance, hinting design potential for breeding. This highlighted a module triggering 2',3'-cAMP to improve chilling tolerance in crops.

Structures of human γδ T cell receptor-CD3 complex.

Gamma delta (γδ) T cells, a unique T cell subgroup, are crucial in various immune responses and immunopathology1-3. The γδ T cell receptor (TCR), which is generated by γδ T cells, recognizes a diverse range of antigens independently of the major histocompatibility complex2. The γδ TCR associates with CD3 subunits, initiating T cell activation and holding great potential in immunotherapy4. Here we report the structures of two prototypical human Vγ9Vδ2 and Vγ5Vδ1 TCR-CD3 complexes5,6, revealing two distinct assembly mechanisms that depend on Vγ usage. The Vγ9Vδ2 TCR-CD3 complex is monomeric, with considerable conformational flexibility in the TCRγ-TCRδ extracellular domain and connecting peptides. The length of the connecting peptides regulates the ligand association and T cell activation. A cholesterol-like molecule wedges into the transmembrane region, exerting an inhibitory role in TCR signalling. The Vγ5Vδ1 TCR-CD3 complex displays a dimeric architecture, whereby two protomers nestle back to back through the Vγ5 domains of the TCR extracellular domains. Our biochemical and biophysical assays further corroborate the dimeric structure. Importantly, the dimeric form of the Vγ5Vδ1 TCR is essential for T cell activation. These findings reveal organizing principles of the γδ TCR-CD3 complex, providing insights into the unique properties of γδ TCR and facilitating immunotherapeutic interventions.

Mitochondrial complex I promotes kidney cancer metastasis.

Most kidney cancers are metabolically dysfunctional1-4, but how this dysfunction affects cancer progression in humans is unknown. We infused 13C-labelled nutrients in over 80 patients with kidney cancer during surgical tumour resection. Labelling from [U-13C]glucose varies across subtypes, indicating that the kidney environment alone cannot account for all tumour metabolic reprogramming. Compared with the adjacent kidney, clear cell renal cell carcinomas (ccRCCs) display suppressed labelling of tricarboxylic acid (TCA) cycle intermediates in vivo and in ex vivo organotypic cultures, indicating that suppressed labelling is tissue intrinsic. [1,2-13C]acetate and [U-13C]glutamine infusions in patients, coupled with measurements of respiration in isolated human kidney and tumour mitochondria, reveal lower electron transport chain activity in ccRCCs that contributes to decreased oxidative and enhanced reductive TCA cycle labelling. However, ccRCC metastases unexpectedly have enhanced TCA cycle labelling compared with that of primary ccRCCs, indicating a divergent metabolic program during metastasis in patients. In mice, stimulating respiration or NADH recycling in kidney cancer cells is sufficient to promote metastasis, whereas inhibiting electron transport chain complex I decreases metastasis. These findings in humans and mice indicate that metabolic properties and liabilities evolve during kidney cancer progression, and that mitochondrial function is limiting for metastasis but not growth at the original site.

The Arabidopsis SGN3/GSO1 receptor kinase integrates soil nitrogen status into shoot development.

The Casparian strip is a barrier in the endodermal cell walls of plants that allows the selective uptake of nutrients and water. In the model plant Arabidopsis thaliana, its development and establishment are under the control of a receptor-ligand mechanism termed the Schengen pathway. This pathway facilitates barrier formation and activates downstream compensatory responses in case of dysfunction. However, due to a very tight functional association with the Casparian strip, other potential signaling functions of the Schengen pathway remain obscure. In this work, we created a MYB36-dependent synthetic positive feedback loop that drives Casparian strip formation independently of Schengen-induced signaling. We evaluated this by subjecting plants in which the Schengen pathway has been uncoupled from barrier formation, as well as a number of established barrier-mutant plants, to agar-based and soil conditions that mimic agricultural settings. Under the latter conditions, the Schengen pathway is necessary for the establishment of nitrogen-deficiency responses in shoots. These data highlight Schengen signaling as an essential hub for the adaptive integration of signaling from the rhizosphere to aboveground tissues.

DNA-PKcs and ATM modulate mitochondrial ADP-ATP exchange as an oxidative stress checkpoint mechanism.

DNA-PKcs is a key regulator of DNA double-strand break repair. Apart from its canonical role in the DNA damage response, DNA-PKcs is involved in the cellular response to oxidative stress (OS), but its exact role remains unclear. Here, we report that DNA-PKcs-deficient human cells display depolarized mitochondria membrane potential (MMP) and reoriented metabolism, supporting a role for DNA-PKcs in oxidative phosphorylation (OXPHOS). DNA-PKcs directly interacts with mitochondria proteins ANT2 and VDAC2, and formation of the DNA-PKcs/ANT2/VDAC2 (DAV) complex supports optimal exchange of ADP and ATP across mitochondrial membranes to energize the cell via OXPHOS and to maintain MMP. Moreover, we demonstrate that the DAV complex temporarily dissociates in response to oxidative stress to attenuate ADP-ATP exchange, a rate-limiting step for OXPHOS. Finally, we found that dissociation of the DAV complex is mediated by phosphorylation of DNA-PKcs at its Thr2609 cluster by ATM kinase. Based on these findings, we propose that the coordination between the DAV complex and ATM serves as a novel oxidative stress checkpoint to decrease ROS production from mitochondrial OXPHOS and to hasten cellular recovery from OS.

Analyzing super-enhancer temporal dynamics reveals potential critical enhancers and their gene regulatory networks underlying skeletal muscle development.

Super-enhancers (SEs) govern the expression of genes defining cell identity. However, the dynamic landscape of SEs and their critical constituent enhancers involved in skeletal muscle development remains unclear. In this study, using pig as a model, we employed CUT&Tag to profile the enhancer-associated histone modification marker H3K27ac in skeletal muscle across two prenatal and three postnatal stages and investigated how SEs influence skeletal muscle development. We identified three SE families with distinct temporal dynamics: continuous (Con, 397), transient (TS, 434), and de novo (DN, 756). These SE families are associated with different temporal gene expression trajectories, biological functions, and DNA methylation levels. Notably, several lines of evidence suggest a potential prominent role of Con SEs in regulating porcine muscle development and meat traits. To pinpoint key cis-regulatory units in Con SEs, we developed an integrative approach that leverages information from eRNA annotation, GWAS signals and high-throughput capture STARR-seq experiments. Within Con SEs, we identified 20 candidate critical enhancers with meat and carcass-associated DNA variations that affect enhancer activity and inferred their upstream TFs and downstream target genes. As a proof of concept, we experimentally validated the role of one such enhancer and its potential target gene during myogenesis. Our findings reveal the dynamic regulatory features of SEs in skeletal muscle development and provide a general integrative framework for identifying critical enhancers underlying the formation of complex traits.

A genome-scale deep learning model to predict gene expression changes of genetic perturbations from multiplex biological networks.

Systematic characterization of biological effects to genetic perturbation is essential to the application of molecular biology and biomedicine. However, the experimental exhaustion of genetic perturbations on the genome-wide scale is challenging. Here, we show TranscriptionNet, a deep learning model that integrates multiple biological networks to systematically predict transcriptional profiles to three types of genetic perturbations based on transcriptional profiles induced by genetic perturbations in the L1000 project: RNA interference, clustered regularly interspaced short palindromic repeat, and overexpression. TranscriptionNet performs better than existing approaches in predicting inducible gene expression changes for all three types of genetic perturbations. TranscriptionNet can predict transcriptional profiles for all genes in existing biological networks and increases perturbational gene expression changes for each type of genetic perturbation from a few thousand to 26 945 genes. TranscriptionNet demonstrates strong generalization ability when comparing predicted and true gene expression changes on different external tasks. Overall, TranscriptionNet can systemically predict transcriptional consequences induced by perturbing genes on a genome-wide scale and thus holds promise to systemically detect gene function and enhance drug development and target discovery.

A novel interpretable deep learning-based computational framework designed synthetic enhancers with broad cross-species activity.

Enhancers play a critical role in dynamically regulating spatial-temporal gene expression and establishing cell identity, underscoring the significance of designing them with specific properties for applications in biosynthetic engineering and gene therapy. Despite numerous high-throughput methods facilitating genome-wide enhancer identification, deciphering the sequence determinants of their activity remains challenging. Here, we present the DREAM (DNA cis-Regulatory Elements with controllable Activity design platforM) framework, a novel deep learning-based approach for synthetic enhancer design. Proficient in uncovering subtle and intricate patterns within extensive enhancer screening data, DREAM achieves cutting-edge sequence-based enhancer activity prediction and highlights critical sequence features implicating strong enhancer activity. Leveraging DREAM, we have engineered enhancers that surpass the potency of the strongest enhancer within the Drosophila genome by approximately 3.6-fold. Remarkably, these synthetic enhancers exhibited conserved functionality across species that have diverged more than billion years, indicating that DREAM was able to learn highly conserved enhancer regulatory grammar. Additionally, we designed silencers and cell line-specific enhancers using DREAM, demonstrating its versatility. Overall, our study not only introduces an interpretable approach for enhancer design but also lays out a general framework applicable to the design of other types of cis-regulatory elements.

Generative Adversarial Matrix Completion Network based on Multi-Source Data Fusion for miRNA-Disease Associations Prediction.

Numerous biological studies have shown that considering disease-associated micro RNAs (miRNAs) as potential biomarkers or therapeutic targets offers new avenues for the diagnosis of complex diseases. Computational methods have gradually been introduced to reveal disease-related miRNAs. Considering that previous models have not fused sufficiently diverse similarities, that their inappropriate fusion methods may lead to poor quality of the comprehensive similarity network and that their results are often limited by insufficiently known associations, we propose a computational model called Generative Adversarial Matrix Completion Network based on Multi-source Data Fusion (GAMCNMDF) for miRNA-disease association prediction. We create a diverse network connecting miRNAs and diseases, which is then represented using a matrix. The main task of GAMCNMDF is to complete the matrix and obtain the predicted results. The main innovations of GAMCNMDF are reflected in two aspects: GAMCNMDF integrates diverse data sources and employs a nonlinear fusion approach to update the similarity networks of miRNAs and diseases. Also, some additional information is provided to GAMCNMDF in the form of a 'hint' so that GAMCNMDF can work successfully even when complete data are not available. Compared with other methods, the outcomes of 10-fold cross-validation on two distinct databases validate the superior performance of GAMCNMDF with statistically significant results. It is worth mentioning that we apply GAMCNMDF in the identification of underlying small molecule-related miRNAs, yielding outstanding performance results in this specific domain. In addition, two case studies about two important neoplasms show that GAMCNMDF is a promising prediction method.

Predicting potential small molecule-miRNA associations utilizing truncated schatten p-norm.

MicroRNAs (miRNAs) have significant implications in diverse human diseases and have proven to be effectively targeted by small molecules (SMs) for therapeutic interventions. However, current SM-miRNA association prediction models do not adequately capture SM/miRNA similarity. Matrix completion is an effective method for association prediction, but existing models use nuclear norm instead of rank function, which has some drawbacks. Therefore, we proposed a new approach for predicting SM-miRNA associations by utilizing the truncated schatten p-norm (TSPN). First, the SM/miRNA similarity was preprocessed by incorporating the Gaussian interaction profile kernel similarity method. This identified more SM/miRNA similarities and significantly improved the SM-miRNA prediction accuracy. Next, we constructed a heterogeneous SM-miRNA network by combining biological information from three matrices and represented the network with its adjacency matrix. Finally, we constructed the prediction model by minimizing the truncated schatten p-norm of this adjacency matrix and we developed an efficient iterative algorithmic framework to solve the model. In this framework, we also used a weighted singular value shrinkage algorithm to avoid the problem of excessive singular value shrinkage. The truncated schatten p-norm approximates the rank function more closely than the nuclear norm, so the predictions are more accurate. We performed four different cross-validation experiments on two separate datasets, and TSPN outperformed various most advanced methods. In addition, public literature confirms a large number of predictive associations of TSPN in four case studies. Therefore, TSPN is a reliable model for SM-miRNA association prediction.

Enhanced protein-metabolite correlation analysis: To investigate the association between Staphylococcus aureus mastitis and metabolic immune pathways.

Mastitis is a disease characterized by congestion, swelling, and inflammation of the mammary gland and usually caused by infection with pathogenic microorganisms. Furthermore, the development of mastitis is closely linked to the exogenous pathway of the gastrointestinal tract. However, the regulatory mechanisms governing the gut-metabolism-mammary axis remain incompletely understood. The present study revealed alterations in the gut microbiota of mastitis rats characterized by an increased abundance of the Proteobacteria phylum. Plasma analysis revealed significantly higher levels of L-isoleucine and cholic acid along with 7-ketodeoxycholic acid. Mammary tissue showed elevated levels of arachidonic acid metabolites and norlithocholic acid. Proteomic analysis showed increased levels of IFIH1, Tnfaip8l2, IRGM, and IRF5 in mastitis rats, which suggests that mastitis triggers an inflammatory response and immune stress. Follistatin (Fst) and progesterone receptor (Pgr) were significantly downregulated, raising the risk of breast cancer. Extracellular matrix (ECM) receptors and focal adhesion signaling pathways were downregulated, while blood-milk barrier integrity was disrupted. Analysis of protein-metabolic network regulation revealed that necroptosis, protein digestion and absorption, and arachidonic acid metabolism were the principal regulatory pathways involved in the development of mastitis. In short, the onset of mastitis leads to changes in the microbiota and alterations in the metabolic profiles of various biological samples, including colonic contents, plasma, and mammary tissue. Key manifestations include disturbances in bile acid metabolism, amino acid metabolism, and arachidonic acid metabolism. At the same time, the integrity of the blood-milk barrier is compromised while inflammation is promoted, thereby reducing cell adhesion in the mammary glands. These findings contribute to a more comprehensive understanding of the metabolic status of mastitis and provide new insights into its impact on the immune system.

mPFC DUSP1 mediates adolescent cocaine exposure-induced higher sensitivity to drug in adulthood.

Adolescent cocaine abuse increases the risk for developing addiction in later life, but the underlying molecular mechanism remains poorly understood. Here, we establish adolescent cocaine-exposed (ACE) male mouse models. A subthreshold dose of cocaine (sdC) treatment, insufficient to produce conditioned place preference (CPP) in adolescent mice, induces CPP in ACE mice during adulthood, along with more activated CaMKII-positive neurons, higher dual specificity protein kinase phosphatase-1 (Dusp1) mRNA, lower DUSP1 activity, and lower DUSP1 expression in CaMKII-positive neurons in the medial prefrontal cortex (mPFC). Overexpressing DUSP1 in CaMKII-positive neurons suppresses neuron activity and blocks sdC-induced CPP in ACE mice during adulthood. On the contrary, depleting DUSP1 in CaMKII-positive neurons activates more neurons and further enhances sdC-induced behavior in ACE mice during adulthood. Also, ERK1/2 might be a downstream signal of DUSP1 in the process. Our findings reveal a role of mPFC DUSP1 in ACE-induced higher sensitivity to the drug in adult mice. DUSP1 might be a potential pharmacological target to predict or treat the susceptibility to addictive drugs caused by adolescent substance use.

Regulation of YAP Promotor Accessibility in Endothelial Mechanotransduction.

BACKGROUND: Endothelial cells are constantly exposed to mechanical forces in the form of fluid shear stress, extracellular stiffness, and cyclic strain. The mechanoresponsive activity of YAP (Yes-associated protein) and its role in vascular development are well described; however, whether changes to transcription or epigenetic regulation of YAP are involved in these processes remains unanswered. Furthermore, how mechanical forces are transduced to the nucleus to drive transcriptional reprogramming in endothelial cells is poorly understood. The YAP target gene, AmotL2 (angiomotin-like 2), is a junctional mechanotransducer that connects cell-cell junctions to the nuclear membrane via the actin cytoskeleton. METHODS: We applied mechanical manipulations including shear flow, stretching, and substrate stiffness to endothelial cells to investigate the role of mechanical forces in modulating YAP transcription. Using in vitro and in vivo endothelial depletion of AmotL2, we assess nuclear morphology, chromatin organization (using transposase-accessible chromatin sequencing), and whole-mount immunofluorescent staining of the aorta to determine the regulation and functionality of YAP. Finally, we use genetic and chemical inhibition to uncouple the nuclear-cytoskeletal connection to investigate the role of this pathway on YAP transcription. RESULTS: Our results reveal that mechanical forces sensed at cell-cell junctions by the YAP target gene AmotL2 are directly involved in changes in global chromatin accessibility and activity of the histone methyltransferase EZH2, leading to modulation of YAP promotor activity. Functionally, shear stress-induced proliferation of endothelial cells in vivo was reliant on AmotL2 and YAP/TAZ (transcriptional coactivator with PDZ-binding motif) expression. Mechanistically, uncoupling of the nuclear-cytoskeletal connection from junctions and focal adhesions led to altered nuclear morphology, chromatin accessibility, and YAP promotor activity. CONCLUSIONS: Our findings reveal a role for AmotL2 and nuclear-cytoskeletal force transmission in modulating the epigenetic and transcriptional regulation of YAP to maintain a mechano-enforced positive feedback loop of vascular homeostasis. These findings may offer an explanation as to the proinflammatory phenotype that leads to aneurysm formation observed in AmotL2 endothelial deletion models.

Pathological polarizations from microglia to astrocyte contributes to spatial memory deficit in methamphetamine abstinence mice.

Previously, we found that dCA1 A1-like polarization of astrocytes contributes a lot to the spatial memory deficit in methamphetamine abstinence mice. However, the underlying mechanism remains unclear, resulting in a lack of promising therapeutic targets. Here, we found that methamphetamine abstinence mice exhibited an increased M1-like microglia and A1-like astrocytes, together with elevated levels of interleukin 1α and tumor necrosis factor α in dCA1. In vitro, the M1-like BV2 microglia cell medium, containing high levels of Interleukin 1α and tumor necrosis factor α, elevated A1-like polarization of astrocytes, which weakened their capacity for glutamate clearance. Locally suppressing dCA1 M1-like microglia activation with minocycline administration attenuated A1-like polarization of astrocytes, ameliorated dCA1 neurotoxicity, and, most importantly, rescued spatial memory in methamphetamine abstinence mice. The effective time window of minocycline treatment on spatial memory is the methamphetamine exposure period, rather than the long-term methamphetamine abstinence.

Ginkgo seed shell provides a unique model for bioinspired design.

Natural structural materials typically feature complex hierarchical anisotropic architectures, resulting in excellent damage tolerance. Such highly anisotropic structures, however, also provide an easy path for crack propagation, often leading to catastrophic fracture as evidenced, for example, by wood splitting. Here, we describe the weakly anisotropic structure of Ginkgo biloba (ginkgo) seed shell, which has excellent crack resistance in different directions. Ginkgo seed shell is composed of tightly packed polygonal sclereids with cell walls in which the cellulose microfibrils are oriented in a helicoidal pattern. We found that the sclereids contain distinct pits, special fine tubes like a "screw fastener," that interlock the helicoidal cell walls together. As a result, ginkgo seed shell demonstrates crack resistance in all directions, exhibiting specific fracture toughness that can rival other highly anisotropic natural materials, such as wood, bone, insect cuticle, and nacre. In situ characterization reveals ginkgo's unique toughening mechanism: pit-guided crack propagation. This mechanism forces the crack to depart from the weak compound middle lamella and enter into the sclereid, where the helicoidal cell wall significantly inhibits crack growth by the cleavage and breakage of the fibril-based cell walls. Ginkgo's toughening mechanism could provide guidelines for a new bioinspired strategy for the design of high-performance bulk materials.

Let-7 underlies metformin-induced inhibition of hepatic glucose production.

SignificanceA clear mechanistic understanding of metformin's antidiabetic effects is lacking. This is because suprapharmacological concentrations of metformin have been used in most studies. Using mouse models and human primary hepatocytes, we show that metformin, at clinically relevant doses, suppresses hepatic glucose production by activating a conserved regulatory pathway encompassing let-7, TET3, and a fetal isoform of hepatocyte nuclear factor 4 alpha (HNF4α). We demonstrate that metformin no longer has potent antidiabetic actions in a liver-specific let-7 loss-of-function mouse model and that hepatic delivery of let-7 ameliorates hyperglycemia and improves glucose homeostasis. Our results thus reveal an important role of the hepatic let-7/TET3/HNF4α axis in mediating the therapeutic effects of metformin and suggest that targeting this axis may be a potential therapeutic for diabetes.