Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells.

It has long been the dream of biologists to map gene expression at the single-cell level. With such data one might track heterogeneous cell sub-populations, and infer regulatory relationships between genes and pathways. Recently, RNA sequencing has achieved single-cell resolution. What is limiting is an effective way to routinely isolate and process large numbers of individual cells for quantitative in-depth sequencing. We have developed a high-throughput droplet-microfluidic approach for barcoding the RNA from thousands of individual cells for subsequent analysis by next-generation sequencing. The method shows a surprisingly low noise profile and is readily adaptable to other sequencing-based assays. We analyzed mouse embryonic stem cells, revealing in detail the population structure and the heterogeneous onset of differentiation after leukemia inhibitory factor (LIF) withdrawal. The reproducibility of these high-throughput single-cell data allowed us to deconstruct cell populations and infer gene expression relationships. VIDEO ABSTRACT.

A noncanonical Frizzled2 pathway regulates epithelial-mesenchymal transition and metastasis.

Wnt signaling plays a critical role in embryonic development, and genetic aberrations in this network have been broadly implicated in colorectal cancer. We find that the Wnt receptor Frizzled2 (Fzd2) and its ligands Wnt5a/b are elevated in metastatic liver, lung, colon, and breast cancer cell lines and in high-grade tumors and that their expression correlates with markers of epithelial-mesenchymal transition (EMT). Pharmacologic and genetic perturbations reveal that Fzd2 drives EMT and cell migration through a previously unrecognized, noncanonical pathway that includes Fyn and Stat3. A gene signature regulated by this pathway predicts metastasis and overall survival in patients. We have developed an antibody to Fzd2 that reduces cell migration and invasion and inhibits tumor growth and metastasis in xenografts. We propose that targeting this pathway could provide benefit for patients with tumors expressing high levels of Fzd2 and Wnt5a/b.

A time-resolved single-cell roadmap of the logic driving anterior neural crest diversification from neural border to migration stages.

Neural crest cells exemplify cellular diversification from a multipotent progenitor population. However, the full sequence of early molecular choices orchestrating the emergence of neural crest heterogeneity from the embryonic ectoderm remains elusive. Gene-regulatory-networks (GRN) govern early development and cell specification toward definitive neural crest. Here, we combine ultradense single-cell transcriptomes with machine-learning and large-scale transcriptomic and epigenomic experimental validation of selected trajectories, to provide the general principles and highlight specific features of the GRN underlying neural crest fate diversification from induction to early migration stages using Xenopus frog embryos as a model. During gastrulation, a transient neural border zone state precedes the choice between neural crest and placodes which includes multiple converging gene programs. During neurulation, transcription factor connectome, and bifurcation analyses demonstrate the early emergence of neural crest fates at the neural plate stage, alongside an unbiased multipotent-like lineage persisting until epithelial-mesenchymal transition stage. We also decipher circuits driving cranial and vagal neural crest formation and provide a broadly applicable high-throughput validation strategy for investigating single-cell transcriptomes in vertebrate GRNs in development, evolution, and disease.

Open Genes-a new comprehensive database of human genes associated with aging and longevity.

The Open Genes database was created to enhance and simplify the search for potential aging therapy targets. We collected data on 2402 genes associated with aging and developed convenient tools for searching and comparing gene features. A comprehensive description of genes has been provided, including lifespan-extending interventions, age-related changes, longevity associations, gene evolution, associations with diseases and hallmarks of aging, and functions of gene products. For each experiment, we presented the necessary structured data for evaluating the experiment's quality and interpreting the study's findings. Our goal was to stay objective and precise while connecting a particular gene to human aging. We distinguished six types of studies and 12 criteria for adding genes to our database. Genes were classified according to the confidence level of the link between the gene and aging. All the data collected in a database are provided both by an API and a user interface. The database is publicly available on a website at https://open-genes.org/.

A new atlas to study embryonic cell types in Xenopus.

This paper introduces a single-cell atlas for pivotal developmental stages in Xenopus, encompassing gastrulation, neurulation, and early tailbud. Notably surpassing its predecessors, the new atlas enhances gene mapping, read counts, and gene/cell type nomenclature. Leveraging the latest Xenopus tropicalis genome version, alongside advanced alignment pipelines and machine learning for cell type assignment, this release maintains consistency with previous cell type annotations while rectifying nomenclature issues. Employing an unbiased approach for cell type assignment proves especially apt for embryonic contexts, given the considerable number of non-terminally differentiated cell types. An alternative cell type attribution here adopts a fuzzy, non-deterministic stance, capturing the transient nature of early embryo progenitor cells by presenting an ensemble of types in superposition. The value of the new resource is emphasized through numerous examples, with a focus on previously unexplored germ cell populations where we uncover novel transcription onset features. Offering interactive exploration via a user-friendly web portal and facilitating complete data downloads, this atlas serves as a comprehensive and accessible reference.

Computationally enhanced quantitative phase microscopy reveals autonomous oscillations in mammalian cell growth.

The fine balance of growth and division is a fundamental property of the physiology of cells, and one of the least understood. Its study has been thwarted by difficulties in the accurate measurement of cell size and the even greater challenges of measuring growth of a single cell over time. We address these limitations by demonstrating a computationally enhanced methodology for quantitative phase microscopy for adherent cells, using improved image processing algorithms and automated cell-tracking software. Accuracy has been improved more than twofold and this improvement is sufficient to establish the dynamics of cell growth and adherence to simple growth laws. It is also sufficient to reveal unknown features of cell growth, previously unmeasurable. With these methodological and analytical improvements, in several cell lines we document a remarkable oscillation in growth rate, occurring throughout the cell cycle, coupled to cell division or birth yet independent of cell cycle progression. We expect that further exploration with this advanced tool will provide a better understanding of growth rate regulation in mammalian cells.

Inverse Lansing Effect: Maternal Age and Provisioning Affecting Daughters' Longevity and Male Offspring Production.

AbstractMaternal age effects on offspring life history are known in a variety of organisms, with offspring of older mothers typically having lower life expectancy (the Lansing effect). However, there is no consensus on the generality and mechanisms of this pattern. We tested predictions of the Lansing effect in several Daphnia magna clones and observed clone-specific magnitude and direction of the maternal age effect on offspring longevity. We also report ambidirectional, genotype-specific effects of maternal age on the propensity of daughters to produce male offspring. Focusing on two clones with contrasting life histories, we demonstrate that maternal age effects can be explained by lipid provisioning of embryos by mothers of different ages. Individuals from a single-generation maternal age reversal treatment showed intermediate life span and intermediate lipid content at birth. In the clone characterized by the "inverse Lansing effect," neonates produced by older mothers showed higher mitochondrial membrane potential in neural tissues than their counterparts born to younger mothers. We conclude that an inverse Lansing effect is possible and hypothesize that it may be caused by age-specific maternal lipid provisioning creating a calorically restricted environment during embryonic development, which in turn reduces fecundity and increases life span in offspring.

Lack of age-related respiratory changes in Daphnia.

Aging is a multifaceted process of accumulation of damage and waste in cells and tissues; age-related changes in mitochondria and in respiratory metabolism have the focus of aging research for decades. Studies of aging in nematodes, flies and mammals all revealed age-related decline in respiratory functions, with somewhat controversial causative role. Here we investigated age-related changes in respiration rates, lactate/pyruvate ratio, a commonly used proxy for NADH/NAD+ balance, and mitochondrial membrane potential in 4 genotypes of an emerging model organism for aging research, a cyclic parthenogen Daphnia magna. We show that total body weight-adjusted respiration rate decreased with age, although this decrease was small in magnitude and could be fully accounted for by the decrease in locomotion and feeding activity. Neither total respiration normalized by protein content, nor basal respiration rate measured in anaesthetized animals decreased with age. Lactate/pyruvate ratio and mitochondrial membrane potential (∆Ψmt) showed no age-related changes, with possible exceptions of ∆Ψmt in epipodites (excretory and gas exchange organs) in which ∆Ψmt decreased with age and in the optical lobe of the brain, in which ∆Ψmt showed a maximum at middle age. We conclude that actuarial senescence in Daphnia is not caused by a decline in respiratory metabolism and discuss possible mechanisms of maintaining mitochondrial healthspan throughout the lifespan.

Hemichordate genomes and deuterostome origins.

Acorn worms, also known as enteropneust (literally, 'gut-breathing') hemichordates, are marine invertebrates that share features with echinoderms and chordates. Together, these three phyla comprise the deuterostomes. Here we report the draft genome sequences of two acorn worms, Saccoglossus kowalevskii and Ptychodera flava. By comparing them with diverse bilaterian genomes, we identify shared traits that were probably inherited from the last common deuterostome ancestor, and then explore evolutionary trajectories leading from this ancestor to hemichordates, echinoderms and chordates. The hemichordate genomes exhibit extensive conserved synteny with amphioxus and other bilaterians, and deeply conserved non-coding sequences that are candidates for conserved gene-regulatory elements. Notably, hemichordates possess a deuterostome-specific genomic cluster of four ordered transcription factor genes, the expression of which is associated with the development of pharyngeal 'gill' slits, the foremost morphological innovation of early deuterostomes, and is probably central to their filter-feeding lifestyle. Comparative analysis reveals numerous deuterostome-specific gene novelties, including genes found in deuterostomes and marine microbes, but not other animals. The putative functions of these genes can be linked to physiological, metabolic and developmental specializations of the filter-feeding ancestor.

Bayesian Confidence Intervals for Multiplexed Proteomics Integrate Ion-statistics with Peptide Quantification Concordance.

Multiplexed proteomics has emerged as a powerful tool to measure relative protein expression levels across multiple conditions. The relative protein abundances are inferred by comparing the signals generated by isobaric tags, which encode the samples' origins. Intuitively, the trust associated with a protein measurement depends on the similarity of ratios from the protein's peptides and the signal-strength of these measurements. However, typically the average peptide ratio is reported as the estimate of relative protein abundance, which is only the most likely ratio with a very naive model. Moreover, there is no sense on the confidence in these measurements. Here, we present a mathematically rigorous approach that integrates peptide signal strengths and peptide-measurement agreement into an estimation of the true protein ratio and the associated confidence (BACIQ). The main advantages of BACIQ are: (1) It removes the need to threshold reported peptide signal based on an arbitrary cut-off, thereby reporting more measurements from a given experiment; (2) Confidence can be assigned without replicates; (3) For repeated experiments BACIQ provides confidence intervals for the union, not the intersection, of quantified proteins; (4) For repeated experiments, BACIQ confidence intervals are more predictive than confidence intervals based on protein measurement agreement. To demonstrate the power of BACIQ we reanalyzed previously published data on subcellular protein movement on treatment with an Exportin-1 inhibiting drug. We detect ∼2× more highly significant movers, down to subcellular localization changes of ∼1%. Thus, our method drastically increases the value obtainable from quantitative proteomics experiments, helping researchers to interpret their data and prioritize resources. To make our approach easily accessible we distribute it via a Python/Stan package.

A cell type annotation Jamboree-Revival of а communal science forum.

Cell Atlases are currently being constructed for human tissues as well as several model organisms. New technologies make creation of vast datasets in many species possible, but the value of such data crucially depends on the quality of annotation. The tools of annotating single cell data and creating knowledge representations comparable across organisms have been lagging. We argue that successfully creating Cell Atlases will require a revival of a boot-camp style forum for communal annotation combined with an intensive learning workshop, dubbed a "Jamboree". We report on our experience of successfully developing a structure and curriculum and running such a Jamboree for Xenopus Embryonic Cell Types at the Janelia Farms campus of the Howard Hughes Medical Institute.

Proteomics of phosphorylation and protein dynamics during fertilization and meiotic exit in the Xenopus egg.

Fertilization releases the meiotic arrest and initiates the events that prepare the egg for the ensuing developmental program. Protein degradation and phosphorylation are known to regulate protein activity during this process. However, the full extent of protein loss and phosphoregulation is still unknown. We examined absolute protein and phosphosite dynamics of the fertilization response by mass spectrometry-based proteomics in electroactivated eggs. To do this, we developed an approach for calculating the stoichiometry of phosphosites from multiplexed proteomics that is compatible with dynamic, stable, and multisite phosphorylation. Overall, the data suggest that degradation is limited to a few low-abundance proteins. However, this degradation promotes extensive dephosphorylation that occurs over a wide range of abundances during meiotic exit. We also show that eggs release a large amount of protein into the medium just after fertilization, most likely related to the blocks to polyspermy. Concomitantly, there is a substantial increase in phosphorylation likely tied to calcium-activated kinases. We identify putative degradation targets and components of the slow block to polyspermy. The analytical approaches demonstrated here are broadly applicable to studies of dynamic biological systems.

Quantitative Proteomics Reveals Remodeling of Protein Repertoire Across Life Phases of Daphnia pulex.

Although the microcrustacean Daphnia is becoming an organism of choice for proteomic studies, protein expression across its life cycle have not been fully characterized. Proteomes of adult females, juveniles, asexually produced embryos, and the ephippia-resting stages containing sexually produced diapausing freezing- and desiccation-resistant embryos are analyzed. Overall, proteins with known molecular functions are more likely to be detected than proteins with no detectable orthology. Similarly, proteins with stronger gene model support in two independent genome assemblies can be detected, than those without such support. This suggests that the proteomics pipeline can be applied to verify hypothesized proteins, even given questionable reference gene models. In particular, upregulation of vitellogenins and downregulation of actins and myosins in embryos of both types, relative to juveniles and adults, and overrepresentation of cell-cycle related proteins in the developing embryos, relative to diapausing embryos and adults, are observed. Upregulation of small heat-shock proteins and peroxidases, as well as overrepresentation of stress-response proteins in the ephippium relative to the asexually produced non-diapausing embryos, is found. The ephippium also shows upregulation of three trehalose-synthesis proteins and downregulation of a trehalose hydrolase, consistent with the role of trehalose in protection against freezing and desiccation.

Luteinizing Hormone is an effective replacement for hCG to induce ovulation in Xenopus.

Injection of human Chorionic Gonadotropin (hCG) directly into the dorsal lymph sac of Xenopus is a commonly used protocol for induction of ovulation, but recent shortages in the stocks of commercially available hCG as well as lack of a well tested alternative have resulted in frustrating experimental delays in laboratories that predominantly use Xenopus in their research. Mammalian Luteinizing Hormones (LH) share structural similarity, functional equivalency, and bind the same receptor as hCG; this suggests that LH may serve as a good alternative to hCG for promoting ovulation in Xenopus. LH has been found to induce maturation of Xenopus oocytes in vitro, but whether it can be used to induce ovulation in vivo has not been examined. Here we compared the ability of four mammalian LH proteins, bovine (bLH), human (hLH), ovine (oLH), porcine (pLH), to induce ovulation in Xenopus when injected into the dorsal lymph sac of sexually mature females. We find that both ovine and human LH, but not bovine or porcine, are good substitutes for hCG for induction of ovulation in WT and J strain Xenopus laevis and Xenopus tropicalis.

Transcriptomic insights into genetic diversity of protein-coding genes in X. laevis.

We characterize the genetic diversity of Xenopus laevis strains using RNA-seq data and allele-specific analysis. This data provides a catalogue of coding variation, which can be used for improving the genomic sequence, as well as for better sequence alignment, probe design, and proteomic analysis. In addition, we paint a broad picture of the genetic landscape of the species by functionally annotating different classes of mutations with a well-established prediction tool (PolyPhen-2). Further, we specifically compare the variation in the progeny of four crosses: inbred genomic (J)-strain, outbred albino (B)-strain, and two hybrid crosses of J and B strains. We identify a subset of mutations specific to the B strain, which allows us to investigate the selection pressures affecting duplicated genes in this allotetraploid. From these crosses we find the ratio of non-synonymous to synonymous mutations is lower in duplicated genes, which suggests that they are under greater purifying selection. Surprisingly, we also find that function-altering ("damaging") mutations constitute a greater fraction of the non-synonymous variants in this group, which suggests a role for subfunctionalization in coding variation affecting duplicated genes.

Genome sequencing reveals insights into physiology and longevity of the naked mole rat.

The naked mole rat (Heterocephalus glaber) is a strictly subterranean, extraordinarily long-lived eusocial mammal. Although it is the size of a mouse, its maximum lifespan exceeds 30 years, making this animal the longest-living rodent. Naked mole rats show negligible senescence, no age-related increase in mortality, and high fecundity until death. In addition to delayed ageing, they are resistant to both spontaneous cancer and experimentally induced tumorigenesis. Naked mole rats pose a challenge to the theories that link ageing, cancer and redox homeostasis. Although characterized by significant oxidative stress, the naked mole rat proteome does not show age-related susceptibility to oxidative damage or increased ubiquitination. Naked mole rats naturally reside in large colonies with a single breeding female, the 'queen', who suppresses the sexual maturity of her subordinates. They also live in full darkness, at low oxygen and high carbon dioxide concentrations, and are unable to sustain thermogenesis nor feel certain types of pain. Here we report the sequencing and analysis of the naked mole rat genome, which reveals unique genome features and molecular adaptations consistent with cancer resistance, poikilothermy, hairlessness and insensitivity to low oxygen, and altered visual function, circadian rythms and taste sensing. This information provides insights into the naked mole rat's exceptional longevity and ability to live in hostile conditions, in the dark and at low oxygen. The extreme traits of the naked mole rat, together with the reported genome and transcriptome information, offer opportunities for understanding ageing and advancing other areas of biological and biomedical research.

Exploiting polypharmacology for drug target deconvolution.

Polypharmacology (action of drugs against multiple targets) represents a tempting avenue for new drug development; unfortunately, methods capable of exploiting the known polypharmacology of drugs for target deconvolution are lacking. Here, we present an ensemble approach using elastic net regularization combined with mRNA expression profiling and previously characterized data on a large set of kinase inhibitors to identify kinases that are important for epithelial and mesenchymal cell migration. By profiling a selected optimal set of 32 kinase inhibitors in a panel against six cell lines, we identified cell type-specific kinases that regulate cell migration. Our discovery of several informative kinases with a previously uncharacterized role in cell migration (such as Mst and Taok family of MAPK kinases in mesenchymal cells) may represent novel targets that warrant further investigation. Target deconvolution using our ensemble approach has the potential to aid in the rational design of more potent but less toxic drug combinations.

A public resource facilitating clinical use of genomes.

Rapid advances in DNA sequencing promise to enable new diagnostics and individualized therapies. Achieving personalized medicine, however, will require extensive research on highly reidentifiable, integrated datasets of genomic and health information. To assist with this, participants in the Personal Genome Project choose to forgo privacy via our institutional review board- approved "open consent" process. The contribution of public data and samples facilitates both scientific discovery and standardization of methods. We present our findings after enrollment of more than 1,800 participants, including whole-genome sequencing of 10 pilot participant genomes (the PGP-10). We introduce the Genome-Environment-Trait Evidence (GET-Evidence) system. This tool automatically processes genomes and prioritizes both published and novel variants for interpretation. In the process of reviewing the presumed healthy PGP-10 genomes, we find numerous literature references implying serious disease. Although it is sometimes impossible to rule out a late-onset effect, stringent evidence requirements can address the high rate of incidental findings. To that end we develop a peer production system for recording and organizing variant evaluations according to standard evidence guidelines, creating a public forum for reaching consensus on interpretation of clinically relevant variants. Genome analysis becomes a two-step process: using a prioritized list to record variant evaluations, then automatically sorting reviewed variants using these annotations. Genome data, health and trait information, participant samples, and variant interpretations are all shared in the public domain-we invite others to review our results using our participant samples and contribute to our interpretations. We offer our public resource and methods to further personalized medical research.

Techniques for Rapidly Sampling Six Crucial Organs in Adult Xenopus.

Xenopus has been a powerful model organism for understanding vertebrate development and disease for over a hundred years. While experimental analysis and dissection techniques of the embryo have been well documented, descriptions of adult Xenopus structures and organs, together with techniques for working with adults, have not been updated to take into consideration the requirements of such modern approaches as quantitative proteomics and single-cell transcriptomics. The cell-type and gene-centric perspectives require contrasting observations in embryonic stages to those in adult tissues. The organs of the larva undergo significant changes in their overall structure, morphology, and anatomical location all along the larval to adult transition, most notably during massive metamorphosis remodeling. Establishing robust standards for organ identification and dissection is crucial to ensure datasets resulting from studies performed at different laboratories can be consistent. The present protocol identifies six of the organs in the adult Xenopus, demonstrating methods for dissection and sampling of the heart ventricle, liver, fat body, pancreas, paired kidney, and skin of the adult Xenopus. Depending on the preservation methods, the dissected organs can be used for quantitative proteomics, single cell/nuclei transcriptomics, in situ hybridization, immunohistochemistry, histology, etc. This protocol aims to standardize tissue sampling and facilitate multi-lab investigations of the adult organ systems.

On standardization of controls in lifespan studies.

The search for interventions to slow down and even reverse aging is a burgeoning field. The literature cites hundreds of supposedly beneficial pharmacological and genetic interventions in model organisms: mice, rats, flies and worms, where research into physiology is routinely accompanied by lifespan data. However, when experimental animals from one article live as long as controls from another article, comparing the results of interventions across studies can yield misleading outcomes. Theoretically, all lifespan data are ripe for re-analysis: we could contrast the molecular targets and pathways across studies and help focus the further search for interventions. Alas, the results of most longevity studies are difficult to compare. This is in part because there are no clear, universally accepted standards for conducting such experiments or even for reporting such data. The situation is worsened by the fact that the authors often do not describe experimental conditions completely. As a result, works on longevity make up a set of precedents, each of which might be interesting in its own right, yet incoherent and incomparable at least for the reason that in a general context, it may indicate, for example, not prolonging the life of an average organism, but compensating for any genetic abnormalities of a particular sample or inappropriate living conditions. Here we point out specific issues and propose solutions for quality control by checking both inter- and intra-study consistency of lifespan data.