Limited access to antigen drives generation of early B cell memory while restraining the plasmablast response.

Cell fate decisions during early B cell activation determine the outcome of responses to pathogens and vaccines. We examined the early B cell response to T-dependent antigen in mice by single-cell RNA sequencing. Early after immunization, a homogeneous population of activated precursors (APs) gave rise to a transient wave of plasmablasts (PBs), followed a day later by the emergence of germinal center B cells (GCBCs). Most APs rapidly exited the cell cycle, giving rise to non-GC-derived early memory B cells (eMBCs) that retained an AP-like transcriptional profile. Rapid decline of antigen availability controlled these events; provision of excess antigen precluded cell cycle exit and induced a new wave of PBs. Fate mapping revealed a prominent contribution of eMBCs to the MBC pool. Quiescent cells with an MBC phenotype dominated the early response to immunization in primates. A reservoir of APs/eMBCs may enable rapid readjustment of the immune response when failure to contain a threat is manifested by increased antigen availability.

A radical switch in clonality reveals a stem cell niche in the epiphyseal growth plate.

Longitudinal bone growth in children is sustained by growth plates, narrow discs of cartilage that provide a continuous supply of chondrocytes for endochondral ossification1. However, it remains unknown how this supply is maintained throughout childhood growth. Chondroprogenitors in the resting zone are thought to be gradually consumed as they supply cells for longitudinal growth1,2, but this model has never been proved. Here, using clonal genetic tracing with multicolour reporters and functional perturbations, we demonstrate that longitudinal growth during the fetal and neonatal periods involves depletion of chondroprogenitors, whereas later in life, coinciding with the formation of the secondary ossification centre, chondroprogenitors acquire the capacity for self-renewal, resulting in the formation of large, stable monoclonal columns of chondrocytes. Simultaneously, chondroprogenitors begin to express stem cell markers and undergo symmetric cell division. Regulation of the pool of self-renewing progenitors involves the hedgehog and mammalian target of rapamycin complex 1 (mTORC1) signalling pathways. Our findings indicate that a stem cell niche develops postnatally in the epiphyseal growth plate, which provides a continuous supply of chondrocytes over a prolonged period.

Schwann cell precursors contribute to skeletal formation during embryonic development in mice and zebrafish.

Immature multipotent embryonic peripheral glial cells, the Schwann cell precursors (SCPs), differentiate into melanocytes, parasympathetic neurons, chromaffin cells, and dental mesenchymal populations. Here, genetic lineage tracing revealed that, during murine embryonic development, some SCPs detach from nerve fibers to become mesenchymal cells, which differentiate further into chondrocytes and mature osteocytes. This occurred only during embryonic development, producing numerous craniofacial and trunk skeletal elements, without contributing to development of the appendicular skeleton. Formation of chondrocytes from SCPs also occurred in zebrafish, indicating evolutionary conservation. Our findings reveal multipotency of SCPs, providing a developmental link between the nervous system and skeleton.

Recurrent Translocations in Topoisomerase Inhibitor-Related Leukemia Are Determined by the Features of DNA Breaks Rather Than by the Proximity of the Translocating Genes.

Topoisomerase inhibitors are widely used in cancer chemotherapy. However, one of the potential long-term adverse effects of such therapy is acute leukemia. A key feature of such therapy-induced acute myeloid leukemia (t-AML) is recurrent chromosomal translocations involving AML1 (RUNX1) or MLL (KMT2A) genes. The formation of chromosomal translocation depends on the spatial proximity of translocation partners and the mobility of the DNA ends. It is unclear which of these two factors might be decisive for recurrent t-AML translocations. Here, we used fluorescence in situ hybridization (FISH) and chromosome conformation capture followed by sequencing (4C-seq) to investigate double-strand DNA break formation and the mobility of broken ends upon etoposide treatment, as well as contacts between translocation partner genes. We detected the separation of the parts of the broken AML1 gene, as well as the increased mobility of these separated parts. 4C-seq analysis showed no evident contacts of AML1 and MLL with loci, implicated in recurrent t-AML translocations, either before or after etoposide treatment. We suggest that separation of the break ends and their increased non-targeted mobility-but not spatial predisposition of the rearrangement partners-plays a major role in the formation of these translocations.

CpG traffic lights are markers of regulatory regions in human genome.

BACKGROUND: DNA methylation is involved in the regulation of gene expression. Although bisulfite-sequencing based methods profile DNA methylation at a single CpG resolution, methylation levels are usually averaged over genomic regions in the downstream bioinformatic analysis. RESULTS: We demonstrate that on the genome level a single CpG methylation can serve as a more accurate predictor of gene expression than an average promoter / gene body methylation. We define CpG traffic lights (CpG TL) as CpG dinucleotides with a significant correlation between methylation and expression of a gene nearby. CpG TL are enriched in all regulatory regions. Among all promoters, CpG TL are especially enriched in poised ones, suggesting involvement of DNA methylation in their regulation. Yet, binding of only a handful of transcription factors, such as NRF1, ETS, STAT and IRF-family members, could be regulated by direct methylation of transcription factor binding sites (TFBS) or its close proximity. For the majority of TF, an alternative scenario is more likely: methylation and inactivation of the whole regulatory element indirectly represses functional TF binding with a CpG TL being a reliable marker of such inactivation. CONCLUSIONS: CpG TL provide a promising insight into mechanisms of enhancer activity and gene regulation linking methylation of single CpG to gene expression. CpG TL methylation can be used as reliable markers of enhancer activity and gene expression in applications, e.g. in clinic where measuring DNA methylation is easier compared to directly measuring gene expression due to more stable nature of DNA.

Genome-Wide DNA Methylation Profiling Reveals Epigenetic Adaptation of Stickleback to Marine and Freshwater Conditions.

The three-spined stickleback (Gasterosteus aculeatus) represents a convenient model to study microevolution-adaptation to a freshwater environment. Although genetic adaptations to freshwater environments are well-studied, epigenetic adaptations have attracted little attention. In this work, we investigated the role of DNA methylation in the adaptation of the marine stickleback population to freshwater conditions. DNA methylation profiling was performed in marine and freshwater populations of sticklebacks, as well as in marine sticklebacks placed into a freshwater environment and freshwater sticklebacks placed into seawater. We showed that the DNA methylation profile after placing a marine stickleback into fresh water partially converged to that of a freshwater stickleback. For six genes including ATP4A ion pump and NELL1, believed to be involved in skeletal ossification, we demonstrated similar changes in DNA methylation in both evolutionary and short-term adaptation. This suggested that an immediate epigenetic response to freshwater conditions can be maintained in freshwater population. Interestingly, we observed enhanced epigenetic plasticity in freshwater sticklebacks that may serve as a compensatory regulatory mechanism for the lack of genetic variation in the freshwater population. For the first time, we demonstrated that genes encoding ion channels KCND3, CACNA1FB, and ATP4A were differentially methylated between the marine and the freshwater populations. Other genes encoding ion channels were previously reported to be under selection in freshwater populations. Nevertheless, the genes that harbor genetic and epigenetic changes were not the same, suggesting that epigenetic adaptation is a complementary mechanism to selection of genetic variants favorable for freshwater environment.

Deep Learning Applications for Predicting Pharmacological Properties of Drugs and Drug Repurposing Using Transcriptomic Data.

Deep learning is rapidly advancing many areas of science and technology with multiple success stories in image, text, voice and video recognition, robotics, and autonomous driving. In this paper we demonstrate how deep neural networks (DNN) trained on large transcriptional response data sets can classify various drugs to therapeutic categories solely based on their transcriptional profiles. We used the perturbation samples of 678 drugs across A549, MCF-7, and PC-3 cell lines from the LINCS Project and linked those to 12 therapeutic use categories derived from MeSH. To train the DNN, we utilized both gene level transcriptomic data and transcriptomic data processed using a pathway activation scoring algorithm, for a pooled data set of samples perturbed with different concentrations of the drug for 6 and 24 hours. In both pathway and gene level classification, DNN achieved high classification accuracy and convincingly outperformed the support vector machine (SVM) model on every multiclass classification problem, however, models based on pathway level data performed significantly better. For the first time we demonstrate a deep learning neural net trained on transcriptomic data to recognize pharmacological properties of multiple drugs across different biological systems and conditions. We also propose using deep neural net confusion matrices for drug repositioning. This work is a proof of principle for applying deep learning to drug discovery and development.

Non-invasive prenatal diagnostics of aneuploidy using next-generation DNA sequencing technologies, and clinical considerations.

Rapidly developing next-generation sequencing (NGS) technologies produce a large amount of data across the whole human genome and allow a large number of DNA samples to be analyzed simultaneously. Screening cell-free fetal DNA (cffDNA) obtained from maternal blood using NGS technologies has provided new opportunities for non-invasive prenatal diagnosis (NIPD) of fetal aneuploidies. One of the major challenges to the analysis of fetal abnormalities is the development of accurate and reliable algorithms capable of analyzing large numbers of short sequence reads. Several such algorithms have recently been developed. Here, we provide a review of recent NGS-based NIPD methods as well as the available algorithms for short-read sequence analysis. We furthermore introduce the practical application of these algorithms for the detection of different types of fetal aneuploidies, and compare the performance, cost and complexity of each approach for clinical deployment. Our review identifies several main technologies and trends in NGS-based NIPD. The main considerations for clinical development for NIPD and screening tests using DNA sequencing are: accuracy, intellectual property, cost and the ability to screen for a wide range of chromosomal abnormalities and genetic defects. The cost of the diagnostic test depends on the sequencing method, diagnostic algorithm and volume of the tests. If the cost of sequencing equipment and reagents remains at or around current levels, targeted approaches for sequencing-based aneuploidy testing and SNP-based methods are preferred.

STRA6 is essential for induction of vascular smooth muscle lineages in human embryonic cardiac outflow tract development.

AIMS: Retinoic acid (RA) signalling is essential for heart development, and dysregulation of the RA signalling can cause several types of cardiac outflow tract (OFT) defects, the most frequent congenital heart disease (CHD) in humans. Matthew-Wood syndrome is caused by inactivating mutations of a transmembrane protein gene STRA6 that transports vitamin A (retinol) from extracellular into intracellular spaces. This syndrome shows a broad spectrum of malformations including CHD, although murine Stra6-null neonates did not exhibit overt heart defects. Thus, the detailed mechanisms by which STRA6 mutations could lead to cardiac malformations in humans remain unclear. Here, we investigated the role of STRA6 in the context of human cardiogenesis and CHD. METHODS AND RESULTS: To gain molecular signatures in species-specific cardiac development, we first compared single-cell RNA sequencing (RNA-seq) datasets, uniquely obtained from human and murine embryonic hearts. We found that while STRA6 mRNA was much less frequently expressed in murine embryonic heart cells derived from the Mesp1+ lineage tracing mice (Mesp1Cre/+; Rosa26tdTomato), it was expressed predominantly in the OFT region-specific heart progenitors in human developing hearts. Next, we revealed that STRA6-knockout human embryonic stem cells (hESCs) could differentiate into cardiomyocytes similarly to wild-type hESCs, but could not differentiate properly into mesodermal nor neural crest cell-derived smooth muscle cells (SMCs) in vitro. This is supported by the population RNA-seq data showing down-regulation of the SMC-related genes in the STRA6-knockout hESC-derived cells. Further, through machinery assays, we identified the previously unrecognized interaction between RA nuclear receptors RARα/RXRα and TBX1, an OFT-specific cardiogenic transcription factor, which would likely act downstream to STRA6-mediated RA signalling in human cardiogenesis. CONCLUSION: Our study highlights the critical role of human-specific STRA6 progenitors for proper induction of vascular SMCs that is essential for normal OFT formation. Thus, these results shed light on novel and human-specific CHD mechanisms, driven by STRA6 mutations.

A branching model of lineage differentiation underpinning the neurogenic potential of enteric glia.

Glial cells have been proposed as a source of neural progenitors, but the mechanisms underpinning the neurogenic potential of adult glia are not known. Using single cell transcriptomic profiling, we show that enteric glial cells represent a cell state attained by autonomic neural crest cells as they transition along a linear differentiation trajectory that allows them to retain neurogenic potential while acquiring mature glial functions. Key neurogenic loci in early enteric nervous system progenitors remain in open chromatin configuration in mature enteric glia, thus facilitating neuronal differentiation under appropriate conditions. Molecular profiling and gene targeting of enteric glial cells in a cell culture model of enteric neurogenesis and a gut injury model demonstrate that neuronal differentiation of glia is driven by transcriptional programs employed in vivo by early progenitors. Our work provides mechanistic insight into the regulatory landscape underpinning the development of intestinal neural circuits and generates a platform for advancing glial cells as therapeutic agents for the treatment of neural deficits.

A previously uncharacterized Factor Associated with Metabolism and Energy (FAME/C14orf105/CCDC198/1700011H14Rik) is related to evolutionary adaptation, energy balance, and kidney physiology.

In this study we use comparative genomics to uncover a gene with uncharacterized function (1700011H14Rik/C14orf105/CCDC198), which we hereby name FAME (Factor Associated with Metabolism and Energy). We observe that FAME shows an unusually high evolutionary divergence in birds and mammals. Through the comparison of single nucleotide polymorphisms, we identify gene flow of FAME from Neandertals into modern humans. We conduct knockout experiments on animals and observe altered body weight and decreased energy expenditure in Fame knockout animals, corresponding to genome-wide association studies linking FAME with higher body mass index in humans. Gene expression and subcellular localization analyses reveal that FAME is a membrane-bound protein enriched in the kidneys. Although the gene knockout results in structurally normal kidneys, we detect higher albumin in urine and lowered ferritin in the blood. Through experimental validation, we confirm interactions between FAME and ferritin and show co-localization in vesicular and plasma membranes.

Bhlhe40 function in activated B and TFH cells restrains the GC reaction and prevents lymphomagenesis.

The generation of high-affinity antibodies against pathogens and vaccines requires the germinal center (GC) reaction, which relies on a complex interplay between specialized effector B and CD4 T lymphocytes, the GC B cells and T follicular helper (TFH) cells. Intriguingly, several positive key regulators of the GC reaction are common for both cell types. Here, we report that the transcription factor Bhlhe40 is a crucial cell-intrinsic negative regulator affecting both the B and T cell sides of the GC reaction. In activated CD4 T cells, Bhlhe40 was required to restrain proliferation, thus limiting the number of TFH cells. In B cells, Bhlhe40 executed its function in the first days after immunization by selectively restricting the generation of the earliest GC B cells but not of early memory B cells or plasmablasts. Bhlhe40-deficient mice with progressing age succumbed to a B cell lymphoma characterized by the accumulation of monoclonal GC B-like cells and polyclonal TFH cells in various tissues.

An IDH-independent mechanism of DNA hypermethylation upon VHL inactivation in cancer.

Hypermethylation of tumour suppressors and other aberrations of DNA methylation in tumours play a significant role in cancer progression. DNA methylation can be affected by various environmental conditions, including hypoxia. The response to hypoxia is mainly achieved through activation of the transcriptional program associated with HIF1A transcription factor. Inactivation of Von Hippel-Lindau Tumour Suppressor gene (VHL) by genetic or epigenetic events, which also induces aberrant activation of HIF1A, is the most common driver event for renal cancer. With whole-genome bisulphite sequencing and LC-MS, we demonstrated that VHL inactivation induced global genome hypermethylation in human kidney cancer cells under normoxic conditions. This effect was reverted by exogenous expression of wild-type VHL. We showed that global genome hypermethylation in VHL mutants can be explained by transcriptional changes in MDH and L2HGDH genes that cause the accumulation of 2-hydroxyglutarate - a metabolite that inhibits DNA demethylation by TET enzymes. Unlike the known cases of DNA hypermethylation in cancer, 2-hydroxyglutarate was accumulated in the cells with the wild-type isocitrate dehydrogenases.

Serotonin limits generation of chromaffin cells during adrenal organ development.

Adrenal glands are the major organs releasing catecholamines and regulating our stress response. The mechanisms balancing generation of adrenergic chromaffin cells and protecting against neuroblastoma tumors are still enigmatic. Here we revealed that serotonin (5HT) controls the numbers of chromaffin cells by acting upon their immediate progenitor "bridge" cells via 5-hydroxytryptamine receptor 3A (HTR3A), and the aggressive HTR3Ahigh human neuroblastoma cell lines reduce proliferation in response to HTR3A-specific agonists. In embryos (in vivo), the physiological increase of 5HT caused a prolongation of the cell cycle in "bridge" progenitors leading to a smaller chromaffin population and changing the balance of hormones and behavioral patterns in adulthood. These behavioral effects and smaller adrenals were mirrored in the progeny of pregnant female mice subjected to experimental stress, suggesting a maternal-fetal link that controls developmental adaptations. Finally, these results corresponded to a size-distribution of adrenals found in wild rodents with different coping strategies.

Correction: Targeting focal adhesion kinase overcomes erlotinib resistance in smoke induced lung cancer by altering phosphorylation of epidermal growth factor receptor.

[This corrects the article DOI: 10.18632/oncoscience.395.].

Single-cell transcriptomics of human embryos identifies multiple sympathoblast lineages with potential implications for neuroblastoma origin.

Characterization of the progression of cellular states during human embryogenesis can provide insights into the origin of pediatric diseases. We examined the transcriptional states of neural crest- and mesoderm-derived lineages differentiating into adrenal glands, kidneys, endothelium and hematopoietic tissue between post-conception weeks 6 and 14 of human development. Our results reveal transitions connecting the intermediate mesoderm and progenitors of organ primordia, the hematopoietic system and endothelial subtypes. Unexpectedly, by using a combination of single-cell transcriptomics and lineage tracing, we found that intra-adrenal sympathoblasts at that stage are directly derived from nerve-associated Schwann cell precursors, similarly to local chromaffin cells, whereas the majority of extra-adrenal sympathoblasts arise from the migratory neural crest. In humans, this process persists during several weeks of development within the large intra-adrenal ganglia-like structures, which may also serve as reservoirs of originating cells in neuroblastoma.

VHL inactivation without hypoxia is sufficient to achieve genome hypermethylation.

VHL inactivation is a key oncogenic event for renal carcinomas. In normoxia, VHL suppresses HIF1a-mediated transcriptional response, which is characteristic to hypoxia. It has previously been shown that hypoxic conditions inhibit TET-dependent hydroxymethylation of cytosines and cause DNA hypermethylation at gene promoters. In this work, we performed VHL inactivation by CRISPR/Cas9 and studied its effects on gene expression and DNA methylation. We showed that even without hypoxia, VHL inactivation leads to hypermethylation of the genome. Hypermethylated cytosines were evenly distributed throughout the genome with a slight preference for AP-1 (JUN and FOS) binding sites. Hypermethylated cytosines tended to be enriched within the binding sites of transcription factors that showed increased gene expression after VHL inactivation. We also observed promoter hypermethylation associated with decreased gene expression for several regulators of transcription and DNA methylation including SALL3.

Bifunctional immune checkpoint-targeted antibody-ligand traps that simultaneously disable TGFß enhance the efficacy of cancer immunotherapy.

A majority of cancers fail to respond to immunotherapy with antibodies targeting immune checkpoints, such as cytotoxic T-lymphocyte antigen-4 (CTLA-4) or programmed death-1 (PD-1)/PD-1 ligand (PD-L1). Cancers frequently express transforming growth factor-ß (TGFß), which drives immune dysfunction in the tumor microenvironment by inducing regulatory T cells (Tregs) and inhibiting CD8+ and TH1 cells. To address this therapeutic challenge, we invent bifunctional antibody-ligand traps (Y-traps) comprising an antibody targeting CTLA-4 or PD-L1 fused to a TGFß receptor II ectodomain sequence that simultaneously disables autocrine/paracrine TGFß in the target cell microenvironment (a-CTLA4-TGFßRIIecd and a-PDL1-TGFßRIIecd). a-CTLA4-TGFßRIIecd is more effective in reducing tumor-infiltrating Tregs and inhibiting tumor progression compared with CTLA-4 antibody (Ipilimumab). Likewise, a-PDL1-TGFßRIIecd exhibits superior antitumor efficacy compared with PD-L1 antibodies (Atezolizumab or Avelumab). Our data demonstrate that Y-traps counteract TGFß-mediated differentiation of Tregs and immune tolerance, thereby providing a potentially more effective immunotherapeutic strategy against cancers that are resistant to current immune checkpoint inhibitors.

Use of deep neural network ensembles to identify embryonic-fetal transition markers: repression of COX7A1 in embryonic and cancer cells.

Here we present the application of deep neural network (DNN) ensembles trained on transcriptomic data to identify the novel markers associated with the mammalian embryonic-fetal transition (EFT). Molecular markers of this process could provide important insights into regulatory mechanisms of normal development, epimorphic tissue regeneration and cancer. Subsequent analysis of the most significant genes behind the DNNs classifier on an independent dataset of adult-derived and human embryonic stem cell (hESC)-derived progenitor cell lines led to the identification of COX7A1 gene as a potential EFT marker. COX7A1, encoding a cytochrome C oxidase subunit, was up-regulated in post-EFT murine and human cells including adult stem cells, but was not expressed in pre-EFT pluripotent embryonic stem cells or their in vitro-derived progeny. COX7A1 expression level was observed to be undetectable or low in multiple sarcoma and carcinoma cell lines as compared to normal controls. The knockout of the gene in mice led to a marked glycolytic shift reminiscent of the Warburg effect that occurs in cancer cells. The DNN approach facilitated the elucidation of a potentially new biomarker of cancer and pre-EFT cells, the embryo-onco phenotype, which may potentially be used as a target for controlling the embryonic-fetal transition.

Targeting focal adhesion kinase overcomes erlotinib resistance in smoke induced lung cancer by altering phosphorylation of epidermal growth factor receptor.

EGFR-based targeted therapies have shown limited success in smokers. Identification of alternate signaling mechanism(s) leading to TKI resistance in smokers is critically important. We observed increased resistance to erlotinib in H358 NSCLC (non-small cell lung carcinoma) cells chronically exposed to cigarette smoke (H358-S) compared to parental cells. SILAC-based mass-spectrometry approach was used to study altered signaling in H358-S cell line. Importantly, among the top phosphosites in H358-S cells we observed hyperphosphorylation of EGFR (Y1197) and non-receptor tyrosine kinase FAK (Y576/577). Supporting these observations, a transcriptomic-based pathway activation analysis of TCGA NSCLC datasets revealed that FAK and EGFR internalization pathways were significantly upregulated in smoking patients, compared to the never-smokers and were associated with elevated PI3K signaling and lower level of caspase cascade and E-cadherin pathways activation. We show that inhibition of FAK led to decreased cellular proliferation and invasive ability of the smoke-exposed cells, and restored their dependency on EGFR signaling. Our data suggests that activation of focal adhesion pathway significantly contributes to erlotinib resistance, and that FAK is a potential therapeutic target for management of erlotinib resistance in smoke-induced NSCLC.