The 50 most-cited articles in eLife, 2016

Being taller is associated with enhanced longevity, and higher education and earnings. We reanalysed 1472 population-based studies, with measurement of height on more than 18.6 million participants to estimate mean height for people born between 1896 and 1996 in 200 countries. The largest gain in adult height over the past century has occurred in South Korean women and Iranian men, who became 20.2 cm (95% credible interval 17.5-22.7) and 16.5 cm (13.3-19.7) taller, respectively. In contrast, there was little change in adult height in some sub-Saharan African countries and in South Asia over the century of analysis. The tallest people over these 100 years are men born in the Netherlands in the last quarter of 20th century, whose average heights surpassed 182.5 cm, and the shortest were women born in Guatemala in 1896 (140.3 cm; 135.8-144.8). The height differential between the tallest and shortest populations was 19-20 cm a century ago, and has remained the same for women and increased for men a century later despite substantial changes in the ranking of countries.

Genetically encoded calcium indicators (GECIs) allow measurement of activity in large populations of neurons and in small neuronal compartments, over times of milliseconds to months. Although GFP-based GECIs are widely used for in vivo neurophysiology, GECIs with red-shifted excitation and emission spectra have advantages for in vivo imaging because of reduced scattering and absorption in tissue, and a consequent reduction in phototoxicity. However, current red GECIs are inferior to the state-of-the-art GFP-based GCaMP6 indicators for detecting and quantifying neural activity. Here we present improved red GECIs based on mRuby (jRCaMP1a, b) and mApple (jRGECO1a), with sensitivity comparable to GCaMP6. We characterized the performance of the new red GECIs in cultured neurons and in mouse, Drosophila, zebrafish and C. elegans in vivo. Red GECIs facilitate deep-tissue imaging, dual-color imaging together with GFP-based reporters, and the use of optogenetics in combination with calcium imaging.

Mutations in Park8, encoding for the multidomain Leucine-rich repeat kinase 2 (LRRK2) protein, comprise the predominant genetic cause of Parkinson's disease (PD). G2019S, the most common amino acid substitution activates the kinase two- to threefold. This has motivated the development of LRRK2 kinase inhibitors; however, poor consensus on physiological LRRK2 substrates has hampered clinical development of such therapeutics. We employ a combination of phosphoproteomics, genetics, and pharmacology to unambiguously identify a subset of Rab GTPases as key LRRK2 substrates. LRRK2 directly phosphorylates these both in vivo and in vitro on an evolutionary conserved residue in the switch II domain. Pathogenic LRRK2 variants mapping to different functional domains increase phosphorylation of Rabs and this strongly decreases their affinity to regulatory proteins including Rab GDP dissociation inhibitors (GDIs). Our findings uncover a key class of bona-fide LRRK2 substrates and a novel regulatory mechanism of Rabs that connects them to PD.

Cancer-associated fibroblasts (CAFs) are a major cellular component of tumor microenvironment in most solid cancers. Altered cellular metabolism is a hallmark of cancer, and much of the published literature has focused on neoplastic cell-autonomous processes for these adaptations. We demonstrate that exosomes secreted by patient-derived CAFs can strikingly reprogram the metabolic machinery following their uptake by cancer cells. We find that CAF-derived exosomes (CDEs) inhibit mitochondrial oxidative phosphorylation, thereby increasing glycolysis and glutamine-dependent reductive carboxylation in cancer cells. Through 13C-labeled isotope labeling experiments we elucidate that exosomes supply amino acids to nutrient-deprived cancer cells in a mechanism similar to macropinocytosis, albeit without the previously described dependence on oncogenic-Kras signaling. Using intra-exosomal metabolomics, we provide compelling evidence that CDEs contain intact metabolites, including amino acids, lipids, and TCA-cycle intermediates that are avidly utilized by cancer cells for central carbon…

By reaching near-atomic resolution for a wide range of specimens, single-particle cryo-EM structure determination is transforming structural biology. However, the necessary calculations come at large computational costs, which has introduced a bottleneck that is currently limiting throughput and the development of new methods. Here, we present an implementation of the RELION image processing software that uses graphics processors (GPUs) to address the most computationally intensive steps of its cryo-EM structure determination workflow. Both image classification and high-resolution refinement have been accelerated more than an order-of-magnitude, and template-based particle selection has been accelerated well over two orders-of-magnitude on desktop hardware. Memory requirements on GPUs have been reduced to fit widely available hardware, and we show that the use of single precision arithmetic does not adversely affect results. This enables high-resolution cryo-EM structure determination in a matter of days on a single workstation.

We recently found that nucleosomes directly block access of CRISPR/Cas9 to DNA (Horlbeck et al., 2016). Here, we build on this observation with a comprehensive algorithm that incorporates chromatin, position, and sequence features to accurately predict highly effective single guide RNAs (sgRNAs) for targeting nuclease-dead Cas9-mediated transcriptional repression (CRISPRi) and activation (CRISPRa). We use this algorithm to design next-generation genome-scale CRISPRi and CRISPRa libraries targeting human and mouse genomes. A CRISPRi screen for essential genes in K562 cells demonstrates that the large majority of sgRNAs are highly active. We also find CRISPRi does not exhibit any detectable non-specific toxicity recently observed with CRISPR nuclease approaches. Precision-recall analysis shows that we detect over 90% of essential genes with minimal false positives using a compact 5 sgRNA/gene library. Our results establish CRISPRi and CRISPRa as premier tools for loss- or gain-of-function studies and provide a general strategy for identifying Cas9 target sites.

Exercise induces beneficial responses in the brain, which is accompanied by an increase in BDNF, a trophic factor associated with cognitive improvement and the alleviation of depression and anxiety. However, the exact mechanisms whereby physical exercise produces an induction in brain Bdnf gene expression are not well understood. While pharmacological doses of HDAC inhibitors exert positive effects on Bdnf gene transcription, the inhibitors represent small molecules that do not occur in vivo. Here, we report that an endogenous molecule released after exercise is capable of inducing key promoters of the Mus musculus Bdnf gene. The metabolite β-hydroxybutyrate, which increases after prolonged exercise, induces the activities of Bdnf promoters, particularly promoter I, which is activity-dependent. We have discovered that the action of β-hydroxybutyrate is specifically upon HDAC2 and HDAC3, which act upon selective Bdnf promoters. Moreover, the effects upon hippocampal Bdnf expression were observed after direct ventricular application of β-hydroxybutyrate. Electrophysiological…

Stress granules are non-membrane bound RNA-protein (RNP) assemblies that form when translation initiation is limited and contain a biphasic structure with stable core structures surrounded by a less concentrated shell. The order of assembly and disassembly of these two structures remains unknown. Time course analysis of granule assembly suggests that core formation is an early event in granule assembly. Stress granule disassembly is also a stepwise process with shell dissipation followed by core clearance. Perturbations that alter liquid-liquid phase separations (LLPS) driven by intrinsically disordered protein regions (IDR) of RNA binding proteins in vitro have the opposite effect on stress granule assembly in vivo. Taken together, these observations argue that stress granules assemble through a multistep process initiated by stable assembly of untranslated mRNPs into core structures, which could provide sufficient high local concentrations to allow for a localized LLPS driven by IDRs on RNA binding proteins.

Open access, open data, open source and other open scholarship practices are growing in popularity and necessity. However, widespread adoption of these practices has not yet been achieved. One reason is that researchers are uncertain about how sharing their work will affect their careers. We review literature demonstrating that open research is associated with increases in citations, media attention, potential collaborators, job opportunities and funding opportunities. These findings are evidence that open research practices bring significant benefits to researchers relative to more traditional closed practices.

How metabolism is reprogrammed during neuronal differentiation is unknown. We found that the loss of hexokinase (HK2) and lactate dehydrogenase (LDHA) expression, together with a switch in pyruvate kinase gene splicing from PKM2 to PKM1, marks the transition from aerobic glycolysis in neural progenitor cells (NPC) to neuronal oxidative phosphorylation. The protein levels of c-MYC and N-MYC, transcriptional activators of the HK2 and LDHA genes, decrease dramatically. Constitutive expression of HK2 and LDHA during differentiation leads to neuronal cell death, indicating that the shut-off aerobic glycolysis is essential for neuronal survival. The metabolic regulators PGC-1α and ERRγ increase significantly upon neuronal differentiation to sustain the transcription of metabolic and mitochondrial genes, whose levels are unchanged compared to NPCs, revealing distinct transcriptional regulation of metabolic genes in the proliferation and post-mitotic differentiation states. Mitochondrial mass increases proportionally with neuronal mass growth, indicating an unknown mechanism linking…

Exosomes are small vesicles that are secreted from metazoan cells and may convey selected membrane proteins and small RNAs to target cells for the control of cell migration, development and metastasis. To study the mechanisms of RNA packaging into exosomes, we devised a purification scheme based on the membrane marker CD63 to isolate a single exosome species secreted from HEK293T cells. Using immunoisolated CD63-containing exosomes we identified a set of miRNAs that are highly enriched with respect to their cellular levels. To explore the biochemical requirements for exosome biogenesis and RNA packaging, we devised a cell-free reaction that recapitulates the species-selective enclosure of miR-223 in isolated membranes supplemented with cytosol. We found that the RNA-binding protein Y-box protein I (YBX1) binds to and is required for the sorting of miR-223 in the cell-free reaction. Furthermore, YBX1 serves an important role in the secretion of miRNAs in exosomes by HEK293T cells.

Abstract Imaging is used to map activity across populations of neurons. Microscopes with cellular resolution have small (< 1 millimeter) fields of view and cannot simultaneously image activity distributed across multiple brain areas. Typical large field of view microscopes do not resolve single cells, especially in the axial dimension. We developed a 2-photon random access mesoscope (2p-RAM) that allows high-resolution imaging anywhere within a volume spanning multiple brain areas (Ø 5 mm × 1 mm cylinder). 2p-RAM resolution is near diffraction limited (lateral, 0.66 μm, axial 4.09 μm at the center; excitation wavelength = 970 nm; numerical aperture = 0.6) over a large range of excitation wavelengths. A fast threedimensional scanning system allows efficient sampling of neural activity in arbitrary regions of interest across the entire imaging volume. We illustrate the use of the 2p-RAM by imaging neural activity in multiple, non-contiguous brain areas in transgenic mice expressing protein calcium sensors.

Subcellular localization critically influences protein function, and cells control protein localization to regulate biological processes. We have developed and applied Dynamic Organellar Maps, a proteomic method that allows global mapping of protein translocation events. We initially used maps statically to generate a database with localization and absolute copy number information for over 8700 proteins from HeLa cells, approaching comprehensive coverage. All major organelles were resolved, with exceptional prediction accuracy (estimated at >92%). Combining spatial and abundance information yielded an unprecedented quantitative view of HeLa cell anatomy and organellar composition, at the protein level. We subsequently demonstrated the dynamic capabilities of the approach by capturing translocation events following EGF stimulation, which we integrated into a quantitative model. Dynamic Organellar Maps enable the proteome-wide analysis of physiological protein movements, without requiring any reagents specific to the investigated process, and will thus be widely applicable in cell…

Neurons in higher cortical areas, such as the prefrontal cortex, are often tuned to a variety of sensory and motor variables, and are therefore said to display mixed selectivity. This complexity of single neuron responses can obscure what information these areas represent and how it is represented. Here we demonstrate the advantages of a new dimensionality reduction technique, demixed principal component analysis (dPCA), that decomposes population activity into a few components. In addition to systematically capturing the majority of the variance of the data, dPCA also exposes the dependence of the neural representation on task parameters such as stimuli, decisions, or rewards. To illustrate our method we reanalyze population data from four datasets comprising different species, different cortical areas and different experimental tasks. In each case, dPCA provides a concise way of visualizing the data that summarizes the task-dependent features of the population response in a single figure.

Prominent theories suggest that compulsive behaviors, characteristic of obsessive-compulsive disorder and addiction, are driven by shared deficits in goal-directed control, which confers vulnerability for developing rigid habits. However, recent studies have shown that deficient goal-directed control accompanies several disorders, including those without an obvious compulsive element. Reasoning that this lack of clinical specificity might reflect broader issues with psychiatric diagnostic categories, we investigated whether a dimensional approach would better delineate the clinical manifestations of goal-directed deficits. Using large-scale online assessment of psychiatric symptoms and neurocognitive performance in two independent general-population samples, we found that deficits in goal-directed control were most strongly associated with a symptom dimension comprising compulsive behavior and intrusive thought. This association was highly specific when compared to other non-compulsive aspects of psychopathology. These data showcase a powerful new methodology and highlight the…

The striatum integrates excitatory inputs from the cortex and the thalamus to control diverse functions. Although the striatum is thought to consist of sensorimotor, associative and limbic domains, their precise demarcations and whether additional functional subdivisions exist remain unclear. How striatal inputs are differentially segregated into each domain is also poorly understood. This study presents a comprehensive map of the excitatory inputs to the mouse striatum. The input patterns reveal boundaries between the known striatal domains. The most posterior striatum likely represents the 4th functional subdivision, and the dorsomedial striatum integrates highly heterogeneous, multimodal inputs. The complete thalamo-cortico-striatal loop is also presented, which reveals that the thalamic subregions innervated by the basal ganglia preferentially interconnect with motor-related cortical areas. Optogenetic experiments show the subregion-specific heterogeneity in the synaptic properties of striatal inputs from both the cortex and the thalamus. This projectome will guide functional…

Mutations, deletions, and changes in copy number of mitochondrial DNA (mtDNA), are observed throughout cancers. Here, we survey mtDNA copy number variation across 22 tumor types profiled by The Cancer Genome Atlas project. We observe a tendency for some cancers, especially of the bladder, breast, and kidney, to be depleted of mtDNA, relative to matched normal tissue. Analysis of genetic context reveals an association between incidence of several somatic alterations, including IDH1 mutations in gliomas, and mtDNA content. In some but not all cancer types, mtDNA content is correlated with the expression of respiratory genes, and anti-correlated to the expression of immune response and cell-cycle genes. In tandem with immunohistochemical evidence, we find that some tumors may compensate for mtDNA depletion to sustain levels of respiratory proteins. Our results highlight the extent of mtDNA copy number variation in tumors and point to related therapeutic opportunities.

Methylation of adenosines at the N(6) position (m(6)A) is a dynamic and abundant epitranscriptomic mark that regulates critical aspects of eukaryotic RNA metabolism in numerous biological processes. The RNA methyltransferases METTL3 and METTL14 are components of a multisubunit m(6)A writer complex whose enzymatic activity is substantially higher than the activities of METTL3 or METTL14 alone. The molecular mechanism underpinning this synergistic effect is poorly understood. Here we report the crystal structure of the catalytic core of the human m(6)A writer complex comprising METTL3 and METTL14. The structure reveals the heterodimeric architecture of the complex and donor substrate binding by METTL3. Structure-guided mutagenesis indicates that METTL3 is the catalytic subunit of the complex, whereas METTL14 has a degenerate active site and plays non-catalytic roles in maintaining complex integrity and substrate RNA binding. These studies illuminate the molecular mechanism and evolutionary history of eukaryotic m(6)A modification in post-transcriptional genome regulation.

The nucleolus is a membrane-less organelle formed through liquid-liquid phase separation of its components from the surrounding nucleoplasm. Here, we show that nucleophosmin (NPM1) integrates within the nucleolus via a multi-modal mechanism involving multivalent interactions with proteins containing arginine-rich linear motifs (R-motifs) and ribosomal RNA (rRNA). Importantly, these R-motifs are found in canonical nucleolar localization signals. Based on a novel combination of biophysical approaches, we propose a model for the molecular organization within liquid-like droplets formed by the N-terminal domain of NPM1 and R-motif peptides, thus providing insights into the structural organization of the nucleolus. We identify multivalency of acidic tracts and folded nucleic acid binding domains, mediated by N-terminal domain oligomerization, as structural features required for phase separation of NPM1 with other nucleolar components in vitro and for localization within mammalian nucleoli. We propose that one mechanism of nucleolar localization involves phase separation of proteins…

Alternative splicing (AS) can critically affect gene function and disease, yet mapping splicing variations remains a challenge. Here, we propose a new approach to define and quantify mRNA splicing in units of local splicing variations (LSVs). LSVs capture previously defined types of alternative splicing as well as more complex transcript variations. Building the first genome wide map of LSVs from twelve mouse tissues, we find complex LSVs constitute over 30% of tissue dependent transcript variations and affect specific protein families. We show the prevalence of complex LSVs is conserved in humans and identify hundreds of LSVs that are specific to brain subregions or altered in Alzheimer's patients. Amongst those are novel isoforms in the Camk2 family and a novel poison exon in Ptbp1, a key splice factor in neurogenesis. We anticipate the approach presented here will advance the ability to relate tissue-specific splice variation to genetic variation, phenotype, and disease.

Cells can enter into a dormant state when faced with unfavorable conditions. However, how cells enter into and recover from this state is still poorly understood. Here, we study dormancy in different eukaryotic organisms and find it to be associated with a significant decrease in the mobility of organelles and foreign tracer particles. We show that this reduced mobility is caused by an influx of protons and a marked acidification of the cytoplasm, which leads to widespread macromolecular assembly of proteins and triggers a transition of the cytoplasm to a solid-like state with increased mechanical stability. We further demonstrate that this transition is required for cellular survival under conditions of starvation. Our findings have broad implications for understanding alternative physiological states, such as quiescence and dormancy, and create a new view of the cytoplasm as an adaptable fluid that can reversibly transition into a protective solid-like state.

The FDA approved drug rapamycin increases lifespan in rodents and delays age-related dysfunction in rodents and humans. Nevertheless, important questions remain regarding the optimal dose, duration, and mechanisms of action in the context of healthy aging. Here we show that 3 months of rapamycin treatment is sufficient to increase life expectancy by up to 60% and improve measures of healthspan in middle-aged mice. This transient treatment is also associated with a remodeling of the microbiome, including dramatically increased prevalence of segmented filamentous bacteria in the small intestine. We also define a dose in female mice that does not extend lifespan, but is associated with a striking shift in cancer prevalence toward aggressive hematopoietic cancers and away from non-hematopoietic malignancies. These data suggest that a short-term rapamycin treatment late in life has persistent effects that can robustly delay aging, influence cancer prevalence, and modulate the microbiome.

A process of global importance in carbon cycling is the remineralization of algae biomass by heterotrophic bacteria, most notably during massive marine algae blooms. Such blooms can trigger secondary blooms of planktonic bacteria that consist of swift successions of distinct bacterial clades, most prominently members of the Flavobacteriia, Gammaproteobacteria and the alphaproteobacterial Roseobacter clade. We investigated such successions during spring phytoplankton blooms in the southern North Sea (German Bight) for four consecutive years. Dense sampling and high-resolution taxonomic analyses allowed the detection of recurring patterns down to the genus level. Metagenome analyses also revealed recurrent patterns at the functional level, in particular with respect to algal polysaccharide degradation genes. We, therefore, hypothesize that even though there is substantial inter-annual variation between spring phytoplankton blooms, the accompanying succession of bacterial clades is largely governed by deterministic principles such as substrate-induced forcing.

Macrophage activation/polarization to distinct functional states is critically supported by metabolic shifts. How polarizing signals coordinate metabolic and functional reprogramming, and the potential implications for control of macrophage activation, remains poorly understood. Here we show that IL-4 signaling co-opts the Akt-mTORC1 pathway to regulate Acly, a key enzyme in Ac-CoA synthesis, leading to increased histone acetylation and M2 gene induction. Only a subset of M2 genes is controlled in this way, including those regulating cellular proliferation and chemokine production. Moreover, metabolic signals impinge on the Akt-mTORC1 axis for such control of M2 activation. We propose that Akt-mTORC1 signaling calibrates metabolic state to energetically demanding aspects of M2 activation, which may define a new role for metabolism in supporting macrophage activation.

Cryo-EM has revealed the structures of many challenging yet exciting macromolecular assemblies at near-atomic resolution (3-4.5Å), providing biological phenomena with molecular descriptions. However, at these resolutions, accurately positioning individual atoms remains challenging and error-prone. Manually refining thousands of amino acids - typical in a macromolecular assembly - is tedious and time-consuming. We present an automated method that can improve the atomic details in models that are manually built in near-atomic-resolution cryo-EM maps. Applying the method to three systems recently solved by cryo-EM, we are able to improve model geometry while maintaining the fit-to-density. Backbone placement errors are automatically detected and corrected, and the refinement shows a large radius of convergence. The results demonstrate that the method is amenable to structures with symmetry, of very large size, and containing RNA as well as covalently bound ligands. The method should streamline the cryo-EM structure determination process, providing accurate and unbiased atomic structure…

Mitochondrial dysfunction is associated with a spectrum of human disorders, ranging from rare, inborn errors of metabolism to common, age-associated diseases such as neurodegeneration. How these lesions give rise to diverse pathology is not well understood, partly because their proximal consequences have not been well-studied in mammalian cells. Here we provide two lines of evidence that mitochondrial respiratory chain dysfunction leads to alterations in one-carbon metabolism pathways. First, using hypothesis-generating metabolic, proteomic, and transcriptional profiling, followed by confirmatory experiments, we report that mitochondrial DNA depletion leads to an ATF4-mediated increase in serine biosynthesis and transsulfuration. Second, we show that lesioning the respiratory chain impairs mitochondrial production of formate from serine, and that in some cells, respiratory chain inhibition leads to growth defects upon serine withdrawal that are rescuable with purine or formate supplementation. Our work underscores the connection between the respiratory chain and one-carbon…

Clarifying gene expression in narrowly defined neuronal populations can provide insight into cellular identity, computation, and functionality. Here, we used next-generation RNA sequencing (RNA-seq) to produce a quantitative, whole genome characterization of gene expression for the major excitatory neuronal classes of the hippocampus; namely, granule cells and mossy cells of the dentate gyrus, and pyramidal cells of areas CA3, CA2, and CA1. Moreover, for the canonical cell classes of the trisynaptic loop, we profiled transcriptomes at both dorsal and ventral poles, producing a cell-class- and region-specific transcriptional description for these populations. This dataset clarifies the transcriptional properties and identities of lesser-known cell classes, and moreover reveals unexpected variation in the trisynaptic loop across the dorsal-ventral axis. We have created a public resource, Hipposeq (http://hipposeq.janelia.org), which provides analysis and visualization of these data and will act as a roadmap relating molecules to cells, circuits, and computation in the hippocampus.

The prokaryotic CRISPR (clustered regularly interspaced palindromic repeats)-associated protein, Cas9, has been widely adopted as a tool for editing, imaging, and regulating eukaryotic genomes. However, our understanding of how to select single-guide RNAs (sgRNAs) that mediate efficient Cas9 activity is incomplete, as we lack insight into how chromatin impacts Cas9 targeting. To address this gap, we analyzed large-scale genetic screens performed in human cell lines using either nuclease-active or nuclease-dead Cas9 (dCas9). We observed that highly active sgRNAs for Cas9 and dCas9 were found almost exclusively in regions of low nucleosome occupancy. In vitro experiments demonstrated that nucleosomes in fact directly impede Cas9 binding and cleavage, while chromatin remodeling can restore Cas9 access. Our results reveal a critical role of eukaryotic chromatin in dictating the targeting specificity of this transplanted bacterial enzyme, and provide rules for selecting Cas9 target sites distinct from and complementary to those based on sequence properties.

The structure of axonal arbors controls how signals from individual neurons are routed within the mammalian brain. However, the arbors of very few long-range projection neurons have been reconstructed in their entirety, as axons with diameters as small as 100 nm arborize in target regions dispersed over many millimeters of tissue. We introduce a platform for high-resolution, three-dimensional fluorescence imaging of complete tissue volumes that enables the visualization and reconstruction of long-range axonal arbors. This platform relies on a high-speed two-photon microscope integrated with a tissue vibratome and a suite of computational tools for large-scale image data. We demonstrate the power of this approach by reconstructing the axonal arbors of multiple neurons in the motor cortex across a single mouse brain.

and human cells, we show here that seipin, an ER protein implicated in LD biology, mediates a discrete step in LD formation-the conversion of small, nascent LDs to larger, mature LDs. Seipin forms discrete and dynamic foci in the ER that interact with nascent LDs to enable their growth. In the absence of seipin, numerous small, nascent LDs accumulate near the ER and most often fail to grow. Those that do grow prematurely acquire lipid synthesis enzymes and undergo expansion, eventually leading to the giant LDs characteristic of seipin deficiency. Our studies identify a discrete step of LD formation, namely the conversion of nascent LDs to mature LDs, and define a molecular role for seipin in this process, most likely by acting at ER-LD contact sites to enable lipid transfer to nascent LDs.

Inducible epigenetic changes in eukaryotes are believed to enable rapid adaptation to environmental fluctuations. We have found distinct regions of the Arabidopsis genome that are susceptible to DNA (de)methylation in response to hyperosmotic stress. The stress-induced epigenetic changes are associated with conditionally heritable adaptive phenotypic stress responses. However, these stress responses are primarily transmitted to the next generation through the female lineage due to widespread DNA glycosylase activity in the male germline, and extensively reset in the absence of stress. Using the CNI1/ATL31 locus as an example, we demonstrate that epigenetically targeted sequences function as distantly-acting control elements of antisense long non-coding RNAs, which in turn regulate targeted gene expression in response to stress. Collectively, our findings reveal that plants use a highly dynamic maternal 'short-term stress memory' with which to respond to adverse external conditions. This transient memory relies on the DNA methylation machinery and associated transcriptional changes…

Zika virus was discovered in Uganda in 1947 and is transmitted by Aedes mosquitoes, which also act as vectors for dengue and chikungunya viruses throughout much of the tropical world. In 2007, an outbreak in the Federated States of Micronesia sparked public health concern. In 2013, the virus began to spread across other parts of Oceania and in 2015, a large outbreak in Latin America began in Brazil. Possible associations with microcephaly and Guillain-Barré syndrome observed in this outbreak have raised concerns about continued global spread of Zika virus, prompting its declaration as a Public Health Emergency of International Concern by the World Health Organization. We conducted species distribution modelling to map environmental suitability for Zika. We show a large portion of tropical and sub-tropical regions globally have suitable environmental conditions with over 2.17 billion people inhabiting these areas.

Antigen Ki-67 is a nuclear protein expressed in proliferating mammalian cells. It is widely used in cancer histopathology but its functions remain unclear. Here, we show that Ki-67 controls heterochromatin organisation. Altering Ki-67 expression levels did not significantly affect cell proliferation in vivo. Ki-67 mutant mice developed normally and cells lacking Ki-67 proliferated efficiently. Conversely, upregulation of Ki-67 expression in differentiated tissues did not prevent cell cycle arrest. Ki-67 interactors included proteins involved in nucleolar processes and chromatin regulators. Ki-67 depletion disrupted nucleologenesis but did not inhibit pre-rRNA processing. In contrast, it altered gene expression. Ki-67 silencing also had wide-ranging effects on chromatin organisation, disrupting heterochromatin compaction and long-range genomic interactions. Trimethylation of histone H3K9 and H4K20 was relocalised within the nucleus. Finally, overexpression of human or Xenopus Ki-67 induced ectopic heterochromatin formation. Altogether, our results suggest that Ki-67 expression in…

Although it is clear that trisomy 21 causes Down syndrome, the molecular events acting downstream of the trisomy remain ill defined. Using complementary genomics analyses, we identified the interferon pathway as the major signaling cascade consistently activated by trisomy 21 in human cells. Transcriptome analysis revealed that trisomy 21 activates the interferon transcriptional response in fibroblast and lymphoblastoid cell lines, as well as circulating monocytes and T cells. Trisomy 21 cells show increased induction of interferon-stimulated genes and decreased expression of ribosomal proteins and translation factors. An shRNA screen determined that the interferon-activated kinases JAK1 and TYK2 suppress proliferation of trisomy 21 fibroblasts, and this defect is rescued by pharmacological JAK inhibition. Therefore, we propose that interferon activation, likely via increased gene dosage of the four interferon receptors encoded on chromosome 21, contributes to many of the clinical impacts of trisomy 21, and that interferon antagonists could have therapeutic benefits.

Viruses interact with hundreds to thousands of proteins in mammals, yet adaptation against viruses has only been studied in a few proteins specialized in antiviral defense. Whether adaptation to viruses typically involves only specialized antiviral proteins or affects a broad array of virus-interacting proteins is unknown. Here, we analyze adaptation in ~1300 virus-interacting proteins manually curated from a set of 9900 proteins conserved in all sequenced mammalian genomes. We show that viruses (i) use the more evolutionarily constrained proteins within the cellular functions they interact with and that (ii) despite this high constraint, virus-interacting proteins account for a high proportion of all protein adaptation in humans and other mammals. Adaptation is elevated in virus-interacting proteins across all functional categories, including both immune and non-immune functions. We conservatively estimate that viruses have driven close to 30% of all adaptive amino acid changes in the part of the human proteome conserved within mammals. Our results suggest that viruses are one of…

The Drosophila genome contains >13000 protein-coding genes, the majority of which remain poorly investigated. Important reasons include the lack of antibodies or reporter constructs to visualise these proteins. Here, we present a genome-wide fosmid library of 10000 GFP-tagged clones, comprising tagged genes and most of their regulatory information. For 880 tagged proteins, we created transgenic lines, and for a total of 207 lines, we assessed protein expression and localisation in ovaries, embryos, pupae or adults by stainings and live imaging approaches. Importantly, we visualised many proteins at endogenous expression levels and found a large fraction of them localising to subcellular compartments. By applying genetic complementation tests, we estimate that about two-thirds of the tagged proteins are functional. Moreover, these tagged proteins enable interaction proteomics from developing pupae and adult flies. Taken together, this resource will boost systematic analysis of protein expression and localisation in various cellular and developmental contexts.

Asia is considered an important source of influenza A virus (IAV) pandemics, owing to large, diverse viral reservoirs in poultry and swine. However, the zoonotic origins of the 2009 A/H1N1 influenza pandemic virus (pdmH1N1) remain unclear, due to conflicting evidence from swine and humans. There is strong evidence that the first human outbreak of pdmH1N1 occurred in Mexico in early 2009. However, no related swine viruses have been detected in Mexico or any part of the Americas, and to date the most closely related ancestor viruses were identified in Asian swine. Here, we use 58 new whole-genome sequences from IAVs collected in Mexican swine to establish that the swine virus responsible for the 2009 pandemic evolved in central Mexico. This finding highlights how the 2009 pandemic arose from a region not considered a pandemic risk, owing to an expansion of IAV diversity in swine resulting from long-distance live swine trade.

A major constraint on the evolution of large body sizes in animals is an increased risk of developing cancer. There is no correlation, however, between body size and cancer risk. This lack of correlation is often referred to as 'Peto's Paradox'. Here, we show that the elephant genome encodes 20 copies of the tumor suppressor gene TP53 and that the increase in TP53 copy number occurred coincident with the evolution of large body sizes, the evolution of extreme sensitivity to genotoxic stress, and a hyperactive TP53 signaling pathway in the elephant (Proboscidean) lineage. Furthermore, we show that several of the TP53 retrogenes (TP53RTGs) are transcribed and likely translated. While TP53RTGs do not appear to directly function as transcription factors, they do contribute to the enhanced sensitivity of elephant cells to DNA damage and the induction of apoptosis by regulating activity of the TP53 signaling pathway. These results suggest that an increase in the copy number of TP53 may have played a direct role in the evolution of very large body sizes and the resolution of Peto's paradox…

A defining feature of mitochondria is their maternal mode of inheritance. However, little is understood about the cellular mechanism through which paternal mitochondria, delivered from sperm, are eliminated from early mammalian embryos. Autophagy has been implicated in nematodes, but whether this mechanism is conserved in mammals has been disputed. Here, we show that cultured mouse fibroblasts and pre-implantation embryos use a common pathway for elimination of mitochondria. Both situations utilize mitophagy, in which mitochondria are sequestered by autophagosomes and delivered to lysosomes for degradation. The E3 ubiquitin ligases PARKIN and MUL1 play redundant roles in elimination of paternal mitochondria. The process is associated with depolarization of paternal mitochondria and additionally requires the mitochondrial outer membrane protein FIS1, the autophagy adaptor P62, and PINK1 kinase. Our results indicate that strict maternal transmission of mitochondria relies on mitophagy and uncover a collaboration between MUL1 and PARKIN in this process.

Transposable elements (TEs) are powerful motors of genome evolution yet a comprehensive assessment of recent transposition activity at the species level is lacking for most organisms. Here, using genome sequencing data for 211 Arabidopsis thaliana accessions taken from across the globe, we identify thousands of recent transposition events involving half of the 326 TE families annotated in this plant species. We further show that the composition and activity of the 'mobilome' vary extensively between accessions in relation to climate and genetic factors. Moreover, TEs insert equally throughout the genome and are rapidly purged by natural selection from gene-rich regions because they frequently affect genes, in multiple ways. Remarkably, loci controlling adaptive responses to the environment are the most frequent transposition targets observed. These findings demonstrate the pervasive, species-wide impact that a rich mobilome can have and the importance of transposition as a recurrent generator of large-effect alleles.

Restoration of touch after hand amputation is a desirable feature of ideal prostheses. Here, we show that texture discrimination can be artificially provided in human subjects by implementing a neuromorphic real-time mechano-neuro-transduction (MNT), which emulates to some extent the firing dynamics of SA1 cutaneous afferents. The MNT process was used to modulate the temporal pattern of electrical spikes delivered to the human median nerve via percutaneous microstimulation in four intact subjects and via implanted intrafascicular stimulation in one transradial amputee. Both approaches allowed the subjects to reliably discriminate spatial coarseness of surfaces as confirmed also by a hybrid neural model of the median nerve. Moreover, MNT-evoked EEG activity showed physiologically plausible responses that were superimposable in time and topography to the ones elicited by a natural mechanical tactile stimulation. These findings can open up novel opportunities for sensory restoration in the next generation of neuro-prosthetic hands.

Neuronal circuit mapping using electron microscopy demands laborious proofreading or reconciliation of multiple independent reconstructions. Here, we describe new methods to apply quantitative arbor and network context to iteratively proofread and reconstruct circuits and create anatomically enriched wiring diagrams. We measured the morphological underpinnings of connectivity in new and existing reconstructions of Drosophila sensorimotor (larva) and visual (adult) systems. Synaptic inputs were preferentially located on numerous small, microtubule-free 'twigs' which branch off a single microtubule-containing 'backbone'. Omission of individual twigs accounted for 96% of errors. However, the synapses of highly connected neurons were distributed across multiple twigs. Thus, the robustness of a strong connection to detailed twig anatomy was associated with robustness to reconstruction error. By comparing iterative reconstruction to the consensus of multiple reconstructions, we show that our method overcomes the need for redundant effort through the discovery and application of…

In vitro selection of antibodies allows to obtain highly functional binders, rapidly and at lower cost. Here, we describe the first fully synthetic phage display library of humanized llama single domain antibody (NaLi-H1: Nanobody Library Humanized 1). Based on a humanized synthetic single domain antibody (hs2dAb) scaffold optimized for intracellular stability, the highly diverse library provides high affinity binders without animal immunization. NaLi-H1 was screened following several selection schemes against various targets (Fluorescent proteins, actin, tubulin, p53, HP1). Conformation antibodies against active RHO GTPase were also obtained. Selected hs2dAb were used in various immunoassays and were often found to be functional intrabodies, enabling tracking or inhibition of endogenous targets. Functionalization of intrabodies allowed specific protein knockdown in living cells. Finally, direct selection against the surface of tumor cells produced hs2dAb directed against tumor-specific antigens further highlighting the potential use of this library for therapeutic applications.

Demyelination in MS disrupts nerve signals and contributes to axon degeneration. While remyelination promises to restore lost function, it remains unclear whether remyelination will prevent axonal loss. Inflammatory demyelination is accompanied by significant neuronal loss in the experimental autoimmune encephalomyelitis (EAE) mouse model and evidence for remyelination in this model is complicated by ongoing inflammation, degeneration and possible remyelination. Demonstrating the functional significance of remyelination necessitates selectively altering the timing of remyelination relative to inflammation and degeneration. We demonstrate accelerated remyelination after EAE induction by direct lineage analysis and hypothesize that newly formed myelin remains stable at the height of inflammation due in part to the absence of MOG expression in immature myelin. Oligodendroglial-specific genetic ablation of the M1 muscarinic receptor, a potent negative regulator of oligodendrocyte differentiation and myelination, results in accelerated remyelination, preventing axonal loss and improving…

Somatic mutations in the estrogen receptor alpha (ERα) gene (ESR1), especially Y537S and D538G, have been linked to acquired resistance to endocrine therapies. Cell-based studies demonstrated that these mutants confer ERα constitutive activity and antiestrogen resistance and suggest that ligand-binding domain dysfunction leads to endocrine therapy resistance. Here, we integrate biophysical and structural biology data to reveal how these mutations lead to a constitutively active and antiestrogen-resistant ERα. We show that these mutant ERs recruit coactivator in the absence of hormone while their affinities for estrogen agonist (estradiol) and antagonist (4-hydroxytamoxifen) are reduced. Further, they confer antiestrogen resistance by altering the conformational dynamics of the loop connecting Helix 11 and Helix 12 in the ligand-binding domain of ERα, which leads to a stabilized agonist state and an altered antagonist state that resists inhibition.

Little is known about the excess mortality caused by multidrug-resistant (MDR) bacterial infection in low- and middle-income countries (LMICs). We retrospectively obtained microbiology laboratory and hospital databases of nine public hospitals in northeast Thailand from 2004 to 2010, and linked these with the national death registry to obtain the 30-day mortality outcome. The 30-day mortality in those with MDR community-acquired bacteraemia, healthcare-associated bacteraemia, and hospital-acquired bacteraemia were 35% (549/1555), 49% (247/500), and 53% (640/1198), respectively. We estimate that 19,122 of 45,209 (43%) deaths in patients with hospital-acquired infection due to MDR bacteria in Thailand in 2010 represented excess mortality caused by MDR. We demonstrate that national statistics on the epidemiology and burden of MDR in LMICs could be improved by integrating information from readily available databases. The prevalence and mortality attributable to MDR in Thailand are high. This is likely to reflect the situation in other LMICs.

The mechanism by which chromatids and chromosomes are segregated during mitosis and meiosis is a major puzzle of biology and biophysics. Using polymer simulations of chromosome dynamics, we show that a single mechanism of loop extrusion by condensins can robustly compact, segregate and disentangle chromosomes, arriving at individualized chromatids with morphology observed in vivo. Our model resolves the paradox of topological simplification concomitant with chromosome 'condensation', and explains how enzymes a few nanometers in size are able to control chromosome geometry and topology at micron length scales. We suggest that loop extrusion is a universal mechanism of genome folding that mediates functional interactions during interphase and compacts chromosomes during mitosis.

Long lasting pyrethroid treated bednets are the most important tool for preventing malaria. Pyrethroid resistant Anopheline mosquitoes are now ubiquitous in Africa, though the public health impact remains unclear, impeding the deployment of more expensive nets. Meta-analyses of bioassay studies and experimental hut trials are used to characterise how pyrethroid resistance changes the efficacy of standard bednets, and those containing the synergist piperonyl butoxide (PBO), and assess its impact on malaria control. New bednets provide substantial personal protection until high levels of resistance, though protection may wane faster against more resistant mosquito populations as nets age. Transmission dynamics models indicate that even low levels of resistance would increase the incidence of malaria due to reduced mosquito mortality and lower overall community protection over the life-time of the net. Switching to PBO bednets could avert up to 0.5 clinical cases per person per year in some resistance scenarios.

In the absence of salient sensory cues to guide behavior, animals must still execute sequences of motor actions in order to forage and explore. How such successive motor actions are coordinated to form global locomotion trajectories is unknown. We mapped the structure of larval zebrafish swim trajectories in homogeneous environments and found that trajectories were characterized by alternating sequences of repeated turns to the left and to the right. Using whole-brain light-sheet imaging, we identified activity relating to the behavior in specific neural populations that we termed the anterior rhombencephalic turning region (ARTR). ARTR perturbations biased swim direction and reduced the dependence of turn direction on turn history, indicating that the ARTR is part of a network generating the temporal correlations in turn direction. We also find suggestive evidence for ARTR mutual inhibition and ARTR projections to premotor neurons. Finally, simulations suggest the observed turn sequences may underlie efficient exploration of local environments.

Gene-environment interactions impact the development of neuropsychiatric disorders, but the relative contributions are unclear. Here, we identify gut microbiota as sufficient to induce depressive-like behaviors in genetically distinct mouse strains. Daily gavage of vehicle (dH2O) in nonobese diabetic (NOD) mice induced a social avoidance behavior that was not observed in C57BL/6 mice. This was not observed in NOD animals with depleted microbiota via oral administration of antibiotics. Transfer of intestinal microbiota, including members of the Clostridiales, Lachnospiraceae and Ruminococcaceae, from vehicle-gavaged NOD donors to microbiota-depleted C57BL/6 recipients was sufficient to induce social avoidance and change gene expression and myelination in the prefrontal cortex. Metabolomic analysis identified increased cresol levels in these mice, and exposure of cultured oligodendrocytes to this metabolite prevented myelin gene expression and differentiation. Our results thus demonstrate that the gut microbiota modifies the synthesis of key metabolites affecting gene expression in…