1. Scope and analytical basis of the review
This review combined large-scale scientometric mapping with a critical interpretation of the biological literature on fish responses to environmental change. Scientometric indicators were used to identify patterns of scholarly activity and association; they were not interpreted as direct evidence of biological causation, research quality, or policy influence17,18.
Bibliographic records were retrieved from the multidisciplinary citation database specified in Table 1. The database was selected because its standardized bibliographic records, cited-reference metadata, author affiliations, and historical coverage support longitudinal citation, collaboration, co-citation, and burst analyses19. Alternative databases may provide broader coverage of particular disciplines or regional journals, but differences in record structure, cited-reference standardization, and retrospective coverage can affect the reproducibility of large-scale citation-network analyses.
The final search was conducted on 12 December 2025 using the Topic field, which covered article titles, abstracts, author keywords, and database-generated keywords. The Boolean query was fish* AND (adapt* OR acclimati* OR toleran* OR plasticit* OR adjust* OR modif*).
Wildcards captured grammatical variants. Adaptation, acclimation, acclimatization, phenotypic plasticity, and tolerance were included as related retrieval concepts but were not treated as biologically identical. The broad stems adjust* and modif* increased search sensitivity but may also have retrieved peripheral or non-biological uses; this limitation is addressed in Section 7.
The initial search identified 100,516 records. Restriction to research articles retained 85,918 records and excluded 14,598 non-article document types. Restriction to English-language publications retained 84,452 records and excluded 1,466 non-English articles. Exclusion of 76 records assigned to the incomplete 2026 publication year left 84,376 records after database filtering. Duplicate detection removed 599 records, leaving 83,777 unique records. A further 163 records were excluded because complete cited-reference information required for citation-based analysis was unavailable. The final analytical dataset, therefore, comprised 83,614 research articles and 4,272,627 cited references (Figure 1). No eligible records were identified before 1975.
The annual publication series used in Section 3 was generated from the 84,376 records retained after database filtering. The final network, co-citation, keyword, and citation-burst analyses used the cleaned dataset of 83,614 articles. This distinction is maintained because annual publication counts were generated before duplicate removal and cited-reference completeness assessment.
Exported records were combined into a unified corpus before analysis. Duplicate records were removed before network construction. Metadata were reviewed for obvious inconsistencies in author names, institutional affiliations, keyword forms, and cited references. Where equivalence could be established unambiguously, spelling variants, abbreviations, singular and plural forms, and synonymous expressions were consolidated. Methodological and software-related citations were distinguished from biological, ecological, and conceptual references during interpretation because widely used analytical references can accumulate citations across many fields.
Network analyses covered January 1975 through December 2025 using one-year time slices. Term sources comprised Title, Abstract, Author, Keywords (DE), and Keywords Plus (ID). Link strength was calculated using cosine similarity, with the scope restricted to within slices. The Top N selection criterion retained the 50 most cited or occurring items from each annual slice. Pathfinder pruning and pruning of sliced networks were enabled, whereas minimum spanning tree and pruning of the merged network were disabled. The saved project configuration reported LRF = 2.5, L/N = 10, LBY = 5, and e = 1.0. Node type was selected according to the analysis; the verified project screenshot shows the institution selected for the institutional collaboration network. The analytical components included annual publication analysis; country, institution, and author co-authorship analysis; reference and journal co-citation analysis; keyword co-occurrence and burst analysis; and cluster and timeline visualization. Co-citation was interpreted as intellectual association rather than evidence that one publication caused a subsequent development20. Collaboration centrality was interpreted as structural position within the observed network rather than a direct measure of scientific quality, national strategy, or policy influence21. Cluster labels were generated algorithmically, and cluster cohesion was assessed using silhouette values22.
Annual publication growth was examined using yearly publication counts from the 84,376-record filtered corpus. A log-linear regression model was fitted using the natural logarithm of annual publication count as the response variable and publication year as the predictor. The 2025 count was excluded from model estimation because the search date may have produced incomplete annual and indexing coverage. Compound annual growth rates were calculated descriptively for 1975-2009, 2010-2018, and 2019-2024. A fixed-knot segmented log-linear model assessed whether slopes changed at 2010 and 201923. Neither slope change was statistically significant, so the periods were retained as post-analysis descriptive phases based on the combined publication and keyword patterns rather than treated as statistically confirmed growth regimes.
2. Global collaboration, institutional capacity, and knowledge equity
Country- and institution-level co-authorship analyses revealed a geographically broad but structurally uneven research landscape. Publication output was concentrated among a relatively small number of countries and organizations, while the strongest bridging positions were primarily occupied by contributors from North America and Western Europe. Publication count and betweenness centrality were interpreted separately because they represent different dimensions of activity: publication count reflects productivity, whereas betweenness centrality indicates the extent to which a node connects otherwise weakly connected parts of a network21. The geographical distribution of publications, country-level co-authorship network, leading institutions, and institutional collaboration network are presented in Figure 2.
The United States was the largest contributor, with 23,990 publications, and occupied the strongest bridging position in the country-level network, with a betweenness centrality of 0.37. China ranked second in publication output with 10,841 publications, but had a centrality value of 0.01. France and Canada contributed 4,565 and 7,005 publications, respectively, but had comparatively high centrality values of 0.26 and 0.15. England, Belgium, and Denmark also occupied bridging positions. These findings show that high productivity does not necessarily correspond to a strong bridging position and that moderate-output countries can contribute substantially to network connectivity.
Malaysia contributed 573 publications, Indonesia 414, and Nigeria 84, but none had a centrality value exceeding 0.01. These values indicate limited bridging roles relative to the most connected countries, but they do not demonstrate an absence of international collaboration. A country may participate in bilateral or regional partnerships while still having low betweenness centrality when those collaborations do not connect otherwise separate parts of the global network. The overall pattern is consistent with a core-periphery structure reported in broader analyses of international scientific collaboration24.
The network results do not establish why particular countries occupy central or peripheral positions. Low centrality cannot demonstrate insufficient funding, political barriers, domestic research orientation, limited expertise, or exclusion from international partnerships. Similarly, high centrality does not prove that a country determines theoretical frameworks, methodological standards, or research agendas. Evaluating those mechanisms would require evidence on research expenditure, grant leadership, corresponding authorship, project coordination, data ownership, and long-term institutional partnerships.
Institutional output was similarly concentrated. The University of California system contributed 2,747 publications, followed by the French National Centre for Scientific Research with 2,344. The Chinese Academy of Sciences and the United States Department of the Interior each contributed 1,541 publications. Other major contributors included the National Oceanic and Atmospheric Administration, the State University System of Florida, the United States Geological Survey, Harvard University, and the Spanish National Research Council. The Chinese Academy of Sciences is identified separately as an East Asian national research organization and is not grouped with United States institutions.
The prominence of universities, national research organizations, and government agencies demonstrates that fish-adaptation research is conducted across academic, environmental-monitoring, and resource-management settings. However, co-authorship data alone cannot establish whether these institutions translated research into management guidelines, thermal benchmarks, habitat assessments, or fisheries policies.
The peripheral position of many countries and institutions from Southeast Asia, Africa, South Asia, and parts of Latin America has potential implications for knowledge equity. Global research expenditure and infrastructure remain unevenly distributed, even as international scientific collaboration has increased25. Equitable participation requires more than inclusion as a sampling location or data provider; it involves meaningful roles in research design, project leadership, corresponding authorship, data interpretation, infrastructure development, and management application. The present co-authorship network cannot determine whether these dimensions were equitable.
The geographical pattern may also affect ecological representation. Climate-related changes in fisheries productivity and species distributions vary strongly among regions, with important consequences for low-latitude fisheries and communities dependent on aquatic resources26. Nevertheless, underrepresentation of species, habitats, and life stages cannot be demonstrated from centrality alone and requires direct evidence-gap mapping. The results support co-led international projects, reciprocal postgraduate training, shared experimental and computational infrastructure, transparent data-governance agreements, and sustained North-South and South-South partnerships.
3. Temporal growth and conceptual evolution of fish-adaptation research
Annual publication output and keyword patterns revealed substantial long-term growth and conceptual reorganization between 1975 and 2025. For interpretation, the field was organized into a foundational physiological phase from 1975 to 2009, a climate-oriented expansion phase from 2010 to 2018, and a high-output molecular and integrative phase from 2019 to 2025. These periods were defined post-analysis by combining publication and keyword patterns and were not imposed before analysis. The keyword co-occurrence structure across the three periods is presented in Figure 3.
Foundational physiological phase, 1975–2009
Annual publication output increased from 23 records in 1975 to 1,880 in 2009. This period contributed 22,554 records, representing approximately 26.7% of the annual publication corpus. Its endpoint-based compound annual growth rate was 13.83%, although this value was influenced by the very small publication count at the beginning of the series.
The strongest reported keyword bursts included the historical taxonomic term Salmo gairdneri, trout Salmo gairdneri, plasma, and smoltification. These patterns indicate a strong early concentration on salmonid experimental systems, endocrine stress responses, circulating physiological indicators, osmoregulatory transitions, and life-history changes associated with movement between freshwater and marine environments.
This research provided important foundations for understanding how corticosteroids, metabolic processes, ionic regulation, and integrated physiological systems contribute to stress responses and homeostasis in fish27,28. The prominence of salmonid terminology demonstrates the visibility of salmonid research within the indexed literature rather than the absence of other taxa. Nevertheless, the strong taxonomic concentration raises questions about how well frameworks derived from temperate salmonids apply to tropical, polar, estuarine, and non-salmonid fishes.
Climate-oriented expansion phase, 2010–2018
Annual publication output increased from 2,087 records in 2010 to 3,658 in 2018. This period contributed 26,942 records, equivalent to approximately 31.9% of the annual publication corpus, and had an endpoint-based compound annual growth rate of 7.27%.
The keyword landscape broadened beyond species-specific physiology. Climate change exhibited a strong reported burst during 2016–2018 and was accompanied by morphology, abundance, fish consumption, biodiversity, assemblage, fisheries, and management terminology. This represented an expansion rather than the replacement of physiological research. Mechanistic evidence was increasingly interpreted within broader contexts involving warming, oxygen limitation, acidification, habitat alteration, species redistribution, and fisheries vulnerability2,29.
The increase in research activity occurred during growing international scientific and policy attention to climate change, biodiversity loss, and ecosystem resilience. However, the present analysis did not include annual policy, funding, or climate-event variables and therefore cannot demonstrate that specific agreements, extreme events, or funding initiatives caused the increase in publications.
High-output molecular and integrative phase, 2019–2025
The third phase contributed 34,880 records, representing approximately 41.3% of the annual publication corpus. Annual output increased from 4,309 records in 2019 to a complete-year peak of 5,341 in 2024. The endpoint-based compound annual growth rate from 2019 to 2024 was 4.39%, lower than in the preceding periods. The term high-output phase is therefore more accurate than the acceleration phase. The 2025 count of 5,335 records should be interpreted cautiously because the search was conducted on 12 December 2025 and may not include publications released or indexed during the final weeks of the year.
The strongest verified Phase 3 keyword bursts were rainbow trout, climate change, gene expression, and disease. The prominence of gene expression indicates increased attention to molecular regulation, environmentally induced phenotypic responses, and pathways connecting environmental exposure with organismal performance. Molecular approaches extended rather than replaced the physiological foundations established during earlier periods9,10. The disease burst indicates increasing recognition that fish responses are influenced by interactions between abiotic stress and biotic pressures30.
Despite greater use of molecular and genomic approaches, rainbow trout remained the strongest Phase 3 burst term. Its continued prominence indicates that increased analytical resolution has not necessarily been accompanied by equivalent taxonomic diversification. Established models and aquaculture species provide extensive genomic resources and standardized protocols, but excessive reliance on a limited number of taxa may constrain the wider applicability of adaptation frameworks.
Interpretation of keyword evolution
The temporal keyword trajectory indicates conceptual accumulation and integration rather than simple replacement. Experimental physiology and salmonid research established much of the mechanistic foundation. Climate-oriented research expanded analysis toward populations, communities, ecosystems, fisheries, and human use. More recent molecular, disease, oxidative-stress, and genomic approaches added greater analytical resolution to these physiological and ecological perspectives.
The absence of DNA, hypoxia, or early life stage as an individual major burst term should not be interpreted as evidence that those topics were absent. Relevant studies may be distributed across gene expression, genomics, transcriptomics, genetic variation, metabolism, thermal tolerance, species-specific terms, and environmental DNA. Burst detection identifies unusually rapid increases in attention and does not provide an exhaustive measure of all consistently important concepts17,18. Complete annual or phase-specific frequencies were not available for every keyword, so burst strength was interpreted as evidence of rapidly increasing attention rather than proportional importance, quality, or biological significance.
Statistical assessment of publication growth
A log-linear regression model fitted to annual publication counts from 1975 to 2024 indicated an estimated exponential annual growth rate of 9.02% (95% confidence interval, 8.49%–9.56%; R2 = 0.963; p < 0.001). A fixed-knot segmented log-linear model did not identify statistically significant slope changes at 2010 (p = 0.292) or 2019 (p = 0.676). The three periods should therefore be retained as descriptive post-analysis phases rather than presented as statistically confirmed inflection points23. Annual publication output and the three descriptive periods are shown in Figure 4.
Citation-density and citation-half-life comparisons were excluded because annual total-citation data and a formal citation half-life calculation were not available. Older publications have also had more time to accumulate citations than recent publications. Scientometric indicators should be interpreted as measures of scholarly activity, visibility, and knowledge circulation rather than direct measures of scientific quality31.
4. Citation-burst structure and influential references
Reference citation-burst analysis identifies publications that experience a rapid increase in citation attention during a particular period. Burst strength differs from total citation count: a reference can exhibit a strong burst because its use increased sharply over a short interval, even when it is not among the most cited references across the complete dataset. Citation bursts are therefore useful for identifying periods of concentrated attention and methodological transitions, but should not be interpreted as direct measures of research quality, biological importance, or causal influence17,18.
The timeline and duration of the 17 strongest cited-reference bursts are presented in Table 1. For interpretation, verified records were grouped into methodological and computational references and biological or conceptual references. Unresolved or malformed citation records were not assigned biological interpretations.
Methodological and computational visibility
The burst structure was dominated by methodological and computational references rather than empirical studies focused specifically on fish-adaptation mechanisms. Multiple version-specific records associated with general statistical computing and community ecology analysis appeared separately, increasing the apparent concentration of methodological influence. This pattern indicates that contemporary fish-adaptation research depends on shared analytical infrastructure for statistical modeling, multivariate analysis, genomic processing, and phylogenetic inference. Repeated citation of a method does not mean that it made a biological discovery or contributed to a fish-specific theory.
Prominent methodological references included molecular-evolution frameworks32,33 and a sequence-read preprocessing method34. These tools are used across many biological disciplines and can attract citations from a much wider literature than fish-adaptation research alone. Version-specific records should be retained transparently in the raw analytical output but grouped conceptually in the narrative to avoid treating software versions as independent intellectual advances.
Biological and conceptual references
Verified biological burst references represented complementary dimensions of fish-adaptation science. Evans and colleagues synthesized the multifunctional role of the fish gill in gas exchange, osmoregulation, acid-base regulation, and nitrogenous-waste excretion35. Eliason and colleagues demonstrated population-level differences in cardiorespiratory performance and thermal tolerance among sockeye salmon populations36. Clark and colleagues evaluated the measurement, interpretation, and limitations of aerobic scope in climate-change research11. Dahlke and colleagues identified reproductive adults and embryos as particularly temperature-sensitive life stages5, while Alfonso and colleagues connected rising temperatures with neuroendocrine, metabolic, immune, and oxidative-stress responses37.
Genomic references reflected expansion toward the architecture of adaptive evolution. Research on threespine sticklebacks demonstrated repeated ecological adaptation through shared standing genetic variation and regulatory and coding changes38. Comparative genomic analysis of African cichlids showed how gene duplication, regulatory divergence, transposable elements, microRNAs, and standing variation contributed to adaptive radiation39. A broader freshwater synthesis linked climate change, disease, contaminants, salinization, cumulative stressors, environmental DNA, and conservation responses40.
Methodological citation dominance creates a form of citation bias because analytical references are cited whenever a procedure is used, whereas empirical biological studies may be relevant only to particular taxa, life stages, habitats, or stressors. Journal visibility, accessibility, international collaboration, and alignment with rapidly developing topics can also influence citation bursts. The present analysis did not evaluate submission behavior or publication pressure and therefore cannot determine whether incentives to publish in high-impact journals caused the observed pattern31.
Citation-burst analysis and co-citation clustering answer different questions. Bursts identify references that experienced rapid increases in attention, whereas co-citation clusters identify groups of references that were cited together and represent broader intellectual or thematic communities. The relative absence of habitat-specific, early-life-stage, or multi-stressor studies from the top burst set does not prove that these topics are absent or unimportant; they may be distributed across many context-specific publications.
5. Thematic co-citation structure and biological evolution of fish-adaptation research
Reference co-citation analysis of 4,272,627 cited references identified 26 thematic clusters. A co-citation cluster comprises references frequently cited together and therefore reflects a shared intellectual, methodological, taxonomic, or biological context. Cluster labels were generated algorithmically and interpreted by examining associated cited and citing literature (Figure 5). They should be treated as concise thematic descriptors rather than exact classifications of every publication within a cluster.
The clusters discussed below had silhouette values ranging from 0.860 to 0.996. Silhouette indicates internal cohesion and separation but does not measure research quality, biological importance, or the reliability of every interpretation. Cluster size represents the number of cited-reference nodes grouped together and is not a direct measure of the total number of studies on a topic. Average publication year provides an approximate temporal position rather than the exact origin of a theme.
Model species as historical research anchors
Species-specific clusters formed some of the largest components of the network. Brown trout was represented by Cluster #0 (337 nodes; silhouette 0.914; average year 1996), with literature on growth efficiency, metabolic regulation, temperature-photoperiod interactions, acclimation, and energetic performance. Sockeye salmon formed Cluster #1 (326 nodes; silhouette 0.885; average year 2007), connecting cardiorespiratory performance, aerobic scope, migration ecology, thermal sensitivity, and population-specific vulnerability36. North Sea cod formed Cluster #4 (228 nodes; silhouette 0.919; average year 2004), linking marine warming, redistribution, productivity, environmental variability, and fisheries. Largemouth bass formed Cluster #5 (206 nodes; silhouette 0.895; average year 2020), with research on chronic and extreme warming, upper thermal limits, immune responses, and molecular regulation.
These clusters demonstrate sustained visibility of well-studied fishes, likely reflecting historical datasets, ecological or commercial relevance, established husbandry methods, and repeated use in physiological experiments. Their prominence does not prove that these species determined the theoretical direction of the field. Concentration on temperate and commercially important fishes may limit transferability to tropical, polar, deep-sea, estuarine, short-lived, and other non-model species.
Physiological, behavioral, and cellular response mechanisms
Thermal tolerance formed Cluster #3 (266 nodes; silhouette 0.860; average year 2017), linking critical thermal limits, acclimation capacity, metabolic compensation, aerobic performance, and warming vulnerability. No single threshold captures the complete adaptive capacity of a population because thermal limits, aerobic performance, reproduction, growth, behavior, disease susceptibility, and early-life-stage survival can respond differently2,5,11.
Bioenergetics modeling formed Cluster #10 (153 nodes; silhouette 0.958; average year 1989), connecting energy acquisition and allocation with growth, reproduction, survival, and population consequences. Seawater stage formed Cluster #8 (185 nodes; silhouette 0.936; average year 2001) and is better interpreted as a life-history and osmoregulatory theme involving smoltification, ion regulation, energetic costs, and developmental timing.
Older and smaller clusters represented caudal myotome function, phagocytic function, antifreeze proteins, movement behavior, and predation-risk seasonality. A smaller size indicates fewer co-cited reference nodes under a specific label, not biological unimportance or the absence of a wider literature. Together, the mechanistic clusters show that tissue function, immune responses, energetic balance, locomotion, osmoregulation, and thermal tolerance interact to determine organismal performance.
Climate-driven environmental stressors
Marine heatwaves formed Cluster #2 (325 nodes; average year 2018), one of the largest and most recent clusters. The associated literature addressed physiological disruption, species redistribution, ecosystem change, fisheries effects, and repeated or compound warming events. Marine heatwaves can expose organisms to rapid and extreme temperatures that exceed short-term acclimation capacity6,7. The co-citation map labels Cluster #13 as lake characteristics; therefore, it was not interpreted as a second marine-heatwave cluster.
Ocean acidification formed Cluster #9 (165 nodes; silhouette 0.902; average year 2013), with research on respiratory performance, acid-base regulation, behavior, development, early-life stages, and interactions with warming. Fish can regulate internal acid-base conditions, but this capacity varies among species, stages, exposure durations, and interacting stressors4. Hypoxia did not appear as a comparably dominant labeled cluster despite its importance. This absence does not indicate a lack of hypoxia research; related studies may be distributed across thermal tolerance, gill physiology, metabolism, marine heatwaves, acidification, and species-specific clusters.
Place-based ecological themes and monitoring methods
The Eastern Bering Sea formed Cluster #6 (198 nodes; silhouette 0.896; average year 2010), linking fisheries oceanography, climate variability, redistribution, ecosystem productivity, and population dynamics. Cluster #13 was labeled lake characteristics and connected thermal regimes, habitat structure, seasonal variation, species composition, and changes in ecological dominance. Place-based clusters reflect research communities in which environmental change, population response, fisheries, and monitoring are considered together; they do not imply that every study directly measured biological adaptation.
Environmental DNA formed Cluster #14 (76 nodes; silhouette 0.962; average year 2018), representing a recent methodological theme associated with non-invasive biodiversity assessment, species detection, community monitoring, and distributional change. Environmental DNA is not itself an adaptation mechanism; it is a monitoring approach that can document where populations and communities are changing. Its interpretation remains affected by DNA production, transport, degradation, detection probability, sampling design, and reference-database quality12,13. Its greatest value will arise when integrated with physiological measurements, conventional surveys, fisheries observations, genomic information, and environmental data.
Cross-cluster interpretation and evidence gaps
The co-citation structure reveals layered rather than strictly linear development. Older clusters were associated with tissue function, energetic performance, osmoregulation, behavior, cold adaptation, and established model species. Intermediate clusters connected physiology with migration, fisheries, oceanography, population dynamics, and climate sensitivity. Recent clusters emphasized thermal extremes, molecular regulation, immune responses, and field monitoring. This pattern supports the temporal changes identified through keyword analysis, but temporal proximity does not demonstrate that one cluster caused another to emerge.
The network simultaneously contains species-centered, mechanistic, climate- and ecosystem-oriented, and monitoring communities. The field has not shifted completely from physiology to genomics or from laboratory research to ecosystem science; newer approaches have accumulated around established foundations. Oxidative stress was visible in recent keyword patterns but did not form a major independent co-citation label, suggesting that it may be distributed across thermal stress, metabolism, immunity, disease, pollution, and species-specific literature. This interpretation requires direct examination of the underlying publications.
Important gaps include taxonomic concentration, fragmented life-stage evidence, limited integration of multiple stressors, and incomplete linkage between monitoring and mechanism. Embryos, larvae, juveniles, and reproductive adults may be particularly sensitive, but their evidence is distributed across thermal-tolerance, acidification, reproduction, growth, and species-specific clusters5. More realistic experiments are needed to connect molecular and physiological mechanisms with demographic and ecosystem outcomes.
6. Emerging research frontiers, evidence gaps, and management perspectives
The combined keyword, citation-burst, and co-citation patterns indicate movement toward greater mechanistic resolution, wider ecological application, and closer integration between laboratory and field evidence. Recent themes include gene expression, oxidative stress, disease, thermal extremes, adaptive capacity, environmental DNA, and climate-sensitive life stages. These patterns identify areas of increasing scholarly attention but do not establish that they are biologically more important than themes with weaker bursts or smaller clusters.
Integrating molecular responses with adaptive capacity
The increased prominence of gene expression and related molecular terminology reflects wider use of transcriptomic and genomic approaches to examine regulatory pathways, biomarkers, and population differences. Altered gene expression alone does not demonstrate evolutionary adaptation; it may represent transient physiological response, acclimation, developmental plasticity, or cellular damage. Molecular responses should therefore be connected with thermal tolerance, aerobic capacity, osmoregulation, immunity, growth, behavior, reproduction, survival, and population persistence. Common-environment, reciprocal-exposure, quantitative-genetic, population-genomic, and multigenerational approaches are needed to distinguish environmentally induced plasticity from heritable change8,9.
Multi-omics approaches may strengthen this integration by combining genomics, transcriptomics, proteomics, metabolomics, and physiological measurements. Their value lies not simply in generating more molecular information but in determining whether changes across biological levels converge on consistent mechanisms and predict organismal or demographic outcomes. Oxidative stress should be measured alongside whole-organism performance rather than treated as an independent indicator of adaptive capacity.
Climate, disease, and host-pathogen interactions
The recent disease burst indicates growing attention to interactions between climate stress and biotic challenges. Temperature, oxygen availability, salinity, and other environmental factors may influence pathogen development, transmission, geographic range, and virulence while simultaneously affecting fish immunity, energetic condition, and susceptibility30. Future work should use factorial designs that combine environmental stress with pathogen exposure and should include longitudinal monitoring because acute immune responses may not predict long-term disease or population outcomes.
Sensitive life stages and demographic bottlenecks
Embryos, larvae, juveniles, spawning adults, and non-reproductive adults may differ substantially in thermal tolerance, oxygen demand, osmoregulatory capacity, mobility, and behavioral avoidance. Evidence across hundreds of species identifies embryos and reproductive adults as particularly temperature sensitive5. Future experiments should encompass multiple stages within the same species and measure immediate mortality, development, growth, deformity, behavior, reproductive performance, carry-over effects, recovery, and recruitment. Linking these outcomes with demographic models would clarify whether sensitivity at one stage creates a population bottleneck.
Multiple stressors, hypoxia, and climate extremes
Natural populations rarely experience stressors in isolation. Warming commonly interacts with oxygen limitation, acidification, salinity change, food availability, pollution, pathogens, habitat degradation, and fishing pressure. Combined effects may be additive, antagonistic, or synergistic and depend on exposure sequence, duration, intensity, fluctuation, and recovery. Hypoxia-related studies may be distributed across gill physiology, metabolism, thermal tolerance, marine heatwaves, acidification, disease, and species-specific clusters. Fish responses include behavioral avoidance, ventilatory and cardiovascular adjustment, metabolic suppression, gene regulation, and, in some species, evolved hypoxia tolerance15.
The interaction between warming and hypoxia is particularly important because warming can increase metabolic oxygen demand while reducing oxygen solubility, and coastal eutrophication and stratification can intensify low-oxygen conditions3. Marine-heatwave experiments should reproduce event magnitude, duration, recurrence, season, previous thermal history, and recovery periods rather than relying solely on gradual changes in mean temperature6,7.
Taxonomic, geographical, and habitat representation
The persistence of salmonids, cod, largemouth bass, and other established species indicates the value of model and commercially important fishes but may limit transferability to species with different evolutionary histories and environmental conditions. Future research should broaden representation across tropical, polar, estuarine, deep-sea, freshwater, migratory, short-lived, and data-limited fishes, particularly species supporting small-scale fisheries and local livelihoods.
Collaboration centrality alone cannot establish that ecosystems or species are quantitatively underrepresented. A formal evidence-gap analysis should map study locations, taxa, climatic zones, habitats, life stages, and stressors. Addressing geographical imbalance requires local leadership, corresponding authorship, shared infrastructure, postgraduate training, data ownership, and participation in interpretation and management translation.
Long-term monitoring and environmental DNA
Environmental DNA can improve detection of rare, cryptic, migratory, or widely dispersed species and support repeated sampling where conventional surveys are difficult12,13. It is not a direct measure of physiological or evolutionary adaptation. Changes in detection can indicate altered occurrence, distribution, or community composition, but do not identify the mechanism responsible. Its greatest contribution will arise when integrated with conventional surveys, tagging, fisheries data, physiological measurements, genomic analysis, and environmental observations.
Long-term monitoring is necessary because short experiments may not capture acclimatization, recovery, repeated extremes, interannual variability, or demographic response. Coordinated time series combining temperature, oxygen, salinity, productivity, habitat condition, disease occurrence, distribution, abundance, and biological traits would provide a stronger basis for evaluating resistance, recovery, and long-term adaptive capacity.
Fisheries, aquaculture, and conservation implications
For fisheries, adaptation research can help identify climate-sensitive species, populations, seasons, and life stages. Information on thermal limits, oxygen sensitivity, reproductive windows, distributional shifts, and recovery may support vulnerability assessment, spatial planning, seasonal management, stock assessment, and interpretation of changes in catch or recruitment. Climate-driven redistribution also requires coordination across administrative boundaries because populations and fisheries may shift beyond historical management areas41.
For aquaculture, physiological and molecular evidence can inform site selection, temperature and oxygen management, broodstock choice, disease surveillance, feeding strategies, and strain evaluation. Selection for one favorable trait may create trade-offs in growth, reproduction, disease resistance, or performance under other stressors, so climate-resilient aquaculture should evaluate multiple traits and realistic environmental combinations.
For conservation, adaptation research can assist in identifying thermal refugia, climate-sensitive nursery areas, habitat corridors, vulnerable populations, and locations where restoration could reduce non-climatic stress. Management application must treat uncertainty explicitly. Laboratory thresholds should not be transferred directly into policy without considering population variation, acclimation history, environmental fluctuation, ecological context, and data quality.
Research and management priorities
The combined evidence supports five connected priorities: validate molecular and cellular responses against organismal and demographic outcomes; conduct full-life-cycle and multi-stressor studies under realistic exposure sequences; expand taxonomic, geographical, habitat, and life-stage representation through equitable partnerships; integrate long-term field monitoring, including environmental DNA, with physiology, genomics, fisheries observations, and environmental data; and translate adaptation evidence into decision-relevant indicators while preserving uncertainty and population-specific variation.
7. Limitations of the scientometric evidence
Scientometric analysis provides a broad overview of how a research field is organized and how scholarly attention changes over time. It does not provide a complete inventory of all relevant evidence, directly assess the methodological quality of individual studies, quantify biological effect sizes, or establish biological and policy causation.
Database coverage and document selection
The corpus was retrieved from one multidisciplinary citation index. The selected database provides standardized records and cited-reference metadata suited to longitudinal network analysis19, but no database provides complete coverage. Relevant publications may be indexed only in other international, discipline-specific, or regional databases. Technical reports, governmental documents, theses, books, and other grey literature may be absent. Database coverage also changes over time and may contribute to apparent publication growth.
Only English-language research articles were retained. This improved document consistency but may have reduced representation of research from regions where English is not the primary publication language and may have excluded influential reviews, proceedings, books, technical findings, and locally relevant fisheries literature. The search was completed on 12 December 2025, so 2025 may be incomplete, and 2024 was treated as the latest complete year in growth modeling.
Search breadth and conceptual specificity
The search strategy intentionally included adaptation, acclimation, acclimatization, tolerance, plasticity, adjustment, and modification to capture a broad domain. The stems adjust* and modif* may have retrieved technical, statistical, morphological, or non-adaptive contexts. Conversely, relevant studies that did not use the selected terms in titles, abstracts, or keywords may have been omitted. The analysis was a field-mapping review rather than a conventional systematic review based on full-text eligibility assessment of all 83,614 articles, so peripheral studies may remain in the corpus.
Dataset stages and metadata quality
The annual publication series used 84,376 records retained after database filtering, whereas citation and network analyses used the final 83,614 articles. The 762 records removed after database filtering represented approximately 0.90% of the filtered corpus and were unlikely to alter the overall five-decade trajectory substantially, although small annual differences may remain. Records without completely cited references were excluded from citation-based analyses and may be disproportionately associated with older or less standardized metadata.
Complete disambiguation of authors, institutions, keywords, and references is difficult in a dataset of this size. The same author or institution may appear under multiple names, while different entities may be combined incorrectly. Keyword analysis is affected by synonyms, abbreviations, historical terminology, taxonomic revisions, and singular or plural forms. Cited-reference formatting differences and version-specific software citations can split one reference into several records. The citation-burst output contained unresolved, duplicated, or malformed records, demonstrating that automated reference matching was not fully error-free.
Dependence on analytical parameters
Time-slice length, node-selection criteria, pruning algorithms, burst settings, and clustering procedures influence which nodes and links remain visible. One-year slices provide temporal detail but can increase fragmentation in low-output years. Pruning improves readability but may reduce visibility of weak, specialized, or emerging relationships. Algorithmic cluster labels are thematic summaries, not exact taxonomies. Silhouette values indicate cohesion and separation, not biological importance, study quality, or theoretical validity17,18.
Citation, keyword, and collaboration bias
Citation indicators are affected by publication age, field size, journal visibility, accessibility, collaboration, self-citation, and disciplinary citation practices. Methodological references may receive high counts or strong bursts because they are cited whenever an analytical procedure is used. Citation prominence should therefore not be interpreted as a direct measure of biological importance or quality31. Keyword bursts are also influenced by corpus growth, indexing, author terminology, synonym use, and normalization. The absence of a term from the strongest bursts or largest clusters does not demonstrate that the topic is absent from the literature.
Country and institutional co-authorship networks identify publication activity and structural connectivity, not funding capacity, scientific leadership, project ownership, equity, or policy influence. Betweenness centrality indicates a bridging position within the observed network, but does not show how funding, labor, data ownership, corresponding authorship, or decision-making were distributed. Geographical centrality also does not directly quantify taxonomic or habitat underrepresentation.
Temporal and biological inference
The periods 1975–2009, 2010–2018, and 2019–2025 were retained as interpretive phases. Segmented modeling did not identify statistically significant publication-growth slope changes at 2010 or 2019. Publication growth cannot be attributed directly to climate events, policy agreements, or funding programs because those variables were not included. Similarly, temporal proximity among clusters or citation bursts does not demonstrate causal succession.
This review maps research activity and critically interprets associated biological literature. It is not a formal risk-of-bias assessment, systematic review of experimental outcomes, or meta-analysis of effect sizes. It cannot determine which species, populations, or life stages are most vulnerable solely from publication counts, bursts, or co-citation clusters. Management implications should therefore be treated as research and policy considerations rather than prescriptions directly validated by scientometric analysis.
Despite these limitations, the large and temporally extensive dataset provides a useful overview of the development, collaboration structure, terminology, citation patterns, and thematic organization of fish-adaptation research. Future syntheses could integrate multiple databases, non-English and grey literature, systematic full-text assessment, normalized keyword frequencies, improved entity disambiguation, and meta-analysis of biological outcomes.