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.     

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.     

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.     

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.     

Testosterone and cortisol concentrations vary with reproductive status in wild female red deer

Although hormones are key regulators of many fitness and life history traits, the causes of individual level variation in hormones, particularly in wild systems, remain understudied. Whilst we know that androgen and glucocorticoid levels vary within and among individuals in mammalian populations, how this relates to key reproductive processes such as gestation and lactation, and their effects on a female''s measurable hormone levels are poorly understood in wild systems. Using fecal samples collected from females in a wild red deer population between 2001 and 2013, we explore how fecal androgen (FAM) and cortisol (FCM) metabolite concentrations change with age and season, and how individual differences relate to variation in reproductive state. Both FAM and FCM levels increase toward parturition, although this only affects FCM levels in older females. FCM levels are also higher when females suckle a male rather than a female calf, possibly due to the higher energetic costs of raising a son. This illustrates the importance of accounting for a female''s life history and current reproductive status, as well as temporal variation, when examining individual differences in hormone levels. We discuss these findings in relation to other studies of mammalian systems and in particular to the relatively scarce information on variation in natural levels of hormones in wild populations.     

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.     

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.     

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.     

The C9orf72 protein interacts with Rab1a and the ULK1 complex to regulate initiation of autophagy

A GGGGCC hexanucleotide repeat expansion in the C9orf72 gene is the most common genetic cause of amyotrophic lateral sclerosis and frontotemporal dementia (C9ALS/FTD). C9orf72 encodes two C9orf72 protein isoforms of unclear function. Reduced levels of C9orf72 expression have been reported in C9ALS/FTD patients, and although C9orf72 haploinsufficiency has been proposed to contribute to C9ALS/FTD, its significance is not yet clear. Here, we report that C9orf72 interacts with Rab1a and the Unc‐51‐like kinase 1 (ULK1) autophagy initiation complex. As a Rab1a effector, C9orf72 controls initiation of autophagy by regulating the Rab1a‐dependent trafficking of the ULK1 autophagy initiation complex to the phagophore. Accordingly, reduction of C9orf72 expression in cell lines and primary neurons attenuated autophagy and caused accumulation of p62‐positive puncta reminiscent of the p62 pathology observed in C9ALS/FTD patients. Finally, basal levels of autophagy were markedly reduced in C9ALS/FTD patient‐derived iNeurons. Thus, our data identify C9orf72 as a novel Rab1a effector in the regulation of autophagy and indicate that C9orf72 haploinsufficiency and associated reductions in autophagy might be the underlying cause of C9ALS/FTD‐associated p62 pathology.     

Improving microbial fitness in the mammalian gut by in vivo temporal functional metagenomics

Elucidating functions of commensal microbial genes in the mammalian gut is challenging because many commensals are recalcitrant to laboratory cultivation and genetic manipulation. We present Temporal FUnctional Metagenomics sequencing (TFUMseq), a platform to functionally mine bacterial genomes for genes that contribute to fitness of commensal bacteria in vivo. Our approach uses metagenomic DNA to construct large‐scale heterologous expression libraries that are tracked over time in vivo by deep sequencing and computational methods. To demonstrate our approach, we built a TFUMseq plasmid library using the gut commensal Bacteroides thetaiotaomicron (Bt) and introduced Escherichia coli carrying this library into germfree mice. Population dynamics of library clones revealed Bt genes conferring significant fitness advantages in E. coli over time, including carbohydrate utilization genes, with a Bt galactokinase central to early colonization, and subsequent dominance by a Bt glycoside hydrolase enabling sucrose metabolism coupled with co‐evolution of the plasmid library and E. coli genome driving increased galactose utilization. Our findings highlight the utility of functional metagenomics for engineering commensal bacteria with improved properties, including expanded colonization capabilities in vivo.     

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.     

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.