Search

Comprehensive characterization of cellular heterogeneity and the underlying regulatory landscapes of tissues and organs requires a highly robust and scalable method to acquire matched RNA and chromatin accessibility profiles on the same cells. Here, we describe a single-nucleus chromatin accessibility and mRNA expression sequencing 2 (SNARE-seq2) assay, implemented with cellular combinatorial indexing. This method involves tagmentation within permeabilized and fixed single-nucleus isolates to capture accessible chromatin (AC) regions, followed by the capture and reverse transcription of RNA transcripts. Through combinatorial split pool ligations, cDNA and AC within each single nucleus become appended with a common cell barcode combination. The captured cDNA and AC are then co-amplified before splitting and enrichment into single-nucleus RNA and single-nucleus AC sequencing libraries. This protocol is compatible with both nuclei and whole cells and can be completed in 3.5 d. SNARE-seq2 permits robust generation of high-quality, joint single-cell RNA and AC sequencing libraries from hundreds of thousands of single cells per experiment.

To study the heterogeneity of complex tissues by joint profiling of gene expression and its regulation, we require an accurate and high-throughput method. Here we described improved high-throughput combinatorial indexing-based single-nucleus chromatin accessibility and mRNA expression sequencing 2 (SNARE-Seq2) co-assay. This protocol involves fixing and permeabilizing the nucleus followed by tagmentation, chromatin barcode ligation, reverse transcription, pooling and splitting for the next rounds of cell barcode ligation into cDNA and accessible chromatin (AC) on the same nucleus. The captured cDNA and AC are co-amplified before splitting and enrichment into single-nucleus RNA and single-nucleus AC sequencing libraries. The protocol can also be applied to both nuclei and whole cells to capture mRNA in the cytoplasm. This improvement allows us to generate hundreds of thousands of data set of each assay and can be scaled up to half a million cells from a single experiment. The entire procedure can be complete in 3.5 d for generating joint single-nucleus RNA and single-nucleus ATAC sequencing libraries.

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch–seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations., An examination of motor cortex in humans, marmosets and mice reveals a generally conserved cellular makeup that is likely to extend to many mammalian species, but also differences in gene expression, DNA methylation and chromatin state that lead to species-dependent specializations.

Understanding kidney disease relies upon defining the complexity of cell types and states, their associated molecular profiles, and interactions within tissue neighborhoods. We have applied multiple single-cell or -nucleus assays (>400,000 nuclei/cells) and spatial imaging technologies to a broad spectrum of healthy reference (n = 42) and disease (n = 42) kidneys. This has provided a high resolution cellular atlas of 100 cell types that include rare and novel cell populations. The multi-omic approach provides detailed transcriptomic profiles, epigenomic regulatory factors, and spatial localizations for major cell types spanning the entire kidney. We further identify and define cellular states altered in kidney injury, encompassing cycling, adaptive or maladaptive repair, transitioning and degenerative states affecting several segments. Molecular signatures of these states permitted their localization within injury neighborhoods using spatial transcriptomics, and large-scale 3D imaging analysis of ∼1.2 million neighborhoods provided linkages to active immune responses. These analyses further defined biological pathways relevant to injury niches, including signatures underlying the transition from reference to predicted maladaptive states that were associated with a decline in kidney function during chronic kidney disease. This human kidney cell atlas, including injury cell states and neighborhoods, will be a valuable resource for future studies.

To study the heterogeneity of complex tissues by joint profiling of gene expression and its regulation, we require an accurate and high-throughput method. Here we described improved high-throughput combinatorial indexing-based single-nucleus chromatin accessibility and mRNA expression sequencing 2 (SNARE-Seq2) co-assay. This protocol involves fixing and permeabilizing the nucleus followed by tagmentation, chromatin barcode ligation, reverse transcription, pooling and splitting for the next rounds of cell barcode ligation into cDNA and accessible chromatin (AC) on the same nucleus. The captured cDNA and AC are co-amplified before splitting and enrichment into single-nucleus RNA and single-nucleus AC sequencing libraries. The protocol can also be applied to both nuclei and whole cells to capture mRNA in the cytoplasm. This improvement allows us to generate hundreds of thousands of data set of each assay and can be scaled up to half a million cells from a single experiment. The entire procedure can be complete in 3.5 d for generating joint single-nucleus RNA and single-nucleus ATAC sequencing libraries.

The primary motor cortex (M1) is essential for voluntary fine motor control and is functionally conserved across mammals. Using high-throughput transcriptomic and epigenomic profiling of over 450,000 single nuclei in human, marmoset monkey, and mouse, we demonstrate a broadly conserved cellular makeup of this region, whose similarity mirrors evolutionary distance and is consistent between the transcriptome and epigenome. The core conserved molecular identity of neuronal and non-neuronal types allowed the generation of a cross-species consensus cell type classification and inference of conserved cell type properties across species. Despite overall conservation, many species specializations were apparent, including differences in cell type proportions, gene expression, DNA methylation, and chromatin state. Few cell type marker genes were conserved across species, providing a short list of candidate genes and regulatory mechanisms responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allowed the Patch-seq identification of layer 5 (L5) corticospinal Betz cells in non-human primate and human and characterization of their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell type diversity in M1 across mammals and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

Comprehensive characterization of cellular heterogeneity and the underlying regulatory landscapes of tissues and organs requires a highly robust and scalable method to acquire matched RNA and chromatin accessibility profiles on the same cells. Here, we describe a single-nucleus chromatin accessibility and mRNA expression sequencing 2 (SNARE-seq2) assay, implemented with cellular combinatorial indexing. This method involves tagmentation within permeabilized and fixed single-nucleus isolates to capture accessible chromatin (AC) regions, followed by the capture and reverse transcription of RNA transcripts. Through combinatorial split pool ligations, cDNA and AC within each single nucleus become appended with a common cell barcode combination. The captured cDNA and AC are then co-amplified before splitting and enrichment into single-nucleus RNA and single-nucleus AC sequencing libraries. This protocol is compatible with both nuclei and whole cells and can be completed in 3.5 d. SNARE-seq2 permits robust generation of high-quality, joint single-cell RNA and AC sequencing libraries from hundreds of thousands of single cells per experiment.

Herceptin or trastuzumab (anti-HER2) is a widely used therapeutic antibody to treat breast cancer with 2017 worldwide sales of US$ 7.5 billion. As the patent on Herceptin expired in Europe in July 2014 and in the US in June 2019, its biosimilars developments have been of great interest and some of them have already been approved in many countries. The production of monoclonal antibody is still relied heavily on mammalian cell culture. Nonetheless, for heterologous protein production, Escherichia coli (E. coli) is often considered first if possible as large scale bioprocessing in the bacterial system is simpler and cheaper than the mammalian system. Recently, E. coli strains called SHuffle strains were specifically engineered for expression of multi-disulfide bonded and difficult-to-express proteins. The aim of our study was to investigate the effects of temperature and induction level on full-length trastuzumab production in the cytoplasm of two different types of SHuffle strains at various growth phases in order to identify the suitable conditions. Trastuzumab was produced as ‘cyclonal’ from pMAZ360-cIgG-Herceptin plasmid. Western blot was used to analyze the solubility and relative assembly of full-length trastuzumab produced under each investigated conditions. It was found that the most suitable condition from our study for trastuzumab production in the cytoplasm of the engineered E. coli was the expression at 30°C with 0.1 mM IPTG in SHuffle T7 express, which was developed based on a protease deficient B strain. In the future, these conditions may be used as guidelines for scaling up trastuzumab production in a fermenter for large-scale production.

The reduced graphene oxide (rGO) and poly(L-methionine) film on screen printed carbon electrodes (SPCE) was proposed for tartrazine determination. The rGO was synthesized based on glucose reduction and modified on SPCE by drop casting. Cyclic voltammetry (CV) was applied to electropolymerize poly(L-methionine) film on SPCE. The electrochemical behaviour of tartrazine on rGO-methionine modified SPCE and the effects of pH were characterized by CV. The modified electrode exhibited a linear response to tartrazine at concentrations from 10 to 200 μM with a detection limit of 0.287 μM. Due to the facile preparation, disposability, and low cost, this modified sensor can be a good alternative for tartrazine determination.

Defining cellular and molecular identities within the kidney is necessary to understand its organization and function in health and disease. Here we demonstrate a reproducible method with minimal artifacts for single-nucleus Droplet-based RNA sequencing (snDrop-Seq) that we use to resolve thirty distinct cell populations in human adult kidney. We define molecular transition states along more than ten nephron segments spanning two major kidney regions. We further delineate cell type-specific expression of genes associated with chronic kidney disease, diabetes and hypertension, providing insight into possible targeted therapies. This includes expression of a hypertension-associated mechano-sensory ion channel in mesangial cells, and identification of proximal tubule cell populations defined by pathogenic expression signatures. Our fully optimized, quality-controlled transcriptomic profiling pipeline constitutes a tool for the generation of healthy and diseased molecular atlases applicable to clinical samples., Single-cell studies in solid tissues remain challenging and have benefited from the development of single-nuclei RNA sequencing strategies. Here Lake et al. apply single-nucleus RNA sequencing to human kidney tissues to provide a comprehensive molecular and cellular atlas of the human kidney, with potential implications for the understanding of kidney physiology and disease.

Bisulfite padlock probes (BSPP) are a method for the targeted quantification of DNA methylation in mammalian genomes. They can simultaneously characterize the level of methylcytosine modification in a large number of targeted regions at single-base resolution. A major advantage of BSPP is that it allows the flexible capture of an arbitrary subset of genomic regions (hundreds to hundreds of thousands of genomic loci) in single-tube reactions. Large number of samples can be processed efficiently and converted into multiplexed sequencing libraries with only three enzymatic steps, without the conventional library preparation procedures. BSPP are applicable to clinical studies, screening cell lines, and for quantifying low abundance regions using deep sequencing.

Bisulfite padlock probes (BSPP) are a method for the targeted quantification of DNA methylation in mammalian genomes. They can simultaneously characterize the level of methylcytosine modification in a large number of targeted regions at single-base resolution. A major advantage of BSPP is that it allows the flexible capture of an arbitrary subset of genomic regions (hundreds to hundreds of thousands of genomic loci) in single-tube reactions. Large number of samples can be processed efficiently and converted into multiplexed sequencing libraries with only three enzymatic steps, without the conventional library preparation procedures. BSPP are applicable to clinical studies, screening cell lines, and for quantifying low abundance regions using deep sequencing.

Chemical modifications of DNA have been recognized as key epigenetic mechanisms for maintenance of the cellular state and memory. Such DNA modifications include canonical 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxycytosine (5caC). Recent advances in detection and quantification of DNA modifications have enabled epigenetic variation to be connected to phenotypic consequences on an unprecedented scale. These methods may use chemical or enzymatic DNA treatment, may be targeted or non-targeted and may utilize array-based hybridization or sequencing. Key considerations in the choice of assay are cost, minimum sample input requirements, accuracy and throughput. This Review discusses the principles behind recently developed techniques, compares their respective strengths and limitations and provides general guidelines for selecting appropriate methods for specific experimental contexts.

Adjacent CpG sites in mammalian genomes can be co-methylated due to the processivity of methyltransferases or demethylases. Yet discordant methylation patterns have also been observed, and found related to stochastic or uncoordinated molecular processes. We focused on a systematic search and investigation of regions in the full human genome that exhibit highly coordinated methylation. We defined 147,888 blocks of tightly coupled CpG sites, called methylation haplotype blocks (MHBs) with 61 sets of whole genome bisulfite sequencing (WGBS) data, and further validated with 101 sets of reduced representation bisulfite sequencing (RRBS) data and 637 sets of methylation array data. Using a metric called methylation haplotype load (MHL), we performed tissue-specific methylation analysis at the block level. Subsets of informative blocks were further identified for deconvolution of heterogeneous samples. Finally, we demonstrated quantitative estimation of tumor load and tissue-of-origin mapping in the circulating cell-free DNA of 59 cancer patients using methylation haplotypes.

Yellow head virus (YHV) is an invertebrate nidovirus which causes a severe mortality in cultured Penaeus monodon. The mortality may be prevented by prior treatment of shrimps with YHV-protease dsRNA. Whether the YHV infected shrimp might be cured by the dsRNA remains to be investigated. P. monodon injected with 10(-6) YHV showed a high virus replication and mortality within 2 days. Injection of 25 microg YHV-protease dsRNA at 3, 6, 12 or 24 h post YHV infection showed a strong inhibition of YHV replication up to 12 h. Unrelated dsRNA-GFP showed no inhibition, indicating that the inhibition was nucleic acid sequence specific through RNAi pathway. Shrimp mortality could be prevented at 3h post YHV infection by the dsRNA, but not at 24 h. These results demonstrate that YHV-protease dsRNA gives therapeutic effect and pave the way to develop a cure for YHV-infected shrimps.

The ability of cells to survive and recover from deteriorating effects of endoplasmic reticulum (ER) stress relies on the unfolded protein response (UPR). The signaling pathway of Ire1p mediate mRNA splicing plays diverge role in UPR response in different organisms from yeast to mammals. Here, we report that Ire1p mediated XBP1 mRNA splicing mechanism is extremely conserved and exerts a critical role for modulating Xbp1 protein synthesis in Drosophila melanogaster. This system is operative in Drosophila S2 cells as a prominent mechanism to mediate transcriptional activation of UPR responsive genes during ER stress.

Fanconi anaemia (FA) is a recessive disorder characterized by genomic instability, congenital abnormalities, cancer predisposition and bone marrow (BM) failure. However, the pathogenesis of FA is not fully understood partly due to the limitations of current disease models. Here, we derive integration free-induced pluripotent stem cells (iPSCs) from an FA patient without genetic complementation and report in situ gene correction in FA-iPSCs as well as the generation of isogenic FANCA-deficient human embryonic stem cell (ESC) lines. FA cellular phenotypes are recapitulated in iPSCs/ESCs and their adult stem/progenitor cell derivatives. By using isogenic pathogenic mutation-free controls as well as cellular and genomic tools, our model serves to facilitate the discovery of novel disease features. We validate our model as a drug-screening platform by identifying several compounds that improve hematopoietic differentiation of FA-iPSCs. These compounds are also able to rescue the hematopoietic phenotype of FA patient BM cells. Altres ajuts: Strategic Priority Research Program of the Chinese Academy of Sciences (XDA01020312), National Basic Research Program of China (973 Program,2014CB964600;2014CB910500), NSFC (81271266, 31222039, 81330008, 31201111, 81371342, 81300261, 81300677), Key Research Program of the Chinese Academy of Sciences (KJZD-EW-TZ-L05), Beijing Natural Science Foundation (7141005; 5142016), the Thousand Young Talents program of China, National Laboratory of Biomacromolecules (012kf02, 2013kf05;2013kf11;2014kf02), and State Key Laboratory of Drug Research (SIMM1302KF-17). M.L. and K.S. are supported by CIRM fellowship. N.M was partially supported by La Fundació Privada La Marató de TV3, 121430/31/32. Y.T. was partially supported by an Uehara Memorial Foundation research fellowship. E.N. was partially supported by an F.M. Kirby Foundation postdoctoral fellowship. J.S. was supported by Fundació Marató TV3 (464/C/2012). J.A.B. was supported by grants from La Fundació Privada La Marató de TV3, 121430/31/32. J.C.I.B. was supported by grants from the G. Harold and Leila Y. Mathers Charitable Foundation, The California Institute of Regenerative Medicine, Ellison Medical Foundation, and The Leona M. and Harry B. Helmsley Charitable Trust grant #2012-PG-MED002.

With defined culture protocol, human embryonic stem cells (hESCs) are able to generate cardiomyocytes in vitro, therefore providing a great model for human heart development, and holding great potential for cardiac disease therapies. In this study, we successfully generated a highly pure population of human cardiomyocytes (hCMs) (>95% cTnT+) from hESC line, which enabled us to identify and characterize an hCM-specific signature, at both the gene expression and DNA methylation levels. Gene functional association network and gene-disease network analyses of these hCM-enriched genes provide new insights into the mechanisms of hCM transcriptional regulation, and stand as an informative and rich resource for investigating cardiac gene functions and disease mechanisms. Moreover, we show that cardiac-structural genes and cardiac-transcription factors have distinct epigenetic mechanisms to regulate their gene expression, providing a better understanding of how the epigenetic machinery coordinates to regulate gene expression in different cell types. Electronic supplementary material The online version of this article (doi:10.1007/s13238-013-0016-x) contains supplementary material, which is available to authorized users.

Genetic polymorphisms can shape the global landscape of DNA methylation, by either changing substrates for DNA methyltransferases or altering the DNA binding affinity of cis-regulatory proteins. The interactions between CpG methylation and genetic polymorphisms have been previously investigated by methylation quantitative trait loci (mQTL) and allele-specific methylation (ASM) analysis. However, it remains unclear whether these approaches can effectively and comprehensively identify all genetic variants that contribute to the inter-individual variation of DNA methylation levels. Here we used three independent approaches to systematically investigate the influence of genetic polymorphisms on variability in DNA methylation by characterizing the methylation state of 96 whole blood samples in 52 parent-child trios from 22 nuclear pedigrees. We performed targeted bisulfite sequencing with padlock probes to quantify the absolute DNA methylation levels at a set of 411,800 CpG sites in the human genome. With mid-parent offspring analysis (MPO), we identified 10,593 CpG sites that exhibited heritable methylation patterns, among which 70.1% were SNPs directly present in methylated CpG dinucleotides. We determined the mQTL analysis identified 49.9% of heritable CpG sites for which regulation occurred in a distal cis-regulatory manner, and that ASM analysis was only able to identify 5%. Finally, we identified hundreds of clusters in the human genome for which the degree of variation of CpG methylation, as opposed to whether or not CpG sites were methylated, was associated with genetic polymorphisms, supporting a recent hypothesis on the genetic influence of phenotypic plasticity. These results show that cis-regulatory SNPs identified by mQTL do not comprise the full extent of heritable CpG methylation, and that ASM appears overall unreliable. Overall, the extent of genome-methylome interactions is well beyond what is detectible with the commonly used mQTL and ASM approaches, and is likely to include effects on plasticity.

Metabolism is vital to every aspect of cell function, yet the metabolome of induced pluripotent stem cells (iPSCs) remains largely unexplored. Here we report, using an untargeted metabolomics approach, that human iPSCs share a pluripotent metabolomic signature with embryonic stem cells (ESCs) that is distinct from their parental cells, and that is characterized by changes in metabolites involved in cellular respiration. Examination of cellular bioenergetics corroborated with our metabolomic analysis, and demonstrated that somatic cells convert from an oxidative state to a glycolytic state in pluripotency. Interestingly, the bioenergetics of various somatic cells correlated with their reprogramming efficiencies. We further identified metabolites that differ between iPSCs and ESCs, which revealed novel metabolic pathways that play a critical role in regulating somatic cell reprogramming. Our findings are the first to globally analyze the metabolome of iPSCs, and provide mechanistic insight into a new layer of regulation involved in inducing pluripotency, and in evaluating iPSC and ESC equivalence.

Targeted quantification of DNA methylation allows for interrogation of the most informative loci across many samples quickly and cost-effectively. Here we report improved bisulfite padlock probes (BSPPs) with a design algorithm to generate efficient padlock probes, a library-free protocol that dramatically reduces sample-preparation cost and time and is compatible with automation, and an efficient bioinformatics pipeline to accurately obtain both methylation levels and genotypes from sequencing of bisulfite-converted DNA.

Generation of human induced pluripotent stem cells (hiPSCs) by the expression of specific transcription factors depends on successful epigenetic reprogramming to a pluripotent state. Although hiPSCs and human embryonic stem cells (hESCs) display a similar epigenome, recent reports demonstrated the persistence of specific epigenetic marks from the somatic cell type of origin and aberrant methylation patterns in hiPSCs. However, it remains unknown whether the use of different somatic cell sources, encompassing variable levels of selection pressure during reprogramming, influences the level of epigenetic aberrations in hiPSCs. In this work, we characterized the epigenomic integrity of 17 hiPSC lines derived from six different cell types with varied reprogramming efficiencies. We demonstrate that epigenetic aberrations are a general feature of the hiPSC state and are independent of the somatic cell source. Interestingly, we observe that the reprogramming efficiency of somatic cell lines inversely correlates with the amount of methylation change needed to acquire pluripotency. Additionally, we determine that both shared and line-specific epigenetic aberrations in hiPSCs can directly translate into changes in gene expression in both the pluripotent and differentiated states. Significantly, our analysis of different hiPSC lines from multiple cell types of origin allow us to identify a reprogramming-specific epigenetic signature comprised of nine aberrantly methylated genes that is able to segregate hESC and hiPSC lines regardless of the somatic cell source or differentiation state.

International audience; Combination of stem cell-based approaches with gene-editing technologies represents an attractive strategy for studying human disease and developing therapies. However, gene-editing methodologies described to date for human cells suffer from technical limitations including limited target gene size, low targeting efficiency at transcriptionally inactive loci, and off-target genetic effects that could hamper broad clinical application. To address these limitations, and as a proof of principle, we focused on homologous recombination-based gene correction of multiple mutations on lamin A (LMNA), which are associated with various degenerative diseases. We show that helper-dependent adenoviral vectors (HDAdVs) provide a highly efficient and safe method for correcting mutations in large genomic regions in human induced pluripotent stem cells and can also be effective in adult human mesenchymal stem cells. This type of approach could be used to generate genotype-matched cell lines for disease modeling and drug discovery and potentially also in therapeutics.