In vivo genome-wide analysis of multiple tissues identifies gene regulatory networks, novel functions and downstream regulatory genes for Bapx1 and its co-regulation with Sox9 in the mammalian vertebral column

RET enhancer haplotype-dependent remodeling of the human fetal gut development program

Hirschsprung disease (HSCR) is associated with deficiency of the receptor tyrosine kinase RET, resulting in loss of cells of the enteric nervous system (ENS) during fetal gut development. The major contribution to HSCR risk is from common sequence variants in RET enhancers with additional risk from rare coding variants in many genes. Here, we demonstrate that these RET enhancer variants specifically alter the human fetal gut development program through significant decreases in gene expression of RET, members of the RET-EDNRB gene regulatory network (GRN), other HSCR genes, with an altered transcriptome of 2,382 differentially expressed genes across diverse neuronal and mesenchymal functions. A parsimonious hypothesis for these results is that beyond RET's direct effect on its GRN, it also has a major role in enteric neural crest-derived cell (ENCDC) precursor proliferation, its deficiency reducing ENCDCs with relative expansion of non-ENCDC cells. Thus, genes reducing RET proliferative activity can potentially cause HSCR. One such class is the 23 RET-dependent transcription factors enriched in early gut development. We show that their knockdown in human neuroblastoma SK-N-SH cells reduces RET and/or EDNRB gene expression, expanding the RET-EDNRB GRN. The human embryos we studied had major remodeling of the gut transcriptome but were unlikely to have had HSCR: thus, genetic or epigenetic changes in addition to those in RET are required for aganglionosis.

Intervertebral disc degeneration-Current therapeutic options and challenges

Degeneration of the intervertebral disc (IVD) is a normal part of aging. Due to the spine's declining function and the development of pain, it may affect one's physical health, mental health, and socioeconomic status. Most of the intervertebral disc degeneration (IVDD) therapies today focus on the symptoms of low back pain rather than the underlying etiology or mechanical function of the disc. The deteriorated disc is typically not restored by conservative or surgical therapies that largely focus on correcting symptoms and structural abnormalities. To enhance the clinical outcome and the quality of life of a patient, several therapeutic modalities have been created. In this review, we discuss genetic and environmental causes of IVDD and describe promising modern endogenous and exogenous therapeutic approaches including their applicability and relevance to the degeneration process.

Glis1 and oxaloacetate in nucleus pulposus stromal cell somatic reprogramming and survival

Regenerative medicine aims to repair degenerate tissue through cell refurbishment with minimally invasive procedures. Adipose tissue (FAT)-derived stem or stromal cells are a convenient autologous choice for many regenerative cell therapy approaches. The intervertebral disc (IVD) is a suitable target. Comprised of an inner nucleus pulposus (NP) and an outer annulus fibrosus (AF), the degeneration of the IVD through trauma or aging presents a substantial socio-economic burden worldwide. The avascular nature of the mature NP forces cells to reside in a unique environment with increased lactate levels, conditions that pose a challenge to cell-based therapies. We assessed adipose and IVD tissue-derived stromal cells through in vitro transcriptome analysis in 2D and 3D culture and suggested that the transcription factor Glis1 and metabolite oxaloacetic acid (OAA) could provide NP cells with survival tools for the harsh niche conditions in the IVD.

Development, Pathogenesis, and Regeneration of the Intervertebral Disc: Current and Future Insights Spanning Traditional to Omics Methods

The intervertebral disc (IVD) is the fibrocartilaginous joint located between each vertebral body that confers flexibility and weight bearing capabilities to the spine. The IVD plays an important role in absorbing shock and stress applied to the spine, which helps to protect not only the vertebral bones, but also the brain and the rest of the central nervous system. Degeneration of the IVD is correlated with back pain, which can be debilitating and severely affects quality of life. Indeed, back pain results in substantial socioeconomic losses and healthcare costs globally each year, with about 85% of the world population experiencing back pain at some point in their lifetimes. Currently, therapeutic strategies for treating IVD degeneration are limited, and as such, there is great interest in advancing treatments for back pain. Ideally, treatments for back pain would restore native structure and thereby function to the degenerated IVD. However, the complex developmental origin and tissue composition of the IVD along with the avascular nature of the mature disc makes regeneration of the IVD a uniquely challenging task. Investigators across the field of IVD research have been working to elucidate the mechanisms behind the formation of this multifaceted structure, which may identify new therapeutic targets and inform development of novel regenerative strategies. This review summarizes current knowledge base on IVD development, degeneration, and regenerative strategies taken from traditional genetic approaches and omics studies and discusses the future landscape of investigations in IVD research and advancement of clinical therapies.

Low-pass whole genome bisulfite sequencing of neonatal dried blood spots identifies a role for RUNX1 in Down syndrome DNA methylation profiles

Neonatal dried blood spots (NDBS) are a widely banked sample source that enables retrospective investigation into early life molecular events. Here, we performed low-pass whole genome bisulfite sequencing (WGBS) of 86 NDBS DNA to examine early life Down syndrome (DS) DNA methylation profiles. DS represents an example of genetics shaping epigenetics, as multiple array-based studies have demonstrated that trisomy 21 is characterized by genome-wide alterations to DNA methylation. By assaying over 24 million CpG sites, thousands of genome-wide significant (q < 0.05) differentially methylated regions (DMRs) that distinguished DS from typical development and idiopathic developmental delay were identified. Machine learning feature selection refined these DMRs to 22 loci. The DS DMRs mapped to genes involved in neurodevelopment, metabolism, and transcriptional regulation. Based on comparisons with previous DS methylation studies and reference epigenomes, the hypermethylated DS DMRs were significantly (q < 0.05) enriched across tissues while the hypomethylated DS DMRs were significantly (q < 0.05) enriched for blood-specific chromatin states. A ~28 kb block of hypermethylation was observed on chromosome 21 in the RUNX1 locus, which encodes a hematopoietic transcription factor whose binding motif was the most significantly enriched (q < 0.05) overall and specifically within the hypomethylated DMRs. Finally, we also identified DMRs that distinguished DS NDBS based on the presence or absence of congenital heart disease (CHD). Together, these results not only demonstrate the utility of low-pass WGBS on NDBS samples for epigenome-wide association studies, but also provide new insights into the early life mechanisms of epigenomic dysregulation resulting from trisomy 21.

Simultaneous detection of multiple mRNAs and proteins in bovine IVD cells and tissue with single cell resolution

OBJECTIVES

Interactions of cells with their neighbors and influences by the surrounding extracellular matrix (ECM) is reflected in a cells transcriptome and proteome. In tissues comprised of heterogeneous cell populations or cells depending on ECM signalling cues such as those of the intervertebral disc (IVD), this information is obscured or lost when cells are pooled for the commonly used transcript analysis by quantitative PCR or RNA sequencing. Instead, these cells require means to analyse RNA transcript and protein distribution at a single cell or subcellular level to identify different cell types and functions, without removing them from their surrounding signalling cues.

RESULTS

We developed a simple, sequential protocol combining RNA is situ hybridisation (RISH) and immunohistochemistry (IHC) for the simultaneous analysis of multiple transcripts alongside proteins. This allows one to characterize heterogeneous cell populations at the single cell level in the natural cell environment and signalling context, both in vivo and in vitro. This protocol is demonstrated on cells of the bovine IVD, for transcripts and proteins involved in mechanotransduction, stemness and cell proliferation.

CONCLUSIONS

A simple, sequential protocol combining RISH and IHC is presented that allows for simultaneous information on RNA transcripts and proteins to characterize cells within a heterogeneous cell population and complex signalling environments such as those of the IVD.

Quantitative Single-Cell Transcript Assessment of Biomarkers Supports Cellular Heterogeneity in the Bovine IVD

Severe and chronic low back pain is often associated with intervertebral disc (IVD) degeneration. While imposing a considerable socio-economic burden worldwide, IVD degeneration is also severely impacting on the quality of life of affected individuals. Cell-based regenerative medicine approaches have moved into clinical trials, yet IVD cell identities in the mature disc remain to be fully elucidated and tissue heterogeneity exists, requiring a better characterization of IVD cells. The bovine coccygeal IVD is an accepted research model to study IVD mechano-biology and disc homeostasis. Recently, we identified novel IVD biomarkers in the outer annulus fibrosus (AF) and nucleus pulposus (NP) of the mature bovine coccygeal IVD through RNA in situ hybridization (AP-RISH) and z-proportion test. Here we follow up on Lam1, Thy1, Gli1, Gli3, Noto, Ptprc, Scx, Sox2 and Zscan10 with fluorescent RNA in situ hybridization (FL-RISH) and confocal microscopy. This permits sub-cellular transcript localization and the addition of quantitative single-cell derived values of mRNA expression levels to our previous analysis. Lastly, we used a Gaussian mixture modeling approach for the exploratory analysis of IVD cells. This work complements our earlier cell population proportion-based study, confirms the previously proposed biomarkers and indicates even further heterogeneity of cells in the outer AF and NP of a mature IVD.

Klhl14 Antisense RNA is a Target of Key Skeletogenic Transcription Factors in the Developing Intervertebral Disc

STUDY DESIGN

RNA in situ hybridization (RISH) allows for validation and characterization of the long noncoding (lnc) natural antisense RNA (NAT) Klhl14as in the embryonic murine intervertebral disc (IVD) in the context of loss-of-function mutants for key transcription factors (TFs) in axial skeleton development.

OBJECTIVE

Validation of Klhl14as in the developing murine IVD.

SUMMARY OF BACKGROUND DATA

The IVD is a focus of regenerative medicine; however, processes and signaling cascades resulting in the different cell types in a mature IVD still require clarification in most animals including humans. Technological advances increasingly point to implications of lnc NATs in transcription/translation regulation. Transcriptome data generation and analysis identified a protein encoding transcript and related noncoding antisense transcript as downregulated in embryos devoid of key TFs during axial skeleton development. Here, primarily, the antisense transcript is analyzed in this loss-of-function context.

METHODS

4930426D05Rik and 6330403N15Rik were identified as Klhl14as and sense, respectively, two transcripts downregulated in the vertebral column of midgestation Pax1 and Pax9 mutant mouse embryos. RISH on wildtype and mutant embryos for the TF encoding genes Pax1/Pax9, Sox5/Sox6/Sox9, and Bapx1 was used to further analyze Klhl14as in the developing IVD.

RESULTS

Klhl14as and Klhl14 were the top downregulated transcripts in Pax1; Pax9 E12.5 embryos. Our data demonstrate expression of Klhl14as and sense transcripts in the annulus fibrosus (AF) and notochord of the developing IVD. Klhl14as expression in the inner annulus fibrosus (iAF) seems dependent on the TFs Pax1/Pax9, Sox6, Sox9, and Bapx1.

CONCLUSION

We are the first to suggest a role for the lncRNA Klhl14as in the developing IVD. Our data link Klhl14as to a previously established gene regulatory network during axial skeleton development and contribute further evidence that lnc NATs are involved in crucial gene regulatory networks in eukaryotic cells.

LEVEL OF EVIDENCE

N/A.

Potential biomarkers of the mature intervertebral disc identified at the single cell level

Intervertebral disc (IVD) degeneration and trauma is a major socio-economic burden and the focus of cell-based regenerative medicine approaches. Despite numerous ongoing clinical trials attempting to replace ailing IVD cells with mesenchymal stem cells, a solid understanding of the identity and nature of cells in a healthy mature IVD is still in need of refinement. Although anatomically simple, the IVD is comprised of heterogeneous cell populations. Therefore, methods involving cell pooling for RNA profiling could be misleading. Here, by using RNA in situ hybridization and z proportion test, we have identified potential novel biomarkers through single cell assessment. We quantified the proportion of RNA transcribing cells for 50 genetic loci in the outer annulus fibrosus (AF) and nucleus pulposus (NP) in coccygeal bovine discs isolated from tails of four skeletally mature animals. Our data reconfirm existing data and suggest 10 novel markers such as Lam1 and Thy1 in the outer AF and Gli1, Gli3, Noto, Scx, Ptprc, Sox2, Zscan10 and LOC101904175 in the NP, including pluripotency markers, that indicate stemness potential of IVD cells. These markers could be added to existing biomarker panels for cell type characterization. Furthermore, our data once more demonstrate heterogeneity in cells of the AF and NP, indicating the need for single cell assessment by methods such as RNA in situ hybridization. Our work refines the molecular identity of outer AF and NP cells, which can benefit future regenerative medicine and tissue engineering strategies in humans.

RNA in situ hybridization characterization of non-enzymatic derived bovine intervertebral disc cell lineages suggests progenitor cell potential

Degeneration of the intervertebral disc (IVD) is a meritorious target for therapeutic cell based regenerative medicine approaches, however, controversy over what defines the precise identity of mature IVD cells and lack of single cell based quality control measures is of concern. Bos taurus and human IVDs are histologically more similar than is Mus musculus. The mature bovine IVD is well suited as model system for technology development to be translated into therapeutic cell based regenerative medicine applications. We present a reproducible non-enzymatic protocol to isolate cell progenitor populations of three distinct areas of the mature bovine IVD. Bovine specific RNA probes were validated in situ and employed to assess fate changes, heterogeneity, stem cell characteristics and differentiation potential of the cultures. Quality control measures with single cell resolution like RNA in situ hybridization to assess culture heterogeneity (PISH) followed by optimization of culture conditions could be translated to human IVD cell culture to increase the safety of cell based regenerative medicine.

Genome wide binding (ChIP-Seq) of murine Bapx1 and Sox9 proteins in vivo and in vitro

This work pertains to GEO submission GSE36672, in vivo and in vitro genome wide binding (ChIP-Seq) of Bapx1/Nkx3.2 and Sox9 proteins. We have previously shown that data from a genome wide binding assay combined with transcriptional profiling is an insightful means to divulge the mechanisms directing cell type specification and the generation of tissues and subsequent organs [1]. Our earlier work identified the role of the DNA-binding homeodomain containing protein Bapx1/Nkx3.2 in midgestation murine embryos. Microarray analysis of EGFP-tagged cells (both wildtype and null) was integrated using ChIP-Seq analysis of Bapx1/Nkx3.2 and Sox9 DNA-binding proteins in living tissue.

Revisiting the role of the Gcm transcription factor, from master regulator to Swiss army knife

Master genes are known to induce the differentiation of a multipotent cell into a specific cell type. These molecules are often transcription factors that switch on the regulatory cascade that triggers cell specification. Gcm was first described as the master gene of the glial fate in Drosophila as it induces the differentiation of neuroblasts into glia in the developing nervous system. Later on, Gcm was also shown to regulate the differentiation of blood, tendon and peritracheal cells as well as that of neuronal subsets. Thus, the glial master gene is used in at least 4 additional systems to promote differentiation. To understand the numerous roles of Gcm, we recently reported a genome-wide screen of Gcm direct targets in the Drosophila embryo. This screen provided new insight into the role and mode of action of this powerful transcription factor, notably on the interactions between Gcm and major differentiation pathways such as the Hedgehog, Notch and JAK/STAT. Here, we discuss the mode of action of Gcm in the different systems, we present new tissues that require Gcm and we revise the concept of 'master gene'.

Gene expression profiles of Bapx1 expressing FACS sorted cells from wildtype and Bapx1-EGFP null mouse embryos

The data described in this article refers to Chatterjee et al. (2015) “In vivo genome-wide analysis of multiple tissues identifies gene regulatory networks, novel functions and downstream regulatory genes for Bapx1 and its co-regulation with Sox9 in the mammalian vertebral column” (GEO GSE35649) [1]. Transcriptional profiling combined with genome wide binding data is a powerful tool to elucidate the molecular mechanism behind vertebrate organogenesis. It also helps to uncover multiple roles of a single gene in different organs. In the above mentioned report we reveal the function of the homeobox gene Bapx1 during the embryogenesis of five distinct organs (vertebral column, spleen, gut, forelimb and hindlimb) at a relevant developmental stage (E12.5), microarray analysis of isolated wildtype and mutant cells in is compared in conjunction with ChIP-Seq analysis. We also analyzed the development of the vertebral column by comparing microarray and ChIP-Seq data for Bapx1 with similarly generated data sets for Sox9 to generate a gene regulatory network controlling various facets of the organogenesis.

Plate in situ hybridization (PISH) as a time and cost effective RNA expression assay to study phenotypic heterogeneity in a population of cultured murine cells at single cell resolution

OBJECTIVES

Regenerative medicine approaches using reprogrammed or transdifferentiated cells require efficient single cell expression profiling to analyze culture homogeneity for quality control and recipients' safety.

RESULTS

While antigen-antibody based systems have been developed for several proteins, probing at the mRNA level allows for more flexibility, faster adaption to the ever increasing new data from next generation sequencing and increased specificity, especially for genes of conserved gene families.

CONCLUSIONS

We developed a time and cost effective expression profiling assay for monolayer cell culture in 96-well plates based on RNA in situ hybridization, termed PISH, at single cell resolution.