First-in-class MKK4 inhibitors enhance liver regeneration and prevent liver failure.

Diminished hepatocyte regeneration is a key feature of acute and chronic liver diseases and after extended liver resections, resulting in the inability to maintain or restore a sufficient functional liver mass. Therapies to restore hepatocyte regeneration are lacking, making liver transplantation the only curative option for end-stage liver disease. Here, we report on the structure-based development and characterization (nuclear magnetic resonance [NMR] spectroscopy) of first-in-class small molecule inhibitors of the dual-specificity kinase MKK4 (MKK4i). MKK4i increased liver regeneration upon hepatectomy in murine and porcine models, allowed for survival of pigs in a lethal 85% hepatectomy model, and showed antisteatotic and antifibrotic effects in liver disease mouse models. A first-in-human phase I trial (European Union Drug Regulating Authorities Clinical Trials [EudraCT] 2021-000193-28) with the clinical candidate HRX215 was conducted and revealed excellent safety and pharmacokinetics. Clinical trials to probe HRX215 for prevention/treatment of liver failure after extensive oncological liver resections or after transplantation of small grafts are warranted.

Tracking cell-type-specific temporal dynamics in human and mouse brains.

Progenitor cells are critical in preserving organismal homeostasis, yet their diversity and dynamics in the aged brain remain underexplored. We introduced TrackerSci, a single-cell genomic method that combines newborn cell labeling and combinatorial indexing to characterize the transcriptome and chromatin landscape of proliferating progenitor cells in vivo. Using TrackerSci, we investigated the dynamics of newborn cells in mouse brains across various ages and in a mouse model of Alzheimer's disease. Our dataset revealed diverse progenitor cell types in the brain and their epigenetic signatures. We further quantified aging-associated shifts in cell-type-specific proliferation and differentiation and deciphered the associated molecular programs. Extending our study to the progenitor cells in the aged human brain, we identified conserved genetic signatures across species and pinpointed region-specific cellular dynamics, such as the reduced oligodendrogenesis in the cerebellum. We anticipate that TrackerSci will be broadly applicable to unveil cell-type-specific temporal dynamics in diverse systems.

Systematic discovery of mutation-directed neo-protein-protein interactions in cancer.

Comprehensive sequencing of patient tumors reveals genomic mutations across tumor types that enable tumorigenesis and progression. A subset of oncogenic driver mutations results in neomorphic activity where the mutant protein mediates functions not engaged by the parental molecule. Here, we identify prevalent variant-enabled neomorph-protein-protein interactions (neoPPI) with a quantitative high-throughput differential screening (qHT-dS) platform. The coupling of highly sensitive BRET biosensors with miniaturized coexpression in an ultra-HTS format allows large-scale monitoring of the interactions of wild-type and mutant variant counterparts with a library of cancer-associated proteins in live cells. The screening of 17,792 interactions with 2,172,864 data points revealed a landscape of gain of interactions encompassing both oncogenic and tumor suppressor mutations. For example, the recurrent BRAF V600E lesion mediates KEAP1 neoPPI, rewiring a BRAFV600E/KEAP1 signaling axis and creating collateral vulnerability to NQO1 substrates, offering a combination therapeutic strategy. Thus, cancer genomic alterations can create neo-interactions, informing variant-directed therapeutic approaches for precision medicine.

Host-Specific Evolutionary and Transmission Dynamics Shape the Functional Diversification of Staphylococcus epidermidis in Human Skin.

Metagenomic inferences of bacterial strain diversity and infectious disease transmission studies largely assume a dominant, within-individual haplotype. We hypothesize that within-individual bacterial population diversity is critical for homeostasis of a healthy microbiome and infection risk. We characterized the evolutionary trajectory and functional distribution of Staphylococcus epidermidis-a keystone skin microbe and opportunistic pathogen. Analyzing 1,482 S. epidermidis genomes from 5 healthy individuals, we found that skin S. epidermidis isolates coalesce into multiple founder lineages rather than a single colonizer. Transmission events, natural selection, and pervasive horizontal gene transfer result in population admixture within skin sites and dissemination of antibiotic resistance genes within-individual. We provide experimental evidence for how admixture can modulate virulence and metabolism. Leveraging data on the contextual microbiome, we assess how interspecies interactions can shape genetic diversity and mobile gene elements. Our study provides insights into how within-individual evolution of human skin microbes shapes their functional diversification.

LECT2, a Ligand for Tie1, Plays a Crucial Role in Liver Fibrogenesis.

Liver fibrosis is a very common condition seen in millions of patients with various liver diseases, and yet no effective treatments are available owing to poorly characterized molecular pathogenesis. Here, we show that leukocyte cell-derived chemotaxin 2 (LECT2) is a functional ligand of Tie1, a poorly characterized endothelial cell (EC)-specific orphan receptor. Upon binding to Tie1, LECT2 interrupts Tie1/Tie2 heterodimerization, facilitates Tie2/Tie2 homodimerization, activates PPAR signaling, and inhibits the migration and tube formations of EC. In vivo studies showed that LECT2 overexpression inhibits portal angiogenesis, promotes sinusoid capillarization, and worsens fibrosis, whereas these changes were reversed in Lect2-KO mice. Adeno-associated viral vector serotype 9 (AAV9)-LECT2 small hairpin RNA (shRNA) treatment significantly attenuates fibrosis. Upregulation of LECT2 is associated with advanced human liver fibrosis staging. We concluded that targeting LECT2/Tie1 signaling may represent a potential therapeutic target for liver fibrosis, and serum LECT2 level may be a potential biomarker for the screening and diagnosis of liver fibrosis.

Therapeutic Ligands Antagonize Estrogen Receptor Function by Impairing Its Mobility.

Estrogen receptor-positive (ER+) breast cancers frequently remain dependent on ER signaling even after acquiring resistance to endocrine agents, prompting the development of optimized ER antagonists. Fulvestrant is unique among approved ER therapeutics due to its capacity for full ER antagonism, thought to be achieved through ER degradation. The clinical potential of fulvestrant is limited by poor physicochemical features, spurring attempts to generate ER degraders with improved drug-like properties. We show that optimization of ER degradation does not guarantee full ER antagonism in breast cancer cells; ER "degraders" exhibit a spectrum of transcriptional activities and anti-proliferative potential. Mechanistically, we find that fulvestrant-like antagonists suppress ER transcriptional activity not by ER elimination, but by markedly slowing the intra-nuclear mobility of ER. Increased ER turnover occurs as a consequence of ER immobilization. These findings provide proof-of-concept that small molecule perturbation of transcription factor mobility may enable therapeutic targeting of this challenging target class.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans.

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

FinnGen provides genetic insights from a well-phenotyped isolated population.

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

Thermal-expansion offset for high-performance fuel cell cathodes.

One challenge for the commercial development of solid oxide fuel cells as efficient energy-conversion devices is thermo-mechanical instability. Large internal-strain gradients caused by the mismatch in thermal expansion behaviour between different fuel cell components are the main cause of this instability, which can lead to cell degradation, delamination or fracture1-4. Here we demonstrate an approach to realizing full thermo-mechanical compatibility between the cathode and other cell components by introducing a thermal-expansion offset. We use reactive sintering to combine a cobalt-based perovskite with high electrochemical activity and large thermal-expansion coefficient with a negative-thermal-expansion material, thus forming a composite electrode with a thermal-expansion behaviour that is well matched to that of the electrolyte. A new interphase is formed because of the limited reaction between the two materials in the composite during the calcination process, which also creates A-site deficiencies in the perovskite. As a result, the composite shows both high activity and excellent stability. The introduction of reactive negative-thermal-expansion components may provide a general strategy for the development of fully compatible and highly active electrodes for solid oxide fuel cells.

The power of genetic diversity in genome-wide association studies of lipids.

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

SENP1-Sirt3 Signaling Controls Mitochondrial Protein Acetylation and Metabolism.

Sirt3, as a major mitochondrial nicotinamide adenine dinucleotide (NAD)-dependent deacetylase, is required for mitochondrial metabolic adaption to various stresses. However, how to regulate Sirt3 activity responding to metabolic stress remains largely unknown. Here, we report Sirt3 as a SUMOylated protein in mitochondria. SUMOylation suppresses Sirt3 catalytic activity. SUMOylation-deficient Sirt3 shows elevated deacetylation on mitochondrial proteins and increased fatty acid oxidation. During fasting, SUMO-specific protease SENP1 is accumulated in mitochondria and quickly de-SUMOylates and activates Sirt3. SENP1 deficiency results in hyper-SUMOylation of Sirt3 and hyper-acetylation of mitochondrial proteins, which reduces mitochondrial metabolic adaption responding to fasting. Furthermore, we find that fasting induces SENP1 translocation into mitochondria to activate Sirt3. The studies on mice show that Sirt3 SUMOylation mutation reduces fat mass and antagonizes high-fat diet (HFD)-induced obesity via increasing oxidative phosphorylation and energy expenditure. Our results reveal that SENP1-Sirt3 signaling modulates Sirt3 activation and mitochondrial metabolism during metabolic stress.

Boosted hydrogen evolution kinetics of heteroatom-doped carbons with isolated Zn as an accelerant.

Carbon-based single-atom catalysts, a promising candidate in electrocatalysis, offer insights into electron-donating effects of metal center on adjacent atoms. Herein, we present a practical strategy to rationally design a model catalyst with a single zinc (Zn) atom coordinated with nitrogen and sulfur atoms in a multilevel carbon matrix. The Zn site exhibits an atomic interface configuration of ZnN4S1, where Zn's electron injection effect enables thermal-neutral hydrogen adsorption on neighboring atoms, pushing the activity boundaries of carbon electrocatalysts toward electrochemical hydrogen evolution to an unprecedented level. Experimental and theoretical analyses confirm the low-barrier Volmer-Tafel mechanism of proton reduction, while the multishell hollow structures facilitate the hydrogen evolution even at high current intensities. This work provides insights for understanding the actual active species during hydrogen evolution reaction and paves the way for designing high-performance electrocatalysts.

Hidden hotspots of amphibian biodiversity in China.

Identifying and protecting hotspots of endemism and species richness is crucial for mitigating the global biodiversity crisis. However, our understanding of spatial diversity patterns is far from complete, which severely limits our ability to conserve biodiversity hotspots. Here, we report a comprehensive analysis of amphibian species diversity in China, one of the most species-rich countries on Earth. Our study combines 20 y of field surveys with new molecular analyses of 521 described species and also identifies 100 potential cryptic species. We identify 10 hotspots of amphibian diversity in China, each with exceptional species richness and endemism and with exceptional phylogenetic diversity and phylogenetic endemism (based on a new time-calibrated, species-level phylogeny for Chinese amphibians). These 10 hotspots encompass 59.6% of China's described amphibian species, 49.0% of cryptic species, and 55.6% of species endemic to China. Only four of these 10 hotspots correspond to previously recognized biodiversity hotspots. The six new hotspots include the Nanling Mountains and other mountain ranges in South China. Among the 186 species in the six new hotspots, only 9.7% are well covered by protected areas and most (88.2%) are exposed to high human impacts. Five of the six new hotspots are under very high human pressure and are in urgent need of protection. We also find that patterns of richness in cryptic species are significantly related to those in described species but are not identical.

A multi-ancestry genome-wide association study in type 1 diabetes.

Type 1 diabetes (T1D) is an autoimmune disease caused by destruction of the pancreatic ß-cells. Genome-wide association (GWAS) and fine mapping studies have been conducted mainly in European ancestry (EUR) populations. We performed a multi-ancestry GWAS to identify SNPs and HLA alleles associated with T1D risk and age at onset. EUR families (N = 3223), and unrelated individuals of African (AFR, N = 891) and admixed (Hispanic/Latino) ancestry (AMR, N = 308) were genotyped using the Illumina HumanCoreExome BeadArray, with imputation to the TOPMed reference panel. The Multi-Ethnic HLA reference panel was utilized to impute HLA alleles and amino acid residues. Logistic mixed models (T1D risk) and frailty models (age at onset) were used for analysis. In GWAS meta-analysis, seven loci were associated with T1D risk at genome-wide significance: PTPN22, HLA-DQA1, IL2RA, RNLS, INS, IKZF4-RPS26-ERBB3, and SH2B3, with four associated with T1D age at onset (PTPN22, HLA-DQB1, INS, and ERBB3). AFR and AMR meta-analysis revealed NRP1 as associated with T1D risk and age at onset, although NRP1 variants were not associated in EUR ancestry. In contrast, the PTPN22 variant was significantly associated with risk only in EUR ancestry. HLA alleles and haplotypes most significantly associated with T1D risk in AFR and AMR ancestry differed from that seen in EUR ancestry; in addition, the HLA-DRB1*08:02-DQA1*04:01-DQB1*04:02 haplotype was 'protective' in AMR while HLA-DRB1*08:01-DQA1*04:01-DQB1*04:02 haplotype was 'risk' in EUR ancestry, differing only at HLA-DRB1*08. These results suggest that much larger sample sizes in non-EUR populations are required to capture novel loci associated with T1D risk.

Scalable mixed model methods for set-based association studies on large-scale categorical data analysis and its application to exome-sequencing data in UK Biobank.

The ongoing release of large-scale sequencing data in the UK Biobank allows for the identification of associations between rare variants and complex traits. SAIGE-GENE+ is a valid approach to conducting set-based association tests for quantitative and binary traits. However, for ordinal categorical phenotypes, applying SAIGE-GENE+ with treating the trait as quantitative or binarizing the trait can cause inflated type I error rates or power loss. In this study, we propose a scalable and accurate method for rare-variant association tests, POLMM-GENE, in which we used a proportional odds logistic mixed model to characterize ordinal categorical phenotypes while adjusting for sample relatedness. POLMM-GENE fully utilizes the categorical nature of phenotypes and thus can well control type I error rates while remaining powerful. In the analyses of UK Biobank 450k whole-exome-sequencing data for five ordinal categorical traits, POLMM-GENE identified 54 gene-phenotype associations.

Structural basis for transcription activation by the nitrate-responsive regulator NarL.

Transcription activation is a crucial step of regulation during transcription initiation and a classic check point in response to different stimuli and stress factors. The Escherichia coli NarL is a nitrate-responsive global transcription factor that controls the expression of nearly 100 genes. However, the molecular mechanism of NarL-mediated transcription activation is not well defined. Here we present a cryo-EM structure of NarL-dependent transcription activation complex (TAC) assembled on the yeaR promoter at 3.2 Å resolution. Our structure shows that the NarL dimer binds at the -43.5 site of the promoter DNA with its C-terminal domain (CTD) not only binding to the DNA but also making interactions with RNA polymerase subunit alpha CTD (αCTD). The key role of these NarL-mediated interactions in transcription activation was further confirmed by in vivo and in vitro transcription assays. Additionally, the NarL dimer binds DNA in a different plane from that observed in the structure of class II TACs. Unlike the canonical class II activation mechanism, NarL does not interact with σ4, while RNAP αCTD is bound to DNA on the opposite side of NarL. Our findings provide a structural basis for detailed mechanistic understanding of NarL-dependent transcription activation on yeaR promoter and reveal a potentially novel mechanism of transcription activation.

CLCF1 signaling restrains thermogenesis and disrupts metabolic homeostasis by inhibiting mitochondrial biogenesis in brown adipocytes.

Great progress has been made in identifying positive regulators that activate adipocyte thermogenesis, but negative regulatory signaling of thermogenesis remains poorly understood. Here, we found that cardiotrophin-like cytokine factor 1 (CLCF1) signaling led to loss of brown fat identity, which impaired thermogenic capacity. CLCF1 levels decreased during thermogenic stimulation but were considerably increased in obesity. Adipocyte-specific CLCF1 transgenic (CLCF1-ATG) mice showed impaired energy expenditure and severe cold intolerance. Elevated CLCF1 triggered whitening of brown adipose tissue by suppressing mitochondrial biogenesis. Mechanistically, CLCF1 bound and activated ciliary neurotrophic factor receptor (CNTFR) and augmented signal transducer and activator of transcription 3 (STAT3) signaling. STAT3 transcriptionally inhibited both peroxisome proliferator-activated receptor-γ coactivator (PGC) 1α and 1ß, which thereafter restrained mitochondrial biogenesis in adipocytes. Inhibition of CNTFR or STAT3 could diminish the inhibitory effects of CLCF1 on mitochondrial biogenesis and thermogenesis. As a result, CLCF1-TG mice were predisposed to develop metabolic dysfunction even without external metabolic stress. Our findings revealed a brake signal on nonshivering thermogenesis and suggested that targeting this pathway could be used to restore brown fat activity and systemic metabolic homeostasis in obesity.

The mechano-chemical circuit drives skin organoid self-organization.

Stem cells in organoids self-organize into tissue patterns with unknown mechanisms. Here, we use skin organoids to analyze this process. Cell behavior videos show that the morphological transformation from multiple spheroidal units with morphogenesis competence (CMU) to planar skin is characterized by two abrupt cell motility-increasing events before calming down. The self-organizing processes are controlled by a morphogenetic module composed of molecular sensors, modulators, and executers. Increasing dermal stiffness provides the initial driving force (driver) which activates Yap1 (sensor) in epidermal cysts. Notch signaling (modulator 1) in epidermal cyst tunes the threshold of Yap1 activation. Activated Yap1 induces Wnts and MMPs (epidermal executers) in basal cells to facilitate cellular flows, allowing epidermal cells to protrude out from the CMU. Dermal cell-expressed Rock (dermal executer) generates a stiff force bridge between two CMU and accelerates tissue mixing via activating Laminin and ß1-integrin. Thus, this self-organizing coalescence process is controlled by a mechano-chemical circuit. Beyond skin, self-organization in organoids may use similar mechano-chemical circuit structures.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids.

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Aminopeptidase MNP-1 triggers intestine protease production by activating daf-16 nuclear location to degrade pore-forming toxins in Caenorhabditis elegans.

Pore-forming toxins (PFTs) are effective tools for pathogens infection. By disrupting epithelial barriers and killing immune cells, PFTs promotes the colonization and reproduction of pathogenic microorganisms in their host. In turn, the host triggers defense responses, such as endocytosis, exocytosis, or autophagy. Bacillus thuringiensis (Bt) bacteria produce PFT, known as crystal proteins (Cry) which damage the intestinal cells of insects or nematodes, eventually killing them. In insects, aminopeptidase N (APN) has been shown to act as an important receptor for Cry toxins. Here, using the nematode Caenorhabditis elegans as model, an extensive screening of APN gene family was performed to analyze the potential role of these proteins in the mode of action of Cry5Ba against the nematode. We found that one APN, MNP-1, participate in the toxin defense response, since the mnp-1(ok2434) mutant showed a Cry5Ba hypersensitive phenotype. Gene expression analysis in mnp-1(ok2434) mutant revealed the involvement of two protease genes, F19C6.4 and R03G8.6, that participate in Cry5Ba degradation. Finally, analysis of the transduction pathway involved in F19C6.4 and R03G8.6 expression revealed that upon Cry5Ba exposure, the worms up regulated both protease genes through the activation of the FOXO transcription factor DAF-16, which was translocated into the nucleus. The nuclear location of DAF-16 was found to be dependent on mnp-1 under Cry5Ba treatment. Our work provides evidence of new host responses against PFTs produced by an enteric pathogenic bacterium, resulting in activation of host intestinal proteases that degrade the PFT in the intestine.