Hoxa2 and Hoxb2 control dorsoventral patterns of neuronal development in the rostral hindbrain

Higher order interaction analysis quantifies coordination in the epigenome revealing novel biological relationships in Kabuki syndrome

Complex direct and indirect relationships between multiple variables, termed higher order interactions (HOIs), are characteristics of all natural systems. Traditional differential and network analyses fail to account for the omic datasets richness and miss HOIs. We investigated peripheral blood DNA methylation data from Kabuki syndrome type 1 (KS1) and control individuals, identified 2,002 differentially methylated points (DMPs), and inferred 17 differentially methylated regions, which represent only 189 DMPs. We applied hypergraph models to measure HOIs on all the CpGs and revealed differences in the coordination of DMPs with lower entropy and higher coordination of the peripheral epigenome in KS1 implying reduced network complexity. Hypergraphs also capture epigenomic trans-relationships, and identify biologically relevant pathways that escape the standard analyses. These findings construct the basis of a suitable model for the analysis of organization in the epigenome in rare diseases, which can be applied to investigate mechanism in big data.

Single-cell transcriptomics reveals the molecular basis of human iPS cell differentiation into ectodermal ocular lineages

The generation of a self-formed, ectodermal, autonomous multi-zone (SEAM) from human induced pluripotent stem cells (hiPSCs) offers a unique perspective to study the dynamics of ocular cell differentiation over time. Here, by utilising single-cell transcriptomics, we have (i) identified, (ii) molecularly characterised and (iii) ascertained the developmental trajectories of ectodermally-derived ocular cell populations which emerge within SEAMs as they form. Our analysis reveals interdependency between tissues of the early eye and delineates the sequential formation and maturation of distinct cell types over a 12-week period. We demonstrate a progression from pluripotency through to tissue specification and differentiation which encompasses both surface ectodermal and neuroectodermal ocular lineages and the generation of iPSC-derived components of the developing cornea, conjunctiva, lens, and retina. Our findings not only advance the understanding of ocular development in a stem cell-based system of human origin, but also establish a robust methodological paradigm for exploring cellular and molecular dynamics during SEAM formation at single-cell resolution and highlight the potential of hiPSC-derived systems as powerful platforms for modelling human eye development and disease.

Dissecting the genetic and causal relationship between sleep-related traits and common brain disorders

BACKGROUND

There is a profound connection between abnormal sleep patterns and brain disorders, suggesting a shared influential association. However, the shared genetic basis and potential causal relationships between sleep-related traits and brain disorders are yet to be fully elucidated.

METHODS

Utilizing linkage disequilibrium score regression (LDSC) and bidirectional two-sample univariable Mendelian Randomization (UVMR) analyses with large-scale GWAS datasets, we investigated the genetic correlations and causal associations across six sleep traits and 24 prevalent brain disorders. Additionally, a multivariable Mendelian Randomization (MVMR) analysis evaluated the cumulative effects of various sleep traits on each brain disorder, complemented by genetic loci characterization to pinpoint pertinent genes and pathways.

RESULTS

LDSC analysis identified significant genetic correlations in 66 out of 144 (45.8 %) pairs between sleep-related traits and brain disorders, with the most pronounced correlations observed in psychiatric disorders (66 %, 48/72). UVMR analysis identified 29 causal relationships (FDR<0.05) between sleep traits and brain disorders, with 19 associations newly discovered according to our knowledge. Notably, major depression, attention-deficit/hyperactivity disorder, bipolar disorder, cannabis use disorder, and anorexia nervosa showed bidirectional causal relations with sleep traits, especially insomnia's marked influence on major depression (IVW beta 0.468, FDR = 5.24E-09). MVMR analysis revealed a nuanced interplay among various sleep traits and their impact on brain disorders. Genetic loci characterization underscored potential genes, such as HOXB2, while further enrichment analyses illuminated the importance of synaptic processes in these relationships.

CONCLUSIONS

This study provides compelling evidence for the causal relationships and shared genetic backgrounds between common sleep-related traits and brain disorders.

Molecular mechanisms reconstruction from single-cell multi-omics data with HuMMuS

MOTIVATION

The molecular identity of a cell results from a complex interplay between heterogeneous molecular layers. Recent advances in single-cell sequencing technologies have opened the possibility to measure such molecular layers of regulation.

RESULTS

Here, we present HuMMuS, a new method for inferring regulatory mechanisms from single-cell multi-omics data. Differently from the state-of-the-art, HuMMuS captures cooperation between biological macromolecules and can easily include additional layers of molecular regulation. We benchmarked HuMMuS with respect to the state-of-the-art on both paired and unpaired multi-omics datasets. Our results proved the improvements provided by HuMMuS in terms of transcription factor (TF) targets, TF binding motifs and regulatory regions prediction. Finally, once applied to snmC-seq, scATAC-seq and scRNA-seq data from mouse brain cortex, HuMMuS enabled to accurately cluster scRNA profiles and to identify potential driver TFs.

AVAILABILITY AND IMPLEMENTATION

HuMMuS is available at https://github.com/cantinilab/HuMMuS.

Transcriptomic profiling across human serotonin neuron differentiation via the FEV reporter system

BACKGROUND

The detailed transcriptomic profiles during human serotonin neuron (SN) differentiation remain elusive. The establishment of a reporter system based on SN terminal selector holds promise to produce highly-purified cells with an early serotonergic fate and help elucidate the molecular events during human SN development process.

METHODS

A fifth Ewing variant (FEV)-EGFP reporter system was established by CRISPR/Cas9 technology to indicate SN since postmitotic stage. FACS was performed to purify SN from the heterogeneous cell populations. RNA-sequencing analysis was performed for cells at four key stages of differentiation (pluripotent stem cells, serotonergic neural progenitors, purified postmitotic SN and purifed mature SN) to explore the transcriptomic dynamics during SN differentiation.

RESULTS

We found that human serotonergic fate specification may commence as early as day 21 of differentiation from human pluripotent stem cells. Furthermore, the transcriptional factors ZIC1, HOXA2 and MSX2 were identified as the hub genes responsible for orchestrating serotonergic fate determination.

CONCLUSIONS

For the first time, we exposed the developmental transcriptomic profiles of human SN via FEV reporter system, which will further our understanding for the development process of human SN.

Developmental patterns of extracellular matrix molecules in the embryonic and postnatal mouse hindbrain

Normal brain development requires continuous communication between developing neurons and their environment filled by a complex network referred to as extracellular matrix (ECM). The ECM is divided into distinct families of molecules including hyaluronic acid, proteoglycans, glycoproteins such as tenascins, and link proteins. In this study, we characterize the temporal and spatial distribution of the extracellular matrix molecules in the embryonic and postnatal mouse hindbrain by using antibodies and lectin histochemistry. In the embryo, hyaluronan and neurocan were found in high amounts until the time of birth whereas versican and tenascin-R were detected in lower intensities during the whole embryonic period. After birth, both hyaluronic acid and neurocan still produced intense staining in almost all areas of the hindbrain, while tenascin-R labeling showed a continuous increase during postnatal development. The reaction with WFA and aggrecan was revealed first 4th postnatal day (P4) with low staining intensities, while HAPLN was detected two weeks after birth (P14). The perineuronal net appeared first around the facial and vestibular neurons at P4 with hyaluronic acid cytochemistry. One week after birth aggrecan, neurocan, tenascin-R, and WFA were also accumulated around the neurons located in several hindbrain nuclei, but HAPLN1 was detected on the second postnatal week. Our results provide further evidence that many extracellular macromolecules that will be incorporated into the perineuronal net are already expressed at embryonic and early postnatal stages of development to control differentiation, migration, and synaptogenesis of neurons. In late postnatal period, the experience-driven neuronal activity induces formation of perineuronal net to stabilize synaptic connections.

Transcriptional dynamics of delaminating neuroblasts in the mouse otic vesicle

An abundance of research has recently highlighted the susceptibility of cochleovestibular ganglion (CVG) neurons to noise damage and aging in the adult cochlea, resulting in hearing deficits. Furthering our understanding of the transcriptional cascades that contribute to CVG development may provide insight into how these cells can be regenerated to treat inner ear dysfunction. Here we perform a high-depth single-cell RNA sequencing analysis of the E10.5 otic vesicle and its surrounding tissues, including CVG precursor neuroblasts and emerging CVG neurons. Clustering and trajectory analysis of otic-lineage cells reveals otic markers and the changes in gene expression that occur from neuroblast delamination toward the development of the CVG. This dataset provides a valuable resource for further identifying the mechanisms associated with CVG development from neurosensory competent cells within the otic vesicle.

Crotonylation of GAPDH regulates human embryonic stem cell endodermal lineage differentiation and metabolic switch

BACKGROUND

Post-translational modifications of proteins are crucial to the regulation of their activity and function. As a newly discovered acylation modification, crotonylation of non-histone proteins remains largely unexplored, particularly in human embryonic stem cells (hESCs).

METHODS

We investigated the role of crotonylation in hESC differentiation by introduce crotonate into the culture medium of GFP tagged LTR7 primed H9 cell and extended pluripotent stem cell lines. RNA-seq assay was used to determine the hESC transcriptional features. Through morphological changes, qPCR of pluripotent and germ layer-specific gene markers and flow cytometry analysis, we determined that the induced crotonylation resulted in hESC differentiating into the endodermal lineage. We performed targeted metabolomic analysis and seahorse metabolic measurement to investigate the metabolism features after crotonate induction. Then high-resolution tandem mass spectrometry (LC-MS/MS) revealed the target proteins in hESCs. In addition, the role of crotonylated glycolytic enzymes (GAPDH and ENOA) was evaluated by in vitro crotonylation and enzymatic activity assays. Finally, we used knocked-down hESCs by shRNA, wild GAPDH and GAPDH mutants to explore potential role of GAPDH crotonylation in regulating human embryonic stem cell differentiation and metabolic switch.

RESULT

We found that induced crotonylation in hESCs resulted in hESCs of different pluripotency states differentiating into the endodermal lineage. Increased protein crotonylation in hESCs was accompanied by transcriptomic shifts and decreased glycolysis. Large-scale crotonylation profiling of non-histone proteins revealed that metabolic enzymes were major targets of inducible crotonylation in hESCs. We further discovered GAPDH as a key glycolytic enzyme regulated by crotonylation during endodermal differentiation from hESCs.

CONCLUSIONS

Crotonylation of GAPDH decreased its enzymatic activity thereby leading to reduced glycolysis during endodermal differentiation from hESCs.

Cannabichromene Induces Neuronal Differentiation in NSC-34 Cells: Insights from Transcriptomic Analysis

Phytocannabinoids, with their variety of beneficial effects, represent a valid group of substances that could be employed as neurogenesis-enhancers or neuronal differentiation inducers. We focused our attention on the neuronal-related potential of cannabichromene (CBC) when administered to undifferentiated NSC-34 for 24 h. Transcriptomic analysis showed an upregulation of several neuronal markers, such as Neurod1 and Tubb3, as well as indicators of neuronal differentiation process progression, such as Pax6. An in-depth investigation of the processes involved in neuronal differentiation indicates positive cytoskeleton remodeling by upregulation of Cfl2 and Tubg1, and active differentiation-targeted transcriptional program, suggested by Phox2b and Hes1. After 48 h of treatment, the markers previously examined in the transcriptomic analysis are still overexpressed, like Ache and Hes1, indicating that the differentiation process is still in progress. The lack of GFAP protein suggests that no astroglial differentiation is taking place, and it is reasonable to indicate the neuronal one as the ongoing one. These results indicate CBC as a potential neuronal differentiation inducer for NSC-34 cells.

Hox genes in development and beyond

Hox genes encode evolutionarily conserved transcription factors that are essential for the proper development of bilaterian organisms. Hox genes are unique because they are spatially and temporally regulated during development in a manner that is dictated by their tightly linked genomic organization. Although their genetic function during embryonic development has been interrogated, less is known about how these transcription factors regulate downstream genes to direct morphogenetic events. Moreover, the continued expression and function of Hox genes at postnatal and adult stages highlights crucial roles for these genes throughout the life of an organism. Here, we provide an overview of Hox genes, highlighting their evolutionary history, their unique genomic organization and how this impacts the regulation of their expression, what is known about their protein structure, and their deployment in development and beyond.

Modular derivation of diverse, regionally discrete human posterior CNS neurons enables discovery of transcriptomic patterns

Our inability to derive the neuronal diversity that comprises the posterior central nervous system (pCNS) using human pluripotent stem cells (hPSCs) poses an impediment to understanding human neurodevelopment and disease in the hindbrain and spinal cord. Here, we establish a modular, monolayer differentiation paradigm that recapitulates both rostrocaudal (R/C) and dorsoventral (D/V) patterning, enabling derivation of diverse pCNS neurons with discrete regional specificity. First, neuromesodermal progenitors (NMPs) with discrete HOX profiles are converted to pCNS progenitors (pCNSPs). Then, by tuning D/V signaling, pCNSPs are directed to locomotor or somatosensory neurons. Expansive single-cell RNA-sequencing (scRNA-seq) analysis coupled with a novel computational pipeline allowed us to detect hundreds of transcriptional markers within region-specific phenotypes, enabling discovery of gene expression patterns across R/C and D/V developmental axes. These findings highlight the potential of these resources to advance a mechanistic understanding of pCNS development, enhance in vitro models, and inform therapeutic strategies.

Dynamic regulatory module networks for inference of cell type-specific transcriptional networks

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Sex-Stratified Single-Cell RNA-Seq Analysis Identifies Sex-Specific and Cell Type-Specific Transcriptional Responses in Alzheimer's Disease Across Two Brain Regions

Alzheimer’s disease (AD) is a pervasive neurodegenerative disorder that disproportionately affects women. Since neural anatomy and disease pathophysiology differ by sex, investigating sex-specific mechanisms in AD pathophysiology can inform new therapeutic approaches for both sexes. Previous bulk human brain RNA sequencing studies have revealed sex differences in dysregulated molecular pathways related to energy production, neuronal function, and immune response; however, the sex differences in disease mechanisms are yet to be examined comprehensively on a single-cell level. We leveraged nearly 74,000 cells from human prefrontal and entorhinal cortex samples from the first two publicly available single-cell RNA sequencing AD datasets to perform a case versus control sex-stratified differential gene expression analysis and pathway network enrichment in a cell type-specific manner for each brain region. Our examination at the single-cell level revealed sex differences in AD prominently in glial cells of the prefrontal cortex. In the entorhinal cortex, we observed the same genes and networks to be perturbed in opposing directions between sexes in AD relative to healthy state. Our findings contribute to growing evidence of sex differences in AD-related transcriptomic changes, which can fuel the development of therapies that may prove more effective at reversing AD pathophysiology.

Association of peripheral blood DNA methylation level with Alzheimer's disease progression

Background

Identifying biomarkers associated with Alzheimer’s disease (AD) progression may enable patient enrichment and improve clinical trial designs. Epigenome-wide association studies have revealed correlations between DNA methylation at cytosine-phosphate-guanine (CpG) sites and AD pathology and diagnosis. Here, we report relationships between peripheral blood DNA methylation profiles measured using Infinium® MethylationEPIC BeadChip and AD progression in participants from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort.

Results

The rate of cognitive decline from initial DNA sampling visit to subsequent visits was estimated by the slopes of the modified Preclinical Alzheimer Cognitive Composite (mPACC; mPACC_digit and mPACC_trailsB) and Clinical Dementia Rating Scale Sum of Boxes (CDR-SB) plots using robust linear regression in cognitively normal (CN) participants and patients with mild cognitive impairment (MCI), respectively. In addition, diagnosis conversion status was assessed using a dichotomized endpoint. Two CpG sites were significantly associated with the slope of mPACC in CN participants (P < 5.79 × 10⁻⁸ [Bonferroni correction threshold]); cg00386386 was associated with the slope of mPACC_digit, and cg09422696 annotated to RP11-661A12.5 was associated with the slope of CDR-SB. No significant CpG sites associated with diagnosis conversion status were identified. Genes involved in cognition and learning were enriched. A total of 19, 13, and 5 differentially methylated regions (DMRs) associated with the slopes of mPACC_trailsB, mPACC_digit, and CDR-SB, respectively, were identified by both comb-p and DMRcate algorithms; these included DMRs annotated to HOXA4. Furthermore, 5 and 19 DMRs were associated with conversion status in CN and MCI participants, respectively. The most significant DMR was annotated to the AD-associated gene PM20D1 (chr1: 205,818,956 to 205,820,014 [13 probes], Sidak-corrected P = 7.74 × 10⁻²⁴), which was associated with both the slope of CDR-SB and the MCI conversion status.

Conclusion

Candidate CpG sites and regions in peripheral blood were identified as associated with the rate of cognitive decline in participants in the ADNI cohort. While we did not identify a single CpG site with sufficient clinical utility to be used by itself due to the observed effect size, a biosignature composed of DNA methylation changes may have utility as a prognostic biomarker for AD progression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01179-2.

In Silico Analysis to Explore Lineage-Independent and -Dependent Transcriptional Programs Associated with the Process of Endothelial and Neural Differentiation of Human Induced Pluripotent Stem Cells

Despite a major interest in understanding how the endothelial cell phenotype is established, the underlying molecular basis of this process is not yet fully understood. We have previously reported the generation of induced pluripotent stem cells (iPS) from human umbilical vein endothelial cells and differentiation of the resulting HiPS back to endothelial cells (Ec-Diff), as well as neural (Nn-Diff) cell lineage that contained both neurons and astrocytes. Furthermore, the identities of these cell lineages were established by gene array analysis. Here, we explored the same arrays to gain insight into the gene alteration processes that accompany the establishment of endothelial vs. non-endothelial neural cell phenotypes. We compared the expression of genes that code for transcription factors and epigenetic regulators when HiPS is differentiated into these endothelial and non-endothelial lineages. Our in silico analyses have identified cohorts of genes that are similarly up- or downregulated in both lineages, as well as those that exhibit lineage-specific alterations. Based on these results, we propose that genes that are similarly altered in both lineages participate in priming the stem cell for differentiation in a lineage-independent manner, whereas those that are differentially altered in endothelial compared to neural cells participate in a lineage-specific differentiation process. Specific GATA family members and their cofactors and epigenetic regulators (DNMT3B, PRDM14, HELLS) with a major role in regulating DNA methylation were among participants in priming HiPS for lineage-independent differentiation. In addition, we identified distinct cohorts of transcription factors and epigenetic regulators whose alterations correlated specifically with the establishment of endothelial vs. non-endothelial neural lineages.

Targeted chromatin conformation analysis identifies novel distal neural enhancers of ZEB2 in pluripotent stem cell differentiation

The transcription factor zinc finger E-box binding protein 2 (ZEB2) controls embryonic and adult cell fate decisions and cellular maturation in many stem/progenitor cell types. Defects in these processes in specific cell types underlie several aspects of Mowat–Wilson syndrome (MOWS), which is caused by ZEB2 haplo-insufficiency. Human ZEB2, like mouse Zeb2, is located on chromosome 2 downstream of a ±3.5 Mb-long gene-desert, lacking any protein-coding gene. Using temporal targeted chromatin capture (T2C), we show major chromatin structural changes based on mapping in-cis proximities between the ZEB2 promoter and this gene desert during neural differentiation of human-induced pluripotent stem cells, including at early neuroprogenitor cell (NPC)/rosette state, where ZEB2 mRNA levels increase significantly. Combining T2C with histone-3 acetylation mapping, we identified three novel candidate enhancers about 500 kb upstream of the ZEB2 transcription start site. Functional luciferase-based assays in heterologous cells and NPCs reveal co-operation between these three enhancers. This study is the first to document in-cis Regulatory Elements located in ZEB2’s gene desert. The results further show the usability of T2C for future studies of ZEB2 REs in differentiation and maturation of multiple cell types and the molecular characterization of newly identified MOWS patients that lack mutations in ZEB2 protein-coding exons.

Anterior Hox Genes and the Process of Cephalization

During evolution, bilateral animals have experienced a progressive process of cephalization with the anterior concentration of nervous tissue, sensory organs and the appearance of dedicated feeding structures surrounding the mouth. Cephalization has been achieved by the specialization of the unsegmented anterior end of the body (the acron) and the sequential recruitment to the head of adjacent anterior segments. Here we review the key developmental contribution of Hox1-5 genes to the formation of cephalic structures in vertebrates and arthropods and discuss how this evolved. The appearance of Hox cephalic genes preceded the evolution of a highly specialized head in both groups, indicating that Hox gene involvement in the control of cephalic structures was acquired independently during the evolution of vertebrates and invertebrates to regulate the genes required for head innovation.

Segmentation and patterning of the vertebrate hindbrain

During early development, the hindbrain is sub-divided into rhombomeres that underlie the organisation of neurons and adjacent craniofacial tissues. A gene regulatory network of signals and transcription factors establish and pattern segments with a distinct anteroposterior identity. Initially, the borders of segmental gene expression are imprecise, but then become sharply defined, and specialised boundary cells form. In this Review, we summarise key aspects of the conserved regulatory cascade that underlies the formation of hindbrain segments. We describe how the pattern is sharpened and stabilised through the dynamic regulation of cell identity, acting in parallel with cell segregation. Finally, we discuss evidence that boundary cells have roles in local patterning, and act as a site of neurogenesis within the hindbrain.

Neural tube patterning: from a minimal model for rostrocaudal patterning toward an integrated 3D model

Rostrocaudal patterning of the neural tube is a defining event in vertebrate brain development. This process is driven by morphogen gradients which specify the fate of neural progenitor cells, leading to the partitioning of the tube. Although this is extensively studied experimentally, an integrated view of the genetic circuitry is lacking. Here, we present a minimal gene regulatory model for rostrocaudal patterning, whose tristable topology was determined in a data-driven way. Using this model, we identified the repression of hindbrain fate as promising strategy for the improvement of current protocols for the generation of dopaminergic neurons. Furthermore, we combined our model with an established minimal model for dorsoventral patterning on a realistic 3D neural tube and found that key features of neural tube patterning could be recapitulated. Doing so, we demonstrate how data and models from different sources can be combined to simulate complex in vivo processes.

Postmitotic Hoxa5 Expression Specifies Pontine Neuron Positional Identity and Input Connectivity of Cortical Afferent Subsets

The mammalian precerebellar pontine nucleus (PN) has a main role in relaying cortical information to the cerebellum. The molecular determinants establishing ordered connectivity patterns between cortical afferents and precerebellar neurons are largely unknown. We show that expression of Hox5 transcription factors is induced in specific subsets of postmitotic PN neurons at migration onset. Hox5 induction is achieved by response to retinoic acid signaling, resulting in Jmjd3-dependent derepression of Polycomb chromatin and 3D conformational changes. Hoxa5 drives neurons to settle posteriorly in the PN, where they are monosynaptically targeted by cortical neuron subsets mainly carrying limb somatosensation. Furthermore, Hoxa5 postmigratory ectopic expression in PN neurons is sufficient to attract cortical somatosensory inputs regardless of position and avoid visual afferents. Transcriptome analysis further suggests that Hoxa5 is involved in circuit formation. Thus, Hoxa5 coordinates postmitotic specification, migration, settling position, and sub-circuit assembly of PN neuron subsets in the cortico-cerebellar pathway.

Estrogen suppresses HOXB2 expression via ERα in breast cancer cells

The expression of HOXB2, a homeobox transcription factor, is altered in a variety of solid tumors. Using an in vivo screen to identify regulators of breast tumor growth in murine mammary fat pads, Boimel and co-workers recently identified HOXB2 as a tumor suppressor. However, the mechanistic underpinnings of its role in breast cancer is not understood. Given the emerging interaction of estrogen-regulated gene expression and altered HOX gene expression network in the pathophysiology of breast cancer, this study addressed the relationship between estrogen signaling and HOXB2 expression. Using a mouse model and human breast cancer cell lines, we show that estrogen suppresses HOXB2 expression. Suppression of HOXB2 by PPT, a known ERα agonist, in MCF-7 and T47D cells indicated the involvement of ERα, which was confirmed by siRNA-mediated ERα knockdown experiments. In-silico analysis of the upstream promoter region revealed the presence of three putative EREs. Chromatin immunoprecipitation experiments showed that upon estrogen binding, ERα engaged with EREs in the 5' upstream region of HOXB2 in MCF-7 and T47D cells. Future investigations should address the implications of estrogen-mediated suppression on the proposed tumor suppressor function of HOXB2.

Gbx2 Is Required for the Migration and Survival of a Subpopulation of Trigeminal Cranial Neural Crest Cells

The development of key structures within the mature vertebrate hindbrain requires the migration of neural crest (NC) cells and motor neurons to their appropriate target sites. Functional analyses in multiple species have revealed a requirement for the transcription factor gastrulation-brain-homeobox 2 (Gbx2) in NC cell migration and positioning of motor neurons in the developing hindbrain. In addition, loss of Gbx2 function studies in mutant mouse embryos, Gbx2^neo, demonstrate a requirement for Gbx2 for the development of NC-derived sensory neurons and axons constituting the mandibular branch of the trigeminal nerve (CNV). Our recent GBX2 target gene identification study identified multiple genes required for the migration and survival of NC cells (e.g., Robo1, Slit3, Nrp1). In this report, we performed loss-of-function analyses using Gbx2^neo mutant embryos, to improve our understanding of the molecular and genetic mechanisms regulated by Gbx2 during anterior hindbrain and CNV development. Analysis of Tbx20 expression in the hindbrain of Gbx2^neo homozygotes revealed a severely truncated rhombomere (r)2. Our data also provide evidence demonstrating a requirement for Gbx2 in the temporal regulation of Krox20 expression in r3. Lastly, we show that Gbx2 is required for the expression of Nrp1 in a subpopulation of trigeminal NC cells, and correct migration and survival of cranial NC cells that populate the trigeminal ganglion. Taken together, these findings provide additional insight into molecular and genetic mechanisms regulated by Gbx2 that underlie NC migration, trigeminal ganglion assembly, and, more broadly, anterior hindbrain development.

Transcriptomic and epigenomic dynamics associated with development of human iPSC-derived GABAergic interneurons

GABAergic interneurons (GINs) are a heterogeneous population of inhibitory neurons that collectively contribute to the maintenance of normal neuronal excitability and network activity. Identification of the genetic regulatory elements and transcription factors that contribute toward GIN function may provide new insight into the pathways underlying proper GIN activity while also indicating potential therapeutic targets for GIN-associated disorders, such as schizophrenia and epilepsy. In this study, we examined the temporal changes in gene expression and chromatin accessibility during GIN development by performing transcriptomic and epigenomic analyses on human induced pluripotent stem cell-derived neurons at 22, 50 and 78 days (D) post-differentiation. We observed 13 221 differentially accessible regions (DARs) of chromatin that associate with temporal changes in gene expression at D78 and D50, relative to D22. We also classified families of transcription factors that are increasingly enriched at DARs during differentiation, indicating regulatory networks that likely drive GIN development. Collectively, these data provide a resource for examining the molecular networks regulating GIN functionality.

Retinoids and developmental neurotoxicity: Utilizing toxicogenomics to enhance adverse outcome pathways and testing strategies

The use of genomic approaches in toxicological studies has greatly increased our ability to define the molecular profiles of environmental chemicals associated with developmental neurotoxicity (DNT). Integration of these approaches with adverse outcome pathways (AOPs), a framework that translates environmental exposures to adverse developmental phenotypes, can potentially inform DNT testing strategies. Here, using retinoic acid (RA) as a case example, we demonstrate that the integration of toxicogenomic profiles into the AOP framework can be used to establish a paradigm for chemical testing. RA is a critical regulatory signaling molecule involved in multiple aspects of mammalian central nervous system (CNS) development, including hindbrain formation/patterning and neuronal differentiation, and imbalances in RA signaling pathways are linked with DNT. While the mechanisms remain unresolved, environmental chemicals can cause DNT by disrupting the RA signaling pathway. First, we reviewed literature evidence of RA and other retinoid exposures and DNT to define a provisional AOP related to imbalances in RA embryonic bioavailability and hindbrain development. Next, by integrating toxicogenomic datasets, we defined a relevant transcriptomic signature associated with RA-induced developmental neurotoxicity (RA-DNT) in human and rodent models that was tested against zebrafish model data, demonstrating potential for integration into an AOP framework. Finally, we demonstrated how these approaches may be systematically utilized to identify chemical hazards by testing the RA-DNT signature against azoles, a proposed class of compounds that alters RA-signaling. The provisional AOP from this study can be expanded in the future to better define DNT biomarkers relevant to RA signaling and toxicity.

Homeobox B2 is a potential prognostic biomarker of glioblastoma

SUMMARY OBJECTIVES HOXB2 is a new prognostic indicator for lung cancer. But it is unclear whether HOXB2 holds an effect in glioblastoma (GBM) progression. The purpose of this article was to probe the influences of HOXB2 on GBM pathogenesis. METHODS HOXB2 expression level and prognostic power in GBM patients were analyzed. Then the mRNA and protein expression levels of HOXB2 in GBM cell lines were tested by qRT-PCR and western blotting. Cell proliferation, invasion, and migration were determined by CCK8 and transwell assay, severally. The protein levels of PI3K/AKT-pathway associated proteins were analyzed by western blotting. RESULTS The results indicated that HOXB2 was distinctly overexpressed in GBM patients and high expression of HOXB2 was related to a poor prognosis. Moreover, the expression of HOXB2 was higher in all GBM cell lines U251, U-87MG, GOS-3 than that in HEB cells (normal control). Meanwhile, decreased expression of p-PI3K and p-AKT were identified after HOXB2 knockdown. CONCLUSIONS These data demonstrated that HOXB2 had a vital role in GBM progression and could serve as a promising target for GBM treatment.

Single-Cell RNA-seq Identifies Cell Diversity in Embryonic Salivary Glands

Branching organs, including the salivary and mammary glands, lung, and kidney, arise as epithelial buds that are morphologically very similar. However, the mesenchyme is known to guide epithelial morphogenesis and to help govern cell fate and eventual organ specificity. We performed single-cell transcriptome analyses of 14,441 cells from embryonic day 12 submandibular and parotid salivary glands to characterize their molecular identities during bud initiation. The mesenchymal cells were considerably more heterogeneous by clustering analysis than the epithelial cells. Nonetheless, distinct clusters were evident among even the epithelial cells, where unique molecular markers separated presumptive bud and duct cells. Mesenchymal cells formed separate, well-defined clusters specific to each gland. Neuronal and muscle cells of the 2 glands in particular showed different markers and localization patterns. Several gland-specific genes were characteristic of different rhombomeres. A muscle cluster was prominent in the parotid, which was not myoepithelial or vascular smooth muscle. Instead, the muscle cluster expressed genes that mediate skeletal muscle differentiation and function. Striated muscle was indeed found later in development surrounding the parotid gland. Distinct spatial localization patterns of neuronal and muscle cells in embryonic stages appear to foreshadow later differences in adult organ function. These findings demonstrate that the establishment of transcriptional identities emerges early in development, primarily in the mesenchyme of developing salivary glands. We present the first comprehensive description of molecular signatures that define specific cellular landmarks for the bud initiation stage, when the neural crest-derived ectomesenchyme predominates in the salivary mesenchyme that immediately surrounds the budding epithelium. We also provide the first transcriptome data for the largely understudied embryonic parotid gland as compared with the submandibular gland, focusing on the mesenchymal cell populations.

Hindbrain induction and patterning during early vertebrate development

The hindbrain is a key relay hub of the central nervous system (CNS), linking the bilaterally symmetric half-sides of lower and upper CNS centers via an extensive network of neural pathways. Dedicated neural assemblies within the hindbrain control many physiological processes, including respiration, blood pressure, motor coordination and different sensations. During early development, the hindbrain forms metameric segmented units known as rhombomeres along the antero-posterior (AP) axis of the nervous system. These compartmentalized units are highly conserved during vertebrate evolution and act as the template for adult brainstem structure and function. TALE and HOX homeodomain family transcription factors play a key role in the initial induction of the hindbrain and its specification into rhombomeric cell fate identities along the AP axis. Signaling pathways, such as canonical-Wnt, FGF and retinoic acid, play multiple roles to initially induce the hindbrain and regulate Hox gene-family expression to control rhombomeric identity. Additional transcription factors including Krox20, Kreisler and others act both upstream and downstream to Hox genes, modulating their expression and protein activity. In this review, we will examine the earliest embryonic signaling pathways that induce the hindbrain and subsequent rhombomeric segmentation via Hox and other gene expression. We will examine how these signaling pathways and transcription factors interact to activate downstream targets that organize the segmented AP pattern of the embryonic vertebrate hindbrain.

Hox5 genes direct elastin network formation during alveologenesis by regulating myofibroblast adhesion

Hox5 genes (Hoxa5, Hoxb5, Hoxc5) are exclusively expressed in the lung mesenchyme during embryogenesis, and the most severe phenotypes result from constitutive loss of function of all three genes. Because Hox5 triple null mutants exhibit perinatal lethality, the contribution of this paralogous group to postembryonic lung development is unknown. Intriguingly, expression of all three Hox5 genes peaks during the first 2 weeks after birth, reaching levels far exceeding those measured at embryonic stages, and surviving Hoxa5 single and Hox5 AabbCc compound mutants exhibit defects in the localization of alveolar myofibroblasts. To define the contribution of the entire Hox5 paralogous group to this process, we generated an Hoxa5 conditional allele to use with our existing null alleles for Hoxb5 and Hoxc5 Postnatally, mesenchymal deletion of Hoxa5 in an Hoxb5/Hoxc5 double-mutant background results in severe alveolar simplification. The elastin network required for alveolar formation is dramatically disrupted in Hox5 triple mutants, while the basal lamina, interstitial matrix, and fibronectin are normal. Alveolar myofibroblasts remain Pdgfrα+/SMA+ double positive and present in normal numbers, indicating that the irregular elastin network is not due to fibroblast differentiation defects. Rather, we observe that SMA+ myofibroblasts of Hox5 triple mutants are morphologically abnormal both in vivo and in vitro with highly reduced adherence to fibronectin. This loss of adhesion is a result of loss of the integrin heterodimer Itga5b1 in mutant fibroblasts. Collectively, these data show an important role for Hox5 genes in lung fibroblast adhesion necessary for proper elastin network formation during alveologenesis.

Cellular and Molecular Underpinnings of Neuronal Assembly in the Central Auditory System during Mouse Development

During development, the organization of the auditory system into distinct functional subcircuits depends on the spatially and temporally ordered sequence of neuronal specification, differentiation, migration and connectivity. Regional patterning along the antero-posterior axis and neuronal subtype specification along the dorso-ventral axis intersect to determine proper neuronal fate and assembly of rhombomere-specific auditory subcircuits. By taking advantage of the increasing number of transgenic mouse lines, recent studies have expanded the knowledge of developmental mechanisms involved in the formation and refinement of the auditory system. Here, we summarize several findings dealing with the molecular and cellular mechanisms that underlie the assembly of central auditory subcircuits during mouse development, focusing primarily on the rhombomeric and dorso-ventral origin of auditory nuclei and their associated molecular genetic pathways.

Characterization of chromosomal abnormalities in pregnancy losses reveals critical genes and loci for human early development

Detailed characterization of chromosomal abnormalities, a common cause for congenital abnormalities and pregnancy loss, is critical for elucidating genes for human fetal development. Here, 2,186 product-of-conception samples were tested for copy-number variations (CNVs) at two clinical diagnostic centers using whole-genome sequencing and high-resolution chromosomal microarray analysis. We developed a new gene discovery approach to predict potential developmental genes and identified 275 candidate genes from CNVs detected from both datasets. Based on Mouse Genome Informatics (MGI) and Zebrafish model organism database (ZFIN), 75% of identified genes could lead to developmental defects when mutated. Genes involved in embryonic development, gene transcription, and regulation of biological processes were significantly enriched. Especially, transcription factors and gene families sharing specific protein domains predominated, which included known developmental genes such as HOX, NKX homeodomain genes, and helix-loop-helix containing HAND2, NEUROG2, and NEUROD1 as well as potential novel developmental genes. We observed that developmental genes were denser in certain chromosomal regions, enabling identification of 31 potential genomic loci with clustered genes associated with development.

The vertebrate Hox gene regulatory network for hindbrain segmentation: Evolution and diversification: Coupling of a Hox gene regulatory network to hindbrain segmentation is an ancient trait originating at the base of vertebrates

Hindbrain development is orchestrated by a vertebrate gene regulatory network that generates segmental patterning along the anterior-posterior axis via Hox genes. Here, we review analyses of vertebrate and invertebrate chordate models that inform upon the evolutionary origin and diversification of this network. Evidence from the sea lamprey reveals that the hindbrain regulatory network generates rhombomeric compartments with segmental Hox expression and an underlying Hox code. We infer that this basal feature was present in ancestral vertebrates and, as an evolutionarily constrained developmental state, is fundamentally important for patterning of the vertebrate hindbrain across diverse lineages. Despite the common ground plan, vertebrates exhibit neuroanatomical diversity in lineage-specific patterns, with different vertebrates revealing variations of Hox expression in the hindbrain that could underlie this diversification. Invertebrate chordates lack hindbrain segmentation but exhibit some conserved aspects of this network, with retinoic acid signaling playing a role in establishing nested domains of Hox expression.

Functional and Comparative Genomics of Hoxa2 Gene cis-Regulatory Elements: Evidence for Evolutionary Modification of Ancestral Core Element Activity

Hoxa2 is an evolutionarily conserved developmental regulatory gene that functions to specify rhombomere (r) and pharyngeal arch (PA) identities throughout the Osteichthyes. Japanese medaka (Oryzias latipes) hoxa2a, like orthologous Hoxa2 genes from other osteichthyans, is expressed during embryogenesis in r2–7 and PA2-7, whereas the paralogous medaka pseudogene, ψhoxa2b, is expressed in noncanonical Hoxa2 domains, including the pectoral fin buds. To understand the evolution of cis-regulatory element (CRE) control of gene expression, we conducted eGFP reporter gene expression studies with extensive functional mapping of several conserved CREs upstream of medaka hoxa2a and ψhoxa2b in transient and stable-line transgenic medaka embryos. The CREs tested were previously shown to contribute to directing mouse Hoxa2 gene expression in r3, r5, and PA2-4. Our results reveal the presence of sequence elements embedded in the medaka hoxa2a and ψhoxa2b upstream enhancer regions (UERs) that mediate expression in r4 and the PAs (hoxa2a r4/CNCC element) or in r3–7 and the PAs ψhoxa2b r3–7/CNCC element), respectively. Further, these elements were shown to be highly conserved among osteichthyans, which suggests that the r4 specifying element embedded in the UER of Hoxa2 is a deeply rooted rhombomere specifying element and the activity of this element has been modified by the evolution of flanking sequences that redirect its activity to alternative developmental compartments.

Hox Proteins Act as Transcriptional Guarantors to Ensure Terminal Differentiation

Cell differentiation usually occurs with high fidelity, but the expression of many transcription factors is variable. Using the touch receptor neurons (TRNs) in C. elegans, we found that the Hox proteins CEH-13/lab and EGL-5/Abd-B overcome this variability by facilitating the activation of the common TRN fate determinant mec-3 in the anterior and posterior TRNs, respectively. CEH-13 and EGL-5 increase the probability of mec-3 transcriptional activation by the POU-homeodomain transcription factor UNC-86 using the same Hox/Pbx binding site. Mutation of ceh-13 and egl-5 resulted in an incomplete (∼40%) loss of the TRN fate in respective TRNs, which correlates with quantitative mRNA measurements showing two distinct modes (all or none) of mec-3 transcription. Therefore, Hox proteins act as transcriptional "guarantors" in order to ensure reliable and robust gene expression during terminal neuronal differentiation. Guarantors do not activate gene expression by themselves but promote full activation of target genes regulated by other transcription factors.

Molecular Analysis of the HOXA2-Dependent Degradation of RCHY1

The homeodomain transcription factor Hoxa2 interacts with the RING-finger type E3 ubiquitin ligase RCHY1 and induces its proteasomal degradation. In this work, we dissected this non-transcriptional activity of Hoxa2 at the molecular level. The Hoxa2-mediated decay of RCHY1 involves both the 19S and 20S proteasome complexes. It relies on both the Hoxa2 homeodomain and C-terminal moiety although no single deletion in the Hoxa2 sequence could disrupt the RCHY1 interaction. That the Hoxa2 homeodomain alone could mediate RCHY1 binding is consistent with the shared ability all the Hox proteins we tested to interact with RCHY1. Nonetheless, the ability to induce RCHY1 degradation although critically relying on the homeodomain is not common to all Hox proteins. This identifies the homeodomain as necessary but not sufficient for what appears to be an almost generic Hox protein activity. Finally we provide evidence that the Hoxa2-induced degradation of RCHY1 is evolutionarily conserved among vertebrates. These data therefore support the hypothesis that the molecular and functional interaction between Hox proteins and RCHY1 is an ancestral Hox property.

Molecular dissection of segment formation in the developing hindbrain

Although many components of the genetic pathways that provide positional information during embryogenesis have been identified, it remains unclear how these signals are integrated to specify discrete tissue territories. Here, we investigate the molecular mechanisms underlying the formation of one of the hindbrain segments, rhombomere (r) 3, specified by the expression of the gene krox20. Dissecting krox20 transcriptional regulation has identified several input pathways: Hox paralogous 1 (PG1) factors, which both directly activate krox20 and indirectly repress it via Nlz factors, and the molecular components of an Fgf-dependent effector pathway. These different inputs are channelled through a single initiator enhancer element to shape krox20 initial transcriptional response: Hox PG1 and Nlz factors define the anterior-posterior extent of the enhancer's domain of activity, whereas Fgf signalling modulates the magnitude of activity in a spatially uniform manner. Final positioning of r3 boundaries requires interpretation of this initial pattern by a krox20 positive-feedback loop, orchestrated by another enhancer. Overall, this study shows how positional information provided by different patterning mechanisms is integrated through a gene regulatory network involving two cis-acting elements operating on the same gene, thus offering a comprehensive view of the delimitation of a territory.

The gene regulatory networks underlying formation of the auditory hindbrain

Development and evolution of auditory hindbrain nuclei are two major unsolved issues in hearing research. Recent characterization of transgenic mice identified the rhombomeric origins of mammalian auditory nuclei and unraveled genes involved in their formation. Here, we provide an overview on these data by assembling them into rhombomere-specific gene regulatory networks (GRNs), as they underlie developmental and evolutionary processes. To explore evolutionary mechanisms, we compare the GRNs operating in the mammalian auditory hindbrain with data available from the inner ear and other vertebrate groups. Finally, we propose that the availability of genomic sequences from all major vertebrate taxa and novel genetic techniques for non-model organisms provide an unprecedented opportunity to investigate development and evolution of the auditory hindbrain by comparative molecular approaches. The dissection of the molecular mechanisms leading to auditory structures will also provide an important framework for auditory processing disorders, a clinical problem difficult to tackle so far. These data will, therefore, foster basic and clinical hearing research alike.

Aberrantly hypermethylated Homeobox A2 derepresses metalloproteinase-9 through TBP and promotes invasion in Nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is notorious for its high invasiveness and metastatic ability. In this study, we identified a differential hypermethylated transcription repressor, Homeobox A2 (HOXA2), which may render NPC cells invasive and metastatic. Aberrant hypermethylation of HOXA2 led to low RNA expression in NPC tumors and cells. Addition of methylation inhibitor 5'Aza restored HOXA2 RNA expression in NPC cells. Methylated HOXA2 promoter reduces the binding affinity of the transcriptional co-activator p300, causing transcriptional repression of HOXA2. In NPC cells, re-expression of ectopic HOXA2 was correlated with decreased invasive ability and reduced metalloproteinase MMP-9 RNA and protein expression. Promoter, ChIP and DNA-pull down assays indicated that HOXA2 competes with the transcription activator, TATA-box binding protein (TBP) for a recognition sequence near the MMP-9 transcription start site, and suppresses MMP-9 transcription. Thus, HOXA2 acts as a suppressor or TBP-antagonist to inhibit MMP-9 expression; while methylation-mediated inactivation of HOXA2 in NPC derepresses MMP-9 production and increases invasion of NPC cells. In NPC plasma samples, increased plasma EBV copy number was correlated with increased in cell-free HOXA2 hypermethylation and elevated MMP-9 levels. Plasma EBV DNA and methylated cell-free HOXA2 can be used as biomarkers for monitoring NPC treatment.

Hox proteins drive cell segregation and non-autonomous apical remodelling during hindbrain segmentation

Hox genes encode a conserved family of homeodomain transcription factors regulating development along the major body axis. During embryogenesis, Hox proteins are expressed in segment-specific patterns and control numerous different segment-specific cell fates. It has been unclear, however, whether Hox proteins drive the epithelial cell segregation mechanism that is thought to initiate the segmentation process. Here, we investigate the role of vertebrate Hox proteins during the partitioning of the developing hindbrain into lineage-restricted units called rhombomeres. Loss-of-function mutants and ectopic expression assays reveal that Hoxb4 and its paralogue Hoxd4 are necessary and sufficient for cell segregation, and for the most caudal rhombomere boundary (r6/r7). Hox4 proteins regulate Eph/ephrins and other cell-surface proteins, and can function in a non-cell-autonomous manner to induce apical cell enlargement on both sides of their expression border. Similarly, other Hox proteins expressed at more rostral rhombomere interfaces can also regulate Eph/ephrins, induce apical remodelling and drive cell segregation in ectopic expression assays. However, Krox20, a key segmentation factor expressed in odd rhombomeres (r3 and r5), can largely override Hox proteins at the level of regulation of a cell surface target, Epha4. This study suggests that most, if not all, Hox proteins share a common potential to induce cell segregation but in some contexts this is masked or modulated by other transcription factors.

Hox genes and region-specific sensorimotor circuit formation in the hindbrain and spinal cord

Homeobox (Hox) genes were originally discovered in the fruit fly Drosophila, where they function through a conserved homeodomain as transcriptional regulators to control embryonic morphogenesis. In vertebrates, 39 Hox genes have been identified and like their Drosophila counterparts they are organized within chromosomal clusters. Hox genes interact with various cofactors, such as the TALE homeodomain proteins, in recognition of consensus sequences within regulatory elements of their target genes. In vertebrates, Hox genes display spatially restricted patterns of expression within the developing hindbrain and spinal cord, and are considered crucial determinants of segmental identity and cell specification along the anterioposterior and dorsoventral axes of the embryo. Here, we review their later roles in the assembly of neuronal circuitry, in stereotypic neuronal migration, axon pathfinding, and topographic connectivity. Importantly, we will put some emphasis on how their early-segmented expression patterns can influence the formation of complex vital hindbrain and spinal cord circuitries.

Hox genes: choreographers in neural development, architects of circuit organization

The neural circuits governing vital behaviors, such as respiration and locomotion, are comprised of discrete neuronal populations residing within the brainstem and spinal cord. Work over the past decade has provided a fairly comprehensive understanding of the developmental pathways that determine the identity of major neuronal classes within the neural tube. However, the steps through which neurons acquire the subtype diversities necessary for their incorporation into a particular circuit are still poorly defined. Studies on the specification of motor neurons indicate that the large family of Hox transcription factors has a key role in generating the subtypes required for selective muscle innervation. There is also emerging evidence that Hox genes function in multiple neuronal classes to shape synaptic specificity during development, suggesting a broader role in circuit assembly. This Review highlights the functions and mechanisms of Hox gene networks and their multifaceted roles during neuronal specification and connectivity.

Adaptive evolution of the Hox gene family for development in bats and dolphins

Bats and cetaceans (i.e., whales, dolphins, porpoises) are two kinds of mammals with unique locomotive styles and occupy novel niches. Bats are the only mammals capable of sustained flight in the sky, while cetaceans have returned to the aquatic environment and are specialized for swimming. Associated with these novel adaptations to their environment, various development changes have occurred to their body plans and associated structures. Given the importance of Hox genes in many aspects of embryonic development, we conducted an analysis of the coding regions of all Hox gene family members from bats (represented by Pteropus vampyrus, Pteropus alecto, Myotis lucifugus and Myotis davidii) and cetaceans (represented by Tursiops truncatus) for adaptive evolution using the available draft genome sequences. Differences in the selective pressures acting on many Hox genes in bats and cetaceans were found compared to other mammals. Positive selection, however, was not found to act on any of the Hox genes in the common ancestor of bats and only upon Hoxb9 in cetaceans. PCR amplification data from additional bat and cetacean species, and application of the branch-site test 2, showed that the Hoxb2 gene within bats had significant evidence of positive selection. Thus, our study, with genomic and newly sequenced Hox genes, identifies two candidate Hox genes that may be closely linked with developmental changes in bats and cetaceans, such as those associated with the pancreatic, neuronal, thymus shape and forelimb. In addition, the difference in our results from the genome-wide scan and newly sequenced data reveals that great care must be taken in interpreting results from draft genome data from a limited number of species, and deep genetic sampling of a particular clade is a powerful tool for generating complementary data to address this limitation.

Opposing roles for Hoxa2 and Hoxb2 in hindbrain oligodendrocyte patterning

Oligodendrocytes are the myelin-forming cells of the vertebrate CNS. Little is known about the molecular control of region-specific oligodendrocyte development. Here, we show that oligodendrogenesis in the mouse rostral hindbrain, which is organized in a metameric series of rhombomere-derived (rd) territories, follows a rhombomere-specific pattern, with extensive production of oligodendrocytes in the pontine territory (r4d) and delayed and reduced oligodendrocyte production in the prepontine region (r2d, r3d). We demonstrate that segmental organization of oligodendrocytes is controlled by Hox genes, namely Hoxa2 and Hoxb2. Specifically, Hoxa2 loss of function induced a dorsoventral enlargement of the Olig2/Nkx2.2-expressing oligodendrocyte progenitor domain, whereas conditional Hoxa2 overexpression in the Olig2(+) domain inhibited oligodendrogenesis throughout the brain. In contrast, Hoxb2 deletion resulted in a reduction of the pontine oligodendrogenic domain. Compound Hoxa2(-/-)/Hoxb2(-/-) mutant mice displayed the phenotype of Hoxb2(-/-) mutants in territories coexpressing Hoxa2 and Hoxb2 (rd3, rd4), indicating that Hoxb2 antagonizes Hoxa2 during rostral hindbrain oligodendrogenesis. This study provides the first in vivo evidence that Hox genes determine oligodendrocyte regional identity in the mammalian brain.