EvoDevo

Silk glands have been found in two groups of amphipods: the Corophiida and the Ampeliscidae. The silk glands in Ampeliscidae, however, have yet to be examined in detail. Here we report, for the first time, the morphology and distribution of pereopodal glands in the Ampeliscidae, in non-thread producing Synopiidae, and in the Paragammaropsidae. In the Ampeliscidae we found two gland types distributed throughout all pereopods which have the ability to create threads. Pereopods three and four have additional silk extrusion morphology at the tip of the dactylus in which silk is transformed into semi-cylindrical threads used for building domiciles. Synopiid outgroup species have one of the gland types but lack silk extrusion morphology. Using ancestral state reconstruction analysis, we find that glands in the Synopiidae are likely ancestral and hypothesize that silk glands in Ampeliscidae are derived from these ancestral glands. Silk-spinning pereopods in the Paragammaropsidae had similarities with both Corophiida and Ampeliscidae but had distinctions. Ampeliscidae silk-spinning systems bear surprising resemblance to the Corophiida which presents one to reconsider the taxonomic placement of Ampeliscidae and the origins of silk-spinning in amphipods. This is the first comprehensive study on the glandular systems of Ampeliscidae, Synopiidae, and Paragammaropsidae using advanced microscopy, providing pertinent morphological data to the study of arthropod silk gland evolution and complex traits.

Acropora spp. are the dominant reef-builders of the Indo-Pacific but are also amongst the most stress-sensitive corals. For these reasons, Acropora spp. have become the most studied of corals, two species (A. digitifera and A. millepora) often essentially serving as the basis for understanding molecular responses and processes across the sub-order Refertina and corals in general. The early development of these species has been well-characterised in terms of morphology and gene expression but as yet we have a limited understanding of how transcription is regulated during development. In "higher" animals (bilaterians) microRNAs (miRNAs) are critical regulators of gene expression but until now their involvement in coral development has not been investigated. Building on the existing developmental data for Acropora spp., we catalogued microRNAs (miRNAs) expressed during the early development of Acropora digitifera and profiled their expression in 21 stages from unfertilised eggs to 24h after treatment with a natural settlement cue (CCA chips). 157 miRNAs were recognised, many of which ([~]60%) were novel. These fell into three distinct groups, corresponding to three distinct developmental phases: (1) those present in eggs through to gastrulation (2) a larvally expressed group and (3) those expressed following settlement induction. Exposure of competent larvae to a natural settlement inducer resulted in major changes in the miRNA profile within 10 minutes, indicating that miRNAs may be particularly important in mediating the larva/polyp transition but are also likely to play important regulatory roles throughout early coral development in addition to possible roles in disease resistance.

In many oviparous animals, egg yolk is the sole source of nutrition until feeding begins, and carbohydrates are present in only small amounts in the yolk. Glucose plays an important role in the developmental processes of various animals. In addition, gluconeogenesis has been reported to occur in the yolk syncytial layer (YSL) of cartilaginous fish and teleosts. In contrast, the role of gluconeogenesis in tetrapods remains unclear. In this study, we used Xenopus tropicalis, an anuran amphibian, which lacks YSL, and therefore provide an opportunity to examine the evolutionary conservation of gluconeogenic mechanisms among vertebrates. In X. tropicalis, liquid chromatography/mass spectrometry revealed that glucose levels increased before liver formation. Subsequent tracer experiments using 13C-labeled metabolic substrates detected gluconeogenesis activity from glycerol and lactate. Expression analyses showed that gluconeogenic genes are expressed in the epidermis and endoderm. Consistently, G0 knockout of fbp1, a key gluconeogenic gene, resulted in a significant reduction in glucose levels, affecting brain development. These findings first demonstrate that gluconeogenesis supports development of X. tropicalis. To the best of our knowledge, gluconeogenesis in developing epidermis has not been reported, highlighting previously unrecognized diversity in tissue-specific metabolism during vertebrate development. Comparative analyses across species will provide further insights into the evolution and functional significance of embryonic gluconeogenesis and nutrient metabolism.

The wings and legs of insects are both appendages that develop along a proximal-distal (PD) axis and likely share many underlying patterning mechanisms. Comparisons between the two appendages have been detailed in Drosophila melanogaster, a derived insect where both traits develop initially as imaginal discs. Here, we visualise and compare the expression of prominent PD developmental genes in the embryonic legs and larval wings of a lepidopteran species, Bicyclus anynana. We examine the domains of twelve leg gap genes that subdivide legs into distal, medial, or proximal domains, three morphogens, and two genes that refine the PD axis after initial specification. Our results reveal high spatial congruence in the order of PD gene expression between the two appendages. Notably, we observe a distinct loss of the medial domain in the Bicyclus wing compared to the leg, providing evidence for the evolutionary re-patterning of these structures. Comparisons with Drosophila further highlight conserved versus lineage-specific regulatory architectures. These findings suggest a deeply conserved PD patterning logic across Holometabola, while pointing to divergent mechanisms that likely facilitated the morphological innovation of butterfly wings.

The morphogenesis of complex vertebrate appendages requires precise regulation of growth, governed by distinct positional identities. The zebrafish caudal fin achieves a symmetrical, forked morphology through the regional specialization of the bony rays: peripheral rays are composed of relatively long, thick segments; while the central rays are made up of shorter, thinner segments, and their overall length is restricted. This length differential establishes the definitive forked shape of the organ. We asked whether these regional morphological differences reflect distinct underlying positional identities. Transcriptomic profiling of intact tissues from adult wild-type zebrafish suggested that central rays possess unique expression profiles, distinct from those of peripheral rays. We previously identified a treatment during embryogenesis that allows excess growth in the central rays, creating a truncate fin shape in adults-we asked whether this novel fin shape was caused by a peripheralization of the central rays. Indeed, the central rays of truncate fins were not only longer, but were composed of longer and thicker individual segments, reminiscent of peripheral rays. Further, gene expression in the central regions of truncate backgrounds showed signatures of peripheral identity. During development of the truncate phenotype, peripheral markers became expressed in more central domains of the growing truncate caudal fin, and in the supportive endoskeleton, the central hypural diastema was lost from the earliest stages. Ultimately, our results demonstrate how adult morphologies may be altered by shifts in positional identities. These findings clarify the anatomical patterning and molecular profiles that underlie regional specialization during caudal fin development.

Calpains constitute an ancient, extensive family of calcium-dependent cysteine proteases found in some bacteria and most eukaryotes. They are involved in a wide variety of developmental and cellular processes and are implicated in major human diseases, yet it remains to be seen if they have a common core function explaining their widespread and varied presence across taxa. Beyond their core CysPc catalytic domain, calpains contain diverse domain combinations and can be either cytosolic or membrane bound. Here we hypothesize a general role for both cytosolic and transmembrane calpains in cellular cytokinesis through positional anchoring and organization of microtubules (MTs). We propose that during plant cell division, the singular transmembrane calpain DEK1 localizes and organizes the array of cortical MTs from the microtubule organizing center (MTOC) to establish the location of the preprophase band and/or the site of cell plate formation according to the positional activation of DEK1 proteins in the nuclear membrane. Similarly, during cell division in animals, their calpains may be involved in setting the point of membrane invagination via their association with membrane-bound proteins. This proposition adds to the current picture of animal MTOC/centrosome function and suggests how a calcium peak during the initial cytokinetic furrowing might be transmitted. We discuss this novel mechanistic model for calpain activity in the context of data from the animal and plant literature, as well as of our novel discovery here of calpain sequences in both brown and red algal genomes. Finally, we speculate that the ancestral role of calpains in early eukaryotes, before the split into the major eukaryotic supergroups, may have been to facilitate the formation and function of MT arrays in flagella and cilia. From this origin, calpains may have developed new functions in eukaryote cell division processes by anchoring centrosomes/MTOC to set the cell division orientations that are especially important for complex multicellularity.

How morphology forms during development and changes during evolution remain major questions in biology. In vertebrates, teeth have long served as model systems to address these questions. In threespine stickleback fish (Gasterosteus aculeatus), repeated and convergent increases in pharyngeal tooth number in derived freshwater sticklebacks occur, suggesting increased tooth number is adaptive in freshwater environments, likely due to a diet of larger prey in freshwater. Whether changes in oral tooth patterning also occur in freshwater sticklebacks was unknown. Here we describe oral tooth number and patterning in a dense developmental time course of lab-reared ancestral marine and derived freshwater fish. We address three major questions. First, is the spatial sequence of early oral tooth formation invariant as we previously described for the pharyngeal dentition? Second, is oral tooth patterning in the upper and lower jaw sexually dimorphic, and if so, when during development does this dimorphism arise? Third, have freshwater fish evolved increases in oral tooth number? We find that (1) unlike the pharyngeal dentition, the oral jaw early spatial sequence is variable, especially in the lower jaw (2) sexual dimorphism in both oral jaws arises at the late juvenile stage with males having more teeth and (3) freshwater fish have evolved more oral teeth similar to the evolved tooth gain in the pharyngeal jaw. Together our morphological descriptions advance the stickleback oral jaw as a model system to study how morphology forms during development and evolves in nature.

How gene expression patterns change spatially as the embryo transitions from simple to complex structures remains a major developmental biology question. Recently developed imaging-based spatial transcriptomics (ST) enable mapping expression of multiple gene at a single-cell resolution. Although Xenopus is a key model in embryology there is no established ST pipeline, and commercially available techniques face many challenges (sample preparation, probe design, cell segmentation). Furthermore, the highly diverse cell shapes and sizes across developmental stages and between different tissues represent major hurdles to accurately defining cells. Here, we describe an optimized workflow for ST in blastula-to-tailbud-stage frog embryos using Merscope, commercial MERFISH (Multiplexed Error-Robust Fluorescence In Situ Hybridization) originally designed for standard mammalian tissues. With stringent quality control and tailored computational pipelines, we optimize this technology for robust, semi-quantitative profiling of spatial transcriptomic landscapes in non-mammalian embryos. Reliable tissue preservation and cell-segmentation enable high-resolution mapping of gene expression during the development of a complex multi-tissue organization. This versatile strategy applies broadly to various dynamic systems, from embryos of various model organisms to complex and heterogeneous organs in mammals. Summary statementThis Single-cell Spatial Transcriptomics pipeline and reference atlas in Xenopus - a model organism in embryology - overcome technical challenges and resolve dynamic changes in patterning during development.

AO_SCPLOWBSTRACTC_SCPLOWEnvironmental heterogeneity across freshwater systems often promotes phenotypic variation, yet disentangling environmentally induced variation from heritable differentiation remains a central goal in evolutionary ecology. We investigated the geographic distribution and morphological differentiation, and heritability of shell traits among populations of the freshwater lymnaeid snail Pectinidens diaphanus in Patagonia. Extensive field surveys across 196 freshwater sites revealed that the species occupies a broad range of lentic and lotic habitats and constitutes the only lymnaeid inhabiting southern Patagonia. While reproductive anatomical structures were conserved across populations, shell shape differed markedly among populations from contrasting habitat types, with population identity explaining nearly 50% of total shape variation. Populations from hydrologically unstable habitats (ponds and streams) exhibited more elongated shells and relatively smaller apertures, a pattern consistent with functional responses to hydroperiod variability and desiccation risk. To assess the heritability of this differentiation, we conducted a common-garden experiment across two generations. Shell shape differences between permanent- (lagoon) and temporary- (pond) habitat-derived populations persisted into the G2 generation reared under standardized laboratory conditions, indicating that the observed variation is not solely a response to local environmental conditions but includes a heritable component. Together, our findings demonstrate that P. diaphanus constitutes the sole lymnaeid across southern Patagonia, occupying a broader range than previously documented, and that populations show heritable shell differentiation potentially associated with contrasting freshwater habitats. By integrating large-scale biogeographic surveys with morphometric and experimental approaches, this study provides new insight into how habitat variation may contribute to ecological and evolutionary differentiation in freshwater gastropods.

Arthropods exhibit an exceptional diversity of life histories, where developmental modes involve moulting stage progressions with changes ranging from the bare minimal to the dramatically transformative. While this variability drives many research questions aiming to understand evolutionary and developmental underpinnings of life history differences, it can complicate comparative analyses across taxa. However, this can be approached by applying a framework that defines metamorphosis as a post-embryonic stage progression characterised by substantial changes in morphology and adaptive landscape. Employing this framework with a phylogenomic dataset spanning 26 orders and encompassing four independently arising metamorphic lineages, we explore gene repertoire evolutionary dynamics potentially associated with metamorphosis in Pancrustacea. The approach contrasts gene family evolutionary dynamics inferred to have occurred in the last common ancestors of the metamorphic Insecta, Copepoda, Eucarida, and Thecostraca, with those of their sister lineages, as well as of descendent and ancestral nodes. The results reveal that the metamorphosis ancestors are characterised by an elevated number of gene family births and expansions. Expanded gene families share a set of commonly enriched biological processes across all metamorphosis ancestors, suggesting functional convergence by independent evolution of distinct gene families involved in embryonic and post-embryonic development and nervous system differentiation. Evolutionary modelling further highlights a subset of these families exhibiting signatures of adaptive, lineage-specific gene family size increases associated with metamorphic development. These families include genes implicated in neural and sensory development, segmentation, and moulting. These findings support a model of the evolution of pancrustacean metamorphosis where distinct gene families from a common functional toolkit expand and are co-opted into facilitating transitions to multi-phasic life cycles. This reframes the role of moulting in arthropod diversification to be recognised as an important reservoir of genetic change that can potentiate truly remarkable life history transitions.

The presence of multiple discrete color patterns within a species has captivated evolutionary biologists for more than a century, especially when such polymorphism is confined to one sex. The brown anole Anolis sagrei exhibits a female-limited polymorphism in dorsal patterning, which is controlled by allelic variation at the autosomal gene CCDC170. Here, we present and test a threshold model that can explain why the polymorphism is female-limited. We hypothesize that allelic variation at the CCDC170 locus affects only female color pattern because this gene is co-expressed with its neighboring gene ESR1, highly expressed in female, but not male, embryos. By manipulating embryonic estradiol levels, we show that genetic males can be induced to express the polymorphism according to allelic variation at the CCDC170 locus, which is naturally masked by low expression levels of this gene. Inversely, treating genetic females with fadrozole, which depletes estradiol, leads to monomorphic patterns irrespective of genotype, as for natural males. Using RT-qPCR, we demonstrate that these effects are accompanied by a direct influence of estradiol and fadrozole on gene expression levels of CCDC170 and ESR1, thereby validating the threshold model. Our results suggest that the CCDC170-ESR1-locus is part of a mechanistic link between the morph-determining and the sex differentiation systems and provide a causal explanation for the developmental origin of a sex-limited color polymorphism.

The diversity of body shapes is one of the most prominent features of phenotypic variation in mammals. Yet, mammalian body shapes are poorly quantified and the underlying components contributing to its diversity as well as its relationship to other components of the skeleton are rarely tested. Here, we use lagomorphs (hares, rabbits and pikas) as a model system to (1) investigate which components of the skeleton contributed the most to body shape diversity, (2) examine the relationships between body shape and relative limb lengths, and (3) test how body size, ecotype, burrowing behavior, and locomotor mode influenced variation in lagomorph body shape and appendicular morphology. We quantified the body shape and functional proxies of the appendicular skeleton in 40 lagomorph species from osteological specimens held at museum collections. Using phylogenetic comparative methods, we found the relative length of the ribs and elongation or shortening of the thoracic and lumbar regions contributed the most to body shape evolution across lagomorphs. Second, we found that only leporids (hares and rabbits) exhibited a significant relationship between limb length and body shape, where more elongate species exhibit relatively shorter forelimbs and hindlimbs. Lastly, we found that models incorporating body size were the best predictors of lagomorph body shape and the majority of the appendicular traits, whereas models incorporating burrowing behavior and locomotor mode were largely poor fits. Broadly, these results indicate that larger lagomorphs tend to exhibit more robust body shapes with longer, more gracile forelimbs, whereas smaller lagomorphs tend to exhibit more elongate body shapes with shorter, more robust forelimbs. Overall, this work contributes to the growing understanding of mammalian body shape evolution and demonstrates the importance of not omitting body size in ecomorphological analyses.

Understanding how anatomical structures evolve requires disentangling the roles of integration and modularity in shaping morphological variation. The vertebral column, a serially repeated and regionally differentiated structure, provides a powerful system for investigating these processes. Here, we examine how vertebral morphology evolves in relation to whole-body elongation across the adaptive radiation of Lake Malawi cichlid fishes. We tested for evolutionary integration between the precaudal and caudal domains, as well as assessed the contributions of vertebral count, centrum shape, and intervertebral spacing on body elongation. We find strong evolutionary integration between precaudal and caudal vertebral shape, with both vertebral shapes varying along shared axes of multivariate shape change. Despite this, precaudal and caudal vertebral counts evolve independently, indicating a decoupling between the evolution of identity and morphology. Whole-body elongation is significantly associated with coordinated changes in vertebral and rib morphology, including proportional increases in centrum size, posterior displacement of neural and haemal spines, and increased rib curvature. In contrast, centrum elongation and intervertebral spacing do not independently explain body elongation beyond vertebral counts. These results demonstrate that body elongation in cichlids necessitates integrated, multivariate changes in axial morphology. Our findings highlight the importance of morphological integration in facilitating coordinated evolutionary responses in anatomical systems.

Immunoglobulin genes are a central component of jawed-vertebrate adaptive immunity. A previous study showed that the blunt-snouted clingfish Gouania willdenowi lacks immunoglobulin genes and T-cell receptor gamma/delta loci, while retaining T-cell receptor alpha/beta genes, MHC genes, and RAG1 /RAG2. Here I extend that observation to the family Gobiesocidae using all seven chromosome-level Gobiesocidae genome assemblies currently available. Manual tblastn and synteny-guided searches found no convincing immunoglobulin heavy-chain or light-chain loci in G. willdenowi, Gouania pigra, Gobiesox punctulatus, Apletodon dentatus, Lepadogaster candolii, Lepadogaster purpurea, or Diplecogaster bimaculata. Thus, the absence of antibody genes is best interpreted as a root-level character of clingfishes. The latest seven-species screen of 40 additional immune-associated genes shifts the broader interpretation in the same direction: the B-cell/adaptive core genes CD79A, CD79B, CIITA, TNFRSF13B, and TNFSF13B lack strong tblastn support in all sampled Gobiesocidae, and 37 of the 40 tested targets show an all-zero binary pattern at the presence threshold. Only IL21R.1, TYROBP, and TNFRSF11A show strong hits in one or more species. I therefore interpret the principal immune-gene erosion as occurring at or near the Gobiesocidae root rather than as a recent Gouania-specific process, while keeping weak, paralog-sensitive, and patchy loci provisional. RAG2 comparisons show a shared Gobiesocidae PHD-domain C-to-S replacement in the zinc-binding motif, with apparently intact RAG2 coding sequence. A family-wide TRG/TRD screen did not recover TRGV V segments or accepted TRDC constant-region exons, but it did detect TRGC-like constant exons in several genomes. These TRGC-like sequences are probably not canonical TRG constant exons without further validation, so I treat the gamma/delta system as eroded or rearranged rather than as a complete root-level loss equivalent to the Ig loss. The RAG2 variant provides a plausible molecular context for antigen-receptor remodeling, but it is not evidence that RAG genes are pseudogenized, because TCR alpha/beta, MHC genes, and RAG1 /RAG2 are retained. Gobiesocidae are therefore best described as a vertebrate family with ancestral loss of canonical immunoglobulin genes and associated root-level erosion of B-cell and immune-related genes, not as a lineage lacking adaptive immunity in its entirety. HighlightsO_LISeven chromosome-level Gobiesocidae genomes lack convincing canonical IgH and IgL loci. C_LIO_LIThe strongest non-Ig losses map to the B-cell/adaptive core: CD79A, CD79B, CIITA, TNFRSF13B, and TNFSF13B. C_LIO_LITCR alpha/beta, MHC genes, and RAG1 /RAG2 are retained, so Gobiesocidae should not be described as lacking adaptive immunity in full. C_LIO_LIA shared Gobiesocidae RAG2 PHD-domain C-to-S variant provides candidate molecular context for antigen-receptor remodeling. C_LI

Mammalian limb development is a complex system involving several signaling centers and coordinated cell behaviors to sculpt a functioning limb capable of the diverse locomotory strategies that mammals exhibit. To investigate the changes in development that facilitate the generation of the wide array of limb phenotypes across mammals, we take a correlation network approach to investigate the developing limbs of mice, bats, and opossums, which represent typical limb development, a novel limb phenotype, and a shift in developmental timing, respectively. Using transcriptomic data of early limb development across these taxa, we build module correlation networks and identify a difference in network connectivity and the distribution of limb development genes across bat limb development. We identify a unique signature of increased modularity in the bat forelimb that is not detected in mouse or opossum. This modularity is not associated with increased specialization of limb development modules, but rather is marked by target limb development genes being spread evenly across several modules. The opossum, with its standard phenotype but altered developmental timing, does not show a difference in modularity relative to mouse. This work points toward the benefit of a network-minded approach to transcriptomic networks, which reveals developmental modularity and potential gene targets for exploration of developmental system evolution.

Sex ratio distorters (SRDs) are heritable elements that bias offspring sex ratios to enhance their transmission. In the terrestrial isopod Armadillidium vulgare, feminization of genetic males can occur through vertical transmission of the sex ratio distorter known as the f-element, as well as through infection by Wolbachia, a maternally inherited bacterial endosymbiont that can alter host reproduction. Previous studies have focused on the distribution of SRDs and their associations with mitochondrial haplotypes in native European populations, but these patterns are poorly understood in the United States. In this study, we sampled A. vulgare in 12 states, screening individuals for Wolbachia infection, the presence of the f-element, and mitochondrial haplotypes. We found that Wolbachia shows a heterogeneous distribution across populations and haplotypes, in contrast with stronger associations in Europe. The f-element occurred in lower overall frequencies but showed a strong association with mitochondrial haplotype VI. These results indicate that patterns associated with SRD differ from those observed in Europe and suggest that multiple introductions and population mixing have shaped these distributions.

Despite recent explorations of tissue-level epithelial morphogenesis, how the embryonic brain wall, an epithelial derivative unique in its extensive cellular stratification coupled with epithelial-cell heightening to [~]0.5 mm, achieves thickening is unknown. Furthermore, the role of the inner curvature of the wall in this process remains unclear. Thus, whether the apically concave dorsal cerebral (pallial) wall thickens inward, not only outward as previously thought, was examined in culture and in vivo. The pallial wall in the midembryonic period, but not the earlier pallial wall, thickened inward at a pace equivalent to that of neuronal accumulation, which was found via stress-release tests to proceed in a compressive manner. The inhibition of actomyosin-mediated contraction of the inner/apical surface prevented the pallium from thickening, while inducing apical-surface buckling, suggesting the necessity of this inward thickening to actively avoid overcrowding. In contrast, inward bulging of the apically convex ganglionic eminence occurs more passively via pushing of the low-actomyosin apical surface by the constituent cells.

Mosquitoes are classified into two subfamilies, each monophyletic, and typically considered to both be ancient, having diverged more than 100 million years ago based on previous divergence analyses. A recent publication challenged this view with phylogenomic results primarily from the third codon position and UCEs. Utilizing alternative fossil placement and these phylogenomic data, these authors find that the Culicidae and Chaoboridae diverged in the lower Cretaceous, and that one mosquito subfamily, the Anophelinae, is nested within the Culicinae. These results are in stark contrast to previous results from diverse data sources, ranging from other genomic data, to morphology, to fossils. Here, we briefly detail the substantial evidence that supports two monophyletic subfamilies of extant mosquitoes, along with fossil evidence that supports the ancient divergence of these lineages.

Organisms have evolved a remarkable diversity of reproductive strategies in response to environmental variations and selective pressures. Although most vertebrates do reproduce biparentally, rare alternative modes such as selfing (self-fertilization) and different forms of parthenogenesis exist, but remain poorly characterized. Here, we investigated an unusual reproductive event in the normally biparental cichlid fish Cyphotilapia frontosa, in which a female produced offspring in the absence of a male. Using whole-genome sequencing data, we analyzed whether reproduction occurred via selfing or parthenogenesis by comparing patterns of heterozygosity with those from a wild, genetically diverse C. frontosa family collected in Lake Tanganyika and a closely related inbred Ctenochromis benthicola family. The uniparental family exhibited reduced genetic diversity, elevated relatedness, and genome-wide patterns of homozygosity distinct from those expected under parthenogenesis or inbreeding, but consistent with self-fertilization. Our study provides rare genomic evidence of selfing in a vertebrate and suggests that such alternative reproductive modes may be overlooked rather than truly absent. These findings contribute to a broader understanding of how alternative reproductive strategies evolve in vertebrate lineages. SignificanceThe overwhelming majority of vertebrates reproduce sexually, requiring a male and a female to produce genetically distinct offspring. Yet, rare alternative modes involving only a single parent such as asexual parthenogenesis ("virgin birth") or self-fertilization challenge this paradigm. Among these, selfing is exceptionally uncommon and poorly studied in vertebrates. Here, we unveiled - based on genomic analyses - the reproductive strategy of a member of the extraordinarily diverse cichlid fish radiation in Lake Tanganyika that reproduced in captivity in the absence of a male. By comparing patterns of genome-wide heterozygosity with both wild and inbred reference families, we identified a rare case of selfing. This finding adds to the limited records of selfing in vertebrates and expands current understanding of reproductive diversity, highlighting the power of whole-genome sequencing to distinguish among alternative reproductive mechanisms.

Idiopathic scoliosis is a common spinal disorder characterized by progressive three-dimensional curvature of unknown cause. Although biomechanical imbalance has long been proposed to contribute to scoliosis, the early physiological states that precede curvature onset remain poorly understood. Here, we investigated this problem using zebrafish uts2r3 mutants, which develop fully penetrant juvenile-onset spinal curvature following disruption of urotensin signaling. Transcriptomic analysis before curvature revealed altered expression of muscle-associated genes, suggesting that Uts2r3 influences axial muscle development or function. However, immunofluorescence, birefringence imaging, and quantitative analysis of myotome morphology showed that mutants lack overt muscle architectural defects or dystrophic pathology. By contrast, direct measurements of isolated larval trunks revealed pre-curvature biomechanical abnormalities: namely, uts2r3 mutants generated reduced active force following electrical stimulation while also exhibiting increased passive resistance to stretch. These findings identify urotensin signaling as a regulator of axial tissue biomechanics during growth and suggest that scoliosis-like curvature can arise from an early imbalance between active force generation and passive tissue stiffness. SignificanceSpinal curvature is common, but the biological events that cause the spine to bend during growth remain poorly understood. Animal models, especially zebrafish, make it possible to study these events before curvature begins. Zebrafish lacking urotensin signaling develop spinal curves that arise during juvenile growth, similar to idiopathic scoliosis in humans. Here, we demonstrate that zebrafish lacking the urotensin pathway receptor Uts2r3 develop an abnormal biomechanical state prior to curve onset. Their axial tissues generate less active force when contracting and, at the same time, show increased passive resistance to stretch--an unexpected combination that reveals a distinct pre-curvature biomechanical state. These findings suggest that spinal curvature can arise from an early imbalance in tissue mechanics during growth and identify urotensin signaling as a pathway that helps preserve spinal morphology through a biomechanical mechanism.