Identification of novel cerebellar developmental transcriptional regulators with motif activity analysis.

BACKGROUND: The work of the FANTOM5 Consortium has brought forth a new level of understanding of the regulation of gene transcription and the cellular processes involved in creating diversity of cell types. In this study, we extended the analysis of the FANTOM5 Cap Analysis of Gene Expression (CAGE) transcriptome data to focus on understanding the genetic regulators involved in mouse cerebellar development. RESULTS: We used the HeliScopeCAGE library sequencing on cerebellar samples over 8 embryonic and 4 early postnatal times. This study showcases temporal expression pattern changes during cerebellar development. Through a bioinformatics analysis that focused on transcription factors, their promoters and binding sites, we identified genes that appear as strong candidates for involvement in cerebellar development. We selected several candidate transcriptional regulators for validation experiments including qRT-PCR and shRNA transcript knockdown. We observed marked and reproducible developmental defects in Atf4, Rfx3, and Scrt2 knockdown embryos, which support the role of these genes in cerebellar development. CONCLUSIONS: The successful identification of these novel gene regulators in cerebellar development demonstrates that the FANTOM5 cerebellum time series is a high-quality transcriptome database for functional investigation of gene regulatory networks in cerebellar development.

Correction to: Relatively frequent switching of transcription start sites during cerebellar development.

CORRECTION: The authors of the original article [1] would like to recognize the critical contribution of core members of the FANTOM5 Consortium, who played the critical role of HeliScopeCAGE sequencing experiments, quality control of tag reads and processing of the raw sequencing data.

Discovery of Transcription Factors Novel to Mouse Cerebellar Granule Cell Development Through Laser-Capture Microdissection.

Laser-capture microdissection was used to isolate external germinal layer tissue from three developmental periods of mouse cerebellar development: embryonic days 13, 15, and 18. The cerebellar granule cell-enriched mRNA library was generated with next-generation sequencing using the Helicos technology. Our objective was to discover transcriptional regulators that could be important for the development of cerebellar granule cells-the most numerous neuron in the central nervous system. Through differential expression analysis, we have identified 82 differentially expressed transcription factors (TFs) from a total of 1311 differentially expressed genes. In addition, with TF-binding sequence analysis, we have identified 46 TF candidates that could be key regulators responsible for the variation in the granule cell transcriptome between developmental stages. Altogether, we identified 125 potential TFs (82 from differential expression analysis, 46 from motif analysis with 3 overlaps in the two sets). From this gene set, 37 TFs are considered novel due to the lack of previous knowledge about their roles in cerebellar development. The results from transcriptome-wide analyses were validated with existing online databases, qRT-PCR, and in situ hybridization. This study provides an initial insight into the TFs of cerebellar granule cells that might be important for development and provide valuable information for further functional studies on these transcriptional regulators.

A Novel and Multivalent Role of Pax6 in Cerebellar Development.

UNLABELLED: Pax6 is a prominent gene in brain development. The deletion of Pax6 results in devastated development of eye, olfactory bulb, and cortex. However, it has been reported that the Pax6-null Sey cerebellum only has minor defects involving granule cells despite Pax6 being expressed throughout cerebellar development. The present work has uncovered a requirement of Pax6 in the development of all rhombic lip (RL) lineages. A significant downregulation of Tbr1 and Tbr2 expression is found in the Sey cerebellum, these are cell-specific markers of cerebellar nuclear (CN) neurons and unipolar brush cells (UBCs), respectively. The examination of Tbr1 and Lmx1a immunolabeling and Nissl staining confirmed the loss of CN neurons from the Sey cerebellum. CN neuron progenitors are produced in the mutant but there is an enhanced death of these neurons as shown by increased presence of caspase-3-positive cells. These data indicate that Pax6 regulates the survival of CN neuron progenitors. Furthermore, the analysis of experimental mouse chimeras suggests a cell-extrinsic role of Pax6 in CN neuron survival. For UBCs, using Tbr2 immunolabeling, these cells are significantly reduced in the Sey cerebellum. The loss of UBCs in the mutant is due partly to cell death in the RL and also to the reduced production of progenitors from the RL. These results demonstrate a critical role for Pax6 in regulating the generation and survival of UBCs. This and previous work from our laboratory demonstrate a seminal role of Pax6 in the development of all cerebellar glutamatergic neurons. SIGNIFICANCE STATEMENT: Pax6 is a key molecule in development. Pax6 is best known as the master control gene in eye development with mutations causing aniridia in humans. Pax6 also plays important developmental roles in the cortex and olfactory bulb. During cerebellar development, Pax6 is robustly expressed in the germinal zone of all glutamatergic neurons [cerebellar nuclear (CN) neurons, granule cells, and unipolar brush cells (UBCs)]. Past work has not found abnormalities in the CN and UBC populations. Our study reveals that the Pax6-null mutation dramatically affects these cells and identifies Pax6 as a key regulator of cell survival in CN neurons and of cell production in UBCs. The present study shows how Pax6 is key to the development of glutamatergic cells in the cerebellum.

Relatively frequent switching of transcription start sites during cerebellar development.

BACKGROUND: Alternative transcription start site (TSS) usage plays important roles in transcriptional control of mammalian gene expression. The growing interest in alternative TSSs and their role in genome diversification spawned many single-gene studies on differential usages of tissue-specific or temporal-specific alternative TSSs. However, exploration of the switching usage of alternative TSS usage on a genomic level, especially in the central nervous system, is largely lacking. RESULTS: In this study, We have prepared a unique set of time-course data for the developing cerebellum, as part of the FANTOM5 consortium ( http://fantom.gsc.riken.jp/5/ ) that uses their innovative capturing of 5' ends of all transcripts followed by Helicos next generation sequencing. We analyzed the usage of all transcription start sites (TSSs) at each time point during cerebellar development that provided information on multiple RNA isoforms that emerged from the same gene. We developed a mathematical method that systematically compares the expression of different TSSs of a gene to identify temporal crossover and non-crossover switching events. We identified 48,489 novel TSS switching events in 5433 genes during cerebellar development. This includes 9767 crossover TSS switching events in 1511 genes, where the dominant TSS shifts over time. CONCLUSIONS: We observed a relatively high prevalence of TSS switching in cerebellar development where the resulting temporally-specific gene transcripts and protein products can play important regulatory and functional roles.

CbGRiTS: cerebellar gene regulation in time and space.

The mammalian CNS is one of the most complex biological systems to understand at the molecular level. The temporal information from time series transcriptome analysis can serve as a potent source of associative information between developmental processes and regulatory genes. Here, we introduce a new transcriptome database called, Cerebellar Gene Regulation in Time and Space (CbGRiTS). This dataset is populated with transcriptome data across embryonic and postnatal development from two standard mouse strains, C57BL/6J and DBA/2J, several recombinant inbred lines and cerebellar mutant strains. Users can evaluate expression profiles across cerebellar development in a deep time series with graphical interfaces for data exploration and link-out to anatomical expression databases. We present three analytical approaches that take advantage of specific aspects of the time series for transcriptome analysis. We demonstrate the use of CbGRiTS dataset as a community resource to explore patterns of gene expression and develop hypotheses concerning gene regulatory networks in brain development.

Wls provides a new compartmental view of the rhombic lip in mouse cerebellar development.

Math1 is the defining molecule of the cerebellar rhombic lip and Pax6 is downstream in the Math1 pathway. In the present study, we discover that Wntless (Wls) is a novel molecular marker of the cells in the interior face of the rhombic lip throughout normal mouse cerebellar development. Wls expression is found complementary to the expression of Math1 and Pax6, which are localized to the exterior face of the rhombic lip. To determine the interaction between these molecules, we examine the loss-of-Math1 or loss-of-Pax6 in the cerebellum, i.e., the Math1-null and Pax6-null (Sey) mutant cerebella. The presence of Wls-positive cells in the Math1-null rhombic lip indicates that Wls expression is independent of Math1. In the Sey mutant cerebellum, there is an expansion of Wls-expressing cells into regions that are normally colonized by Pax6-expressing cells. The ectopic expression of Wls in the Pax6-null cerebellum suggests a negative interaction between Wls-expressing cells and Pax6-positive cells. These findings suggest that the rhombic lip is dynamically patterned by the expression of Wls, Math1, and Pax6. We also examine five rhombic lip cell markers (Wls, Math1, Pax6, Lmx1a, and Tbr2) to identify four molecularly distinct compartments in the rhombic lip during cerebellar development. The existence of spatial compartmentation in the rhombic lip and the interplay between Wls, Math1, and Pax6 in the rhombic lip provides novel views of early cerebellar development.

CAGE-defined promoter regions of the genes implicated in Rett Syndrome.

BACKGROUND: Mutations in three functionally diverse genes cause Rett Syndrome. Although the functions of Forkhead box G1 (FOXG1), Methyl CpG binding protein 2 (MECP2) and Cyclin-dependent kinase-like 5 (CDKL5) have been studied individually, not much is known about their relation to each other with respect to expression levels and regulatory regions. Here we analyzed data from hundreds of mouse and human samples included in the FANTOM5 project, to identify transcript initiation sites, expression levels, expression correlations and regulatory regions of the three genes. RESULTS: Our investigations reveal the predominantly used transcription start sites (TSSs) for each gene including novel transcription start sites for FOXG1. We show that FOXG1 expression is poorly correlated with the expression of MECP2 and CDKL5. We identify promoter shapes for each TSS, the predicted location of enhancers for each gene and the common transcription factors likely to regulate the three genes. Our data imply Polycomb Repressive Complex 2 (PRC2) mediated silencing of Foxg1 in cerebellum. CONCLUSIONS: Our analyses provide a comprehensive picture of the regulatory regions of the three genes involved in Rett Syndrome.

Analytical and Clinical Validation of the Oncomine Dx Target Test to Assess HER2 Mutation Status in Tumor Tissue Samples From Patients With Non-Small Cell Lung Cancer Treated With Trastuzumab Deruxtecan in the DESTINY-Lung01 and DESTINY-Lung02 Studies.

CONTEXT.­: Trastuzumab deruxtecan (T-DXd), a human epidermal growth factor receptor 2 (HER2)-targeted therapy, has demonstrated durable anticancer activity in patients with advanced, metastatic HER2 (also known as ERBB2)-mutant (HER2m) non-small cell lung cancer (NSCLC) who have limited treatment options and poor prognosis. OBJECTIVE.­: To analytically validate and assess the clinical utility of the Oncomine Dx Target (ODxT) Test as a companion diagnostic to identify patients with HER2m NSCLC. DESIGN.­: Tumor samples from patients in DESTINY-Lung01 and DESTINY-Lung02 were retrospectively analyzed alongside commercially procured samples using the ODxT test and compared to the assays used for screening in these clinical trials. RESULTS.­: Positive percent agreement (PPA) and negative percent agreement (NPA) between the ODxT Test and TruSight Tumor 170 assay when testing DESTINY-Lung01 and commercially procured samples met prespecified thresholds (PPA and NPA ≥90%) for analytical accuracy (100% and 99.1%). The ODxT Test results were highly concordant with clinical trial assays (CTAs) used in DESTINY-Lung01 (PPA and NPA, 98.0% and 100%) and DESTINY-Lung02 (PPA and NPA, 96.7% and 100%) to identify activating HER2 mutations in tumor samples. Confirmed objective response rates were similar between patients with HER2m tumors identified by the ODxT Test and by CTAs in DESTINY-Lung01 (58.3% and 54.9%) and DESTINY-Lung02 (53.6% and 53.8%). Response duration was 12.0 and 9.3 months for patients identified by the ODxT Test and CTAs, respectively, in DESTINY-Lung01. CONCLUSIONS.­: The ODxT Test detected HER2 mutations in NSCLC with high analytical and clinical accuracy and identified HER2m populations with response rates similar to populations identified by CTAs, supporting clinical utility of the ODxT Test to inform treatment decisions for HER2m NSCLC.

Spatial and temporal requirements for huntingtin (Htt) in neuronal migration and survival during brain development.

Huntington's disease (HD), caused by an expanded triplet repeat in the huntingtin (Htt) gene, results in extensive neuropathology, but study of the Htt gene in CNS development through gene knockout is problematic as the knockout leads to embryonic lethality in mice. Here, we report that the knockdown of Htt expression in neuroepithelial cells of neocortex results in disturbed cell migration, reduced proliferation, and increased cell death that is relatively specific to early neural development. In the cerebellum, however, Htt knockdown results in cell death but not perturbed migration. The cell death phenotype in cortex can be partially reversed with co-knockdown of Casp9, indicating that mitochondria-mediated cell apoptotic processes are involved in the neuronal death. The timing of knockdown during early development is also an important variable. These results indicate a spatial and temporal requirement for Htt expression in neural development. Although it is uncertain whether the loss of wild-type huntingtin function contributes to pathogenesis in Huntington's disease, these results clearly contraindicate the use of nonspecific knockdown of Htt as a therapeutic measure in HD, particularly in utero.

Genome-wide microarray comparison reveals downstream genes of Pax6 in the developing mouse cerebellum.

The Pax6 transcription factor is expressed in cerebellar granule cells and when mutated, as in the Sey/Sey mouse, produces granule cells with disturbed survival and migration and with defects in neurite extension. The impact of Pax6 on other genes in the context of cerebellar development has not been identified. In this study, we performed transcriptome comparisons between wildtype and Pax6-null whole cerebellar tissue at embryonic day (E) 13.5, 15.5 and 18.5 using Affymetrix arrays (U74Av2). Statistical analyses identified 136 differentially regulated transcripts (FDR 0.05, 1.2-fold change cutoff) over time in Pax6-null cerebellar tissue. In parallel we examined the Math1-null granuloprival cerebellum and identified 228 down-regulated transcripts (FDR 0.05, 1.2-fold change cutoff). The intersection of these two microarray datasets produced a total of 21 differentially regulated transcripts. For a subset of the identified transcripts, we used qRT-PCR to validate the microarray data and demonstrated the expression in the rhombic lip lineage and differential expression in Pax6-null cerebellum with in situ hybridisation analysis. The candidate genes identified in this way represent direct or indirect Pax6-downstream genes involved in cerebellar development.

Conserved and differential gene interactions in dynamical biological systems.

MOTIVATION: While biological systems operated from a common genome can be conserved in various ways, they can also manifest highly diverse dynamics and functions. This is because the same set of genes can interact differentially across specific molecular contexts. For example, differential gene interactions give rise to various stages of morphogenesis during cerebellar development. However, after over a decade of efforts toward reverse engineering biological networks from high-throughput omic data, gene networks of most organisms remain sketchy. This hindrance has motivated us to develop comparative modeling to highlight conserved and differential gene interactions across experimental conditions, without reconstructing complete gene networks first. RESULTS: We established a comparative dynamical system modeling (CDSM) approach to identify conserved and differential interactions across molecular contexts. In CDSM, interactions are represented by ordinary differential equations and compared across conditions through statistical heterogeneity and homogeneity tests. CDSM demonstrated a consistent superiority over differential correlation and reconstruct-then-compare in simulation studies. We exploited CDSM to elucidate gene interactions important for cellular processes poorly understood during mouse cerebellar development. We generated hypotheses on 66 differential genetic interactions involved in expansion of the external granule layer. These interactions are implicated in cell cycle, differentiation, apoptosis and morphogenesis. Additional 1639 differential interactions among gene clusters were also identified when we compared gene interactions during the presence of Rhombic lip versus the presence of distinct internal granule layer. Moreover, compared with differential correlation and reconstruct-then-compare, CDSM makes fewer assumptions on data and thus is applicable to a wider range of biological assays. AVAILABILITY: Source code in C++ and R is available for non-commercial organizations upon request from the corresponding author. The cerebellum gene expression dataset used in this article is available upon request from the Goldowitz lab (dang@cmmt.ubc.ca, http://grits.dglab.org/). CONTACT: joemsong@cs.nmsu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Role of nitric oxide in classical conditioning of siphon withdrawal in Aplysia.

Nitric oxide (NO) is thought to be involved in several forms of learning in vivo and synaptic plasticity in vitro, but very little is known about the role of NO during physiological forms of plasticity that occur during learning. We addressed that question in a simplified preparation of the Aplysia siphon-withdrawal reflex. We first used in situ hybridization to show that the identified L29 facilitator neurons express NO synthase. Furthermore, exogenous NO produced facilitation of sensory-motor neuron EPSPs, and an inhibitor of NO synthase or an NO scavenger blocked behavioral conditioning. Application of the scavenger to the ganglion or injection into a sensory neuron blocked facilitation of the EPSP and changes in the sensory-neuron membrane properties during conditioning. Injection of the scavenger into the motor neuron reduced facilitation without affecting sensory neuron membrane properties, and injection of an inhibitor of NO synthase had no effect. Postsynaptic injection of an inhibitor of exocytosis had effects similar to injection of the scavenger. However, changes in the shape of the EPSP during conditioning were not consistent with postsynaptic AMPA-like receptor insertion but were mimicked by presynaptic spike broadening. These results suggest that NO makes an important contribution during conditioning and acts directly in both the sensory and motor neurons to affect different processes of facilitation at the synapses between them. In addition, they suggest that NO does not come from either the sensory or motor neurons but rather comes from another source, perhaps the L29 interneurons.

An integrated expression atlas of miRNAs and their promoters in human and mouse.

MicroRNAs (miRNAs) are short non-coding RNAs with key roles in cellular regulation. As part of the fifth edition of the Functional Annotation of Mammalian Genome (FANTOM5) project, we created an integrated expression atlas of miRNAs and their promoters by deep-sequencing 492 short RNA (sRNA) libraries, with matching Cap Analysis Gene Expression (CAGE) data, from 396 human and 47 mouse RNA samples. Promoters were identified for 1,357 human and 804 mouse miRNAs and showed strong sequence conservation between species. We also found that primary and mature miRNA expression levels were correlated, allowing us to use the primary miRNA measurements as a proxy for mature miRNA levels in a total of 1,829 human and 1,029 mouse CAGE libraries. We thus provide a broad atlas of miRNA expression and promoters in primary mammalian cells, establishing a foundation for detailed analysis of miRNA expression patterns and transcriptional control regions.

Terahertz pulsed spectroscopy of human Basal cell carcinoma.

Good contrast is seen between normal tissue and regions of tumor in terahertz pulsed imaging of basal cell carcinoma (BCC). To date, the source of contrast at terahertz frequencies is not well understood. In this paper we present results of a spectroscopy study comparing the terahertz properties (absorption coefficient and refractive index) of excised normal human skin and BCC. Both the absorption coefficient and refractive index were higher for skin that contained BCC. The difference was statistically significant over the range 0.2 to 2.0 THz (6.6 cm(-1) to 66.6 cm(-1)) for absorption coefficient and 0.25 to 0.90 THz (8.3 cm(-1) to 30 cm(-1)) for refractive index. The maximum difference for absorption was at 0.5 THz(16.7 cm(-1)). These changes are consistent with higher water content. These results account for the contrast seen in terahertz images of BCC and explain why parameters relating to the reflected terahertz pulse provide information about the lateral spread of the tumor. Knowing the properties of the tissue over the terahertz frequency range will enable the use of mathematical models to improve understanding of the terahertz response of normal and diseased tissue.

Molecular characterization of NMDA-like receptors in Aplysia and Lymnaea: relevance to memory mechanisms.

The N-methyl-D-aspartate (NMDA) receptor belongs to the group of ionotropic glutamate receptors and has been implicated in synaptic plasticity, memory acquisition, and learning in both vertebrates and invertebrates, including molluscs. However, the molecular identity of NMDA-type receptors in molluscs remains unknown. Here, we cloned two NMDA-type receptors from the sea slug Aplysia californica, AcNR1-1 and AcNR1-2, as well as their homologs from the freshwater pulmonate snail Lymnaea stagnalis, LsNR1-1 and LsNR1-2. The cloned receptors contain a signal peptide, two extracellular segments with predicted binding sites for glycine and glutamate, three recognized transmembrane regions, and a fourth hydrophobic domain that makes a hairpin turn to form a pore-like structure. Phylogenetic analysis suggests that both the AcNR1s and LsNR1s belong to the NR1 subgroup of ionotrophic glutamate receptors. Our in situ hybridization data indicate highly abundant, but predominantly neuron-specific expression of molluscan NR1-type receptors in all central ganglia, including identified motor neurons in the buccal and abdominal ganglia as well as groups of mechanosensory cells. AcNR1 transcripts were detected extrasynaptically in the neurites of metacerebral cells of Aplysia. The widespread distribution of AcNR1 and LsNR1 transcripts also implies diverse functions, including their involvement in the organization of feeding, locomotory, and defensive behaviors.

Simulating the response of terahertz radiation to basal cell carcinoma using ex vivo spectroscopy measurements.

Studies of basal cell carcinoma using terahertz pulsed imaging have revealed a significant difference between regions of tumor and healthy tissue. These differences are manifested in the reflected pulse due to what is thought to be changes in refractive index and absorption. We present measurements of the refractive index and absorption coefficient of excised normal tissue and basal cell carcinoma using terahertz (THz) transmission spectroscopy. We extract Debye parameters from these data and enter them into a finite difference time domain simulation to predict the shape of the waveforms reflected off the normal tissue and basal cell carcinoma and compare them with published in vivo data. Simulating the interaction of terahertz radiation with normal and cancerous tissue is a key step toward understanding the origin of contrast in terahertz images of skin cancer.