Climate change induces shifts in coastal Baltic Sea surface water microorganism stress and photosynthesis gene expression.

The world's oceans are challenged by climate change linked warming with typically highly populated coastal areas being particularly susceptible to these effects. Many studies of climate change on the marine environment use large, short-term temperature manipulations that neglect factors such as long-term adaptation and seasonal cycles. In this study, a Baltic Sea 'heated' bay influenced by thermal discharge since the 1970s from a nuclear reactor (in relation to an unaffected nearby 'control' bay) was used to investigate how elevated temperature impacts surface water microbial communities and activities. 16S rRNA gene amplicon based microbial diversity and population structure showed no difference in alpha diversity in surface water microbial communities, while the beta diversity showed a dissimilarity between the bays. Amplicon sequencing variant relative abundances between the bays showed statistically higher values for, e.g., Ilumatobacteraceae and Burkholderiaceae in the heated and control bays, respectively. RNA transcript-derived activities followed a similar pattern in alpha and beta diversity with no effect on Shannon's H diversity but a significant difference in the beta diversity between the bays. The RNA data further showed more elevated transcript counts assigned to stress related genes in the heated bay that included heat shock protein genes dnaKJ, the co-chaperonin groS, and the nucleotide exchange factor heat shock protein grpE. The RNA data also showed elevated oxidative phosphorylation transcripts in the heated (e.g., atpHG) compared to control (e.g., atpAEFB) bay. Furthermore, genes related to photosynthesis had generally higher transcript numbers in the control bay, such as photosystem I (psaAC) and II genes (psbABCEH). These increased stress gene responses in the heated bay will likely have additional cascading effects on marine carbon cycling and ecosystem services.

Baltic Sea coastal sediment-bound eukaryotes have increased year-round activities under predicted climate change related warming.

Climate change related warming is a serious environmental problem attributed to anthropogenic activities, causing ocean water temperatures to rise in the coastal marine ecosystem since the last century. This particularly affects benthic microbial communities, which are crucial for biogeochemical cycles. While bacterial communities have received considerable scientific attention, the benthic eukaryotic community response to climate change remains relatively overlooked. In this study, sediments were sampled from a heated (average 5°C increase over the whole year for over 50 years) and a control (contemporary conditions) Baltic Sea bay during four different seasons across a year. RNA transcript counts were then used to investigate eukaryotic community changes under long-term warming. The composition of active species in the heated and control bay sediment eukaryotic communities differed, which was mainly attributed to salinity and temperature. The family level RNA transcript alpha diversity in the heated bay was higher during May but lower in November, compared with the control bay, suggesting altered seasonal activity patterns and dynamics. In addition, structures of the active eukaryotic communities varied between the two bays during the same season. Hence, this study revealed that long-term warming can change seasonality in eukaryotic diversity patterns. Relative abundances and transcript expression comparisons between bays suggested that some taxa that now have lower mRNA transcripts numbers could be favored by future warming. Furthermore, long-term warming can lead to a more active metabolism in these communities throughout the year, such as higher transcript numbers associated with diatom energy production and protein synthesis in the heated bay during winter. In all, these data can help predict how future global warming will affect the ecology and metabolism of eukaryotic community in coastal sediments.

Genetic structure, UV-vision, wing coloration and size coincide with colour polymorphism in Fabriciana adippe butterflies.

Colour polymorphisms have long served as model systems in evolutionary studies and continue to inform about processes involved in the origin and dynamics of biodiversity. Modern sequencing tools allow for evaluating whether phenotypic differences between morphs reflect genetic differentiation rather than developmental plasticity, and for investigating whether polymorphisms represent intermediate stages of diversification towards speciation. We investigated phenotypic and genetic differentiation between two colour morphs of the butterfly Fabriciana adippe using a combination of ddRAD-sequencing and comparisons of body size, colour patterns and optical properties of bright wing spots. The silvery-spotted adippe form had larger and darker wings and reflected UV light, while the yellow cleodoxa form displayed more green scales and reflected very little UV, showcasing that they constitute distinct and alternative integrated phenotypes. Genomic analyses revealed genetic structuring according to source population, and to colour morph, suggesting that the phenotypic differentiation reflects evolutionary modifications. We report 17 outlier loci associated with colour morph, including ultraviolet-sensitive visual pigment (UVRh1), which is associated with intraspecific communication and mate choice in butterflies. Together with the demonstration that the wings of the adippe (but essentially not the cleodoxa) morph reflect UV light, that UV reflectance is higher in females than males and that morphs differ in wing size, this suggests that these colour morphs might represent genetically integrated phenotypes, possibly adapted to different microhabitats. We propose that non-random mating might contribute to the differentiation and maintenance of the polymorphism.

Effects of environmental translocation and host characteristics on skin microbiomes of sun-basking fish.

Variation in the composition of skin-associated microbiomes has been attributed to host species, geographical location and habitat, but the role of intraspecific phenotypic variation among host individuals remains elusive. We explored if and how host environment and different phenotypic traits were associated with microbiome composition. We conducted repeated sampling of dorsal and ventral skin microbiomes of carp individuals (Cyprinus carpio) before and after translocation from laboratory conditions to a semi-natural environment. Both alpha and beta diversity of skin-associated microbiomes increased substantially within and among individuals following translocation, particularly on dorsal body sites. The variation in microbiome composition among hosts was significantly associated with body site, sun-basking, habitat switch and growth, but not temperature gain while basking, sex, personality nor colour morph. We suggest that the overall increase in the alpha and beta diversity estimates among hosts were induced by individuals expressing greater variation in behaviours and thus exposure to potential colonizers in the pond environment compared with the laboratory. Our results exemplify how biological diversity at one level of organization (phenotypic variation among and within fish host individuals) together with the external environment impacts biological diversity at a higher hierarchical level of organization (richness and composition of fish-associated microbial communities).

Anthropogenic Influence on Moth Populations: A Comparative Study in Southern Sweden.

As moths are vital components of ecosystems and serve as important bioindicators, understanding the dynamics of their communities and the factors influencing these dynamics, such as anthropogenic impacts, is crucial to understand the ecological processes. Our study focuses on two provinces in southern Sweden, Västergötland and Småland, where we used province records from 1974 to 2019 in combination with light traps (in 2020) to record the presence and abundance of moth species, subsequently assessing species traits to determine potential associations with their presence in anthropogenically modified landscapes. This study design provides a unique opportunity to assess temporal changes in moth communities and their responses to shifts in environmental conditions, including anthropogenic impacts. Across the Västergötland and Småland provinces in Sweden, we recorded 776 moth taxa belonging to fourteen different taxonomic families of mainly Macroheterocera. We captured 44% and 28% of the total moth species known from these provinces in our traps in Borås (Västergötland) and Kalmar (Småland), respectively. In 2020, the species richness and abundance were higher in Borås than in Kalmar, while the Shannon and Simpson diversity indices revealed a higher species diversity in Kalmar. Between 1974 and 2019, the colonisation rates of the provinces increased faster in Småland. Ninety-three species were found to have colonised these provinces since 1974, showing that species richness increased over the study period. We reveal significant associations between the probability of a species being present in the traps and distinct traits compared to a provincial species pool. Traits over-represented in the traps included species with a high variation in colour patterns, generalist habitat preferences, extended flight periods, lower host plant specificity, and overwintering primarily as eggs. Our findings underscore the ongoing ecological filtering that favours certain species-specific traits. This study sheds light on the roles of climate change and anthropogenic impacts in shaping moth biodiversity, offers key insights into the ecological processes involved, and can guide future conservation efforts.

Century-long butterfly range expansions in northern Europe depend on climate, land use and species traits.

Climate change is an important driver of range shifts and community composition changes. Still, little is known about how the responses are influenced by the combination of land use, species interactions and species traits. We integrate climate and distributional data for 131 butterfly species in Sweden and Finland and show that cumulative species richness has increased with increasing temperature over the past 120 years. Average provincial species richness increased by 64% (range 15-229%), from 46 to 70. The rate and direction of range expansions have not matched the temperature changes, in part because colonisations have been modified by other climatic variables, land use and vary according to species characteristics representing ecological generalisation and species interactions. Results emphasise the role of a broad ecological filtering, whereby a mismatch between environmental conditions and species preferences limit the ability to disperse and establish populations in emerging climates and novel areas, with potentially widespread implications for ecosystem functioning.

Long-term warming modulates diversity, vertical structuring of microbial communities, and sulfate reduction in coastal Baltic Sea sediments.

Coastal waters such as those found in the Baltic Sea already suffer from anthropogenic related problems including increased algal blooming and hypoxia while ongoing and future climate change will likely worsen these effects. Microbial communities in sediments play a crucial role in the marine energy- and nutrient cycling, and how they are affected by climate change and shape the environment in the future is of great interest. The aims of this study were to investigate potential effects of prolonged warming on microbial community composition and nutrient cycling including sulfate reduction in surface (∼0.5 cm) to deeper sediments (∼ 24 cm). To investigate this, 16S rRNA gene amplicon sequencing was performed, and sulfate concentrations were measured and compared between sediments in a heated bay (which has been used as a cooling water outlet from a nearby nuclear power plant for approximately 50 years) and a nearby but unaffected control bay. The results showed variation in overall microbial diversity according to sediment depth and higher sulfate flux in the heated bay compared to the control bay. A difference in vertical community structure reflected increased relative abundances of sulfur oxidizing- and sulfate reducing bacteria along with a higher proportion of archaea, such as Bathyarchaeota, in the heated compared to the control bay. This was particularly evident closer to the sediment surface, indicating a compression of geochemical zones in the heated bay. These results corroborate findings in previous studies and additionally point to an amplified effect of prolonged warming deeper in the sediment, which could result in elevated concentrations of toxic compounds and greenhouse gases closer to the sediment surface.

Climate change-related warming reduces thermal sensitivity and modifies metabolic activity of coastal benthic bacterial communities.

Besides long-term average temperature increases, climate change is projected to result in a higher frequency of marine heatwaves. Coastal zones are some of the most productive and vulnerable ecosystems, with many stretches already under anthropogenic pressure. Microorganisms in coastal areas are central to marine energy and nutrient cycling and therefore, it is important to understand how climate change will alter these ecosystems. Using a long-term heated bay (warmed for 50 years) in comparison with an unaffected adjacent control bay and an experimental short-term thermal (9 days at 6-35 °C) incubation experiment, this study provides new insights into how coastal benthic water and surface sediment bacterial communities respond to temperature change. Benthic bacterial communities in the two bays reacted differently to temperature increases with productivity in the heated bay having a broader thermal tolerance compared with that in the control bay. Furthermore, the transcriptional analysis showed that the heated bay benthic bacteria had higher transcript numbers related to energy metabolism and stress compared to the control bay, while short-term elevated temperatures in the control bay incubation experiment induced a transcript response resembling that observed in the heated bay field conditions. In contrast, a reciprocal response was not observed for the heated bay community RNA transcripts exposed to lower temperatures indicating a potential tipping point in community response may have been reached. In summary, long-term warming modulates the performance, productivity, and resilience of bacterial communities in response to warming.

Associations of 16-Year Population Dynamics in Range-Expanding Moths with Temperature and Years since Establishment.

Parallel to the widespread decline of plants and animals, there is also an ongoing expansion of many species, which is especially pronounced in certain taxonomic groups and in northern latitudes. In order to inform an improved understanding of population dynamics in range-expanding taxa, we studied species richness, abundance and population growth in a sample of 25,138 individuals representing 107 range-expanding moth species at three light-trap sites in southeastern Sweden over 16 years (from 2005 to 2020) in relation to temperature and years since colonisation. Species richness and average abundance across range-expanding moths increased significantly over time, indicating a continuous influx of species expanding their ranges northward. Furthermore, average abundance and population growth increased significantly with increasing average ambient air temperature during the recording year, and average abundance also increased significantly with increasing temperature during the previous year. In general, population growth increased between years (growth rate > 1), although the population growth rate decreased significantly in association with years since colonisation. These findings highlight that, in contrast to several other studies in different parts of the world, species richness and abundance have increased in southeastern Sweden, partly because the warming climate enables range-expanding moths to realise their capacity for rapid distribution shifts and population growth. This may lead to fast and dramatic changes in community composition, with consequences for species interactions and the functioning of ecosystems. These findings are also of applied relevance for agriculture and forestry in that they can help to forecast the impacts of future invasive pest species.

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment.

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large-scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD-sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species-specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

Long-Term Warming of Baltic Sea Coastal Waters Affects Bacterial Communities in Bottom Water and Sediments Differently.

Coastal marine ecosystems are some of the most diverse natural habitats while being highly vulnerable in the face of climate change. The combination of anthropogenic influence from land and ongoing climate change will likely have severe effects on the environment, but the precise response remains uncertain. This study compared an unaffected "control" Baltic Sea bay to a "heated" bay that has undergone artificial warming from cooling water release from a nuclear power plant for ~50 years. This heated the water in a similar degree to IPCC SSP5-8.5 predictions by 2100 as natural systems to study temperature-related climate change effects. Bottom water and surface sediment bacterial communities and their biogeochemical processes were investigated to test how future coastal water warming alters microbial communities; shifts seasonal patterns, such as increased algae blooming; and influences nutrient and energy cycling, including elevated respiration rates. 16S rRNA gene amplicon sequencing and geochemical parameters demonstrated that heated bay bottom water bacterial communities were influenced by increased average temperatures across changing seasons, resulting in an overall Shannon's H diversity loss and shifts in relative abundances. In contrast, Shannon's diversity increased in the heated surface sediments. The results also suggested a trend toward smaller-sized microorganisms within the heated bay bottom waters, with a 30% increased relative abundance of small size picocyanobacteria in the summer (June). Furthermore, bacterial communities in the heated bay surface sediment displayed little seasonal variability but did show potential changes of long-term increased average temperature in the interplay with related effects on bottom waters. Finally, heated bay metabolic gene predictions from the 16S rRNA gene sequences suggested raised anaerobic processes closer to the sediment-water interface. In conclusion, climate change will likely alter microbial seasonality and diversity, leading to prolonged and increased algae blooming and elevated respiration rates within coastal waters.

Differences in phenology, daily timing of activity, and associations of temperature utilization with survival in three threatened butterflies.

We used observational data collected during a mark-recapture study that generated a total of 7503 captures of 6108 unique individuals representing three endangered butterfly species to quantify inter-and intraindividual variation in temperature utilization and examine how activity patterns vary according to season, time of day, and ambient temperature. The Marsh Fritillary, the Apollo, and the Large Blue differed in utilized temperatures and phenology. Their daily activity patterns responded differently to temperature, in part depending on whether they were active in the beginning, middle or end of the season, in part reflecting interindividual variation and intraindividual flexibility, and in part owing to differences in ecology, morphology, and colouration. Activity temperatures varied over the season, and the Apollo and the Large Blue were primarily active at the highest available ambient temperatures (on the warmest days and during the warmest part of the day). The Marsh Fritillary was active early in the season and decreased activity during the highest temperatures. The relationship between individual lifespan and the average temperature was qualitatively different in the three species pointing to species-specific selection. Lifespan increased with an increasing range of utilized temperatures in all species, possibly reflecting that intra-individual flexibility comes with a general survival benefit.

Weakened resilience of benthic microbial communities in the face of climate change.

Increased ocean temperature associated with climate change is especially intensified in coastal areas and its influence on microbial communities and biogeochemical cycling is poorly understood. In this study, we sampled a Baltic Sea bay that has undergone 50 years of warmer temperatures similar to RCP5-8.5 predictions due to cooling water release from a nuclear power plant. The system demonstrated reduced oxygen concentrations, decreased anaerobic electron acceptors, and higher rates of sulfate reduction. Chemical analyses, 16S rRNA gene amplicons, and RNA transcripts all supported sediment anaerobic reactions occurring closer to the sediment-water interface. This resulted in higher microbial diversities and raised sulfate reduction and methanogenesis transcripts, also supporting increased production of toxic sulfide and the greenhouse gas methane closer to the sediment surface, with possible release to oxygen deficient waters. RNA transcripts supported prolonged periods of cyanobacterial bloom that may result in increased climate change related coastal anoxia. Finally, while metatranscriptomics suggested increased energy production in the heated bay, a large number of stress transcripts indicated the communities had not adapted to the increased temperature and had weakened resilience. The results point to a potential feedback loop, whereby increased temperatures may amplify negative effects at the base of coastal biochemical cycling.

Climate shapes the spatiotemporal variation in color morph diversity and composition across the distribution range of Chrysomela lapponica leaf beetle.

Color polymorphism offers rich opportunities for studying the eco-evolutionary mechanisms that drive the adaptations of local populations to heterogeneous and changing environments. We explored the color morph diversity and composition in a Chrysomela lapponica leaf beetle across its entire distribution range to test the hypothesis that environmental and climatic variables shape spatiotemporal variation in the phenotypic structure of a polymorphic species. We obtained information on 13 617 specimens of this beetle from museums, private collections, and websites. These specimens (collected from 1830-2020) originated from 959 localities spanning 33° latitude, 178° longitude, and 4200 m altitude. We classified the beetles into five color morphs and searched for environmental factors that could explain the variation in the level of polymorphism (quantified by the Shannon diversity index) and in the relative frequencies of individual color morphs. The highest level of polymorphism was found at high latitudes and altitudes. The color morphs differed in their climatic requirements; composition of colour morphs was independent of the geographic distance that separated populations but changed with collection year, longitude, mean July temperature and between-year temperature fluctuations. The proportion of melanic beetles, in line with the thermal melanism hypothesis, increased with increasing latitude and altitude and decreased with increasing climate seasonality. Melanic morph frequencies also declined during the past century, but only at high latitudes and altitudes where recent climate warming was especially strong. The observed patterns suggest that color polymorphism is especially advantageous for populations inhabiting unpredictable environments, presumably due to the different climatic requirements of coexisting color morphs.

Drivers of neutral and adaptive differentiation in pike (Esox lucius) populations from contrasting environments.

Understanding how eco-evolutionary processes and environmental factors drive population differentiation and adaptation are key challenges in evolutionary biology of relevance for biodiversity protection. Differentiation requires at least partial reproductive separation, which may result from different modes of isolation such as geographic isolation (allopatry) or isolation by distance (IBD), resistance (IBR), and environment (IBE). Despite that multiple modes might jointly influence differentiation, studies that compare the relative contributions are scarce. Using RADseq, we analyse neutral and adaptive genetic diversity and structure in 11 pike (Esox lucius) populations from contrasting environments along a latitudinal gradient (54.9-63.6°N), to investigate the relative effects of IBD, IBE and IBR, and to assess whether the effects differ between neutral and adaptive variation, or across structural levels. Patterns of neutral and adaptive variation differed, probably reflecting that they have been differently affected by stochastic and deterministic processes. The importance of the different modes of isolation differed between neutral and adaptive diversity, yet were consistent across structural levels. Neutral variation was influenced by interactions among all three modes of isolation, with IBR (seascape features) playing a central role, wheares adaptive variation was mainly influenced by IBE (environmental conditions). Taken together, this and previous studies suggest that it is common that multiple modes of isolation interactively shape patterns of genetic variation, and that their relative contributions differ among systems. To enable identification of general patterns and understand how various factors influence the relative contributions, it is important that several modes are simultaneously investigated in additional populations, species and environmental settings.

Size, connectivity and edge effects of stream habitats explain spatio-temporal variation in brown trout (Salmo trutta) density.

Ecological theory postulates that the size and isolation of habitat patches impact the colonization/extinction dynamics that determine community species richness and population persistence. Given the key role of lotic habitats for life-history completion in rheophilic fish, evaluating how the distribution of swift-flowing habitats affects the abundance and dynamics of subpopulations is essential. Using extensive electrofishing data, we show that merging island biogeography with meta-population theory, where lotic habitats are considered as islands in a lentic matrix, can explain spatio-temporal variation in occurrence and density of brown trout (Salmo trutta). Subpopulations in larger and less isolated lotic habitat patches had higher average densities and smaller between-year density fluctuations. Larger lotic habitat patches also had a lower predicted risk of excessive zero-catches, indicative of lower extinction risk. Trout density further increased with distance from the edge of adjacent lentic habitats with predator (Esox lucius) presence, suggesting that edge- and matrix-related mortality contributes to the observed patterns. These results can inform the prioritization of sites for habitat restoration, dam removal and reintroduction by highlighting the role of suitable habitat size and connectivity in population abundance and stability for riverine fish populations.

Comparisons of historical Dutch commons inform about the long-term dynamics of social-ecological systems.

Human societies and natural ecosystems are under threat by growing populations, overexploitation of natural resources and climate change. This calls for more sustainable utilization of resources based on past experiences and insights from many different disciplines. Interdisciplinary approaches to studies of historical commons have potential to identify drivers of change and keys to success in the past, and offer advice about the management and use of shared resources in contemporary and future systems. We address these issues by applying an ecological perspective to historical data on social-ecological systems. We perform comparisons and time series analyses for nine successful Dutch commons for which high-resolution data on the regulatory activities and use of shared resources is available for on average 380 years (range 236 to 568) during the period 1300 to 1972. Within commons, institutional developments were oscillating, with periods of intense regulatory activity being separated by periods of low activity, and with the dynamics of regulations being largely independent across commons. Ecological theory posits that species that occupy similar niches should show correlated responses to environmental challenges; however, commons using more similar resources did not have more parallel or similar institutional developments. One notable exception was that sanctioning was more frequent in commons that directed more regulatory activities towards non-renewable subsoil resources, whereas there was no association between sanctioning and the use of renewable resources. This might indicate that commoners were aware of potential resource depletion and attempted to influence freeriding by actively trying to solve the underlying social dilemmas. Sanctioning regulations were more frequent during the first than during the second part of a common's life, indicating that while sanctioning might have been important for the establishment of commons it was not key to the long-term persistence of historical commons.

How well do genetic markers inform about responses to intraspecific admixture? A comparative analysis of microsatellites and RADseq.

BACKGROUND: Fitness consequences of intraspecific genetic admixture can vary from positive to negative depending on the genetic composition of the populations and environmental conditions. Because admixture has potential to influence the success of management and conservation efforts, genetic similarity has been suggested to be used as a proxy to predict the outcome. Studies utilizing microsatellites (a neutral marker) to investigate associations between genetic distance and admixture effects show conflicting results. Marker types that yield information on genome-wide and/or adaptive variation might be more useful for predicting responses to inter-population hybridization. In this study we utilized published data for three populations of pike (Esox lucius) to investigate associations between offspring performance (hatching success) and parental genetic similarity in experimentally purebred and admixed families, based on neutral (microsatellites), genome-wide neutral (RADseq SNPs), and adaptive (SNPs under selection) markers. RESULTS: Estimated similarity varied among the markers, likely reflecting differences in their inherent properties, but was consistently higher in purebred than admixed families. A significant interaction between marker type and admixture treatment reflected that neutral SNPs yielded higher estimates than adaptive SNPs for admixed families whereas no difference was found for purebred families, which indicates that neutral similarity was not reflective of adaptive similarity. When all samples were pooled, no association between similarity and performance was found for any marker. For microsatellites, similarity was positively correlated with hatching success in purebred families, whereas no association was found in admixed families; however, the direction of the effect differed between the population combinations. CONCLUSIONS: The results strengthen the notion that, as of today, there is no proxy that can reliably predicted the outcome of admixture. This emphasizes the need of further studies to advance knowledge that can shed light on how to safeguard against negative consequences of admixture, and thereby inform management and promote conservation of biological diversity.

Intra-population variation in reproductive timing covaries with thermal plasticity of offspring performance in perch Perca fluviatilis.

Life history theory posits that organisms should time their reproduction to coincide with environmental conditions that maximize their fitness. Population-level comparisons have contributed important insights on the adaptive value of reproductive timing and its association to environmental variation. Yet, despite its central role to ecology and evolution, the causes and consequences of variation in reproductive timing among individuals within populations are poorly understood in vertebrates other than birds. Using a combination of observational field studies and a split-brood experiment, we investigated whether differences in breeding time were associated with changes in hatching success, reproductive allocation and reaction norms linking offspring performance to temperature within an anadromous Baltic Sea population of perch Perca fluviatilis. Field observations revealed substantial variation in reproductive timing, with the breeding period lasting almost 2 months and occurring in temperatures ranging from 10 to 21℃. The hatching success of perch decreased as the reproductive season progressed. At the same time, the reproductive allocation strategy changed over the season, late breeders (the offspring of which were introduced into a high resource environment and increased predation pressure) produced more and smaller eggs that resulted in smaller larvae, compared with early breeders. The split-brood experiment in which eggs were incubated in different temperatures (10, 12, 15, 18°C) showed that differences in reproductive timing were associated with a change in the shape of the reaction norm linking offspring performance to water temperature indicative of adaptive phenotypic plasticity, with the offspring of early breeders performing best in low temperatures and the offspring of late breeders performing best in high temperatures. The seasonal changes in reproductive traits and the shape of the thermal performance suggest time-dependent adaptive differences among individuals within the population. Management actions aimed at preserving and restoring variation in the timing of reproductive events will thus likely also influence variation in associated life history traits and thermal performance curves, which could safeguard populations against environmental challenges and changes associated with exploitation and global warming.

Fish Skin Microbiomes Are Highly Variable Among Individuals and Populations but Not Within Individuals.

Fish skin-associated microbial communities are highly variable among populations and species and can impact host fitness. Still, the sources of variation in microbiome composition, and particularly how they vary among and within host individuals, have rarely been investigated. To tackle this issue, we explored patterns of variation in fish skin microbiomes across different spatial scales. We conducted replicate sampling of dorsal and ventral body sites of perch (Perca fluviatilis) from two populations and characterized the variation of fish skin-associated microbial communities with 16S rRNA gene metabarcoding. Results showed a high similarity of microbiome samples taken from the left and right side of the same fish individuals, suggesting that fish skin microbiomes can be reliably assessed and characterized even using a single sample from a specific body site. The microbiome composition of fish skin differed markedly from the bacterioplankton communities in the surrounding water and was highly variable among individuals. No ASV was present in all samples, and the most prevalent phyla, Actinobacteria, Bacteroidetes, and Proteobacteria, varied in relative abundance among fish hosts. Microbiome composition was both individual- and population specific, with most of the variation explained by individual host. At the individual level, we found no diversification in microbiome composition between dorsal and ventral body sites, but the degree of intra-individual heterogeneity varied among individuals. To identify how genetic and phenotypic characteristics of fish hosts impact the rate and nature of intra-individual temporal dynamics of the skin microbiome, and thereby contribute to the host-specific patterns documented here, remains an important task for future research.