Developmental roles of glomerular epithelial protein-1 in mice molar morphogenesis.

Glomerular epithelial protein-1 (Glepp1), a R3 subtype family of receptor-type protein tyrosine phosphatases, plays important role in the activation of Src family kinases and regulates cellular processes such as cell proliferation, differentiation, and apoptosis. In this study, we firstly examined the functional evaluation of Glepp1 in tooth development and morphogenesis. The precise expression level and developmental function of Glepp1 were examined by RT-qPCR, in situ hybridization, and loss and gain of functional study using a range of in vitro organ cultivation methods. Expression of Glepp1 was detected in the developing tooth germs in cap and bell stage of tooth development. Knocking down Glepp1 at E13 for 2 days showed the altered expression levels of tooth development-related signaling molecules, including Bmps, Dspp, Fgf4, Lef1, and Shh. Moreover, transient knock down of Glepp1 revealed alterations in cellular physiology, examined by the localization patterns of Ki67 and E-cadherin. Similarly, knocking down of Glepp1 showed disrupted enamel rod and interrod formation in 3-week renal transplanted teeth. In addition, due to attrition of odontoblastic layers, the expression signals of Dspp and the localization of NESTIN were almost not detected after knock down of Glepp1; however, their expressions were increased after Glepp1 overexpression. Thus, our results suggested that Glepp1 plays modulating roles during odontogenesis by regulating the expression levels of signaling molecules and cellular events to achieve the proper structural formation of hard tissue matrices in mice molar development.

RBFOX2-regulated TEAD1 alternative splicing plays a pivotal role in Hippo-YAP signaling.

Alternative pre-mRNA splicing is key to proteome diversity; however, the biological roles of alternative splicing (AS) in signaling pathways remain elusive. Here, we focus on TEA domain transcription factor 1 (TEAD1), a YAP binding factor in the Hippo signaling pathway. Public database analyses showed that expression of YAP-TEAD target genes negatively correlated with the expression of a TEAD1 isoform lacking exon 6 (TEAD1ΔE6) but did not correlate with overall TEAD1 expression. We confirmed that the transcriptional activity and oncogenic properties of the full-length TEAD1 isoform were greater than those of TEAD1ΔE6, with the difference in transcription related to YAP interaction. Furthermore, we showed that RNA-binding Fox-1 homolog 2 (RBFOX2) promoted the inclusion of TEAD1 exon 6 via binding to the conserved GCAUG element in the downstream intron. These results suggest a regulatory mechanism of RBFOX2-mediated TEAD1 AS and provide insight into AS-specific modulation of signaling pathways.

The KRAB zinc finger protein ZFP809 is required to initiate epigenetic silencing of endogenous retroviruses.

Retroviruses have been invading mammalian germlines for millions of years, accumulating in the form of endogenous retroviruses (ERVs) that account for nearly one-tenth of the mouse and human genomes. ERVs are epigenetically silenced during development, yet the cellular factors recognizing ERVs in a sequence-specific manner remain elusive. Here we demonstrate that ZFP809, a member of the Krüppel-associated box zinc finger protein (KRAB-ZFP) family, initiates the silencing of ERVs in a sequence-specific manner via recruitment of heterochromatin-inducing complexes. ZFP809 knockout mice display highly elevated levels of ZFP809-targeted ERVs in somatic tissues. ERV reactivation is accompanied by an epigenetic shift from repressive to active histone modifications but only slight destabilization of DNA methylation. Importantly, using conditional alleles and rescue experiments, we demonstrate that ZFP809 is required to initiate ERV silencing during embryonic development but becomes largely dispensable in somatic tissues. Finally, we show that the DNA-binding specificity of ZFP809 is evolutionarily conserved in the Muroidea superfamily of rodents and predates the endogenization of retroviruses presently targeted by ZFP809 in Mus musculus. In sum, these data provide compelling evidence that ZFP809 evolved to recognize foreign DNA and establish histone modification-based epigenetic silencing of ERVs.

Effects of Mitochondrial Transplantation on Transcriptomics in a Polymicrobial Sepsis Model.

Previously, we demonstrated that mitochondrial transplantation has beneficial effects in a polymicrobial sepsis model. However, the mechanism has not been fully investigated. Mitochondria have their own genes, and genomic changes in sepsis are an important issue in terms of pathophysiology, biomarkers, and therapeutic targets. To investigate the changes in transcriptomic features after mitochondrial transplantation in a polymicrobial sepsis model, we used a rat model of fecal slurry polymicrobial sepsis. Total RNA from splenocytes of sham-operated (SHAM, n = 10), sepsis-induced (SEPSIS, n = 7), and sepsis receiving mitochondrial transplantation (SEPSIS + MT, n = 8) samples was extracted and we conducted a comparative transcriptome-wide analysis between three groups. We also confirmed these results with qPCR. In terms of percentage of mitochondrial mapped reads, the SEPSIS + MT group had a significantly higher mapping ratio than the others. RT1-M2 and Cbln2 were identified as highly expressed in SEPSIS + MT compared with SEPSIS. Using SHAM expression levels as another control variable, we further identified six genes (Fxyd4, Apex2l1, Kctd4, 7SK, SNORD94, and SNORA53) that were highly expressed after sepsis induction and observed that their expression levels were attenuated by mitochondrial transplantation. Changes in transcriptomic features were identified after mitochondrial transplantation in sepsis. This might provide a hint for exploring the mechanism of mitochondrial transplantation in sepsis.

A machine learning model for classifying G-protein-coupled receptors as agonists or antagonists.

BACKGROUND: G-protein coupled receptors (GPCRs) sense and transmit extracellular signals into the intracellular machinery by regulating G proteins. GPCR malfunctions are associated with a variety of signaling-related diseases, including cancer and diabetes; at least a third of the marketed drugs target GPCRs. Thus, characterization of their signaling and regulatory mechanisms is crucial for the development of effective drugs. RESULTS: In this study, we developed a machine learning model to identify GPCR agonists and antagonists. We designed two-step prediction models: the first model identified the ligands binding to GPCRs and the second model classified the ligands as agonists or antagonists. Using 990 selected subset features from 5270 molecular descriptors calculated from 4590 ligands deposited in two drug databases, our model classified non-ligands, agonists, and antagonists of GPCRs, and achieved an area under the ROC curve (AUC) of 0.795, sensitivity of 0.716, specificity of 0.744, and accuracy of 0.733. In addition, we verified that 70% (44 out of 63) of FDA-approved GPCR-targeting drugs were correctly classified into their respective groups. CONCLUSIONS: Studies of ligand-GPCR interaction recognition are important for the characterization of drug action mechanisms. Our GPCR-ligand interaction prediction model can be employed in the pharmaceutical sciences for the efficient virtual screening of putative GPCR-binding agonists and antagonists.

Pulmonary fibrosis model using micro-CT analyzable human PSC-derived alveolar organoids containing alveolar macrophage-like cells.

Human lung organoids (hLOs) are useful for disease modelling and drug screening. However, a lack of immune cells in hLOs limits the recapitulation of in vivo cellular physiology. Here, we generated hLOs containing alveolar macrophage (AMφ)-like cells derived from pluripotent stem cells (PSC). To bridge hLOs with advanced human lung high-resolution X-ray computed tomography (CT), we acquired quantitative micro-CT images. Three hLO types were observed during differentiation. Among them, alveolar hLOs highly expressed not only lung epithelial cell markers but also AMφ-specific markers. Furthermore, CD68+ AMφ-like cells were spatially organized on the luminal epithelial surface of alveolar hLOs. Bleomycin-treated alveolar hLOs showed upregulated expression of fibrosis-related markers and extracellular matrix deposits in the alveolar sacs. Alveolar hLOs also showed structural alterations such as excessive tissue fraction under bleomycin treatment. Therefore, we suggest that micro-CT analyzable PSC-derived alveolar hLOs are a promising in vitro model to predict lung toxicity manifestations, including fibrosis.

Identification of Novel Genes for Cell Fusion during Osteoclast Formation.

Osteoclasts are derived from hematopoietic stem cells. Monocyte preosteoclasts obtain resorbing activity via cell-cell fusion to generate multinucleated cells. However, the mechanisms and molecules involved in the fusion process are poorly understood. In this study, we performed RNA sequencing with single nucleated cells (SNCs) and multinucleated cells (MNCs) to identify the fusion-specific genes. The SNCs and MNCs were isolated under the same conditions during osteoclastogenesis with the receptor activator of nuclear factor-κB ligand (RANKL) administration. Based on this analysis, the expression of seven genes was found to be significantly increased in MNCs but decreased in SNCs, compared to that in bone marrow-derived macrophages (BMMs). We then generated knockout macrophage cell lines using a CRISPR-Cas9 genome-editing tool to examine their function during osteoclastogenesis. Calcrl-, Marco-, or Ube3a-deficient cells could not develop multinucleated giant osteoclasts upon RANKL stimulation. However, Tmem26-deficient cells fused more efficiently than control cells. Our findings demonstrate that Calcrl, Marco, and Ube3a are novel determinants of osteoclastogenesis, especially with respect to cell fusion, and highlight potential targets for osteoporosis therapy.

Comparability of reference-based and reference-free transcriptome analysis approaches at the gene expression level.

BACKGROUND: Lately, high-throughput RNA sequencing has been extensively used to elucidate the transcriptome landscape and dynamics of cell types of different species. In particular, for most non-model organisms lacking complete reference genomes with high-quality annotation of genetic information, reference-free (RF) de novo transcriptome analyses, rather than reference-based (RB) approaches, are widely used, and RF analyses have substantially contributed toward understanding the mechanisms regulating key biological processes and functions. To date, numerous bioinformatics studies have been conducted for assessing the workflow, production rate, and completeness of transcriptome assemblies within and between RF and RB datasets. However, the degree of consistency and variability of results obtained by analyzing gene expression levels through these two different approaches have not been adequately documented. RESULTS: In the present study, we evaluated the differences in expression profiles obtained with RF and RB approaches and revealed that the former tends to be satisfactorily replaced by the latter with respect to transcriptome repertoires, as well as from a gene expression quantification perspective. In addition, we urge cautious interpretation of these findings. Several genes that are lowly expressed, have long coding sequences, or belong to large gene families must be validated carefully, whenever gene expression levels are calculated using the RF method. CONCLUSIONS: Our empirical results indicate important contributions toward addressing transcriptome-related biological questions in non-model organisms.

PIC-Me: paralogs and isoforms classifier based on machine-learning approaches.

BACKGROUND: Paralogs formed through gene duplication and isoforms formed through alternative splicing have been important processes for increasing protein diversity and maintaining cellular homeostasis. Despite their recognized importance and the advent of large-scale genomic and transcriptomic analyses, paradoxically, accurate annotations of all gene loci to allow the identification of paralogs and isoforms remain surprisingly incomplete. In particular, the global analysis of the transcriptome of a non-model organism for which there is no reference genome is especially challenging. RESULTS: To reliably discriminate between the paralogs and isoforms in RNA-seq data, we redefined the pre-existing sequence features (sequence similarity, inverse count of consecutive identical or non-identical blocks, and match-mismatch fraction) previously derived from full-length cDNAs and EST sequences and described newly discovered genomic and transcriptomic features (twilight zone of protein sequence alignment and expression level difference). In addition, the effectiveness and relevance of the proposed features were verified with two widely used support vector machine (SVM) and random forest (RF) models. From nine RNA-seq datasets, all AUC (area under the curve) scores of ROC (receiver operating characteristic) curves were over 0.9 in the RF model and significantly higher than those in the SVM model. CONCLUSIONS: In this study, using an RF model with five proposed RNA-seq features, we implemented our method called Paralogs and Isoforms Classifier based on Machine-learning approaches (PIC-Me) and showed that it outperformed an existing method. Finally, we envision that our tool will be a valuable computational resource for the genomics community to help with gene annotation and will aid in comparative transcriptomics and evolutionary genomics studies, especially those on non-model organisms.

Is Phylotranscriptomics as Reliable as Phylogenomics?

Phylogenomics, the study of phylogenetic relationships among taxa based on their genome sequences, has emerged as the preferred phylogenetic method because of the wealth of phylogenetic information contained in genome sequences. Genome sequencing, however, can be prohibitively expensive, especially for taxa with huge genomes and when many taxa need sequencing. Consequently, the less costly phylotranscriptomics has seen an increased use in recent years. Phylotranscriptomics reconstructs phylogenies using DNA sequences derived from transcriptomes, which are often orders of magnitude smaller than genomes. However, in the absence of corresponding genome sequences, comparative analyses of transcriptomes can be challenging and it is unclear whether phylotranscriptomics is as reliable as phylogenomics. Here, we respectively compare the phylogenomic and phylotranscriptomic trees of 22 mammals and 15 plants that have both sequenced nuclear genomes and publicly available RNA sequencing data from multiple tissues. We found that phylotranscriptomic analysis can be sensitive to orthologous gene identification. When a rigorous method for identifying orthologs is employed, phylogenomic and phylotranscriptomic trees are virtually identical to each other, regardless of the tissue of origin of the transcriptomes and whether the same tissue is used across species. These findings validate phylotranscriptomics, brighten its prospect, and illustrate the criticality of reliable ortholog detection in such practices.

Biochemical activity of magnesium ions on human osteoblast migration.

Magnesium is well known as a biodegradable biomaterial that has been reported to promote bone remodeling in several studies; however, the underlying biological mechanism remains unclear. In the present study, the role of magnesium ions in the migration of U-2 OS cells, which are osteoblast-like cell lines, was investigated. Magnesium treatment did not significantly alter the global transcriptome of U-2 OS cells, but increased the protein expression level of SNAI2, an epithelial-mesenchymal transition (EMT) marker. In addition, it was confirmed that the junctional site localization of Zona-occludens 1 (ZO-1), a representative tight junction protein, was destroyed by magnesium treatment; furthermore, it was determined that cytoplasmic localization increased, and alkaline phosphatase (ALP) activity increased. The obtained results on the mechanism by which magnesium is involved in osteoblast migration, which is important for fracture healing, will contribute to the understanding of the bone-formation process in patients with osteoporosis and musculoskeletal injury.

Identification of Differentially Expressed Genes Associated with Extracellular Matrix Degradation and Inflammatory Regulation in Calcific Tendinopathy Using RNA Sequencing.

Calcific tendinopathy (CT), developed due to calcium hydroxyapatite deposition in the rotator cuff tendon, mostly affects women in their 40 s and 50 s and causes severe shoulder pain. However, the molecular basis of its pathogenesis and appropriate treatment methods are largely unknown. In this study, we identified 202 differentially expressed genes (DEGs) between calcific and adjacent normal tendon tissues of rotator cuff using RNA sequencing-based transcriptome analysis. The DEGs were highly enriched in extracellular matrix (ECM) degradation and inflammation-related processes. Further, matrix metalloproteinase 9 (MMP9) and matrix metalloproteinase 13 (MMP13), two of the enzymes associated with ECM degradation, were found to be highly upregulated 25.85- and 19.40-fold, respectively, in the calcific tendon tissues compared to the adjacent normal tendon tissues. Histopathological analyses indicated collagen degradation and macrophage infiltration at the sites of calcific deposit in the rotator cuff tendon. Our study acts as a foundation that may help in better understanding of the pathogenesis associated with CT, and thus in better management of the disease.

COVID-19 Patients Upregulate Toll-like Receptor 4-mediated Inflammatory Signaling That Mimics Bacterial Sepsis.

BACKGROUND: Observational studies of the ongoing coronavirus disease 2019 (COVID-19) outbreak suggest that a 'cytokine storm' is involved in the pathogenesis of severe illness. However, the molecular mechanisms underlying the altered pathological inflammation in COVID-19 are largely unknown. We report here that toll-like receptor (TLR) 4-mediated inflammatory signaling molecules are upregulated in peripheral blood mononuclear cells (PBMCs) from COVID-19 patients, compared with healthy controls (HC). METHODS: A total of 48 subjects including 28 COVID-19 patients (8 severe/critical vs. 20 mild/moderate cases) admitted to Chungnam National University Hospital, and age/sex-matched 20 HC were enrolled in this study. PBMCs from the subjects were processed for nCounter Human Immunology gene expression assay to analyze the immune related transcriptome profiles. Recombinant proteins of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) were used to stimulate the PBMCs and monocyte-derived macrophages, and real-time polymerase chain reaction was performed to quantify the mRNA expressions of the pro-inflammatory cytokines/chemokines. RESULTS: Among the most highly increased inflammatory mediators in severe/critically ill patients, S100A9, an alarmin and TLR4 ligand, was found as a noteworthy biomarker, because it inversely correlated with the serum albumin levels. We also observed that recombinant S2 and nucleocapsid proteins of SARS-CoV-2 significantly increased pro-inflammatory cytokines/chemokines and S100A9 in human primary PBMCs. CONCLUSION: These data support a link between TLR4 signaling and pathological inflammation during COVID-19 and contribute to develop therapeutic approaches through targeting TLR4-mediated inflammation.

Conventionally used reference genes are not outstanding for normalization of gene expression in human cancer research.

BACKGROUND: The selection of reference genes is essential for quantifying gene expression. Theoretically they should be expressed stably and not regulated by experimental or pathological conditions. However, identification and validation of reference genes for human cancer research are still being regarded as a critical point, because cancerous tissues often represent genetic instability and heterogeneity. Recent pan-cancer studies have demonstrated the importance of the appropriate selection of reference genes for use as internal controls for the normalization of gene expression; however, no stably expressed, consensus reference genes valid for a range of different human cancers have yet been identified. RESULTS: In the present study, we used large-scale cancer gene expression datasets from The Cancer Genome Atlas (TCGA) database, which contains 10,028 (9,364 cancerous and 664 normal) samples from 32 different cancer types, to confirm that the expression of the most commonly used reference genes is not consistent across a range of cancer types. Furthermore, we identified 38 novel candidate reference genes for the normalization of gene expression, independent of cancer type. These genes were found to be highly expressed and highly connected to relevant gene networks, and to be enriched in transcription-translation regulation processes. The expression stability of the newly identified reference genes across 29 cancerous and matched normal tissues were validated via quantitative reverse transcription PCR (RT-qPCR). CONCLUSIONS: We reveal that most commonly used reference genes in current cancer studies cannot be appropriate to serve as representative control genes for quantifying cancer-related gene expression levels, and propose in this study three potential reference genes (HNRNPL, PCBP1, and RER1) to be the most stably expressed across various cancerous and normal human tissues.

Population Genomic Analysis Reveals Contrasting Demographic Changes of Two Closely Related Dolphin Species in the Last Glacial.

Population genomic data can be used to infer historical effective population sizes (Ne), which help study the impact of past climate changes on biodiversity. Previous genome sequencing of one individual of the common bottlenose dolphin Tursiops truncatus revealed an unusual, sharp rise in Ne during the last glacial, raising questions about the reliability, generality, underlying cause, and biological implication of this finding. Here we first verify this result by additional sampling of T. truncatus. We then sequence and analyze the genomes of its close relative, the Indo-Pacific bottlenose dolphin T. aduncus. The two species exhibit contrasting demographic changes in the last glacial, likely through actual changes in population size and/or alterations in the level of gene flow among populations. Our findings suggest that even closely related species can have drastically different responses to climatic changes, making predicting the fate of individual species in the ongoing global warming a serious challenge.

Identification and sequencing of the gene encoding DNA methyltransferase 3 (DNMT3) from sea cucumber, Apostichopus japonicus.

The sea cucumber Apostichopus japonicus is well known as a traditional tonic food and as a commercially important cultured aquatic species. This species produces saponins, and has a high potential to cope with environmental stress, such as aestivation, organ regeneration, and wound healing. Recently, several studies have shown that cellular reprogramming and the physiological responses of the sea cucumber to environmental changes, including aestivation, are potentially mediated by epigenetic DNA methylation. The DNA methyltransferase (DNMT)1 and DNMT3 genes are independent participants in the maintenance and de novo methylation of specific sequences. Sea urchin (Strongylocentrotus purpuratus) and starfish (Asterina pectinifera), which belong to the same phylum as A. japonicus, have both DNMT1 and DNMT3 genes. However, it was previously reported that DNMT1 is present, but DNMT3 is absent, in A. japonicus. In the present study, we sequenced the full-length cDNA of the A. japonicus DNMT3 gene. The newly sequenced DNMT3 gene comprises three major conserved domains (Pro-Trp-Trp-Pro (PWWP), plant homeodomain (PHD), and S-adenosylmethionine-dependent methyltransferase (AdoMet-MTase)), indicating that the DNMT3 possibly has de novo DNA methylation catalytic activity. Gene structure and phylogenetic analysis showed that sea cucumber DNMT3 is evolutionarily conserved in the Echinodermata. Next, we demonstrated the conservation of DNMT3 gene expression in sea cucumber and starfish belong to same phylum, echinoderm. Using reverse transcription-polymerase chain reaction, sea cucumber DNMT3 mRNA was detected in testis tissue, but not in other tissues tested, including the respiratory tree, muscle, tentacle, intestine, and ovary. This is inconsistent with previous reports, which showed the expression of DNMT3 in ovary, but not in testis of the starfish A. pectinifera, indicating the tissue- and species-specific expression of DNMT3 gene. Although further studies are needed to clarify the epigenetic regulatory mechanisms of DNMT3 and its application to the aquaculture industry, our findings may provide insights into the sea cucumber biology.

In silico prediction of potential chemical reactions mediated by human enzymes.

BACKGROUND: Administered drugs are often converted into an ineffective or activated form by enzymes in our body. Conventional in silico prediction approaches focused on therapeutically important enzymes such as CYP450. However, there are more than thousands of different cellular enzymes that potentially convert administered drug into other forms. RESULT: We developed an in silico model to predict which of human enzymes including metabolic enzymes as well as CYP450 family can catalyze a given chemical compound. The prediction is based on the chemical and physical similarity between known enzyme substrates and a query chemical compound. Our in silico model was developed using multiple linear regression and the model showed high performance (AUC = 0.896) despite of the large number of enzymes. When evaluated on a test dataset, it also showed significantly high performance (AUC = 0.746). Interestingly, evaluation with literature data showed that our model can be used to predict not only enzymatic reactions but also drug conversion and enzyme inhibition. CONCLUSION: Our model was able to predict enzymatic reactions of a query molecule with a high accuracy. This may foster to discover new metabolic routes and to accelerate the computational development of drug candidates by enabling the prediction of the potential conversion of administered drugs into active or inactive forms.

Intra and Interspecific Variations of Gene Expression Levels in Yeast Are Largely Neutral: (Nei Lecture, SMBE 2016, Gold Coast).

It is commonly, although not universally, accepted that most intra and interspecific genome sequence variations are more or less neutral, whereas a large fraction of organism-level phenotypic variations are adaptive. Gene expression levels are molecular phenotypes that bridge the gap between genotypes and corresponding organism-level phenotypes. Yet, it is unknown whether natural variations in gene expression levels are mostly neutral or adaptive. Here we address this fundamental question by genome-wide profiling and comparison of gene expression levels in nine yeast strains belonging to three closely related Saccharomyces species and originating from five different ecological environments. We find that the transcriptome-based clustering of the nine strains approximates the genome sequence-based phylogeny irrespective of their ecological environments. Remarkably, only ∼0.5% of genes exhibit similar expression levels among strains from a common ecological environment, no greater than that among strains with comparable phylogenetic relationships but different environments. These and other observations strongly suggest that most intra and interspecific variations in yeast gene expression levels result from the accumulation of random mutations rather than environmental adaptations. This finding has profound implications for understanding the driving force of gene expression evolution, genetic basis of phenotypic adaptation, and general role of stochasticity in evolution.

Histone and RNA-binding protein interaction creates crosstalk network for regulation of alternative splicing.

Alternative splicing is an essential process in eukaryotes, as it increases the complexity of gene expression by generating multiple proteins from a single pre-mRNA. However, information on the regulatory mechanisms for alternative splicing is lacking, because splicing occurs over a short period via the transient interactions of proteins within functional complexes of the spliceosome. Here, we investigated in detail the molecular mechanisms connecting alternative splicing with epigenetic mechanisms. We identified interactions between histone proteins and splicing factors such as Rbfox2, Rbfox3, and splicing factor proline and glutamine rich protein (SFPQ) by in vivo crosslinking and immunoprecipitation. Furthermore, we confirmed that splicing factors were bound to specific modified residues of histone proteins. Additionally, changes in histone methylation due to histone methyltransferase inhibitor treatment notably affected alternative splicing in selected genes. Therefore, we suggested that there may be crosstalk mechanisms connecting histone modifications and RNA-binding proteins that increase the local concentration of RNA-binding proteins in alternative exon loci of nucleosomes by binding specific modified histone proteins, leading to alternative splicing. This crosstalk mechanism may play a major role in epigenetic processes such as histone modification and the regulation of alternative splicing.

Functional Diversification of Motor Neuron-specific Isl1 Enhancers during Evolution.

Functional diversification of motor neurons has occurred in order to selectively control the movements of different body parts including head, trunk and limbs. Here we report that transcription of Isl1, a major gene necessary for motor neuron identity, is controlled by two enhancers, CREST1 (E1) and CREST2 (E2) that allow selective gene expression of Isl1 in motor neurons. Introduction of GFP reporters into the chick neural tube revealed that E1 is active in hindbrain motor neurons and spinal cord motor neurons, whereas E2 is active in the lateral motor column (LMC) of the spinal cord, which controls the limb muscles. Genome-wide ChIP-Seq analysis combined with reporter assays showed that Phox2 and the Isl1-Lhx3 complex bind to E1 and drive hindbrain and spinal cord-specific expression of Isl1, respectively. Interestingly, Lhx3 alone was sufficient to activate E1, and this may contribute to the initiation of Isl1 expression when progenitors have just developed into motor neurons. E2 was induced by onecut 1 (OC-1) factor that permits Isl1 expression in LMCm neurons. Interestingly, the core region of E1 has been conserved in evolution, even in the lamprey, a jawless vertebrate with primitive motor neurons. All E1 sequences from lamprey to mouse responded equally well to Phox2a and the Isl1-Lhx3 complex. Conversely, E2, the enhancer for limb-innervating motor neurons, was only found in tetrapod animals. This suggests that evolutionarily-conserved enhancers permit the diversification of motor neurons.