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.     

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.     

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.     

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.     

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.     

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.     

Tree species richness decreases while species evenness increases with disturbance frequency in a natural boreal forest landscape

Understanding species diversity and disturbance relationships is important for biodiversity conservation in disturbance‐driven boreal forests. Species richness and evenness may respond differently with stand development following fire. Furthermore, few studies have simultaneously accounted for the influences of climate and local site conditions on species diversity. Using forest inventory data, we examined the relationships between species richness, Shannon''s index, evenness, and time since last stand‐replacing fire (TSF) in a large landscape of disturbance‐driven boreal forest. TSF has negative effect on species richness and Shannon''s index, and a positive effect on species evenness. Path analysis revealed that the environmental variables affect richness and Shannon''s index only through their effects on TSF while affecting evenness directly as well as through their effects on TSF. Synthesis and applications. Our results demonstrate that species richness and Shannon''s index decrease while species evenness increases with TSF in a boreal forest landscape. Furthermore, we show that disturbance frequency, local site conditions, and climate simultaneously influence tree species diversity through complex direct and indirect effects in the studied boreal forest.     

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.     

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.     

Gruffi: an algorithm for computational removal of stressed cells from brain organoid transcriptomic datasets

Organoids enable in vitro modeling of complex developmental processes and disease pathologies. Like most 3D cultures, organoids lack sufficient oxygen supply and therefore experience cellular stress. These negative effects are particularly prominent in complex models, such as brain organoids, and can affect lineage commitment. Here, we analyze brain organoid and fetal single‐cell RNA sequencing (scRNAseq) data from published and new datasets, totaling about 190,000 cells. We identify a unique stress signature in the data from all organoid samples, but not in fetal samples. We demonstrate that cell stress is limited to a defined subpopulation of cells that is unique to organoids and does not affect neuronal specification or maturation. We have developed a computational algorithm, Gruffi, which uses granular functional filtering to identify and remove stressed cells from any organoid scRNAseq dataset in an unbiased manner. We validated our method using six additional datasets from different organoid protocols and early brains, and show its usefulness to other organoid systems including retinal organoids. Our data show that the adverse effects of cell stress can be corrected by bioinformatic analysis for improved delineation of developmental trajectories and resemblance to in vivo data.     

Summer temperature can predict the distribution of wild yeast populations

The wine yeast, Saccharomyces cerevisiae, is the best understood microbial eukaryote at the molecular and cellular level, yet its natural geographic distribution is unknown. Here we report the results of a field survey for S. cerevisiae,S. paradoxus and other budding yeast on oak trees in Europe. We show that yeast species differ in their geographic distributions, and investigated which ecological variables can predict the isolation rate of S. paradoxus, the most abundant species. We find a positive association between trunk girth and S. paradoxus abundance suggesting that older trees harbor more yeast. S. paradoxus isolation frequency is also associated with summer temperature, showing highest isolation rates at intermediate temperatures. Using our statistical model, we estimated a range of summer temperatures at which we expect high S. paradoxus isolation rates, and show that the geographic distribution predicted by this optimum temperature range is consistent with the worldwide distribution of sites where S. paradoxus has been isolated. Using laboratory estimates of optimal growth temperatures for S. cerevisiae relative to S. paradoxus, we also estimated an optimum range of summer temperatures for S. cerevisiae. The geographic distribution of these optimum temperatures is consistent with the locations where wild S. cerevisiae have been reported, and can explain why only human‐associated S. cerevisiae strains are isolated at northernmost latitudes. Our results provide a starting point for targeted isolation of S. cerevisiae from natural habitats, which could lead to a better understanding of climate associations and natural history in this important model microbe.     

Population genetic analyses are consistent with the introduction of Ceramium secundatum (Ceramiaceae,Rhodophyta) to Narragansett Bay,Rhode Island,USA

During ongoing DNA barcode (COI‐5P) surveys of the macroalgal flora along the northwest Atlantic coast, we discovered a population of Ceramium secundatum in Narragansett Bay, Rhode Island, USA. This species is regarded as common and widespread in the northeast Atlantic, ranging from Norway to Morocco, but until now has not been reported from the western Atlantic. Several lines of evidence suggest that C. secundatum may be introduced to Narragansett Bay: (1) despite extensive collecting, specimens have only been obtained from a limited geographic range in the northwest Atlantic; (2) three other nonindigenous seaweed species are reportedly introduced in this region, which is thought to be a consequence of shipping; and (3) this species is introduced to South Africa and New Zealand. To investigate this suspected introduction, we applied population genetic analyses (using the cox2‐3 spacer) to compare the Narragansett Bay C. secundatum population to native populations in the Republic of Ireland and the United Kingdom. Collectively, analyses of biogeographical and molecular data indicate that C. secundatum is likely introduced to Narragansett Bay. The implications of this discovery are discussed.