Developmental hourglass and heterochronic shifts in fin and limb development

Comparative Hox genes expression within the dimorphic annelid Streblospio benedicti reveals patterning variation during development

Hox genes are transcriptional regulators that elicit cell positional identity along the anterior-posterior region of the body plan across different lineages of Metazoan. Comparison of Hox gene expression across distinct species reveals their evolutionary conservation; however, their gains and losses in different lineages can correlate with body plan modifications and morphological novelty. We compare the expression of 11 Hox genes found within Streblospio benedicti, a marine annelid that produces two types of offspring with distinct developmental and morphological features. For these two distinct larval types, we compare Hox gene expression through ontogeny using hybridization chain reaction (HCR) probes for in situ hybridization and RNA-seq data. We find that Hox gene expression patterning for both types is typically similar at equivalent developmental stages. However, some Hox genes have spatial or temporal differences between the larval types that are associated with morphological and life-history differences. This is the first comparison of developmental divergence in Hox gene expression within a single species and these changes reveal how body plan differences may arise in larval evolution.

Splicing is dynamically regulated during limb development

Modifications to highly conserved developmental gene regulatory networks are thought to underlie morphological diversification in evolution and contribute to human congenital malformations. Relationships between gene expression and morphology have been extensively investigated in the limb, where most of the evidence for alterations to gene regulation in development consists of pre-transcriptional mechanisms that affect expression levels, such as epigenetic alterations to regulatory sequences and changes to cis-regulatory elements. Here we report evidence that alternative splicing (AS), a post-transcriptional process that modifies and diversifies mRNA transcripts, is dynamic during limb development in two mammalian species. We evaluated AS patterns in mouse (Mus musculus) and opossum (Monodelphis domestica) across the three key limb developmental stages: the ridge, bud, and paddle. Our data show that splicing patterns are dynamic over developmental time and suggest differences between the two mammalian taxa. Additionally, multiple key limb development genes, including Fgf8, are differentially spliced across the three stages in both species, with expression levels of the conserved splice variants, Fgf8a and Fgf8b, changing across developmental time. Our data demonstrates that AS is a critical mediator of mRNA diversity in limb development and provides an additional mechanism for evolutionary tweaking of gene dosage.

Paired fins in vertebrate evolution and ontogeny

The origin of paired appendages became one of the most important adaptations of vertebrates, allowing them to lead active lifestyles and explore a wide range of ecological niches. The basic form of paired appendages in evolution is the fins of fishes. The problem of paired appendages has attracted the attention of researchers for more than 150 years. During this time, a number of theories have been proposed, mainly based on morphological data, two of which, the Balfour-Thacher-Mivart lateral fold theory and Gegenbaur's gill arch theory, have not lost their relevance. So far, however, none of the proposed ideas has been supported by decisive evidence. The study of the evolutionary history of the appearance and development of paired appendages lies at the intersection of several disciplines and involves the synthesis of paleontological, morphological, embryological, and genetic data. In this review, we attempt to summarize and discuss the results accumulated in these fields and to analyze the theories put forward regarding the prerequisites and mechanisms that gave rise to paired fins and limbs in vertebrates.

The contributions of DNA accessibility and transcription factor occupancy to enhancer activity during cellular differentiation

During gene regulation, DNA accessibility is thought to limit the availability of transcription factor (TF) binding sites, while TFs can increase DNA accessibility to recruit additional factors that upregulate gene expression. Given this interplay, the causative regulatory events in the modulation of gene expression remain unknown for the vast majority of genes. We utilized deeply sequenced ATAC-Seq data and site-specific knock-in reporter genes to investigate the relationship between the binding-site resolution dynamics of DNA accessibility and the expression dynamics of the enhancers of Cebpa during macrophage-neutrophil differentiation. While the enhancers upregulate reporter expression during the earliest stages of differentiation, there is little corresponding increase in their total accessibility. Conversely, total accessibility peaks during the last stages of differentiation without any increase in enhancer activity. The accessibility of positions neighboring C/EBP-family TF binding sites, which indicates TF occupancy, does increase significantly during early differentiation, showing that the early upregulation of enhancer activity is driven by TF binding. These results imply that a generalized increase in DNA accessibility is not sufficient, and binding by enhancer-specific TFs is necessary, for the upregulation of gene expression. Additionally, high-coverage ATAC-Seq combined with time-series expression data can infer the sequence of regulatory events at binding-site resolution.

Comparative Hox genes expression within the dimorphic annelid Streblospio benedicti reveals patterning variation during development

Hox genes are transcriptional regulators that elicit cell positional identity along the anterior-posterior region of the body plan across different lineages of Metazoan. Comparison of Hox gene expression across distinct species reveals their evolutionary conservation, however their gains and losses in different lineages can correlate with body plan modifications and morphological novelty. We compare the expression of eleven Hox genes found within Streblospio benedicti, a marine annelid that produces two types of offspring with distinct developmental and morphological features. For these two distinct larval types, we compare Hox gene expression through ontogeny using HCR (hybridization chain reaction) probes for in-situ hybridization and RNA-seq data. We find that Hox gene expression patterning for both types is typically similar at equivalent developmental stages. However, some Hox genes have spatial or temporal differences between the larval types that are associated with morphological and life-history differences. This is the first comparison of developmental divergence in Hox genes expression within a single species and these changes reveal how body plan differences may arise in larval evolution.

The contributions of DNA accessibility and transcription factor occupancy to enhancer activity during cellular differentiation

The upregulation of gene expression by enhancers depends upon the interplay between the binding of sequence-specific transcription factors (TFs) and DNA accessibility. DNA accessibility is thought to limit the ability of TFs to bind to their sites, while TFs can increase accessibility to recruit additional factors that upregulate gene expression. Given this interplay, the causative regulatory events underlying the modulation of gene expression during cellular differentiation remain unknown for the vast majority of genes. We investigated the binding-site resolution dynamics of DNA accessibility and the expression dynamics of the enhancers of an important neutrophil gene, Cebpa, during macrophage-neutrophil differentiation. Reporter genes were integrated in a site-specific manner in PUER cells, which are progenitors that can be differentiated into neutrophils or macrophages in vitro by activating the pan-leukocyte TF PU.1. Time series data show that two enhancers upregulate reporter expression during the first 48 hours of neutrophil differentiation. Surprisingly, there is little or no increase in the total accessibility, measured by ATAC-Seq, of the enhancers during the same time period. Conversely, total accessibility peaks 96 hrs after PU.1 activation-consistent with its role as a pioneer-but the enhancers do not upregulate gene expression. Combining deeply sequenced ATAC-Seq data with a new bias-correction method allowed the profiling of accessibility at single-nucleotide resolution and revealed protected regions in the enhancers that match all previously characterized TF binding sites and ChIP-Seq data. Although the accessibility of most positions does not change during early differentiation, that of positions neighboring TF binding sites, an indicator of TF occupancy, did increase significantly. The localized accessibility changes are limited to nucleotides neighboring C/EBP-family TF binding sites, showing that the upregulation of enhancer activity during early differentiation is driven by C/EBP-family TF binding. These results show that increasing the total accessibility of enhancers is not sufficient for upregulating their activity and other events such as TF binding are necessary for upregulation. Also, TF binding can cause upregulation without a perceptible increase in total accessibility. Finally, this study demonstrates the feasibility of comprehensively mapping individual TF binding sites as footprints using high coverage ATAC-Seq and inferring the sequence of events in gene regulation by combining with time-series gene expression data.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

Identification of Transcription Factor Networks during Mouse Hindlimb Development

Mammalian hindlimb development involves a variety of cells and the regulation of spatiotemporal molecular events, but regulatory networks of transcription factors contributing to hindlimb morphogenesis are not well understood. Here, we identified transcription factor networks during mouse hindlimb morphology establishment through transcriptome analysis. We used four stages of embryonic hindlimb transcription profiles acquired from the Gene Expression Omnibus database (GSE30138), including E10.5, E11.5, E12.5 and E13.5, to construct a gene network using Weighted Gene Co-expression Network Analysis (WGCNA), and defined seven stage-associated modules. After filtering 7625 hub genes, we further prioritized 555 transcription factors with AnimalTFDB3.0. Gene ontology enrichment showed that transcription factors of different modules were enriched in muscle tissue development, connective tissue development, embryonic organ development, skeletal system morphogenesis, pattern specification process and urogenital system development separately. Six regulatory networks were constructed with key transcription factors, which contribute to the development of different tissues. Knockdown of four transcription factors from regulatory networks, including Sox9, Twist1, Snai2 and Klf4, showed that the expression of limb-development-related genes was also inhibited, which indicated the crucial role of transcription factor networks in hindlimb development.

Squalomix: shark and ray genome analysis consortium and its data sharing platform

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.

Early shape divergence of developmental trajectories in the jaw of galeomorph sharks

Background

The onset of morphological differences between related groups can be tracked at early stages during embryological development. This is expressed in functional traits that start with minor variations, but eventually diverge to defined specific morphologies. Several processes during this period, like proliferation, remodelling, and apoptosis for instance, can account for the variability observed between related groups. Morphological divergence through development is often associated with the hourglass model, in which early stages display higher variability and reach a conserved point with reduced variability from which divergence occurs again to the final phenotype.

Results

Here we explored the patterns of developmental shape changes in the lower jaw of two shark species, the bamboo shark (Chiloscyllium punctatum) and the catshark (Scyliorhinus canicula). These two species present marked differences in their foraging behaviour, which is reflected in their adult jaw morphology. By tracing the developmental sequence of the cartilage condensation, we identified the onset of cartilage for both species at around stage 31. Other structures that developed later without a noticeable anlage were the labial cartilages, which appear at around stage 33. We observed that the lower jaw displays striking differences in shape from the earliest moments, without any overlap in shape through the compared stages.

Conclusions

The differences observed are also reflected in the functional variation in feeding mechanism between both species. Likewise, the trajectory analysis shows that the main differences are in the magnitude of the shape change through time. Both species follow a unique trajectory, which is explained by the timing between stages.

Spatial regulation by multiple Gremlin1 enhancers provides digit development with cis-regulatory robustness and evolutionary plasticity

Precise cis-regulatory control of gene expression is essential for normal embryogenesis and tissue development. The BMP antagonist Gremlin1 (Grem1) is a key node in the signalling system that coordinately controls limb bud development. Here, we use mouse reverse genetics to identify the enhancers in the Grem1 genomic landscape and the underlying cis-regulatory logics that orchestrate the spatio-temporal Grem1 expression dynamics during limb bud development. We establish that transcript levels are controlled in an additive manner while spatial regulation requires synergistic interactions among multiple enhancers. Disrupting these interactions shows that altered spatial regulation rather than reduced Grem1 transcript levels prefigures digit fusions and loss. Two of the enhancers are evolutionary ancient and highly conserved from basal fishes to mammals. Analysing these enhancers from different species reveal the substantial spatial plasticity in Grem1 regulation in tetrapods and basal fishes, which provides insights into the fin-to-limb transition and evolutionary diversification of pentadactyl limbs.

The BMP antagonist Gremlin1 balances BMP and SHH signalling, endowing limb bud development with robustness. Here, the authors identify enhancers controlling Grem1 levels in an additive, and spatial regulation in a synergistic manner, providing digit patterning with cis-regulatory robustness and evolutionary plasticity.

New Species Can Broaden Myelin Research: Suitability of Little Skate, Leucoraja erinacea

Although myelinated nervous systems are shared among 60,000 jawed vertebrates, studies aimed at understanding myelination have focused more and more on mice and zebrafish. To obtain a broader understanding of the myelination process, we examined the little skate, Leucoraja erinacea. The reasons behind initiating studies at this time include: the desire to study a species belonging to an out group of other jawed vertebrates; using a species with embryos accessible throughout development; the availability of genome sequences; and the likelihood that mammalian antibodies recognize homologs in the chosen species. We report that the morphological features of myelination in a skate hatchling, a stage that supports complex behavioral repertoires needed for survival, are highly similar in terms of: appearances of myelinating oligodendrocytes (CNS) and Schwann cells (PNS); the way their levels of myelination conform to axon caliber; and their identity in terms of nodal and paranodal specializations. These features provide a core for further studies to determine: axon-myelinating cell communication; the structures of the proteins and lipids upon which myelinated fibers are formed; the pathways used to transport these molecules to sites of myelin assembly and maintenance; and the gene regulatory networks that control their expressions.