Mutations in E3 ubiquitin ligase Parkin have been linked to familial Parkinson's disease. Accumulating evidence suggests that Parkin is a tumor suppressor, but the underlying mechanism is poorly understood. Here we show that Parkin is an E3 ubiquitin ligase for hypoxia-inducible factor 1α (HIF-1α). Parkin interacts with HIF-1α and promotes HIF-1α degradation through ubiquitination, which in turn inhibits metastasis of breast cancer cells. Parkin downregulation in breast cancer cells promotes metastasis, which can be inhibited by targeting HIF-1α with RNA interference or the small-molecule inhibitor YC-1. We further identify lysine 477 (K477) of HIF-1α as a major ubiquitination site for Parkin. K477R HIF-1α mutation and specific cancer-associated Parkin mutations largely abolish the functions of Parkin to ubiquitinate HIF-1α and inhibit cancer metastasis. Importantly, Parkin expression is inversely correlated with HIF-1α expression and metastasis in breast cancer. Our results reveal an important mechanism for Parkin in tumor suppression and HIF-1α regulation.

Axonemal dyneins are tethered to doublet microtubules inside cilia to drive ciliary beating, a process critical for cellular motility and extracellular fluid flow. Axonemal dyneins are evolutionarily and biochemically distinct from cytoplasmic dyneins that transport cargo, and the mechanisms regulating their localization and function are poorly understood. Here, we report a single-particle cryo-EM reconstruction of a three-headed axonemal dynein natively bound to doublet microtubules isolated from cilia. The slanted conformation of the axonemal dynein causes interaction of its motor domains with the neighboring dynein complex. Our structure shows how a heterotrimeric docking complex specifically localizes the linear array of axonemal dyneins to the doublet microtubule by directly interacting with the heavy chains. Our structural analysis establishes the arrangement of conserved heavy, intermediate and light chain subunits, and provides a framework to understand the roles of individual subunits and the interactions between dyneins during ciliary waveform generation.

Silymarin has been used for improving hepatic damage and lipid disorders, but its action mechanism remains to be clarified. Here, we investigate the contributions of the gut microbiota to the improvement of liver lipid metabolism by silymarin. We find i) strong and significant microbial shifts upon silymarin but not silibinin treatment; ii) over 60% variations of liver fat are explained by silymarin-induced bacterial B12 production in male rats but not in male germ-free mice; iii) fecal microbiota transplantation confirms their protective roles against liver fat accumulation; iv) upregulation of one-carbon metabolism and fatty acid degradation pathways are observed based on the liver transcriptome analyses; and v) in humans the delta changes of serum B12 associate negatively with the fluctuations of serum triglycerides. Overall, we reveal a mechanism of action underpinning the lipid-lowering effect of silymarin via the gut microbiota and its vitamin B12 producing capabilities.

Araliaceae species produce various classes of triterpene and triterpenoid saponins, such as the oleanane-type triterpenoids in Aralia species and dammarane-type saponins in Panax, valued for their medicinal properties. The lack of genome sequences of Panax relatives has hindered mechanistic insight into the divergence of triterpene saponins in Araliaceae. Here, we report a chromosome-level genome of Aralia elata with a total length of 1.05 Gb. The loss of 12 exons in the dammarenediol synthase (DDS)-encoding gene in A. elata after divergence from Panax might have caused the lack of dammarane-type saponin production, and a complementation assay shows that overexpression of the PgDDS gene from Panax ginseng in callus of A. elata recovers the accumulation of dammarane-type saponins. Tandem duplication events of triterpene biosynthetic genes are common in the A. elata genome, especially for AeCYP72As, AeCSLMs, and AeUGT73s, which function as tailoring enzymes of oleanane-type saponins and aralosides. More than 13 aralosides are de novo synthesized in Saccharomyces cerevisiae by overexpression of these genes in combination. This study sheds light on the diversity of saponins biosynthetic pathway in Araliaceae and will facilitate heterologous bioproduction of aralosides.

Pluripotent stem cells hold great promise in regenerative medicine and developmental biology studies. Mitochondrial metabolites, including tricarboxylic acid (TCA) cycle intermediates, have been reported to play critical roles in pluripotency. Here we show that TCA cycle enzymes including Pdha1, Pcb, Aco2, Cs, Idh3a, Ogdh, Sdha and Mdh2 are translocated to the nucleus during somatic cell reprogramming, primed-to-naive transition and totipotency acquisition. The nuclear-localized TCA cycle enzymes Pdha1, Pcb, Aco2, Cs, Idh3a promote somatic cell reprogramming and primed-to-naive transition. In addition, nuclear-localized TCA cycle enzymes, particularly nuclear-targeted Pdha1, facilitate the 2-cell program in pluripotent stem cells. Mechanistically, nuclear Pdha1 increases the acetyl-CoA and metabolite pool in the nucleus, leading to chromatin remodeling at pluripotency genes by enhancing histone H3 acetylation. Our results reveal an important role of mitochondrial TCA cycle enzymes in the epigenetic regulation of pluripotency that constitutes a mitochondria-to-nucleus retrograde signaling mode in different states of pluripotent acquisition.

The dual functions of TMEM16F as Ca2+-activated ion channel and lipid scramblase raise intriguing questions regarding their molecular basis. Intrigued by the ability of the FDA-approved drug niclosamide to inhibit TMEM16F-dependent syncytia formation induced by SARS-CoV-2, we examined cryo-EM structures of TMEM16F with or without bound niclosamide or 1PBC, a known blocker of TMEM16A Ca2+-activated Cl- channel. Here, we report evidence for a lipid scrambling pathway along a groove harboring a lipid trail outside the ion permeation pore. This groove contains the binding pocket for niclosamide and 1PBC. Mutations of two residues in this groove specifically affect lipid scrambling. Whereas mutations of some residues in the binding pocket of niclosamide and 1PBC reduce their inhibition of TMEM16F-mediated Ca2+ influx and PS exposure, other mutations preferentially affect the ability of niclosamide and/or 1PBC to inhibit TMEM16F-mediated PS exposure, providing further support for separate pathways for ion permeation and lipid scrambling.

Pancreatic ductal adenocarcinoma (PDAC) is a highly metastatic disease refractory to all targeted and immune therapies. However, our understanding of PDAC microenvironment especially the metastatic microenvironment is very limited partly due to the inaccessibility to metastatic tumor tissues. Here, we present the single-cell transcriptomic landscape of synchronously resected PDAC primary tumors and matched liver metastases. We perform comparative analysis on both cellular composition and functional phenotype between primary and metastatic tumors. Tumor cells exhibit distinct transcriptomic profile in liver metastasis with clearly defined evolutionary routes from cancer cells in primary tumor. We also identify specific subtypes of stromal and immune cells critical to the formation of the pro-tumor microenvironment in metastatic lesions, including RGS5+ cancer-associated fibroblasts, CCL18+ lipid-associated macrophages, S100A8+ neutrophils and FOXP3+ regulatory T cells. Cellular interactome analysis further reveals that the lack of tumor-immune cell interaction in metastatic tissues contributes to the formation of the immunosuppressive microenvironment. Our study provides a comprehensive characterization of the transcriptional landscape of PDAC liver metastasis.

Patients affected by neurofibromatosis type 1 (NF1) frequently show muscle weakness with unknown etiology. Here we show that, in mice, Neurofibromin 1 (Nf1) is not required in muscle fibers, but specifically in early postnatal myogenic progenitors (MPs), where Nf1 loss led to cell cycle exit and differentiation blockade, depleting the MP pool resulting in reduced myonuclear accretion as well as reduced muscle stem cell numbers. This was caused by precocious induction of stem cell quiescence coupled to metabolic reprogramming of MPs impinging on glycolytic shutdown, which was conserved in muscle fibers. We show that a Mek/Erk/NOS pathway hypersensitizes Nf1-deficient MPs to Notch signaling, consequently, early postnatal Notch pathway inhibition ameliorated premature quiescence, metabolic reprogramming and muscle growth. This reveals an unexpected role of Ras/Mek/Erk signaling supporting postnatal MP quiescence in concert with Notch signaling, which is controlled by Nf1 safeguarding coordinated muscle growth and muscle stem cell pool establishment. Furthermore, our data suggest transmission of metabolic reprogramming across cellular differentiation, affecting fiber metabolism and function in NF1.

Understanding the molecular mechanisms regulating the maintenance and destruction of intervertebral disc may lead to the development of new therapies for intervertebral disc degeneration (IDD). Here we present evidence from miRNA microarray analyses of clinical data sets along with in vitro and in vivo experiments that miR-141 is a key regulator of IDD. Gain- and loss-of-function studies show that miR-141 drives IDD by inducing nucleus pulposus (NP) apoptosis. Furthermore, miR-141 KO in mice attenuated spontaneous and surgically induced IDD. Mechanistically, miR-141 promotes IDD development by targeting and depleting SIRT1, a negative regulator of NF-κB pathway. Therapeutically, upregulation or downregulation of miR-141 by nanoparticle delivery in IDD model aggravated or alleviated experimental IDD, respectively. Our findings reveal a novel mechanism by which miR-141, in part, promotes IDD progression by interacting with SIRT1/NF-κB pathway. Blockade of miR-141 in vivo may serve as a potential therapeutic approach in the treatment of IDD.

Protein-protein interactions featuring intricate binding epitopes remain challenging targets for synthetic inhibitors. Interactions of NEMO, a scaffolding protein central to NF-κB signaling, exemplify this challenge. Various regulators are known to interact with different coiled coil regions of NEMO, but the topological complexity of this protein has limited inhibitor design. We undertook a comprehensive effort to block the interaction between vFLIP, a Kaposi's sarcoma herpesviral oncoprotein, and NEMO using small molecule screening and rational design. Our efforts reveal that a tertiary protein structure mimic of NEMO is necessary for potent inhibition. The rationally designed mimic engages vFLIP directly causing complex disruption, protein degradation and suppression of NF-κB signaling in primary effusion lymphoma (PEL). NEMO mimic treatment induces cell death and delays tumor growth in a PEL xenograft model. Our studies with this inhibitor reveal the critical nexus of signaling complex stability in the regulation of NF-κB by a viral oncoprotein.

LIN28A is important in somatic reprogramming and pluripotency regulation. Although previous studies addressed that LIN28A can repress let-7 microRNA maturation in the cytoplasm, few focused on its role within the nucleus. Here, we show that the nucleolus-localized LIN28A protein undergoes liquid-liquid phase separation (LLPS) in mouse embryonic stem cells (mESCs) and in vitro. The RNA binding domains (RBD) and intrinsically disordered regions (IDR) of LIN28A contribute to LIN28A and the other nucleolar proteins' phase-separated condensate establishment. S120A, S200A and R192G mutations in the IDR result in subcellular mislocalization of LIN28A and abnormal nucleolar phase separation. Moreover, we find that the naive-to-primed pluripotency state conversion and the reprogramming are associated with dynamic nucleolar remodeling, which depends on LIN28A's phase separation capacity, because the LIN28A IDR point mutations abolish its role in regulating nucleolus and in these cell fate decision processes, and an exogenous IDR rescues it. These findings shed light on the nucleolar function in pluripotent stem cell states and on a non-canonical RNA-independent role of LIN28A in phase separation and cell fate decisions.

Date palm (Phoenix dactylifera L.) is a cultivated woody plant species with agricultural and economic importance. Here we report a genome assembly for an elite variety (Khalas), which is 605.4 Mb in size and covers >90% of the genome (~671 Mb) and >96% of its genes (~41,660 genes). Genomic sequence analysis demonstrates that P. dactylifera experienced a clear genome-wide duplication after either ancient whole genome duplications or massive segmental duplications. Genetic diversity analysis indicates that its stress resistance and sugar metabolism-related genes tend to be enriched in the chromosomal regions where the density of single-nucleotide polymorphisms is relatively low. Using transcriptomic data, we also illustrate the date palm's unique sugar metabolism that underlies fruit development and ripening. Our large-scale genomic and transcriptomic data pave the way for further genomic studies not only on P. dactylifera but also other Arecaceae plants.

Access to experimental X-ray diffraction image data is fundamental for validation and reproduction of macromolecular models and indispensable for development of structural biology processing methods. Here, we established a diffraction data publication and dissemination system, Structural Biology Data Grid (SBDG; data.sbgrid.org), to preserve primary experimental data sets that support scientific publications. Data sets are accessible to researchers through a community driven data grid, which facilitates global data access. Our analysis of a pilot collection of crystallographic data sets demonstrates that the information archived by SBDG is sufficient to reprocess data to statistics that meet or exceed the quality of the original published structures. SBDG has extended its services to the entire community and is used to develop support for other types of biomedical data sets. It is anticipated that access to the experimental data sets will enhance the paradigm shift in the community towards a much more dynamic body of continuously improving data analysis.

Understanding and control of structures and rates involved in protein ligand binding are essential for drug design. Unfortunately, atomistic molecular dynamics (MD) simulations cannot directly sample the excessively long residence and rearrangement times of tightly binding complexes. Here we exploit the recently developed multi-ensemble Markov model framework to compute full protein-peptide kinetics of the oncoprotein fragment 25-109Mdm2 and the nano-molar inhibitor peptide PMI. Using this system, we report, for the first time, direct estimates of kinetics beyond the seconds timescale using simulations of an all-atom MD model, with high accuracy and precision. These results only require explicit simulations on the sub-milliseconds timescale and are tested against existing mutagenesis data and our own experimental measurements of the dissociation and association rates. The full kinetic model reveals an overall downhill but rugged binding funnel with multiple pathways. The overall strong binding arises from a variety of conformations with different hydrophobic contact surfaces that interconvert on the milliseconds timescale.

NLRP1 and CARD8 are related cytosolic sensors that upon activation form supramolecular signalling complexes known as canonical inflammasomes, resulting in caspase-1 activation, cytokine maturation and/or pyroptotic cell death. NLRP1 and CARD8 use their C-terminal (CT) fragments containing a caspase recruitment domain (CARD) and the UPA (conserved in UNC5, PIDD, and ankyrins) subdomain for self-oligomerization, which in turn form the platform to recruit the inflammasome adaptor ASC (apoptosis-associated speck-like protein containing a CARD) or caspase-1, respectively. Here, we report cryo-EM structures of NLRP1-CT and CARD8-CT assemblies, in which the respective CARDs form central helical filaments that are promoted by oligomerized, but flexibly linked, UPAs surrounding the filaments. Through biochemical and cellular approaches, we demonstrate that the UPA itself reduces the threshold needed for NLRP1-CT and CARD8-CT filament formation and signalling. Structural analyses provide insights on the mode of ASC recruitment by NLRP1-CT and the contrasting direct recruitment of caspase-1 by CARD8-CT. We also discover that subunits in the central NLRP1CARD filament dimerize with additional exterior CARDs, which roughly doubles its thickness and is unique among all known CARD filaments. Finally, we engineer and determine the structure of an ASCCARD-caspase-1CARD octamer, which suggests that ASC uses opposing surfaces for NLRP1, versus caspase-1, recruitment. Together these structures capture the architecture and specificity of the active NLRP1 and CARD8 inflammasomes in addition to key heteromeric CARD-CARD interactions governing inflammasome signalling.

The epsin family of endocytic adapter proteins are widely expressed, and interact with both proteins and lipids to regulate a variety of cell functions. However, the role of epsins in atherosclerosis is poorly understood. Here, we show that deletion of endothelial epsin proteins reduces inflammation and attenuates atherosclerosis using both cell culture and mouse models of this disease. In atherogenic cholesterol-treated murine aortic endothelial cells, epsins interact with the ubiquitinated endoplasmic reticulum protein inositol 1,4,5-trisphosphate receptor type 1 (IP3R1), which triggers proteasomal degradation of this calcium release channel. Epsins potentiate its degradation via this interaction. Genetic reduction of endothelial IP3R1 accelerates atherosclerosis, whereas deletion of endothelial epsins stabilizes IP3R1 and mitigates inflammation. Reduction of IP3R1 in epsin-deficient mice restores atherosclerotic progression. Taken together, epsin-mediated degradation of IP3R1 represents a previously undiscovered biological role for epsin proteins and may provide new therapeutic targets for the treatment of atherosclerosis and other diseases.

Developing injectable antibacterial and conductive shape memory hemostatic with high blood absorption and fast recovery for irregularly shaped and noncompressible hemorrhage remains a challenge. Here we report injectable antibacterial conductive cryogels based on carbon nanotube (CNT) and glycidyl methacrylate functionalized quaternized chitosan for lethal noncompressible hemorrhage hemostasis and wound healing. These cryogels present robust mechanical strength, rapid blood-triggered shape recovery and absorption speed, and high blood uptake capacity. Moreover, cryogels show better blood-clotting ability, higher blood cell and platelet adhesion and activation than gelatin sponge and gauze. Cryogel with 4 mg/mL CNT (QCSG/CNT4) shows better hemostatic capability than gauze and gelatin hemostatic sponge in mouse-liver injury model and mouse-tail amputation model, and better wound healing performance than Tegaderm™ film. Importantly, QCSG/CNT4 presents excellent hemostatic performance in rabbit liver defect lethal noncompressible hemorrhage model and even better hemostatic ability than Combat Gauze in standardized circular liver bleeding model.

Computational cell type identification is a fundamental step in single-cell omics data analysis. Supervised celltyping methods have gained increasing popularity in single-cell RNA-seq data because of the superior performance and the availability of high-quality reference datasets. Recent technological advances in profiling chromatin accessibility at single-cell resolution (scATAC-seq) have brought new insights to the understanding of epigenetic heterogeneity. With continuous accumulation of scATAC-seq datasets, supervised celltyping method specifically designed for scATAC-seq is in urgent need. Here we develop Cellcano, a computational method based on a two-round supervised learning algorithm to identify cell types from scATAC-seq data. The method alleviates the distributional shift between reference and target data and improves the prediction performance. After systematically benchmarking Cellcano on 50 well-designed celltyping tasks from various datasets, we show that Cellcano is accurate, robust, and computationally efficient. Cellcano is well-documented and freely available at https://marvinquiet.github.io/Cellcano/ .

Stereocilia are actin-based cell protrusions of inner ear hair cells and are indispensable for mechanotransduction. Ankle links connect the ankle region of developing stereocilia, playing an essential role in stereocilia development. WHRN, PDZD7, ADGRV1 and USH2A have been identified to form the so-called ankle link complex (ALC); however, the detailed mechanism underlying the temporal emergence and degeneration of ankle links remains elusive. Here we show that WHRN and PDZD7 orchestrate ADGRV1 and USH2A to assemble the ALC through liquid-liquid phase separation (LLPS). Disruption of the ALC multivalency for LLPS largely abolishes the distribution of WHRN at the ankle region of stereocilia. Interestingly, high concentration of ADGRV1 inhibits LLPS, providing a potential mechanism for ALC disassembly. Moreover, certain deafness mutations of ALC genes weaken the multivalent interactions of ALC and impair LLPS. In conclusion, our study demonstrates that LLPS mediates ALC formation, providing essential clues for understanding the pathogenesis of deafness.

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.