To what extent is digit patterning a Turing System?

Building precise, robust patterns and structures from an initially homogeneous state is fundamental to developmental biology. Digit patterning is a representative example of a periodic pattern in development. Previous studies have shown that a reaction–diffusion (Turing) system, in which diffusible activators and inhibitors interact, is the most likely explanation of how the spatial pattern of the digits is formed. Although self-organisation mechanisms such as the Turing system successfully recapitulate many aspects of digit patterning, critical questions remain regarding its timing and behaviour. First I addressed the question of timing, or how long reaction-diffusion plays a role in the developing digits. I perturbed the digit patterning process of embryonic limbs by inserting beads that contain morphogens involved in the reaction-diffusion mechanism. Then I quantified the degree of pattern change, or plasticity of the patterning, from limbs harvested at different developmental timing throughout the digit patterning stage. For quantification, I developed a custom image analytic pipeline that extracts relevant topology and represents the difference between perturbed and unperturbed patterns. Modelling the plasticity profile over the digit patterning process, through extensive interplay of experiments and modelling, revealed that plasticity during digit patterning decreases in a sigmoidal manner. Transcriptomics analysis that matches with the sigmoidal decrease observed in expression patterns further identified gene candidates that could be critical to the digit patterning. Further, the timing of reaction-diffusion is discussed in the context of the tissue movements, revealing that Sox9 digit patterning happens significantly earlier than cell density changes. The second part aims at improving our understanding about which pathways and components of the pathways are involved in the digit forming Turing network. Previously identified digit patterning Turing network, such as BSW model, abstracts the entire Wnt and Bmp signalling pathways’ activities into each node. Thus there is insufficient knowledge on the mechanistic role of Wnt signalling mediated Sox9 repression. To further clarify detailed mechanisms of the Turing network, I used an unbiased screening approach to systematically perturb digit patterning using small molecule inhibitors, ligands, and peptides at different doses in systems such as limb culture and micromass. Out of multiple steps critical to Wnt signalling, including Wnt production, Wnt receptor interaction, Wnt canonical pathway cytosolic interactions, and Wnt canonical pathway transcriptional interactions, I identified that inhibition of Wnt production and Wnt transcriptional component inhibition category most effectively disrupt digit patterning. I also identified candidate ligands such as sFRP1 and Dkk1 as potential extracellular Wnt inhibitors that effectively change digit patterning upon application. These results provide the first quantitative insight into the duration of the reaction-diffusion based mechanism in a biological system, and how a screening approach complemented with data driven modelling can complement and clarify workings of a reaction diffusion based system. Further work in improving our knowledge on the Turing system with tissue growth, cell movements, and ectodermal-mesenchymal interaction will eventually allow generation of a complete organogenesis simulation model

Regulating Retinoic Acid Availability during Development and Regeneration: The Role of the CYP26 Enzymes.

This review focuses on the role of the Cytochrome p450 subfamily 26 (CYP26) retinoic acid (RA) degrading enzymes during development and regeneration. Cyp26 enzymes, along with retinoic acid synthesising enzymes, are absolutely required for RA homeostasis in these processes by regulating availability of RA for receptor binding and signalling. Cyp26 enzymes are necessary to generate RA gradients and to protect specific tissues from RA signalling. Disruption of RA homeostasis leads to a wide variety of embryonic defects affecting many tissues. Here, the function of CYP26 enzymes is discussed in the context of the RA signalling pathway, enzymatic structure and biochemistry, human genetic disease, and function in development and regeneration as elucidated from animal model studies

Emerging properties of animal gene regulatory networks

Gene regulatory networks (GRNs) provide system level explanations of developmental and physiological functions in the terms of the genomic regulatory code. Depending on their developmental functions, GRNs differ in their degree of hierarchy, and also in the types of modular sub-circuit of which they are composed, although there is a commonly employed sub-circuit repertoire. Mathematical modelling of some types of GRN sub-circuit has deepened biological understanding of the functions they mediate. The structural organization of various kinds of GRN reflects their roles in the life process, and causally illuminates both developmental and evolutionary process

Computational Stem Cell Biology: Open Questions and Guiding Principles

Computational biology is enabling an explosive growth in our understanding of stem cells and our ability to use them for disease modeling, regenerative medicine, and drug discovery. We discuss four topics that exemplify applications of computation to stem cell biology: cell typing, lineage tracing, trajectory inference, and regulatory networks. We use these examples to articulate principles that have guided computational biology broadly and call for renewed attention to these principles as computation becomes increasingly important in stem cell biology. We also discuss important challenges for this field with the hope that it will inspire more to join this exciting area

Sox9 Expression in Amniotes: Species-Specific Differences in the Formation of Digits

In tetrapods the digit pattern has evolved to adapt to distinct locomotive strategies. The number of digits varies between species or even between hindlimb and forelimb within the same species. These facts illustrate the plasticity of embryonic limb autopods. Sox9 is a precocious marker of skeletal differentiation of limb mesenchymal cells. Its pattern of expression in the developing limb has been widely studied and reflects the activity of signaling cascades responsible for skeletogenesis. In this assay we stress previously overlooked differences in the pattern of expression of Sox9 in limbs of avian, mouse and turtle embryos which may reflect signaling differences associated with distinct limb skeletal morphologies observed in these species. Furthermore, we show that Sox9 gene expression is higher and maintained in the interdigital region in species with webbed digits in comparison with free digit animals

Analysis of craniofacial defects in Six1/Eya1-associated Branchio-Oto-Renal Syndrome

Poster Session I - Morphogenesis: 205/B10117th ISDB 2013 cum 72nd Annual Meeting of the Society for Developmental Biology, 7th Latin American Society of Developmental Biology Meeting and 11th Congreso de la Sociedad Mexicana de Biologia del Desarrollo.Branchio-Oto-Renal (BOR) syndrome patients exhibit craniofacial and renal anomalies as well as deafness. BOR syndrome is caused by mutations in Six1 or Eya1, both of which regulate cell proliferation and differentiation. The molecular mechanism underlying the craniofacial and branchial arch (BA) defects in BOR syndrome is unclear. We have found that Hoxb3 is up-regulated in the second branchial arch (BA2) of Six1-/- mutants. Moreover, Hoxb3 over-expression in transgenic mice leads to BA abnormalities which are similar to the BA defects in Six1-/- or Eya1-/- mutants, suggesting a regulatory relationship among Six1, Eya1 and Hoxb3 genes. The aim of this study is to investigate the molecular mechanism underlying abnormal BA development in BOR syndrome using Six1 and Eya1 mutant mice. Two potential Six1 binding sites were identified on the Hoxb3 gene. In vitro and in vivo Chromatin IP assays showed that Six1 could directly bind to one of the sites specifically. Furthermore, using a chick in ovo luciferase assay we showed that Six1 could suppress gene expression through one of the specific binding sites. On the other hand, in Six1-/- mutants, we found that the Notch ligand Jag1 was up-regulated in BA2. Similarly, in Hoxb3 transgenic mice, ectopic expression of Jag1 could be also detected in BA2. To investigate the activation of Notch signaling pathway, we found that Notch intracellular domain (NICD), a direct indicator of Notch pathway activation, was up-regulated in BAs of Six1-/-; Eya1-/- double mutants. Our results indicate that Hoxb3 and Notch signaling pathway are involved in mediating the craniofacial defects of Six1/Eya1-associated Branchio-Oto-Renal Syndrome.postprin

Sox10 regulates enteric neural crest cell migration in the developing gut

Concurrent Sessions 1: 1.3 - Organs to organisms: Models of Human Diseases: abstract no. 1417th ISDB 2013 cum 72nd Annual Meeting of the Society for Developmental Biology, VII Latin American Society of Developmental Biology Meeting and XI Congreso de la Sociedad Mexicana de Biologia del Desarrollo. The Conference's web site is located at http://www.inb.unam.mx/isdb/Sox10 is a HMG-domain containing transcription factor which plays important roles in neural crest cell survival and differentiation. Mutations of Sox10 have been identified in patients with Waardenburg-Hirschsprung syndrome, who suffer from deafness, pigmentation defects and intestinal aganglionosis. Enteric neural crest cells (ENCCs) with Sox10 mutation undergo premature differentiation and fail to colonize the distal hindgut. It is unclear, however, whether Sox10 plays a role in the migration of ENCCs. To visualize the migration behaviour of mutant ENCCs, we generated a Sox10NGFP mouse model where EGFP is fused to the N-terminal domain of Sox10. Using time-lapse imaging, we found that ENCCs in Sox10NGFP/+ mutants displays lower migration speed and altered trajectories compared to normal controls. This behaviour was cell-autonomous, as shown by organotypic grafting of Sox10NGFP/+ gut segments onto control guts and vice versa. ENCCs encounter different extracellular matrix (ECM) molecules along the developing gut. We performed gut explant culture on various ECM and found that Sox10NGFP/+ ENCCs tend to form aggregates, particularly on fibronectin. Time-lapse imaging of single cells in gut explant culture indicated that the tightly-packed Sox10 mutant cells failed to exhibit contact inhibition of locomotion. We determined the expression of adhesion molecule families by qPCR analysis, and found integrin expression unaffected while L1-cam and selected cadherins were altered, suggesting that Sox10 mutation affects cell adhesion properties of ENCCs. Our findings identify a de novo role of Sox10 in regulating the migration behaviour of ENCCs, which has important implications for the treatment of Hirschsprung disease.postprin

ClockOME: searching for oscillatory genes in early vertebrate development

Embryo development is a dynamic process regulated in space and time. Cells must integrate biochemical and mechanical signals to generate fully functional organisms, where oscillatory gene expression plays a key role. The embryo molecular clock (EMC) is the best known genetic oscillator active in embryo segmentation, involving genes from the Notch, FGF, and WNT pathways. However, the list of cyclic genes is still incomplete mostly due to the challenges involved with studying periodic systems. Recently, such studies have become more feasible with the development of pseudo-time ordering algorithms that search for candidate oscillatory genes using large transcriptomics datasets sampled without explicit time measurements. This study aims at finding candidate oscillatory genes - ClockOME - active in early chick embryo development. Two Gallus gallus microarray transcriptomics datasets from Presomitic mesoderm (PSM), and one dataset from limb segmentation were gathered from GEO and ArrayExpress. To normalize these data from different experiments, an RData package - FrozenChicken - was developed to apply a frozen Robust MultiArray (fRMA) normalization to the data. Next the datasets were processed with Oscope (a pseudo-time ordering algorithm) to search for candidate periodic genes clustered by similar oscillatory behaviour. The clusters of predicted oscillators were then subject to functional enrichment and interaction network analyses to highlight the biological functions associated with these genes. Oscope predicted three clusters of oscillators: two in PSM (106 and 32 genes), and one in Limb (162 genes). Overall, the genes are associated with regulatory, morphological, and developmental processes. Mesp2, a gene involved with the EMC, was found in this dataset, validating the approach, however, the majority of genes are novel oscillatory candidates, associated with chromatin and transcriptional regulation, as well as protein and oxygen metabolism. The list of candidate oscillators represents a valuable resource for guided experimental validation to discover additional members of the chick EMC. Six genes have been proposed for high-priority experimental validation: SRC, PTCH1, NOTCH2, YAP1, KDR, CTR9.O desenvolvimento embrionário é um processo dinâmico que envolve alterações moleculares no espaço e no tempo. As células embrionárias são constantemente expostas a estímulos bioquímicos e mecânicos, e respondem ao ambiente em que se encontram alterando o seu programa genético. Quando corretamente integradas, estas respostas celulares culminam com o desenvolvimento bem-sucedido de um organismo funcional. Assim, a embriogénese envolve processos moleculares estritamente regulados, sendo a expressão oscilatória de genes uma das formas possíveis para a regulação do comportamento das células ao longo do tempo. O relógio molecular embrionário é um conhecido oscilador genético, e está envolvido na segmentação do tecido paraxial embrionário. O conceito de relógio molecular foi inicialmente proposto em 1976 por Cooke e Zeeman, ao qual chamaram o modelo Clock and Wavefront (Relógio e Frente de Onda)1. Este modelo foi concebido para descrever teoricamente a formação rítmica de sómitos em ambos os lados da mesoderme paraxial (PSM) nos vertebrados, e baseia-se na existência de osciladores genéticos que regulam esse processo de segmentação da PSM ao longo do tempo. Para além do relógio, como diz o nome, o modelo inclui a existência de uma frente de onda, que determina espacialmente o comportamento das células presentes na mesoderme pré-somítica (PSM). Assim, os dois mecanismos guiam a diferenciação das células da PSM, que consequentemente sofrem transformações genéticas que precedem a formação dos sómitos. A base deste relógio molecular consiste na expressão periódica de genes que fazem parte das vias moleculares Notch, FGF e WNT. Contudo, a lista de genes envolvidos no relógio embrionário ainda não se encontra completa, facto este que se deve principalmente às dificuldades experimentais relacionadas com o estudo de sistemas periódicos quando não se conhece de antemão a periodicidade/ritmo da expressão dos genes envolvidos. Com o advento de novas técnicas de transcriptómica que permitem o estudo dos valores de expressão de todos os genes simultaneamente, nomeadamente usando Microarrays, ou mais recentemente através de métodos de sequenciação, como RNA-sequencing ou Single-Cell RNA-sequencing, surge a oportunidade de procurar alargar a lista de genes com expressão oscilatória. Porém, estes métodos implicam a extração do RNA das células amostradas resultando na morte celular. Assim, este processamento inviabiliza o estudo das mesmas células ao longo do tempo, originando dados moleculares estáticos, isto é, os níveis de expressão obtidos representam uma única amostra temporal. Para o estudo de processos periódicos, seria então necessário fazer uma série temporal amostrando diferentes indivíduos ao longo do tempo de desenvolvimento, aumentando grandemente o número de amostras biológicas necessárias para resolver o ciclo de oscilação para cada gene estudado. Assim, sem informação temporal medida explicitamente, a expressão oscilatória de genes pode apenas ser estudada usando modelos matemáticos apropriados, nomeadamente através da aplicação de algoritmos de ordenação pseudo-temporal. Estes métodos ordenam as amostras ao longo do tempo de uma oscilação de forma a obter o padrão do comportamento cíclico para todos os genes cuja expressão oscila concomitantemente. Torna-se assim possível, bioinformaticamente, inferir o potencial oscilatório de genes medidos por estas técnicas de transcriptómica, sem informação temporal explícita. Deste modo, o objetivo deste estudo é encontrar novos genes oscilatórios, a que coletivamente chamamos ClockOME, que estão ativos durante as primeiras etapas do desenvolvimento embrionário (somitogénese) da galinha, nos tecidos da mesoderme présomítica (PSM), e no membro superior (Limb); tecidos estes onde o relógio molecular foi descrito, atuando como regulador temporal das alterações genéticas subjacentes. Para tal, recolheu-se 3 conjuntos de dados (datasets) de transcriptómica obtidos por microarray de dois repositórios de dados públicos: GEO (da instituição americana NCBI) e ArrayExpress (da instituição europeia EMBL-EBI). Dois datasets continham dados de mesoderme paraxial (PSM) – tecido onde ocorre a somitogénese; e um dataset de dados de obtidos do membro superior do embrião de galinha. Com o objetivo de normalizar os três datasets de forma a torná-los comparáveis (uma vez que são oriundos de processos experimentais diferentes), foi desenvolvido um pacote de R denominado “FrozenChicken: Promoting the meta-analysis of chicken microarray data” (publicado em 2021) (https://doi.org/10.1101/2021.02.25.432894). Este pacote contém dados sumarizados de 472 datasets de microarrays de embriões de galinha, tornando possível a normalização por fRMA (frozen Robust MultiArray) de microarrays de Gallus gallus. Após normalização e controlo de qualidade dos valores de expressão genética, os dados da PSM e do membro foram processados com o Oscope (algoritmo de ordenação pseudo-temporal), com o propósito de prever genes oscilatórios. Este algoritmo avalia todas as combinações de pares de genes, agrupando aqueles que apresentem padrões de expressão semelhantes, ou seja, cujos valores de expressão ao longo das amostras seguem trajetórias semelhantes, indiciando um período de oscilação potencialmente semelhante. Os clusters de genes previstos pelo Oscope foram posteriormente submetidos a uma análise de enriquecimento funcional e a uma análise de interações funcionais, com o intuito de perceber o seu potencial papel biológico, e funções moleculares subjacentes. O Oscope reportou três listas de genes potencialmente oscilatórios: dois grupos foram encontrados a partir dos dados da PSM (com 106 e 32 genes cada) e o terceiro grupo de 162 genes foi encontrado nos dados do membro superior. No total, a lista de genes que denominamos ClockOME é composta por 296 genes potencialmente oscilatórios, envolvidos em diversos mecanismos regulatórios importantes para o desenvolvimento embrionário e para a morfogénese. A maioria dos genes presentes nesta lista não estão descritos na literatura como sendo oscilatórios (novel candidates), representando, portanto, uma mais-valia para a comunidade científica que estuda o relógio molecular embrionário. Estes genes parecem estar associados a funções como remodelação da cromatina, regulação da transcrição, metabolismo proteico e metabolismo do oxigénio, sendo, portanto, bons candidatos para futura validação experimental. Notavelmente, o Oscope identificou com sucesso o Mesp2, um gene oscilatório bem descrito na literatura, mostrando assim a validade e o potencial desta abordagem teórica. Em suma, este trabalho produziu uma lista de 296 genes potencialmente oscilatórios. Com base na sua novidade e na função molecular anotada, foi proposta uma lista de seis genes candidatos de particular relevância para validação experimental no futuro próximo, nomeadamente: SRC, PTCH1, NOTCH2, YAP1, KDR, CTR9. Assim, as listas resultantes do trabalho desta tese poderão agora guiar futuras experiências laboratoriais capazes de adicionar novos interactores moleculares ao atual modelo do relógio molecular embrionário

Understanding axial progenitor biology in vivo and in vitro

The generation of the components that make up the embryonic body axis, such as the spinal cord and vertebral column, takes place in an anterior-to-posterior (head-to-tail) direction. This process is driven by the coordinated production of various cell types from a pool of posteriorly-located axial progenitors. Here, we review the key features of this process and the biology of axial progenitors, including neuromesodermal progenitors, the common precursors of the spinal cord and trunk musculature. We discuss recent developments in the in vitro production of axial progenitors and their potential implications in disease modelling and regenerative medicine