More stories

  • in

    Fin whale movement patterns from the Mediterranean sea reveal summer habitat use in the North Atlantic and potential feeding behavior

    AbstractUnderstanding the movements of highly mobile marine species is essential for identifying critical habitats and developing effective conservation measures. Although occasional movements of fin whales (Balaenoptera physalus) between the Mediterranean Sea and the North Atlantic Ocean have been reported, the extent to which these movements represent a recurrent component of Mediterranean fin whale ecology remains unclear. Here, we present satellite telemetry tracks from seven fin whales (2022–2025) with five years of supporting evidence from systematic land-based observations (2021–2025) in the western Ibiza Channel to investigate movement patterns, behavior, and habitat use beyond the North-Western Mediterranean Sea. All tracked whales left the Mediterranean Sea passing through the Strait of Gibraltar and subsequently moved northward along the Atlantic coast of the Iberian Peninsula, providing direct evidence of cross-basin movements from the Mediterranean to the North Atlantic. Behavioral analyses using two independent movement-modeling approaches identified a common pattern of directed transit followed by area-restricted search (ARS) behavior—commonly associated with foraging—in several Atlantic areas, including waters off the Iberian Peninsula and the Bay of Biscay. The occurrence of prolonged ARS behavior in these areas is consistent with the use of Atlantic feeding habitats. The similarity of movement trajectories across individuals, years, tagging locations, and data streams indicates that seasonal use of Atlantic waters may represent a regular component of the movement ecology of at least some fin whales in the Mediterranean Sea. Given that these movements span multiple jurisdictions and areas of intense maritime activity, conservation measures should consider both important feeding habitats and the migratory routes that connect them.

    Explore related subjects
    Discover the latest articles and news in related subjects.

    Ecology

    Ocean sciences

    AcknowledgementsThis work was supported by the CaboRorcual project, co-funded by the Universitat Politècnica de València (UPV) with the support of the Biodiversity Foundation of the Ministry for the Ecological Transition and the Demographic Challenge (MITECO) and project MysticMed with the support of the Biodiversity Foundation of MITECO, through co-funding from the European Regional Development Fund (ERDF). Marina del Portet supported the project by providing berthing facilities and access to its infrastructure. We are particularly grateful to WWF Protecting Whales and Dolphins Initiative and SPA/RAC who provided financial aid to support SP and VP on the field and to purchase satellite transmitters. The Malizia Explorer team kindly offered logistical support by hosting the tagging crew onboard the Malizia Explorer in 2025. Terra X also provided support during the 2025 fieldwork. The Tethys Telemetry Lab also supported data collection over the years. We are grateful to P. Sánchez, M. Giovanni, A. Almagro and A. Hernández for their support in fin whale detection and assistance in locating whales for satellite tagging during the 2023-2025 field seasons. We thank Associació EUCRANTE and the Ajuntament de Dénia for supporting systematic land-based surveys from the coastal observatory between 2021 and 2025, with special thanks to Meritxell Maimi for her dedication to sighting efforts and volunteer coordination. We also thank Toni Martínez (Servicio Ambiental, Marino y de Pesca, Ajuntament de Dénia) for providing opportunistic sighting records. We are very grateful to Tiago Marques and Giuseppe Notarbartolo di Sciara for reviewing an earlier version of the manuscript and greatly improving the quality of the paper. Nino Pierantonio provided significant support with the generation of the code for the interactive map in the Supplementary Information, publicly hosted on GitHub. The authors thank the Consejo Superior de Investigaciones Científicas (CSIC) for its agreement with Scientific Reports, as well as the Duke Global Publications Fund, which supported the open access publication of this article.FundingEJB and BFT received funding from the UPV project CABORORCUAL, El cabo de los rorcuales: estudio de la presencia y origen del rorcual común en la plataforma y talud del Cabo de la Nao-Canal de Ibiza, with the support of the Biodiversity Foundation of the Ministry for the Ecological Transition and the Demographic Challenge (MITECO). EJB and VG received funding from the Universitat Politècnica de València (UPV), for the projects: 1) MYSTICMED, Misticetos del Mediterráneo: distribución, migración y conservación del rorcual común (Balaenoptera physalus) en aguas españolas, with the support of the Biodiversity Foundation (MITECO), through co-funding from the European Regional Development Fund (ERDF) (Reference 171); and 2) MEGAMED, Conservación de Megafauna Marina (CE 20230598). SP and VP received funding from WWF Protecting Whales and Dolphins Initiative and SPA/RAC, while the Malizia Team and the Tethys Telemetry Lab provided in-kind support.Author informationAuthor notesViola Panigada and Blanca Feliu-Tena contributed equally to this work.Authors and AffiliationsDuke University Marine Laboratory, Nicholas School of the Environment, Duke University, NC, Beaufort, USAViola Panigada & Douglas P. NowacekPratt School of Engineering, Duke University, Durham, NC, USADouglas P. NowacekTethys Research Institute, Viale G.B. Gadio 2, , 20121, Milano, ItalyViola Panigada & Simone PanigadaCentro Oceanográfico de Murcia (COMU-IEO), CSIC, C/El Varadero 1, Lo Pagan, 30740, Murcia, SpainBlanca Feliu-TenaInstituto de Investigación para la Gestión Integrada de Zonas Costeras, Universitat Politècnica de València, C/ Paranimf 1, 46730, Gandia, SpainBlanca Feliu-Tena, Eduardo J. Belda & Pablo SantonjaEDMAKTUB Association, Caldes d’Estrac, SpainEduard Degollada & Beatriu TortInstituto Investigación de Ciencia y Tecnología Animal, Universitat Politècnica de València, Camino de Vera S/N, 46022, Valencia, SpainVictor GallegoAuthorsViola PanigadaView author publicationsSearch author on:PubMed Google ScholarBlanca Feliu-TenaView author publicationsSearch author on:PubMed Google ScholarEduardo J. BeldaView author publicationsSearch author on:PubMed Google ScholarEduard DegolladaView author publicationsSearch author on:PubMed Google ScholarVictor GallegoView author publicationsSearch author on:PubMed Google ScholarDouglas P. NowacekView author publicationsSearch author on:PubMed Google ScholarPablo SantonjaView author publicationsSearch author on:PubMed Google ScholarBeatriu TortView author publicationsSearch author on:PubMed Google ScholarSimone PanigadaView author publicationsSearch author on:PubMed Google ScholarCorresponding authorsCorrespondence to
    Viola Panigada or Blanca Feliu-Tena.Ethics declarations

    Competing interests
    The track of fin whale PTT 232686 was partially published in V. Panigada et al. (2024), where analyses were limited to movements and behavior within the North-Western Mediterranean Sea. The transit of this individual from the Mediterranean Sea to the North Atlantic Ocean, together with the associated analyses presented here, has not been published previously. The authors declare no competing interests.

    Additional informationPublisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary InformationBelow is the link to the electronic supplementary material.Supplementary Material 1 (download PDF )Rights and permissions
    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
    Reprints and permissionsAbout this articleCite this articlePanigada, V., Feliu-Tena, B., Belda, E.J. et al. Fin whale movement patterns from the Mediterranean sea reveal summer habitat use in the North Atlantic and potential feeding behavior.
    Sci Rep (2026). https://doi.org/10.1038/s41598-026-74862-0Download citationReceived: 07 June 2026Accepted: 01 October 2026Published: 08 October 2026DOI: https://doi.org/10.1038/s41598-026-74862-0Share this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy shareable link to clipboard
    Provided by the Springer Nature SharedIt content-sharing initiative
    Keywords
    Balaenoptera physalus
    MigrationArea-restricted searchBiologgingState-space modelConservation More

  • in

    Pathbreak: A Biodiversity-Food-Governance Game for experiential reimagination of social-environmental futures

    AbstractSocial dilemmas are central to addressing biodiversity loss, yet tools to make them tangible remain critically scarce. Pathbreak: A Biodiversity-Food-Governance Game provides an engaging and hands-on experience to explore social dilemmas between short- and long-term, individual and collective interests, alongside those pertaining to trust, elections, democracy, intersectionality, and power. The platform offers researchers and practitioners a unique arena to study and foster collective learning and reimagine alternative social-environmental futures.

    Explore related subjects
    Discover the latest articles and news in related subjects.

    Ecology

    Environmental social sciences

    Environmental studies

    Geography

    Social sciences

    Addressing biodiversity loss is a significant governance and social learning challengeThere is compelling evidence of unprecedented and alarming rates of biodiversity loss1,2 and the overall trends posing a critical threat for humans at local-to-global scales3,4. However, addressing biodiversity loss presents a significant governance challenge. Effective biodiversity governance necessitates collaboration and coordination among involved actors to manage shared resources to meet often-competing needs such as producing food and maintaining habitat for diverse species. Achieving this requires overcoming the absence of direct or immediately observable causal relationships between actions and outcomes. It is a particular challenge to prioritize biodiversity in decision making, especially in light of tangible and short-term benefits from so-called free ride, when one benefits from the resource system without contributing to its longer-term maintenance5,6,7. The Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES)1 identified human activities such as land use change, direct exploitation, and pollution as direct drivers of biodiversity loss. However, observing how, for example, change in land use leads to biodiversity loss is not easy due to delayed and geographically dispersed impacts4. Thus, both economic measures (e.g., internalizing externalities—unaccounted costs to other actors and the environment) and regulatory interventions (e.g., developing regulations that create incentives preventing externalities) have practical limitations in addressing biodiversity loss. The latest IPBES assessment2 recognizes this and calls for fundamental shifts in views, practices, and structures to address the alarming rates of biodiversity loss and particularly accelerate food systems transformations8.Moreover, because “biodiversity” encompasses the totality and plurality of all living beings, the term entails particularly diverse perceptions, complicating consensus on definitions, issues, and solutions9. Under these conditions, social dilemmas, where interests of specific actors clash with interests of the broader system, have a tendency of creating even more room for free ride as preventing, detecting, and compensating for it is extremely difficult. Further, looking at available choices to any individual actor from a historical-institutional perspective10,11,12, it is evident that our actions are bounded by path dependencies (existing norms, infrastructure and policies, networks, etc.) that make it even more difficult to break out of established paths and make new choices. For example, if actors (individuals, communities, societies) have unsustainable consumption patterns with impact on biodiversity in other parts of the globe, breaking path dependencies (existing norms, infrastructure and policies, networks, etc.) and setting up a new governance system that takes into account the costs and either prevents or compensates free ride, especially in a centralized manner and top-down, is indeed a tremendous challenge. Hence, internalization of norms where actors develop their own intrinsic motivation compatible with preventing actions that lead to broader undesirable outcomes13,14,15 becomes particularly relevant for facilitating change in such complex systems.Games are known as a promising participatory intervention for social and experiential learning about complex social-ecological processes and internalization of norms16,17. Well-designed and implemented games can start transformations from within or bottom-up, that is individual actors might adjust their beliefs or develop new beliefs about issues at hand, which then can affect other actions and interactions they are involved in. A number of games that particularly focus on learning about social dilemmas and prevention of free ride have emerged recently, where relatively well-studied social dilemmas (e.g., public goods or common pool resource dilemmas) are framed in specific social-ecological contexts13,18,19. Generally, many studies show that context matters when generic games are stylized or framed for a certain social-ecological context, or generally when games focus on facilitating learning situated within specific time, place, situations, and in relation to specific people20,21,22,23. Likewise, developing a game for facilitating social learning about biodiversity is both a theoretical and practical challenge, as decisions must be made on myriad features of a game that can influence its outcomes both within and beyond the game, ranging from those on game narrative and experience to those on game rules and importance of player attributes19,22. Games that effectively facilitate learning about complex subjects, particularly by addressing the multiple layers of underlying social dilemmas to overcome the free ride problem, are relatively rare in the world of experiential games17,24,25,26 and virtually non-existent at the intersection of biodiversity, food, and governance.Social dilemmas at the heart of addressing biodiversity lossHere, we aim to fill this gap by sharing our lessons of developing and piloting Pathbreak: A Biodiversity-Food-Governance Game (www.pathbreak.eu). We present Pathbreak as a proof-of-concept boundary object that operationalizes layered biodiversity–food–governance dilemmas and elicits participant reflection through facilitated debriefing. We focus on the importance of integrating multiple direct and indirect layers of dilemmas for social learning. We reflect on to what extent the new game succeeds in bringing these dilemmas to the attention of the participants using an in-depth analysis of debriefing sessions from five pilots implemented in a three-year Horizon Europe project dedicated to societal dimensions of addressing biodiversity loss. Our specific research questions are: How does Pathbreak operationalize layered biodiversity–food–governance dilemmas? What kinds of first-order and second-order dilemmas emerge in participants’ debriefing reflections? How do participants connect the game experience to real-world biodiversity governance challenges?A social dilemma is any group situation where one actor is, at least temporarily, better off when this actor does not cooperate while others in the group do; however, if all actors defect, the payoff is substantially lower than if they cooperate27,28. As this structure of incentives generates tension between individual and group rationality, such problems were widely analyzed by researchers. The investigations included not only the economic framework of the dilemma29,30,31 but also, among others, the role of monitoring, sanctioning, and institutions32,33,34,35, communication36,37,38,39, individual values, mental models, and social norms40,41,42,43, as well as group identity44,45,46 that moderate the cooperation and competition in such situations.Research shows that experiential learning through social dilemmas can be particularly beneficial for learning about social processes when facing dilemmas, especially in social-ecological contexts. Bush et al.24 highlight that “a degree of firsthand discovery and novelty (e.g., violation of basic assumptions)” when players experience the dilemma situations and the consequences of their decisions may facilitate conceptual understanding of complex phenomena. Bush et al. further note that this is in line with studies on both decision-making in naturalistic settings47,48,49 and educational theory50,51. Generally, ecological issues are deeply rooted in social dilemmas, where individuals are inclined or forced to prioritize personal interests over societal consequences52. Despite being the core drivers of major social-ecological challenges worldwide5,6, this connection often goes unrecognized by many53 and also within the latest biodiversity assessments1,2. However, understanding this fundamental link between social dilemmas and environmental degradation is essential for fostering effective societal cooperation and shaping public policies that promote human well-being and global sustainability5,54. Thus, there is a strong need for enhancing public awareness of social-ecological systems and their inherent dilemmas55,56.In the biodiversity domain the need for social learning games is well documented57 especially for ones based on social dilemmas, as well as social-political processes around them. In contextualizing generic social dilemma games (e.g., prisoner dilemma, public goods, common pool resources games), a caution is necessary about designing choices players have to make in the games. The complexity often arises from social interactions and group processes, beyond the underlying formal model of the game that might focus on decisions between sets of limited choices. As Levit and List20 put it, “in experiments, the researcher creates a set of rules governing the interactions, chooses the wording of instructions, and defines the set of actions the subject is allowed to take. In stark contrast, in naturally occurring environments, the choice set often is almost limitless and institutions arise endogenously”.Based on these theoretical considerations and building on what Ostrom5,6 described as first-order and second-order social dilemmas, within the domain of biodiversity-food-governance, we identify at least the following three sets of factors that define the key social dilemmas central to addressing biodiversity loss.Direct or first-order biodiversity-food dilemmaAlthough this dilemma is not linear, biodiversity loss is frequently associated with agricultural expansion or intensification, which typically requires more land or intensifies existing land use58,59,60. Biodiversity loss, in turn, negatively impacts productivity of ecosystems, affecting the overall resilience of social-ecological systems—their ability to respond to changing needs and future uncertainties. Increasing numbers of extreme weather events in the last decades are often closely associated with the declining diversity and health of the ecosystems, habitats, and species, which otherwise provide often-vital ecosystem services and perform functions to stabilize climate and to ensure longer-term sustainability of (non-human) nature’s contributions to people61,62,63. At the same time, some literature emphasizes the unavoidable need for high yields in the context of global food shortages due to the detrimental consequences of reduced yields, both for the impoverished communities and countries and in terms of risks associated with further agricultural expansion to meet newly generated import demands64. This first-order dilemma presents producers with a choice between lower and higher yields (through expansion, intensification, and their mixtures), while consumers face the questions of food affordability and/or make choices between higher or lower levels of consumption. Each decision must balance the individual and shorter-term payoff with the group and longer-term payoff.Indirect or second-order governance and political dilemmasThe governance and political dilemmas focus on how much energy should be invested in communicating, developing new rules, and generally being active beyond one’s own private sphere with an attempt to influence others in a given community. Here, depending on the individual circumstances, the consequences of the dilemma can be substantially different, therefore the term dilemmas is used in plural. It is well known that effective governance in complex social-ecological systems that can address both the equitability and sustainability considerations in the system and ensure participants avoid the worst-case scenarios, requires effective coordination52. Individuals in positions of power often invest significant effort to attain their positions and derive privileges from their authority, including pursuing policies that serve public and specific interests. They also endure public scrutiny and bear the risks associated with power misuse or changes in power dynamics, including retaliation. One can also be politically active without getting in a formal position of power, through mobilizing, rallying behind a certain cause or candidate, expressing or withholding endorsement, etc65. All of these options also have potential consequences consisting of risks and benefits, including the prospect of being on the losing side, enduring disappointing policies enacted by the supported candidate, being rejected in communities that supported the other candidate(s), and so on. Further, the coordination of providing public goods and maintaining common pool resource systems is often costly in larger scales66. When decisions have to be made that change how society prioritizes various objectives—for example, limiting or expanding available land for agriculture, raising or increasing tax for use of public goods or common pool resources, dedicating resources to various issues—one can think of the immediate decision maker, for instance, a governor, or entire society where historical events led to specific norms and values in that society and specific leaders being elected and setting specific political agenda. Considering the state of biodiversity and food as a result of such broader historical social-political developments, it becomes evident that any transformative change in the biodiversity and food domains would necessitate transformative change in the governance domain.Place-based intersectionality dimensionsIntersectionality dimensions such as age, gender, education, income, cultural background, religion, disability combined with specific place are known to make opportunities and barriers asymmetric, explain systematically different perceptions of the same phenomenon (e.g., the dilemmas above), and result in fundamentally different individual interactions, sometimes unlike what is predicted by formal theory or models that report mean values or “typical” behavior67,68,69. Thus, contextualization of the game for addressing biodiversity loss might have some general dilemmas, for example of deciding based on individual versus group interest, yet for a very specific person with very specific background in a specific place these dilemmas can mean substantially different trade-offs. For some a public goods game might mean an exercise in economics, while others might be imagining a very specific school that is being constructed in their neighborhood. For some it could mean selecting food over biodiversity for livelihood and avoiding extreme poverty already in times of hardship, for others it could mean selecting food over biodiversity as a matter of maximizing profit and preserving a long-term agricultural family tradition. For some, speaking out might mean an integral element in the battle for power and nurturing democracy, for others it might mean risking their immigration status or attracting personal attacks while already being disadvantaged.Looking at these dilemmas and factors, it becomes evident that addressing the indirect drivers and underlying causes of biodiversity loss goes well beyond the biodiversity domain. Experiential learning is considered particularly useful for such complex and multiple layers of dilemmas. According to David Kolb, experiential learning is “the process whereby knowledge is created through the transformation of experience. […] Knowledge results from the combination of grasping and transforming experience”70. One should note that Kolb uses the term “transforming experience”—that is when experience turns into knowledge through reflection and realization of what was experienced. Kolb70 suggests a four-stage experiential learning model where knowledge is continuously created through concrete experience, reflective observation, abstract conceptualization, and active experimentation. In using games for facilitating learning, active testing in the game is often linked to direct experience, reflected upon, and conceptualized in debriefings, which scholars consider as the most important point of learning from the use of games and simulations71,72,73,74. One can distinguish cognitive (focus on knowing) and non-cognitive learning such as affective (focus on emotional understanding of others and the situation), motivational, relational, normative, sociocultural (focus on understanding or building relationships and trust, norms of interaction, etc.)19,75,76. Overall, both David Kolb as the pioneer in experiential learning theory and many who have worked on it over the last decades explain that experiential learning is a transformative process where one readjusts one’s own knowledge and relational system to integrate the new experience.Thus, transformation in learning for us has two important dimensions—it is a process that effectively transforms the simulated experience into actionable knowledge70, and one that in a complex social-environmental context such as in addressing biodiversity loss facilitates learning that goes beyond any single level of social change2,9. Our experiences of developing and piloting the Pathbreak game in PLANET4B (see background to the project in the “Methods” section) confirm these multi-faceted societal dimensions of biodiversity on the one hand and games and learning on the other.Pathbreak experiential learning platform: Integrating layered dilemmasBefore introducing the Pathbreak framework and how the layered dilemmas are integrated, we should note that the game elements are simultaneously both the result of “translating” the theory of social dilemmas discussed earlier into a tangible and practical framework that allows to experience these social dilemmas and the platform for understanding further insights from experiencing it.The conceptPathbreak: A Biodiversity-Food-Governance Game immerses participants in a simulated community where they can explore and actively test70 the complex interconnections between biodiversity, agricultural practices, and policy-making. The game is called Pathbreak as the participants will face the challenge of breaking the path consisting of unsustainable production and consumption patterns and requiring reconciliation of one’s own, various sub-groups’, and broader system’s interests (resilience of the social-ecological system). The activity runs for approximately 90–120 min, including test rounds, the main gameplay, and a debriefing session. It is facilitated through a combination of physical game boards, role cards, choice cards, and an online Web Application for data input and estimation of points (see Fig. 1, the detailed game materials and facilitator’s manual available at https://doi.org/10.5281/zenodo.17476375).Fig. 1Full size imageA recommended table layout in Pathbreak: A Biodiversity-Food-Governance Game.The interactive setup and embedding the key dilemma in goal-settingIn the game, one facilitator and seven players—three producers (representing farms) and four consumers (representing households)—sit around a table and make a number of decisions (Fig. 1) (more than one participant can play as one player). At the beginning participants are informed that they can set their own goal for the game. They are also informed that if by the end of the game the state of biodiversity is not in red, one person with the highest points will receive a symbolic reward (e.g., sticker, origami, pen).Structured and unstructured decisionsThey make decisions that have consequences for their individual livelihood (individual points) and the state of biodiversity. They can communicate and make decisions about governance that can result in policy change (virtually limitless set of decisions except for the limited time, such as to what extent to be politically active and how).Selected policiesThere are two specific policies in the game. Land Policy applies to producers and can reduce or increase available land for production. The default policy in the first round is high land availability for production. Tax Policy applies to consumers and can reduce or increase net income for spending. The default policy is low tax/high disposable net income.Direct first order social dilemmaIn each round players representing farmers decide on the intensity of production (high or low). Consumers can choose between high and low consumption (Fig. 2). They receive higher or lower individual points accordingly. Points are estimated through a market-based supply and demand model (see “Methods”) that is made available to the participants live via the Web Application.Fig. 2: Two examples of what the main game board could look like after placing the cards on the board.Full size imageA Policies are high land availability for producers and low taxes/more net income available to spend for households; and B Policies are low land availability for producers and low taxes/more net income available to spend for households. Policy cards are framed as “Tax & Investments” and “Rewilding Area: Diverse Trees, Bushes and Flowers”.Impact on and from biodiversityThere are three states of biodiversity: high (green), reduced (yellow), and low (red) (Fig. 3). Any combination of decisions except one will worsen the state of biodiversity or keep it at low if the players have already reached it. The exception is when the Land Policy is low land availability and all of the three farms produce with low intensity. Two rounds of coordinated action are needed to improve the state of biodiversity and one round of uncoordinated action for worsening it (or keeping it at low/red).Fig. 3Full size imageThe three states of biodiversity.The state of biodiversity has an impact on the likelihood of events that can take place after each round (based on the works revealing the importance of biodiversity for increasing the resistance of ecosystem productivity to climate extremes, see earlier discussion). Events—extreme weather events—penalize all players by deducting substantial points (−80, see the Methods section for calibration details). When the state of biodiversity is high, the likelihood of events is 0%, when reduced 30%, and when low 80%. Players draw from a stack of cards after each round with the color corresponding to the state of biodiversity in that round (yellow or red).Indirect second order social dilemmasBetween the rounds, there are community council meetings where players can communicate about anything they would like.Then players make decisions about the Governor. First, they are asked whether anyone would like to nominate themselves or another player to become the Governor who can then change the policies. Nominees have the opportunity to make a brief statement about why others should vote for them. Once the statements have been made, the players vote for candidates. In case of a tie, if there are more than two candidates, the second round of voting takes place only with the two frontrunners. The frontrunners can make new statements. If there is still a tie, the Governor’s office remains unoccupied and policies remain unchanged.If there is a Governor, the Governor makes decisions on two policies: Land Policy and Tax Policy—whether to keep or change each of these.Each elected Governor can serve for the maximum of three terms (three rounds in the game). Once there is a sitting Governor, they can be removed from the office sooner through a referendum. When there is a sitting Governor, the facilitator asks players whether there is a need for a referendum. If the majority votes for referendum, the facilitator asks players to vote for keeping or removing the current Governor before proceeding with the possible new election.A structured intersectionality dimensionOne of the consumer players is a migrant without citizenship and therefore is not eligible to vote in the referendums and elections, nor to run as a candidate. The facilitator informs players about this during the first round of the game after asking who would like to run for the office of the Governor. However, the facilitator reveals in round five that the migrant’s application for citizenship has just been approved and now this player can participate in referendums and elections with the same rights and responsibilities as the others.An unstructured intersectionality dimensionA further (optional) feature is the role cards. After the initial two test rounds and before the actual game begins, each player receives a role card with a short profile. These cards are distributed randomly. They contain information on their more general role (consumer 1, consumer 2, consumer 3, consumer 4, producer 1, producer 2, producer 3) and profile. The profile of each player was generated to make sure that some of the key arguments from real-life social dilemmas (link to direct experience70) relevant to the broader topics of the game are available to the players (livelihood and personal background, views towards biodiversity and agriculture, the role of government and policies). However, the facilitator informs that this information is optional and confidential, and it is up to the players to what extent to share and build on the provided profile information. This is done to ensure that the players can still actively bring in their own experiences and thoughts, not necessarily being limited by the pre-defined profile.Sanctions to sharpen second-order dilemmas during the gameOne further option that becomes available to the Governor later in the game is sanctioning any given one player alongside the decisions on the two policies. If the state of biodiversity is low at the beginning of the round and remains low after players make their decisions in the current round, the Facilitator informs the players that the Governor has emergency powers to sanction someone. If the Governor uses these powers, the sanctioned player loses points (30 points).Number of roundsThe game continues for seven rounds (based on the modeled calibration provided in the Methods section). Although it is apparent from the game materials that there will be seven rounds, it is not announced by the Facilitator, to avoid making it salient and nudging players to maximize profit in the last round. This information is neither concealed, and the facilitator can inform about the number of rounds, should the players ask.Debriefing and reflectionsAfter the game there is debriefing that focuses on reflections on the game experience, potential real-world parallels, and conceptualization of the experiential learning as a whole – thus completing the Kolb’s70 learning cycle.The overall gameflow for each round is depicted in Fig. 4 and also placed on the table for the players as an orientation.Fig. 4Full size imageGameflow for each round in the game.Pathbreak and experiential learning cycleThus, pedagogically the game is structurally aligned with Kolb’s experiential learning cycle. The gameplay as a whole serves as the Concrete Experience, immersing participants in the visceral tension of managing the first order and second order dilemmas. Because the game is played in repeated rounds, it naturally embeds Active Experimentation directly into the gameplay loop. Participants continuously test and adapt their strategies through both structured experimentation—such as deploying choice cards, voting, policy decisions, and use of sanctions—and less structured social experimentation, such as adjusting their level of engagement during community council discussions, navigating elections, and deciding whether to actually follow through on stated intentions from previous rounds and potentially face corresponding consequences. Following the gameplay, the debriefing protocol shifts the focus to Reflective Observation, providing a dedicated space for players to articulate their feelings (positive, negative, neutral), review their in-game actions, and critically assess the social dilemmas they faced. Finally, the debriefing guides participants toward Abstract Conceptualization, allowing them to map their localized gameplay experiences and iterative strategies onto complexities of real-world biodiversity governance and broader theoretical thoughts on how society can or should address challenges such as biodiversity loss.Real-world complexity and simplification in PathbreakThe game is intentionally designed around a core dilemma: the trade-off between short-term gains from high production and consumption and long-term sustainability linked to biodiversity. To make this dilemma salient and experientially accessible, several aspects of real-world systems are deliberately simplified. First, ecological dynamics are represented through an aggregated biodiversity indicator, rather than multiple species, habitats, or ecosystem services. Second, players are given a limited set of policy and decision options, each with relatively strong and transparent effects on system outcomes. Third, actor roles and socio-economic diversity are simplified, focusing on a small number of asymmetric roles rather than a full range of stakeholders and institutions. By reducing the number of variables and decision pathways, the game provides clear and immediate feedback, thereby shortening the experiential learning loop (in the sense of Kolb). Within the limited number of rounds, participants can observe how repeated choices lead to system degradation, attempt to reverse these trends through cooperation, and experience how such recovery efforts may be undermined by individual incentives, role asymmetries, and probabilistic disasters.Despite its abstraction, the game is designed to reflect several well-established real-world dynamics in biodiversity–food–governance systems. At its core, the model captures a “business-as-usual” trajectory (e.g., refs. 1,2), in which increasing production and consumption generate short-term benefits while progressively eroding biodiversity. This reflects widely documented trade-offs in agricultural intensification. The game operationalizes this through a feedback structure in which declining biodiversity increases the probability of adverse events, representing real-world risks such as ecosystem service loss, climate-related shocks, or system instability. The game also embeds multiple layers of social dilemmas. Sustainable outcomes require coordinated efforts across actors, as unilateral actions (e.g., one farmer reducing production) are insufficient to shift the system trajectory. This mirrors real-world challenges of environmental governance, where individual short-term incentives often conflict with collective long-term benefits. Additionally, the model includes elements of social dynamics, such as asymmetries in roles and influence over decision-making processes, as well as the possibility that individual “defection” can undermine both ecological outcomes and social trust. Together, these mechanisms correspond to key real-world phenomena, including externalities in agricultural systems, tragedy-of-the-commons dynamics (collapse in the absence of communication and coordinated collective action), delayed ecological feedback, and the overall importance of governance structures in managing coupled human–natural systems.We should note that pathways such as agroecology, regenerative agriculture, or biodiversity-compatible intensification are not explicitly modeled in the current game design. This is an intentional design choice rather than an omission. Incorporating these approaches would require representing multiple alternative production systems, heterogeneous landscapes, and more complex policy instruments, which would substantially increase the number of rules, decision points, and possible system trajectories. Such an expansion would risk obscuring the core learning objective and reducing the clarity of feedback within the limited time frame of the activity. The current design prioritizes portability, ease of facilitation, and the ability to complete multiple rounds, allowing participants to observe system dynamics and potentially shift trajectories within a single session. Importantly, the absence of these pathways in the formal model does not mean they are excluded from the learning process. On the contrary, the game is designed to create a conceptual space in which the need for such approaches becomes apparent. During the debriefing phase, participants are encouraged to reflect on the model’s simplifications. In this way, the game functions as a foundational abstraction and a boundary object, while the broader activity integrates more complex and context-specific knowledge and experience, shared by the participants themselves.Proof of concept: results from piloting and debriefingExperiencing the direct biodiversity-food dilemma in PathbreakDebriefing analyses demonstrate that the game’s architecture operationalizes direct social dilemmas as intended. Participants consistently identified the structural friction between individual profit-seeking (via intensified production or consumption) and the collective objective of sustaining biodiversity to mitigate extreme weather events. The complex relationship between biodiversity and food was the most frequently discussed theme in the debriefings (Fig. 5). It seemed also clear that participants experienced several layers of individual versus group dilemma—in the sense that there were nominally two groups of producers and consumers initially, but during the game players would realize that there were often points of clash and cooperation also within various sub-groups. They reflected on the dilemmas in personal terms, that is “I and others”, as well as highlighting broader clashes in “individual and group” or “short-term and long-term” interests.Fig. 5Full size imageFrequency and distribution of themes emerging during the post-game debriefings.The way participants described their experiences also allows us to distinguish roughly between two types of social learning: cognitive and non-cognitive. Cognitive learning about the dilemma manifested itself as participants tried to make sense of the dilemma, for example by trying to calculate or figure out how the economic model behind the game would predetermine the game outcomes. Non-cognitive learning about the dilemma manifested itself in the way the participants expressed their frustrations by emphasizing the feelings, emotions, the role of trust or lack of it, and how they felt they could or could not identify with the group or the others, highlighting the behavioral dimensions of economic decision-making77.There were also a number of valuable critical perspectives voiced by the participants. One of the key critical points was related to how the game’s design forces players to choose between profit (individual points) and biodiversity, at least implicitly neglecting the possibility of achieving both, which might in turn reinforce the image of the current unsustainable consumption and production patterns. This is an important critique for at least two reasons. First, it is true, as noted earlier, there are non-linear relationships between biodiversity and food, and these should not be neglected. However, that this discussion emerged serves the purpose of Pathbreak—to experience prevailing unsustainable paradigms and talk about possible ways out. Second, the participant’s perception that the game forces players to choose between profit and biodiversity is rather a disappointing news for the developers as participants in the game face myriad invisible choices well beyond biodiversity and food: to be active or not, to communicate or not and to what extent, to run for the political office or not, to enact various policies or not, to support specific candidates or not, to bring specific issues to the attention of the community council or not—the list goes on. Of course, the fact that the participant did not notice that there were also these non-biodiversity and non-formal choices at all times might signal the game’s weakness in facilitating learning. Although, one could argue it is a good simulation of contextual reality where actors might not realize the presence of these rather invisible higher-level political-institutional choices11 (see more on this next). Further, while at least one participant perceived the game extraordinarily difficult, particularly in terms of how difficult it was to recover from extreme weather events, another participant found the game limited in its parameters. However, both of these observations are in line with the intended design of the game.Uncovering the indirect governance and political dilemmas in PathbreakThe pilots and debriefings showed that the participants questioned dilemmas they encountered beyond biodiversity and food. They discussed the role of communication, building relationships, intersectionality—privilege or lack of it originating from one’s background, power of the market, difference between being a consumer and an active citizen, and the value and limitations of elections and democracy. They also talked about the complexity of how policies are made, and some examples of how changes made by a new government, even if beneficial in some sense, might not be immediately apparent due to the time required for implementation and the legacy of previous administrations. Hence, they also discussed how political campaigns did not necessarily focus on what could be actually delivered. One participant formulated it: “…at the end of the day, it was more about agreeing between the people. I mean, it showed that that improvement comes from agreement [process], in between the people, among the people, that power” (Debriefing 3, lines 1672–1674, p. 48).While reflecting on these, often-underestimated, rather abstract elements of the common-pool resource games, the debriefings also showed that the participants made real-life references across all themes (Fig. 6).Fig. 6Full size imageFrequency and distribution of themes by references to real-life and game mechanics.The participants highlighted the game’s potential for improving policy-making. This would be in line with the literature that highlights the value of games for facilitating stakeholder engagement78. One participant stated: “I think it would be really useful for people making local land management decisions. I think where, so I live in an area of sort of lots of high conservation interest, there’s a lot of conflict between the local community groups, councils and people like the National Trust and NGOs and big landowners. And I think when they’re having their big decision making, two-day workshops, it might be nice to play that game right at the beginning. That’s a sort of war game thing. And then they go into their negotiations maybe afterwards”. (Debriefing 1, lines 530–536, p. 16).There were also those who were connecting their game experience with real life, while, to some extent ironically, being frustrated that the game was not realistic enough: “There were some elements that, uh, frustrates me, which I struggle to connect with, let’s say, real life. Other than the fun, and joy that comes with it. There’s no question of this. But, for example, I question myself, why do we start from green? Uh, because it’s frustrating that actually in real life we’re not, realistically speaking, in a green area, right? It feels like we are in a red area somewhere. Um, and of course it does, gives the, all the way back, gives the impression of a long road ahead to reach the green space, so to speak. But still, the frustration comes when we’re not there. We’re told, all the way, red and so on. So, these elements give me some frustration to connect with real life. Not, again, not touching upon the fun part and then everything else that comes with it, which is also great”. (Debriefing 1, lines 33–41, p. 3).Intersectionality and identity in governanceA defining feature of the Pathbreak simulation is the deliberate introduction of power asymmetries through role cards. In particular, the inclusion of a consumer role representing a migrant without voting rights served as a primary driver for discussions regarding political legitimacy and environmental justice. In several pilot sessions, this role acted as a “governance bottleneck”, where the excluded player could influence community dialog but remained legally disenfranchised during formal referendums and elections.Participants reacted to this legal barrier with a mix of frustration and strategic adaptation. The eventual granting of citizenship and voting rights in the fifth round often functioned as a critical turning point, shifting the community’s focus from mere economic survival to the broader ethics of democratic participation in environmental decision-making. This transition illustrated that addressing biodiversity loss is not solely a matter of scientific or economic consensus, but is deeply intertwined with the inclusion of marginalized voices and the dismantling of institutional barriers.We have also observed that participants often started the game making references to the profile information and gradually merged it with the knowledge from their own backgrounds. In addition, while the sample of 50 participants across five pilot studies is informative for a proof-of-concept, the heterogeneity of the sample (e.g., experts, students, community members) raises important questions. Differences in prior knowledge, and “cultural capital” may have influenced participants’ engagement with the game and the resulting learning outcomes. One should be aware of the self-selection bias – all participants were already either open for experimenting with such topics and activities or embedded in a professional or broader social environment where trying a game like this was not something completely out of ordinary. It would be useful to test the game with audiences with little familiarity with the topic and across cultural contexts. Yet, the fact that all games revealed the multiple layers of the embedded dilemmas, despite the heterogeneity, gives a good indication that these are genuine and appropriate dilemmas of biodiversity-food-governance for consideration in multi-actor settings with diverse backgrounds.Lessons specific to tackling biodiversity lossBased on the analysis of the debriefings, the overall impression in terms of lessons specific to tackling biodiversity was related to the problem’s complexity and frustration coming from the realization of this complexity, which is in line with the literature that highlights games being particularly useful for facilitating learning about complex subjects79.“I think for me […] I use that word frightening. […] it made me think about the bigger picture and […] again, not so much for me, but […] about the next generation […] and how or what we can potentially do to influence things going forward. So, if nothing else, the game certainly made me think more about what I can do going forward and maybe a desire to do more quickly”. (Debriefing 1, lines: 695–700, p. 21).Some participants shared thoughts that indicated they were still hopeful, despite the complexity and frustration, what seemed to be rather an endorsement for the game and relatively deep discussion they had just experienced.“Yeah. So, then I think it creates […], kind of sets the ground for more ideas, like just a brainstorming session and discussion session. Just information sharing, especially with the […] group. That would be nice to get kids from different places to exchange on a topic that’s common to all of us. So that’s really nice, I think, about this game”. (Debriefing 1, lines 744–747, p. 21).Looking at these two illustrative quotes, one can also note that the participants thought the game could be useful for different age groups in terms of both making one think about future generations and helping younger generations to learn about layered dimensions of the topic.Learning across contexts with different emphasisWe have also witnessed some similarities and differences across the debriefings. Figure 7 shows dominance of one or another group of themes. The game and the debriefing generally led to a frame that consists of social dilemmas, political processes, complexity of biodiversity. The figure also shows that it still allows participants to develop their emphasis which can vary considerably. While experts tend to talk about the game as a tool and provide more feedback on its design, students talk more about politics and governance, social dilemmas, and specific biodiversity lessons, which is in line with their generally active learning phase of life. The learning community members, some of whom played the game twice, talked more about the key social dilemma. The external group focused on the key social dilemma most, then politics and governance of it, and then lessons specific to tackling biodiversity loss. It should also be noted here that while the transcripts from the debriefings with the experts were more challenging to code as they seem to be in constant search of topics at higher level and therefore not necessarily address themes directly, the students’ group was easiest to code and analyze. We also see that some of the discussions were very detailed and free-flowing, while others were shorter and more focused, once again indicating that the design of the game had both the structure to talk about specific dilemmas and flexibility to bring in the perspectives of the participants.Fig. 7: Debriefing focus diagrams for five pilot studies illustrating thematic composition and dynamics.Full size imageEach figure represents the changing themes over time, where the start of the debriefing is at the top left and the end is at the bottom right.DiscussionAddressing biodiversity loss requires grappling with complex layers of social dilemmas that link individual action, collective outcomes, and institutional arrangements. Pathbreak: A Biodiversity–Food–Governance Game addresses this challenge by providing an experiential setting where participants confront the tensions between short-term individual gains and long-term collective well-being. We acknowledge that the presented work has been piloted in a small number of cases, yet we argue that it serves as a proof of concept, and its value lies in the framework and the qualitative depth of the boundary object, rather than statistical generalizability. In line with Nobel Laureate Elinor Ostrom’s seminal works on governing the commons, our findings show that participants frequently experienced and debated both first-order dilemmas (production/consumption vs. biodiversity) and second-order dilemmas (governance, trust, and coordination). These multilayered dilemmas are rarely addressed simultaneously in biodiversity education and learning tools, yet they lie at the heart of biodiversity decline1,2.A central lesson from the pilots is that players did not merely interpret the game in ecological or economic terms, but also engaged with broader political and normative questions: power asymmetries, governance legitimacy, and the difficulty of sustaining collective agreements. This echoes Ostrom’s insight that overcoming social dilemmas depends on collectively crafted rules, communication, and the internalization of reciprocity and trust5. Several participants highlighted that “improvement comes from agreement between the people”, resonating with Ostrom’s emphasis on social capital as a condition for cooperation capable of overcoming dilemmas. That the game sparked similar reflections across diverse groups—from students to practitioners and experts—suggests that it has a broad potential to support social learning and dialog on biodiversity-food-governance dilemmas.At the same time, our pilots highlighted the importance of being aware of critical design limitations when using experiential games: some simplification or abstraction is necessary, but one should avoid oversimplification. Our observations and experiences show that the simplified ecological dynamics and the associated hysteresis effect embedded in the game, and as observed in many ecosystems80, have triggered reactions similar to what one would expect in the real world. Systems that pass a tipping point do not necessarily recover when pressure is released. That takes additional effort and/or time. Our experiences in the pilots are in line with the previous collective action games with tipping points21,81,82,83,84,85,86,87,88,89,90. However, most other studies do not include events happening dependent on the state of the ecosystem (but see Janssen et al.66), and Pathbreak extends the current literature by offering a social-dilemmas game on biodiversity with a dynamic ecological feedback mechanism. The ecological dynamics in the game is simple compared to real ecosystems also because the goal is not to create a model with similar complexity as reality. As we have seen in our pilots as well, some participants might find it challenging to understand the dynamic processes even within this simplified environment. For participants, these simple dynamics are typically a challenge to comprehend within a short amount of time. In fact, research has shown that humans, including graduate students at MIT, have difficulty managing dynamic systems91,92, and some scholars have proposed that the tragedy of the commons might be caused by people’s inability to understand dynamical system instead of selfish behavior93. To avoid that we measure confusion instead of intended decisions, a simplified ecological model that captures relevant dynamics for our purposes is therefore a necessity. This is still without adding the dynamics stemming from the multiple layers of social dilemmas. The pilots demonstrate that the game functions as a sandbox for social interaction which enriches the existing simplicity in the ecological system with social complexity. In contrast to other studies where all players have the same role, this study demonstrates how integrating asymmetric roles, power imbalances, economic inequality, and intersectionality dimensions such as migration in a stylized form can generate virtually limitless scenarios of cooperation and conflict that can fundamentally re-shape how players address the challenges of the underlying dynamic ecological model.Finally, our pilots show the importance of seeing the game not as an end in itself but as a part of a broader activity. Its primary purpose should be to provide a shared experience that informs subsequent discussion and debriefing. In line with Kolb’s experiential learning cycle, our pilots confirm that the debriefing phase enables participants to connect the model to real-world contexts, reflect on its limitations, and engage in learning and exploration of often-limitless scenarios. This stresses the importance of careful debriefing where perceived limitations can be critically discussed in depth. We should further note, given that the game developers also facilitated the sessions and conducted the primary analysis, recognizing and mitigating insider bias was a methodological priority for this study. The developers’ deep familiarity with the game’s mechanics was essential for accurately mapping participant actions to the intended social dilemmas and keeping in mind that the game experience should be addressed in the subsequent debriefing. Several additional safeguards were implemented to prevent this positionality from unduly influencing the results. During data collection, facilitators minimized their interference by adhering to standardized scripts for the introduction and rules explanation, maintaining a pre-defined and observational role during actual gameplay. To counter the inherent power dynamic during the subsequent debriefings, facilitators explicitly encouraged feedback, including critical concerns and recommendations for future improvements, empowering participants as expert evaluators to disrupt the tendency for polite, confirmatory responses. Furthermore, to ensure analytical distance, a triangulation strategy was employed during the coding process: the coding team always included a researcher who was not part of the pilot but acted as “a critical friend” during weekly meetings. This way we ensured that the primary coder’s interpretations are actively challenged, and that the final themes—especially regarding deviant cases where participants challenged the game’s premise—remained strictly grounded in the empirical data rather than the developers’ original design intentions. Thus, future iterations could explore game experience with facilitators not involved in game development. Similarly, one could imagine exploring more nuanced pathways where biodiversity-compatible production strategies are viable, but it would be important to preserve the clarity of the dilemma structure. Likewise, it should be acknowledged here that the presented study captures immediate post-game reactions. Further research is needed to generate evidence regarding the long-term retention of learning or its translation into real-world change in views and actions.Importantly, our findings from the pilots point to the contextual nature of learning. As path dependence and bounded rationality literature notes11,12, actors interpret dilemmas through the lens of their histories, institutions, and identities. In Pathbreak, participants drew on personal experiences of land conflicts, consumer choices, or political activism, illustrating that the same formal dilemma can carry different meanings across settings. This underscores once again the need to treat Pathbreak not as a “model of reality”, but as a boundary object—a tool for eliciting situated knowledge and facilitating dialog in multi-actor settings.By integrating direct biodiversity–food dilemmas with indirect governance and intersectionality dimensions, Pathbreak extends the repertoire of social dilemma games beyond classic public goods or common-pool resource framings. The design advances the agenda set by Elinor Ostrom and others to move beyond stylized economic models and to embrace the complexity of real-world governance of shared resources, environment, and place.Taken together, the pilots reveal that Pathbreak shows promise in fostering both cognitive and affective learning, prompting reflection on trade-offs, power, and pathways for collective action. Once again, while the present study is limited to five pilots, it builds on decades of expertise from academia, policy, and practice. The in-depth analysis demonstrates the game’s potential to contribute to biodiversity-food-governance debates by making abstract dilemmas tangible. Future research could test the game in decision-making contexts with practitioners and policymakers with wide-range of backgrounds, and explore evaluation of learning outcomes across cultural contexts and social groups.To move from experiential learning to systemic impact, we propose a “Theory of Change” where Pathbreak can potentially serve as a catalyst for better informed policy negotiation. Our findings suggest the game could be beneficial not only for students and experts studying social-ecological challenges, but also as a precursor to formal stakeholder workshops as it creates points of references to key social dilemmas behind addressing a complex social-ecological challenge. By providing a low-risk environment to test the social legitimacy of interventions and multiple scenarios in repeated rounds, Pathbreak could offer a practical tool for stakeholders tasked with developing and implementing international goals, such as the targets set in the Kunming-Montreal Global Biodiversity Framework or to be developed for post-2030 Sustainable Development Goals. The game’s ability to expose “invisible” institutional strengths and barriers makes it potentially applicable for navigating the localized land-use conflicts inherent in scaling up protected and rewilded areas. By immersing actors in the layered dilemmas of trust and reciprocity before negotiations begin, Pathbreak could potentially de-escalate “us-versus-them” mentalities common in land-use conflicts. Specifically, beyond research and education, we recommend testing the tool further by integrating it into:

    Local land and environmental management decisions: Facilitating dialog between community groups, municipal councils, and NGOs to reveal the “invisible” political-institutional choices that shape land use and environmental governance.

    Transdisciplinary science communication: Bridging the gap between STEM-focused biodiversity data and the social-political realities of implementation.

    Policy simulation laboratories: Using the game as a low-risk environment for policymakers to test the social legitimacy of new regulations, such as land-use taxes or rewilding mandates, before their real-world rollout.

    Overall, Pathbreak illustrates the promise of social dilemma games as laboratories for experiencing practically limitless scenarios of how one can address social-environmental challenges and reimagining what is possible. It operationalizes Ostrom’s insight that people can achieve outcomes “better than rational” when reciprocity, trust, and dialog are fostered5. This is particularly fitting as Ostrom and her colleagues’ advances in institutional scholarship served as the foundation for the landmark IPBES global assessment1 that stressed the importance of addressing indirect drivers of biodiversity loss. More broadly, embedding such tools within broader participatory processes—involving citizens, scientists, and policymakers—can be illuminating for assessing their broader contextualized impacts94. It may likewise help bridge the persistent gap between biodiversity/STEM science and the societal transformations urgently needed to halt and reverse biodiversity loss while democratizing the science itself and making it more equitable95,96.MethodsBackground to the projectThe Project PLANET4B, funded by the European Union’s Horizon Europe research and innovation programme, grant agreement No 101082212 (2022-2025), as detailed on www.planet4b.eu, was an initiative focused on understanding and facilitating social change to foster improved biodiversity prioritization in decision-making at individual, community, and institutional levels (see also Soliev et al.9). PLANET4B is an acronym that stands for understanding Plural values, intersectionality, Leverage points, Attitudes, Norms, behavior and social lEarning in Transformation for Biodiversity decision making. By integrating interdisciplinary research methods, the project aimed to identify key institutional and behavioral drivers and develop innovative strategies that promote responsible actions across various sectors and better take into account potential social and environmental externalities. As a transdisciplinary project, PLANET4B also collaborated with diverse stakeholders, including policymakers, businesses, and communities, to analyze and create practical solutions that support processes towards a more biodiversity-rich and socially equitable future. Experiential learning has been one of many closely analyzed, explored, adopted, developed and piloted methods in PLANET4B.Game calibration: theory, the goal, and the contextIn Pathbreak, players face the stylized version of the direct social dilemma. We designed it in a way that short-term profit maximization leads to high production and consumption decisions. In the long term, if producers follow this path, it leads to negative feedback from the environment—a higher susceptibility to extreme weather events. They provide penalties for the community when biodiversity is reduced or low. To emphasize the dilemma, we decided to set a sustainable level of production at the lowest production decision.From the perspective of a single round, the value for the production is matched with the consumption. In theory, an imbalance between the two can create a shortage or surplus that will benefit either producers or consumers through the increase or decrease of the value of one unit of food in the market. Still, it affects the group payoff for the specific type of players. Within the group, the payoffs are distributed proportionally, according to the decision of each player. It means higher individual production/consumption, when all other elements are kept stable, leads to a higher individual outcome. Therefore, each player has an incentive to increase their production or consumption. If all decisions are maximized no player has an incentive to change their action, thus creating Nash equilibrium. This way, rational players should not create artificial shortages or surpluses. Finally, from the perspective of a single round, Nash equilibrium is not affected by the occurrence of an event. Even if there is a probability of a penalty, it does not depend on the actions in the current round. Therefore, maximizing production or consumption increases the expected payoff even in the face of potential events that can reduce payoff (extreme weather events).To calibrate long-term solutions, we used a set of four strategies as a benchmark (Fig. 8). The first two are derived from the general idea of social dilemma. Myopic maximization means that all players follow the single-round Nash equilibrium, and the policies are set for maximum production/consumption. On the other hand, a sustainable strategy means that all players choose in unison sustainable production and consumption along with the matching policies. Alongside these two strategies, we added two more, to have a better understanding of the potential results. The first one shows a short recovery loop, where players maximize their production/consumption when biodiversity is high and switch to sustainable production/consumption when it is reduced, allowing it to recover. The second one follows a long recovery loop, where players maximize their actions until the biodiversity is low. Then, they switch to sustainable production/consumption until it recovers to a high level. At this point, they start maximizing their actions again. Therefore, players create a long cycle of maximizing production and recovery from the low biodiversity.Fig. 8Full size imageCalibration of accumulated payoff and biodiversity for benchmark strategies.After establishing these benchmark strategies, we were able to adjust the parameters linked to events (probability and value of the penalty) so that the sustainable strategy would yield the highest long-term payoff for the participants, while the myopic maximization the lowest one, with recovery loop strategies in the middle. Furthermore, we had to take into account that the policies are initially set up for maximizing production/consumption and can be only changed after the players’ decisions. Rounds of community council meetings are organized to discuss anything players would like to discuss. This is then followed by elections or referendum that determines which players serve as a governor and can make decisions for the subsequent round of policies: land (high or low land availability for agriculture, with respectively no or some areas allocated for rewilding that is expected to increase biodiversity in line with Pereira et al.97) and tax (high or low tax for households to cover the externalities, see also Carattini et al.98 for carbon taxes). The first round starts with policies maximizing production and consumption to increase ecological validity of contextualization. Therefore, in the first round, from the perspective of individual points it is optimal to maximize the payoff, only switching to a predefined strategy from the second round. Finally, as the game is planned for seven rounds it is also optimal for individual points to make the maximum decision in the final round. Payoff values represent expected payoff, taking into account penalties from the events weighted by the probability of event.Based on the social dilemma, when all players follow the sustainable strategy, they should achieve the highest individual payoff by the end of the game, while short-term maximization strategy leads to the lowest payoff by the end of the game. The strategy that maximizes production and consumption when biodiversity is high provides an outcome similar to but slightly lower than the sustainable one. Still, it has to be emphasized that these results present the expected payoff, taking into account the probability of events and penalties linked to them. In a single game, the results may differ due to the probabilistic nature of these events.As mentioned before, the first round has fixed decisions that generate a balanced market and a group payoff of 600. If players reduce production and consumption to a sustainable level, they can accumulate a payoff of up to 1968 in the long term. One-round results decrease but there is no significant risk of disasters. Myopic optimization results in an expected negative group payoff of −444, taking into account the probabilities of the penalties. Individual results may vary due to the imbalance in the market or the randomness of the disaster events.Pilots, analysis of debriefings, and reflectionsThe game was continuously pre-tested by the authors and their teams. The piloting followed a purposive sampling99, where we were interested to test the earlier versions of the game with specific audiences—we wanted to ensure it resonates with experts in biodiversity and participants without expertise in biodiversity. Here we report on the results of our in-depth content analysis of debriefing transcripts from the five pilots with a total of 50 participants conducted during the development of the game (more than one participant can play as one player). The pilots included (1) two groups of experts that represented the consortium partners with expertise on biodiversity-related topics in the context of different sectors (agriculture, education, trade, finance, fashion), intersectionality dimensions (age, gender, ethnicity, disability, religion), and geographic focus (Austria, Belgium, Brazil, Czech Republic, Germany, Hungary, Italy, Netherlands, Norway, Switzerland, UK) that met at the workshop conducted in Berlin on 18 January 2024 (n = 26; 16 female, 10 male; information on age was not collected, but it is a mixed and professionally-active age adult group of approximately 25–60); (2) sociology students studying structural change in a society as part of the course “SO2 Economy, State and Society” at the Martin Luther University Halle-Wittenberg that met in Halle on 2 June 2024 (n = 9; 5 female, 3 male, 1 non-binary, information on age was not collected, but it is a group of bachelor students mid-way through their studies in sociology and political sciences); and (3) a group of young people, including both prior participants in biodiversity-related activities of the learning community of the PLANET4B case study on urban youth and nature and newcomers primarily drawn by interest in the game itself that met in Erfurt on 9 September 2024 (n = 8; 1 female, 7 male; one 19-year old, four university students with migration background in their mid-20s, one international volunteer in their mid-20s, two young adults in their late-20s, one with family migration history and one unemployed after their bachelor studies; three of the participants played the earlier version of the game); and (4) an external group of audience that had no direct connection to the project, who attended the game piloting in Hamburg on 17 September 2024 (n = 7; 5 female, 2 male; information on age was not collected, one senior lecturer at a university in Hamburg, and a group of practitioners in the field of urban planning and architecture, two of whom were also students) (see the example photos in Fig. 9).Fig. 9Full size imageMoments captured from the game pilots with the project-related learning community members (left) and participants beyond the project (right).An iterative deductive and inductive analytical process was employed to explore patterns emerging from the qualitative data following Standards for Reporting Qualitative Research100. With the informed consents of the participants, the debriefings were recorded and transcribed using AI assisted transcription of sound and three rounds of subsequent manual correction. The transcripts were then analyzed using in-depth content analysis99,101, where initial codes were developed based on the theories discussed earlier beginning with open coding informed by close engagement with the transcripts, and iteratively adjusted based on the empirical evidence. Analysis was conducted using MAXQDA qualitative data analysis software. One researcher undertook the primary coding, and regular weekly meetings with two supervisors were conducted to enable collective reflection on emerging interpretations until we reached the theoretical saturation. An interpretive agreement was reached by first anchoring our initial codes in the discussed theories of social dilemmas. When the primary coder encountered participant statements that did not fit this deductive framework or shared critical feedback about the game’s premise and mechanics, these instances were treated as “critical cases”. These cases were actively sought out in subsequent iterations and brought to the weekly supervisory meetings, where the team collectively agreed on new, inductive codes (such as “affective-relational focus” or “frustration about the challenge”) to capture the added value from the empirical data. This iterative regrouping ensured that initially novel participant viewpoints drove the evolution of the coding framework rather than being forced into pre-existing categories. A framework with a hierarchy of codes (Table 1; see also the Supplementary Information for illustrative examples of quotations) evolved throughout the analysis, with new themes and sub-themes prompting iterative revisions and regrouping of codes to ensure alignment with the developing understanding of the data.Table 1 The coding framework for debriefing analysisFull size tableEthical statementEthical approval for this project was obtained from the German Association for Experimental Economic Research e.V. (Institutional Review Board Certificate No. Yuh2H9y1) and from the Ethics Commission of the Martin Luther University Halle-Wittenberg (Reference Approval No. 2308SIU). The authors have obtained written informed consent from the participants before joining the study as well as to publish images.

    Data availability

    The pseudonymized data from the debriefing sessions that were recorded, transcribed, and analyzed for this research is available at https://doi.org/10.5281/zenodo.17476103. The fully coded datasets are not publicly available due to notes that contain information that can lead to identification of individuals but are available from the corresponding author on reasonable request.
    ReferencesIPBES. Global Assessment Report on Biodiversity and Ecosystem Services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES Secretariat, 2019).IPBES. The Thematic Assessment Report on the Transformative Change and Biodiversity of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES Secretariat, 2024).Cardinale, B. J. et al. Biodiversity loss and its impact on humanity. Nature 486, 59–67 (2012).Article 
    CAS 
    PubMed 

    Google Scholar 
    Keck, F. et al. The global human impact on biodiversity. Nature 641, 395–400 (2025).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Ostrom, E. A behavioral approach to the rational choice theory of collective action: presidential address, American Political Science Association, 1997. Am. Polit. Sci. Rev. 92, 1–22 (1998).Article 

    Google Scholar 
    Ostrom, E. Polycentric systems for coping with collective action and global environmental change. Glob. Environ. Change 20, 550–527 (2010).Article 

    Google Scholar 
    Dasgupta, P. The Economics of Biodiversity (Cambridge University Press, 2024).Fesenfeld, L. P., Candel, J. & Gaupp, F. Governance principles for accelerating food systems transformation in the European Union. Nat. Food 4, 826–829 (2023).Article 
    PubMed 

    Google Scholar 
    Soliev, I., Zolyomi, A. & Franklin, A. Reflexivity-Situatedness Matrix: an intervention centered analytical framework for enabling transformations through transdisciplinary partnerships. Earth Steward. 2, e70020 https://doi.org/10.1002/eas2.70020 (2025).Article 

    Google Scholar 
    Williamson, O. E. The new institutional economics: taking stock, looking ahead. J. Econ. Lit. 38, 595–613 (2000).Article 

    Google Scholar 
    Meyer, J. W. World society, institutional theories, and the actor. Annu. Rev. Sociol. 36, 1–20 (2010).Article 

    Google Scholar 
    North, D. C. Institutional change: a framework of analysis. In Social Rules 189–201 (Routledge, 2018).Meinzen-Dick, R. et al. Playing games to save water: collective action games for groundwater management in Andhra Pradesh, India. World Dev. 107, 40–53 https://doi.org/10.5751/ES-15836-300125 (2018).Article 

    Google Scholar 
    ElDidi, H. et al. Getting ahead of the game: experiential learning for groundwater governance in Ethiopia. Int. J. Commons 18, 66–81 https://doi.org/10.5334/ijc.1316 (2024).Article 

    Google Scholar 
    Bell, A. R. et al. Understanding pastoralist adaptations to drought via games and choice experiments: field testing among Borana communities. Ecol. Soc. 30, 25 (2025).Article 

    Google Scholar 
    Lukosch, H. K., Bekebrede, G., Kurapati, S. & Lukosch, S. G. A scientific foundation of simulation games for the analysis and design of complex systems. Simul. Gaming 49, 279–314 (2018).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Janssen, M. A., Falk, T., Meinzen-Dick, R. & Vollan, B. Using games for social learning to promote self-governance. Curr. Opin. Environ. Sustain. 62, 101289 https://doi.org/10.1016/j.cosust.2023.101289 (2023).Article 

    Google Scholar 
    García-Barrios, L., García-Barrios, R., Waterman, A. & Cruz-Morales, J. Social dilemmas and individual/group coordination strategies in a complex rural land-use game. Int. J. Commons 5, 364–387 https://doi.org/10.18352/ijc.289 (2011).Article 

    Google Scholar 
    Falk, T. et al. Games for experiential learning: triggering collective changes in commons management. Ecol. Soc. 28, 30 https://doi.org/10.5751/ES-13862-280130 (2023).Article 

    Google Scholar 
    Levitt, S. D. & List, J. A. What do laboratory experiments measuring social preferences reveal about the real world? J. Econ. Perspect. 21, 153–174 (2007).Article 

    Google Scholar 
    Castillo, D., Bousquet, F., Janssen, M. A., Worrapimphong, K. & Cardenas, J. C. Context matters to explain field experiments: results from Colombian and Thai fishing villages. Ecol. Econ. 70, 1609–1620 (2011).Article 

    Google Scholar 
    Alekseev, A., Charness, G. & Gneezy, U. Experimental methods: when and why contextual instructions are important. J. Econ. Behav. Organ. 134, 48–59 (2017).Article 

    Google Scholar 
    Naar, N. Gaming anthropology: The problem of external validity and the challenge of interpreting experimental games. Am. Anthropol. 122, 784–798 (2020).Article 

    Google Scholar 
    Bush, J., DeCaro, M. S. & DeCaro, D. A. Playing a social dilemma game as an exploratory learning activity before instruction improves conceptual understanding. J. Exp. Psychol. Appl. 29, 725–746 (2023).Article 
    PubMed 

    Google Scholar 
    Hordijk, D. J., Andreotti, F. & van Zanten, H. H. Food system games for sustainability transformation—a review. Glob. Food Secur. 45, 100864 https://doi.org/10.1016/j.gfs.2025.100864 (2025).Article 

    Google Scholar 
    Motlagh, M., Horcea-Milcu, A. I. & König, B. Discovering the potential of serious games for transformative sustainability research. Discov. Sustain. 6, 30 https://doi.org/10.1007/s43621-024-00756-8 (2025).Article 

    Google Scholar 
    Olson, M. The Logic of Collective Action: Public Goods and the Theory of Groups (Harvard University Press, 1974).Dawes, R. M. Formal models of dilemmas in social decision-making. In Hum. Judg. Decis. Process. 87–107 (Academic Press, 1975). https://doi.org/10.1016/B978-0-12-397250-7.50010-6.Bohm, P. Estimating demand for public goods: an experiment. Eur. Econ. Rev. 3, 111–130 (1972).Article 

    Google Scholar 
    Ostrom, E., Gardner, R. & Walker, J. Rules, Games, and Common-Pool Resources (University of Michigan Press, 1994).Fischbacher, U., Gächter, S. & Fehr, E. Are people conditionally cooperative? Evidence from a public goods experiment. Econ. Lett. 71, 397–404 (2001).Article 

    Google Scholar 
    Tenbrunsel, A. E. & Messick, D. M. Sanctioning systems, decision frames, and cooperation. Admin. Sci. Q. 44, 684–707 (1999).Article 

    Google Scholar 
    Ostrom, E. The value-added of laboratory experiments for the study of institutions and collective action. J. Econ. Behav. Organ. 61, 149–163 (2006).Article 

    Google Scholar 
    Denant-Boemont, L., Masclet, D. & Noussair, C. N. Punishment, counter-punishment and sanction enforcement in a social dilemma experiment. J. Econ. Theory 33, 145–167 (2007).Article 

    Google Scholar 
    Cardenas, J. C., Janssen, M. & Bousquet, F. Dynamics of rules and resources: three new field experiments on water, forests and fisheries. In Handbook on Experimental Economics and the Environment 319–345 (Edward Elgar Publishing, 2013).Ostrom, E., Walker, J. & Gardner, R. Covenants, with and without a sword: self-governance is possible. Am. Polit. Sci. Rev. 86, 404–417 (1992).Article 

    Google Scholar 
    Fehr, E. & Gächter, S. Cooperation and punishment in public goods experiments. Am. Econ. Rev. 90, 980–994 (2000).Article 

    Google Scholar 
    Brosig, J. & Weimann, J. The effect of communication media on cooperation. Ger. Econ. Rev. 4, 217–241 (2003).Article 

    Google Scholar 
    Bochet, O., Page, T. & Putterman, L. Communication and punishment in voluntary contribution experiments. J. Econ. Behav. Organ. 60, 11–26 (2006).Article 

    Google Scholar 
    Biel, A. & Thøgersen, J. Activation of social norms in social dilemmas: a review of the evidence and reflections on the implications for environmental behaviour. J. Econ. Psychol. 28, 93–112 (2007).Article 

    Google Scholar 
    Jones, N. A., Ross, H., Lynam, T., Perez, P. & Leitch, A. Mental models: an interdisciplinary synthesis of theory and methods. Ecol. Soc. 16, 46 https://doi.org/10.5751/ES-03802-160146 (2011).Article 

    Google Scholar 
    Sagiv, L., Sverdlik, N. & Schwarz, N. To compete or to cooperate? Values’ impact on perception and action in social dilemma games. Eur. J. Soc. Psychol. 41, 64–77 (2011).Article 

    Google Scholar 
    Vesely, S., Klöckner, C. A. & Brick, C. Pro-environmental behavior as a signal of cooperativeness: evidence from a social dilemma experiment. J. Environ. Psychol. 67, 101362 https://doi.org/10.1016/j.jenvp.2019.101362 (2020).Article 

    Google Scholar 
    Insko, C. A., Schopler, J., Hoyle, R. H., Dardis, G. J. & Graetz, K. A. Individual-group discontinuity as a function of fear and greed. J. Pers. Soc. Psychol. 58, 68–79 (1990).Article 

    Google Scholar 
    Bornstein, G. The free-rider problem in intergroup conflicts over step-level and continuous public goods. J. Pers. Soc. Psychol. 62, 597–606 (1992).Article 

    Google Scholar 
    Dawes, R. M. & Messick, D. M. Social dilemmas. Int. J. Psychol. 35, 111–116 (2000).Article 

    Google Scholar 
    Klein, G. A. Sources of Power: How People Make Decisions (MIT Press, 1998).Ostrom, E. Governing the Commons: The Evolution of Institutions for Collective Action (Cambridge University Press, 1990).Ostrom, E. Understanding Institutional Diversity (Princeton University Press, 2005).Loibl, K., Roll, I. & Rummel, N. Towards a theory of when and how problem solving followed by instruction supports learning. Educ. Psychol. Rev. 29, 693–715 (2017).Article 

    Google Scholar 
    Schwartz, D. L. & Bransford, J. D. A time for telling. Cogn. Instr. 16, 475–522 (1998).Article 

    Google Scholar 
    Parks, C. D., Joireman, J. & Van Lange, P. A. M. Cooperation, trust, and antagonism: how public goods are promoted. Psychol. Sci. Public Interest 14, 119–165 (2013).Article 
    PubMed 

    Google Scholar 
    DeCaro, D. A., Arnol, C. A. T., Boama, E. F. & Garmestani, A. S. Understanding and applying principles of social cognition and decision making in adaptive environmental governance. Ecol. Soc. 22, 33 https://doi.org/10.5751/ES-09154-220133 (2017).Article 

    Google Scholar 
    Folke, C. et al. Resilience and sustainable development: building adaptive capacity in a world of transformations. Ambio 31, 437–440 (2002).Article 
    PubMed 

    Google Scholar 
    Baird, J., Plummer, R., Haug, C. & Huitema, D. Learning effects of interactive decision-making processes for climate change adaptation. Glob. Environ. Change 27, 51–63 (2014).Article 

    Google Scholar 
    Blackmore, C. What kinds of knowledge, knowing and learning are required for addressing resource dilemmas? A theoretical overview. Environ. Sci. Policy 10, 512–525 (2007).Article 

    Google Scholar 
    Parreño, M. A. et al. How puzzles are shaping our understanding of biodiversity: a call for more research into biodiversity representation in educational games. GAIA 31, 139–145 (2022).Article 

    Google Scholar 
    Phalan, B., Onial, M., Balmford, A. & Green, R. E. Reconciling food production and biodiversity conservation: land sharing and land sparing compared. Science 333, 1289–1291 (2011).Article 
    CAS 
    PubMed 

    Google Scholar 
    Tscharntke, T. et al. Global food security, biodiversity conservation and the future of agricultural intensification. Biol. Conserv. 151, 53–59 (2012).Article 

    Google Scholar 
    Soliev, I. Social-technological dimensions of land use in agriculture: applying and advancing the IPBES framework for transformative change assessment in the examples of land sharing and land sparing. Encycl. Agric. Food Syst. 6, 143–160 (2026).
    Google Scholar 
    Isbell, F. et al. Biodiversity increases the resistance of ecosystem productivity to climate extremes. Nature 526, 574–577 (2015).Article 
    CAS 
    PubMed 

    Google Scholar 
    Loreau, M. et al. Biodiversity as insurance: from concept to measurement and application. Biol. Rev. 96, 2333–2354 (2021).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Mahecha, M. D. et al. Biodiversity loss and climate extremes—study the feedbacks. Nature 612, 30–32 (2022).Article 
    PubMed 

    Google Scholar 
    Balmford, A., Bateman, I. J., Eyres, A., Swinfield, T. & Ball, T. S. Sustainable high-yield farming is essential for bending the curve of biodiversity loss. Philos. Trans. R. Soc. B 380, 20230216 (2025).Article 

    Google Scholar 
    Stern, P. C. New environmental theories: toward a coherent theory of environmentally significant behavior. J. Soc. Issues 56, 407–424 (2000).Article 

    Google Scholar 
    Janssen, M. A., Gharavi, L. & Yichao, M. Keeping up shared infrastructure on a Port of Mars: an experimental study. Int. J. Commons 14, 404–417 https://doi.org/10.5334/ijc.1017 (2020).Article 

    Google Scholar 
    Crenshaw, K. Demarginalizing the intersection of race and sex: a Black feminist critique of antidiscrimination doctrine, feminist theory and antiracist politics. Univ. Chic. Leg. Forum 1989, 139–167 (1989).
    Google Scholar 
    Mangelsdorf, M., Pregernig, M. & Kuni, V. (Bio-) diversity, gender, and intersectionality. FZG Freib. Z. Geschlechterstudien 22, 3–4 (2016).
    Google Scholar 
    Gonzalez, C. Racial capitalism, climate justice, and climate displacement. OSLS 11, 108–147 (2020).Article 

    Google Scholar 
    Kolb, D. A. Experiential Learning: Experience as the Source of Learning and Development (Pearson Education, 2015).Thatcher, D. C. Promoting learning through games and simulations. Simul. Gaming 21, 262–273 (1990).Article 

    Google Scholar 
    Lederman, L. C. Debriefing: toward a systematic assessment of theory and practice. Simul. Gaming 23, 145–160 (1992).Article 

    Google Scholar 
    Crookall, D. Serious games, debriefing, and simulation/gaming as a discipline. Simul. Gaming 41, 898–920 (2010).Article 

    Google Scholar 
    Kriz, W. C. A systemic-constructivist approach to the facilitation and debriefing of simulations and games. Simul. Gaming 41, 663–680 (2010).Article 

    Google Scholar 
    Wouters, P., Van Nimwegen, C., Van Oostendorp, H. & Van Der Spek, E. D. A meta-analysis of the cognitive and motivational effects of serious games. J. Educ. Psychol. 105, 249–265 (2013).Article 

    Google Scholar 
    Plass, J. L., Homer, B. D. & Kinzer, C. K. Foundations of game-based learning. Educ. Psychol. 50, 258–283 (2015).Article 

    Google Scholar 
    Loewenstein, G. Experimental economics from the vantage-point of behavioural economics. Econ. J. 109, 25–34 (1999).Article 

    Google Scholar 
    Garcia, C. A. et al. Strategy games to improve environmental policymaking. Nat. Sustain. 5, 464–471 (2022).Article 

    Google Scholar 
    Dernat, S., Grillot, M., Andreotti, F. & Martel, G. A sustainable game changer? Systematic review of serious games used for agriculture and research agenda. Agric. Syst. 222, 104178 (2025).Article 

    Google Scholar 
    Scheffer, M., Carpenter, S., Foley, J. A., Folke, C. & Walker, B. Catastrophic shifts in ecosystems. Nature 413, 591–596 (2001).Article 
    CAS 
    PubMed 

    Google Scholar 
    Milinski, M., Sommerfeld, R. D., Krambeck, H. J., Reed, F. A. & Marotzke, J. The collective-risk social dilemma and the prevention of simulated dangerous climate change. Proc. Natl. Acad. Sci. USA 105, 2291–2294 (2008).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Rocha, J. C., Schill, C., Saavedra-Díaz, L. M., Moreno, R. D. P. & Maldonado, J. H. Cooperation in the face of thresholds, risk, and uncertainty: experimental evidence in fisher communities from Colombia. PLoS ONE 15, e0242363 https://doi.org/10.1371/journal.pone.0242363 (2020).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Au, W. Criticality and environmental uncertainty in step-level public goods dilemmas. Group Dyn. 8, 40–61 (2004).Article 

    Google Scholar 
    Suleiman, R., Rapoport, A. & Budescu, D. V. Fixed position and property rights in sequential resource dilemmas under uncertainty. Acta Psychol. 93, 229–245 (1996).Article 

    Google Scholar 
    Maas, A., Goemans, C., Manning, D., Kroll, S. & Brown, T. Dilemmas, coordination and defection: How uncertain tipping points induce common pool resource destruction. Games Econ. Behav. 104, 760–774 (2017).Article 

    Google Scholar 
    Dannenberg, A., Löschel, A., Paolacci, G., Reif, C. & Tavoni, A. On the provision of public goods with probabilistic and ambiguous thresholds. Environ. Resour. Econ. 61, 365–383 (2015).Article 

    Google Scholar 
    Schill, C., Lindahl, T. & Crépin, A. S. Collective action and the risk of ecosystem regime shifts: insights from a laboratory experiment. Ecol. Soc. 20 https://doi.org/10.5751/ES-07318-200148 (2015).Ahsanuzzaman, Palm-Forster, L. H. & Suter, J. F. Experimental evidence of common pool resource use in the presence of uncertainty. J. Econ. Behav. Organ. 194, 139–160 (2022).Article 

    Google Scholar 
    Schill, C. & Rocha, J. C. Sustaining local commons in the face of uncertain ecological thresholds: Evidence from a framed field experiment with Colombian small-scale fishers. Ecol. Econ. 207, 107605 (2023).Article 

    Google Scholar 
    Hembach-Stunden, K., Vorlaufer, T. & Engel, S. Threshold ambiguity and sustainable resource management: a lab experiment. Ecol. Econ. 226, 108353 (2024).Article 

    Google Scholar 
    Sweeney, L. B. & Sterman, J. D. Bathtub dynamics: initial results of a systems thinking inventory. Syst. Dyn. Rev. 16, 249–286 (2000).Article 

    Google Scholar 
    Sterman, J. D. Risk communication on climate: mental models and mass balance. Science 322, 532–533 (2008).Article 
    CAS 
    PubMed 

    Google Scholar 
    Moxnes, E. Not only the tragedy of the commons: misperceptions of feedback and policies for sustainable development. Syst. Dyn. Rev. 16, 325–348 (2000).Article 

    Google Scholar 
    Rodela, R. & Speelman, E. N. Serious games in natural resource management: steps toward assessment of their contextualized impacts. Curr. Opin. Environ. Sustain. 65, 101375 (2023).Article 

    Google Scholar 
    Genua-Olmedo, A. et al. Breaking barriers: ten essential steps to achieve gender equality in academia through scientific societies. npj Biodivers. 4, 37 https://doi.org/10.1038/s44185-025-00105-6 (2025).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Bonet, R. J. A call to re-collectivize science and make it more equitable. npj Biodivers. 5, 2 https://doi.org/10.1038/s44185-025-00115-4 (2026).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Pereira, H. M., Hines, J. & Fernández, N. Conservation: meeting biodiversity targets through rewilding. Curr. Biol. 34, R918–R921 (2024).Article 
    CAS 
    PubMed 

    Google Scholar 
    Carattini, S., Carvalho, M. & Fankhauser, S. Overcoming public resistance to carbon taxes. Wiley Interdiscip. Rev. Clim. Change 9, e531 https://doi.org/10.1002/wcc.531 (2018).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Bryman, A. Social Research Methods (Oxford University Press, 2016).O’Brien, B. C., Harris, I. B., Beckman, T. J., Reed, D. A. & Cook, D. A. Standards for reporting qualitative research: a synthesis of recommendations. Acad. Med. 89, 1245–1251 (2014).Article 
    PubMed 

    Google Scholar 
    Mayring, P. Qualitative Content Analysis: Theoretical Foundation, Basic Procedures and Software Solution (Beltz, 2014).Download referencesAcknowledgementsThe initial concept of this game has emerged as a result of the first author’s participation in the games workshop organized by Leticia Merino, Ruth Meinzen-Dick, and Kate Ashbrook at the 2014 International Association for the Study of the Commons Conference in Edmonton. This has sharpened during the research visit of the first author to the Center for Behavior, Institutions and the Environment at Arizona State University, directed by the last author. Many conversations with various colleagues during these times have indirectly contributed to the initial concept. We thank the PLANET4B consortium, particularly as part of Tasks 2.1 and 2.2, for their valuable inputs throughout the project to understand and debate interventions, as well as their willingness to experiment with them in practice. We thank Torsten Wähler, Alex Franklin, Agnes Zolyomi, Sven Grüner, and Lindy Binder for the continuous and helpful discussions related to various aspects of the game. We thank Paul Wongsam for support in developing professional illustrations. We also thank Vladislav Artiukhov, Mariana Matoso, Zafar Saydaliev, Kevin Verduci, the students of the SO2 course at the Martin Luther University Halle-Wittenberg, the participants of the PLANET4B learning community meetings at the Culture Goes Europe e.V., and others beyond the project team for participating in the pilots. The support from Annabel Eckmann and Julian Jäger in transcribing and coding the debriefing sessions is gratefully acknowledged. We are also grateful to numerous participants at the 2023 Nairobi Conference of the International Association for the Study of the Commons, 2023 Radboud Conference on Earth System Governance, 2024 Berlin Workshop for Facilitators on Creative Methods, 2025 Brussels Biodiversity Week, and 2026 Brussels Transformative Change Cluster Projects’ Final Conference for the valuable feedback on the earlier works related to this manuscript. This work has partly evolved as part of the PLANET4B project funded by the European Union’s Horizon Europe research and innovation programme under grant agreement No 101082212. It was additionally funded by the UK Research and Innovation (UKRI) under the UK government’s Horizon Europe funding guarantee and the Swiss State Secretariat for Education, Research and Innovation (SERI).FundingOpen Access funding enabled and organized by Projekt DEAL.Author informationAuthors and AffiliationsDepartment of Environmental Sociology, Institute of Sociology, Faculty of Philosophy I: Social Sciences and Historical Cultural Studies, Martin Luther University Halle-Wittenberg, Halle (Saale), GermanyIlkhom SolievSocial Learning and Environmental Governance (SLEG) Lab, Centre for Interdisciplinary Regional Studies (ZIRS), Martin Luther University Halle-Wittenberg, Halle (Saale), GermanyIlkhom Soliev & Michał PająkGerman Centre for Integrative Biodiversity Research (iDiv), Leipzig, GermanyIlkhom SolievWroclaw University of Economics and Business, Wrocław, PolandMichał PająkCulture Goes Europe (CGE) Erfurt e.V., Erfurt, GermanyMaryna BykovaCenter for Behavior, Institutions and the Environment (CBIE), Arizona State University, Tempe, AZ, USAMarco A. JanssenAuthorsIlkhom SolievView author publicationsSearch author on:PubMed Google ScholarMichał PająkView author publicationsSearch author on:PubMed Google ScholarMaryna BykovaView author publicationsSearch author on:PubMed Google ScholarMarco A. JanssenView author publicationsSearch author on:PubMed Google ScholarContributionsI.S. conceived the study, conceptualized the model, led the design of the game, undertook formal analysis, secured funding, conducted investigations, developed methodology, provided resources, supervised the overall project, validated results, created visualizations, and wrote the original draft. I.S. also led reviewing, editing, and revising the manuscript. M.P. contributed to the design of the game, carried out formal analysis and methodological work, contributed to validation and visualizations, contributed to drafting of the original manuscript, and participated in reviewing, editing, and revision. M.B. performed analyses, supervised pilot work, validated findings, and created visualizations, as well as reviewed, edited, and revised the manuscript. M.A.J. contributed to the conceptualization and methodological design of the study, and participated in reviewing, editing, and revision of the manuscript. All authors reviewed the manuscript.Corresponding authorCorrespondence to
    Ilkhom Soliev.Ethics declarations

    Competing interests
    The authors declare no competing interests.

    Additional informationPublisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary information260907_Supplementary_Information_Pathbreak. (download PDF )Rights and permissions
    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
    Reprints and permissionsAbout this articleCite this articleSoliev, I., Pająk, M., Bykova, M. et al. Pathbreak: A Biodiversity-Food-Governance Game for experiential reimagination of social-environmental futures.
    npj biodivers 5, 35 (2026). https://doi.org/10.1038/s44185-026-00160-7Download citationReceived: 06 March 2026Accepted: 18 September 2026Published: 08 October 2026Version of record: 08 October 2026DOI: https://doi.org/10.1038/s44185-026-00160-7Share this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy shareable link to clipboard
    Provided by the Springer Nature SharedIt content-sharing initiative More

  • in

    Global patterns are key to understanding macroevolution

    AbstractRegional diversification patterns are increasingly being used to investigate spatial heterogeneity in biodiversity evolution and its drivers. However, estimating and interpreting regional diversification rates is complicated by the fact that regional diversity fluctuations reflect multiple biological processes and may be more sensitive to temporal sampling bias than global patterns. Here, we assess the biological fidelity and sampling sensitivity of regional diversification rates using fossil occurrences of marine invertebrates. We show that immigrations and extirpations account for up to 76% of apparent plate-level generic turnover. Moreover, regional diversity fluctuations are more strongly shaped by preservation, sampling, and tectono-sedimentary history than global patterns. In widely distributed groups, global analyses integrate data across regions and environments, dampen region-specific biases, and recover a higher proportion of shared biological signal. Global analyses therefore remain indispensable for inferring macroevolutionary dynamics over deep time.

    Explore related subjects
    Discover the latest articles and news in related subjects.

    Biodiversity

    Palaeontology

    IntroductionAn accurate reconstruction of deep-time diversification patterns is a prerequisite for identifying their underlying drivers. Over the past half century, analyses of the global fossil record have enabled the identification of important macroevolutionary patterns, such as the major diversification episodes, the three Great Evolutionary Faunas, and the ‘Big Five’ Phanerozoic mass extinctions1,2,3,4,5,6. These discoveries formed the foundation for subsequent research into the mechanisms underlying macroevolutionary dynamics7,8,9,10,11. Nevertheless, diversity and the processes driving diversification are, and have always been, spatially structured12,13,14,15,16,17. Global patterns, which represent aggregated signals across regions, may therefore obscure processes operating at regional scales. Moreover, a growing number of studies have shown that fossil sampling is spatially heterogeneous, further complicating the interpretation of global patterns13,14,16,17. These concerns have stimulated growing interest in the spatial structure of macroevolutionary patterns and placed increasing emphasis on analyses conducted at regional scales16,17,18,19,20. In some studies, the validity of global patterns has even been questioned16,17.However, estimating and interpreting diversification rates at the regional scale is inherently more complex and potentially more sensitive to bias than at global scale. First, while global patterns necessarily reflect global diversification dynamics (i.e., origination and extinction), regional first and last appearances reflect both diversification dynamics and regional range dynamics (i.e., immigration and extirpation). This distinction is important because environmental and ecological drivers of range shifts are not usually the same as those of origination and global extinction21, but most available methods do not explicitly separate these processes22,23. Second, both diversification and range dynamics are further filtered by preservation and sampling. Numerous studies have demonstrated strong correlations between fossil diversity and rock availability24,25,26,27,28,29,30,31,32,33,34, and the initiation or truncation of sedimentary successions may generate flawed originations and extinctions in the fossil record35,36,37,38,39,40. Such effects may be amplified in regional datasets, which typically encompass fewer depositional settings and are more strongly influenced by idiosyncratic tectonic, sedimentary, and sampling histories. Indeed, the explanatory power of the rock record on biodiversity is greater within some regions than at larger scales27,34. Although these issues are well recognised, the extent to which regional diversification rates are generally affected by range dynamics, preservation and sampling biases has not previously been quantified.Here, we assess the biological significance of regional diversification rates derived from fossil occurrence data. We focus on three major groups of marine invertebrates (Bivalvia, Brachiopoda, Cephalopoda), which have relatively complete global fossil records41,42. In this study, regional diversification rates refer to rates calculated from first and last appearances of genera at regional scale, as in some previous studies16,17,43,44. Such rates may conflate true diversification dynamics with range dynamics. We used tectonic plates to delineate geographic regions (see Material and Methods). In these regions, we quantified the proportion of first and last appearances that may reflect range expansion or shrinkage as opposed to in situ origination or global extinction. We then evaluated whether apparent regional diversity fluctuations coincided with changes in the availability of fossil-bearing localities. Using marine bivalves as an example, we further illustrated how apparent regional diversification rates can be shaped by these processes, and how the global dataset can mitigate some of these region-specific biases. Our aim is not to diminish the value of regional studies, but to clarify the biases that should be considered when estimating and interpreting regional diversification patterns.ResultsA large proportion of regional generic turnover reflects range dynamicsTo investigate how many apparent regional generic turnovers represent ‘true’ originations (i.e., in situ originations) and extinctions (i.e., global extinctions), or range shifts (i.e., immigrations and extirpations), we compared the first and last appearances of taxa at the regional and global scales. The first and last appearances of genera on a given plate are defined as ‘regional events’ regardless of plate dimensions. Similar events have been treated as ‘regional originations and extinctions’ in some previous studies and used to estimate regional diversification rates16,17,43,44. Comparison between the regional and global datasets shows that many regional appearances occurred at least one stage after the global first appearance of the same genus, whereas many regional disappearances occurred at least one stage before its global last appearance (Fig. 1; Table S1). These regional events are here referred to as immigrations and extirpations, although they may also reflect artefacts of preservation and sampling.Fig. 1: Many regional appearances and disappearances reflect range shifts rather than origination and extinction.Full size imageA–C numbers of regional first appearances and immigrations for bivalves, brachiopods, and cephalopods, respectively. D–F numbers of regional last appearances and extirpations for the three groups. Each dot represents data from a specific tectonic plate in a specific stage. For example, 44 bivalve genera last appeared in South China during the Changhsingian, of which 11 were recorded on other plates in younger stages; the coordinate of this data point is therefore (44, 11) in D. Least-squares fits were applied to illustrate the proportion of such events (blue lines; intercept constrained to zero). However, their proportions are underestimated (see text). Data from post-mass-extinction intervals in A–C and from mass extinction intervals in D–F are highlighted with coloured circles; grey circles denote all other intervals. Grey dashed lines have a slope of one. Rhu, Rhuddanian; Fam, Famennian; Ind, Induan; Het, Hettangian; Dan, Danian; Hir, Hirnantian; Fra, Frasnian; Cha, Changhsingian; Rha, Rhaetian; Maa, Maastrichtian.A substantial proportion of regional first and last appearances represent immigrations and extirpations, respectively, ranging from 23% to 76% across clades, indicating that range dynamics are important components of apparent regional turnover (Fig. 1). Because range shifts were identified by comparison with the global dataset, their proportions are partly influenced by the degree of faunal endemism: plates containing higher proportions of endemic taxa tend to exhibit lower proportions of immigrations and extirpations (Fig. S5). These proportions also showed weak and often negative correlations with plate size (Spearman’s ρ = −0.21 to 0.01 across clades; Fig. S6), suggesting that larger regions are more likely to capture the global origination and extinction times of taxa. In addition, the proportion of extirpations is lower during intervals of high global extinction intensity (Spearman’s ρ = −0.48 to −0.25), indicating that regional disappearances are more likely to coincide with global extinctions during major extinction episodes (Fig. 1, S4; Table S1). Similarly, the proportion of immigrations tends to be lower following mass extinctions (Spearman’s ρ = −0.37 to −0.35; Fig. 1, S4; Table S1).The estimated contribution of immigrations and extirpations to regional diversity turnover varies markedly among clades, whether measured as weighted ratios across all plate-stage combinations or as slopes of the fitted relationships in Fig. 1. Bivalves show the highest values, with both metrics ranging from 0.72 to 0.76. Cephalopods show the lowest values, ranging from 0.23 to 0.28, whereas brachiopods are intermediate, ranging from 0.47 to 0.55 (Table S1). To explore the factors underlying these differences, we calculated genus durations, sampling rates, and geographic ranges for the three focal groups (see Methods for details of these metrics). There are pronounced differences in duration and sampling rate among groups (Fig. S11). Bivalves have the longest genus durations (median: 16.67 Myr) and the lowest sampling rates (0.78 occurrences Myr-1), whereas cephalopods have the shortest durations (5.82 Myr) and the highest sampling rates (1.14 occurrences Myr-1); brachiopods show intermediate values (11.52 Myr and 0.85 occurrences Myr-1, respectively). Geographic ranges are similar among the three groups, although bivalves tend to have slightly larger ranges (3,623, 3,008, and 2,981 km for bivalves, brachiopods, and cephalopods, respectively). The long durations, low sampling rates and wide geographical ranges of bivalves imply that a bivalve genus is more likely to disappear from one plate while persisting elsewhere, or to have originated earlier outside the focal plate, thereby elevating the proportions of immigration and extirpations. The high evolutionary rates45, readily identifiable morphologies and high dispersal abilities of cephalopods have contributed to their widespread use as index fossils. Consequently, the regional first and last appearances of cephalopods matche best their global stratigraphical ranges.The proportions of immigrations and extirpations reported here should be regarded as conservative estimates. New taxa originate in geographically restricted areas and may subsequently expand into other regions. As the temporal resolution of PBDB data is generally limited to stage levels, the exact place and time of origin cannot reliably be traced and the first appearance of a genus may be recorded on multiple plates within the same stage. Such cases were not counted as immigrations, even though some of them must represent later regional arrivals within that stage. This underestimation may be particularly pronounced for cephalopods, which have short durations but broad geographic distributions, suggesting rapid dispersal after origination. The same issue also contributes to the lower proportion of immigrations in the aftermath of mass extinctions, when newly originated taxa rapidly migrated across regions46,47. It is possible that regional extinctions occurred synchronously, but some extirpations cannot be detected because of limited temporal resolution.The high proportions of immigration and extirpation contributing to regional generic turnover are noteworthy. Origination and extinction rates are not biologically equivalent to immigration and extirpation rates, and the environmental or ecological drivers of these processes may differ (e.g., refs. 21,48,49,50; see below). Therefore, driver-seeking analyses based on apparent regional turnover (i.e., raw regional appearances and disappearances) require caution. Failure to distinguish diversification from range dynamics may lead to misleading inferences about the processes shaping regional turnover.Sampling and sedimentary history influence apparent turnover signalsWe assessed the influence of sampling on regional diversity patterns by quantifying how many regional events could be associated with changes in the availability of fossil-bearing localities. Regional events were classified as ‘sampling-associated’ when generic appearances or disappearances coincided with the emergence or loss of sampled fossil localities in equal-area grid cells (Fig. S3). Importantly, this classification does not imply that these events were caused solely by preservation or sampling bias. Changes in the availability of fossil-bearing localities may reflect tectono-sedimentary evolution and facies shifts that simultaneously alter both the distribution of suitable habitats and the probability of fossil preservation. For example, when a depositional environment becomes unsuitable for a taxon, the taxon may be locally extirpated or track its preferred habitat into another area that is not sampled, while the original fossil-bearing facies also disappears. Ecological habitat loss and loss of the sampling window can therefore generate the same observed pattern and may represent a common-cause effect39. Because these mechanisms cannot be distinguished reliably with our data, we grouped such events operationally as ‘sampling-associated’, while recognising that some may reflect genuine ecologically driven range shifts or extirpations.Our results show that a substantial fraction of regional events is indeed sampling-associated (Fig. 2). Notably, the proportions of sampling-associated events are relatively high across the three groups (range of weighted ratios and slopes of fitted lines: 0.60–0.80; Table S2), indicating that even the cephalopod fossil record, despite its relatively high sampling intensity, is strongly influenced by local environmental shifts and sampling bias. The proportions of these events also show a weak negative correlation with plate size (Spearman’s ρ = −0.30 to −0.03; Fig. S7), suggesting that larger regions tend to buffer taxon losses or gains associated with shifts in fossil localities. The global dataset is also affected by sampling bias, as shown by earlier studies24,25,26,38,39, but for clades analysed here, the proportion of sampling-associated originations and extinctions is slightly lower (0.51–0.74; Table S3), and extreme proportions approaching one are less frequent (Fig. S8). Sensitivity analyses using larger grid cells result in lower numbers and proportions of sampling-related events, but the overall patterns remain unchanged (Fig. S9, S10; Table S4, S5).Fig. 2: Regional appearances and disappearances are frequently associated with sampling bias.Full size imageA–C total regional first appearances versus first appearances in newly sampled cells. D–F total regional last appearances versus last appearances in subsequently unsampled cells. Each dot represents a certain plate of a certain stage. For example, in D, 44 bivalve genera last appeared in South China in the Changhsingian stage, and the cells of five genus occurrences were not sampled in the next stage (Induan); the coordinate of this data point is therefore (44, 5). Least-squares fits were applied to illustrate the proportion of such events (blue lines; intercept constrained to zero). Data from post-mass-extinction intervals in A–C and from mass extinction intervals in D–F are highlighted with coloured circles; grey circles denote all other intervals. Grey dashed lines have a slope of one. The edge length of hexagonal equal area grid cells is ~123 km. Abbreviations as in Fig. 1.Sampling rates and sampling completeness in the global occurrence datasets are also greater than in regional datasets. Although global datasets contain more occurrences overall, their sampling rates and completeness are not by design greater than those of regional subsets. When a clade exhibits some degree of biogeographical differentiation, a well-sampled region can have high sampling rates and completeness, whereas global estimates may be reduced by poorly sampled regions. Sampling metrics varied greatly through time (Fig. 3). Median regional sampling rates of our three focal groups (0.26–0.27) are much lower than those of the global dataset (0.65–0.71) (Table S6). Global datasets also exhibit higher median sampling completeness (0.81–0.86) than that on individual plates (0.66–0.71) (Table S6). Although some plates have higher sampling rates or completeness than the global datasets in certain stages, overall, the global datasets are expectedly more complete and more continuous than regional datasets (one-sided Mann-Whitney test: p < 0.003, except for cephalopod sampling completeness, p = 0.053; Fig. 3). The advantage of global integration suggests that the analysed clades have relatively broad geographic distributions and many taxa are shared among regions.Fig. 3: Global sampling is generally more complete than regional sampling.Full size imageA–C sampling rate (number of occurrences per million years) for bivalves, brachiopods, and cephalopods. Sampling rates are log-transformed; regional sampling rates of zero (i.e., intervals with no fossil occurrences) are omitted. D–F three-timer sampling completeness. In all panels, open circles represent regional values and solid dots represent values of the global dataset. Although some regions show higher sampling rates or completeness than the global dataset in particular stages, the global dataset generally exhibits better sampling than regional datasets overall. This difference is supported by one-sided Mann-Whitney tests: p < 0.003 in A–E and p = 0.053 in F. Sample sizes for the global dataset are n = 81, 81, 81, 79, 79, and 79 in A–F, respectively; sample sizes for regional datasets are n = 1694, 1309, 1420, 743, 489, and 353 in A–F, respectively. See Table S6 for numerical comparisons.Bivalves showcase the discrepancy between regional and global diversification ratesDiscrepancies between plate-level and global diversification rates in marine bivalves illustrate how inferred regional diversification dynamics can be influenced by range dynamics, regional sedimentary change, and local sampling bias (Fig. 4). The global extinction rate of bivalves is similar to rates reported for other marine groups3,4,51,52, with pronounced peaks during the mass extinctions (Fig. 4A). However, regional disappearance rates vary wildly and often exhibit high values during intervals when the global extinction rate was low (Fig. 4A). After separating extirpations from final extinctions, regional extinction rates show less intense fluctuations during background times, further confirming that many regional disappearances represent extirpations. Nevertheless, several erratic peaks remained evident (Fig. 4B), indicating the spatial heterogeneity of observed extinction rates.Fig. 4: Regional extinction patterns can differ substantially from the global composite because of range dynamics and sampling variation.Full size imageA global extinction rates and regional disappearance rates in marine bivalves from the Silurian to the Miocene. Blue lines represent disappearance rates for individual plates, based on regional last appearances without distinguishing final extinction from extirpation; black line denotes the global extinction rate. B comparison of global and regional extinction rates. Red lines represent regional extinction rates; extirpations were not included in the calculation. Mass-extinction stages are highlighted by yellow bars in A, B. C regional disappearance rates (blue solid lines) and global extinction rates (black dashed lines) from the Capitanian to Hettangian for the Siberian craton, South China, and Adriatic plate, respectively. D regional extinction rates (red solid lines) and extirpation rates (purple dashed lines) for the three plates. Grey bars in C and D indicate the number of fossil-bearing cells on each plate. Abbreviations in A and B: O, Ordovician; S, Silurian; D, Devonian; C, Carboniferous; P, Permian; Tr, Triassic; J, Jurassic; K, Cretaceous; Pg, Paleogene; Ng, Neogene. Abbreviations in C and D: W, Wuchiapingian; C, Changhsingian; I, Induan; O, Olenekian; A, Anisian; L, Ladinian; Car, Carnian; Nor, Norian; R, Rhaetian; Reg., regional; Glo., global; disapp., disappearance; extin, extinction; extir, extirpation.Three relatively well-sampled regions exhibit very different patterns of apparent disappearance rates during the Triassic (Fig. 4C). On the Siberian craton, regional extinction rates show clear peaks associated with both the end-Permian and end-Triassic mass extinctions, but also exhibit elevated values during the Wuchiapingian, Olenekian, and Ladinian stages (Fig. 4C). In South China, regional disappearance rates record a strong signal of the end-Permian but little evidence for the end-Triassic mass extinction, with elevated rates in the Anisian and Norian stages (Fig. 4C). The rates in the Adriatic (Apulian) plate capture both the end-Permian and end-Triassic mass extinctions, but show a prominent Carnian peak that is not prominent in the global pattern (Fig. 4C).Decomposing regional disappearances into extirpations and extinctions, together with changes in the number of sampled fossil localities, helps explain these discrepancies (Fig. 4D). Both the end-Permian and end-Triassic mass extinctions were driven by severe environmental perturbations and had global-scale impacts53,54,55. However, the decline in sampled fossil localities in the Changhsingian of Siberia and the Late Triassic marine regression in South China shifted a substantial proportion of regional disappearances to stages preceding the corresponding global mass extinctions (Fig. 4D). Consequently, many extirpations predated the mass-extinction intervals. The Carnian Pluvial Episode also represented a global environmental perturbation that is thought to have profoundly affected animal evolution56. Carnian extinction rates exhibit pronounced spatial heterogeneity17 (Fig. 4C, D), and this pattern partly reflects contrasting sampling regimes. The high apparent diversity and extinction rate on the Adriatic Plate reflect oversampling of the Cassian Formation, a ‘liberation Lagerstätte’57, whereas the lower values in Siberia and South China likely result from the limited number of sampled fossil localities and the concentration of available occurrences in low-diversity shales. Other regional extinction and extirpation spikes, such as those in the Olenekian and Ladinian of Siberia and the Anisian of South China, also closely correspond to reductions in sampling (Fig. 4C, D).These examples show how regional diversification patterns can be distorted by a combination of range dynamics, sampling bias, and tectono-sedimentary evolution, as also illustrated by a recent study on brachiopods58. By integrating occurrences across plates and environments, global analyses dampen region-specific biases that confound local records, thereby recovering a clearer and more temporally consistent macroevolutionary signal.Implications for estimating diversification dynamics across spatial scalesBoth global and regional perspectives have distinct advantages and limitations. Diversification dynamics observed at different spatial and temporal scales may reflect different evolutionary processes13. Local and regional datasets can have relatively high temporal resolution, allowing their dynamics to capture fine-scale controls on diversity, as well as evidence for biotic interactions and local environmental conditions15,59,60. By contrast, global patterns are more likely to reflect large-scale mechanisms regulating diversity, such as climatic change and continental breakup10,61,62,63,64. Moreover, regional diversity fluctuations contain signals of range dynamics and can therefore provide insights into the formation of spatial diversity gradients and biotic invasions12,21,49,50. Global patterns, in contrast, are less affected by immigration and extirpation because these processes occur only within the global range of a taxon. The most appropriate spatial scale depends on the scientific question, and integrating global and regional perspectives can provide a more complete understanding of evolutionary drivers. However, given the current incompleteness of fossil data, regional patterns remain difficult to interpret, constraining the macroevolutionary information that can be inferred from them.Currently, estimating regional diversification rates remains challenging for both biological and methodological reasons. First, fluctuations in regional diversity reflect the combined effects of diversification and range dynamics. These processes can have different temporal-spatial patterns48,49,50 (Fig. 4), and should be distinguished when assessing their drivers. The most widely used analytical methods, including phylogenetic tree-based and fossil occurrence-based approaches, were primarily developed for processes operating at global scales23 (but see refs. 65,66). Although diversification and range dynamics can sometimes be separated analytically48,49,50, the low sample sizes of fossil occurrences in many regional datasets reduce the statistical power and reliability of regional diversification rate estimates22.Second, regional tectono-sedimentary history and fossil sampling complicate the detection of diversification drivers. Accounting for these processes is not straightforward because the stratigraphic distribution of fossils is non-random and strongly controlled by local environments and their sea-level-induced shifts35,40, and no single sampling proxy can fully capture these biases34,67,68. A range of analytical tools has been developed to mitigate temporal sampling bias in diversification rate estimation4,69,70,71,72,73,74, but these cannot yet fully accommodate the full suite of biases. These approaches can reduce the risk of inferring spurious rate shifts, or Type I errors, in which an apparent origination or extinction pulse is detected even though no such biological event occurred. For example, a decline in fossil sampling within an interval may produce an apparent extinction pulse, whereas sampling-corrected methods may weaken this signal or express the associated uncertainty through broad confidence intervals. However, these methods are less able to resolve Type II errors arising from false absences, in which genuine origination or extinction events are not detected because fossil records are sparse or discontinuous. For instance, the end-Triassic extinction is barely detected in South China because marine regression greatly reduced the availability of the marine fossil record beforehand (Fig. 4D). Consequently, the extent to which apparent regional diversification dynamics reflect intrinsic biological drivers remains uncertain.Global patterns are also subject to temporal and spatial sampling bias13,14,24,25,39,75,76,77,78. Nevertheless, by integrating information from multiple incomplete regional records, global datasets generally record more true originations and extinctions and can partly mitigate the effects of low sample size and local sedimentary shifts encountered in regional studies, and as a result, shared large-scale signals may be retained while some asynchronous regional noise is dampened. However, the performance of global diversification rate estimates depends on the biogeographic structure of the clade under study. When interregional exchange is frequent and biogeographical differentiation is relatively weak (e.g., marine animals with high dispersal ability like foraminifers50,79, particularly in globally warm climates64), larger spatial coverage enhances fossil record completeness, improves estimates of true origination and extinction times, and increases the reliability of diversification-rate estimates. Conversely, in systems characterised by strong provincialism and regionally distinct diversification histories (e.g., terrestrial animals75,76,80), aggregating regional datasets can provide little improvement in sampling completeness. In such cases, global analyses may capture only the largest-scale events shared across regions (e.g., mass extinctions). In reality, no clade is completely uniform in biogeography, and the spatial coverage of fossil data may vary through time13,17. Different spatial-standardisation methods have been developed to reduce apparent rate variation caused by changes in spatial coverage, allowing fairer comparisons of diversity patterns among regions and time bins14,17,19,75,76. However, these methods cannot account for sedimentary and sampling biases that shape regional diversity fluctuations. Moreover, by substantially reducing sample sizes, spatial standardisation may increase the risk of Type II errors.In summary, global and regional approaches are complementary, but their inferential power is asymmetrical under the constraints of the fossil record. At regional scales, reconstructing diversification dynamics and identifying their drivers remain limited by uneven sampling, sedimentary biases, and the overprinting effects of range dynamics. Analyses must navigate inherent trade-offs between spatial detail and temporal robustness, and between capturing local variability and recovering generalisable macroevolutionary patterns. For many clades, regional diversification signals remain too incomplete to support strong inferences. Integrative global analyses, however, can partially overcome these limitations by averaging across regions and environments, thereby isolating the shared macroevolutionary signal. Advancing beyond this asymmetry will require both increased fossil sampling and the development of analytical approaches that explicitly model and separate the interacting effects of diversification, dispersal, and sampling. Such advances are essential for moving from descriptive patterns toward a mechanistic understanding of biodiversity dynamics through deep time.MethodsFossil occurrence dataPost-Cambrian fossil occurrences of Bivalvia, Brachiopoda, and Cephalopoda were downloaded from the Paleobiology Database (PBDB) on 26/05/2026. The URLs used to download the datasets are provided in the Supplementary Material. Terrestrial records indicated by the ‘environment’ field in PBDB were discarded. Previous studies have shown that macroecological and macroevolutionary patterns may vary among taxonomic levels81,82. Given the greater completeness of genus-level data relative to species-level data83,84, and their widespread use in palaeontological studies, all analyses were carried out at the genus level. Subgeneric names were removed and only genus names were retained. Occurrences that cannot be assigned to a genus and those with ambiguous generic assignments (i.e., primary generic resolution noted by cf., aff., ?, informal, ex.gr., “”) were excluded. The PBDB is known to contain taxonomic errors that may affect estimates of diversification rates85,86,87. We filtered the PBDB data using the pacmacro_ranges function in fossilbrush R package87 (v1.0.6) to reduce the influence of taxonomic errors and wastebasket taxa on rate estimation. We applied the default settings and discarded occurrences falling outside the 90% density intervals of genus-level stratigraphic ranges. Although this procedure cannot eliminate all errors, random errors are unlikely to generate large-scale systematic biases85, particularly because our analyses compare regional and global patterns derived from the same underlying dataset.We assigned fossil occurrences to stratigraphic stages using the divDyn R package52 (v0.8.3). Palaeo-coordinates of occurrences were reconstructed with the rgplates R package88,89 (v0.6.1) using the PALEOMAP rotation model90. We used the ‘geoplates2’ field in the PBDB data as an indicator of regions, which records the plate IDs used in the GPlates software88 under the PALEOMAP model90 (Fig. S1). Plate IDs in GPlates define groups of geological features that share the same reconstructed motion through time88. Although they do not correspond exactly to physical tectonic plates, they may represent the motion of features located on those physical plates. Plates were employed as spatial units because: (1) they allow regional fossil records to be consistently tracked through time, which is important when estimating regional diversification dynamics16,17,22; (2) fauna on different plates tend to share a lower compositional similarity and marine invertebrate biogeographic regions often coincide broadly with major plates43,91; and (3) the sizes of fossil-bearing areas on plates are comparable to those of regional units analysed in some previous studies16,19 (Fig. S1). Occurrences that could not be assigned to a stratigraphic stage or plate were discarded. After data cleaning, a total of 157,231, 130,418, and 89,936 occurrences in 2472, 3324, and 3679 genera were retained for the bivalve, brachiopod, and cephalopod datasets, respectively. All subsequent analyses were performed on these occurrences without accounting for the collections to which they belonged.Regional eventsRegional events were defined as first and last appearances of genera on a given plate. These events were then classified in two downstream analyses.In the first analysis, we quantified the extent to which regional events reflect range dynamics (immigrations and extirpations) rather than diversification dynamics (origination and extinction rates). To do so, we compared the timing of each regional event with the corresponding global first and last appearances of the same genus. Following other fossil-based studies12,48,92, global first and last appearances were used as approximations of taxonomic origination and extinction, although these dates are affected by preservation and sampling. If a genus first appeared on a plate later than its global origination, that regional appearance was interpreted as an immigration into the focal plate. Similarly, regional disappearances earlier than final global extinctions were treated as regional extirpations. It should be noted that we referred to these events as immigrations and extirpations for simplicity, but they could also be the result of preservation and sampling gaps. The numbers of immigrations and extirpations were calculated for each plate-stage combination. We then quantified their proportions relative to all regional events using two approaches: (1) the weighted ratio, defined as the total number of immigrations and extirpations across all plates and stages divided by the total number of regional events, and (2) the slope of a linear regression between the number of immigrations and extirpations and the number of regional events on each plate and stage. The intercept of the linear fit was forced to zero because neither immigrations nor extirpations can occur in the absence of regional events. To reduce noise associated with small sample sizes, only plate-stage combinations containing more than 100 occurrences were considered in the analyses. We first classified regional events using occurrences from all plates, and then retained only those plate-stage combinations that met this threshold when calculating the proportions of different event types. Consequently, 26,480 of 157,231 bivalve occurrences, 20,178 of 130,418 brachiopod occurrences, and 20,267 of 89,936 cephalopod occurrences were excluded from the calculation of these proportions (Fig. S2). However, the classification of regional events was performed without this filter. The first and last appearances of a genus are sometimes recorded on multiple plates within the same stage. In such cases, we cannot determine which occurrences represent the earliest or latest records, and these regional events were not counted as range shifts. This treatment will artificially bias estimates of range shifts downward.In the second analysis, we assessed the impact of fossil sampling on the estimates of regional diversification rates by testing the correspondence between regional events and the changes in the availability of fossil localities. Changes in the spatial distribution of fossil localities were measured using grid-cell occupancy of fossil occurrences. Based on their palaeocoordinates, fossil occurrences were assigned to equal-area hexagonal cells with edge lengths of ~123 km, using the icosa R package93 (v0.12.0). A regional first appearance was classified as sampling-associated if the genus first appeared in cells that were not sampled in the preceding stage on the same plate. A regional last appearance was considered sampling-associated if all cells occupied by a genus in one stage were unsampled in the subsequent stage and the genus disappeared from the plate. An illustration of our approach is provided in Figure S3. Although plates drift through time, the palaeogeographical positions of fossil localities do not change substantially between adjacent geological stages (ca. 5 Myr in duration) and were thus assumed to be constant. The number of sampling-associated regional events was calculated for every plate and stage. Again, the proportions of sampling-associated events were quantified using the same weighted ratio and linear-fit approaches described above. The threshold of 100 occurrences was also applied, thereby excluding data points from poorly sampled plate-stage combinations with only a few cells (Fig. S2). We also investigated the relationship between the proportion of sampling-associated events and plate size. Plate size was quantified for each time interval as the maximum great-circle distance between occurrences on a given plate, using the divDyn package52. Because the spatial extent of sedimentary areas varies through time, the estimated maximum great-circle distance of each plate also varies through time.To compare the prevalence of sampling-associated events at regional and global scales, we also conducted the same analysis on the global dataset using the same method. Because the identification of sampling-associated events depends on grid cell size, we also performed a sensitivity analysis using larger cells with edge lengths of ~247 km.Region- and genus-specific sampling patternsTo assess the quality of the fossil record, we calculated sampling rates and sampling completeness for individual plates and the global dataset. Sampling rate was measured as the number of fossil occurrences per million years. For each fossil dataset and stage, we calculated the sampling rate of every genus by dividing its number of occurrences by the duration of the stage. The median across genera was then used to represent the sampling rate of that dataset in that stage. As a proxy for sampling completeness, we used the three-timer completeness metric4, calculated using the divDyn package52. The metric is defined as (frac{{t}_{3}}{{t}_{3}+{t}_{{part}}}), where t3 is the number of taxa sampled in three consecutive bins and tpart is the number of taxa sampled in the bins immediately before and after, but not within, the focal bin. A low value indicates that many taxa expected to occur in the focal bin were not sampled there, suggesting poor sampling completeness. The two metrics were calculated for each plate-stage combination and for the global dataset in each stage, and the results were plotted in Fig. 3. We further calculated the median and its 95% bootstrap interval, based on 100 bootstrap replicates, for each metric to provide overall estimates for the regional and global datasets (Table S6).We also calculated genus-specific sampling rate, duration, and geographical range for genera in each animal group. Because the exact ages of fossil occurrences are uncertain and long stages may lead to overestimation of taxon durations, we randomly sampled occurrence ages within stage boundaries and calculated the interval between first and last appearances. This procedure was repeated 100 times for each genus, and the median value was retained as the duration. Taxonomic sampling rate was then calculated as the total number of occurrences of the genus across its duration divided by that duration. Geographic range size was quantified as the maximum great-circle distance of all occurrences of the genus in all stages, and calculated using the divDyn package52. After calculating these metrics for each genus, we estimated group-level medians for bivalves, brachiopods, and cephalopods. Uncertainty in these medians was assessed using bootstrap resampling. For each group, genera were resampled with replacement 100 times, and the median value of each metric was recalculated for each bootstrap replicate. We report the median of these 100 bootstrap medians, together with the 2.5% and 97.5% quantiles as the 95% bootstrap interval.Contrast between regional and global ratesTo assess differences between global and regional diversification rates, we analysed the occurrences of marine bivalves, a well-established model system for macroevolution92,94,95. We focused on extinction because global mass extinction events provide a reference framework for comparing regional and global signals.For the global dataset, we calculated per-capita extinction rates96 using the divDyn package52. For regional datasets, we calculated three analogous rates: (1) regional disappearance rate, in which all regional last appearances were considered without distinguishing final extinction from extirpation; (2) regional extinction rate, in which only regional last appearances that coincided with the global final extinction of the genus were considered22; and (3) regional extirpation rate, in which only regional last appearances that preceded the global final extinction of the genus were considered. For example, the regional per-capita extirpation rate in a given stage was calculated as (-{mathrm{ln}}(frac{{N}_{{through}}}{{N}_{{extirpation}}+{N}_{{through}}})), where Nthrough is the number of genera that have occurrences on the focal plate before and after the focal stage, and Nextirpation is the number of genera that occurred on the focal plate before the focal stage, had their last occurrence on that plate within the focal stage, but had a last occurrence on other plates after the focal stage. We highlighted the rate variation from the Late Permian to earliest Jurassic, which encompasses two of the five mass extinctions3 (end-Permian mass extinction and end-Triassic mass extinction), on three plates with the most abundant fossil occurrences during this interval: Siberian craton, South China plate, and Adriatic (Apulian) plate in southern Europe. For each of the three regions, we also calculated the number of fossil-bearing hexagonal cells (with edge lengths of ~123 km) in each stage.Statistics and reproducibilityAll analyses were conducted in R97 (v4.5.3). After data cleaning, the datasets comprised 157,231 occurrences representing 2,472 bivalve genera, 130,418 occurrences representing 3,324 brachiopod genera, and 89,936 occurrences representing 3,679 cephalopod genera. For calculations of regional event proportions, only plate-stage combinations containing more than 100 occurrences were included to reduce noise associated with small sample sizes. One-sided Mann-Whitney U tests were used to determine whether global sampling rates and three-timer completeness were greater than their regional counterparts. These metrics were also summarised using medians, with 95% bootstrap intervals calculated from 100 bootstrap replicates. Spearman’s rank correlation tests were used to assess the relationships of the plate-level proportion of range shifts with the proportion of endemic taxa, global origination and extinction rates, and plate size. The same test was used to examine the relationship between the proportion of sampling-associated events and plate size. Boxplots were used to compare genus durations, geographical ranges, and taxonomic sampling rates among bivalves, brachiopods, and cephalopods. For genus-level duration estimates, occurrence ages were randomly sampled within stratigraphic stage boundaries 100 times, and the median estimate was retained. Uncertainty in group-level medians of genus duration, geographical range, and taxonomic sampling rate was assessed by resampling genera with replacement 100 times and reporting the median bootstrap estimate together with the 2.5% and 97.5% quantiles. All figures and tables can be reproduced using the data and code provided. In addition to the packages described above, packages used for data manipulation and visualisation included chronosphere98(v0.6.1), scales99 (v1.4.0), dplyr100 (v1.2.1), and stringr101 (v1.6.0).Reporting summaryFurther information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

    Data availability

    Fossil occurrence data were downloaded from the Paleobiology Database (paleobiodb.org). All datasets analysed in this study are available on Zenodo102.
    Code availability

    The code used for data analysis and figure generation is available on Zenodo102.
    ReferencesSepkoski, J. J. Jr., Bambach, R. K., Raup, D. M. & Valentine, J. W. Phanerozoic marine diversity and the fossil record. Nature 293, 435–437 (1981).Article 

    Google Scholar 
    Sepkoski, J. J. Jr. A factor analytic description of the Phanerozoic marine fossil record. Paleobiology 7, 36–53 (1981).
    Google Scholar 
    Raup, D. M. & Sepkoski, J. J. Jr. Mass extinctions in the marine fossil record. Science 215, 1501–1503 (1982).PubMed 

    Google Scholar 
    Alroy, J. Dynamics of origination and extinction in the marine fossil record. Proc. Natl. Acad. Sci. USA 105, 11536–11542 (2008).PubMed 
    PubMed Central 

    Google Scholar 
    Alroy, J. et al. Phanerozoic trends in the global diversity of marine invertebrates. Science 321, 97–100 (2008).PubMed 

    Google Scholar 
    Alroy, J. The shifting balance of diversity among major marine animal groups. Science 329, 1191–1194 (2010).PubMed 

    Google Scholar 
    Cárdenas, A. L. & Harries, P. J. Effect of nutrient availability on marine origination rates throughout the Phanerozoic eon. Nat. Geosci. 3, 430–434 (2010).Article 

    Google Scholar 
    Hannisdal, B. & Peters, S. E. Phanerozoic earth system evolution and marine biodiversity. Science 334, 1121–1124 (2011).PubMed 

    Google Scholar 
    Mayhew, P. J., Bell, M. A., Benton, T. G. & McGowan, A. J. Biodiversity tracks temperature over time. Proc. Natl. Acad. Sci. USA 109, 15141–15145 (2012).PubMed 
    PubMed Central 

    Google Scholar 
    Zaffos, A., Finnegan, S. & Peters, S. E. Plate tectonic regulation of global marine animal diversity. Proc. Natl. Acad. Sci. USA 114, 5653–5658 (2017).PubMed 
    PubMed Central 

    Google Scholar 
    Malanoski, C. M., Farnsworth, A., Lunt, D. J., Valdes, P. J. & Saupe, E. E. Climate change is an important predictor of extinction risk on macroevolutionary timescales. Science 383, 1130–1134 (2024).PubMed 

    Google Scholar 
    Jablonski, D., Roy, K. & Valentine, J. W. Out of the tropics: evolutionary dynamics of the latitudinal diversity gradient. Science 314, 102–106 (2006).PubMed 

    Google Scholar 
    Benson, R. B. J., Butler, R., Close, R. A., Saupe, E. & Rabosky, D. L. Biodiversity across space and time in the fossil record. Curr. Biol. 31, R1225–R1236 (2021).PubMed 

    Google Scholar 
    Close, R. A., Benson, R. B. J., Saupe, E. E., Clapham, M. E. & Butler, R. J. The spatial structure of Phanerozoic marine animal diversity. Science 368, 420–424 (2020).PubMed 

    Google Scholar 
    Close, R. A., Benson, R. B. J., Kiessling, W. & Saupe, E. E. Reefal regions were biodiversity hotspots throughout the Phanerozoic. Sci. Adv. 11, eadv9793 (2025).PubMed 
    PubMed Central 

    Google Scholar 
    Flannery-Sutherland, J. T. et al. Late Cretaceous ammonoids show that drivers of diversification are regionally heterogeneous. Nat. Commun. 15, 5382 (2024).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Flannery-Sutherland, J. T., Silvestro, D. & Benton, M. J. Global diversity dynamics in the fossil record are regionally heterogeneous. Nat. Commun. 13, 2751 (2022).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Vermeij, G. J. & Leighton, L. R. Does global diversity mean anything? Paleobiology 29, 3–7 (2003).
    Google Scholar 
    Antell, G. T., Benson, R. B. J. & Saupe, E. E. Spatial standardization of taxon occurrence data—a call to action. Paleobiology 50, 177–193 (2024).
    Google Scholar 
    Dillon, E. M. et al. Challenges and directions in analytical paleobiology. Paleobiology 49, 377–393 (2023).PubMed 
    PubMed Central 

    Google Scholar 
    Saupe, E. E. Explanations for latitudinal diversity gradients must invoke rate variation. Proc. Natl. Acad. Sci. USA 120, e2306220120 (2023).PubMed 
    PubMed Central 

    Google Scholar 
    Allen, B. J. et al. Estimating spatial variation in origination and extinction in deep time: a case study using the Permian–Triassic marine invertebrate fossil record. Paleobiology 49, 509–526 (2023).
    Google Scholar 
    Allen, B., Warnock, R. & Dunhill, A. M. “A history of the world imperfectly kept”: Will we ever know how biodiversity has changed over deep time? Preprint at https://doi.org/10.32942/X2DD1V (2025).Raup, D. M. Taxonomic diversity during the Phanerozoic: the increase in the number of marine species since the Paleozoic may be more apparent than real. Science 177, 1065–1071 (1972).PubMed 

    Google Scholar 
    Raup, D. M. Species diversity in the Phanerozoic: an interpretation. Paleobiology 2, 289–297 (1976).
    Google Scholar 
    Peters, S. E. & Foote, M. Biodiversity in the Phanerozoic: a reinterpretation. Paleobiology 27, 583–601 (2001).
    Google Scholar 
    Crampton, J. S. et al. Estimating the rock volume bias in paleobiodiversity studies. Science 301, 358–360 (2003).PubMed 

    Google Scholar 
    Smith, A. B. & McGowan, A. J. Cyclicity in the fossil record mirrors rock outcrop area. Biol. Lett. 1, 443–445 (2005).PubMed 
    PubMed Central 

    Google Scholar 
    Smith, A. B. & McGowan, A. J. The shape of the Phanerozoic marine palaeodiversity curve: how much can be predicted from the sedimentary rock record of Western Europe? Palaeontology 50, 765–774 (2007).
    Google Scholar 
    Barrett, P. M., McGowan, A. J. & Page, V. Dinosaur diversity and the rock record. Proc. R. Soc. B: Biol. Sci. 276, 2667–2674 (2009).
    Google Scholar 
    Wall, P. D., Ivany, L. C. & Wilkinson, B. H. Revisiting Raup: exploring the influence of outcrop area on diversity in light of modern sample-standardization techniques. Paleobiology 35, 146–167 (2009).
    Google Scholar 
    Butler, R. J., Barrett, P. M., Nowbath, S. & Upchurch, P. Estimating the effects of sampling biases on pterosaur diversity patterns: implications for hypotheses of bird/pterosaur competitive replacement. Paleobiology 35, 432–446 (2009).
    Google Scholar 
    Benson, R. B. J., Butler, R. J., Lindgren, J. & Smith, A. S. Mesozoic marine tetrapod diversity: mass extinctions and temporal heterogeneity in geological megabiases affecting vertebrates. Proc. R. Soc. B: Biol. Sci. 277, 829–834 (2010).
    Google Scholar 
    Na, L., Li, Q., Krause, C., Zhu, M. & Kiessling, W. Revisiting the Phanerozoic rock–diversity relationship. Geol. Mag. 160, 1853–1862 (2023).
    Google Scholar 
    Holland, S. M. The stratigraphic distribution of fossils. Paleobiology 21, 92–109 (1995).
    Google Scholar 
    Smith, A. B., Gale, A. S. & Monks, N. E. Sea-level change and rock-record bias in the Cretaceous: a problem for extinction and biodiversity studies. Paleobiology 27, 241–253 (2001).
    Google Scholar 
    Smith, A. B. Large–scale heterogeneity of the fossil record: implications for Phanerozoic biodiversity studies. Philos. Trans. R. Soc. Lond. Ser. B: Biol. Sci. 356, 351–367 (2001).
    Google Scholar 
    Peters, S. E. & Foote, M. Determinants of extinction in the fossil record. Nature 416, 420–424 (2002).Article 
    PubMed 

    Google Scholar 
    Peters, S. E. Geologic constraints on the macroevolutionary history of marine animals. Proc. Natl. Acad. Sci. USA 102, 12326–12331 (2005).PubMed 
    PubMed Central 

    Google Scholar 
    Patzkowsky, M. E. & Holland, S. M. Stratigraphic Paleobiology: Understanding the Distribution of Fossil Taxa in Time and Space. (University of Chicago Press, 2012).Foote, M. & Raup, D. M. Fossil preservation and the stratigraphic ranges of taxa. Paleobiology 22, 121–140 (1996).PubMed 

    Google Scholar 
    Foote, M. & Sepkoski, J. J. Jr. Absolute measures of the completeness of the fossil record. Nature 398, 415–417 (1999).Article 
    PubMed 

    Google Scholar 
    Guo, Z., Flannery-Sutherland, J. T., Benton, M. J. & Chen, Z. Q. Bayesian analyses indicate bivalves did not drive the downfall of brachiopods following the Permian-Triassic mass extinction. Nat. Commun. 14, 5566 (2023).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Tian, S. Y. et al. Cenozoic history of the tropical marine biodiversity hotspot. Nature 632, 343–349 (2024).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Stanley, S. M. An analysis of the history of marine animal diversity. Paleobiology 33, 1–55 (2007).
    Google Scholar 
    Button, D. J., Lloyd, G. T., Ezcurra, M. D. & Butler, R. J. Mass extinctions drove increased global faunal cosmopolitanism on the supercontinent Pangaea. Nat. Commun. 8, 733 (2017).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Dai, X. & Song, H. Toward an understanding of cosmopolitanism in deep time: a case study of ammonoids from the middle Permian to the Middle Triassic. Paleobiology 46, 533–549 (2020).
    Google Scholar 
    Powell, M. G., Moore, B. R. & Smith, T. J. Origination, extinction, invasion, and extirpation components of the brachiopod latitudinal biodiversity gradient through the Phanerozoic Eon. Paleobiology 41, 330–341 (2015).
    Google Scholar 
    Carrillo, J. D. et al. Disproportionate extinction of South American mammals drove the asymmetry of the Great American Biotic Interchange. Proc. Natl. Acad. Sci. USA 117, 26281–26287 (2020).PubMed 
    PubMed Central 

    Google Scholar 
    Fenton, I. S., Aze, T., Farnsworth, A., Valdes, P. & Saupe, E. E. Origination of the modern-style diversity gradient 15 million years ago. Nature 614, 708–712 (2023).Article 
    PubMed 

    Google Scholar 
    Raup, D. M. & Boyajian, G. E. Patterns of generic extinction in the fossil record. Paleobiology 14, 109–125 (1988).PubMed 

    Google Scholar 
    Kocsis, ÁT., Reddin, C. J., Alroy, J. & Kiessling, W. The R package divDyn for quantifying diversity dynamics using fossil sampling data. Methods Ecol. Evolution 10, 735–743 (2019).
    Google Scholar 
    Chen, Z. Q. & Benton, M. J. The timing and pattern of biotic recovery following the end-Permian mass extinction. Nat. Geosci. 5, 375–383 (2012).Article 

    Google Scholar 
    Dal Corso, J. et al. Environmental crises at the Permian–Triassic mass extinction. Nat. Rev. Earth Environ. 3, 197–214 (2022).Article 

    Google Scholar 
    Schoepfer, S. D., Algeo, T. J., van de Schootbrugge, B. & Whiteside, J. H. The Triassic–Jurassic transition – A review of environmental change at the dawn of modern life. Earth-Sci. Rev. 232, 104099 (2022).
    Google Scholar 
    Dal Corso, J. et al. Extinction and dawn of the modern world in the Carnian (Late Triassic). Sci. Adv. 6, eaba0099 (2020).
    Google Scholar 
    Roden, V. J. et al. Fossil liberation: a model to explain high biodiversity in the Triassic Cassian Formation. Palaeontology 63, 85–102 (2020).
    Google Scholar 
    Viaretti, M., Muttoni, G. & Angiolini, L. Late Palaeozoic brachiopod diversity: global vs regional trends. Evol. Earth 4, 100145 (2026).
    Google Scholar 
    Antell, G. T., Kiessling, W., Aberhan, M. & Saupe, E. E. Marine biodiversity and geographic distributions are independent on large scales. Curr. Biol. 30, 115–121 (2020).PubMed 

    Google Scholar 
    Lidgard, S., Di Martino, E., Zágoršek, K. & Liow, L. H. When fossil clades ‘compete’: local dominance, global diversification dynamics and causation. Proc. R. Soc. B: Biol. Sci. 288, 20211632 (2021).
    Google Scholar 
    Valentine, J. W. & Moores, E. M. Plate-tectonic regulation of faunal diversity and sea level: a model. Nature 228, 657–659 (1970).Article 
    PubMed 

    Google Scholar 
    Valentine, J. W., Foin, T. C. & Peart, D. A provincial model of Phanerozoic marine diversity. Paleobiology 4, 55–66 (1978).
    Google Scholar 
    Cermeño, P. et al. Post-extinction recovery of the Phanerozoic oceans and biodiversity hotspots. Nature 607, 507–511 (2022).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Kocsis, ÁT., Reddin, C. J., Scotese, C. R., Valdes, P. J. & Kiessling, W. Increase in marine provinciality over the last 250 million years governed more by climate change than plate tectonics. Proc. R. Soc. B: Biol. Sci. 288, 20211342 (2021).
    Google Scholar 
    Silvestro, D. et al. Fossil biogeography: a new model to infer dispersal, extinction and sampling from palaeontological data. Philos. Trans. R. Soc. B: Biol. Sci. 371, 20150225 (2016).
    Google Scholar 
    Hauffe, T., Pires, M. M., Quental, T. B., Wilke, T. & Silvestro, D. A quantitative framework to infer the effect of traits, diversity and environment on dispersal and extinction rates from fossils. Methods Ecol. Evolution 13, 1201–1213 (2022).
    Google Scholar 
    Dunhill, A. M., Benton, M. J., Twitchett, R. J. & Newell, A. J. Testing the fossil record: Sampling proxies and scaling in the British Triassic–Jurassic. Palaeogeography, Palaeoclimatology, Palaeoecology 404, 1–11 (2014).Dunhill, A. M., Hannisdal, B. & Benton, M. J. Disentangling rock record bias and common-cause from redundancy in the British fossil record. Nat. Commun. 5, 4818 (2014).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Alroy, J. A more precise speciation and extinction rate estimator. Paleobiology 41, 633–639 (2015).
    Google Scholar 
    Alroy, J. Accurate and precise estimates of origination and extinction rates. Paleobiology 40, 374–397 (2014).
    Google Scholar 
    Liow, L. H. & Nichols, J. D. Estimating rates and probabilities of origination and extinction using taxonomic occurrence data: capture-mark-recapture (CMR) approaches. Paleontological Soc. Pap. 16, 81–94 (2010).
    Google Scholar 
    Silvestro, D., Salamin, N., Antonelli, A. & Meyer, X. Improved estimation of macroevolutionary rates from fossil data using a Bayesian framework. Paleobiology 45, 546–570 (2019).
    Google Scholar 
    Silvestro, D., Salamin, N. & Schnitzler, J. PyRate: a new program to estimate speciation and extinction rates from incomplete fossil data. Methods Ecol. Evolution 5, 1126–1131 (2014).
    Google Scholar 
    Silvestro, D., Schnitzler, J., Liow, L. H., Antonelli, A. & Salamin, N. Bayesian estimation of speciation and extinction from incomplete fossil occurrence data. Syst. Biol. 63, 349–367 (2014).PubMed 
    PubMed Central 

    Google Scholar 
    Close, R. A. et al. The apparent exponential radiation of Phanerozoic land vertebrates is an artefact of spatial sampling biases. Proc. R. Soc. B: Biol. Sci. 287, 20200372 (2020).
    Google Scholar 
    Close, R. A., Benson, R. B. J., Upchurch, P. & Butler, R. J. Controlling for the species-area effect supports constrained long-term Mesozoic terrestrial vertebrate diversification. Nat. Commun. 8, 15381 (2017).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Vilhena, D. A. & Smith, A. B. Spatial bias in the marine fossil record. PLoS One 8, e74470 (2013).PubMed 
    PubMed Central 

    Google Scholar 
    Phillipi, D. B., Czekanski-Moir, J. & Ivany, L. C. Global Phanerozoic biodiversity—can variation be explained by spatial sampling intensity? Paleobiology 50, 661–672 (2024).
    Google Scholar 
    Ezard, T. H., Aze, T., Pearson, P. N. & Purvis, A. Interplay between changing climate and species’ ecology drives macroevolutionary dynamics. Science 332, 349–351 (2011).PubMed 

    Google Scholar 
    Mannion, P. D. The spatiotemporal distribution of Mesozoic dinosaur diversity. Biol. Lett. 20, 20240443 (2024).PubMed 
    PubMed Central 

    Google Scholar 
    Tomašových, A. & Jablonski, D. Decoupling of latitudinal gradients in species and genus geographic range size: a signature of clade range expansion. Glob. Ecol. Biogeogr. 26, 288–303 (2017).
    Google Scholar 
    Hendricks, J. R., Saupe, E. E., Myers, C. E., Hermsen, E. J. & Allmon, W. D. The generification of the fossil record. Paleobiology 40, 511–528 (2014).
    Google Scholar 
    Forey, P. L., Fortey, R. A., Kenrick, P. & Smith, A. B. Taxonomy and fossils: a critical appraisal. Philos. Trans. R. Soc. Lond. Ser. B: Biol. Sci. 359, 639–653 (2004).
    Google Scholar 
    Raup, D. M. Biases in the fossil record of species and genera. Bull. Carnegie Mus. Nat. Hist. 13, 85–91 (1979).
    Google Scholar 
    Wagner, P. J., Aberhan, M., Hendy, A. & Kiessling, W. The effects of taxonomic standardization on sampling-standardized estimates of historical diversity. Proc. R. Soc. Lond., Ser. B: Biol. Sci. 274, 439–444 (2007).
    Google Scholar 
    Prothero, D. Garbage in, garbage out: The effect of immature taxonomy on database compilations of North American fossil mammals. N. Mex. Mus. Nat. Hist. Sci. Bull. 67, 257–264 (2015).
    Google Scholar 
    Flannery-Sutherland, J. T., Raja, N. B., Kocsis, Á. T. & Kiessling, W. fossilbrush: An R package for automated detection and resolution of anomalies in palaeontological occurrence data. Methods Ecol. Evolution 13, 2404–2418 (2022).
    Google Scholar 
    Müller, R. D. et al. GPlates: Building a virtual Earth through deep time. Geochem., Geophysics, Geosystems 19, 2243–2261 (2018).
    Google Scholar 
    Kocsis, Á. T., Raja, N. B., Williams, S. & Dowding, E. M. rgplates: R interface for the GPlates Web Service and Desktop Application https://doi.org/10.5281/zenodo.8093990 (2024).Scotese, C. R. PALEOMAP PaleoAtlas for GPlates and the PaleoData Plotter Program www.earthbyte.org/paleomap-paleoatlas-for-gplates (2016).Harper, D. A. et al. Biodiversity, biogeography and phylogeography of Ordovician rhynchonelliform brachiopods. Geol. Soc., Lond., Mem. 38, 127–144 (2013).
    Google Scholar 
    Jablonski, D., Roy, K., Valentine, J. W., Price, R. M. & Anderson, P. S. The impact of the pull of the recent on the history of marine diversity. Science 300, 1133–1135 (2003).PubMed 

    Google Scholar 
    Kocsis, ÁT. icosa: Global Triangular and Penta-Hexagonal Grids Based on Tessellated Icosahedra https://CRAN.R-project.org/package=icosa (2021).Miller, A. I. & Sepkoski, J. J. Jr. Modeling bivalve diversification: the effect of interaction on a macroevolutionary system. Paleobiology 14, 364–369 (1988).PubMed 

    Google Scholar 
    Foote, M., Edie, S. M. & Jablonski, D. Ecological structure of diversity-dependent diversification in Phanerozoic marine bivalves. Biol. Lett. 20, 20230475 (2024).PubMed 
    PubMed Central 

    Google Scholar 
    Foote, M. Origination and extinction components of taxonomic diversity: general problems. Paleobiology 26, 74–102 (2000).
    Google Scholar 
    R. Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria https://www.R-project.org (2026).Kocsis, Á. T. & Raja, N. B. chronosphere: Evolving Earth System Variables https://doi.org/10.5281/zenodo.3525481 (2023).Wickham, H., Pedersen, T. & Seidel, D. scales: Scale Functions for Visualization https://doi.org/10.32614/CRAN.package.scales (2025).Wickham, H., François, R., Henry, L., Müller, K. & Vaughan, D. dplyr: A Grammar of Data Manipulation https://doi.org/10.32614/CRAN.package.dplyr (2026).Wickham, H. stringr: Simple, Consistent Wrappers for Common String Operations https://doi.org/10.32614/CRAN.package.stringr (2025).Guo, Z. Supplementary data for “Global patterns are key to understanding macroevolution”. Zenodo https://doi.org/10.5281/zenodo.21812688 (2026).Download referencesAcknowledgementsWe are grateful to the contributions of the PBDB. This is PBDB publication 563.FundingZ.G. was supported by the Alexander von Humboldt Foundation. Open Access funding enabled and organized by Projekt DEAL.Author informationAuthors and AffiliationsGeoZentrum Nordbayern, Department of Geography and Geosciences, Friedrich-Alexander-Universität Erlangen-Nürnberg, Loewenichstraße 28, Erlangen, GermanyZhen Guo & Wolfgang KiesslingAuthorsZhen GuoView author publicationsSearch author on:PubMed Google ScholarWolfgang KiesslingView author publicationsSearch author on:PubMed Google ScholarContributionsConceptualization: Z.G., W.K. Formal analysis: Z.G. Writing: Z.G., W.K.Corresponding authorCorrespondence to
    Zhen Guo.Ethics declarations

    Competing interests
    The authors declare no competing interests.

    Peer review

    Peer review information
    Communications Biology thanks all reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Dr Katie (E) Davis and Dr David Favero.

    Additional informationPublisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary informationSupplementary Figs. and Tables (download PDF )Reporting Summary (download PDF )Rights and permissions
    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
    Reprints and permissionsAbout this articleCite this articleGuo, Z., Kiessling, W. Global patterns are key to understanding macroevolution.
    Commun Biol 9, 1313 (2026). https://doi.org/10.1038/s42003-026-11118-7Download citationReceived: 23 February 2026Accepted: 25 September 2026Published: 08 October 2026Version of record: 08 October 2026DOI: https://doi.org/10.1038/s42003-026-11118-7Share this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy shareable link to clipboard
    Provided by the Springer Nature SharedIt content-sharing initiative More

  • in

    Different dietary compositions shape reproductive physiology in Litopenaeus vannamei broodstock

    AbstractMaintaining the reproductive performance of Litopenaeus vannamei broodstock toward the end of the breeding season remains a major challenge for commercial hatcheries, primarily due to physiological exhaustion. This study evaluated whether diversified dietary formulations could extend the breeding period and enhance larval quality. Broodstock (males: 50 ± 2 g; females: 69 ± 5 g) were allocated to four dietary treatments comprising worms (Nereis spp., ragworm), squid (Loligo forbesii), and mussels (Mytilus galloprovincialis). Group 1: 80% worms and 20% squid; Group 2: 80% worms and 20% mussels; Group 3: 50% worms, 25% mussels, and 25% squid; and Group 4: equal proportions (25% each) of worms, mussels, squid, and chicken liver. Over a 30-day trial, reproductive indices (HSI, GSI, OMI), hemolymph biochemistry, hepatopancreatic antioxidant status, and the expression of key maturation-related genes (vih, vtg, iag, igfbp1) were assessed, alongside histological examination of the hepatopancreas and gonads. Group 3 exhibited the highest reproductive efficiency, producing significantly greater egg output and larval yield compared with the worm-dominated control (Group 1). Molecular analyses showed a pronounced downregulation of vitellogenesis-inhibiting hormone (vih) in the eyestalks, accompanied by upregulation of vitellogenin (vtg) in the hepatopancreas of Groups 3 and 4. In males, expression of the insulin-like androgenic gland hormone (iag) peaked in Group 4, corresponding with enhanced histological structure of the vas deferens. Group 4 also exhibited enhanced hepatopancreatic antioxidant status, characterized by increased superoxide dismutase activity and reduced lipid peroxidation. Histological observations showed that diversified diets promoted healthy hepatopancreatic structure and sustained vitellogenic progression in the ovaries. Overall, reducing reliance on polychaetes and incorporating a balanced mix of mussels, squid, and chicken liver provides an effective nutritional strategy to rejuvenate aging broodstock. These findings offer a practical nutritional approach for hatcheries to enhance seed production and potentially prolong broodstock productivity, with the observed reproductive responses associated with changes in neuroendocrine-related marker gene expression.

    Explore related subjects
    Discover the latest articles and news in related subjects.

    Ecology

    Physiology

    Zoology

    IntroductionOver the past several years, the global aquaculture industry has experienced significant progress, with whiteleg shrimp farming now accounting for more than 80% of total global penaeid shrimp production1,2. This growth has been driven largely by the combined influence of several key factors, most notably the adoption of improved strains of Pacific white shrimp (Litopenaeus vannamei). These advancements have resulted in higher survival rates, faster growth, reduced production costs, and greater resistance to disease3,4.As a consequence, demand for shrimp has risen substantially across both domestic and international markets, strengthening the industry’s role as a highly profitable and strategically important economic sector5,6. Over the past two decades, the cultivation of Pacific white shrimp has expanded rapidly, with the species being introduced to an increasing number of countries and its global production footprint continuing to grow2. This widespread adoption is largely due to its strong adaptability to a range of aquaculture environments, along with its relatively high tolerance to diseases and environmental stressors compared with other shrimp species7,8,9.This sustained growth has, in turn, intensified interest among researchers and investors in improving broodstock quality to maximize production. Commercial hatcheries rely on consistent broodstock performance to supply high-quality postlarvae for grow-out operations. However, maintaining reproductive performance over extended breeding seasons remains a major challenge10. As broodstock age and undergo repeated spawning cycles, physiological exhaustion, depleted energy reserves, oxidative stress, and endocrine dysregulation can reduce fecundity, egg quality, and larval viability—ultimately lowering seed output and increasing production costs11,12.Nutrition plays a central role in reproductive success in crustaceans. It provides the substrates required for gametogenesis, supports the synthesis of lipoproteins and yolk precursors, modulates antioxidant defenses, and interacts with neuroendocrine pathways that regulate maturation13,14. In penaeid shrimp, the hepatopancreas serves as the primary metabolic organ, storing and converting dietary lipids and proteins into vitellogenin and other compounds essential for oocyte development15. Consequently, diet strongly influences key reproductive indices—hepatosomatic index (HSI), gonadosomatic index (GSI), and ovary maturation index (OMI)—as well as molecular regulators such as vitellogenin (vtg), insulin-like androgenic gland hormone (iag), and neuropeptides including vitellogenesis-inhibiting hormone (vih), produced in the eyestalk X-organ–sinus gland complex16,17,18.Although formulated broodstock diets have been developed and are increasingly incorporated into maturation protocols, commercial shrimp hatcheries continue to rely heavily on fresh or frozen marine ingredients19,20. Traditionally, polychaete worms have been widely used as a live or fresh feed for broodstock due to their high protein and essential fatty acid content, as well as their ability to stimulate spawning and improve larval quality21,22. However, prolonged reliance on a single feed source can be limiting. Monodiets often fail to provide a balanced supply of long-chain polyunsaturated fatty acids (LC-PUFAs), sterols, phospholipids, micronutrients, and other conditionally essential compounds required for sustained vitellogenesis and antioxidant defense. This imbalance can lead to physiological depletion and reduced reproductive output, necessitating dietary supplementation23,24.Among the fresh feed ingredients available for broodstock nutrition, mussels, squid, and chicken liver are of particular interest because their distinct nutrient profiles may address different physiological requirements associated with gonadal maturation and repeated spawning. Mussels are rich in phospholipids and micronutrients that may support hepatopancreatic function and vitellogenesis, whereas squid provides highly digestible protein and docosahexaenoic acid (DHA), both of which contribute to egg membrane integrity and larval neural development25,26,27. Chicken liver is a concentrated source of cholesterol, retinol, and B vitamins, which play important roles in steroidogenesis and energy metabolism during successive spawning cycles28,29. Chicken liver and enzymatically processed chicken liver powder have also been evaluated as alternatives to fishmeal in diets for juveniles Pacific white shrimp, with studies reporting improvements in growth and immune-related parameters30,31. Compared with diets based on a single fresh feed ingredient, combining these nutritionally distinct sources may provide a more balanced supply of essential proteins, lipids, vitamins, and micronutrients. This complementary nutritional profile provided the rationale for evaluating mussels, squid, and chicken liver as dietary components for sustaining broodstock reproductive performance throughout an extended breeding period.Beyond macronutrient supply, diet also influences oxidative status and the ability of broodstock to withstand oxidative stress associated with intense metabolic activity during gametogenesis and spawning32. Antioxidant enzymes such as superoxide dismutase (SOD), catalase (CAT), and glutathione peroxidase (GPx) play a critical role in neutralizing reactive oxygen species (ROS)33,34. These ROS may be generated during yolk formation and embryogenesis32,35. Insufficient antioxidant protection can lead to lipid peroxidation in oocyte membranes and reduced larval viability32,35. Diets enriched with antioxidant vitamins (A and E), selenium, and related cofactors therefore have strong potential to preserve gamete integrity and extend the reproductive lifespan of broodstock36. Consistent evidence in shrimp and other crustaceans links enhanced dietary antioxidant capacity with improved egg quality and larval survival36,37.Reproductive regulation in crustaceans is governed by a tightly integrated endocrine axis that links environmental and nutritional cues to the molecular processes controlling maturation. The eyestalk-derived neuropeptide vih acts as a key inhibitor of vitellogenesis; reduced vih activity permits increased synthesis of vitellogenin, the primary yolk precursor17. In contrast, the androgenic gland hormone (iag) drives male sexual differentiation and supports reproductive function, including sperm maturation and mating success38,39. Nutritional inputs can influence the expression of these regulatory pathways, modulating hormones such as estradiol, vitellogenin, luteinizing hormone, follicle-stimulating hormone, and testosterone40.Consequently, optimizing fresh-feed combinations remains highly relevant to commercial broodstock management. The present study evaluated whether reducing reliance on polychaetes and introducing diversified diets containing different proportions of worms, mussels, squid, and chicken liver could extend the effective breeding period of L. vannamei broodstock toward the end of the reproductive season.Chicken liver was included as an exploratory, nutrient-dense animal ingredient to increase dietary diversity, given that terrestrial animal tissues have previously been used in shrimp maturation diets as supplementary sources of cholesterol, phospholipids, and micronutrients41. The 80:20 formulations maintained a consistent worm-based reference while allowing the effects of individual supplementary ingredients to be compared. In contrast, the 50:25:25 and equal 25% formulations represented intermediate and extensive levels of worm replacement, respectively. These treatments followed a stepwise ingredient-substitution approach commonly used in broodstock nutrition trials. Thus, the selected proportions established a controlled gradient of dietary diversification and ingredient replacement rather than representing previously established nutritional optima41.We hypothesized that mixed diets providing complementary sources of protein, long-chain polyunsaturated fatty acids, sterols, phospholipids, and antioxidants would: (1) sustain or improve the hepatosomatic index, gonadosomatic index, and overall reproductive performance; (2) maintain or improve hemolymph biochemical parameters related to lipid metabolism, energy status, and reproductive physiology; (3) enhance antioxidant enzyme activity while reducing lipid peroxidation; (4) shift neuroendocrine gene expression toward reproductive activation, as indicated by the downregulation of vih, upregulation of vtg, and modulation of iag and igfbp1 in males; and (5) support hepatopancreatic recovery and continued vitellogenic progression, as evidenced by ovarian histology, relative to a worm-dominated control diet. Demonstrating these effects would provide hatcheries with practical dietary strategies for sustaining the reproductive capacity of aging broodstock, increasing seed production, and reducing the operational constraints associated with physiological exhaustion during extended breeding cycles.Materials and methodsBroodstockAt the end of the hatchery breeding season, Whiteleg shrimp broodstock (L. vannamei; synonym: Penaeus vannamei) sourced from Kona Bay, USA, were selected to ensure comparable baseline reproductive performance among experimental groups. Selection was based on previous spawning records and similar reproductive maturity, the latter being assessed from external signs of ovarian development visible through the dorsal exoskeleton. Male broodstock had a mean body weight of 50± 2 g and a mean total length of 17 ± 2 cm, whereas females had a mean body weight of 69 ± 5 g and a mean total length of 19 ± 2 cm. The experiment was initiated at the end of the breeding season to evaluate whether different feeding strategies could prolong the reproductive performance of the broodstock beyond the normal breeding period.SamplingBefore tissue sampling, shrimp were anesthetized by immersion in a clove-oil solution at 450 mg/L until spontaneous movement ceased. The shrimp were subsequently transferred to an ice slurry to maintain insensibility, and lymph was collected prior to dissection. Terminal dissection and tissue collection were then performed immediately while the animals remained insensible, and the collected tissues were promptly processed for the subsequent analyses. The samples were collected at the beginning of the experiment (n = 3 per treatment) and at the end of the experiment (day 30; three per replicate, n = 9 per group).Initial measurementsAt the beginning of the study, initial samples were collected to establish baseline measurements for comparison with the final experimental results and to evaluate the effects of the dietary treatments on broodstock performance and reproductive outcomes. The broodstock used for the initial assessment were collected after completion of the final hatchery spawning cycle; therefore, the females were considered to be in a post-spawning condition. The body length and weight of each shrimp were recorded, and hemolymph samples were collected into citrate–EDTA anticoagulant solution for biochemical analysis. Hepatopancreas samples were collected from both males and females, while testes from males and ovaries from females were obtained for histological examination. In addition, eyestalks and hepatopancreas tissues from females and androgenic glands from males were collected, immediately preserved in liquid nitrogen, and subsequently used for RNA extraction.These analyses were conducted to assess the health status and overall growth of the broodstock. As well as, the measurements provided critical baseline data, reflecting both the condition of the stock at the end of the breeding season and the initial state of the broodstock used in the study. This information is essential for understanding the metabolic and health status of the shrimp under the experimental conditions.Experimental setupFollowing the initial investigation, shrimp were divided into four groups, each with three replicates. Each group consisted of 18 females and 18 males.The experiment was conducted in black concrete tanks measuring 6.3 × 6.3 × 1.2 m (length × width × height). Each tank had a central drainage outlet and separate seawater and freshwater inlets, allowing the salinity to be maintained at 30 ppt. Each tank was divided into three fully separated sections using impermeable partitions, resulting in three independent experimental units, each measuring 2.0 × 6.3 m. These units served as replicates for the respective treatment groups. Water was supplied independently to each experimental unit, thereby preventing water exchange among compartments. Before the experiment began, all tanks and equipment were disinfected with formalin at a concentration of 3 ppm.After adjusting the salinity, the water underwent a series of purification and sterilization processes to ensure high quality. This included filtration through a sand filter and a protein skimmer to remove organic compounds, followed by ozone treatment and UV sterilization, and finally passing through a cotton filter. Additionally, heaters and aerators were installed on all sides to maintain proper aeration, keeping the water temperature at 29 °C.Every day at 1:00 am, the tanks were cleaned and the water was changed completely by 100%. Light cycle was at 12 h of light and 12 h of darkness.Water quality was monitored periodically, with the following parameters recorded: salinity (30 ± 0.00 ppt), temperature (29 °C), pH (8.14 ± 0.02), ammonia nitrogen (0.12 ± 0.02 mg/L), nitrite (0.02 ± 0.01 mg/L), and total alkalinity (110 ± 0.00 mg/L).Weekly water samples were collected. Isolation and enumeration of Vibrio spp. were performed using thiosulfate–citrate–bile salts–sucrose (TCBS) agar. Microbial analyses included assessment of yellow-, green-, fluorescent-, and black colony-forming Vibrio species. These analyses were conducted to evaluate water quality and monitor shrimp health, ensuring optimal rearing conditions throughout the experiment.Feeding protocolsThe experimental design followed standard feeding practices commonly employed in shrimp hatcheries. During the acclimation period, the basal feeding ration was set at approximately 25% of the shrimp biomass, consisting of 80% worms (Nereis spp., ragworms) distributed over three daily meals and 20% squid (Loligo forbesii) provided in a single meal. This feeding regime is consistent with practices commonly reported for shrimp broodstock management in commercial hatcheries42,43,44. Throughout the experimental period, all feed offered was completely consumed by the shrimp.Four feeding treatments were evaluated: Group 1 (control): 80% worms and 20% squid (calamari); Group 2: 80% worms and 20% mussels (Mytilus galloprovincialis); Group 3: 50% worms, 25% mussels, and 25% squid (calamari); and Group 4: 25% worms, 25% mussels, 25% squid (calamari), and 25% chicken liver. All ingredient proportions were calculated on a wet-weight basis. The experimental feeding period lasted 30 days.Fresh squid was cleaned by removing the entrails and eyes and then cut into approximately 1 cm pieces. All feed ingredients were sterilized using povidone-iodine prior to use. The feed components were stored at − 20 °C until feeding, except for the worms, which were offered live.Composite analysis of feed and fatty acid analysis of feed sourcesFeed samples from each ingredient were chemically analyzed for moisture, protein, fat, and ash content using standard methods established by AOAC45. The crude protein was determined using Kjeldahl method, while the crude fat was measured by the Soxhlet extraction. Moisture content was estimated by calculating the weight loss observed after heating samples in an oven at 60 °C for two days until a constant weight was achieved. The ash contents were determined by analyzing the remaining content after heating samples in a muffle furnace at 550 °C for two hours.The fatty acid composition of the feed sources was determined based on data reported in a concurrent companion study (currently unpublished; Supplementary Table S1). Briefly, fat was extracted from each feed sample and subjected to fatty acid methyl ester (FAME) preparation following the method described by Danish and Nizami46.MatingOnce feeding stabilized on the basal diets—approximately two weeks into acclimation—the experimental diets were administered at full rations for 30 days. Upon observing signs of gonadal maturation in all treatments, we proceeded with the mating process and its preparations.To ensure that only fully mature females were introduced into the male tank for mating, ovarian development was assessed visually through the dorsal exoskeleton. Females were considered ready for mating when the ovary was markedly enlarged and turgid, extended from the cephalothoracic region along the dorsal abdomen, exhibited well-developed anterior and lateral lobes, and showed a distinct dark-orange to orange-brown coloration. These features corresponded to Stage 4, indicating full ovarian maturation and advanced yolk deposition. Females meeting these criteria were transferred to the male tank at 10:00 a.m., and the number introduced was recorded for accurate tracking.The male tank is continuously observed throughout the day to record key mating behaviors, such as circular swimming patterns47. At 7:00 PM, females are carefully examined for evidence of mating success. If sperm deposition in the seminal receptacle is confirmed, mated females are promptly transferred to the hatchery tank for egg incubation, as rapid transfer enhances fertilization success.Females that do not mate are returned to their original tanks to continue their reproductive cycle, allowing them to re-enter mating conditions without disruption. This meticulous monitoring process ensures that only fertilized females are moved to the hatchery.These steps were repeated six times for each group, involving six separate matings per treatment, all conducted on the same day. This design enabled a comprehensive assessment of each group’s impact on mating behavior and reproductive success. By repeating the procedures, random variation was minimized, enhancing result accuracy and improving the reliability of the study’s data.Reproductive performance assessmentThe hepatosomatic index (HSI) and gonadosomatic index (GSI) in females were calculated using the following formulas:$$:HSI=left(frac{text{Hepatopancreas}text{:weight}}{text{Total:body:weight}}right)times:100$$$$:GSI=left(frac{text{Gonad:weight}}{text{Total:body:weight}}right)times:100$$Additional reproductive parameters were evaluated separately for each experimental replicate. Females ready for mating were identified based on the external signs of advanced ovarian maturation described above and were recorded as the number of females reaching the appropriate mating stage before pairing. Fertilized females were defined as those that spawned eggs with the presence of an externally attached spermatophore on the female thelycum. Fertilization rate was calculated as follows: fertilization rate (%) = (number of fertilized females/number of females ready for mating) × 100.Egg production was estimated using a volumetric counting method. After spawning, the eggs were uniformly distributed within a known total volume of seawater, and three subsamples of known volume were collected and counted using a graduated glass counting vessel. The mean number of eggs per unit volume was then extrapolated to the total seawater volume to estimate the total number of eggs produced in each replicate. Hatched larvae were quantified at the nauplius stage using the same volumetric counting procedure. Egg production per female was calculated by dividing the total number of eggs by the number of spawning females, while larval production per female was calculated by dividing the total number of larvae by the number of spawning females. Total egg and larval production and their corresponding per-female values were expressed as ×10⁶.In addition, the ovarian maturation stages of female shrimp were recorded daily for 30 days across four dietary treatments to assess the duration of ovarian maturation at each stage. Ovarian-stage progression was evaluated by visually examining the ovaries through the dorsal exoskeleton. Females were classified into five developmental stages, numbered 0–4, based on ovarian visibility, size, coloration, and development of the anterior and lateral lobes. The staging system was adapted from previously published macroscopic ovarian maturity criteria for L. vannamei48,49, which have been validated against histological observations49. Stage 0 represented the absence of visible ovarian development. At Stage 1, the ovary appeared as a very thin, whitish to translucent line along the dorsal region and could be clearly distinguished from the digestive tract. Stage 2 was characterized by an increase in ovarian size and a change in coloration from whitish-translucent to pale yellow. At Stage 3, the ovary became more prominent and enlarged, with a distinct orange coloration and clearly developing anterior and lateral lobes, indicating advanced ovarian development and increased yolk accumulation. Stage 4 represented full ovarian maturation, with a markedly enlarged ovary extending from the cephalothoracic region along the dorsal abdomen, well-developed anterior and lateral lobes, and a dark-orange to orange-brown coloration, indicating advanced yolk deposition and readiness for spawning. The number of females assigned to each stage was recorded for each treatment, and ovarian-stage distribution was expressed as the percentage of examined females in each stage.Hemolymph biochemical parameters and hepatopancreatic oxidative stress biomarkersEstradiol (E2) concentrations were measured using a Siemens IMMULITE 2000/XPi immunoassay analyzer (Siemens Co., Germany). Cholesterol, low-density lipoprotein cholesterol (LDL), and triglyceride concentrations were determined colorimetrically using commercial assay kits (Biodiagnostic Co., Egypt) at wavelengths of 500, 550, and 505 nm, respectively. vLDL was estimated using the Friedewald equation: vLDL = triglycerides/550.Hepatopancreatic antioxidant and oxidative-stress indicators, including superoxide dismutase (SOD), glutathione peroxidase (GPx), malondialdehyde (MDA), and catalase (CAT) levels, were determined colorimetrically using commercial assay kits (Biodiagnostic Co., Egypt) at wavelengths of 560, 340, 534, and 510 nm, respectively. MDA was measured as an indicator of lipid peroxidation.Gene expression analysisRNA was extracted from eyestalks and hepatopancreas of females, as well as androgenic glands of males, using the TRIzol reagent (Invitrogen™ Life Technologies, US). The quantity and quality of the RNA were assessed using a NanoPhotometer® NP80 (Implen, Germany) and 1% agarose gel electrophoresis.cDNA was synthesized from 2 µg of RNA using a cDNA synthesis kit (Applied Biotechnology Co., Egypt). Real-time PCR was then performed using an Applied Biosystems 7900HT real-time PCR instrument with qPCR GreenMaster (Jena Bioscience, Germany). The cycling conditions were as follows: initial denaturation at 94 °C for 2 min, followed by 40 cycles of 94 °C for 30 s (denaturation), a gene-specific annealing temperature for 30 s, and 72 °C for 30 s (extension). A dissociation curve analysis was conducted to assess the specificity of the PCR products. The primers used are listed in Table 1.Table 1 Primers used for real-time PCR.Full size tableThe expression of vitellogenesis-inhibiting hormone (vih) was analyzed in the eyestalks, while vitellogenin (vtg) was examined in the hepatopancreas of females. In males, the expression of insulin-like androgenic gland hormone (iag) and insulin-like growth factor-binding protein-related protein 1 (igfbp1) was assessed in the androgenic glands. Gene expression levels were normalized to elongation factor 1-alpha (ef1α), which was used as the internal reference gene. Because no treatment group was designated as a calibrator, normalized relative gene expression was calculated as 2−ΔCt, where ΔCt = Ct(target gene) − Ct(reference gene). Higher 2−ΔCt values indicate greater normalized transcript abundance34.Histological examinationThe gonads and hepatopancreas of adult male and female shrimp were fixed in Davidson’s solution for 24 h, then dehydrated in a graded ethanol series (70%, 80%, 90%, and 100%), cleared in xylene, embedded in paraffin, and sectioned at a thickness of 4 μm. The tissue slides were stained with hematoxylin and eosin (H&E) and examined under a light microscope (OLYMPUS CKX53, Germany) equipped with a camera (OLYMPUS EP50, Germany) for histological evaluation. Images were captured for documentation53,54.Statistical analysisData were analyzed using GraphPad Prism software (version 4.0; GraphPad Software, San Diego, CA, USA; https://www.graphpad.com/). For each measured parameter, the dietary treatment was considered the independent factor, and differences among the four treatment groups were evaluated using one-way analysis of variance (ANOVA). The assumptions of normality and homogeneity of variance were assessed before analysis using the Shapiro–Wilk test and Levene’s test, respectively. When the overall ANOVA indicated a significant treatment effect, Tukey’s honestly significant difference multiple-comparison test was used to perform all pairwise comparisons among treatment means while controlling the family-wise Type I error rate. Results are presented as mean ± standard deviation (SD), with statistical significance defined as p < 0.05.In addition, ovarian maturation stages of female shrimp were recorded daily for 30 days across four dietary treatments (Groups 1–4). Each individual was classified into one of five ordinal developmental stages (Stages 0–4). Because maturation stage represents an ordered categorical variable, the effects of diet and time were analyzed using ordinal logistic regression (cumulative link model)55. The model was defined as:$$:logitleft[Pleft(Yge:jright)right]={alpha:}_{j}+{beta:}_{1}:left(Groupright)+{beta:}_{2}:left(Dayright)$$where Y represents the ovarian stage and αj denotes the threshold intercepts between stages. Counts of females at each stage were used as frequency weights56.To quantify overall ovarian development, an ovary maturity index (OMI) was calculated as the weighted mean maturity stage based on the distribution of females across stages, an approach conceptually similar to weighted maturity indices derived from the distribution of gonadal developmental stages57,58. Such indices provide a quantitative representation of reproductive development and complement traditional measures such as the gonadosomatic index59,60. The OMI was calculated as follows:$$:text{O}text{M}text{I}:=:frac{sum:(itimes{n}_{i})}{sum:{n}_{i}}$$where i represents the maturation stage (0–4) and ni is the number of females at that stage. These statistical analyses were performed using the statsmodels ordinal package in Python61, with significance set at α = 0.05.ResultsProximate composition and fatty acid profile of the dietsThe formulation in Group 4, consisting of equal parts (25%) worms, mussels, squid (calamari), and chicken liver, had the highest protein percentage, followed by Group 3, which contained 50% worms, 25% mussels, and 25% squid (calamari) (Table 2). The lowest protein percentage was found in Group 2, which included 80% worms and 20% mussels. Group 4 also had the highest crude fat percentage, while the other formulations showed similar levels. Ash content was highest in Group 2, followed by Group 1, Group 3, and Group 4 (Table 2).Table 2 Composite analysis of different diet formulations using worms, squid (calamari), mussels, and chicken liver.Full size tableThe diet formulation used in this study was identical to that described in a concurrent (unpublished) study; the complete fatty acid profile is provided in Supplementary Table S1 for reference. The fatty acid composition of the experimental diets varied significantly, reflecting differences in the proportions of marine and terrestrial ingredients.Group 1, which served as the control diet dominated by polychaetes, was characterized by higher levels of elaidic acid (approximately 1.2% higher than in Group 3), along with notable concentrations of heptadecanoic acid and cis-11,14,17-eicosatrienoic acid. Similar to Group 2, its profile reflects a traditional worm-based diet that, although rich in certain marine lipids, lacked the diversity observed in the higher-performing treatments. Group 2 contained the highest concentrations of palmitic acid, palmitoleic acid, and cis-5,8,11,14,17-eicosapentaenoic acid (EPA).In contrast, the more diversified diets—particularly Group 4—were distinguished by significantly higher levels of linoleic acid (ω−6), oleic acid, and cis-13,16-docosadienoic acid. Additionally, erucic acid, tetracosanoic acid, arachidic acid and cis-11-eicosenoic acid were markedly elevated in both Groups 3 and 4 compared with the worm-dominated diets (Groups 1 and 2). The diets of Groups 2 and 3 contained approximately 1% higher levels of α-linolenic acid (ω−3) than those of Groups 1 and 4. In contrast, the diet of Group 1 contained substantially lower levels of cis-4,7,10,13,16,19-docosahexaenoic acid (DHA) compared with the diets of the other groups.Fatty acids such as butyric, hexanoic, dodecanoic, and behenic acids were present only at negligible levels across all formulations.Female fecundityReproductive performance of shrimp females varied among the four dietary treatments (Group 1–Group 4) (Table 3). Diet significantly affected both the hepatosomatic index (HSI) and gonadosomatic index (GSI), with the highest values recorded in Group 3.Table 3 Female fecundity parameters of female shrimp fed four different diets.Full size tableDiet also significantly influenced the number of females available for reproduction (p = 0.043). Similarly, the number of fertilized females differed significantly among treatments (p = 0.019), with the highest values observed in Group 3 (3.00 ± 0.63) and Group 4 (2.83 ± 1.17). In contrast, fertilization rate did not differ significantly among dietary treatments (p = 0.062).Egg production varied markedly among treatments. Females fed the Group 3 diet produced the highest number of eggs, followed by Group 4, whereas the lowest egg production was observed in Group 2. Statistical analysis confirmed a significant effect of diet on total egg production (p = 0.013). Likewise, the number of eggs per female differed significantly among groups (p = 0.015). Group 3 exhibited a higher eggs-per-female ratio than Group 2, while Groups 1 and 4 showed intermediate values that did not differ significantly from the other treatments.Larval production followed a similar pattern. The highest larval output was recorded in Group 3 (10.34 ± 7.20 × 10⁶ larvae), which was significantly higher than that observed in Groups 1 and 2 (Table 3). Diet also had a significant effect on larval production (p = 0.0077). In addition, the larvae-per-female ratio varied significantly among treatments (p = 0.0015). Group 3 produced the highest number of larvae per female, significantly exceeding the values recorded for Groups 1 and 2. Larval production in Group 4 was intermediate and did not differ significantly from the other treatments.Ovarian development increased significantly over the 30-day experimental period across all treatment groups (β = 0.132, p < 0.001). Dietary treatment significantly influenced the maturation trajectory and the rate at which females transitioned from immature to mature stages (Fig. 1).Fig. 1Full size imageTemporal progression of the Ovary Maturity Index (OMI) in female shrimp fed four experimental diets over a 30-day period. The OMI was calculated as the weighted average of ovarian maturation stages (0–4), with higher values indicating more advanced maturation. Data points represent the daily mean OMI for each treatment group. Groups 3 (green) and 4 (red-orange) exhibited significantly higher maturation rates compared with the Group 1 control (blue) (p < 0.001). Group 2 (orange) showed the lowest overall maturation rate during the first half of the experimental period.Group 3 exhibited the strongest overall performance throughout the trial, with a significantly higher probability of progressing to advanced maturation stages compared with the Group 1 control (β = 0.658, p < 0.001). This response was characterized by an early transition out of Stage 0, as reflected in the distribution of ovarian developmental stages (Fig. 2). Similarly, Group 4 also outperformed the control group (β = 0.503, p < 0.001), showing a marked increase in reproductive readiness during the latter half of the experiment.Fig. 2Full size imageDistribution of ovarian developmental stages in female shrimp across four dietary treatments (Groups 1–4) during a 30-day feeding trial. Stacked bars show the percentage of individuals at each developmental stage: Stage 0 (dark purple; immature), Stage 1 (blue, early development), Stage 2 (teal, vitellogenesis), Stage 3 (green, late vitellogenesis), and Stage 4 (yellow; fully mature/spawning ready). A rapid transition from immature to advanced maturation stages (Stages 3–4) is most evident in Group 3 beginning around Day 13. By Day 30, Groups 3 and 4 displayed the highest proportion of reproductively ready females compared with the control (Group 1) and Group 2.In contrast, Group 2 showed a tendency toward slower maturation during the first two weeks of the study (β = −0.187), although it was not significantly different from the Group 1 control by the end of the trial (p = 0.132). Ovary Maturity Index (OMI) values calculated at Day 20 further supported these findings, with Groups 3 and 4 reaching an index of 2.50, substantially higher than the Group 1 baseline (OMI = 1.38) at the same time point (Fig. 1).Analysis of the stage distribution indicated that by the end of the trial (Day 30), Group 4 achieved the highest final maturity index (OMI = 2.56), while Group 3 maintained the most consistent proportion of Stage 4 (fully mature) females during the final 10 days of the experiment (Fig. 2).Hemolymph biochemistry and hepatopancreatic oxidative stress biomarkersFemalesIn female shrimp, hemolymph biochemical parameters exhibited variable responses among the experimental groups (Table 4). Estradiol (E2) concentrations were markedly lower than the initial values recorded before the experimental diets were introduced. However, Group 4 showed a slight but significant increase compared with the other treatment groups. Cholesterol levels did not differ significantly among the groups or from the initial measurements, with the exception of Group 2, which demonstrated a noticeable increase. Low-density lipoprotein (LDL) concentrations were higher in Groups 2 and 3; however, they remained lower than the initial measurements. In contrast, very-low-density lipoprotein (vLDL) reached its highest level in Group 4, followed by Group 3. Triglycerides displayed a trend similar to that observed for vLDL.Table 4 Hematological parameters and hepatopancreatic oxidative stress biomarkers in female shrimp fed four different diets.Full size tableRegarding hepatopancreatic oxidative stress biomarkers, superoxide dismutase (SOD) activity was significantly higher in Group 4, while no significant differences were detected among the remaining groups or compared with the initial values (Table 4). Glutathione peroxidase (GPx) activity was highest in Groups 2 and 3 and moderate in Group 4. Malondialdehyde (MDA) levels were markedly elevated in Group 2. In contrast, catalase activity did not differ significantly among the experimental groups (Table 4).MalesIn male shrimp, hemolymph biochemical profiles showed a different pattern compared with females (Table 5). E2 levels were lowest in Group 2 and highest in Group 3, while the remaining groups did not differ significantly from each other or from the initial measurements. Cholesterol and LDL concentrations were highest in Group 4 and lowest in Group 3. In contrast, vLDL levels peaked in Group 2, showed moderate reductions in Group 4, and reached their lowest values in Group 3. Triglyceride levels differed significantly among male groups and relative to the initial measurements, with the highest values in Group 2 and the lowest in Group 3. Group 4 also showed relatively high triglyceride levels; however, these did not differ significantly from those in Group 1.Table 5 Hematological parameters and hepatopancreatic oxidative stress biomarkers in male shrimp fed four different diets.Full size tableHepatopancreatic oxidative stress biomarkers varied among the groups. SOD activity was substantially reduced in Group 2 and was moderately lower in Group 3 compared with the other groups and the initial values (Table 5). Similarly, GPx activity was reduced in Groups 2 and 3 relative to the other treatments. MDA levels did not show marked variation among the male groups, although the lowest values were recorded in Group 1. In contrast, catalase activity increased in Groups 2 and 3 and reached its highest level in Group 4 (Table 5).Expression patterns of reproduction-related genesFemalesIn females, the expression of vitellogenesis-inhibiting hormone (vih) in the eyestalks showed a gradual downregulation from Group 1 to Group 4 compared with the initial stage prior to the introduction of the experimental diets (Fig. 3; p < 0.0001; Supplementary File 2). In contrast, vitellogenin (vtg) expression in the hepatopancreas increased progressively from Group 2 to Group 4, with no significant difference observed between the initial stage and Group 1 (Fig. 3; p < 0.0001; Supplementary File 2).Fig. 3Full size imageExpression of vitellogenesis-inhibiting hormone (vih) in the eyestalks and vitellogenin (vtg) in the hepatopancreas of female broodstock shrimp (L. vannamei) fed different experimental diets. Gene expression levels were normalized to elongation factor 1-alpha (ef1α) as the internal reference gene. Results are presented as mean ± SD (n = 9 per group). Relative expression was calculated using the ΔCt method and expressed as 2−ΔCt for clarity. The overall ANOVA indicated a significant effect of dietary treatment (P < 0.0001). Tukey-adjusted P values for the pairwise comparisons are shown between groups.MalesIn the androgenic glands, the expression levels of insulin-like androgenic gland hormone (iag) and insulin-like growth factor-binding protein–related protein 1 (igfbp1) showed a significant increase from Group 1 to Group 4 compared with the initial measurements. However, no significant difference was observed between Groups 3 and 4 (Fig. 4; Supplementary File 2).Fig. 4Full size imageExpression of insulin-like androgenic gland hormone (iag) and insulin-like growth factor-binding protein-related protein 1 (igfbp1) in the androgenic glands of male broodstock shrimp (L. vannamei) fed different experimental diets. Gene expression levels were normalized to elongation factor 1-alpha (ef1α) as the internal reference gene. Results are presented as mean ± SD (n = 9 per group). Relative expression was calculated using the ΔCt method and expressed as 2−ΔCt for clarity. The overall ANOVA indicated a significant effect of dietary treatment (P < 0.0001). Tukey-adjusted P values for the pairwise comparisons are shown between groups.Histological results of hepatopancreas, ovaries and vas deferensHepatopancreasAt the initial stage, the hepatopancreas showed preserved tubular architecture with distinct lumina and identifiable B, R, and F cells. B cells contained a prominent cytoplasmic vacuole, R cells displayed multiple smaller vacuoles, and F cells appeared comparatively basophilic and less vacuolated. No obvious epithelial disruption or cellular debris was observed within the tubule lumina (Fig. 5A).Fig. 5Full size imagePhotomicrographs of histological cross-sections of the hepatopancreas from experimental adult L. vannamei females (H&E staining). L, lumen of hepatopancreatic tubules; B, B-cells; F, F-cells; R, R-cells; E, E-cells.In Group 1, the hepatopancreatic epithelium exhibited prominent cytoplasmic vacuolation, particularly in regions containing R and B cells, together with irregular variation in lumen width (Fig. 5B). Group 2 showed generally rounded or oval tubules with distinct, frequently star-shaped lumina and numerous epithelial vacuoles. No clear epithelial necrosis or cellular sloughing was evident (Fig. 5C).Group 3 exhibited extensive epithelial vacuolation, variable epithelial thickness, and irregular tubule-lumen morphology, although overt epithelial destruction was not apparent (Fig. 5D). In Group 4, the tubules displayed well-defined star-shaped lumina surrounded by a relatively compact epithelial lining. Cytoplasmic vacuolation appeared less extensive than in Groups 1 and 3, and the overall tubular architecture was comparatively well preserved (Fig. 5E).In males at the initial stage, the hepatopancreas showed generally preserved tubule architecture with distinct lumina and identifiable B, F, and R cells. The epithelial cells displayed prominent cytoplasmic vacuolation, while the epithelial lining remained generally continuous, with no obvious cellular debris or extensive epithelial disruption within the lumina (Fig. 6A).Fig. 6Full size imagePhotomicrographs of histological cross-sections of the hepatopancreas from experimental adult L. vannamei males (H&E staining). L, lumen of hepatopancreatic tubules; B, B-cells; F, F-cells; R, R-cells; E, E-cells.In Group 1, the hepatopancreatic tubules retained their general organization and continuous epithelial lining. Prominent cytoplasmic vacuolation was observed in the epithelial cells, together with variation in lumen shape and width. However, an increase in vacuolation relative to the initial stage could not be established from the representative sections alone (Fig. 6B). Group 2 exhibited generally preserved tubule profiles with distinct, frequently irregular or star-shaped lumina. Numerous cytoplasmic vacuoles of variable size were present within the epithelium, but no clear evidence of extensive epithelial disruption or cellular sloughing was observed (Fig. 6C).In Group 3, the hepatopancreas showed heterogeneous tubule morphology, with prominent epithelial vacuolation and variation in lumen size and shape. Several tubules exhibited enlarged or irregular lumina, although definite epithelial thinning could not be confirmed consistently throughout the section (Fig. 6D). In Group 4, the tubules displayed relatively well-defined, frequently star-shaped lumina surrounded by a generally continuous epithelial lining. Cytoplasmic vacuolation remained evident, while the overall tubular architecture appeared comparatively preserved. A definite reduction in vacuolation relative to the other groups would require quantitative histomorphometric assessment (Fig. 6E).OvariesAt the initial stage, the ovarian tissue was composed predominantly of small primary-growth or previtellogenic oocytes, with some larger developing oocytes. The oocytes exhibited distinct nuclei and were supported by ovarian stroma containing visible blood vessels. No overt tissue disruption or degenerative changes were evident in the examined fields (Fig. 7A, B).Fig. 7Full size imagePhotomicrographs of histological cross-sections of ovaries from experimental L. vannamei females, showing multiple stages of oocyte maturation (H&E staining). P, perinucleolar stage; Pg, primary growth ovarian follicle; Lc, lipid and cortical alveoli stage; Vf, vitellogenic follicle; Ev, early vitellogenic stage; Lv, late vitellogenic stage; Pv, post-vitellogenic stage; Ca, cortical alveoli; post-ovulatory follicle (asterisk); yolk granules (arrowhead); N, nucleolus.In Group 1, the ovary remained dominated by small primary-growth and previtellogenic oocytes, together with a limited number of larger developing oocytes. Overall, the ovarian organization was generally preserved, but clear advancement of oocyte maturation relative to the initial stage was not evident (Fig. 7C, D). Group 2 displayed a greater occurrence of enlarged, vitellogenic oocytes with more abundant eosinophilic cytoplasm and distinct nuclei, while smaller primary-growth oocytes remained present. These features indicate progression of ovarian development relative to the initial stage and Group 1 (Fig. 7E, F).In Group 3, numerous large vitellogenic oocytes containing abundant yolk material were observed. The presence of prominent nuclei and cortical structures was consistent with a comparatively advanced stage of oocyte development (Fig. 7G, H). Group 4 contained a mixture of primary-growth, vitellogenic, and advanced or postvitellogenic oocytes, indicating continued ovarian maturation (Fig. 7I, J).vas deferensAt the initial stage, the vas deferens exhibits a healthy and fully functional organization, characterized by clear, unobstructed lumina, a delicate reticular framework, and well-defined basement membrane, submucosa, epithelial, and muscle layers. Cellular density is moderate, with evenly distributed nuclei. The epithelium consists of columnar cells with eosinophilic cytoplasm and prominent apical secretory material or brush border (arrowhead), indicating normal reproductive activity. The lumen is well defined and contains a dense sperm mass (Fig. 8A, B).Fig. 8Full size imagePhotomicrographs of histological cross-sections of the vas deferens from experimental adult L. vannamei males (H&E staining). Vd, vas deferens; SM, sperm mass; EPT, epithelial lining; MS, muscle layer.In Group 1, relative to the initial stage, the vas deferens shows slightly dilated lumina, subtle thickening of the submucosa, and minor architectural disturbances. Overall tissue organization remains largely preserved; however, the epithelial layer appears mildly reduced in height compared with the initial stage (Fig. 8C). The sperm mass remains compact. Fewer apical secretory granules (arrowhead) and a slightly thinner epithelium are evident (Fig. 8D). In Group 2, more pronounced alterations are observed compared with the initial stage. The lumina are well defined but variable in size, and cellular density is increased, with densely packed nuclei suggestive of epithelial hyperplasia. Epithelial cells display heightened secretory activity and irregular mucosal folding, indicating a proliferative response. Intense apical secretory material is evident (arrowhead; Fig. 8E, F).In Group 3, the vas deferens exhibits slightly constricted lumina, a reduction in sperm mass volume, and altered mucosal folding (Fig. 8G). Focal epithelial cell vacuolization is observed (arrowhead), and although the muscle layer remains present, overall mucosal integrity appears compromised (Fig. 8H). In Group 4, the overall architecture shows a moderately tall epithelium with abundant apical secretory material, resembling Group 2 but to a lesser extent. Columnar epithelial cells retain visible secretory granules (arrowhead; Fig. 8I, J). The sperm mass is present but appears less dense than that observed at the initial stage.DiscussionThe reproductive success of L. vannamei (whiteleg shrimp) is closely linked to broodstock nutrition because the diet provides the biochemical precursors required for vitellogenesis and spermatogenesis. In the present study, diversified feeding strategies were evaluated as a means of extending reproductive capacity toward the end of the commercial breeding season, a period typically characterized by physiological exhaustion and declining fecundity.The results demonstrate that dietary composition, particularly the inclusion and relative proportions of worms, mussels, squid, and chicken liver, substantially influenced reproductive performance. Females fed the Group 3 diet, consisting of 50% worms, 25% mussels, and 25% squid, exhibited the highest fecundity, egg production, larval output, and rate of ovarian maturation. The Group 4 diet, containing equal proportions of worms, mussels, squid, and chicken liver, also significantly improved reproductive readiness and maturity indices compared with the control and Group 2 diets. Proximate composition analysis further showed that Group 4 had the highest overall protein and lipid contents, indicating that nutrient-dense, diversified diets may help meet the considerable energetic and nutritional demands of gametogenesis.Collectively, these findings indicate that balanced combinations of marine-derived protein sources, with complementary lipid and fatty-acid profiles, can synergistically enhance reproductive performance, endocrine activity, and tissue integrity in broodstock, consistent with previous reports62,63. Although polychaetes (worms) are traditionally favored due to their high levels of highly unsaturated fatty acids (HUFAs)21,22,64,65, worm-dominated diets alone (Groups 1 and 2) did not yield optimal outcomes during the late breeding season. Instead, incorporating additional marine ingredients such as mussels and squid appears to improve the balance of essential fatty acids—particularly omega-3 to omega-6 ratios and levels of EPA and DHA—while also supplying key micronutrients, including sterols and phospholipids, that support vitellogenesis and ovarian maturation63. The superior performance of Group 3 suggests that maintaining a relatively high proportion of worms, while complementing them with other marine protein sources, optimizes reproductive output without introducing potential imbalances associated with terrestrial ingredients such as chicken liver, which is comparatively rich in polyunsaturated fatty acids and minerals28. Overall, these results support the concept that diversified broodstock diets can prolong reproductive efficiency and sustain spawning performance toward the end of the production cycle, in line with earlier findings in penaeid shrimp16,66,67.Interestingly, Group 4 exhibited significantly higher levels of linoleic and oleic acids. Beyond serving as energy sources, these fatty acids act as precursors for eicosanoid synthesis68,69,70, which regulates key reproductive processes in crustaceans. Notably, the elevated EPA levels observed in Group 2 did not translate into superior fecundity, underscoring that a balanced “cocktail” of lipid classes and amino acids—such as that provided by diversified protein sources—is more effective than a high concentration of any single fatty acid24,40.Proximate and fatty-acid analyses confirmed substantial variation among diets in palmitic, palmitoleic, linoleic, α-linolenic, EPA, and DHA contents, which likely underpin the observed reproductive differences. In penaeids, vitellogenesis is highly sensitive to lipid quality: long-chain polyunsaturated fatty acids (LC-PUFAs), phospholipids, and sterols are essential for oocyte membrane formation, yolk lipoprotein synthesis, and embryonic development13,71. Dietary provision of EPA and DHA has been consistently linked to improved egg quality and larval survival24,27. The combined presence of high-quality protein and a favorable fatty-acid profile in Groups 3 and 4 likely enhanced vitellogenin synthesis in the hepatopancreas and its transport to developing oocytes.The hepatosomatic index (HSI) and gonadosomatic index (GSI) are reliable indicators of nutritional status and reproductive maturation15,59. GSI increases progressively with ovarian development, and both indices are closely linked during vitellogenesis15,59. The peak HSI and GSI values observed in Group 3 indicate efficient mobilization of nutrients from the hepatopancreas to the ovaries, consistent with the elevated egg production and larval output recorded in this group.In line with these observations, the ovary maturity index (OMI) and ordinal logistic regression provided a high-resolution temporal profile of maturation. Groups 3 and 4 showed a rapid transition from immature (Stage 0) to advanced vitellogenic stages (Stages 3 and 4) by the second week of the trial. By Day 20, both groups reached an OMI of 2.50, significantly exceeding the control (OMI = 1.38). These findings suggest that diversified marine protein sources enhance nutrient availability for vitellogenesis and ovarian maturation, helping to overcome reproductive stagnation in aging broodstock and extend the effective spawning window14,72,73.Molecular responses provided additional insight into the dietary effects on female reproduction. Normalized vih transcript abundance in the eyestalk decreased from Group 2 to Group 4, whereas hepatopancreatic vtg abundance showed the opposite pattern. This inverse association is consistent with the established Vih–Vtg regulatory axis in L. vannamei, in which Vih produced by the eyestalk and brain suppresses hepatopancreatic vitellogenin synthesis17,59. Nutritional signals can modulate this axis through energy status, lipid-derived mediators, and nutrient-sensing pathways (e.g., insulin-like signaling), thereby relieving inhibition and promoting vitellogenesis63. Therefore, the lower vih and higher vtg transcript abundance in Groups 3 and 4 suggests reduced inhibitory signaling and greater capacity for yolk-precursor synthesis. Nutrient availability may also regulate vitellogenesis through metabolic signaling pathways72.The treatment comparisons indicate that molecular activation did not translate linearly into reproductive output. Group 2, which received worms and mussels, exhibited the lowest egg and larval production, consistent with its comparatively weaker vitellogenic response. Group 1, which received worms and squid, produced intermediate reproductive results. Group 3, which combined worms, mussels, and squid, achieved the highest GSI, number of fertilized females, and egg and larval production, suggesting that this combination provided the most effective nutritional balance for ovarian maturation and subsequent reproductive performance. Although Group 4 showed a strong molecular response, including lower vih and higher vtg transcript abundance, its reproductive outputs remained below those of Group 3. This discrepancy indicates that increased vtg transcript abundance alone does not ensure maximum fecundity and that nutrient balance, ovarian developmental stage, spawning timing, and post-transcriptional regulation may also influence reproductive performance.This response is likely mediated by diet-induced stimulation of the neuroendocrine system. For example, increased estradiol (E2) levels in Group 4 females coincided with higher vtg expression, consistent with the known role of E2 in promoting yolk protein synthesis in the hepatopancreas63,74,75. The ability of the Group 4 diet to modulate this axis late in the breeding season highlights its potential as a functional broodstock feed, rather than merely a source of energy.Male reproductive performance, although often overlooked, is equally dependent on nutrition. The androgenic gland regulates spermatogenesis and male differentiation through the secretion of insulin-like androgenic gland hormone (iag), which interacts with binding proteins such as igfbp18,76. The observed upregulation of iag and igfbp1 from Group 1 to Group 4 indicates that diet also enhances male reproductive physiology. Elevated iag expression is associated with increased spermatogenic activity, androgenic gland function, and the development of secondary sexual characteristics39. The higher expression levels in Groups 3 and 4 therefore suggest improved male condition and potentially enhanced mating success or sperm quality, contributing to the increased fertilization rates and larval outputs observed. The absence of significant differences between these two groups in male gene expression further suggests that, once a nutritional threshold is reached, activation of male reproductive pathways may plateau.Elevated iag expression in the dietary groups suggests that dietary composition influenced androgenic-gland activity at the transcriptional level. Given the established role of IAG in male sexual differentiation and reproductive regulation in crustaceans, this response may reflect modulation of male reproductive physiology77,78. This molecular response was accompanied by generally preserved vas deferens architecture across the dietary groups, including Groups 2 and 4, with continuous epithelial linings, visible apical epithelial material, and intraluminal sperm masses. These features suggest maintenance of vas deferens structure and epithelial function. The inclusion of chicken liver in the Group 4 diet—providing cholesterol, fat-soluble vitamins, and other nutrients—may have supplied substrates supporting androgenic gland hormone production and protein synthesis28,79.Hemolymph lipid indicators (cholesterol, LDL, vLDL, and triglycerides) varied with both diet and sex, providing insight into metabolic partitioning during reproductive preparation. Very-low-density lipoprotein (vLDL) plays a central role in yolk protein assembly and lipid transport to the ovaries in decapods29,80. Lipids more broadly are critical to reproductive success: triglycerides support molting and early embryonic development; phospholipids influence spawning frequency; EPA and DHA contribute to fecundity and hatching success; and total HUFAs are linked to improved hatch rates37. Dietary cholesterol has also been shown to stimulate ovarian development and vitellogenin production in L. vannamei23,67, as well as increase mating rates23. Similarly, supplementation with krill oil (rich in EPA and DHA) and phospholipids modulates genes involved in cholesterol metabolism, arachidonic acid pathways, and ovarian steroidogenesis16,48.In this context, the elevated vLDL and triglyceride levels observed in Group 4 females are consistent with active lipid mobilization for oogenesis. However, Group 3 achieved superior egg and larval production despite slightly lower vLDL levels, suggesting more efficient lipid utilization and vitellogenesis rather than simple lipid accumulation. Differences in lipid profiles between males and females further underscore sex-specific allocation strategies, with females prioritizing gonadal development and males allocating more resources toward somatic maintenance81,82.Reproduction imposes substantial metabolic demands that can lead to oxidative stress through the accumulation of reactive oxygen species (ROS)35,83. Antioxidant enzyme activities (SOD, GPx, and catalase) and lipid peroxidation levels (MDA) revealed diet-dependent patterns in the hepatopancreas. In females, Group 4 exhibited significantly elevated SOD activity, while GPx peaked in Groups 2 and 3. Notably, MDA—an indicator of lipid peroxidation—was highest in Group 2. Increased antioxidant activity in broodstock likely reflects heightened metabolic rates and the associated risk of oxidative damage during vitellogenesis84,85. Group 3, which yielded the highest larval output, maintained moderate GPx activity alongside lower MDA levels than Group 2, indicating reduced oxidative damage coupled with effective reproductive investment. In contrast, the elevated MDA in Group 2 suggests increased lipid peroxidation, which may compromise oocyte quality and help explain its lower reproductive performance despite certain favorable fatty-acid characteristics32. In males, reduced SOD and GPx activities in Groups 2 and 3 may reflect shifts in metabolic allocation toward gametogenesis86.Maintaining oxidative balance is therefore essential to protect lipid-rich yolk reserves from peroxidative damage. Diets that combine beneficial fatty acids with adequate antioxidant components—such as tocopherols and selenium—are likely to enhance broodstock performance and egg quality36,87. The diversified ingredients in Group 4 likely contributed a broader spectrum of exogenous antioxidants, including selenium and vitamin E from liver and mussels25,88, thereby supporting cellular protection and physiological resilience36,89. A modest increase in catalase activity observed in Group 4 males further supports the role of this diet in maintaining physiological homeostasis during peak reproductive activity, as moderate elevations in catalase activity may confer protective benefits against oxidative stress90.Histological examination of the hepatopancreas, ovaries, and vas deferens provided morphological support for the biochemical and reproductive responses. In Group 3, the hepatopancreas showed prominent epithelial vacuolation and variation in tubule and lumen morphology, suggesting altered metabolic activity during reproductive development and sustained vitellogenin synthesis80,91. Group 4 exhibited comparatively preserved hepatopancreatic architecture, with well-defined lumina and compact epithelial organization, indicating maintenance of hepatopancreatic function during reproductive maturation92.Group 2 displayed extensive epithelial vacuolation and irregular tubule morphology, consistent with its elevated MDA levels and lower reproductive performance. These alterations may impair the digestive and biosynthetic functions required for vitellogenesis and ovarian development93,94. Ovarian histology showed advanced oocyte development in Groups 3 and 4, with numerous vitellogenic oocytes. The sustained occurrence of Stage 4 females in Group 3 was consistent with its higher egg and larval production. This reproductive pattern may be associated with the nutritional value of feed ingredients such as chicken liver, which supplies bioavailable amino acids, polyunsaturated fatty acids, vitamins, and essential minerals that may support ovarian development and reproductive metabolism28.The hepatopancreas, as the central metabolic organ in shrimp, is particularly sensitive to nutritional status. Group 1 exhibited prominent epithelial vacuolation and variation in tubule and lumen morphology—may reflect nutrient depletion typical of late-season broodstock92,95. In contrast, Group 4 showed comparatively preserved tubular architecture, with well-defined star-shaped lumina and less extensive vacuolation, more closely resembling a healthy baseline condition.Finally, vas deferens histology revealed clear diet-dependent differences in epithelial activity and sperm density. Group 2 showed epithelial hyperplasia and intense apical secretion, potentially indicative of compensatory or stress-related responses96,97, In contrast, Group 4 exhibited a more balanced secretory profile, likely supporting efficient sperm maturation and transfer, and ultimately contributing to improved reproductive success97,98,99.When considered together with the observed variation in fatty-acid composition, the proximate-composition results demonstrate that the diets differed in both nutrient quantity and lipid quality, in addition to the number and combination of ingredients. Although Groups 1–3 contained comparable amounts of crude fat (8.22–8.25%), similar total lipid concentrations do not imply nutritional equivalence because diets with comparable lipid contents may differ substantially in their proportions of saturated, monounsaturated, and polyunsaturated fatty acids, including essential fatty acids such as EPA and DHA. Such differences can influence fatty-acid digestibility, tissue fatty-acid deposition, lipid metabolism, growth, and physiological condition100,101,102,103. The variability in fatty-acid composition observed among the present treatments therefore provides additional evidence that their biological effects cannot be interpreted solely in relation to dietary diversity.The treatments also differed in crude-protein concentration. Groups 1 and 2 both contained 80% worms, but substitution of squid with mussels was accompanied by a reduction in crude protein from 60.44% to 54.62% and an increase in ash from 9.97% to 10.36%. Consequently, any biological difference between these treatments may reflect differences in protein supply, amino-acid composition, digestibility, palatability, fatty-acid profile, or mineral composition rather than the replacement ingredient alone. Previous ingredient-evaluation studies have demonstrated that animal protein sources can differ in nutrient digestibility, utilization, and growth responses even when diets are formulated to contain comparable protein and energy concentrations104.Group 3, which contained three animal ingredients, had the second-highest crude-protein concentration (64.08%) and lower ash than Groups 1 and 2, whereas Group 4, containing four ingredients, had the highest crude protein (66.74%) and crude fat (9.29%) and the lowest moisture and ash contents. Dietary protein and lipid concentrations, their ratio, and lipid source can independently or interactively affect feed intake, growth, feed efficiency, body composition, fatty-acid deposition, blood biochemical indices, antioxidant status, and nitrogen excretion101,105,106,107.The strong reproductive performance observed in Groups 3 and 4 suggests that their dietary compositions supplied complementary nutrients that supported high-quality gamete production. In many invertebrates, lipids derived from mollusks, such as mussels, and cephalopods, such as squid, contribute to the deposition of DHA and other essential phospholipids in the egg yolk, thereby supporting embryonic development and larval viability108. Nevertheless, the responses differed between the two dietary treatments. Group 3 appeared to provide a more balanced reproductive response across sexes, with particularly favorable outcomes in males, whereas Group 4 produced comparatively stronger responses in females. These findings suggest that male and female broodstock may have distinct nutritional requirements for optimal reproductive performance.Because broodstock were assessed during the late phase of the reproductive season, further studies spanning the full reproductive cycle are needed to fully elucidate the long-term effects of these feed formulations on reproductive performance in L. vannamei.ConclusionThis study demonstrates that a heavy reliance on polychaetes alone is insufficient for sustaining reproductive performance toward the end of the breeding season in L. vannamei. In contrast, diversified feeding strategies—particularly the 50:25:25 worm/mussel/squid formulation (Group 3) and the four-component diet including chicken liver (Group 4)—substantially enhance reproductive outcomes. Group 3 was most effective in maximizing egg production and larval output, while Group 4 conferred pronounced molecular and physiological advantages, including improved gene expression profiles and antioxidant capacity. Overall, implementing diversified broodstock diets offers a practical and effective strategy for hatcheries to maintain high seed production and reproductive efficiency, even as broodstock approach the later stages of the production cycle.

    Data availability

    The majority of the data supporting the findings of this study are included within the manuscript. Other relevant datasets can be obtained from the corresponding author upon reasonable request.
    ReferencesFawzy, S., Mohammed, S., Abo-Al-Ela, H. G., Helal, M. A. & El-Kassas, S. Economic and productive impacts of dietary turnip (Brassica rapa) leaf powder on Nile tilapia (Oreochromis niloticus) fingerlings. Anim. Feed Sci. Technol. 333, 116622. https://doi.org/10.1016/j.anifeedsci.2025.116622 (2026).Article 
    CAS 

    Google Scholar 
    FAO. The State of World Fisheries and Aquaculture 2024 – Blue Transformation in action. Rome, Italy2024. https://doi.org/10.4060/cd0683enCastillo-Juárez, H., Campos-Montes, G. R., Caballero-Zamora, A. & Montaldo, H. H. Genetic improvement of Pacific white shrimp [Penaeus (Litopenaeus) vannamei]: perspectives for genomic selection. Front. Genet. 6, 93. https://doi.org/10.3389/fgene.2015.00093 (2015).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Abo-Al-Ela, H. G. The emerging regulatory roles of noncoding RNAs in immune function of fish: MicroRNAs versus long noncoding RNAs. Mol. Genet. Genomics 296(4), 765–781. https://doi.org/10.1007/s00438-021-01786-x (2021).Article 
    CAS 
    PubMed 

    Google Scholar 
    Ho, T. Q. & Do, H.-L. Production performance and risk in shrimp farming: The role of social trust belief. Aquaculture 612, 743167. https://doi.org/10.1016/j.aquaculture.2025.743167 (2026).Article 

    Google Scholar 
    Kittichotsatsawat, Y., Wattanutchariya, W., Jongjareonrak, A. & Seesuriyachan, P. Enhancing manufacturing operations within the supply chain for sustainable frozen shrimp production. Sustainability 17(6), 2412. https://doi.org/10.3390/su17062412 (2025).Article 
    ADS 

    Google Scholar 
    Said, M. M., Zaki, F. M., Taw, N. & Snyder, G. S. Light limitation alters water quality, biofloc composition, zooplankton community, and performance of the whiteleg shrimp (Litopenaeus vannamei) reared with intensive biofloc. Egypt. J. Aquat. Biol. Fish. 28 (3), 359–381. https://doi.org/10.21608/ejabf.2024.355970 (2024).Article 

    Google Scholar 
    Abo-Al-Ela, H. G., Mahdi, S., Angthong, P. & Rungrassamee, W. Probiotic modulation of key immune macromolecules in shrimp. Microb. Pathog. 203, 107463. https://doi.org/10.1016/j.micpath.2025.107463 (2025).Article 
    CAS 
    PubMed 

    Google Scholar 
    Zhou, L. et al. Shrimp shapes a nitrite tolerance trait via regulating autophagy and apoptosis. Int. J. Mol. Sci. 26(4), 1641. https://doi.org/10.3390/ijms26041641 (2025).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Emerenciano, M. G. C. et al. Intensification of penaeid shrimp culture: An applied review of advances in production systems, nutrition and breeding. Animals : an open access journal from MDPI 12(3), 236. (2022). https://doi.org/10.3390/ani12030236Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Galal, A. M., Said, M. M. & Sharaf, S. M. Reproductive studies on preparation of whiteleg shrimp Litopenaeus vannamei broodstock in commercial hatcheries. J. Anim. Poult. Fish. Prod. 11 (1), 35–42. https://doi.org/10.21608/japfp.2022.284091 (2022).Article 

    Google Scholar 
    Huang, Y. et al. Identification of structural variation related to spawn capability of Penaeus vannamei. BMC Genomics 25(1), 934. https://doi.org/10.1186/s12864-024-10863-5 (2024).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Barua, H. et al. Dietary phospholipids and their impact on crustacean physiology: Growth, metabolism, immunity, and beyond. Aquac Nutr. 2025 (1), 8180797. (2025). https://doi.org/10.1155/anu/8180797Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Dildar, T., Cui, W. & Ma, H. Physiology of ovarian development in crustaceans: Interactions among hormones, nutrients, and environmental factors from integrated perspectives. Aquac. Nutr. 1, 4900891. https://doi.org/10.1155/anu/4900891 (2025).Article 
    CAS 

    Google Scholar 
    Ren, Y. et al. Insight of vitellogenesis patterns: A comparative analysis of the differences between the primary and secondary vitellogenesis period in the ovary, hepatopancreas, and muscle of mud crab, Scylla paramamosain. Front. Genet. 13, 965070. https://doi.org/10.3389/fgene.2022.965070 (2022).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Liang, X. et al. Effect of dietary krill oil levels on the regulation of ovary development in Pacific white shrimp (Litopenaeus vannamei) broodstock. Aquaculture 589, 740968. https://doi.org/10.1016/j.aquaculture.2024.740968 (2024).Article 
    CAS 

    Google Scholar 
    Chen, T. et al. Pacific white shrimp (Litopenaeus vannamei) vitellogenesis-inhibiting hormone (VIH) is predominantly expressed in the brain and negatively regulates hepatopancreatic vitellogenin (VTG) gene expression. Biol. Reprod. 90 (3), 47. https://doi.org/10.1095/biolreprod.113.115030 (2014).Article 
    CAS 
    PubMed 

    Google Scholar 
    Yang, G. et al. Insight into the regulatory relationships between the insulin-like androgenic gland hormone gene and the insulin-like androgenic gland hormone-binding protein gene in giant freshwater prawns (Macrobrachium rosenbergii). Int. J. Mol. Sci. 21(12), 4207. https://doi.org/10.3390/ijms21124207 (2020).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Ortiz-Guillén, S., Magaña-Gallegos, E., Arévalo, M., Cuzon, G. & Gaxiola, G. Partial substitution of frozen-fresh food by an experimental diet in Litopenaeus vannamei (Boone, 1931) (Decapoda: Dendrobranchiata: Penaeidae) broodstock. Anim. Reprod. Sci. 242, 107001. https://doi.org/10.1016/j.anireprosci.2022.107001 (2022).Article 
    CAS 
    PubMed 

    Google Scholar 
    Sandeep, K. P., Sivaramakrishnan, T., Ambasankar, K. & Vijayan, K. K. Live Feeds in Shrimp Aquaculture: An Overview. In Aquatic Animal Health Management (eds Kooloth Valappil, R. et al.) 809–825 (Springer Nature Singapore, 2025). https://doi.org/10.1007/978-981-96-7987-4_28.Chapter 

    Google Scholar 
    Yang, D. et al. Gonadal maturation in Litopenaeus vannamei fed on four different polychaetes. Aquac. Rep. 22, 100920. https://doi.org/10.1016/j.aqrep.2021.100920 (2022).Article 

    Google Scholar 
    Li, S. et al. Effects of three polychaete species on growth and reproductive performance, biochemical indices, and histological condition of differing tissues in male broodstock of Pacific white shrimp, Litopenaeus vannamei. Aquaculture 598, 742055. https://doi.org/10.1016/j.aquaculture.2024.742055 (2025).Article 
    CAS 

    Google Scholar 
    Ma, Q., Wang, Z., Xu, H., Wei, Y. & Liang, M. Effects of dietary cholesterol on ovary development and reproductive capacity in Pacific white shrimp broodstock, Litopenaeus vannamei. Aquac. Rep. 38, 102346. https://doi.org/10.1016/j.aqrep.2024.102346 (2024).Article 

    Google Scholar 
    Thiruvasagam, T., Chidambaram, P., Ranjan, A. & Komuhi, N. B. Significance of fatty acids in fish broodstock nutrition. Anim. Reprod. Sci. 268, 107573. https://doi.org/10.1016/j.anireprosci.2024.107573 (2024).Article 
    CAS 
    PubMed 

    Google Scholar 
    Bjerknes, H., Elvevoll, E. O., Alterskjær Sundset, M., Langdal, A. & Eilertsen, K.-E. Farmed blue mussels (Mytilus edulis)—a nutrient-dense resource retaining nutritional value through processing. Front. Nutr. 11, 1443229. https://doi.org/10.3389/fnut.2024.1443229 (2024).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Aubourg, S. P., Trigo, M., Prego, R., Cobelo-García, A. & Medina, I. Nutritional and healthy value of chemical constituents obtained from patagonian squid (Doryteuthis gahi) by-products captured at different seasons. Foods 10(9), 2144. https://doi.org/10.3390/foods10092144 (2021).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Roo, J., Montero, D., Raquel, Q.-P., Monzón-Rivero, C. & Izquierdo López, M. Optimizing artemia enrichment: A low DHA/high EPA protocol for enhanced n-3 LC-HUFA levels to support greater amberjack (Seriola dumerili) larval rearing. Aquac. Nutr. 1, 5548991. https://doi.org/10.1155/2023/5548991 (2023).Article 

    Google Scholar 
    Majewska, D. et al. Comparison of the mineral and fatty acid profiles of ostrich, turkey and broiler chicken livers. Br. Poult. Sci. 57(2), 193–200. https://doi.org/10.1080/00071668.2016.1154136 (2016).Article 
    CAS 
    PubMed 

    Google Scholar 
    Liang, X. et al. Evalution of the synergistic effects of cholesterol and phospholipid on ovarian lipid accumulation, steroidogenesis and vitellogenesis in female broodstock in the Pacific whiteleg shrimp, Litopenaeus vannamei. Aquaculture 613, 743351. https://doi.org/10.1016/j.aquaculture.2025.743351 (2026).Article 
    CAS 

    Google Scholar 
    Wei, H. et al. Growth, nonspecific immunity, intestinal flora, hepatopancreas, and intestinal histological results for Litopenaeus vannamei fed with diets supplement with different animal by-products. Aquac. Rep. 29, 101521. https://doi.org/10.1016/j.aqrep.2023.101521 (2023).Article 

    Google Scholar 
    Lin, H. et al. Chicken liver meal as an effective fishmeal replacement for Litopenaeus vannamei diets: A comprehensive evaluation of growth, non-specific immunity and gut health. Aquac. Rep. 42, 102742. https://doi.org/10.1016/j.aqrep.2025.102742 (2025).Article 

    Google Scholar 
    Khan, M., Mushtaq, M., Usman, M., Rahman, M. A. U. & Quan, G. Oxidative stress-induced cytotoxicity and the role of dietary antioxidants in farm animals: A review. Advances in Redox Research 16, 100138. (2025). https://doi.org/10.1016/j.arres.2025.100138Article 
    CAS 

    Google Scholar 
    Abo-Al-Ela, H. G. et al. Genotoxic effects of 6PPD-quinone on juvenile European seabass (Dicentrarchus labrax): Impacts on DNA replication fidelity and oxidative DNA damage. Aquat. Toxicol. 296, 107830. https://doi.org/10.1016/j.aquatox.2026.107830 (2026).Article 
    CAS 
    PubMed 

    Google Scholar 
    Bedir, E. A., Said, M. M., Al Wakeel, R. A., Nasr-Allah, A. M. & Abo-Al-Ela, H. G. A comparative study of the genetically improved Abbassa Nile tilapia strain (GIANT-G9) and a commercial strain in Egypt: Growth vs. stress response. Vet. Res. Commun. 49(2), 75. https://doi.org/10.1007/s11259-024-10639-y (2025).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Almansa-Ordonez, A., Bellido, R., Vassena, R., Barragan, M. & Zambelli, F. Oxidative stress in reproduction: A mitochondrial perspective. Biology 9(9), 269. https://doi.org/10.3390/biology9090269 (2020).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    El–Dahhar, A. A., Khalifa, G. A., El–Zaeem, S. Y., Abdel-Rahim, M. M. & Mourad, M. M. Impact of organic and nano selenium on sexual maturation, egg, and offspring quality of European sea bass (Dicentrarchus labrax) broodstock. Sci. Rep. 16 (1), 1342. https://doi.org/10.1038/s41598-025-29108-w (2026).Article 
    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Racotta, I. S., Palacios, E. & Ibarra, A. M. Shrimp larval quality in relation to broodstock condition. Aquaculture 227 (1), 107–130. https://doi.org/10.1016/S0044-8486(03)00498-8 (2003).Article 

    Google Scholar 
    Chen, K., Li, S., Xiang, J., Sagi, A. & Li, F. Transcriptome analysis reveals the endocrine regulation on the expression of IAG in Litopenaeus vannamei. J. Mar. Sci. Eng. 9(6), 677. https://doi.org/10.3390/jmse9060677 (2021).Article 

    Google Scholar 
    Levy, T. & Sagi, A. The IAG-switch—a key controlling element in decapod crustacean sex differentiation. Front. Endocrinol. 11, 651. https://doi.org/10.3389/fendo.2020.00651 (2020).Article 

    Google Scholar 
    Chen, Z. et al. Compared to fishmeal, dietary soybean meal improves the reproductive performance of female yellow catfish (Pelteobagrus fulvidraco) broodstock. Aquac. Nutr. 1, 6240803. https://doi.org/10.1155/2023/6240803 (2023).Article 
    CAS 

    Google Scholar 
    Wouters, R., Lavens, P., Nieto, J. & Sorgeloos, P. Penaeid shrimp broodstock nutrition: An updated review on research and development. Aquaculture 202(1), 1–21. https://doi.org/10.1016/S0044-8486(01)00570-1 (2001).Article 

    Google Scholar 
    Li, S. et al. Effect of three polychaetes on growth and reproductive performance, biochemical indices and histology of different tissues in the female Pacific white shrimp, Litopenaeus vannamei broodstock. Anim. Nutr. 21, 49–62. https://doi.org/10.1016/j.aninu.2024.11.020 (2025).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Estante-Superio, E. G. et al. P. Inclusion of live mud polychaete (Marphysa iloiloensis) in the feeding regime improved the hatchery performance of domesticated Indian white shrimp (Penaeus indicus). Reg. Stud. Mar. Sci. 62, 102923. https://doi.org/10.1016/j.rsma.2023.102923 (2023).Article 

    Google Scholar 
    Sittikankaew, K. et al. Transcriptome analyses reveal the synergistic effects of feeding and eyestalk ablation on ovarian maturation in black tiger shrimp. Sci. Rep. 10(1), 3239. https://doi.org/10.1038/s41598-020-60192-2 (2020).Article 
    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    AOAC. Official methods of analysis of AOAC international 18th Ed. (Association of Official Analytical Chemists, 2007).
    Google Scholar 
    Danish, M. & Nizami, M. Complete fatty acid analysis data of flaxseed oil using GC-FID method. Data Brief 23, 103845. https://doi.org/10.1016/j.dib.2019.103845 (2019).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Dall, W., Hill, B. J., Rothlisberg, P. C. & Sharples, D. J. in The Biology of the Penaeidae. Advances in Marine Biology. (eds Blaxter, J. H. S. & Southward, A. J.) (Academic Press, 1991).Liang, X. et al. Effects and mechanism of different phospholipid diets on ovary development in female broodstock pacific white shrimp, Litopenaeus vannamei. Front. Nutr. 9, 830934. https://doi.org/10.3389/fnut.2022.830934 (2022).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Tito, L., Vera, M., Ordinola, E., Alemán, S. & Perea, Á. Histologically validated gonadal maturity scale of Litopenaeus vannamei (Boone, 1931) females, Peru. Bol. Instituto del. Mar. del. Perú. 39 (1), 26–37. https://doi.org/10.53554/boletin.v39i1.396 (2024).Article 

    Google Scholar 
    Friedewald, W. T., Levy, R. I. & Fredrickson, D. S. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin. Chem. 18(6), 499–502. https://doi.org/10.1093/clinchem/18.6.499 (1972).Article 
    CAS 
    PubMed 

    Google Scholar 
    Pang, Y. et al. Characterization and expression analysis of insulin growth factor binding proteins (IGFBPs) in Pacific white shrimp Litopenaeus vannamei. Int. J. Mol. Sci. 22(3), 1056. https://doi.org/10.3390/ijms22031056 (2021).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Zhang, X. et al. Crustacean hyperglycemic hormone (CHH) regulates the ammonia excretion and metabolism in white shrimp, Litopenaeus vannamei under ammonia-N stress. Sci. Total Environ. 723, 138128. https://doi.org/10.1016/j.scitotenv.2020.138128 (2020).Article 
    CAS 
    PubMed 

    Google Scholar 
    Bell, T. A. & Lightner, D. V. A handbook of normal penaeid shrimp histology (World Aquaculture Society, 1988).
    Google Scholar 
    Campos-Ramos, R., Garza-Torres, R., Guerrero-Tortolero, D. A., Maeda-Martínez, A. M. & Obregón-Barboza, H. Environmental sex determination, external sex differentiation and structure of the androgenic gland in the Pacific white shrimp Litopenaeus vannamei (Boone). Aquac. Res. 37 (15), 1583–1593. https://doi.org/10.1111/j.1365-2109.2006.01604.x (2006).Article 

    Google Scholar 
    de Rooij, M., Breemer, L., Woestenburg, D. & Busing, F. Logistic multidimensional data analysis for ordinal response variables using a cumulative link function. Psychometrika 90(3), 833–869. https://doi.org/10.1017/psy.2025.10 (2025).Article 
    MathSciNet 

    Google Scholar 
    Hosmer, D. W. J., Lemeshow, S. & Sturdivant, R. X. Applied Logistic Regression 3rd ed. (John Wiley & Sons, Inc., 2013). https://doi.org/10.1002/9781118548387.Book 

    Google Scholar 
    Muhd-Farouk, H. et al. Development of ovarian maturations in orange mud crab, Scylla olivacea (Herbst, 1796) through induction of eyestalk ablation and methyl farnesoate. Arab. J. Basic. Appl. Sci. 26 (1), 171–181. https://doi.org/10.1080/25765299.2019.1588197 (2019).Article 

    Google Scholar 
    Foglietta, L. M., Camejo, M., Gallardo, L. & Herrera, F. C. A maturity index for holothurians exhibiting asynchronous development of gonad tubules. J. Exp. Mar. Biol. Ecol. 303(1), 19–30. https://doi.org/10.1016/j.jembe.2003.10.019 (2004).Article 

    Google Scholar 
    Lu, Y. et al. Vitellogenin and its receptor characterize exogenous vitellogenesis mode in the Pacific white shrimp Litopenaeus vannamei. Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology 311, 111947. (2026). https://doi.org/10.1016/j.cbpa.2025.111947Article 
    CAS 
    PubMed 

    Google Scholar 
    Quesada-Ávila, I. et al. Serotonin plus spiperone induce ovarian maturation, spawning and potential release of ovarian maturation inducing pheromones in the Pacific white shrimp Litopenaeus vannamei (Penaeidae). Aquac. Rep. 37, 102257. https://doi.org/10.1016/j.aqrep.2024.102257 (2024).Article 

    Google Scholar 
    Seabold, S., Perktold, J. & Statsmodels. Econometric and statistical modeling with python. Proceedings of the 9th Python in Science Conference, SciPy 2010; Austin, Texas. (2010). https://doi.org/10.25080/Majora-92bf1922-011Song, J. et al. The dietary lipid requirement for ovarian maturation and health in female giant river prawn, Macrobrachium rosenbergii broodstock. Aquac. Nutr. 1, 7462841. https://doi.org/10.1155/2024/7462841 (2024).Article 
    CAS 

    Google Scholar 
    Huang, W. et al. Effects of dietary n-3HUFA levels on growth, antioxidant capacity, and ovarian development in Litopenaeus vannamei broodstock shrimp. Aquac. Rep. 46, 103317. https://doi.org/10.1016/j.aqrep.2025.103317 (2026).Article 

    Google Scholar 
    Wang, H. et al. Growth and nutritional composition of the polychaete Hediste diversicolor (OF Müller, 1776) cultivated on waste from land-based salmon smolt aquaculture. Aquaculture 502, 232–241. https://doi.org/10.1016/j.aquaculture.2018.12.047 (2019).Article 

    Google Scholar 
    Jerónimo, D. et al. Unravelling the fatty acid profiles of different polychaete species cultured under integrated multi-trophic aquaculture (IMTA). Sci. Rep. 11(1), 10812. https://doi.org/10.1038/s41598-021-90185-8 (2021).Article 
    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Liang, X. et al. Growth, health, and gut microbiota of female Pacific white shrimp, Litopenaeus vannamei broodstock fed different phospholipid sources. Antioxidants 11(6), 1143. https://doi.org/10.3390/antiox11061143 (2022).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Liang, X. et al. Positive effects of optimal dietary cholesterol levels on the ovary development and health of female Pacific white shrimp, Litopenaeus vannamei broodstock. Aquaculture 577, 739987. https://doi.org/10.1016/j.aquaculture.2023.739987 (2023).Article 
    CAS 

    Google Scholar 
    Jones, P. J. H. & Lichtenstein, A. H. Chapter 4 – Lipids. In Present Knowledge in Nutrition (eds Marriott, B. P. et al.) 51–69 (Academic, 2020). . https://doi.org/10.1016/B978-0-323-66162-1.00004-4Chapter 

    Google Scholar 
    Zhou, L. & Nilsson, Å. Sources of eicosanoid precursor fatty acid pools in tissues. J. Lipid Res. 42 (10), 1521–1542. https://doi.org/10.1016/S0022-2275(20)32206-9 (2001).Article 
    CAS 
    PubMed 

    Google Scholar 
    Mori, T. A. & Hodgson, J. M. Fatty acids: Health Effects of Omega-6 Polyunsaturated Fatty Acids. In: Caballero B, editor. Encyclopedia of Human Nutrition (Third Edition). Waltham: Academic Press; pp. 209–214. (2013). https://doi.org/10.1016/B978-0-12-375083-9.00100-8Yepiz-Plascencia, G., Vargas-Albores, F. & Higuera-Ciapara, I. Penaeid shrimp hemolymph lipoproteins. Aquaculture 191 (1), 177–189. https://doi.org/10.1016/S0044-8486(00)00427-0 (2000).Article 
    CAS 

    Google Scholar 
    Zhang, X., Yin, Y., Fan, H., Zhou, Q. & Jiao, L. Arginine promoted ovarian development in Pacific white shrimp Litopenaeus vannamei via the NO-sGC-cGMP and TORC1 signaling pathways. Animals 14(13), 1986. https://doi.org/10.3390/ani14131986 (2024).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Liu, Y. et al. Changes revealed by histology, physiology, enteric microbiota, and transcriptome after ovarian ripening of shrimp, Litopenaeus vannamei, through using fresh baits support the innovation of artificial ripening of shrimps. Aquac. Rep. 39, 102494. https://doi.org/10.1016/j.aqrep.2024.102494 (2024).Article 

    Google Scholar 
    Pan, J., Liu, M., Chen, T., Cheng, Y. & Wu, X. Immunolocalization and changes of 17beta-estradiol during ovarian development of Chinese mitten crab Eriocheir Sinensis. Cell Tissue Res. 373 (2), 509–520. https://doi.org/10.1007/s00441-018-2834-x (2018).Article 
    CAS 
    PubMed 

    Google Scholar 
    Zhao, Z., Peng, L., Zhao, Q. & Wang, Z. Genome-wide identification, expression and response to estrogen of vitellogenin gene family in Sichuan bream (Sinibrama taeniatus). Int. J. Mol. Sci. 25(12), 6739. https://doi.org/10.3390/ijms25126739 (2024).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Feng, H. et al. Involvement of insulin-like growth factor binding proteins (IGFBPs) and activation of insulin/IGF-like signaling (IIS)-target of rapamycin (TOR) signaling cascade in pacific white shrimp, Litopenaeus vannamei exposed to acute low-salinity. Aquac. Rep. 30, 101627. https://doi.org/10.1016/j.aqrep.2023.101627 (2023).Article 

    Google Scholar 
    Fu, C. et al. Comparative transcriptome analysis reveals related regulatory mechanisms of androgenic gland in Eriocheir sinensis. Biomed Res. Int. 1, 4956216. https://doi.org/10.1155/2017/4956216 (2017).Article 
    CAS 

    Google Scholar 
    Tan, K. et al. Transcriptomic profiling and novel insights into the effect of AG ablation on gonad development in Macrobrachium rosenbergii. Aquaculture 556, 738224. https://doi.org/10.1016/j.aquaculture.2022.738224 (2022).Article 
    CAS 

    Google Scholar 
    Kumar, V. et al. Metabolism and nutritive role of cholesterol in the growth, gonadal development, and reproduction of crustaceans. Rev. Fish. Sci. Aquac. 26(2), 254–273. https://doi.org/10.1080/23308249.2018.1429384 (2018).Article 

    Google Scholar 
    Li, Z. et al. Transcriptomic analysis reveals yolk accumulation mechanism from the hepatopancreas to ovary in the pacific white shrimp Litopenaeus vannamei. Front. Mar. Sci. 9, 948105. https://doi.org/10.3389/fmars.2022.948105 (2022).Article 

    Google Scholar 
    Martínez-Pita, I., García, F. J. & Pita, M. L. Males and females gonad fatty acids of the sea urchins Paracentrotus lividus and Arbacia lixula (Echinodermata). Helgol. Mar. Res. 64 (2), 135–142. https://doi.org/10.1007/s10152-009-0174-7 (2010).Article 
    ADS 

    Google Scholar 
    Martin-Creuzburg, D., Massier, T. & Wacker, A. Sex-specific differences in essential lipid requirements of Daphnia magna. Front. Ecol. Evol. 6, 89. https://doi.org/10.3389/fevo.2018.00089 (2018).Article 

    Google Scholar 
    Wu, L. T. & Chu, K. H. Characterization of an ovary-specific glutathione peroxidase from the shrimp Metapenaeus ensis and its role in crustacean reproduction. Comparative Biochemistry and Physiology Part B: Biochemistry and Molecular Biology 155(1), 26–33. (2010). https://doi.org/10.1016/j.cbpb.2009.09.005Article 
    CAS 
    PubMed 

    Google Scholar 
    Zanitti, M. et al. Diclofenac activates antioxidant defenses and affects ovarian growth in the freshwater shrimp Macrobrachium borellii. Environ. Toxicol. Chem. 44(6), 1647–1654. https://doi.org/10.1093/etojnl/vgaf070 (2025).Article 
    CAS 
    PubMed 

    Google Scholar 
    Montalvo, G. et al. Effect of vitamin C on the antioxidant and immune response of male white shrimp (Penaeus vannamei) broodstock. Antioxidants 14(8), 988. https://doi.org/10.3390/antiox14080988 (2025).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Ladisa, C., Ma, Y. & Habibi, H. R. Metabolic changes during growth and reproductive phases in the liver of female goldfish (Carassius auratus). Front. Cell. Dev. Biol. 10, 834688. https://doi.org/10.3389/fcell.2022.834688 (2022).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Montalvo, G. et al. Immune gene expression and antioxidant response to vitamin E enriched diets for males Litopenaeus vannamei breeder (Boone, 1931). Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology 268, 111187. (2022). https://doi.org/10.1016/j.cbpa.2022.111187Article 
    CAS 
    PubMed 

    Google Scholar 
    Pilarczyk, B. et al. Selenium content in selected products of animal origin and estimation of the degree of cover daily Se requirement in Poland. Int. J. Food Sci. Technol. 45(1), 186–191. https://doi.org/10.1111/j.1365-2621.2009.02120.x (2010).Article 
    CAS 

    Google Scholar 
    AbdEl-Kader, M. F. et al. Dietary selenium nanoparticles positively modulate the growth and immunity of seabream (Sparus aurata) fingerlings exposed to low salinity stress and Vibrio parahaemolyticus challenge. Aquaculture 576, 739893. https://doi.org/10.1016/j.aquaculture.2023.739893 (2023).Article 
    CAS 

    Google Scholar 
    Abo-Al-Ela, H. G., Siddique, M., Sharawy, Z. Z. & Tahoun, A.-A. Growth and immune effects of standardized Ashwagandha extract in whiteleg shrimp, Litopenaeus vannamei, post-larvae. Aquaculture 617, 743758. https://doi.org/10.1016/j.aquaculture.2026.743758 (2026).Article 
    CAS 

    Google Scholar 
    Du, L. et al. Single-cell RNA sequencing reveals the heterogeneity of hepatopancreas cells and their association with gonadal development in the red swamp crayfish Procambarus clarkii. Aquaculture 601, 742311. https://doi.org/10.1016/j.aquaculture.2025.742311 (2025).Article 
    CAS 

    Google Scholar 
    Yuanlae, S. et al. Ecytonucleospora hepatopenaei causes lipid droplet depletion and imbalanced lipid metabolism in Penaeus vannamei. Sci. Rep. 15(1), 37198. https://doi.org/10.1038/s41598-025-21037-y (2025).Article 
    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Arcos, F. G., Ibarra, A. M. & Racotta, L. S. Vitellogenin in hemolymph predicts gonad maturity in adult female Litopenaeus (Penaeus) vannamei shrimp. Aquaculture 316(1), 93–98. https://doi.org/10.1016/j.aquaculture.2011.02.045 (2011).Article 
    CAS 

    Google Scholar 
    Sun, S. et al. Transciptomic and histological analysis of hepatopancreas, muscle and gill tissues of oriental river prawn (Macrobrachium nipponense) in response to chronic hypoxia. BMC Genom. 16 (1), 491. https://doi.org/10.1186/s12864-015-1701-3 (2015).Article 
    CAS 

    Google Scholar 
    Vogt, G. Functional cytology of the hepatopancreas of decapod crustaceans. J. Morphol. 280 (9), 1405–1444. https://doi.org/10.1002/jmor.21040 (2019).Article 
    CAS 
    PubMed 

    Google Scholar 
    Bugnot, A. B. & López Greco, L. S. Structural alterations in the male reproductive system of the freshwater crayfish, Cherax quadricarinatus (Decapoda, Parastacidae). J. Invertebr. Pathol. 102(2), 160–166. https://doi.org/10.1016/j.jip.2009.08.004 (2009).Article 
    PubMed 

    Google Scholar 
    Fransozo, V., Fernandes, A. B., López-Greco, L. S., Zara, F. J. & Santos, D. C. Functional morphology of the male reproductive system of the white shrimp Litopenaeus schmitti (Burkenroad, 1936) (Crustacea, Penaeidea) compared to other Litopenaeus. Invertebrate Reprod. Dev. 60 (3), 161–174. https://doi.org/10.1080/07924259.2016.1174158 (2016).Article 
    CAS 

    Google Scholar 
    Chow, S., Dougherty, M. M., Dougherty, W. J. & Sandifer, P. A. Spermatophore formation in the white shrimps Penaeus Setiferus and P. Vannamei. J. Crustac. Biol. 11 (2), 201–216. https://doi.org/10.1163/193724091X00022 (1991).Article 

    Google Scholar 
    Feng, T., Paterson, B. & Johnston, S. A morphological study of the male reproductive tract, post-testicular acrosome maturation and spermatophore formation in the black tiger prawn (Penaeus monodon). J. Morphol. 279(9), 1290–1300. https://doi.org/10.1002/jmor.20869 (2018).Article 
    PubMed 

    Google Scholar 
    Qiu, H. et al. Dietary lipid sources influence fatty acid composition in tissue of large yellow croaker (Larmichthys crocea) by regulating triacylglycerol synthesis and catabolism at the transcriptional level. PLoS One 12(1), e0169985. https://doi.org/10.1371/journal.pone.0169985 (2017).Article 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Karalazos, V., Bendiksen, E. Å., Dick, J. R., Tocher, D. R. & Gordon Bell, J. Influence of the dietary protein:lipid ratio and fish oil substitution on fatty acid composition and metabolism of Atlantic salmon (Salmo salar) reared at high water temperatures. Br. J. Nutr. 105(7), 1012–1025. https://doi.org/10.1017/S0007114510004605 (2011).Article 
    CAS 
    PubMed 

    Google Scholar 
    Qin, Y. et al. Growth performance, fatty acid composition, and lipid metabolism are altered in groupers (Epinephelus coioides) by dietary fish oil replacement with palm oil. Anim. Nutr. 8, 102–113. https://doi.org/10.1016/j.aninu.2021.04.007 (2022).Article 
    CAS 
    PubMed 

    Google Scholar 
    Bell, J. G., Tocher, D. R., Henderson, R. J., Dick, J. R. & Crampton, V. O. Altered fatty acid compositions in Atlantic salmon (Salmo salar) fed diets containing linseed and rapeseed oils can be partially restored by a subsequent fish oil finishing diet. J. Nutr. 133 (9), 2793–2801. https://doi.org/10.1093/jn/133.9.2793 (2003).Article 
    CAS 
    PubMed 

    Google Scholar 
    Glencross, B. D., Booth, M. & Allan, G. L. A feed is only as good as its ingredients – a review of ingredient evaluation strategies for aquaculture feeds. Aquac Nutr. 13 (1), 17–34. (2007). https://doi.org/10.1111/j.1365-2095.2007.00450.xArticle 
    CAS 

    Google Scholar 
    Sankian, Z., Khosravi, S., Kim, Y.-O. & Lee, S.-M. Effect of dietary protein and lipid level on growth, feed utilization, and muscle composition in golden mandarin fish Siniperca scherzeri. Fish. Aquat. Sci. 20(1), 7. https://doi.org/10.1186/s41240-017-0053-0 (2017).Article 
    CAS 

    Google Scholar 
    Rahimnejad, S. et al. Effects of dietary protein and lipid levels on growth, body composition, blood biochemistry, antioxidant capacity and ammonia excretion of European grayling (Thymallus thymallus). Front. Mar. Sci. 8, 715636. https://doi.org/10.3389/fmars.2021.715636 (2021).Article 

    Google Scholar 
    Mohanta, K. N., Subramanian, S. & Korikanthimath, V. S. Effect of dietary protein and lipid levels on growth, nutrient utilization and whole-body composition of blue gourami, Trichogaster trichopterus fingerlings. J. Anim. Physiol. Anim. Nutr. 97 (1), 126–136. https://doi.org/10.1111/j.1439-0396.2011.01258.x (2013).Article 
    CAS 

    Google Scholar 
    Zuo, R. et al. Effects of dietary lipid sources on the growth, gonad development, fatty acid composition and spawning performance of broodstock, and early larvae quality of sea urchin (Strongylocentrotus intermedius). Front. Mar. Sci. 9, 927116. https://doi.org/10.3389/fmars.2022.927116 (2022).Article 

    Google Scholar 
    Download referencesAuthor informationAuthors and AffiliationsGenetics and Biotechnology, Department of Aquaculture, Faculty of Fish Resources, Suez University, Suez, 43221, EgyptHaitham G. Abo-Al-ElaShrimp Hatchery, National Company for Fisheries and Aquaculture, Delta Aquaculture Sector, Kafr El-Sheikh, Ghalioun, EgyptAmira MosbahDepartment of Aquaculture, Faculty of Fish Resources, Suez University, Suez, 43221, EgyptAl-Azab TahounAuthorsHaitham G. Abo-Al-ElaView author publicationsSearch author on:PubMed Google ScholarAmira MosbahView author publicationsSearch author on:PubMed Google ScholarAl-Azab TahounView author publicationsSearch author on:PubMed Google ScholarContributionsConceptualization: Abo-Al-Ela H.G. and Tahoun A. Data Curation: Abo-Al-Ela H.G. and Mosbah A. Formal analysis: Abo-Al-Ela H.G. and Mosbah A. Investigation: Abo-Al-Ela H.G., Mosbah A., and Tahoun A. Methodology: Abo-Al-Ela H.G., Mosbah A. Writing – original draft: Abo-Al-Ela H.G. Writing – review & editing: Abo-Al-Ela H.G. All authors discussed the results and commented on the manuscript.Corresponding authorCorrespondence to
    Haitham G. Abo-Al-Ela.Ethics declarations

    Competing interests
    The authors declare no competing interests.

    Declaration of Generative AI and AI-Assisted Technologies in the Writing Process
    The authors used ChatGPT to improve the language and readability of this manuscript. Following its use, all authors reviewed and edited the text as necessary and accept full responsibility for the final content of the manuscript.

    Ethics approval
    The experiment was conducted in accordance with both international and national standards for animal research, having received approval from the Faculty of Fish Resources at Suez University, Egypt (Approval No. SUEZ Sci_IRB: 31/12/2025/36). The study adheres to the ARRIVE guidelines (https://arriveguidelines.org/).

    Additional informationPublisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary InformationBelow is the link to the electronic supplementary material.Supplementary Material 1 (download DOCX )Supplementary Material 2 (download XLSX )Rights and permissions
    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
    Reprints and permissionsAbout this articleCite this articleAbo-Al-Ela, H.G., Mosbah, A. & Tahoun, AA. Different dietary compositions shape reproductive physiology in Litopenaeus vannamei broodstock.
    Sci Rep 16, 31483 (2026). https://doi.org/10.1038/s41598-026-71783-wDownload citationReceived: 02 July 2026Accepted: 11 September 2026Published: 08 October 2026Version of record: 08 October 2026DOI: https://doi.org/10.1038/s41598-026-71783-wShare this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy shareable link to clipboard
    Provided by the Springer Nature SharedIt content-sharing initiative
    KeywordsBroodstock nutritionFatty acid compositionGonadal development
    Litopenaeus vannamei
    Reproductive performanceVitellogenesis More

  • in

    Spatial variation shapes gut microbiome composition and functional potential in an Afrotropical fish

    AbstractHost-associated microbiomes are shaped by environmental heterogeneity across spatial gradients. However, the extent to which site-level variation structures gut microbial communities in Afrotropical fishes remains poorly resolved. In this study, we employed high-throughput 16S rRNA gene sequencing to characterize the gut microbiome of the African forktail catfish (Chrysichthys nigrodigitatus) encompassing two contrasting habitat types, with the aim of assessing spatial patterns in microbial composition and predicted functional potential. Across all sites, gut microbial communities were dominated by Firmicutes, Proteobacteria, Bacteroidota, and Fusobacteriota, with a conserved core microbiome shared among populations. Alpha diversity metrics (Observed richness, Shannon, Simpson, and Chao1 indices) did not differ significantly among sites, indicating comparable levels of microbial richness and evenness. In contrast, beta diversity analysis based on Bray–Curtis dissimilarity revealed pronounced site-specific structuring and clear differentiation between habitat types. Differential abundance and indicator taxa analyses identified distinct bacterial assemblages associated with specific habitats, reflecting localized environmental influences. Functional predictions indicated enrichment of pathways related to xenobiotic degradation, oxidative stress response, and aromatic compound metabolism in lagoon habitat. The concordance between taxonomic composition and predicted functional profiles underscores the role of site variation in shaping gut microbiome structure. Collectively, these findings demonstrate that spatial variation is a key determinant of gut microbiome composition and predicted functional capacity in this widely distributed Afrotropical fish species, providing new insights into host-microbiome-environment interactions in tropical aquatic systems.

    Explore related subjects
    Discover the latest articles and news in related subjects.

    Ecology

    Microbiology

    AcknowledgementsWe thank Olaide Bello for assisting with fieldwork and the local fishers across Lagos Lagoon for their support during field activities.FundingThe work was supported by the Faculty Start-Up grant from the Office of the Provost, Howard University, Washington, DC, to LMN.Author informationAuthors and AffiliationsDepartment of Biology, Howard University, Washington, DC, USAOlabisi Tawakalt Atofarati, Segun Olayinka Oladipo, Maame Kobe Asiamah, Amaya Everett & Lotanna Micah NnejiDepartment of Marine Biology, Shantou University, Shantou, ChinaIfeanyi Christopher NnejiFederal College of Wildlife Management, New Bussa, Niger State, NigeriaKehinde Moruf AdelakunAuthorsOlabisi Tawakalt AtofaratiView author publicationsSearch author on:PubMed Google ScholarSegun Olayinka OladipoView author publicationsSearch author on:PubMed Google ScholarIfeanyi Christopher NnejiView author publicationsSearch author on:PubMed Google ScholarKehinde Moruf AdelakunView author publicationsSearch author on:PubMed Google ScholarMaame Kobe AsiamahView author publicationsSearch author on:PubMed Google ScholarAmaya EverettView author publicationsSearch author on:PubMed Google ScholarLotanna Micah NnejiView author publicationsSearch author on:PubMed Google ScholarCorresponding authorCorrespondence to
    Lotanna Micah Nneji.Ethics declarations

    Competing interests
    The authors declare no competing interests.

    Additional informationPublisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary InformationSupplementary Information 1. (download PDF )Supplementary Information 2. (download PDF )Supplementary Information 3. (download DOCX )Rights and permissions
    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
    Reprints and permissionsAbout this articleCite this articleAtofarati, O.T., Oladipo, S.O., Nneji, I.C. et al. Spatial variation shapes gut microbiome composition and functional potential in an Afrotropical fish.
    Sci Rep (2026). https://doi.org/10.1038/s41598-026-71562-7Download citationReceived: 10 May 2026Accepted: 10 September 2026Published: 08 October 2026DOI: https://doi.org/10.1038/s41598-026-71562-7Share this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy shareable link to clipboard
    Provided by the Springer Nature SharedIt content-sharing initiative
    Keywords
    Chrysichthys nigrodigitatus
    16S rRNA geneGut microbiomeFunctional potential More

  • in

    Multilocus phylogeny and historical biogeography of the pill millipede genus Rhopalomeris (Diplopoda, Glomerida) in Thailand

    AbstractSoutheast Asia is widely recognized as a global biodiversity hotspot with a complex geological history and high species richness across multiple taxa. The pill millipede genus Rhopalomeris Verhoeff, 1906 has received increasing taxonomic attention in recent years, with several new species described, particularly from western Southeast Asia, including Thailand and Myanmar. However, its historical biogeography has remained insufficiently understood. In this study, we used an integrative analysis combining morphology, comprehensive molecular phylogeny based on mitochondrial and nuclear genes, and biogeographic approaches to provide new insights into the evolutionary history of Rhopalomeris in Thailand. The multilocus phylogeny indicated the presence of 21 species-level subclades, with nine previously described species and 12 putative species. Divergence time estimation and historical biogeography analysis revealed that Rhopalomeris likely originated in the Tenasserim Range and the inner area of Indochina around the Late Eocene to Early Oligocene (41.39–28.89 Mya, mean age = 34.78 Mya), followed by southward expansion into the Thai-Malay Peninsula and a series of subsequent diversifications shaped by recurrent sea-level fluctuations and multidirectional dispersal. Our findings highlight the dynamic interplay between Earth’s history and biogeographic processes driving the diversity of these Southeast Asian pill millipedes.

    Explore related subjects
    Discover the latest articles and news in related subjects.

    Ecology

    Evolution

    Genetics

    Zoology

    AcknowledgementsThe authors would like to thank members of the (1) Animal Systematics and Molecular Ecology laboratory, Mahidol University, (2) Animal Systematics Research Unit, Chulalongkorn University, (3) Animal Systematics Research, Khon Kaen University, and (4) Leibniz Institute for the Study of Biodiversity Change, Bonn, Germany, for their guidance and assistance both in the field and laboratory works, and D.J. Anderson for grammar checking and comments.FundingThis research project is supported by National Research Council of Thailand (NRCT) (Contract No. N41A661155) to P.S. and E.J., and received funding support from the NSRF via the Research and Innovation Acceleration Agency for Competitiveness and Area Development (RCAD) (Program Management Unit for Frontier Brainpower and Future Industries) (Contract No. B42G670038) to E.J. and N.L.Author informationAuthors and AffiliationsAnimal Systematics and Molecular Ecology Laboratory, Department of Biology, Faculty of Science, Mahidol University, Bangkok, 10400, ThailandPichsinee Sapparojpattana, Warut Siriwut & Ekgachai JeratthitikulProgram of Agriculture, Faculty of Agricultural Production, Maejo University, Chiang Mai, 50290, ThailandNatdanai LikhitrakarnDepartment of Biology, Faculty of Science, Khon Kaen University, Khon Kaen, 40002, ThailandRuttapon SrisonchaiLeibniz Institute for the Study of Biodiversity Change, Bonn, GermanyThomas WesenerAuthorsPichsinee SapparojpattanaView author publicationsSearch author on:PubMed Google ScholarNatdanai LikhitrakarnView author publicationsSearch author on:PubMed Google ScholarWarut SiriwutView author publicationsSearch author on:PubMed Google ScholarRuttapon SrisonchaiView author publicationsSearch author on:PubMed Google ScholarThomas WesenerView author publicationsSearch author on:PubMed Google ScholarEkgachai JeratthitikulView author publicationsSearch author on:PubMed Google ScholarCorresponding authorCorrespondence to
    Ekgachai Jeratthitikul.Ethics declarations

    Competing interests
    The authors declare no competing interests.

    Additional informationPublisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary InformationBelow is the link to the electronic supplementary material.Supplementary Material 1 (download DOCX )Rights and permissions
    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
    Reprints and permissionsAbout this articleCite this articleSapparojpattana, P., Likhitrakarn, N., Siriwut, W. et al. Multilocus phylogeny and historical biogeography of the pill millipede genus Rhopalomeris (Diplopoda, Glomerida) in Thailand.
    Sci Rep (2026). https://doi.org/10.1038/s41598-026-73309-wDownload citationReceived: 15 July 2026Accepted: 22 September 2026Published: 08 October 2026DOI: https://doi.org/10.1038/s41598-026-73309-wShare this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy shareable link to clipboard
    Provided by the Springer Nature SharedIt content-sharing initiative
    KeywordsBiodiversityEvolutionIndochinaMyriapodaPhylogeographySystematics More

  • in

    Contrasting plant trait strategies underpin ecosystem stability in temperate and alpine grasslands

    AbstractPlant functional traits are key drivers of ecosystem stability, yet the extent to which trait-stability relationships vary among different ecosystems remains poorly understood. Here, we provide a comparative analysis of trait-stability relationships across temperate and alpine grasslands. By integrating in situ measurements of plant functional traits to satellite-derived indices of ecosystem temporal stability, resistance, and resilience for the two types of grasslands in China, we find strikingly biome-dependent patterns. In temperate grasslands, where fast-growing, resource-acquisitive annuals are more prevalent, acquisitive strategies enhance ecosystem temporal stability and resistance. In contrast, in alpine grasslands where annuals are less common, these same strategies reduce ecosystem temporal stability and resistance, and perennial species with slow-growing, conservative traits play the stabilizing role. Ecosystem resilience, however, shows consistent positive responses to acquisitive traits in both grassland types. These patterns persist after accounting for the effects of plant diversity, climate and soil properties. Our findings challenge the conventional view that only conservative plant growth strategies promote ecosystem resistance and temporal stability, and instead highlight distinct pathways to stability across ecosystems—much like the many roads that once led to Rome.

    Explore related subjects
    Discover the latest articles and news in related subjects.

    Community ecology

    Ecosystem ecology

    Grassland ecology

    AcknowledgementsWe thank members of the He Lab for their assistance with field sampling and plant trait measurements.FundingN.H. was supported by the National Key R&D Program of China (2022YFF080210102) and the National Natural Science Foundation of China (32588202 and 32430067). L.J. was supported by the US National Science Foundation (DEB-2520925) and the Georgia Institute of Technology.Author informationAuthors and AffiliationsSchool of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USAPu Yan & Lin JiangKey Laboratory of Boreal Forest Ecosystem Conservation and Restoration, National Forestry and Grassland Administration, Northeast Forestry University, Harbin, ChinaNianpeng HeSouth China National Botanical Garden, Guangzhou, ChinaChangjin ChengMinistry of Education Key Laboratory of Ecology and Resource Use of the Mongolian Plateau and Inner Mongolia Key Laboratory of Grassland Ecology, School of Ecology and Environment, Inner Mongolia University, Hohhot, ChinaYonghui WangDepartment of Biodiversity, Ecology and Evolution, Faculty of Biological Science, Complutense University of Madrid, Madrid, SpainEnrique ValenciaCollege of Life and Environmental Sciences, Minzu University of China, Beijing, ChinaCongcong LiuDepartment of Ecology and Evolutionary Biology, University of California, Los Angeles, CA, USALawren SackAuthorsPu YanView author publicationsSearch author on:PubMed Google ScholarNianpeng HeView author publicationsSearch author on:PubMed Google ScholarChangjin ChengView author publicationsSearch author on:PubMed Google ScholarYonghui WangView author publicationsSearch author on:PubMed Google ScholarEnrique ValenciaView author publicationsSearch author on:PubMed Google ScholarCongcong LiuView author publicationsSearch author on:PubMed Google ScholarLawren SackView author publicationsSearch author on:PubMed Google ScholarLin JiangView author publicationsSearch author on:PubMed Google ScholarCorresponding authorsCorrespondence to
    Nianpeng He or Lin Jiang.Ethics declarations

    Competing interests
    The authors declare no competing interests.

    Additional informationPublisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary informationSupplementary Information (download PDF )Reporting Summary (download PDF )Transparent Peer Review File (download PDF )Source dataSource Data (download XLSX )Rights and permissions
    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
    Reprints and permissionsAbout this articleCite this articleYan, P., He, N., Cheng, C. et al. Contrasting plant trait strategies underpin ecosystem stability in temperate and alpine grasslands.
    Nat Commun (2026). https://doi.org/10.1038/s41467-026-78537-2Download citationReceived: 30 September 2025Accepted: 22 September 2026Published: 08 October 2026DOI: https://doi.org/10.1038/s41467-026-78537-2Share this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy shareable link to clipboard
    Provided by the Springer Nature SharedIt content-sharing initiative More

  • in

    The impact of two forms of environmental agency supply on public space use among older adults in nursing homes: a cluster-randomised field experiment

    AbstractEnvironmental control in nursing-home public spaces is usually treated as a matter of physical conditions, although residents often depend on staff to change those conditions. We conceptualise environmental agency supply as the opportunity created by institutional rules and environmental arrangements for residents to act on, or choose among, available conditions. We tested two forms of this supply in a mixed factorial, cluster-randomised field experiment involving 544 residents in 36 nursing homes across three climate zones in China. Direct control access was randomised between facilities, and preset microclimate choice was assigned within facilities as a crossover. Direct control access increased target-space occupancy by 6.36% points despite prompt staff responses in the comparison condition. Preset microclimate choice increased occupancy by 3.88% points relative to near-uniform conditions. Occupancy of other public spaces also increased, and no interaction between the two protocols was detected. These findings show that the organisation of environmental control can itself shape spatial use: authority to act and selectable environmental variation should be treated as design properties of institutional public space alongside physical conditions.

    Explore related subjects
    Discover the latest articles and news in related subjects.

    Ecology

    Environmental sciences

    Environmental social sciences

    FundingThis research received no external funding.Author informationAuthors and AffiliationsSchool of Architecture, Building and Design, Taylor’s University, 47100, Subang Jaya, MalaysiaYu Sun, Myzatul Aishah Binti Kamarazaly & Filzani Illia Binti IbrahimSchool of Art and Design, Jiangxi Institute of Fashion Technology, 103 Lihu Road, Xiangtang Development Zone, Nanchang, 330201, Jiangxi, ChinaYu SunAuthorsYu SunView author publicationsSearch author on:PubMed Google ScholarMyzatul Aishah Binti KamarazalyView author publicationsSearch author on:PubMed Google ScholarFilzani Illia Binti IbrahimView author publicationsSearch author on:PubMed Google ScholarCorresponding authorCorrespondence to
    Yu Sun.Ethics declarations

    Competing interests
    The authors declare no competing interests.

    Ethical considerations
    This study was reviewed and approved by the Human Research Ethics Committee of Taylor’s University (approval number HEC 2025/013). All methods were performed in accordance with the relevant guidelines and regulations, and the study was conducted in accordance with the Declaration of Helsinki. Written informed consent was obtained from all participants before enrolment. No identifying images or personal identifiers of participants are included in the manuscript or the supplementary materials.

    Additional informationPublisher’s noteSpringer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary InformationBelow is the link to the electronic supplementary material.Supplementary Material 1 (download DOCX )Rights and permissions
    Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
    Reprints and permissionsAbout this articleCite this articleSun, Y., Kamarazaly, M.A.B. & Ibrahim, F.I.B. The impact of two forms of environmental agency supply on public space use among older adults in nursing homes: a cluster-randomised field experiment.
    Sci Rep (2026). https://doi.org/10.1038/s41598-026-75075-1Download citationReceived: 24 July 2026Accepted: 05 October 2026Published: 07 October 2026DOI: https://doi.org/10.1038/s41598-026-75075-1Share this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy shareable link to clipboard
    Provided by the Springer Nature SharedIt content-sharing initiative
    KeywordsEnvironmental agency supplyPersonal environmental controlMicroclimate choiceNursing homesPublic space useCluster-randomised field experiment More