Wnt target gene Ascl4 is dispensable for skin appendage development.

The development of skin appendages, including hair follicles, teeth and mammary glands is initiated through the formation of the placode, a local thickening of the epithelium. The Wnt/ß-catenin signaling cascade is an evolutionary conserved pathway with an essential role in placode morphogenesis, but its downstream targets and their exact functions remain ill defined. In this study, we identify Achaete-scute complex-like 4 (Ascl4) as a novel target of the Wnt/ß-catenin pathway and demonstrate its expression pattern in the signaling centers of developing hair follicles and teeth. Ascl transcription factors belong to the superfamily of basic helix-loop-helix transcriptional regulators involved in cell fate determination in many tissues. However, their specific role in the developing skin remains largely unknown. We report that Ascl4 null mice have no overt phenotype. Absence of Ascl4 did not impair hair follicle morphogenesis or hair shaft formation suggesting that it is non-essential for hair follicle development. No tooth or mammary gland abnormalities were detected either. We suggest that other transcription factors may functionally compensate for the absence of Ascl4, but further research is warranted to assess this possibility.

Mesenchyme instructs growth while epithelium directs branching in the mouse mammary gland.

The mammary gland is a unique organ that undergoes dynamic alterations throughout a female's reproductive life, making it an ideal model for developmental, stem cell and cancer biology research. Mammary gland development begins in utero and proceeds via a quiescent bud stage before the initial outgrowth and subsequent branching morphogenesis. How mammary epithelial cells transit from quiescence to an actively proliferating and branching tissue during embryogenesis and, importantly, how the branch pattern is determined remain largely unknown. Here, we provide evidence indicating that epithelial cell proliferation and onset of branching are independent processes, yet partially coordinated by the Eda signaling pathway. Through heterotypic and heterochronic epithelial-mesenchymal recombination experiments between mouse mammary and salivary gland tissues and ex vivo live imaging, we demonstrate that unlike previously concluded, the mode of branching is an intrinsic property of the mammary epithelium whereas the pace of growth and the density of ductal tree are determined by the mesenchyme. Transcriptomic profiling and ex vivo and in vivo functional studies in mice disclose that mesenchymal Wnt/ß-catenin signaling, and in particular IGF-1 downstream of it critically regulate mammary gland growth. These results underscore the general need to carefully deconstruct the different developmental processes producing branched organs.

Vole genomics links determinate and indeterminate growth of teeth.

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Anatomical variations in the cerebral arterial circle of the Saimaa (Pusa hispida saimensis) and Baltic ringed seals (Pusa hispida botnica).

The intracranial arterial vascularization of the Saimaa ringed seals (Pusa hispida saimensis; Nordquist, 1899) and Baltic ringed seals (Pusa hispida botnica; Gmelin, 1788) disclosed patterns of anatomical architecture comparable to that of other pinniped species. Arterial silicone casts on skull scaffolds, and magnetic resonance imaging (MRI) showed that the besides joining the caudal communicating arteries upon entering the cerebral arterial circle, the bilateral internal carotid arteries bifurcated as laterally oriented rostral choroidal arteries and rostral cerebral arteries. The latter arteries almost immediately gave off the laterally oriented middle cerebral arteries. Numerous individual variations were evident in differences in the exact branching sites of bilateral vessels or the size or number of arterial branches. Two Saimaa ringed seals had only a tiny foramen for the left internal carotid artery to enter the intracranial space, and the intracranial part of this vessel was short. It did not reach the cerebral arterial circle. The intracranial part of the right internal carotid artery is bifurcated and also supplied the left side of the cerebral arterial circle. Both specimens had aplasia of the left rostral cerebral artery. The intracranial arterial arrangement of Saimaa and Baltic ringed seals reflects the arterial architecture of this body region in terrestrial mammals with little evidence for aquatic adaptations or changes related to thermoregulation.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord.

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

The developmental basis for scaling of mammalian tooth size.

When evolution leads to differences in body size, organs generally scale along. A well-known example of the tight relationship between organ and body size is the scaling of mammalian molar teeth. To investigate how teeth scale during development and evolution, we compared molar development from initiation through final size in the mouse and the rat. Whereas the linear dimensions of the rat molars are twice that of the mouse molars, their shapes are largely the same. Here, we focus on the first lower molars that are considered the most reliable dental proxy for size-related patterns due to their low within-species variability. We found that scaling of the molars starts early, and that the rat molar is patterned equally as fast but in a larger size than the mouse molar. Using transcriptomics, we discovered that a known regulator of body size, insulin-like growth factor 1 (Igf1), is more highly expressed in the rat molars compared to the mouse molars. Ex vivo and in vivo mouse models demonstrated that modulation of the IGF pathway reproduces several aspects of the observed scaling process. Furthermore, analysis of IGF1-treated mouse molars and computational modeling indicate that IGF signaling scales teeth by simultaneously enhancing growth and by inhibiting the cusp-patterning program, thereby providing a relatively simple mechanism for scaling teeth during development and evolution. Finally, comparative data from shrews to elephants suggest that this scaling mechanism regulates the minimum tooth size possible, as well as the patterning potential of large teeth.

Gene expression detection in developing mouse tissue using in situ hybridization and µCT imaging.

High resolution and noninvasiveness have made soft-tissue X-ray microtomography (µCT) a widely applicable three-dimensional (3D) imaging method in studies of morphology and development. However, scarcity of molecular probes to visualize gene activity with µCT has remained a challenge. Here, we apply horseradish peroxidase-assisted reduction of silver and catalytic gold enhancement of the silver deposit to in situ hybridization in order to detect gene expression in developing tissues with µCT (here called GECT, gene expression CT). We show that GECT detects expression patterns of collagen type II alpha 1 and sonic hedgehog in developing mouse tissues comparably with an alkaline phosphatase-based detection method. After detection, expression patterns are visualized with laboratory µCT, demonstrating that GECT is compatible with varying levels of gene expression and varying sizes of expression regions. Additionally, we show that the method is compatible with prior phosphotungstic acid staining, a conventional contrast staining approach in µCT imaging of soft tissues. Overall, GECT is a method that can be integrated with existing laboratory routines to obtain spatially accurate 3D detection of gene expression.

Transcriptomic landscape of early hair follicle and epidermal development.

Morphogenesis of ectodermal organs, such as hair, tooth, and mammary gland, starts with the formation of local epithelial thickenings, or placodes, but it remains to be determined how distinct cell types and differentiation programs are established during ontogeny. Here, we use bulk and single-cell transcriptomics and pseudotime modeling to address these questions in developing hair follicles and epidermis and produce a comprehensive transcriptomic profile of cellular populations in the hair placode and interplacodal epithelium. We report previously unknown cell populations and marker genes, including early suprabasal and genuine interfollicular basal markers, and propose the identity of suprabasal progenitors. By uncovering four different hair placode cell populations organized in three spatially distinct areas, with fine gene expression gradients between them, we posit early biases in cell fate establishment. This work is accompanied by a readily accessible online tool to stimulate further research on skin appendages and their progenitors.

Toward a universal measure of robustness across model organs and systems.

The development of an individual must be capable of resisting the harmful effects of internal and external perturbations. This capacity, called robustness, can make the difference between normal variation and disease. Some systems and organs are more resilient in their capacity to correct the effects of internal disturbances such as mutations. Similarly, organs and organisms differ in their capacity to be resilient against external disturbances, such as changes in temperature. Furthermore, all developmental systems must be somewhat flexible to permit evolutionary change, and understanding robustness requires a comparative framework. Over the last decades, most research on developmental robustness has been focusing on specific model systems and organs. Hence, we lack tools that would allow cross-species and cross-organ comparisons. Here, we emphasize the need for a uniform framework to experimentally test and quantify robustness across study systems and suggest that the analysis of fluctuating asymmetry might be a powerful proxy to do so. Such a comparative framework will ultimately help to resolve why and how organs of the same and different species differ in their sensitivity to internal (e.g., mutations) and external (e.g., temperature) perturbations and at what level of biological organization buffering capacities exist and therefore create robustness of the developmental system.

Fragmented habitat compensates for the adverse effects of genetic bottleneck.

In the face of the human-caused biodiversity crisis, understanding the theoretical basis of conservation efforts of endangered species and populations has become increasingly important. According to population genetics theory, population subdivision helps organisms retain genetic diversity, crucial for adaptation in a changing environment. Habitat topography is thought to be important for generating and maintaining population subdivision, but empirical cases are needed to test this assumption. We studied Saimaa ringed seals, landlocked in a labyrinthine lake and recovering from a drastic bottleneck, with additional samples from three other ringed seal subspecies. Using whole-genome sequences of 145 seals, we analyzed the distribution of variation and genetic relatedness among the individuals in relation to the habitat shape. Despite a severe history of genetic bottlenecks with prevalent homozygosity in Saimaa ringed seals, we found evidence for the population structure mirroring the subregions of the lake. Our genome-wide analyses showed that the subpopulations had retained unique variation and largely complementary patterns of homozygosity, highlighting the significance of habitat connectivity in conservation biology and the power of genomic tools in understanding its impact. The central role of the population substructure in preserving genetic diversity at the metapopulation level was confirmed by simulations. Integration of genetic analyses in conservation decisions gives hope to Saimaa ringed seals and other endangered species in fragmented habitats.

Chromosomal neighbourhoods allow identification of organ specific changes in gene expression.

Although most genes share their chromosomal neighbourhood with other genes, distribution of genes has not been explored in the context of individual organ development; the common focus of developmental biology studies. Because developmental processes are often associated with initially subtle changes in gene expression, here we explored whether neighbouring genes are informative in the identification of differentially expressed genes. First, we quantified the chromosomal neighbourhood patterns of genes having related functional roles in the mammalian genome. Although the majority of protein coding genes have at least five neighbours within 1 Mb window around each gene, very few of these neighbours regulate development of the same organ. Analyses of transcriptomes of developing mouse molar teeth revealed that whereas expression of genes regulating tooth development changes, their neighbouring genes show no marked changes, irrespective of their level of expression. Finally, we test whether inclusion of gene neighbourhood in the analyses of differential expression could provide additional benefits. For the analyses, we developed an algorithm, called DELocal that identifies differentially expressed genes by comparing their expression changes to changes in adjacent genes in their chromosomal regions. Our results show that DELocal removes detection bias towards large changes in expression, thereby allowing identification of even subtle changes in development. Future studies, including the detection of differential expression, may benefit from, and further characterize the significance of gene-gene neighbour relationships.

System-level analyses of keystone genes required for mammalian tooth development.

When a null mutation of a gene causes a complete developmental arrest, the gene is typically considered essential for life. Yet, in most cases, null mutations have more subtle effects on the phenotype. Here we used the phenotypic severity of mutations as a tool to examine system-level dynamics of gene expression. We classify genes required for the normal development of the mouse molar into different categories that range from essential to subtle modification of the phenotype. Collectively, we call these the developmental keystone genes. Transcriptome profiling using microarray and RNAseq analyses of patterning stage mouse molars show highly elevated expression levels for genes essential for the progression of tooth development, a result reminiscent of essential genes in single-cell organisms. Elevated expression levels of progression genes were also detected in developing rat molars, suggesting evolutionary conservation of this system-level dynamics. Single-cell RNAseq analyses of developing mouse molars reveal that even though the size of the expression domain, measured in the number of cells, is the main driver of organ-level expression, progression genes show high cell-level transcript abundances. Progression genes are also upregulated within their pathways, which themselves are highly expressed. In contrast, a high proportion of the genes required for normal tooth patterning are secreted ligands that are expressed in fewer cells than their receptors and intracellular components. Overall, even though expression patterns of individual genes can be highly different, conserved system-level principles of gene expression can be detected using phenotypically defined gene categories.

Mapping molar shapes on signaling pathways.

A major challenge in evolutionary developmental biology is to understand how genetic mutations underlie phenotypic changes. In principle, selective pressures on the phenotype screen the gene pool of the population. Teeth are an excellent model for understanding evolutionary changes in the genotype-phenotype relationship since they exist throughout vertebrates. Genetically modified mice (mutants) with abnormalities in teeth have been used to explore tooth development. The relationship between signaling pathways and molar shape, however, remains elusive due to the high intrinsic complexity of tooth crowns. This hampers our understanding of the extent to which developmental factors explored in mutants explain developmental and phenotypic variation in natural species that represent the consequence of natural selection. Here we combine a novel morphometric method with two kinds of data mining techniques to extract data sets from the three-dimensional surface models of lower first molars: i) machine learning to maximize classification accuracy of 22 mutants, and ii) phylogenetic signal for 31 Murinae species. Major shape variation among mutants is explained by the number of cusps and cusp distribution on a tooth crown. The distribution of mutant mice in morphospace suggests a nonlinear relationship between the signaling pathways and molar shape variation. Comparative analysis of mutants and wild murines reveals that mutant variation overlaps naturally occurring diversity, including more ancestral and derived morphologies. However, taxa with transverse lophs are not fully covered by mutant variation, suggesting experimentally unexplored developmental factors in the evolutionary radiation of Murines.

Reptile-like physiology in Early Jurassic stem-mammals.

Despite considerable advances in knowledge of the anatomy, ecology and evolution of early mammals, far less is known about their physiology. Evidence is contradictory concerning the timing and fossil groups in which mammalian endothermy arose. To determine the state of metabolic evolution in two of the earliest stem-mammals, the Early Jurassic Morganucodon and Kuehneotherium, we use separate proxies for basal and maximum metabolic rate. Here we report, using synchrotron X-ray tomographic imaging of incremental tooth cementum, that they had maximum lifespans considerably longer than comparably sized living mammals, but similar to those of reptiles, and so they likely had reptilian-level basal metabolic rates. Measurements of femoral nutrient foramina show Morganucodon had blood flow rates intermediate between living mammals and reptiles, suggesting maximum metabolic rates increased evolutionarily before basal metabolic rates. Stem mammals lacked the elevated endothermic metabolism of living mammals, highlighting the mosaic nature of mammalian physiological evolution.

The effect of amoxicillin on dental enamel development in vivo.

The exposure to amoxicillin has been associated with molar incisor hypomineralization. This study aimed to determine if amoxicillin disturbs the enamel mineralization in in vivo experiments. Fifteen pregnant rats were randomly assigned into three groups to received daily phosphatase-buffered saline or amoxicillin as either 100 or 500 mg/kg. Mice received treatment from day 13 of pregnancy to day 40 postnatal. After birth, the offsprings from each litter continued to receive the same treatment according to their respective group. Calcium (Ca) and phosphorus (P) content in the dental hard tissues were analyzed from 60 upper first molars and 60 upper incisors by the complexometric titration method and colorimetric analysis using a spectrophotometer at 680 nm, respectively. Lower incisors were analyzed by X-ray microtomography, it was measured the electron density of lingual and buccal enamel, and the enamel and dentin thickness. Differences in Ca and P content and electron density among the groups were analyzed by one-way ANOVA. There was no significant difference on enamel electron density and thickness among the groups (p > 0.05). However, in incisors, the higher dose of amoxicillin decreased markedly the electron density in some rats. There were no statistically significant differences in Ca (p = 0.180) or P content (p = 0.054), although the higher dose of amoxicillin could affect the enamel in some animals. The amoxicillin did not significantly alter the enamel mineralization and thickness in rats.

gapFinisher: A reliable gap filling pipeline for SSPACE-LongRead scaffolder output.

Unknown sequences, or gaps, are present in many published genomes across public databases. Gap filling is an important finishing step in de novo genome assembly, especially in large genomes. The gap filling problem is nontrivial and while there are many computational tools partially solving the problem, several have shortcomings as to the reliability and correctness of the output, i.e. the gap filled draft genome. SSPACE-LongRead is a scaffolding tool that utilizes long reads from multiple third-generation sequencing platforms in finding links between contigs and combining them. The long reads potentially contain sequence information to fill the gaps created in the scaffolding, but SSPACE-LongRead currently lacks this functionality. We present an automated pipeline called gapFinisher to process SSPACE-LongRead output to fill gaps after the scaffolding. gapFinisher is based on the controlled use of a previously published gap filling tool FGAP and works on all standard Linux/UNIX command lines. We compare the performance of gapFinisher against two other published gap filling tools PBJelly and GMcloser. We conclude that gapFinisher can fill gaps in draft genomes quickly and reliably. In addition, the serial design of gapFinisher makes it scale well from prokaryote genomes to larger genomes with no increase in the computational footprint.

Modeling enamel matrix secretion in mammalian teeth.

The most mineralized tissue of the mammalian body is tooth enamel. Especially in species with thick enamel, three-dimensional (3D) tomography data has shown that the distribution of enamel varies across the occlusal surface of the tooth crown. Differences in enamel thickness among species and within the tooth crown have been used to examine taxonomic affiliations, life history, and functional properties of teeth. Before becoming fully mineralized, enamel matrix is secreted on the top of a dentine template, and it remains to be explored how matrix thickness is spatially regulated. To provide a predictive framework to examine enamel distribution, we introduce a computational model of enamel matrix secretion that maps the dentine topography to the enamel surface topography. Starting from empirical enamel-dentine junctions, enamel matrix deposition is modeled as a diffusion-limited free boundary problem. Using laboratory microCT and synchrotron tomographic data of pig molars that have markedly different dentine and enamel surface topographies, we show how diffusion-limited matrix deposition accounts for both the process of matrix secretion and the final enamel distribution. Simulations reveal how concave and convex dentine features have distinct effects on enamel surface, thereby explaining why the enamel surface is not a straightforward extrapolation of the dentine template. Human and orangutan molar simulations show that even subtle variation in dentine topography can be mapped to the enamel surface features. Mechanistic models of extracellular matrix deposition can be used to predict occlusal morphologies of teeth.

Bracketing phenogenotypic limits of mammalian hybridization.

An increasing number of mammalian species have been shown to have a history of hybridization and introgression based on genetic analyses. Only relatively few fossils, however, preserve genetic material, and morphology must be used to identify the species and determine whether morphologically intermediate fossils could represent hybrids. Because dental and cranial fossils are typically the key body parts studied in mammalian palaeontology, here we bracket the potential for phenotypically extreme hybridizations by examining uniquely preserved cranio-dental material of a captive hybrid between grey and ringed seals. We analysed how distinct these species are genetically and morphologically, how easy it is to identify the hybrids using morphology and whether comparable hybridizations happen in the wild. We show that the genetic distance between these species is more than twice the modern human-Neanderthal distance, but still within that of morphologically similar species pairs known to hybridize. By contrast, morphological and developmental analyses show grey and ringed seals to be highly disparate, and that the hybrid is a predictable intermediate. Genetic analyses of the parent populations reveal introgression in the wild, suggesting that grey-ringed seal hybridization is not limited to captivity. Taken together, we postulate that there is considerable potential for mammalian hybridization between phenotypically disparate taxa.

A microCT Study of Three-Dimensional Patterns of Biomineralization in Pig Molars.

Domestic pig molars provide an interesting system to study the biomineralization process. The large size, thick enamel and complex crown morphology make pig molars relatively similar to human molars. However, compared to human molars, pig molars develop considerably faster. Here we use microCT to image the developing pig molars and to decipher spatial patterns of biomineralization. We used mineral grains to calibrate individual microCT-scans, which allowed an accurate measure of the electron density of the developing molars. The microCT results show that unerupted molars that are morphologically at the same stage of development, can be at markedly different stage of enamel biomineralization. Erupted molars show increased electron density, suggesting that mineralization continues in oral cavity. Yet, our comparisons show that human enamel has slightly higher electron density than pig enamel. These results support the relatively low hardness values and calcium level values that have been reported earlier in literature for pig teeth. The mineral calibration was an efficient method for the microCT-absorption models, allowing a relatively robust way to detect scanning artifacts. In conclusions, whereas thin sections remain the preferred way to analyze enamel features, such as incremental lines and crystal orientation, the microCT allows efficient and non-destructive comparisons between different teeth and species.