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.     

Africanization of a feral honey bee (Apis mellifera) population in South Texas: does a decade make a difference?

The arrival to the United States of the Africanized honey bee, a hybrid between European subspecies and the African subspecies Apis mellifera scutellata, is a remarkable model for the study of biological invasions. This immigration has created an opportunity to study the dynamics of secondary contact of honey bee subspecies from African and European lineages in a feral population in South Texas. An 11‐year survey of this population (1991–2001) showed that mitochondrial haplotype frequencies changed drastically over time from a resident population of eastern and western European maternal ancestry, to a population dominated by the African haplotype. A subsequent study of the nuclear genome showed that the Africanization process included bidirectional gene flow between European and Africanized honey bees, giving rise to a new panmictic mixture of A. m. scutellata‐ and European‐derived genes. In this study, we examined gene flow patterns in the same population 23 years after the first hybridization event occurred. We found 28 active colonies inhabiting 92 tree cavities surveyed in a 5.14 km² area, resulting in a colony density of 5.4 colonies/km². Of these 28 colonies, 25 were of A. m. scutellata maternal ancestry, and three were of western European maternal ancestry. No colonies of eastern European maternal ancestry were detected, although they were present in the earlier samples. Nuclear DNA revealed little change in the introgression of A. m. scutellata‐derived genes into the population compared to previous surveys. Our results suggest this feral population remains an admixed swarm with continued low levels of European ancestry and a greater presence of African‐derived mitochondrial genetic composition.     

MicroRNA‐455 regulates brown adipogenesis via a novel HIF1an‐AMPK‐PGC1α signaling network

Brown adipose tissue (BAT) dissipates chemical energy as heat and can counteract obesity. MicroRNAs are emerging as key regulators in development and disease. Combining microRNA and mRNA microarray profiling followed by bioinformatic analyses, we identified miR‐455 as a new regulator of brown adipogenesis. miR‐455 exhibits a BAT‐specific expression pattern and is induced by cold and the browning inducer BMP7. In vitro gain‐ and loss‐of‐function studies show that miR‐455 regulates brown adipocyte differentiation and thermogenesis. Adipose‐specific miR‐455 transgenic mice display marked browning of subcutaneous white fat upon cold exposure. miR‐455 activates AMPKα1 by targeting HIF1an, and AMPK promotes the brown adipogenic program and mitochondrial biogenesis. Concomitantly, miR‐455 also targets the adipogenic suppressors Runx1t1 and Necdin, initiating adipogenic differentiation. Taken together, the data reveal a novel microRNA‐regulated signaling network that controls brown adipogenesis and may be a potential therapeutic target for human metabolic disorders.     

Reversible HuR‐microRNA binding controls extracellular export of miR‐122 and augments stress response

microRNAs (miRNAs), the tiny but stable regulatory RNAs in metazoan cells, can undergo selective turnover in presence of specific internal and external cues to control cellular response against the changing environment. We have observed reduction in cellular miR‐122 content, due to their accelerated extracellular export in human hepatic cells starved for small metabolites including amino acids. In this context, a new role of human ELAV protein HuR has been identified. HuR, a negative regulator of miRNA function, accelerates extracellular vesicle (EV)‐mediated export of miRNAs in human cells. In stressed cells, HuR replaces miRNPs from target messages and is both necessary and sufficient for the extracellular export of corresponding miRNAs. HuR could reversibly bind miRNAs to replace them from Ago2 and subsequently itself gets freed from bound miRNAs upon ubiquitination. The ubiquitinated form of HuR is predominantly associated with multivesicular bodies (MVB) where HuR‐unbound miRNAs also reside. These MVB‐associated pool of miRNAs get exported out via EVs thereby delimiting cellular miR‐122 level during starvation. Therefore, by modulating extracellular export of miR‐122, HuR could control stress response in starved human hepatic cells.     

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.     

Molecular architecture of the ribosome‐bound Hepatitis C Virus internal ribosomal entry site RNA

Internal ribosomal entry sites (IRESs) are structured cis‐acting RNAs that drive an alternative, cap‐independent translation initiation pathway. They are used by many viruses to hijack the translational machinery of the host cell. IRESs facilitate translation initiation by recruiting and actively manipulating the eukaryotic ribosome using only a subset of canonical initiation factor and IRES transacting factors. Here we present cryo‐EM reconstructions of the ribosome 80S‐ and 40S‐bound Hepatitis C Virus (HCV) IRES. The presence of four subpopulations for the 80S•HCV IRES complex reveals dynamic conformational modes of the complex. At a global resolution of 3.9 Å for the most stable complex, a derived atomic model reveals a complex fold of the IRES RNA and molecular details of its interaction with the ribosome. The comparison of obtained structures explains how a modular architecture facilitates mRNA loading and tRNA binding to the P‐site. This information provides the structural foundation for understanding the mechanism of HCV IRES RNA‐driven translation initiation.     

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.     

Human TLR8 senses UR/URR motifs in bacterial and mitochondrial RNA

Toll‐like receptor (TLR) 13 and TLR2 are the major sensors of Gram‐positive bacteria in mice. TLR13 recognizes Sa19, a specific 23S ribosomal (r) RNA‐derived fragment and bacterial modification of Sa19 ablates binding to TLR13, and to antibiotics such as erythromycin. Similarly, RNase A‐treated Staphylococcus aureus activate human peripheral blood mononuclear cells (PBMCs) only via TLR2, implying single‐stranded (ss) RNA as major stimulant. Here, we identify human TLR8 as functional TLR13 equivalent that promiscuously senses ssRNA. Accordingly, Sa19 and mitochondrial (mt) 16S rRNA sequence‐derived oligoribonucleotides (ORNs) stimulate PBMCs in a MyD88‐dependent manner. These ORNs, as well as S. aureus‐, Escherichia coli‐, and mt‐RNA, also activate differentiated human monocytoid THP‐1 cells, provided they express TLR8. Moreover, Unc93b1 ^−/−‐ and Tlr8 ^−/−‐THP‐1 cells are refractory, while endogenous and ectopically expressed TLR8 confers responsiveness in a UR/URR RNA ligand consensus motif‐dependent manner. If TLR8 function is inhibited by suppression of lysosomal function, antibiotic treatment efficiently blocks bacteria‐driven inflammatory responses in infected human whole blood cultures. Sepsis therapy might thus benefit from interfering with TLR8 function.     

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.     

Little or no gene flow despite F1 hybrids at two interspecific contact zones

Hybridization can create the selective force that promotes assortative mating but hybridization can also select for increased hybrid fitness. Gene flow resulting from hybridization can increase genetic diversity but also reduce distinctiveness. Thus the formation of hybrids has important implications for long‐term species coexistence. This study compares the interaction between the tree wētā Hemideina thoracica and its two neighboring species; H. crassidens and H. trewicki. We examined the ratio of parent and hybrid forms in natural areas of sympatry. Individuals with intermediate phenotype were confirmed as first generation hybrids using nine independent genetic markers. Evidence of gene flow from successful hybridization was sought from the distribution of morphological and genetic characters. Both species pairs appear to be largely retaining their own identity where they live in sympatry, each with a distinct karyotype. Hemideina thoracica and H. trewicki are probably reproductively isolated, with sterile F₁ hybrids. This species pair shows evidence of niche differences with adult size and timing of maturity differing where Hemideina thoracica is sympatric with H. trewicki. In contrast, evidence of a low level of introgression was detected in phenotypes and genotypes where H. thoracica and H. crassidens are sympatric. We found no evidence of size divergence although color traits in combination with hind tibia spines reliably distinguish the two species. This species pair show a bimodal hybrid zone in the absence of assortative mating and possible sexual exclusion by H. thoracica males in the formation of F₁ hybrids.     

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.     

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.     

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.     

Phosphoproteomic screening identifies Rab GTPases as novel downstream targets of PINK1

Mutations in the PTEN‐induced kinase 1 (PINK1) are causative of autosomal recessive Parkinson''s disease (PD). We have previously reported that PINK1 is activated by mitochondrial depolarisation and phosphorylates serine 65 (Ser⁶⁵) of the ubiquitin ligase Parkin and ubiquitin to stimulate Parkin E3 ligase activity. Here, we have employed quantitative phosphoproteomics to search for novel PINK1‐dependent phosphorylation targets in HEK (human embryonic kidney) 293 cells stimulated by mitochondrial depolarisation. This led to the identification of 14,213 phosphosites from 4,499 gene products. Whilst most phosphosites were unaffected, we strikingly observed three members of a sub‐family of Rab GTPases namely Rab8A, 8B and 13 that are all phosphorylated at the highly conserved residue of serine 111 (Ser¹¹¹) in response to PINK1 activation. Using phospho‐specific antibodies raised against Ser¹¹¹ of each of the Rabs, we demonstrate that Rab Ser¹¹¹ phosphorylation occurs specifically in response to PINK1 activation and is abolished in HeLa PINK1 knockout cells and mutant PINK1 PD patient‐derived fibroblasts stimulated by mitochondrial depolarisation. We provide evidence that Rab8A GTPase Ser¹¹¹ phosphorylation is not directly regulated by PINK1 in vitro and demonstrate in cells the time course of Ser¹¹¹ phosphorylation of Rab8A, 8B and 13 is markedly delayed compared to phosphorylation of Parkin at Ser⁶⁵. We further show mechanistically that phosphorylation at Ser¹¹¹ significantly impairs Rab8A activation by its cognate guanine nucleotide exchange factor (GEF), Rabin8 (by using the Ser111Glu phosphorylation mimic). These findings provide the first evidence that PINK1 is able to regulate the phosphorylation of Rab GTPases and indicate that monitoring phosphorylation of Rab8A/8B/13 at Ser¹¹¹ may represent novel biomarkers of PINK1 activity in vivo. Our findings also suggest that disruption of Rab GTPase‐mediated signalling may represent a major mechanism in the neurodegenerative cascade of Parkinson''s disease.     

Population and genetic outcomes 20 years after reintroducing bobcats (Lynx rufus) to Cumberland Island,Georgia USA

In 1988–1989, 32 bobcats Lynx rufus were reintroduced to Cumberland Island (CUIS), Georgia, USA, from which they had previously been extirpated. They were monitored intensively for 3 years immediately post‐reintroduction, but no estimation of the size or genetic diversity of the population had been conducted in over 20 years since reintroduction. We returned to CUIS in 2012 to estimate abundance and effective population size of the present‐day population, as well as to quantify genetic diversity and inbreeding. We amplified 12 nuclear microsatellite loci from DNA isolated from scats to establish genetic profiles to identify individuals. We used spatially explicit capture–recapture population estimation to estimate abundance. From nine unique genetic profiles, we estimate a population size of 14.4 (SE = 3.052) bobcats, with an effective population size (N _e) of 5–8 breeding individuals. This is consistent with predictions of a population viability analysis conducted at the time of reintroduction, which estimated the population would average 12–13 bobcats after 10 years. We identified several pairs of related bobcats (parent‐offspring and full siblings), but ~75% of the pairwise comparisons were typical of unrelated individuals, and only one individual appeared inbred. Despite the small population size and other indications that it has likely experienced a genetic bottleneck, levels of genetic diversity in the CUIS bobcat population remain high compared to other mammalian carnivores. The reintroduction of bobcats to CUIS provides an opportunity to study changes in genetic diversity in an insular population without risk to this common species. Opportunities for natural immigration to the island are limited; therefore, continued monitoring and supplemental bobcat reintroductions could be used to evaluate the effect of different management strategies to maintain genetic diversity and population viability. The successful reintroduction and maintenance of a bobcat population on CUIS illustrates the suitability of translocation as a management tool for re‐establishing felid populations.