A Triassic stem lepidosaur illuminates the origin of lizard-like reptiles.

The early evolution of diapsid reptiles is marked by a deep contrast between our knowledge of the origin and early evolution of archosauromorphs (crocodiles, avian and non-avian dinosaurs) to that of lepidosauromorphs (squamates (lizards, snakes) and sphenodontians (tuataras)). Whereas the former include hundreds of fossil species across various lineages during the Triassic period1, the latter are represented by an extremely patchy early fossil record comprising only a handful of fragmentary fossils, most of which have uncertain phylogenetic affinities and are confined to Europe1-3. Here we report the discovery of a three-dimensionally preserved reptile skull, assigned as Taytalura alcoberi gen. et sp. nov., from the Late Triassic epoch of Argentina that is robustly inferred phylogenetically as the earliest evolving lepidosauromorph, using various data types and optimality criteria. Micro-computed tomography scans of this skull reveal details about the origin of the lepidosaurian skull from early diapsids, suggesting that several traits traditionally associated with sphenodontians in fact originated much earlier in lepidosauromorph evolution. Taytalura suggests that the strongly evolutionarily conserved skull architecture of sphenodontians represents the plesiomorphic condition for all lepidosaurs, that stem and crown lepidosaurs were contemporaries for at least ten million years during the Triassic, and that early lepidosauromorphs had a much broader geographical distribution than has previously been thought.

Handling Logical Character Dependency in Phylogenetic Inference: Extensive Performance Testing of Assumptions and Solutions Using Simulated and Empirical Data.

Logical character dependency is a major conceptual and methodological problem in phylogenetic inference of morphological data sets, as it violates the assumption of character independence that is common to all phylogenetic methods. It is more frequently observed in higher-level phylogenies or in data sets characterizing major evolutionary transitions, as these represent parts of the tree of life where (primary) anatomical characters either originate or disappear entirely. As a result, secondary traits related to these primary characters become "inapplicable" across all sampled taxa in which that character is absent. Various solutions have been explored over the last three decades to handle character dependency, such as alternative character coding schemes and, more recently, new algorithmic implementations. However, the accuracy of the proposed solutions, or the impact of character dependency across distinct optimality criteria, has never been directly tested using standard performance measures. Here, we utilize simple and complex simulated morphological data sets analyzed under different maximum parsimony optimization procedures and Bayesian inference to test the accuracy of various coding and algorithmic solutions to character dependency. This is complemented by empirical analyses using a recoded data set on palaeognathid birds. We find that in small, simulated data sets, absent coding performs better than other popular coding strategies available (contingent and multistate), whereas in more complex simulations (larger data sets controlled for different tree structure and character distribution models) contingent coding is favored more frequently. Under contingent coding, a recently proposed weighting algorithm produces the most accurate results for maximum parsimony. However, Bayesian inference outperforms all parsimony-based solutions to handle character dependency due to fundamental differences in their optimization procedures-a simple alternative that has been long overlooked. Yet, we show that the more primary characters bearing secondary (dependent) traits there are in a data set, the harder it is to estimate the true phylogenetic tree, regardless of the optimality criterion, owing to a considerable expansion of the tree parameter space. [Bayesian inference, character dependency, character coding, distance metrics, morphological phylogenetics, maximum parsimony, performance, phylogenetic accuracy.].

The origin of squamates revealed by a Middle Triassic lizard from the Italian Alps.

Modern squamates (lizards, snakes and amphisbaenians) are the world's most diverse group of tetrapods along with birds 1 and have a long evolutionary history, with the oldest known fossils dating from the Middle Jurassic period-168 million years ago2-4. The evolutionary origin of squamates is contentious because of several issues: (1) a fossil gap of approximately 70 million years exists between the oldest known fossils and their estimated origin5-7; (2) limited sampling of squamates in reptile phylogenies; and (3) conflicts between morphological and molecular hypotheses regarding the origin of crown squamates6,8,9. Here we shed light on these problems by using high-resolution microfocus X-ray computed tomography data from the articulated fossil reptile Megachirella wachtleri (Middle Triassic period, Italian Alps 10 ). We also present a phylogenetic dataset, combining fossils and extant taxa, and morphological and molecular data. We analysed this dataset under different optimality criteria to assess diapsid reptile relationships and the origins of squamates. Our results re-shape the diapsid phylogeny and present evidence that M. wachtleri is the oldest known stem squamate. Megachirella is 75 million years older than the previously known oldest squamate fossils, partially filling the fossil gap in the origin of lizards, and indicates a more gradual acquisition of squamatan features in diapsid evolution than previously thought. For the first time, to our knowledge, morphological and molecular data are in agreement regarding early squamate evolution, with geckoes-and not iguanians-as the earliest crown clade squamates. Divergence time estimates using relaxed combined morphological and molecular clocks show that lepidosaurs and most other diapsids originated before the Permian/Triassic extinction event, indicating that the Triassic was a period of radiation, not origin, for several diapsid lineages.

Evaluating the Performance of Probabilistic Algorithms for Phylogenetic Analysis of Big Morphological Datasets: A Simulation Study.

Reconstructing the tree of life is an essential task in evolutionary biology. It demands accurate phylogenetic inference for both extant and extinct organisms, the latter being almost entirely dependent on morphological data. While parsimony methods have traditionally dominated the field of morphological phylogenetics, a rapidly growing number of studies are now employing probabilistic methods (maximum likelihood and Bayesian inference). The present-day toolkit of probabilistic methods offers varied software with distinct algorithms and assumptions for reaching global optimality. However, benchmark performance assessments of different software packages for the analyses of morphological data, particularly in the era of big data, are still lacking. Here, we test the performance of four major probabilistic software under variable taxonomic sampling and missing data conditions: the Bayesian inference-based programs MrBayes and RevBayes, and the maximum likelihood-based IQ-TREE and RAxML. We evaluated software performance by calculating the distance between inferred and true trees using a variety of metrics, including Robinson-Foulds (RF), Matching Splits (MS), and Kuhner-Felsenstein (KF) distances. Our results show that increased taxonomic sampling improves accuracy, precision, and resolution of reconstructed topologies across all tested probabilistic software applications and all levels of missing data. Under the RF metric, Bayesian inference applications were the most consistent, accurate, and robust to variation in taxonomic sampling in all tested conditions, especially at high levels of missing data, with little difference in performance between the two tested programs. The MS metric favored more resolved topologies that were generally produced by IQ-TREE. Adding more taxa dramatically reduced performance disparities between programs. Importantly, our results suggest that the RF metric penalizes incorrectly resolved nodes (false positives) more severely than the MS metric, which instead tends to penalize polytomies. If false positives are to be avoided in systematics, Bayesian inference should be preferred over maximum likelihood for the analysis of morphological data.

Sphenodontian phylogeny and the impact of model choice in Bayesian morphological clock estimates of divergence times and evolutionary rates.

BACKGROUND: The vast majority of all life that ever existed on earth is now extinct and several aspects of their evolutionary history can only be assessed by using morphological data from the fossil record. Sphenodontian reptiles are a classic example, having an evolutionary history of at least 230 million years, but currently represented by a single living species (Sphenodon punctatus). Hence, it is imperative to improve the development and implementation of probabilistic models to estimate evolutionary trees from morphological data (e.g., morphological clocks), which has direct benefits to understanding relationships and evolutionary patterns for both fossil and living species. However, the impact of model choice on morphology-only datasets has been poorly explored. RESULTS: Here, we investigate the impact of a wide array of model choices on the inference of evolutionary trees and macroevolutionary parameters (divergence times and evolutionary rates) using a new data matrix on sphenodontian reptiles. Specifically, we tested different clock models, clock partitioning, taxon sampling strategies, sampling for ancestors, and variations on the fossilized birth-death (FBD) tree model parameters through time. We find a strong impact on divergence times and background evolutionary rates when applying widely utilized approaches, such as allowing for ancestors in the tree and the inappropriate assumption of diversification parameters being constant through time. We compare those results with previous studies on the impact of model choice to molecular data analysis and provide suggestions for improving the implementation of morphological clocks. Optimal model combinations find the radiation of most major lineages of sphenodontians to be in the Triassic and a gradual but continuous drop in morphological rates of evolution across distinct regions of the phenotype throughout the history of the group. CONCLUSIONS: We provide a new hypothesis of sphenodontian classification, along with detailed macroevolutionary patterns in the evolutionary history of the group. Importantly, we provide suggestions to avoid overestimated divergence times and biased parameter estimates using morphological clocks. Partitioning relaxed clocks offers methodological limitations, but those can be at least partially circumvented to reveal a detailed assessment of rates of evolution across the phenotype and tests of evolutionary mosaicism.

CavVis-A Field-of-View Geometric Algorithm for Protein Cavity Detection.

Several geometric-based methods have been developed for the last two to three decades to detect and identify cavities (i.e., putative binding sites) on proteins, as needed to study protein-ligand interactions and protein docking. This paper introduces a new protein cavity method, called CavVis, which combines voxelization (i.e., a grid of voxels) and an analytic formulation of Gaussian surfaces that approximates the solvent-excluded surface. This method builds upon visibility of points on protein surface to find its cavities. Specifically, the visibility criterion combines three concepts we borrow from computer graphics, the field-of-view of each surface point, voxel ray casting, and back-face culling.

Giant taxon-character matrices II: a response to Laing et al. (2017).

The trend towards big data analyses in evolutionary biology has been observed in phylogenetics via the assembly of giant datasets composed of genomic and phenotypic data. We recently (Simões et al., 2017. Giant taxon-character matrices: Quality of character constructions remains critical regardless of size. Cladistics 33, 198-219) presented a critique of the phylogenetic character concepts used in current morphological datasets, with the caution that giant datasets did not obviate the empirical requirement of rigor in character construction. Laing et al. (2017. Giant taxon-character matrices: The future of morphological systematics. Cladistics, https://doi.org/10.1111/cla.12197) have since argued that we had 'suggested' that large datasets inherently contain flawed characters, and that we had presented a substandard methodology of phylogenetic analysis. Laing et al. concluded by discussing their approach to phylogenetic signal, total evidence and the inevitability of large datasets. We here reply to Laing et al. by reviewing what we actually wrote regarding dataset size, characters and methodology. We show that Laing et al.'s. central premise is unsupported, thus characterizing a Straw Man argument, and deeply misrepresents our original study. In part two, we discuss total evidence and phylogenetic signal issues raised by Laing et al. that are of major consequence to the appropriate construction of large morphological datasets.

New data on the Late Cretaceous lizard Dicothodon bajaensis (Squamata, Borioteiioidea) from Baja California, Mexico reveals an unusual tooth replacement pattern in squamates.

Borioteiioids comprise an extinct family of squamates that inhabited the Northern Hemisphere during the Cretaceous and were characterized by varying patterns of tooth replacement and dental morphology. Understanding the evolution of these tooth replacement patterns has, however, been largely hampered by an extremely fragmentary fossil record. Here we present new information on Dicothodon bajaensisfrom the Campanian of Baja California (Mexico), so far known only from isolated teeth and jaw fragments. Among abundant new materials there are ten maxillae and five dentaries belonging to distinct ontogenetic stages. Whereas juveniles display active tooth replacement, older specimens show no evidence of replacement. Dicothodon bajaensis is therefore inferred to have had arrested tooth replacement later in ontogeny. This provides the first evidence of cessation of tooth replacement during late ontogeny in lizards (living or extinct). This replacement type is also an evolutionary intermediate between the typical lizard replacement pattern, observed in some borioteiioids (e.g., Bicuspidon) and the complete absence of tooth replacement since early ontogeny as observed in other borioteiioids (e.g., Polyglyphanodon).

Giant taxon-character matrices: quality of character constructions remains critical regardless of size.

Giant morphological data matrices are increasingly common in cladistic analyses of vertebrate phylogeny, reporting numbers of characters never seen or expected before. However, the concern for size is usually not followed by an equivalent, if any, concern for character construction/selection criteria. Therefore, the question of whether quantity parallels quality for such influential works remains open. Here, we provide the largest compilation known to us of character construction methods and criteria, as derived from previous studies, and from our own de novo conceptualizations. Problematic character constructions inhibit the capacity of phylogenetic analyses to recover meaningful homology hypotheses and thus accurate clade structures. Upon a revision of two of the currently largest morphological datasets used to test squamate phylogeny, more than one-third of the almost 1000 characters analysed were classified within at least one of our categories of "types" of characters that should be avoided in cladistic investigations. These characters were removed or recoded, and the data matrices re-analysed, resulting in substantial changes in the sister group relationships for squamates, as compared to the original studies. Our results urge caution against certain types of character choices and constructions, also providing a methodological basis upon which problematic characters might be avoided.

Interrelationships among Early Triassic faunas of Western Gondwana and Laurasia as illuminated by a new South American benthosuchid temnospondyl.

The End-Permian Mass Extinction marked a critical turning point in Earth's history, and the biological recovery that followed the crisis led to the emergence of several modern vertebrate and invertebrate taxa. Even considering the importance of the Early Triassic biotic recovery for the evolution of modern faunas and floras, our knowledge of this event is still hindered by the sparse sampling of crucial geological formations. This leaves our understanding of Early Triassic ecosystems fundamentally biased toward productive and historically well-explored geological units. Recent surveys in poorly known Gondwanan localities, such as those within the Sanga do Cabral Formation in southern Brazil, have unveiled insights into Early Triassic terrestrial ecosystems, shedding light on a diverse and previously unknown tetrapod fauna. Here, we report the discovery of a new temnospondyl genus and species in the Lower Triassic Sanga do Cabral Formation. The new taxon can be confidently assigned to the Benthosuchidae, a stereospondyl clade with a distribution previously restricted to the East European Platform. Phylogenetic analysis confirms the relationship of the new genus to the trematosaurian lineage, being closely related to the genus Benthosuchus. Our results raise questions about the biogeographical history of stereospondyls after the End-Permian Mass Extinction and suggest a potential connection between Russian and South American Early Triassic faunas. Further investigations are needed to thoroughly explore the potential dispersal routes that may explain this seemingly unusual biogeographical pattern.

Metal accumulation in female green sea turtles (Chelonia mydas) from Eastern Atlantic affects their egg quality with potential implications for embryonic development.

Sea turtles, with their global distribution and complex life cycle, often accumulate pollutants such as metals and metalloids due to their extended lifespan and feeding habits. However, there are limited studies exploring the impact of metal pollution on the reproductive health of female sea turtles, specifically focusing on the quality of their eggs, which has significant implications for the future generations of these charismatic animals. São Tomé Island, a crucial nesting and feeding habitat for green sea turtles, underscores the urgent need for comprehensive research in this ecologically significant area. This study aimed to investigate whether metals and metalloids in the blood of nesting female green sea turtles induce genotoxic effects in their erythrocytes and affect their egg morphometric characteristics and the composition of related compartments. Additionally, this study aimed to evaluate whether the quality of energetic reserves for embryo development (fatty acids in yolk's polar and neutral lipids) is influenced by the contamination status of their predecessors. Results revealed correlations between Cu and Hg levels and increased "lobed" erythrocytes, while As and Cu negatively influenced shell thickness. In terms of energy reserves, both polar and neutral lipid fractions contained primarily saturated and monounsaturated fatty acids, with prevalent 18:1n-9, 18:0, 16:0, 14:0, and 12:0 fatty acids in yolk samples. The yolk polar fraction was more susceptible to contaminant levels in female sea turtles, showing consistent negative correlations between pollution load index and essential n3 fatty acids, including linolenic, eicosatrienoic, eicosapentaenoic, and docosapentaenoic acids, crucial for embryonic development. These metals accumulation, coupled with the reduced availability of these key fatty acids, may disrupt the eicosanoid and other important pathways, affecting reproductive development. This study reveals a negative correlation between metal contamination in female sea turtles' blood and egg lipid reserves, raising concerns about embryonic development and the species' future generations.

Liver and Plasma Fatty Acid Characterization in Cultured Brown Trout at Distinct Reproductive Stages.

Fatty acids are energy sources, and their profiles are used as biomarkers of metabolic status and physiological changes in fish. Within this context, the main aim of this study was to identify the fatty acids that best discriminate the reproductive status of male and female farmed brown trout. The fatty acid composition in liver and plasma samples from the adults of both sexes was monitored along four distinct reproductive stages, namely the spawning capable (December), regressing (March), regenerating (July), and developing (November) stages. Irrespective of the sex and stage, the most representative fatty acids were palmitic acid (16:0), stearic acid (18:0), oleic acid (18:1 n-9), arachidonic acid (20:4 n-6), eicosapentaenoic acid (EPA, 20:5 n-3), and docosahexaenoic acid (DHA, 22:6 n-3). There were no significant sex differences in fatty acid classes in the liver and plasma. Despite this, there were several changes in individual fatty acid levels between the sexes. In the liver, both males and females showed high monounsaturated fatty acid and low polyunsaturated fatty acid (PUFA) levels during the regressing and regenerating stages. At spawning capable and developing stages, a reverse profile was noted. The plasma profiles were mainly influenced by changes in saturated fatty acids and PUFAs in males and by PUFA in females. Based on the most representative fatty acids, four patterns were established for female plasma samples, one for each reproductive stage. This scenario suggests that female plasma samples are promising for the discrimination of gonadal reproductive status, and this potential can be further explored in aquaculture and environmental monitoring studies.

The affinities of the Late Triassic Cryptovaranoides and the age of crown squamates.

Most living reptile diversity is concentrated in Squamata (lizards, including snakes), which have poorly known origins in space and time. Recently, †Cryptovaranoides microlanius from the Late Triassic of the United Kingdom was described as the oldest crown squamate. If true, this result would push back the origin of all major lizard clades by 30-65 Myr and suggest that divergence times for reptile clades estimated using genomic and morphological data are grossly inaccurate. Here, we use computed tomography scans and expanded phylogenetic datasets to re-evaluate the phylogenetic affinities of †Cryptovaranoides and other putative early squamates. We robustly reject the crown squamate affinities of †Cryptovaranoides, and instead resolve †Cryptovaranoides as a potential member of the bird and crocodylian total clade, Archosauromorpha. Bayesian total evidence dating supports a Jurassic origin of crown squamates, not Triassic as recently suggested. We highlight how features traditionally linked to lepidosaurs are in fact widespread across Triassic reptiles. Our study reaffirms the importance of critically choosing and constructing morphological datasets and appropriate taxon sampling to test the phylogenetic affinities of problematic fossils and calibrate the Tree of Life.

Redescription of Tijubina pontei, an Early Cretaceous lizard (Reptilia; Squamata) from the Crato Formation of Brazil.

The record of Gondwanan Mesozoic lizards is very poor. Among the few species described for this region there is Tijubina pontei, an Early Cretaceous lizard from the Crato Formation (late Aptian) of northeast Brazil. Its description is very brief and lacks most of its diagnostic characters and clear delimitation from other lizard species. Here, a full redescription of the holotype is provided. T. pontei is demonstrated to be a valid species and a new diagnosis is provided with reference to Olindalacerta brasiliensis, a contemporary species of the Crato Formation. It lacks the posteroventral and posterodorsal processes of the dentary; the tibial/fibular length equals the femoral length and its posterior dentary teeth are robust, cylindrically based, unsculptured and bear no cuspids. The systematic position of T. pontei still needs further clarification, but preliminary analyses indicate that it lies in a rather basal position among the Squamata, similarly to O. brasiliensis.

An exceptionally preserved Sphenodon-like sphenodontian reveals deep time conservation of the tuatara skeleton and ontogeny.

Sphenodontian reptiles are an extremely old evolutionary lineage forming the closest relatives to squamates (lizards and snakes) and were globally distributed and more diverse than squamates during the first half of their evolutionary history. However, the majority of their fossils are highly fragmentary, especially within sphenodontines-the group including its single surviving species, Sphenodon punctatus (the tuatara of New Zealand)-thus severely hampering our understanding on the origins of the tuatara. Here, we present a new sphenodontian species from the Early Jurassic of North America (Arizona, USA) represented by a nearly complete articulated skeleton and dozens of upper and lower jaws forming the most complete ontogenetic series in the sphenodontian fossil record. CT-scanning provides plentitude of data that unambiguously place this new taxon as one of the earliest evolving and oldest known sphenodontines. Comparisons with Sphenodon reveal that fundamental patterns of mandibular ontogeny and skeletal architecture in Sphenodon may have originated at least ~190Mya. In combination with recent findings, our results suggest strong morphological stability and an ancient origin of the modern tuatara morphotype.

Successive climate crises in the deep past drove the early evolution and radiation of reptiles.

Climate change-induced mass extinctions provide unique opportunities to explore the impacts of global environmental disturbances on organismal evolution. However, their influence on terrestrial ecosystems remains poorly understood. Here, we provide a new time tree for the early evolution of reptiles and their closest relatives to reconstruct how the Permian-Triassic climatic crises shaped their long-term evolutionary trajectory. By combining rates of phenotypic evolution, mode of selection, body size, and global temperature data, we reveal an intimate association between reptile evolutionary dynamics and climate change in the deep past. We show that the origin and phenotypic radiation of reptiles was not solely driven by ecological opportunity following the end-Permian extinction as previously thought but also the result of multiple adaptive responses to climatic shifts spanning 57 million years.

An insight into the mechanisms underpinning the anti-browning effect of Codium tomentosum on fresh-cut apples.

This work focuses on understanding the action of a novel seaweed extract with anti-browning functionality in fresh-cut apples. Organic fresh-cut apples were coated by immersion in an aqueous Codium tomentosum seaweed extract (0.5 % w/v), packaged under ambient atmospheric conditions in plastic bags, and stored at 4 ˚C for 15 days. Browning-related enzymatic activities, as well as targeted gene expression related to superficial browning, were monitored immediately after coating and followed at five-day intervals, until a final storage period of 15 days. Gene expression was particularly affected one hour after coating application (day 0), with no expression registered for peroxidase (mdPOD) and phenylalanine ammonia-lyase (mdPAL) genes in the coated samples. A reduction in polyphenol oxidase expression levels was also observed. After 15 days of storage, the coated samples developed lower browning levels and presented distinctly lower activities of polyphenol oxidase and peroxidase - the oxidative enzymes predominantly involved in enzymatic browning. The observed post-coating suppression of mdPAL and mdPOD expression, and reduction in mdPPO expression, suggest that the seaweed C. tomentosum extract delays the activation of these genes, and decreases enzymatic activity, which in turn accounts for the coating's anti-browning effect.

Integrative Phylogenetics: Tools for Palaeontologists to Explore the Tree of Life.

The modern era of analytical and quantitative palaeobiology has only just begun, integrating methods such as morphological and molecular phylogenetics and divergence time estimation, as well as phenotypic and molecular rates of evolution. Calibrating the tree of life to geological time is at the nexus of many disparate disciplines, from palaeontology to molecular systematics and from geochronology to comparative genomics. Creating an evolutionary time scale of the major events that shaped biodiversity is key to all of these fields and draws from each of them. Different methodological approaches and data employed in various disciplines have traditionally made collaborative research efforts difficult among these disciplines. However, the development of new methods is bridging the historical gap between fields, providing a holistic perspective on organismal evolutionary history, integrating all of the available evidence from living and fossil species. Because phylogenies with only extant taxa do not contain enough information to either calibrate the tree of life or fully infer macroevolutionary dynamics, phylogenies should preferably include both extant and extinct taxa, which can only be achieved through the inclusion of phenotypic data. This integrative phylogenetic approach provides ample and novel opportunities for evolutionary biologists to benefit from palaeontological data to help establish an evolutionary time scale and to test core macroevolutionary hypotheses about the drivers of biological diversification across various dimensions of organisms.

Sustained high rates of morphological evolution during the rise of tetrapods.

The fish-to-tetrapod transition is one of the most iconic events in vertebrate evolution, yet fundamental questions regarding the dynamics of this transition remain unresolved. Here, we use advances in Bayesian morphological clock modelling to reveal the evolutionary dynamics of early tetrapodomorphs (tetrapods and their closest fish relatives). We show that combining osteological and ichnological calibration data results in major shifts on the time of origin of all major groups of tetrapodomorphs (up to 25 million years) and that low rates of net diversification, not fossilization, explain long ghost lineages in the early tetrapodomorph fossil record. Further, our findings reveal extremely low rates of morphological change for most early tetrapodomorphs, indicating widespread stabilizing selection upon their 'fish' morphotype. This pattern was broken only by elpistostegalians (including early tetrapods), which underwent sustained high rates of morphological evolution for ~30 Myr during the deployment of the tetrapod body plan.

Exposure to the Insecticide Sulfoxaflor Affects Behaviour and Biomarkers Responses of Carcinus maenas (Crustacea: Decapoda).

Sulfoxaflor is an insecticide belonging to the recent sulfoximine class, acting as a nicotinic acetylcholine receptor (nAChRs) agonist. There are few studies regarding sulfoxaflor's toxicity to non-target organisms. The present study aimed to investigate the acute and sub-lethal effects of sulfoxaflor on Carcinus maenas by addressing survival, behaviour (feed intake and motricity), and neuromuscular, detoxification and oxidative stress, and energy metabolism biomarkers. Adult male green crabs were exposed to sulfoxaflor for 96 h and an LC50 of 2.88 mg L-1 was estimated. All biomarker endpoints were sampled after three (T3) and seven (T7) days of exposure and behavioural endpoints were addressed at T3 and day six (T6). Sulfoxaflor affected the feed intake and motricity of C. maenas at T6. From the integrated analysis of endpoints, with the increase in concentrations of sulfoxaflor, after seven days, one can notice a lower detoxification capacity (lower GST), higher LPO levels and effects on behaviour (higher motricity effects and lower feed intake). This integrated approach proved to be valuable in understanding the negative impacts of sulfoxaflor on green crabs, while contributing to the knowledge of this pesticide toxicity to non-target coastal invertebrates.