����������̳

The recent advent of methods for high-throughput single-cell molecular profiling has catalyzed a growing sense in the scientific community that the time is ripe to complete the 150-year-old effort to identify all cell types in the human body. The Human Cell Atlas Project is an international collaborative effort that aims to define all human cell types in terms of distinctive molecular profiles (such as gene expression profiles) and to connect this information with classical cellular descriptions (such as location and morphology). An open comprehensive reference map of the molecular state of cells in healthy human tissues would propel the systematic study of physiological states, developmental trajectories, regulatory circuitry and interactions of cells, and also provide a framework for understanding cellular dysregulation in human disease. Here we describe the idea, its potential utility, early proofs-of-concept, and some design considerations for the Human Cell Atlas, including a commitment to open data, code, and community.

V(D)J recombination is essential for the generation of diverse antigen receptor (AgR) repertoires. In B cells, immunoglobulin kappa (Igκ) light chain recombination follows immunoglobulin heavy chain (Igh) recombination. We recently developed the DNA-based VDJ-seq assay for the unbiased quantitation of Igh VH and DH repertoires. Integration of VDJ-seq data with genome-wide datasets revealed that two chromatin states at the recombination signal sequence (RSS) of VH genes are highly predictive of recombination in mouse pro-B cells. It is unknown whether local chromatin states contribute to Vκ gene choice during Igκ recombination. Here we adapt VDJ-seq to profile the Igκ VκJκ repertoire and present a comprehensive readout in mouse pre-B cells, revealing highly variable Vκ gene usage. Integration with genome-wide datasets for histone modifications, DNase hypersensitivity, transcription factor binding and germline transcription identified PU.1 binding at the RSS, which was unimportant for Igh, as highly predictive of whether a Vκ gene will recombine or not, suggesting that it plays a binary, all-or-nothing role, priming genes for recombination. Thereafter, the frequency with which these genes recombine was shaped both by the presence and level of enrichment of several other chromatin features, including H3K4 methylation and IKAROS binding. Moreover, in contrast to the Igh locus, the chromatin landscape of the promoter, as well as of the RSS, contributes to Vκ gene recombination. Thus, multiple facets of local chromatin features explain much of the variation in Vκ gene usage. Together, these findings reveal shared and divergent roles for epigenetic features and transcription factors in AgR V(D)J recombination and provide avenues for further investigation of chromatin signatures that may underpin V(D)J-mediated chromosomal translocations.

Blood coagulation functions as part of the innate immune system by preventing bacterial invasion, and it is critical to stopping blood loss (hemostasis). Coagulation involves the external membrane surface of activated platelets and leukocytes. Using lipidomic, genetic, biochemical, and mathematical modeling approaches, we found that enzymatically oxidized phospholipids (eoxPLs) generated by the activity of leukocyte or platelet lipoxygenases (LOXs) were required for normal hemostasis and promoted coagulation factor activities in a Ca- and phosphatidylserine (PS)-dependent manner. In wild-type mice, hydroxyeicosatetraenoic acid-phospholipids (HETE-PLs) enhanced coagulation and restored normal hemostasis in clotting-deficient animals genetically lacking p12-LOX or 12/15-LOX activity. Murine platelets generated 22 eoxPL species, all of which were missing in the absence of p12-LOX. Humans with the thrombotic disorder antiphospholipid syndrome (APS) had statistically significantly increased HETE-PLs in platelets and leukocytes, as well as greater HETE-PL immunoreactivity, than healthy controls. HETE-PLs enhanced membrane binding of the serum protein β2GP1 (β2-glycoprotein 1), an event considered central to the autoimmune reactivity responsible for APS symptoms. Correlation network analysis of 47 platelet eoxPL species in platelets from APS and control subjects identified their enzymatic origin and revealed a complex network of regulation, with the abundance of 31 p12-LOX-derived eoxPL molecules substantially increased in APS. In summary, circulating blood cells generate networks of eoxPL molecules, including HETE-PLs, which change membrane properties to enhance blood coagulation and contribute to the excessive clotting and immunoreactivity of patients with APS.

Migration of activated regulatory T (Treg) cells to inflamed tissue is crucial for their immune-modulatory function. While metabolic reprogramming during Treg cell differentiation has been extensively studied, the bioenergetics of Treg cell trafficking remains undefined. We have investigated the metabolic demands of migrating Treg cells in vitro and in vivo. We show that glycolysis was instrumental for their migration and was initiated by pro-migratory stimuli via a PI3K-mTORC2-mediated pathway culminating in induction of the enzyme glucokinase (GCK). Subsequently, GCK promoted cytoskeletal rearrangements by associating with actin. Treg cells lacking this pathway were functionally suppressive but failed to migrate to skin allografts and inhibit rejection. Similarly, human carriers of a loss-of-function GCK regulatory protein gene-leading to increased GCK activity-had reduced numbers of circulating Treg cells. These cells displayed enhanced migratory activity but similar suppressive function, while conventional T cells were unaffected. Thus, GCK-dependent glycolysis regulates Treg cell migration.

Expression of the transcription factors OCT4, SOX2, KLF4, and cMYC (OSKM) reprograms somatic cells into induced pluripotent stem cells (iPSCs). Reprogramming is a slow and inefficient process, suggesting the presence of safeguarding mechanisms that counteract cell fate conversion. One such mechanism is senescence. To identify modulators of reprogramming-induced senescence, we performed a genome-wide shRNA screen in primary human fibroblasts expressing OSKM. In the screen, we identified novel mediators of OSKM-induced senescence and validated previously implicated genes such as CDKN1A We developed an innovative approach that integrates single-cell RNA sequencing (scRNA-seq) with the shRNA screen to investigate the mechanism of action of the identified candidates. Our data unveiled regulation of senescence as a novel way by which mechanistic target of rapamycin (mTOR) influences reprogramming. On one hand, mTOR inhibition blunts the induction of cyclin-dependent kinase (CDK) inhibitors (CDKIs), including p16(INK4a), p21(CIP1), and p15(INK4b), preventing OSKM-induced senescence. On the other hand, inhibition of mTOR blunts the senescence-associated secretory phenotype (SASP), which itself favors reprogramming. These contrasting actions contribute to explain the complex effect that mTOR has on reprogramming. Overall, our study highlights the advantage of combining functional screens with scRNA-seq to accelerate the discovery of pathways controlling complex phenotypes.

There is limited knowledge about the metabolic reprogramming induced by cancer therapies and how this contributes to therapeutic resistance. Here we show that although inhibition of PI3K-AKT-mTOR signaling markedly decreased glycolysis and restrained tumor growth, these signaling and metabolic restrictions triggered autophagy, which supplied the metabolites required for the maintenance of mitochondrial respiration and redox homeostasis. Specifically, we found that survival of cancer cells was critically dependent on phospholipase A2 (PLA2) to mobilize lysophospholipids and free fatty acids to sustain fatty acid oxidation and oxidative phosphorylation. Consistent with this, we observed significantly increased lipid droplets, with subsequent mobilization to mitochondria. These changes were abrogated in cells deficient for the essential autophagy gene Accordingly, inhibition of PLA2 significantly decreased lipid droplets, decreased oxidative phosphorylation, and increased apoptosis. Together, these results describe how treatment-induced autophagy provides nutrients for cancer cell survival and identifies novel cotreatment strategies to override this survival advantage.

Pancreatic cancer is a high mortality form of cancer, with a median survival only six months. There are multiple associated risk factors associated, most importantly type 2 diabetes, obesity, pancreatitis and smoking. The relative rarity of the disease, however, has made it difficult to dissect causative risk factors, especially with related risk factors. A major unanswered question with important therapeutic implications is the effect of immunological responses on pancreatic cancer formation, with data from other cancers suggesting the potential for local immunological responses to either increase cancer development or increase cancer elimination. Due to the rarity and late diagnosis of pancreatic cancer direct epidemiological evidence is lacking, thus necessitating a reliance on animal models. Here we investigated the relationship between pancreatic autoimmunity and cancer by backcrossing the well characterised Ela1-Tag transgenic model of pancreatic cancer onto the pancreatic autoimmune susceptible NOD mouse strain. Through longitudinal magnetic resonance imaging we found that the NOD genetic background delayed the onset of pancreatic tumours and substantially slowed the growth rate of tumours after development. These results suggest that elevated autoimmune surveillance of the pancreas limits tumour formation and growth, identifying pancreatic cancer as a promising target for immune checkpoint blockade therapies that unleash latent autoimmunity.

BioModels serves as a central repository of mathematical models representing biological processes. It offers a platform to make mathematical models easily shareable across the systems modelling community, thereby supporting model reuse. To facilitate hosting a broader range of model formats derived from diverse modelling approaches and tools, a new infrastructure for BioModels has been developed that is available at http://www.ebi.ac.uk/biomodels. This new system allows submitting and sharing of a wide range of models with improved support for formats other than SBML. It also offers a version-control backed environment in which authors and curators can work collaboratively to curate models. This article summarises the features available in the current system and discusses the potential benefit they offer to the users over the previous system. In summary, the new portal broadens the scope of models accepted in BioModels and supports collaborative model curation which is crucial for model reproducibility and sharing.

Erasure of DNA methylation and repressive chromatin marks in the mammalian germline leads to risk of transcriptional activation of transposable elements (TEs). Here, we used mouse embryonic stem cells (ESCs) to identify an endosiRNA-based mechanism involved in suppression of TE transcription. In ESCs with DNA demethylation induced by acute deletion of Dnmt1, we saw an increase in sense transcription at TEs, resulting in an abundance of sense/antisense transcripts leading to high levels of ARGONAUTE2 (AGO2)-bound small RNAs. Inhibition of Dicer or Ago2 expression revealed that small RNAs are involved in an immediate response to demethylation-induced transposon activation, while the deposition of repressive histone marks follows as a chronic response. In vivo, we also found TE-specific endosiRNAs present during primordial germ cell development. Our results suggest that antisense TE transcription is a "trap" that elicits an endosiRNA response to restrain acute transposon activity during epigenetic reprogramming in the mammalian germline.

Nicotinamide adenine dinucleotide (NAD) and nicotinamide adenine dinucleotide phosphate (NADP) are conserved metabolic cofactors that mediate reduction-oxidation (redox) reactions throughout all domains of life. The diversity of synthetic routes and cellular processes involving the transfer of reducing equivalents to and from these cofactors makes the accurate quantitation and metabolic tracing of NAD(H) and NADP(H) of broad interest. However, current analytical techniques typically rely on standard curves that do not incorporate confounding effects of the sample matrix. We utilized the essential requirement of niacin and tryptophan for NAD synthesis in mammalian cells to devise a stable isotope labeling by essential nutrients in cell culture (SILEC) method for efficient labeling of intracellular NAD(H) and NADP(H) pools. Coupling this approach with detection by liquid chromatography-high resolution mass spectrometry (LC-HRMS), we demonstrate the utility of incorporating a [CN]-nicotinamide moiety into a library of NAD-derived metabolites for use as internal standards in matrixed samples. Using a two-label system incorporating [CN]-nicotinamide and [C]-tryptophan, we quantify the relative contribution of salvage and de novo NAD synthesis, respectively, in S. cerevisiae and HepG2 human hepatocellular carcinoma cells under basal conditions. As a further proof-of-principle, we demonstrate an improvement in the linear range for quantification of NAD and apply this method to analysis of NAD(H) in mouse liver. This method demonstrates the generalizability of SILEC, and provides a simple method for generating a library of stable isotope labeled internal standards for quantifying and tracing the metabolism of cellular and tissue NAD(H) and NADP(H).

We achieve automated and reliable annotation of lipid species and their molecular structures in high-throughput data from chromatography-coupled tandem mass spectrometry using decision rule sets embedded in Lipid Data Analyzer (LDA; http://genome.tugraz.at/lda2). Using various low- and high-resolution mass spectrometry instruments with several collision energies, we proved the method's platform independence. We propose that the software's reliability, flexibility, and ability to identify novel lipid molecular species may now render current state-of-the-art lipid libraries obsolete.

The PI3K signaling pathway regulates cell growth and movement and is heavily mutated in cancer. Class I PI3Ks synthesize the lipid messenger PI(3,4,5)P3. PI(3,4,5)P3 can be dephosphorylated by 3- or 5-phosphatases, the latter producing PI(3,4)P2. The PTEN tumor suppressor is thought to function primarily as a PI(3,4,5)P3 3-phosphatase, limiting activation of this pathway. Here we show that PTEN also functions as a PI(3,4)P2 3-phosphatase, both in vitro and in vivo. PTEN is a major PI(3,4)P2 phosphatase in Mcf10a cytosol, and loss of PTEN and INPP4B, a known PI(3,4)P2 4-phosphatase, leads to synergistic accumulation of PI(3,4)P2, which correlated with increased invadopodia in epidermal growth factor (EGF)-stimulated cells. PTEN deletion increased PI(3,4)P2 levels in a mouse model of prostate cancer, and it inversely correlated with PI(3,4)P2 levels across several EGF-stimulated prostate and breast cancer lines. These results point to a role for PI(3,4)P2 in the phenotype caused by loss-of-function mutations or deletions in PTEN.

In vitro histone modification (HM) assays are used to characterize the activity of chromatin-modifying enzymes. These assays provide information regarding the modification sites on histones, the product specificity, and the impact of other histone or nucleotide modifications on enzyme activity. In particular, histone methyltransferase (HMT) assays have been instrumental in elucidating the activity and site specificity of many plant HMT enzymes. In this chapter, we describe a general protocol that can be used to perform HMT assays using different chromatin substrates, detection methods, and enzymes directly purified from plant material or heterologous sources.

A major limitation of embryo epigenotyping by chromatin immunoprecipitation analysis is the reduced amount of sample available from an embryo biopsy. We developed an in vitro system to expand trophectoderm cells from an embryo biopsy to overcome this limitation. This work analyzes whether expanded trophectoderm (EX) is representative of the trophectoderm (TE) methylation or adaptation to culture has altered its epigenome. We took a small biopsy from the trophectoderm (30-40 cells) of in vitro produced bovine-hatched blastocysts and cultured it on fibronectin-treated plates until we obtained ∼4 × 104 cells. The rest of the embryo was allowed to recover its spherical shape and, subsequently, TE and inner cell mass were separated. We examined whether there were DNA methylation differences between TE and EX of three bovine embryos using whole-genome bisulfite sequencing. As a consequence of adaptation to culture, global methylation, including transposable elements, was higher in EX, with 5.3% of quantified regions showing significant methylation differences between TE and EX. Analysis of individual embryos indicated that TE methylation is more similar to its EX counterpart than to TE from other embryos. Interestingly, these similarly methylated regions are enriched in CpG islands, promoters and transcription units near genes involved in biological processes important for embryo development. Our results indicate that EX is representative of the embryo in terms of DNA methylation, thus providing an informative proxy for embryo epigenotyping.

DNA methylation is an important epigenetic modification in many species that is critical for development, and implicated in ageing and many complex diseases, such as cancer. Many cost-effective genome-wide analyses of DNA modifications rely on restriction enzymes capable of digesting genomic DNA at defined sequence motifs. There are hundreds of restriction enzyme families but few are used to date, because no tool is available for the systematic evaluation of restriction enzyme combinations that can enrich for certain sites of interest in a genome. Herein, we present customised Reduced Representation Bisulfite Sequencing (cuRRBS), a novel and easy-to-use computational method that solves this problem. By computing the optimal enzymatic digestions and size selection steps required, cuRRBS generalises the traditional MspI-based Reduced Representation Bisulfite Sequencing (RRBS) protocol to all restriction enzyme combinations. In addition, cuRRBS estimates the fold-reduction in sequencing costs and provides a robustness value for the personalised RRBS protocol, allowing users to tailor the protocol to their experimental needs. Moreover, we show in silico that cuRRBS-defined restriction enzymes consistently out-perform MspI digestion in many biological systems, considering both CpG and CHG contexts. Finally, we have validated the accuracy of cuRRBS predictions for single and double enzyme digestions using two independent experimental datasets.

Mouse embryonic stem cells derived from the epiblast contribute to the somatic lineages and the germline but are excluded from the extra-embryonic tissues that are derived from the trophectoderm and the primitive endoderm upon reintroduction to the blastocyst. Here we report that cultures of expanded potential stem cells can be established from individual eight-cell blastomeres, and by direct conversion of mouse embryonic stem cells and induced pluripotent stem cells. Remarkably, a single expanded potential stem cell can contribute both to the embryo proper and to the trophectoderm lineages in a chimaera assay. Bona fide trophoblast stem cell lines and extra-embryonic endoderm stem cells can be directly derived from expanded potential stem cells in vitro. Molecular analyses of the epigenome and single-cell transcriptome reveal enrichment for blastomere-specific signature and a dynamic DNA methylome in expanded potential stem cells. The generation of mouse expanded potential stem cells highlights the feasibility of establishing expanded potential stem cells for other mammalian species.

Hi-C is a powerful method to investigate genome-wide, higher-order chromatin and chromosome conformations averaged from a population of cells. To expand the potential of Hi-C for single-cell analysis, we developed single-cell Hi-C. Similar to the existing "ensemble" Hi-C method, single-cell Hi-C detects proximity-dependent ligation events between cross-linked and restriction-digested chromatin fragments in cells. A major difference between the single-cell Hi-C and ensemble Hi-C protocol is that the proximity-dependent ligation is carried out in the nucleus. This allows the isolation of individual cells in which nearly the entire Hi-C procedure has been carried out, enabling the production of a Hi-C library and data from individual cells. With this new method, we studied genome conformations and found evidence for conserved topological domain organization from cell to cell, but highly variable interdomain contacts and chromosome folding genome wide. In addition, we found that the single-cell Hi-C protocol provided cleaner results with less technical noise suggesting it could be used to improve the ensemble Hi-C technique.

Single-cell multi-omics has recently emerged as a powerful technology by which different layers of genomic output-and hence cell identity and function-can be recorded simultaneously. Integrating various components of the epigenome into multi-omics measurements allows for studying cellular heterogeneity at different time scales and for discovering new layers of molecular connectivity between the genome and its functional output. Measurements that are increasingly available range from those that identify transcription factor occupancy and initiation of transcription to long-lasting and heritable epigenetic marks such as DNA methylation. Together with techniques in which cell lineage is recorded, this multilayered information will provide insights into a cell's past history and its future potential. This will allow new levels of understanding of cell fate decisions, identity, and function in normal development, physiology, and disease.

The small GTPase RhoA is involved in a variety of fundamental processes in normal tissue. Spatiotemporal control of RhoA is thought to govern mechanosensing, growth, and motility of cells, while its deregulation is associated with disease development. Here, we describe the generation of a RhoA-fluorescence resonance energy transfer (FRET) biosensor mouse and its utility for monitoring real-time activity of RhoA in a variety of native tissues in vivo. We assess changes in RhoA activity during mechanosensing of osteocytes within the bone and during neutrophil migration. We also demonstrate spatiotemporal order of RhoA activity within crypt cells of the small intestine and during different stages of mammary gestation. Subsequently, we reveal co-option of RhoA activity in both invasive breast and pancreatic cancers, and we assess drug targeting in these disease settings, illustrating the potential for utilizing this mouse to study RhoA activity in vivo in real time.

Upon recognition of auto-antigens, thymocytes are negatively selected or diverted to a regulatory T cell (Treg) fate. CCR7 is required for negative selection of auto-reactive thymocytes in the thymic medulla. Here, we describe an unanticipated contribution of CCR7 to intrathymic Treg generation. Ccr7 mice have increased Treg cellularity because of a hematopoietic but non-T cell autonomous CCR7 function. CCR7 expression by thymic dendritic cells (DCs) promotes survival of mature Sirpα DCs. Thus, CCR7 deficiency results in apoptosis of Sirpα DCs, which is counterbalanced by expansion of immature Sirpα DCs that efficiently induce Treg generation. CCR7 deficiency results in enhanced intrathymic generation of Tregs at the neonatal stage and in lymphopenic adults, when Treg differentiation is critical for establishing self-tolerance. Together, these results reveal a complex function for CCR7 in thymic tolerance induction, where CCR7 not only promotes negative selection but also governs intrathymic Treg generation via non-thymocyte intrinsic mechanisms.