EvoDevo

The evolution of visual systems has compelled numerous investigations of developmental processes underlying eye patterning across Bilateria. It is well-established that homologs of the transcription factor Pax6 play a highly conserved role in eye fate specification and are at the top of the retinal determination gene network (RDGN) hierarchy. In insects, the two Pax6 homologs eyeless (ey) and twin of eyeless (toy) are required for the development of the two visual systems broadly found within the phylum (i.e., median and lateral eyes). Curiously, Pax6 homologs do not appear to maintain this function in well-studied chelicerate models, with emphasis on spiders, a lineage of arachnids with great diversity of eye form and acuity. It was recently proposed that the gene Pax2 (shaven; sv) may have subsumed the role of eye fate specification in chelicerates, a hypothesis predicated upon the observation that one of two spider Pax2 copies is strongly expressed in the developing lateral eyes during embryogenesis. However, no functional data are available for any Pax homologs across Chelicerata. We examined the incidence of Pax family genes across Chelicerata, as well as interrogated the expression and function of Pax2 and Pax6 homologs in the daddy-longlegs Phalangium opilio, an arachnid recently discovered to bear a highly plesiomorphic arrangement of visual systems. Here, we show that ey and toy are expressed early in the developing head lobes of P. opilio, whereas sv is not expressed until well after stages when downstream RDGN members (eyes absent and sine oculis) are already activated. Gene silencing of ey, toy, and sv individually had no discernible effect on eye development. By contrast, double knockdown of ey and toy resulted in an array of median eye defects, spanning loss of some cells of the eye to total loss of the median eyes. Gene expression assays also showed that depletion of the two Pax6 copies resulted in failure of the vestigial median and vestigial lateral eyes. These data are consistent with a conserved role for Pax6 homologs in patterning both visual systems and all three eye pairs in the daddy-longlegs. Our results comprise the first functional data for Pax6 genes in any chelicerate and suggest that heterochronic shifts in expression, rather than changes in function, underlie the atypical dynamics of Pax genes in derived arachnid groups such as spiders.

In tunicates, larvaceans represent a fascinating case of evolution, where the chordate body plan has been maintained despite a rapidly evolving genome characterized by strong In contrast to other tunicates, larvaceans keep the chordate body plan during their entire life. They have acquired a highly specialized epithelium in charge of producing the "house", a complex extracellular apparatus used for filter feeding in the plankton. To what extent the house and this epithelium represent true molecular innovations withing chordates is a question for which thorough transcriptomics can bring novel insights. We conducted a developmental profiling of gene expression at the single-cell level in the larvacean Oikopleura dioica. We provide detailed descriptions of cellular transcriptomes associated with the house-synthesizing organ, which permits to define the molecular specifics of epithelial cell territories. We followed their emergence during development, and we identified genes that represent key candidate molecules for regulating the morphogenesis of the house-producing organ. Dynamic changes in gene expression and cell identities during major developmental transitions of the lifecycle illustrate that our dataset effectively allows access to the diversity of O. dioicas cell types in embryos and in adults. The resources presented here constitute critical assets to investigate larvacean biology and evolution for mechanistic and comparative goals.

Animals sense and integrate complex external cues to make developmental decisions that help them better survive and adapt to their natural habitats. Under environmental adversity, nematodes can enter an alternative developmental pathway to form a diapautic and stress-resistant stage, termed the dauer larvae. While dauer formation has been well characterized in Caenorhabditis elegans, how environmental factors influence analogous stages in other nematode species remains largely unexplored. This study examines how symbiotic bacteria, temperature, and pheromones affect the formation of the infective juvenile (IJ), a dauer-like stage, of the insect-parasitic nematode Steinernema hermaphroditum. In contrast to C. elegans, where dauer entry is promoted by heat, IJ development in S. hermaphroditum development is enhanced by reduced temperature. Moreover, the presence and absence of live symbiotic bacterium Xenorhabdus griffiniae functions as an ON-and-OFF switch that regulates the host IJ formation. Crude pheromone extracts from S. hermaphroditum liquid culture do not robustly induce IJ formation in a dose-responsive manner, unlike the potent pheromone-driven dauer entry observed in C. elegans. Nutrient-rich liver-kidney media that mimics host insect environment showed IJ entry induction in a pheromone-dependent manner. These data suggest that external cues, such as temperature, microbial diet, and pheromone, are perceived differently by S. hermaphroditum in comparison to that of C. elegans, reflecting species-specific adaptations to distinct ecological niches and life history strategies.

To achieve a stereotypic lineage, each embryo of Caenorhabditis elegans follows an invariant cell differentiation process arising from a combination of cell polarisation, asymmetric or symmetric divisions, combined with intercellular signalling processes. This pattern of embryonic cell differentiation is driven by regulated segregation of molecules occurring at each cell division, including polarity proteins or cell fate determinants, transcription factors, p-granules and mRNAs. These distribution patterns are coupled with a robust spatio-temporal orchestration of cortical actin dynamics, which also plays a crucial role in these processes. However, compared to other molecular contents, how the actin per se is segregated from the first asymmetric division onward remains poorly understood. This study presents a thorough quantification of the intracellular distribution from the zygote to the 4-cell stage of key actors related to actin polymerisation: two nucleators (a formin and the Arp2/3 complex), a capping protein and E-cadherin. We additionally developed a novel method to assess actin polymerisation capacities from single blastomere extracts. We found that actin-related signatures arise at these early stages and that differential mechanisms of protein segregation and homeostasis occur, depending both on the cell pair and on the protein considered. Notably, if asymmetric divisions correlated with unequal partitioning of actin-related contents in a process linked with embryonic polarity, differences were revealed between AB daughter cells upon their separation. Taken together, these actin-related asymmetric distributions are adding a layer to the complexity of cell fate acquisition mechanisms in the early embryo.

Dauer larvae are a dormant developmental stage in nematodes that is induced by a range of environmental cues. The molecular mechanisms that transduce these cues to regulate dauer entry have been well characterized in Caenorhabditis elegans, whereas those in other nematode species remain unclear. The closest known sibling species of C. elegans, Caenorhabditis inopinata, occupies a distinct ecological niche and shows an extremely low frequency of dauer formation by starvation in laboratory conditions, suggesting that it could serve as a useful comparative model for analyzing dauer-inducing mechanisms. To support such analysis, we generated a fluorescent dauer reporter, Cin-col-183p::mCherry, in C. inopinata based on a previously reported dauer-specific reporter in C. elegans. This reporter showed fluorescence specifically in the pre-dauer and dauer stages, but not in other developmental stages, indicating that it functions as a dauer-specific marker in C. inopinata. Using these marker strains, we compared the responses to high temperature and RNAi-mediated knockdown of insulin/IGF-1 pathway genes (daf-2, age-1, and pdk-1), and found that dauer induction differs mechanistically between C. elegans and C. inopinata. This dauer-specific fluorescent strain will be a useful tool for investigating the diversity of dauer-inducing mechanisms across nematode species. Article SummaryDauer is a dormant developmental stage in nematodes induced by environmental stress. Although its regulation is well studied in Caenorhabditis elegans, the mechanisms in other species remain unclear. Here, we developed a fluorescent dauer reporter, Cin-col-183p::mCherry, in Caenorhabditis inopinata, a close relative of C. elegans. The reporter was specifically expressed in pre-dauer and dauer stages, confirming its usefulness as a dauer marker. Using this strain, we found that responses to high temperature and insulin/IGF-1 pathway gene knockdown differ between C. elegans and C. inopinata. This reporter will help reveal diversity in dauer-inducing mechanisms across nematode species.

The cockroach Blattella germanica possesses panoistic ovaries, in which oocytes lack nurse cells and therefore need to rely on their own transcriptional activity to support embryogenesis. Ovarian development in this species involves the development of a single basal ovarian follicle (BOF) per gonadotropic cycle, a process strictly regulated by endocrine signals, primarily juvenile hormone and ecdysone, which act at both the transcriptional and translational levels. In addition, transcriptional activity in these ovaries is necessary for both regulating and genome protection, and at this level, PIWI-interacting RNAs (piRNAs) play an essential role. Although insect ovaries are known to be particularly rich in piRNAs, their function in ovary maturation is still not well defined. For this purpose, we characterize the piRNA expression dynamics across seven key developmental and reproductive stages, ranging from late nymphal instars to post-vitellogenic adults. piRNA expression in B. germanica shows coordinated fluctuations. Expression remains stable in previtellogenic ovaries, whereas vitellogenic ovaries show pronounced changes. Moreover, vitellogenic ovaries exhibit reduced piRNA diversity due to strong enrichment of a subset of highly expressed piRNAs. Our data show that although piRNAs predominantly map to transposable elements, particularly LINEs, there is a notable increase in gene-derived piRNAs toward the end of the cycle. Our results suggest regulatory roles of piRNAs in modulating both TEs and mRNAs during BOF maturation, likely related to changes in the follicular cell program.

Sex-responsive trait development generates much of the phenotypic variation found in natural populations and diversifies rapidly among closely-related taxa. Furthermore, rather than exhibiting equal sexual dimorphism across all traits, organisms are mosaics of tissues that vary in their degree of dimorphism. Yet, how these mosaic patterns are generated remains largely an open question, as sexually dimorphic traits have typically been studied individually in select model systems. In this study, we compare gene regulatory landscapes across five traits that differ in the degree of morphological sexual dimorphism in the bull-headed dung beetle Onthophagus taurus by assaying tissue-specific gene expression and chromatin accessibility at the onset of pupal development when future adult form is specified. We identify a modest number of pleiotropic regulators associated with sex differences across traits, yet uncover a high degree of sex- and trait-specificity in chromatin architecture within developing tissues. We then confirm the role of the sex determination factor doublesex in the regulation of sex differences through expression of sex-specific isoforms, and uncover trait- and sex-specific sets of Doublesex binding sites likely underpinning context specific sexual dimorphisms. Further, we identify and functionally validate the transcription factor ventral veinless as a regulator of sexually dimorphic development. Our findings suggest that in contrast to doublesex, ventral veinless does not exhibit sex-biased expression, yet exerts its sex-specific regulation via sets of differentially accessible binding sites. This work furthers our understanding of the molecular mechanisms instructing the development of sex differences and provides novel insights illustrating how transcriptional activity and chromatin remodeling interact to generate sexual dimorphism in a trait-specific manner. More generally, our work contributes to a growing body of knowledge on how development integrates cues such as sex determination to enable highly similar genomes to yield diverse phenotypic outcomes.

Medusozoan cnidarians (e.g., jellyfish) metamorphose from a benthic juvenile polyp into a pelagic adult medusa, providing a well-known example of a clade that uses tissue remodeling to create distinct juvenile and adult body plans. Staurozoans (i.e., stalked jellyfish) are an atypical lineage of medusozoans that have lost their medusa stage; thus, their juvenile and adult body plans look remarkably alike. Their limited metamorphosis is characterized by the regression of primary (juvenile) tentacles and the development of secondary (adult) tentacles. In some staurozoan lineages, metamorphosis also involves development of novel adhesive structures (anchors), which are built on top of the regressing primary tentacles. Understanding how cells are partitioned from making juvenile tissues to making adult tissues is important for understanding how animals can make adult structures in the absence of complete metamorphosis. We compared the abundance and distribution of proliferative cells in tissues undergoing regression (primary tentacles) and development (secondary tentacles and anchors) during the juvenile to adult transition in the San Juan Island stalked jellyfish, Haliclystus sanjuanensis. We show that proliferative cells are lost in regressing primary tentacles but are gained in anchors, consistent with a shift in investment from juvenile to adult tissue. Prior to regression, primary and secondary tentacles show similar patterns in their proliferative cell distribution and in the identity of their cnidocytes (stinging cells), indicating that adult tentacles are made by re-deploying a juvenile tentacle program. Finally, we demonstrate that unlike secondary tentacles, primary tentacles cannot regenerate, illustrating that the temporary investment in this tissue is tied to their loss of proliferative cells. Thus, we propose that continued investment in a population of proliferating cells is an important mechanism for segregating temporary tissues (primary tentacles) from long-term tissues (secondary tentacles). These observations of cell dynamics in H. sanjuanensis suggest that temporary investment into juvenile structures may be used to pattern novel adult tissues, providing an important mechanism for diversifying adult body plans.

Sea star wasting disease has caused widespread mortality in the kelp forest predator, the sunflower sea star (Pycnopodia helianthoides). Wild populations have declined by up to 99% in parts of their native range along the western North American coast. In response, a multi-institutional conservation breeding and rearing program has been initiated to support future reintroduction efforts for the species. We split a full-sibling cohort across four larval density treatments (1 larva/ml, 2 larvae/ml, 5 larvae/ml, and 15-20 larvae/ml) to assess the effects on larval settlement, juvenile survival, and juvenile fitness at 12 months old. Stars raised in the highest density treatment displayed a lower settlement rate and were significantly smaller than the other density groups at 12 months old, but showed no significant difference in flip time, a measure of fitness. Additionally, measurements of diameter, weight, and arm count across modern and historical juvenile and adult stars indicate that P. helianthoides experience exponential weight gain as they grow in length, with corresponding asymptotic growth in arm count. These findings will inform best practices for the aquarium propagation of P. helianthoides and will contribute to broader efforts aimed at reestablishing populations in the wild.

Stony corals are the only cnidarians to secrete a robust calcium carbonate skeleton making them the critical keystone species for one of the most diverse marine ecosystems on earth: coral reefs. Despite their importance, very little is known about the genetic interactions that drive the earliest life history stages of corals during which coral-specific traits emerge, such as specification of coral skeletogenic cells. Here, we used a combination of chromatin profiling and gene expression assays to derive the cis-regulatory architectures of early coral development using the emerging model system Astrangia poculata. From this work we found that elements of the cnidarian endomesoderm gene regulatory network have been co-opted into a subnetwork underlying the specification of coral skeleton secreting cells. We further identified the cis-regulatory element responsible for the novel expression of the endomesodermal gene Brachyury in coral skeletogenic cells and demonstrated that this element is capable of reproducing similar expression patterns in Nematostella vectensis, a distantly related species that does not produce a skeleton. These findings support a novel dual role of the endomesodermal GRN in establishing germ layer identity and specifying skeleton secreting cells in stony corals and provide a gene-regulatory framework that underpins the evolution and diversification of stony corals from other cnidarian lineages.

Understanding the reproductive biology of non-model organisms remains challenging due to the limited availability of high-resolution molecular resources. Here, we present a comprehensive single-cell transcriptomic atlas of the adult ovary of Spodoptera litura (S. litura), a highly polyphagous agricultural pest with a polytrophic meroistuc ovary. By integrating single-cell RNA sequencing with cross-species comparison to Drosophila melanogaster (D. melanogaster), we define major germline and somatic cell populations and delineate conserved and species-specific features of ovarian cell composition. To enhance the interpretability and reuse of this dataset, we combine transcriptomic profiling with in situ hybridization to validate cluster-specific molecular markers across ovarian cell types. We further apply RNA interference (RNAi) to assess the contributions of germline-enriched genes (Hsc70-4, Wech, Polo, Path) to ovarian development and fecundity. Trajectory inference, together with SCENIC and CellChat analyses, provides a system-level view of transcriptional regulatory programs and predicted intercellular communication pathways during oogenesis in S. litura. Collectively, this work establishes a valuable resource for studying lepidopteran insect oogenesis, offering a comparative framework for reproductive biology in non-model insects and highlighting potential targets for RNAi-based pest control strategies.

Colonial animals composed of clonally produced units can achieve a high degree of functional integration, challenging the distinction between an individual and a higher-order organism. Siphonophores (Cnidaria: Hydrozoa) exemplify this condition, forming highly organized colonies in which genetically identical zooids are specialized for functions such as locomotion, feeding, and reproduction, and are precisely arranged along a shared stem. All zooids arise from two spatially separated budding zones, the nectosomal and siphosomal growth zones, suggesting that positional information along the stem patterns colony organization at the level of the colony rather than individual zooids. However, the molecular basis of this colony-level axial patterning remains poorly understood. Here, we analyze gene expression along the stem of the siphonophore Agalma okenii using RNA sequencing and in situ hybridization chain reaction (HCR). We show that conserved developmental regulators, including Hox and Wnt pathway genes, exhibit region-specific expression corresponding to distinct budding zones and zooid distributions. These results indicate that canonical axial patterning systems are deployed at the level of the colony axis. Our findings demonstrate that developmental gene networks classically associated with anterior-posterior patterning can operate at a higher level of biological organization, providing a mechanistic framework for the evolution of integrated, superorganism-like body plans in colonial animals.

The temperate rainforests of the Pacific Northwest of North America harbor many endemic taxa whose evolutionary histories have been shaped by major climatic and geologic events. The enigmatic taildropper slugs (genus Prophysaon) are one example, notable for their ability to autonomize their tails to escape predators. Despite extensive work uncovering the evolutionary history of individual lineages, relationships among the nine recognized species of Prophysaon remain poorly understood due to insufficient molecular data. To address this, we collected transcriptomes for six of the nine currently accepted species of Prophysaon. Using these data, we were able to resolve species relationships, calling into question the existing subgeneric classification based on morphology. We also detected undescribed phenotypic diversity within the P. andersonii--P. foliolatum species complex, with molecular data supporting the distinctness of two phenotypically distinct populations from Washington. Finally, our transcriptomic data suggest a moderate role of introgression in shaping the evolutionary history of Prophysaon. Here, we synonymize the subgenus Mimetarion with nominotypical Prophysaon. Future work should further investigate whether the undescribed diversity detected here represents species level differentiation.

Predator-prey interactions are key drivers of behavioural and life-history evolution, yet their mechanisms remain difficult to study in natural contexts. The nematode Pristionchus pacificus is a model predator, but most studies exploring its behaviours use Caenorhabditis elegans as prey, a species that it likely only rarely encountered in nature. Here, we examine predation within nematode communities associated with beetle carcasses, the native necromenic habitat of P. pacificus. We identify Oscheius myriophilus as a cohabiting species, likely representing natural prey. Using predatory assays, automated tracking, and machine-learning-based behavioural analysis, we show that P. pacificus actively kills and consumes O. myriophilus. Strikingly, predation rates are lower than those observed for C. elegans, suggesting partial resistance or reciprocal adaptation in O. myriophilus. Consistent with this, O. myriophilus exhibits a mixed reproductive strategy, with early oviposition followed by ovoviviparity and matricide. As later developmental stages are more resistant to predation, internal hatching may protect offspring while providing maternal resources for development. These findings establish these nematodes as a tractable model for investigating predator-prey interactions and their evolutionary consequences, highlighting how behavioural strategies and life-history traits can co-evolve in natural communities.

Transitions from sexual to asexual reproduction are well-documented across different taxa. However, despite extensive efforts, the regulatory changes underlying the emergence of asexuality remain largely undiscovered in the majority of species studied. Artemia brine shrimp have multiple closely related sexual and obligate parthenogenetic lineages, making them a promising model for addressing this question. While earlier work suggested that asexuals use a modified meiosis, and inferred a likely role for the Z-chromosome in its transmission, no master regulator or genetic changes have been put forward as the root causes for the shift. Here, we generate single-nucleus RNAseq data of the female reproductive system of individuals from the Aibi lake population of Artemia parthenogenetica and its closely related obligate sexual species Artemia sp. Kazakhstan. We identify the germline cell clusters in the female reproductive system and perform differential expression analysis to infer substantial transcriptional differences at genes putatively involved in cell cycle and oocyte development between the meiotic cells of the two species. Additionally, we use whole-genome sequencing of 32 individuals from two backcrossing experiments to narrow down the genomic regions associated with the transmission of asexuality to an 8 megabase region of the Z chromosome. Within the identified regions, two adjacent genes with known functions in oogenesis, ITPR and USP8, show differential expression and genetic differentiation between sexuals and asexuals, making them promising candidate drivers of asexuality in this species. Significance statementWhile most animals reproduce sexually, many do not, and why and how these shifts occur remains an open question. This paper presents a systematic investigation of the molecular changes that underlie the transition from sexual to asexual reproduction in brine shrimp. We combine multiple computational and experimental approaches to look for differences between close sexual and asexual lineages. We find that a subset of meiotic germ cells is regulated differently in the two, and that two important oogenesis genes are the likely drivers of asexuality. This work is unique in providing an in-depth characterization of the combined genetic and regulatory changes underlying this key transition in reproductive modes.

Ancient human dental calculus is one of the richest archives of archaeological biomolecular information, providing direct evidence of diet, oral health, and the oral microbiome. Proteomic analyses of this biological matrix have so far focused mainly on oral microbes and dietary proteins, with milk proteins such as beta-lactoglobulin (BLG) providing the largest corpus of proteomic evidence. Despite the close relation between the various stages of dental calculus formation and mineralization with the dental enamel surface, proteins from the dental enamel matrix have not previously been reported outside of dental enamel tissue. Here we reanalysed 498 ancient dental calculus proteomes from 14 published studies (n=434 individuals) reporting the presence of BLG, spanning from the Neolithic to the Victorian Era and applying different protein extraction protocols (FASP, GASP, SP3 and in-solution digestion). Dental enamel matrix proteins were identified in ten studies (n=37 individuals), with amelogenin being the most frequently detected. Enamel peptides occurred more often in studies that applied SP3, although amelogenin was successfully identified through both SP3 and FASP. Structural proteins, including enamelin, ameloblastin, and MMP20, were also identified. The detection of AMELX and AMELY peptide sequences provided new insights into cases where the sex was previously undetermined. These findings establish dental enamel proteins as a new category of biomolecules detected in dental calculus, broadening its application beyond diet and microbiome studies to possible sex estimation. HighlightsO_LIDental calculus entraps oral microbes along with endogenous and exogenous particles during formation and mineralization C_LIO_LIWe conduct reanalysis of 14 published ancient dental calculus studies (n = 434 individuals) spanning the Neolithic to Victorian Era C_LIO_LIDental enamel proteins AMELX, AMELY, AMBN, COL17A1, ENAM and MMP20 are identified in ancient human dental calculus C_LIO_LIAmelogenin was the most frequently detected enamel protein C_LIO_LIWe expand dental calculus palaeoproteomics beyond diet and oral microbiome to potentially include sex estimation C_LI

Early neural development involves a combination of planar signals from the vertebrate organiser and vertical signals from its derived structures, the prechordal plate and notochord. However, the relative contribution of each structure to neural development is not clear. Here, we isolate anterior tissues from the primitive streak at successively later stages of development, to identify the extent of patterning that can occur prior to, during, and after the formation of the organiser and its later derivatives. Our results show that acquisition of neural identity occurs gradually and that exposure to planar signals from the developing node is necessary for neural plate specification. We also show that planar node-derived signals are required for AP patterning in isolated anterior tissues and give evidence that early neural tissue is of anterior character which subsequently becomes caudalised by signals (in part) from the developing node. However, anterior neural identity is lost without long-term contact with vertical signals from the axial mesendoderm. These results reveal a previously unappreciated level of autonomy in anterior neural development in the absence of node derived tissues. Summary statementCulture of isolated anterior tissues from the chick embryo reveal the roles of planar and vertical organiser signals for neural specification and anteroposterior patterning and maintenance.

Identifying spatial and temporal patterns of connectivity among populations is fundamental to marine ecology, evolutionary biology, and fisheries management. Yet, due to large population sizes and low genetic differentiation among populations, empirical quantification of population connectivity across a species entire range has not been achieved for an open-coast marine organism. Here, we leverage experimental transcriptomics to develop a genotyping-in-thousands by sequencing (GT-seq) panel to support assignment of recruits of the kelp forest gastropod, Kellets whelk (Kelletia kelletii), collected across the species biogeographic range. Over a three-year period, we identified high self-recruitment in the historical range (100%) and low self-recruitment in the expanded range (10.53 - 13.73%). Additionally, self-recruitment within the expanded range generally increased with recruit age, from 27.14% at 0.93 years to 43.40% at 1.93 years, indicating that the locally spawned individuals were more likely to survive to older ages than migrants from the historical range. Together, these results reveal limited self-recruitment in the expanded range and suggest that a post-settlement selective filter contributes to differential survival in a high gene flow marine system.

Dental enamel, the hardest mineralised tissue in the human body, has proven to be an excellent source of ancient proteins, which have been found to survive within dental enamel for at least twenty million years. In archaeological and palaeontological contexts, the enamel proteome is generally considered to be rather small, consisting of about twelve proteins, most of which are unique to enamel. During amelogenesis these proteins undergo in vivo digestion by matrix metalloproteinase 20 (MMP20) and kallikrein 4 (KLK4) as well as serine phosphorylation by family with sequence similarity member 20-C (FAM20C) that alter their characteristics. Gaining knowledge of the previously understudied influence of amelogenesis on the archaeological human dental enamel proteome could benefit various palaeoproteomic analysis, especially in an human evolutionary context. Here we present archaeological dental enamel proteomes and explore protein cleavage patterns and sequence coverage to estimate the effects of in vivo digestion, as well as explore phosphorylation patterns. Additionally, we present a new marker based on phosphorylation to estimate genetic sex.

Phenotypic diversity is often thought to arise from the evolutionary modification of developmental processes. However, developmental processes are tightly coupled in space and time, with each process beginning from conditions set by the one before it. While we know from dynamical systems theory that initial conditions can significantly affect a systems out-come, their importance as a source of phenotypic evolvability has been largely overlooked. Here we show for the first time, that phenotypic evolution can proceed through changes in developmental initial conditions while the underlying developmental process remains conserved. Somitogenesis is the process by which vertebral precursors, known as somites, are periodically patterned in the pre-somitic mesoderm (PSM). Somitic count (total number of somites) is thought to diversify through the evolution of components of somitogenesis such as the tempo of the segmentation clock or the mechanisms driving axial morphogenesis. Using two closely related species of Lake Malawi cichlid fishes that differ in vertebral counts, we show that somite count evolution has happened without changes to somitogenesis itself, but instead, by altering the size of the PSM at the onset of this process. This work will expand what we consider developmental drivers of phenotypic evolution and highlight the importance of comparative studies to understand the diversification of phenotypes.