The gustatory stalk of the Remo flounder exemplifies how complex evolutionary novelties may arise.

The appearance of evolutionary novelties is a central issue in biology. Since Darwin's theory, difficulties in explaining how novel intricate body parts arose have often been used by creationists and other deniers to challenge evolution. Here, we describe the gustatory stalk of the Remo flounder (Oncopterus darwinii), an anatomically and functionally complex organ presumably used as a chemoreceptor probe to detect prey buried in the substrate. We demonstrate that the gustatory stalk is derived from the first dorsal-fin ray, which acquired remarkable modifications in its external morphology, integument, skeleton, muscles, and nerves. Such structural innovations are echoed in both functional and ecological specializations. We reveal that the gustatory stalk arose through the gradual accumulation of changes that evolved at different levels of the phylogenetic tree of ray-finned fishes. At least five preconditions arose in nodes preceding Oncopterus darwinii. This finding constitutes an interesting example of how evolution can deeply remodel body parts to perform entirely new functions. In this case, a trivial support structure primitively used for swimming became a sophisticated sensory tool to uncover hidden prey.

A rare occurrence of reversal in the four-spotted megrim Lepidorhombus boscii (Risso 1810) (Pleuronectiformes: Scophthalmidae) in the central Adriatic Sea.

One reversal dextral (right-sided) specimen occurred in the sample of 1070 individuals of four-spotted megrim Lepidorhombus boscii caught in the central Adriatic Sea from July 2020 to June 2021. This is the first documented occurrence of reversal in this species in the Adriatic Sea. Besides reverse symmetry, the individual showed similar morphological and morphometric characteristics as non-reversed individuals.

Involvement of glycolysis activation in flatfish sexual size dimorphism: Insights from transcriptomic analyses of Platichthys stellatus and Cynoglossus semilaevis.

The starry flounder (Platichthys stellatus), a flatfish cultured at the margins of the North Pacific, displays an obvious female-biased growth advantage, similar to many other fish species. To reveal the molecular mechanism underlying sexual size dimorphism, a comparative transcriptomic analysis of the somatotropic and reproductive axes was conducted. In total, 156, 67, 3434, and 378 differentially expressed genes (DEGs) between female and male samples were obtained in the brain, liver, gonad, and muscle tissues (q < 0.05). These DEGs were significantly enriched for various GO terms, including ion channel activity, protein binding, lipid transporter activity, and glycolytic process. The significantly enriched KEGG pathways included insulin secretion, axon guidance, and glycolysis/gluconeogenesis. In a detailed analysis of DEGs in these significantly enriched pathways, 35 genes showed higher expression levels in female muscle tissues than in male muscle tissues. A protein-protein interaction network further revealed specific interactions involving the glycolysis related-protein enolase (ENO), triosephosphate isomerase (TPI), Bisphosphoglycerate mutase (BPGM), fructose-bisphosphate aldolase (ALDO), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH). Interestingly, the role of glycolysis/gluconeogenesis was supported by an analysis of common DEGs between P. stellatus and Chinese tongue sole (Cynoglossus semilaevis). These results indicate that the activation of glycolysis in female muscle tissues contributes to flatfish sexual size dimorphism.

A comprehensive description of oocyte developmental stages in Pacific halibut, Hippoglossus stenolepis.

Accurate characterization of oocyte development is essential to understanding foundational aspects of reproductive biology and successful management of Pacific halibut (Hippoglossus stenolepis). Here this study provides complete histological descriptions for eight oocyte developmental stages in addition to postovulatory follicles and demonstrates the potential for oocyte size frequency distribution to act as a proxy for ovarian developmental stage and future maturity assessments. Importantly, it provides the first histological evidence that Pacific halibut have a group-synchronous ovarian developmental pattern with determinate fecundity and support for their batch-spawning strategy.

An hypoxia-tolerant flatfish: consequences of sustained stress on the slender sole Lyopsetta exilis (Pleuronectidae) in the context of a changing ocean.

Slender sole Lyopsetta exilis is an abundant groundfish on the continental shelf and inner waters of British Columbia, Canada, where it reaches a maximum standard length of 44 cm. Benthic image surveys coupled with oxygen measurements in Saanich Inlet document a dense population in bottom conditions near anoxia (0.03 ml l-1 oxygen) where diel migrating zooplankton intersect the bottom; we confirm this species is a planktivore, which limits its depth range to the base of the migration layer. In a comparison with slender sole from a nearby well-oxygenated habitat, several probable effects of living in severe hypoxia emerge: both sexes are significantly smaller in Saanich and the sex ratio is male-skewed. Otoliths from the Saanich fish were difficult to read due to many checks, but both sexes were smaller at age with the largest female (20 cm) from the hypoxia zone registering 17 years. Hypoxia appears to have a direct consequence on growth despite good food supply in this productive basin. Hyperventilation, a low metabolic rate and a very low critical oxygen tension help this fish regulate oxygen uptake in severely hypoxic conditions; it will be particularly resilient as the incidence of hypoxia increases on the continental shelf. Data from small-mesh bottom-trawl surveys over four decades reveal an increase in mean annual catch per unit effort in southern regions of the province, including the outer shelf and the Strait of Georgia. The California Cooperative Oceanic Fisheries Investigations (CalCOFI) ichthyoplankton database records a general decline in fish larvae on the Oregon-California shelf since 1990, but slender sole larvae are increasing there, as they are in the Strait of Georgia. We project that the slender sole will gain relative benefits in the future warming, deoxygenated northeast Pacific Ocean.

Pterobothrium sp. parasitando linguado, Xystreurys rasile (Jordan, 1891), coletados no litoral do Estado do Rio de Janeiro, Brasil e sua importância na inspeção do pescado / Pterobothrium sp. parasitizing flounders, Xystreurys rasile (Jordan, 1891), collected off the coast the State of Rio de Janeiro, Brazil and its fishing inspection importance

Alguns parasitos, podem conferir aspecto repugnante ao pescado, sendo condenados pelo Serviço de Inspeção, quando na indústria de beneficiamento. Entre novembro de 2010 e maio de 2018 foram adquiridos 36 espécimes de Xystreurys rasile em mercados de pescados no município do Rio de Janeiro, Brasil. Os peixes foram identificados, medidos, pesados e necropsiados. Os cestoides foram processados segundo as técnicas usuais em helmintologia. Cinco linguados estavam parasitados por 29 cestoides do gênero Pterobothrium, os índices parasitários foram: prevalência (14%), intensidade média (5,8), abundância média (0,8) e amplitude de infecção (2-12) e os sítios de infecção foram estômago, intestino, rim, fígado, serosas do estômago, baço e fígado e cavidade abdominal. Este é o primeiro registro de Pterobothrium sp. parasitando X. rasile.

Skeletal development of the chondrocranium in the tongue sole Cynoglossus semilaevis (Pleuronectiformes: Cynoglossidae).

This study provides a comprehensive description of chondrocranial development before, during and after larval metamorphosis in the tongue sole Cynoglossus semilaevis, a commercially valuable flatfish in China. Samples were collected at regular intervals ranging from 1 to 23 days post hatching (dph). Based on observations of cleared and double-stained specimens and images from sections stained with safranin O-fast green, major morphological events during early development were described. No cartilaginous structure was visible at hatching. From 2 dph onwards, cartilaginous structures such as the trabecular bar and some elements of the mandibular, hyoid and branchial arches appeared. At this time also, cartilaginous structures of the neurocranium started to form. Hypertrophic chondrocytes could be observed in many splanchnocranium elements at 5 dph. The start of ossification was indicated by alizarin red stain visible at 14 dph. At 17 dph, most of the cartilaginous skeleton was ossified. Soon after, the right eye started to migrate and pass through a slit beneath the dorsal-fin base and above the skull. Metamorphosis was complete at 20 dph, at which time the dorsal-fin base cartilage extended onto the anterior region of the head. Meanwhile, extremities of the hyoid and branchial arch elements remained cartilaginous. At 23 dph, endochondral ossification of the splanchnocranium was nearly complete. Unlike previous observations of other Pleuronectiformes, our study indicates that endochondral ossification of C. semilaevis skull cartilage occurs before metamorphosis.

Automatic interpretation of otoliths using deep learning.

The age structure of a fish population has important implications for recruitment processes and population fluctuations, and is a key input to fisheries-assessment models. The current method of determining age structure relies on manually reading age from otoliths, and the process is labor intensive and dependent on specialist expertise. Recent advances in machine learning have provided methods that have been remarkably successful in a variety of settings, with potential to automate analysis that previously required manual curation. Machine learning models have previously been successfully applied to object recognition and similar image analysis tasks. Here we investigate whether deep learning models can also be used for estimating the age of otoliths from images. We adapt a pre-trained convolutional neural network designed for object recognition, to estimate the age of fish from otolith images. The model is trained and validated on a large collection of images of Greenland halibut otoliths. We show that the model works well, and that its precision is comparable to documented precision obtained by human experts. Automating this analysis may help to improve consistency, lower cost, and increase the extent of age estimation. Given that adequate data are available, this method could also be used to estimate age of other species using images of otoliths or fish scales.

Comparative performance of growth, vertebral structure and muscle composition in diploid and triploid Paralichthys olivaceus.

Growth, skeletal structure and muscle composition of cold-shock-induced triploid olive flounder Paralichthys olivaceus were investigated. The average values of total length and total weight of triploids were higher than those of diploids from 5 to 11 months posthatch (mph). The growth difference disappeared after 11 mph. The skeletal structure of flounder at 11 mph was observed by X-ray imaging method. There are four kinds of vertebral deformity including vertebrae fusion, one-sided compression, two-sided compression and vertically shifted. The trunk region (V8-18) and tailing end of the vertebral column were the predominant locations of deformity. In general, the frequencies of vertebral deformities in triploids (60.0%) were higher than those in diploids (33.3%, p < 0.05). Both the number of fish with deformed vertebrae and the average frequencies of deformed vertebrae in triploids were significantly higher than those in diploids (p < 0.05). The muscle tissues of diploid and triploid flounder at 11 mph contain the same types of fatty acid and free amino acid profiles. The number of fatty acids with significant higher contents in diploids and triploids was one and ten, respectively (p < 0.05). The contents of free amino acids showed no difference between triploid and diploid fish.

Effects of organism and substrate size on burial mechanics of English sole, Parophrys vetulus.

Flatfishes use cyclic body undulations to force water into the sediment and fluidize substrate particles, displacing them into the water column. When water velocity decreases, suspended particles settle back onto the fish, hiding it from view. Burial may become more challenging as flatfishes grow because the area to be covered increases exponentially with the second power of length. In addition, particle size is not uniform in naturally occurring substrates, and larger particles require higher water velocities for fluidization. We quantified the effects of organism and particle-size scaling on burial behavior of English sole, Parophrys vetulus We recorded burial events from a size range of individuals (5-32 cm total length, TL), while maintaining constant substrate grain size. Larger fish used lower cycle frequencies and took longer to bury, but overall burial performance was maintained (∼100% coverage). To test the effect of particle size on burial performance, individuals of similar lengths (5.7-8.1 cm TL) were presented with different substrate sizes (0.125-0.710 mm). Particle size did not affect cycle frequency or time to burial, but fish did not achieve 100% coverage with the largest particles because they could not fluidize this substrate. Taken together, these results suggest that both body size and substrate grain size can potentially limit the ability of flatfishes to bury: a very large fish (>150 cm) may move too slowly to fluidize all but the smallest substrate particles and some particles are simply too large for smaller individuals to fluidize.

Population specific sperm production in European flounder Platichthys flesus: Adaptation to salinity at spawning.

Marine teleosts inhabiting the brackish Baltic Sea have adapted to the less saline water with activation of spermatozoa at low salinity hypo-osmotic conditions but with shorter longevity and lower swimming speed that affect the fertilization capacity. Aiming to elucidate if the fertilization capacity may be maintained by increasing the number of spermatozoa produced, testis size for the euryhaline flounder Platichthys flesus with external fertilization was assessed along a salinity gradient; with spawning at a salinity of c. 7, 10-18 and 30-35. Fulton's condition factor K = 0.881 ± 0.085 (mean ± S.D.), 0.833 ± 0.096 and 0.851 ± 0.086, for fish spawning at salinities of c. 7, 10-18 and 30-35, respectively, with no difference between areas, i.e. analysed fish were in similar nutritional condition. A general linear model, with testes dry mass as the dependent variable and somatic mass as covariate resulted in a significant difference between areas-populations with larger testes for P. flesus spawning at a salinity of c. 7 but no difference between fish spawning at a salinity of 10-18 and 30-35. The result suggests that adaptation by increasing the number of spermatozoa produced may be a key mechanism for marine teleosts spawning in areas with low salinities to sustain the fertilization capacity as shown here for the euryhaline P. flesus.

Structural and functional maturation of skin during metamorphosis in the Atlantic halibut (Hippoglossus hippoglossus).

To establish if the developmental changes in the primary barrier and osmoregulatory capacity of Atlantic halibut skin are modified during metamorphosis, histological, histochemical, gene expression and electrophysiological measurements were made. The morphology of the ocular and abocular skin started to diverge during the metamorphic climax and ocular skin appeared thicker and more stratified. Neutral mucins were the main glycoproteins produced by the goblet cells in skin during metamorphosis. Moreover, the number of goblet cells producing neutral mucins increased during metamorphosis and asymmetry in their abundance was observed between ocular and abocular skin. The increase in goblet cell number and their asymmetric abundance in skin was concomitant with the period that thyroid hormones (THs) increase and suggests that they may be under the control of these hormones. Several mucin transcripts were identified in metamorphosing halibut transcriptomes and Muc18 and Muc5AC were characteristic of the body skin. Na+, K+-ATPase positive (NKA) cells were observed in skin of all metamorphic stages but their number significantly decreased with the onset of metamorphosis. No asymmetry was observed between ocular and abocular skin in NKA cells. The morphological changes observed were linked to modified skin barrier function as revealed by modifications in its electrophysiological properties. However, the maturation of the skin functional characteristics preceded structural maturation and occurred at stage 8 prior to the metamorphic climax. Treatment of Atlantic halibut with the THs disrupter methimazole (MMI) affected the number of goblet cells producing neutral mucins and the NKA cells. The present study reveals that the asymmetric development of the skin in Atlantic halibut is TH sensitive and is associated with metamorphosis and that this barrier's functional properties mature earlier and are independent of metamorphosis.

The genome and transcriptome of Japanese flounder provide insights into flatfish asymmetry.

Flatfish have the most extreme asymmetric body morphology of vertebrates. During metamorphosis, one eye migrates to the contralateral side of the skull, and this migration is accompanied by extensive craniofacial transformations and simultaneous development of lopsided body pigmentation. The evolution of this developmental and physiological innovation remains enigmatic. Comparative genomics of two flatfish and transcriptomic analyses during metamorphosis point to a role for thyroid hormone and retinoic acid signaling, as well as phototransduction pathways. We demonstrate that retinoic acid is critical in establishing asymmetric pigmentation and, via cross-talk with thyroid hormones, in modulating eye migration. The unexpected expression of the visual opsins from the phototransduction pathway in the skin translates illumination differences and generates retinoic acid gradients that underlie the generation of asymmetry. Identifying the genetic underpinning of this unique developmental process answers long-standing questions about the evolutionary origin of asymmetry, but it also provides insight into the mechanisms that control body shape in vertebrates.

Production and verification of a 2nd generation clonal group of Japanese flounder, Paralichthys olivaceus.

Clonal fishes are useful tools in biology and aquaculture studies due to their isogenicity. In Japanese flounder (Paralichthys olivaceus), a group of homozygous clones was created by inducing meiogynogenesis in eggs from a mitogynogenetic homozygous diploid. As the clones reached sexual maturity, meiogynogenesis was again induced in order to produce a 2nd generation clonal group of Japanese flounder. After 3 months, there were 611 healthy, surviving individuals. Twenty-four microsatellite markers, that covered all the linkage groups of Japanese flounder, were used to identify the homozygosity of the 2nd generation clones; no heterozygous locus was detected. This indicates that the production of a 2nd generation clonal group of Japanese flounder was successful. Restriction-site DNA associated sequencing at the genomic level also confirmed the homozygosity and clonality of the 2nd generation clonal group. Furthermore, these 2nd generation clones had a small coefficient of variation for body shape indices at 210 days of age and showed a high degree of similarity in body characteristics among individuals. The successful production of 2nd generation clones has laid the foundation for the large-scale production of clonal Japanese flounder.

Key role of scale morphology in flatfishes (Pleuronectiformes) in the ability to keep sand.

Flatfishes bury themselves for camouflage and protection. Whereas species-specific preferences for certain sediments were previously shown, the role of scales in interaction with sediment has not been investigated. Here, scale morphology and sediment friction were examined in four European pleuronectiforms: Limanda limanda, Platichthys flesus, Pleuronectes platessa, and Solea solea. All species had different scale types ranging from cycloid to ctenoid scales. On the blind side, the number of scales is higher and scales have less ctenial spines than on the eye side. The critical angle of sediment sliding (static friction) significantly depended on the grain size and was considerably higher on the eye side. The effect of mucus was excluded by repeated measurements on resin replicas of the skin. Our results demonstrate the impact of scale morphology on sediment interaction and give an insight about the ability of scales to keep sand. Exposed scales and a higher number of ctenial spines on the eye side lead to an increase of friction forces, especially for sediments with a smaller grain size. Our results suggest that the evolution of scales was at least partly driven by their interactions with sediment which confirms the relevance of sediment for the distribution and radiation of Pleuronectiformes.

Characterization of F-spondin in Japanese flounder (Paralichthys olivaceus) and its role in the nervous system development of teleosts.

F-spondin was originally isolated from the developing embryonic floor plate of vertebrates, secreting numerous kinds of neuron-related molecules. The protein performs a positive function in nervous system development, which is attributed to the high conservation of F-spondin protein, an extracellular matrix (ECM) protein in several species. However, its precise function remains unknown, especially in marine fish. In this study, the F-spondin of Japanese flounder (Paralichthys olivaceus). was cloned, and its expression pattern and structural characteristics were analyzed. The 2421bp-long cDNA ORF of PoF-spondin was obtained and divided into 14 exons spread over 61,496bp of the genomic sequence. Phylogenetic analysis showed that PoF-spondin was actually the ortholog of the human spon1 gene and shared high identities with other teleost spon1a genes. Quantitative RT-PCR analysis showed that PoF-spondin was maternally expressed, and transcripts were present from one-cell stage to hatching stage, peaking at tailbud stage. Tissue distribution analysis indicated that PoF-spondin was detectable mainly in the gonads (especially in the ovary) and the brain. Whole mount in situ hybridization analysis revealed that the PoF-spondin transcription distributed throughout the cleavage of the ball in the early stage and expressed at a high level in the floor plate of the trunk at tailbud and pre-hatching stages. Furthermore, the expression of genes related to nervous system development (spon1b, foxo3b, and foxj1a) was significantly increased after the injection of PoF-spondin into the embryos of wild-type zebrafish. Furthermore, PoF-spondin significantly suppressed the expression of the chordamesoderm marker gene ntl, increased the expression of otx2/krox20, ectoderm mark genes, and left the expression of dorsal mesodermal marker gene gsc unaffected at 50% epiboly stage in zebrafish. In short, our results suggest that PoF-spondin functions in the development of the teleost nervous system.

Bayesian Analysis for Genetic Architectures of Body Weights and Morphological Traits Using Distorted Markers in Japanese Flounder Paralichthys olivaceus.

A doubled haploids (DH) population of 160 individuals is constructed by using mitotic gynogenetics for the study of genetic architectures of the body weight and morphological traits in Japanese flounder. Each DH individual is genotyped for 458 SSR markers, 222 of which segregate distortionally. By modifying conditional probabilities of quantitative trait locus (QTL) genotypes on the distorted flanking markers, Bayesian model selection is used to dissect genetic architectures for the traits. As a result, we identify 42 main-effect QTLs on chromosomes 5, 6, 7, 8, 9, 10, 15, 20, 21, 22, and 59 pairs of interacting QTLs. Among these detected QTLs, the largest interacting QTL is between chromosome 6 and chromosome 9 and accounts for 25.196 % of phenotypic variance for body weights and in a similar trend. Also, many QTLs show pleiotropic effects. The QTL on chromosome 9 simultaneously governs seven traits, BL, BH, FL, HL, PFL, HW, and CW. As compared to method using the uncorrected conditional probabilities of QTL genotypes, our method using corrected conditional probabilities can detect more interacting QTLs for the traits.

High-density aggregations of rodlet cells in the gonads of Greenland halibut Reinhardtius hippoglossoides, a deep-water marine flatfish.

Large aggregations of rodlet cells in the gonads of male and female Greenland halibut Reinhardtius hippoglossoides are reported for the first time. These rodlet cells were not arranged epithelially but rather were found throughout the connective tissue between oocytes (females) or within lymphatic spaces between testicular lobules (males). The reason for large aggregations of rodlet cells in the gonads and not other tissues of this species is uncertain.

Genetic analysis for main length ratio associated with morphological traits in Japanese flounder Paralichthys olivaceus.

The objective of this study was to estimate genetic parameters of main length ratios associated with morphological traits in Paralichthys olivaceus. At 180, 240 and 360 days of age, body length (L1 ), head length (L2 ), body depth (L3 ), caudal peduncle depth (L4 ) and caudal-peduncle length (L5 ) were observed in 4757-5224 individuals from 52 full-sib families. Based on these observed traits, the main length-ratio traits including L3 :L1 , L2 :L1 , L5 :L1 and L4 :L3 were constructed. Three single-trait animal models were used to examine the phenotypic variation of length-ratio traits and were compared using the likelihood ratio test. As a result, the optimal model for all the analysed traits was determined as the model with the fixed, additive genetic and full-sib family effects. Using the optimal model, animal model analysis for single trait showed that heritabilities were estimated to be low, ranging from 0·122 to 0·276. The estimates for all length-ratio traits decreased while the days of age increased except for L3 :L1 . Multivariate animal model analysis indicated that within and among days of age, genetic correlations between length-ratio traits were different in the direction and magnitude. Genetic parameter estimation of main length ratios would be helpful for genetic improvement of body shape in P. olivaceus.