miR‐9‐5p suppresses pro‐fibrogenic transformation of fibroblasts and prevents organ fibrosis by targeting NOX4 and TGFBR2

Uncontrolled extracellular matrix (ECM) production by fibroblasts in response to injury contributes to fibrotic diseases, including idiopathic pulmonary fibrosis (IPF). Reactive oxygen species (ROS) generation is involved in the pathogenesis of IPF. Transforming growth factor‐β1 (TGF‐β1) stimulates the production of NADPH oxidase 4 (NOX4)‐dependent ROS, promoting lung fibrosis (LF). Dysregulation of microRNAs (miRNAs) has been shown to contribute to LF. To identify miRNAs involved in redox regulation relevant for IPF, we performed arrays in human lung fibroblasts exposed to ROS. miR‐9‐5p was selected as the best candidate and we demonstrate its inhibitory effect on TGF‐β receptor type II (TGFBR2) and NOX4 expression. Increased expression of miR‐9‐5p abrogates TGF‐β1‐dependent myofibroblast phenotypic transformation. In the mouse model of bleomycin‐induced LF, miR‐9‐5p dramatically reduces fibrogenesis and inhibition of miR‐9‐5p and prevents its anti‐fibrotic effect both in vitro and in vivo. In lung specimens from patients with IPF, high levels of miR‐9‐5p are found. In omentum‐derived mesothelial cells (MCs) from patients subjected to peritoneal dialysis (PD), miR‐9‐5p also inhibits mesothelial to myofibroblast transformation. We propose that TGF‐β1 induces miR‐9‐5p expression as a self‐limiting homeostatic response.     

Hypoxia and loss of PHD2 inactivate stromal fibroblasts to decrease tumour stiffness and metastasis

Cancer‐associated fibroblasts (CAFs) interact with tumour cells and promote growth and metastasis. Here, we show that CAF activation is reversible: chronic hypoxia deactivates CAFs, resulting in the loss of contractile force, reduced remodelling of the surrounding extracellular matrix and, ultimately, impaired CAF‐mediated cancer cell invasion. Hypoxia inhibits prolyl hydroxylase domain protein 2 (PHD2), leading to hypoxia‐inducible factor (HIF)‐1α stabilisation, reduced expression of αSMA and periostin, and reduced myosin II activity. Loss of PHD2 in CAFs phenocopies the effects of hypoxia, which can be prevented by simultaneous depletion of HIF‐1α. Treatment with the PHD inhibitor DMOG in an orthotopic breast cancer model significantly decreases spontaneous metastases to the lungs and liver, associated with decreased tumour stiffness and fibroblast activation. PHD2 depletion in CAFs co‐injected with tumour cells similarly prevents CAF‐induced metastasis to lungs and liver. Our data argue that reversion of CAFs towards a less active state is possible and could have important clinical implications.     

Integrative genomics positions MKRN1 as a novel ribonucleoprotein within the embryonic stem cell gene regulatory network

In embryonic stem cells (ESCs), gene regulatory networks (GRNs) coordinate gene expression to maintain ESC identity; however, the complete repertoire of factors regulating the ESC state is not fully understood. Our previous temporal microarray analysis of ESC commitment identified the E3 ubiquitin ligase protein Makorin‐1 (MKRN1) as a potential novel component of the ESC GRN. Here, using multilayered systems‐level analyses, we compiled a MKRN1‐centered interactome in undifferentiated ESCs at the proteomic and ribonomic level. Proteomic analyses in undifferentiated ESCs revealed that MKRN1 associates with RNA‐binding proteins, and ensuing RIP‐chip analysis determined that MKRN1 associates with mRNAs encoding functionally related proteins including proteins that function during cellular stress. Subsequent biological validation identified MKRN1 as a novel stress granule‐resident protein, although MKRN1 is not required for stress granule formation, or survival of unstressed ESCs. Thus, our unbiased systems‐level analyses support a role for the E3 ligase MKRN1 as a ribonucleoprotein within the ESC GRN.     

Structure and function of the N‐terminal domain of the human mitochondrial calcium uniporter

The mitochondrial calcium uniporter (MCU) is responsible for mitochondrial calcium uptake and homeostasis. It is also a target for the regulation of cellular anti‐/pro‐apoptosis and necrosis by several oncogenes and tumour suppressors. Herein, we report the crystal structure of the MCU N‐terminal domain (NTD) at a resolution of 1.50 Å in a novel fold and the S92A MCU mutant at 2.75 Å resolution; the residue S92 is a predicted CaMKII phosphorylation site. The assembly of the mitochondrial calcium uniporter complex (uniplex) and the interaction with the MCU regulators such as the mitochondrial calcium uptake‐1 and mitochondrial calcium uptake‐2 proteins (MICU1 and MICU2) are not affected by the deletion of MCU NTD. However, the expression of the S92A mutant or a NTD deletion mutant failed to restore mitochondrial Ca²⁺ uptake in a stable MCU knockdown HeLa cell line and exerted dominant‐negative effects in the wild‐type MCU‐expressing cell line. These results suggest that the NTD of MCU is essential for the modulation of MCU function, although it does not affect the uniplex formation.     

Glucagon signalling in the dorsal vagal complex is sufficient and necessary for high‐protein feeding to regulate glucose homeostasis in vivo

High‐protein feeding acutely lowers postprandial glucose concentration compared to low‐protein feeding, despite a dichotomous rise of circulating glucagon levels. The physiological role of this glucagon rise has been largely overlooked. We here first report that glucagon signalling in the dorsal vagal complex (DVC) of the brain is sufficient to lower glucose production by activating a Gcgr–PKA–ERK–K_ATP channel signalling cascade in the DVC of rats in vivo. We further demonstrate that direct blockade of DVC Gcgr signalling negates the acute ability of high‐ vs. low‐protein feeding to reduce plasma glucose concentration, indicating that the elevated circulating glucagon during high‐protein feeding acts in the brain to lower plasma glucose levels. These data revise the physiological role of glucagon and argue that brain glucagon signalling contributes to glucose homeostasis during dietary protein intake.     

A comparative approach to testing hypotheses for the evolution of sex‐biased dispersal in bean beetles

Understanding the selective forces that shape dispersal strategies is a fundamental goal of evolutionary ecology and is increasingly important in changing, human‐altered environments. Sex‐biased dispersal (SBD) is common in dioecious taxa, and understanding variation in the direction and magnitude of SBD across taxa has been a persistent challenge. We took a comparative, laboratory‐based approach using 16 groups (species or strains) of bean beetles (genera Acanthoscelides, Callosobruchus, and Zabrotes, including 10 strains of one species) to test two predictions that emerge from dominant hypotheses for the evolution of SBD: (1) groups that suffer greater costs of inbreeding should exhibit greater SBD in favor of either sex (inbreeding avoidance hypothesis) and (2) groups with stronger local mate competition should exhibit greater male bias in dispersal (kin competition avoidance hypothesis). We used laboratory experiments to quantify SBD in crawling dispersal, the fitness effects of inbreeding, and the degree of polygyny (number of female mates per male), a proxy for local mate competition. While we found that both polygyny and male‐biased dispersal were common across bean beetle groups, consistent with the kin competition avoidance hypothesis, quantitative relationships between trait values did not support the predictions. Across groups, there was no significant association between SBD and effects of inbreeding nor SBD and degree of polygyny, using either raw values or phylogenetically independent contrasts. We discuss possible limitations of our experimental approach for detecting the predicted relationships, as well as reasons why single‐factor hypotheses may be too simplistic to explain the evolution of SBD.     

Systematic discovery of linear binding motifs targeting an ancient protein interaction surface on MAP kinases

Mitogen‐activated protein kinases (MAPK) are broadly used regulators of cellular signaling. However, how these enzymes can be involved in such a broad spectrum of physiological functions is not understood. Systematic discovery of MAPK networks both experimentally and in silico has been hindered because MAPKs bind to other proteins with low affinity and mostly in less‐characterized disordered regions. We used a structurally consistent model on kinase‐docking motif interactions to facilitate the discovery of short functional sites in the structurally flexible and functionally under‐explored part of the human proteome and applied experimental tools specifically tailored to detect low‐affinity protein–protein interactions for their validation in vitro and in cell‐based assays. The combined computational and experimental approach enabled the identification of many novel MAPK‐docking motifs that were elusive for other large‐scale protein–protein interaction screens. The analysis produced an extensive list of independently evolved linear binding motifs from a functionally diverse set of proteins. These all target, with characteristic binding specificity, an ancient protein interaction surface on evolutionarily related but physiologically clearly distinct three MAPKs (JNK, ERK, and p38). This inventory of human protein kinase binding sites was compared with that of other organisms to examine how kinase‐mediated partnerships evolved over time. The analysis suggests that most human MAPK‐binding motifs are surprisingly new evolutionarily inventions and newly found links highlight (previously hidden) roles of MAPKs. We propose that short MAPK‐binding stretches are created in disordered protein segments through a variety of ways and they represent a major resource for ancient signaling enzymes to acquire new regulatory roles.     

Gillespie eco‐evolutionary models (GEMs) reveal the role of heritable trait variation in eco‐evolutionary dynamics

Heritable trait variation is a central and necessary ingredient of evolution. Trait variation also directly affects ecological processes, generating a clear link between evolutionary and ecological dynamics. Despite the changes in variation that occur through selection, drift, mutation, and recombination, current eco‐evolutionary models usually fail to track how variation changes through time. Moreover, eco‐evolutionary models assume fitness functions for each trait and each ecological context, which often do not have empirical validation. We introduce a new type of model, Gillespie eco‐evolutionary models (GEMs), that resolves these concerns by tracking distributions of traits through time as eco‐evolutionary dynamics progress. This is done by allowing change to be driven by the direct fitness consequences of model parameters within the context of the underlying ecological model, without having to assume a particular fitness function. GEMs work by adding a trait distribution component to the standard Gillespie algorithm – an approach that models stochastic systems in nature that are typically approximated through ordinary differential equations. We illustrate GEMs with the Rosenzweig–MacArthur consumer–resource model. We show not only how heritable trait variation fuels trait evolution and influences eco‐evolutionary dynamics, but also how the erosion of variation through time may hinder eco‐evolutionary dynamics in the long run. GEMs can be developed for any parameter in any ordinary differential equation model and, furthermore, can enable modeling of multiple interacting traits at the same time. We expect GEMs will open the door to a new direction in eco‐evolutionary and evolutionary modeling by removing long‐standing modeling barriers, simplifying the link between traits, fitness, and dynamics, and expanding eco‐evolutionary treatment of a greater diversity of ecological interactions. These factors make GEMs much more than a modeling advance, but an important conceptual advance that bridges ecology and evolution through the central concept of heritable trait variation.     

Using Africa's protected area network to estimate the global population of a threatened and declining species: a case study of the Critically Endangered White‐headed Vulture Trigonoceps occipitalis

The White‐headed Vulture Trigonoceps occipitalis (WhV) is uncommon and largely restricted to protected areas across its range in sub‐Saharan Africa. We used the World Database on Protected Areas to identify protected areas (PAs) likely to contain White‐headed Vultures. Vulture occurrence on road transects in Southern, East, and West Africa was adjusted to nests per km² using data from areas with known numbers of nests and corresponding road transect data. Nest density was used to calculate the number of WhV nests within identified PAs and from there extrapolated to estimate the global population. Across a fragmented range, 400 PAs are estimated to contain 1893 WhV nests. Eastern Africa is estimated to contain 721 nests, Central Africa 548 nests, Southern Africa 468 nests, and West Africa 156 nests. Including immature and nonbreeding birds, and accounting for data deficient PAs, the estimated global population is 5475 ‐ 5493 birds. The identified distribution highlights are alarming: over 78% (n = 313) of identified PAs contain fewer than five nests. A further 17% (n = 68) of PAs contain 5 ‐ 20 nests and 4% (n = 14) of identified PAs are estimated to contain >20 nests. Just 1% (n = 5) of PAs are estimated to contain >40 nests; none is located in West Africa. Whilst ranging behavior of WhVs is currently unknown, 35% of PAs large enough to hold >20 nests are isolated by more than 100 km from other PAs. Spatially discrete and unpredictable mortality events such as poisoning pose major threats to small localized vulture populations and will accelerate ongoing local extinctions. Apart from reducing the threat of poisoning events, conservation actions promoting linkages between protected areas should be pursued. Identifying potential areas for assisted re‐establishment via translocation offers the potential to expand the range of this species and alleviate risk.     

Molecular and iridescent feather reflectance data reveal recent genetic diversification and phenotypic differentiation in a cloud forest hummingbird

The present day distribution and spatial genetic diversity of Mesoamerican biota reflects a long history of responses to habitat change. The hummingbird Lampornis amethystinus is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling across the species range using mitochondrial DNA (mtDNA) sequences and nuclear microsatellites jointly analysed with phenotypic and climatic data, we (1) test whether the fragmented distribution is correlated with main evolutionary lineages, (2) assess body size and plumage color differentiation of populations in geographic isolation, and (3) evaluate a set of divergence scenarios and demographic patterns of the hummingbird populations. Analysis of genetic variation revealed four main groups: blue‐throated populations (Sierra Madre del Sur); two groups of amethyst‐throated populations (Trans‐Mexican Volcanic Belt and Sierra Madre Oriental); and populations east of the Isthmus of Tehuantepec (IT) with males showing an amethyst throat. The most basal split is estimated to have originated in the Pleistocene, 2.39–0.57 million years ago (MYA), and corresponded to groups of populations separated by the IT. However, the estimated recent divergence time between blue‐ and amethyst‐throated populations does not correspond to the 2‐MY needed to be in isolation for substantial plumage divergence, likely because structurally iridescent colors are more malleable than others. Results of species distribution modeling and Approximate Bayesian Computation analysis fit a model of lineage divergence west of the Isthmus after the Last Glacial Maximum (LGM), and that the species’ suitable habitat was disjunct during past and current conditions. These results challenge the generality of the contraction/expansion glacial model to cloud forest‐interior species and urges management of cloud forest, a highly vulnerable ecosystem to climate change and currently facing destruction, to prevent further loss of genetic diversity or extinction.     

N‐glycosylation of PD‐1 promotes binding of camrelizumab

PD‐1 is a highly glycosylated inhibitory receptor expressed mainly on T cells. Targeting of PD‐1 with monoclonal antibodies (MAbs) to block the interaction with its ligand PD‐L1 has been successful for the treatment of multiple tumors. However, polymorphisms at N‐glycosylation sites of PD‐1 exist in the human population that might affect antibody binding, and dysregulated glycosylation has been observed in the tumor microenvironment. Here, we demonstrate varied N‐glycan composition in PD‐1, and show that the binding affinity of camrelizumab, a recently approved PD‐1‐specific MAb, to non‐glycosylated PD‐1 proteins from E. coli is substantially decreased compared with glycosylated PD‐1. The structure of the camrelizumab/PD‐1 complex reveals that camrelizumab mainly utilizes its heavy chain to bind to PD‐1, while the light chain sterically inhibits the binding of PD‐L1 to PD‐1. Glycosylation of asparagine 58 (N58) promotes the interaction with camrelizumab, while the efficiency of camrelizumab to inhibit the binding of PD‐L1 is substantially reduced for glycosylation‐deficient PD‐1. These results increase our understanding of how glycosylation affects the activity of PD‐1‐specific MAbs during immune checkpoint therapy.     

Contribution of intertwined loop to membrane association revealed by Zika virus full‐length NS1 structure

The association of Zika virus (ZIKV) infections with microcephaly and neurological diseases has highlighted an emerging public health concern. Here, we report the crystal structure of the full‐length ZIKV nonstructural protein 1 (NS1), a major host‐interaction molecule that functions in flaviviral replication, pathogenesis, and immune evasion. Of note, a long intertwined loop is observed in the wing domain of ZIKV NS1, and forms a hydrophobic “spike”, which can contribute to cellular membrane association. For different flaviviruses, the amino acid sequences of the “spike” are variable but their common characteristic is either hydrophobic or positively charged, which is a beneficial feature for membrane binding. Comparative studies with West Nile and Dengue virus NS1 structures reveal conserved features, but diversified electrostatic characteristics on both inner and outer faces. Our results suggest different mechanisms of flavivirus pathogenesis and should be considered during the development of diagnostic tools.