Integrated pseudotime analysis of human pre-implantation embryo single-cell transcriptomes reveals the dynamics of lineage specification

Understanding lineage specification during human pre-implantation development is a gateway to improving assisted reproductive technologies and stem cell research. Here we employ pseudotime analysis of single-cell RNA sequencing (scRNA-seq) data to reconstruct early mouse and human embryo development. Using time-lapse imaging of annotated embryos, we provide an integrated, ordered, and continuous analysis of transcriptomics changes throughout human development. We reveal that human trophectoderm/inner cell mass transcriptomes diverge at the transition from the B2 to the B3 blastocyst stage, just before blastocyst expansion. We explore the dynamics of the fate markers IFI16 and GATA4 and show that they gradually become mutually exclusive upon establishment of epiblast and primitive endoderm fates, respectively. We also provide evidence that NR2F2 marks trophectoderm maturation, initiating from the polar side, and subsequently spreads to all cells after implantation. Our study pinpoints the precise timing of lineage specification events in the human embryo and identifies transcriptomics hallmarks and cell fate markers.

Light-activated cell identification and sorting (LACIS) for selection of edited clones on a nanofluidic device

Despite improvements in the CRISPR molecular toolbox, identifying and purifying properly edited clones remains slow, laborious, and low-yield. Here, we establish a method to enable clonal isolation, selection, and expansion of properly edited cells, using OptoElectroPositioning technology for single-cell manipulation on a nanofluidic device. Briefly, after electroporation of primary T cells with CXCR4-targeting Cas9 ribonucleoproteins, single T cells are isolated on a chip and expanded into colonies. Phenotypic consequences of editing are rapidly assessed on-chip with cell-surface staining for CXCR4. Furthermore, individual colonies are identified based on their specific genotype. Each colony is split and sequentially exported for on-target sequencing and further off-chip clonal expansion of the validated clones. Using this method, single-clone editing efficiencies, including the rate of mono- and bi-allelic indels or precise nucleotide replacements, can be assessed within 10 days from Cas9 ribonucleoprotein introduction in cells.

Annamaria Mocciaro et al. present LACIS, a method for identifying and selecting gene-edited clones using a microfluidic device. The authors apply LACIS to primary T cells after CRISPR-Cas9 editing of CXCR4 and show that selection of edited clones was possible within 10 days from initiation of gene editing.

Parallel derivation of isogenic human primed and naive induced pluripotent stem cells

Induced pluripotent stem cells (iPSCs) have considerably impacted human developmental biology and regenerative medicine, notably because they circumvent the use of cells of embryonic origin and offer the potential to generate patient-specific pluripotent stem cells. However, conventional reprogramming protocols produce developmentally advanced, or primed, human iPSCs (hiPSCs), restricting their use to post-implantation human development modeling. Hence, there is a need for hiPSCs resembling preimplantation naive epiblast. Here, we develop a method to generate naive hiPSCs directly from somatic cells, using OKMS overexpression and specific culture conditions, further enabling parallel generation of their isogenic primed counterparts. We benchmark naive hiPSCs against human preimplantation epiblast and reveal remarkable concordance in their transcriptome, dependency on mitochondrial respiration and X-chromosome status. Collectively, our results are essential for the understanding of pluripotency regulation throughout preimplantation development and generate new opportunities for disease modeling and regenerative medicine.

A Comparison of mRNA Sequencing with Random Primed and 3'-Directed Libraries

Creating a cDNA library for deep mRNA sequencing (mRNAseq) is generally done by random priming, creating multiple sequencing fragments along each transcript. A 3'-end-focused library approach cannot detect differential splicing, but has potentially higher throughput at a lower cost, along with the ability to improve quantification by using transcript molecule counting with unique molecular identifiers (UMI) that correct PCR bias. Here, we compare an implementation of such a 3'-digital gene expression (3'-DGE) approach with "conventional" random primed mRNAseq. Given our particular datasets on cultured human cardiomyocyte cell lines, we find that, while conventional mRNAseq detects ~15% more genes and needs ~500,000 fewer reads per sample for equivalent statistical power, the resulting differentially expressed genes, biological conclusions, and gene signatures are highly concordant between two techniques. We also find good quantitative agreement at the level of individual genes between two techniques for both read counts and fold changes between given conditions. We conclude that, for high-throughput applications, the potential cost savings associated with 3'-DGE approach are likely a reasonable tradeoff for modest reduction in sensitivity and inability to observe alternative splicing, and should enable many larger scale studies focusing on not only differential expression analysis, but also quantitative transcriptome profiling.

Organism-Level Analysis of Vaccination Reveals Networks of Protection across Tissues

A fundamental challenge in immunology is to decipher the principles governing immune responses at the whole-organism scale. Here, using a comparative infection model, we observe immune signal propagation within and between organs to obtain a dynamic map of immune processes at the organism level. We uncover two inter-organ mechanisms of protective immunity mediated by soluble and cellular factors. First, analyzing ligand-receptor connectivity across tissues reveals that type I IFNs trigger a whole-body antiviral state, protecting the host within hours after skin vaccination. Second, combining parabiosis, single-cell analyses, and gene knockouts, we uncover a multi-organ web of tissue-resident memory T cells that functionally adapt to their environment to stop viral spread across the organism. These results have implications for manipulating tissue-resident memory T cells through vaccination and open up new lines of inquiry for the analysis of immune responses at the organism level.

The leak channel NALCN controls tonic firing and glycolytic sensitivity of substantia nigra pars reticulata neurons

Certain neuron types fire spontaneously at high rates, an ability that is crucial for their function in brain circuits. The spontaneously active GABAergic neurons of the substantia nigra pars reticulata (SNr), a major output of the basal ganglia, provide tonic inhibition of downstream brain areas. A depolarizing 'leak' current supports this firing pattern, but its molecular basis remains poorly understood. To understand how SNr neurons maintain tonic activity, we used single-cell RNA sequencing to determine the transcriptome of individual mouse SNr neurons. We discovered that SNr neurons express the sodium leak channel, NALCN, and that SNr neurons lacking NALCN have impaired spontaneous firing. In addition, NALCN is involved in the modulation of excitability by changes in glycolysis and by activation of muscarinic acetylcholine receptors. Our findings suggest that disruption of NALCN could impair the basal ganglia circuit, which may underlie the severe motor deficits in humans carrying mutations in NALCN.

DOI:http://dx.doi.org/10.7554/eLife.15271.001

Some neurons in the brain produce electrical signals (or “fire”) spontaneously, without receiving any other signals from the senses or from other neurons. This spontaneous activity has a number of important roles. For example, in a part of the brain known as the substantia nigra pars reticulata (SNr), spontaneously active neurons frequently produce electrical signals that reduce electrical activity in other brain areas.

A current of positively charged ions constantly flows into the spontaneously active SNr neurons and enables them to fire constantly. Ions enter neurons through proteins called ion channels that are embedded in the surface of the neuron. Like all proteins, ion channels are made by “transcribing” genes to form molecules of RNA that are then “translated” to produce the basic sequence of the protein.

Lutas et al. have now used single-cell RNA sequencing to study SNr neurons from mice and investigate which ion channel the positive ion current flows through. The RNA sequences revealed that the neurons have the gene for an ion channel known as NALCN. Recordings of the firing rate of neurons in slices of mouse brain showed that SNr neurons without this channel did not fire as often as SNr neurons with the channel. In addition, neurotransmitters (chemicals that alter the ability of neurons to fire) and changes in cell metabolism had less of an effect on the firing rate of SNr neurons that lacked the NALCN channel than they do on normal neurons.

These findings may help explain why people with mutations in the NALCN gene have movement disorders, as the substantia nigra pars reticulata plays an important role in orchestrating complex movements. Future work is now needed to understand how a change in NALCN activity affects the other brain areas that SNr neurons connect to.

DOI:http://dx.doi.org/10.7554/eLife.15271.002

A targeted functional RNA interference screen uncovers glypican 5 as an entry factor for hepatitis B and D viruses

Chronic hepatitis B and D infections are major causes of liver disease and hepatocellular carcinoma worldwide. Efficient therapeutic approaches for cure are absent. Sharing the same envelope proteins, hepatitis B virus and hepatitis delta virus use the sodium/taurocholate cotransporting polypeptide (a bile acid transporter) as a receptor to enter hepatocytes. However, the detailed mechanisms of the viral entry process are still poorly understood. Here, we established a high-throughput infectious cell culture model enabling functional genomics of hepatitis delta virus entry and infection. Using a targeted RNA interference entry screen, we identified glypican 5 as a common host cell entry factor for hepatitis B and delta viruses.

CONCLUSION

These findings advance our understanding of virus cell entry and open new avenues for curative therapies. As glypicans have been shown to play a role in the control of cell division and growth regulation, virus-glypican 5 interactions may also play a role in the pathogenesis of virus-induced liver disease and cancer.

Integrative Analyses of Human Reprogramming Reveal Dynamic Nature of Induced Pluripotency

Induced pluripotency is a promising avenue for disease modeling and therapy, but the molecular principles underlying this process, particularly in human cells, remain poorly understood due to donor-to-donor variability and intercellular heterogeneity. Here, we constructed and characterized a clonal, inducible human reprogramming system that provides a reliable source of cells at any stage of the process. This system enabled integrative transcriptional and epigenomic analysis across the human reprogramming timeline at high resolution. We observed distinct waves of gene network activation, including the ordered re-activation of broad developmental regulators followed by early embryonic patterning genes and culminating in the emergence of a signature reminiscent of pre-implantation stages. Moreover, complementary functional analyses allowed us to identify and validate novel regulators of the reprogramming process. Altogether, this study sheds light on the molecular underpinnings of induced pluripotency in human cells and provides a robust cell platform for further studies. PAPERCLIP.

The evolution of lncRNA repertoires and expression patterns in tetrapods

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

Birth and expression evolution of mammalian microRNA genes

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression, yet their origins and functional evolution in mammals remain little understood due to the lack of appropriate comparative data. Using RNA sequencing, we have generated extensive and comparable miRNA data for five organs in six species that represent all main mammalian lineages and birds (the evolutionary outgroup) with the aim to unravel the evolution of mammalian miRNAs. Our analyses reveal an overall expansion of miRNA repertoires in mammals, with threefold accelerated birth rates of miRNA families in placentals and marsupials, facilitated by the de novo emergence of miRNAs in host gene introns. Generally, our analyses suggest a high rate of miRNA family turnover in mammals with many newly emerged miRNA families being lost soon after their formation. Selectively preserved mammalian miRNA families gradually evolved higher expression levels, as well as altered mature sequences and target gene repertoires, and were apparently mainly recruited to exert regulatory functions in nervous tissues. However, miRNAs that originated on the X chromosome evolved high expression levels and potentially diverse functions during spermatogenesis, including meiosis, through selectively driven duplication-divergence processes. Overall, our study thus provides detailed insights into the birth and evolution of mammalian miRNA genes and the associated selective forces.

The evolution of gene expression levels in mammalian organs

Changes in gene expression are thought to underlie many of the phenotypic differences between species. However, large-scale analyses of gene expression evolution were until recently prevented by technological limitations. Here we report the sequencing of polyadenylated RNA from six organs across ten species that represent all major mammalian lineages (placentals, marsupials and monotremes) and birds (the evolutionary outgroup), with the goal of understanding the dynamics of mammalian transcriptome evolution. We show that the rate of gene expression evolution varies among organs, lineages and chromosomes, owing to differences in selective pressures: transcriptome change was slow in nervous tissues and rapid in testes, slower in rodents than in apes and monotremes, and rapid for the X chromosome right after its formation. Although gene expression evolution in mammals was strongly shaped by purifying selection, we identify numerous potentially selectively driven expression switches, which occurred at different rates across lineages and tissues and which probably contributed to the specific organ biology of various mammals.

Specificity of the coxsackievirus B4 VP4 capsid protein investigated in silico

The Enterovirus genus encompasses several species and various serotypes, like coxsackievirus-B1 (CV-B1) to CV-B6, and many variants. The role of these viruses, especially CV-B4, in the pathogenesis of type 1 diabetes is strongly suspected. It has been reported that antibodies directed towards the region of amino acids 11-30 of the VP4 capsid protein enhance the infection of human peripheral blood mononuclear cells with CV-B4. In order to predict the inter- and intra-serotype specificity of the region 11-30 of CV-B4 VP4, 362 available protein sequences of CV-B1 to -B6, CV-A9, and swine vesicular disease virus (SVDV) have been aligned and levels of homology have been calculated. Serine residue substitutions in this region of VP4 were observed without predictable subsequent modification of conformation or charge. The amino acids 16-24 region was the most variable. The sequence of amino acids 16-24 of the CV-B4E2 VP4 protein was highly homologous to those of other CV-B4 (64.4%) whereas there was no homology with CV-B3 and B5 and very low levels of homology with CV-B1 and B2 (3.3% and 9.9%, respectively). In conclusion, the bioinformatic analysis suggests that the region 16-24 of the VP4 capsid protein is the feature of the specificity of the target of infection-enhancing antibodies directed towards CV-B.