Forest logging is associated with shape of males and gynes in the wood ant Formica aquilonia

Environmental stress can affect individual development and fitness in insects. Forest logging is a serious environmental stress for forest-specialist insects, such as the mound-building wood ant Formica aquilonia Yarrow 1955, which builds its nests into forests and is dependent on the nutrition provided by the forest habitat. We studied whether the logging causes such a strong environmental stress that it would be visible in disproportionate growth and shape of young winged wood ant gynes (‘queens’) and males. We measured head width, thorax width and wing length of 144 gynes and 100 males from 12 nests from six clear-cut areas, 11 nests from unlogged forest stands and 10 nests in clear-cut—forest margins. We observed disproportionate growth of different body parts and allometric growth of wings in both sexes in areas with different logging-induced disturbance. Gynes had larger heads in clear-cut areas and males had narrowest thoraxes in forest clear-cut edges. With an increasing thorax width, the wing length decreases steeper in clear-cuts than in other habitats in gynes and increases steeper in clear-cuts than other habitats in males. It seems that in a strongly disturbed clear-cut environment, the gynes invest the growth of head width. There seems to be a trade-off between the growth of the thorax and wings in gynes but not in males. The altered body shape may be adaptation to clearings caused by storm events of wild fires, but not forest management practices of modern era. Large-headed gynes may be better in a colony take-over, needed for the establishment of new colonies in early succession stage habitats in this temporally social parasitic species. Long-winged males may have better long-range flight ability, and they may thus have better fitness and change to disperse their genes onto new habitats.


Introduction
The current era, Anthropocene, is challenging insects, including social insects, with plenty of anthropogenic environmental changes that alter abiotic and biotic factors, such as temperature and nutrition that are known to have potential to affect growth of insects (e.g., Chown and Nicolson 2004). During the development of an insect, different body parts can grow disproportionately due to pathways that may be co-affected by environmental factors such as temperature and nutrition (e.g., Koyama et al. 2013;Lavine et al. 2015).
This raises a question whether anthropogenic disturbances in nature can affect shape of insects.
Monomorphic insects and other monomorphic animals tend to develop so that the proportional relationships are preserved between small and large individuals. This is called isometric scaling, and it is opposed by allometric scaling when there is any deviation from isometry. An excellent example of colony-level allometry is polymorphic ants with morphologically different worker castes, e.g., leafcutter ants, army ants and fire ants that have large-headed 'soldiers' (Feener et al. 1988;Tschinkel et al. 2003;Araujo and Tschinkel 2010).

3
In Finland, forest clear-cutting is the most common forest renewal method (Vaahtera 2019). For the wood ants of the Formica rufa group, removal of trees means food resource limitation as the main food sources, aphids, are reared in trees and insect prey is also mainly preyed in trees (Lenoir 2003). The decrease in food availability after clear-cutting has reduced body energy reserves of the wood ant Formica aquilonia Yarrow, 1955 (Sorvari andHakkarainen 2009;Sorvari et al. 2011). In addition, both the micro-climate and uppermost soil layers may become drier in clear-cut areas (Keenan and Kimmins 1993), which alter the temperature and humidity regulation of the mound nests of wood ants (Sorvari and Hakkarainen 2009;Sorvari et al. 2016).
Decreased food availability and malfunction in nest temperature and humidity regulation is detrimental to wood ants and is likely causing the observed nest abandonment process in clear-cuts (Sorvari and Hakkarainen 2007a). Not all nests are abandoned soon after clear-cutting; thus, the poor environment can cause physiological stress to growing ant larvae and pupae, which could affect the adult body size and size and allometric scaling of different body parts. Accordingly, a significant decrease in body size measured as head width was found in F. aquilonia workers grown in clear-cut areas (Sorvari and Hakkarainen 2009). However, in gynes (unmated queens), this effect was not found in an earlier study by the authors (Haatanen and Sorvari 2013) raising the question whether the effect of habitat degradation could be measurable in other body parts or in the allometric scaling of the body parts. This could be a sign of ants trying to adapt to the new environmental conditions. The objective of this research was to study (i) whether forest clear-cutting is associated with size of head width, thorax width and wing length, (ii) whether the body parts are developing allometrically, and (iii) whether the possible allometric scaling is associated with the strong decrease in habitat quality in gynes and males of F. aquilonia.

Materials and methods
Formica aquilonia is the most common red wood ant in northern European boreal forests (e.g., Collingwood 1979;Kvamme 1982;Punttila and Kilpeläinen 2009;Sorvari 2021). It is a typical species in sparse to medium density mature forest stands and forest edges from Scotland to Siberia (Collingwood 1979). It is a highly polygynous (multiple queens) and polydomous (multiple nests) species, with large nest mounds containing over a million workers and hundreds of queens (Pamilo 1982).
To test our hypotheses, recently produced 144 gynes and 100 males of F. aquilonia were collected in the spring of 2001 from forests (11 nests, eight forest stands), and clear-cut areas and their forest-clear-cut margins, here forth 'forest edges' (12 nests in clear-cuts and 10 nests in forest edges in six clear-cut areas), in Jyväskylä, Finland (62° 15′, 25° 46′). One to two nests were sampled from each study stand. All study stands, also clear-cuts prior to the logging, were bilberry (Vaccinium myrtillus L.) growth forests dominated by over 70-year-old Norway spruces (Picea abies L.). Logging has occurred 1-2 years prior to the sampling and the clear-cuts were treeless and planted with one-year-old Norway spruce saplings. All forest interior nests were over 50 m distance from the nearest forest edge. Clear-cut nests were all within the borders of logged area. Forest edge nests were 1-15 m away from the clear-cut margins inside the forest.
Individuals were stored in 96% alcohol until characters were measured in the year 2012. The measured characters were thorax width between wing bases, head width above compound eyes, and wing length (TW, HW, WL, respectively). The WL was measured over three wing cells from the junction of veins 1 M and 1Cu to the junction of veins 5RS and R from both left and right forewing. The mean of the left and right wing-length measurement was used as the WL. The measurement accuracy did not allow us to measure fluctuating asymmetry in WL. The measurements were made using microscope photos taken with Olympus SZ40 microscope combined with Deltapix Invenio 3S 3 Mpix digital camera using 12 × (TW), 20 × (WL), and 35 × (HW) magnification, exposure 45 ms and illuminance 2800 lx. The measurements were made using image analysis software (ImageJ) in mm scale.
We adopted similar method than in Tschinkel et al. (2003) and created size-free estimate of shape using a logarithm of the ratio of two measurements (WL/TW, HW/TW) and then regressed these log-ratios against the logarithm of other body size measurements (TW, HW). Head width (HW) and wing length (WL) were regressed as logarithms of ratios with thorax width (TW) or HW using formula log y = b + a[log (x/y)]. If a slope is different from zero by t-test, it indicates that the shape is allometric, i.e., it changes with body size (Mosiman and James 1979;Tschinkel et al. 2003).
Gynes and males of F. aquilonia are morphologically very different in shape, thus, the majority of morphological analyses between sexes is useless. Statistical analyses were performed using linear mixed model analyses for between-habitat comparisons and linear regression models for allometry in statistical software SAS version 9.4 (SAS Institute, Cary, NC, USA). All means presented are estimated marginal means of models accompanied with 95% confidence limits. Nest of origin was used as a random factor in the analyses with Kenward-Roger approximation for degrees of freedom.

Head and thorax width
In gynes, the head width (HW) was similarly associated with thorax width in all three habitat types (interaction term TW X habitat type: F 2, 134 = 1.93, P = 0.15). HW was not significantly associated with thorax width (F 1, 137 = 1.91, P = 0.17), but the HW differed between habitat types (F 2, 16.9 = 4.99, P = 0.020). The gynes from forest nests had significantly smaller HW than the clear-cut gynes and marginally nonsignificantly smaller HW than that of forest edge gynes; forest edge gynes and clear-cut gynes did not differ significantly from each other (Tukey's tests for pairwise differences, forest vs. clear-cut: P = 0.034; Forest vs. forest edge: P = 0.055, clear-cut vs. edge: P = 0.78; Fig. 1).

Allometry
The slopes of log (log 10 ) ratios of HW/TW and WL/TW against TW, and WL/HW against HW, were all negative and significantly different from zero in both gynes and males, indicating that there is allometric scaling of the shape of body parts (Table 1). The association of the slopes of log HW/TW and log WL/HW did not differ between the habitat types, but the slope log WL/TW was differently associated between the habitat types in both sexes (Table 1). Scaled with the growth of thorax width, the wing length grew less in clear-cuts that in the other habitats in gynes, whereas in  Head width (HW) and wing length (WL) were regressed as logarithms of ratios with thorax width (TW) or HW using formula log y = b + a[log (x/y)]. In both sexes, the slopes of both HW and WL differed significantly from zero indicating size-free shape changes against the size variable (TW, or HW). Both in gynes and males, the slopes WL/TW × TW were different between habitat types the males, the wing length grew less in other habitats than in clear-cuts (Fig. 4).

Discussion
Our study showed that forestry-induced habitat degradation is associated with the shape of gynes and males of F. aquilonia but with different morphological reactions. The different allometric relationships between wing length and thorax width found in gynes and males can simply be the result of differences in developmental constraints between the sexes, which likely is due to the haploid (males) and diploid (gynes) genome that already produces clear morphological differences between ant sexes. Similarly, in both sexes the different scaling relationships between wing length and thorax width among clear-cut vs. forest colonies may be the result of developmental constraints on wing growth in resource poor environment.
Previous studies show that the production of sexual offspring may decrease significantly in clear-cut areas and the produced sex ratios are more gyne-biased Hakkarainen 2005, 2007a). Thus, the developmental reactions in gynes may reflect also weaker competition between siblings in smaller broods and/or weaker competition between sexes. The ants may try in this challenging environment produce less gynes but of higher quality. An adaptive colony level strategy in clear-cuts would be producing high quality gynes and males that are able to disperse to more optimal habitats. However, based on this data long-winged gynes and males with large wing musculature (large TW) was not found more often than in other habitats. Instead, large headed gynes were common in clear-cuts, which could in some sense be a sign of quality improvement.
A mated gyne, 'a queen', of F. aquilonia has three options, (i) staying in the natal nest, (ii) disperse to an existing colony of conspecifics and (iii) try to establish a new colony (Rosengren and Pamilo 1983). The dispersal of individual gyne can be done either by flying or often by foot (Rosengren and Pamilo 1983). Formica aquilonia gyne is establishing new colonies via a risky temporal social parasitism route by taking over colonies of Serviformica species (e.g., Formica fusca) (Gösswald 1989;Buschinger 2009;Chernenko et al. 2013). While not shown in this species, big head and mandibles might be an advantage if it needs to fight. However, joining to conspecific colony is less risky than temporal parasitism strategy; thus, in an open area or forest edge area, such as clear-cuts, a mated gyne dispersing by foot likely seeks primarily an opportunity to join an existing conspecific colony.
The head width of queens was significantly larger in clear-cut areas than in forests, which contrasts our previous finding in the same study areas (Haatanen and Sorvari 2013). In addition, our current result differs from a previous result in related species, Formica podzolica, that showed no difference in gyne size between rich and poor-quality habitats (Deslippe and Savolainen 1994). The data for our previous study was collected in 2003, and current data two years before that in year 2001, when the logging disturbance was more recent. All nests in clear-cuts are not in equally poorly Fig. 4 Size-related change of wing length in gynes and males of F. aquilonia, shown in relation to thorax width in clear-cuts, forest edges (clearcut margins) and forest interiors. Predicted values from linear mixed model, clear-cut values with predicted values 95% confidence limits located, nests closer to remaining forests tend to reproduce and survive better than nests located in the centre of clearcut (Sorvari 2013). Therefore, it may be that the samples from the year 2003 were from nests closer to the forest edge. Another possibility is that the ants somehow adapted to the poor environment; however, this ant is highly forest-dwelling species (Punttila and Kilpeläinen 2009) as its populous colonies need rich food resources provided by forests and only its populous colonies in large nest mounds tend to produce gynes and males (Sorvari and Hakkarainen 2005).
Interestingly, the head size of gynes was not significantly associated with thorax width, at least with this range of sizes. Head width has often been used as an estimate of body size in ants, such as Formica sp. (e.g., Deslippe and Savolainen 1994;Sorvari and Hakkarainen 2009;Skaldina et al. 2018), but according to our results, its use as an estimate of gyne size may not be accurate. Thorax width may be another candidate for this, or dry body mass. Here, thorax width was not different in gynes originating in the three different habitats. Similarly, wing length of gynes was similar among the habitat types.
In males, head width increased with thorax width and did not differ between the habitats. This resembles the results of Deslippe and Savolainen (1994) where head width of F. podzolica males did not differ between rich and poor-quality habitats. However, thorax width was slightly smaller in forest edges compared to the other habitat types. Wing length was similar among the habitat types. For a F. aquilonia male, the dispersal ability by flight is essential, as it needs to seek a gyne to mate with. Slightly narrower thorax width, thus possibly weaker wing muscles, in forest edge habitat may affect the dispersal ability of those males. Since forest edge was here the intermediate-disturbance habitat, we have no strong suggestion for the reason behind the narrower thorax size.
Homozygous individuals tend to be more sensitive to environmental variation than heterozygous ones (e.g., Tomkins and Kotiaho 2002). Therefore, haploid, thus homozygous, males should be more sensitive to environmental change. However, in our study both the diploid gynes and haploid (homozygous) males showed morphological responses to habitat alteration.
Head, thorax and wings gets their final size and shape in metamorphosis in pupa; thus, they could be affected by (1) temperature regimes during pupation and (2) nutrition and temperature regimes during larval stage. Nutrition provides elements for different tissues and temperature is typically linked with development speed in insects (Chown and Nicolson 2004). In insects, such as our ants, the proportion of chitin tissues is large especially in wings; also head capsules and mandibles are rich with chitin whereas thorax has larger proportion of other tissues such as flight muscles (Gullan and Granston 2005). As the proportions of different tissue types differ between the studied body parts, it is possible that differences in nutrition and temperatures among habitat types can result in observed differences in body part sizes and shapes.
Both sexes showed allometric scaling of the measured body parts (HW, TW and WL). Interestingly, the allometric scaling of wings was different in clear-cuts than in the less disturbed habitat types in both sexes. Tschinkel et al. (2003) have shown that bodies of ant Solenopsis invicta workers are only slightly allometric, and usually allometry is strongest in the head. Similarly Breed (2002) showed that the allometry in workers of Paraponera clavata is strongest in head shape. In the present case, we studied only three measurements and, based on t and r 2 values (Table 1), the wing length seems to be most allometric in gynes; whereas, in males, allometry seems to be approximately the same among measured body parts. In addition, in gynes, it seems that allometry of wing size and head width is stronger in clear-cuts and its margins (Table 1). In Breed's (2002) data on P. clavata, the magnitude of allometry varied among colonies and it was thought to be due to differences in colony maturity or colony differences in developmental patterns. In our case, all colonies were mature sexual offspring producing colonies, thus, habitat quality is the most evident candidate to cause observed differences in shapes.
While some morphological response to clear-cutting, e.g., gyne head size, may sound adaptive it surely is not an adaptation to a historically recent forest management practice but more likely to storm events and wildfires in evolutionary time scale. However, there is no studies on wood ants' reactions in storm disaster and wildfire areas. Previously in this species clear-cutting has shown to decrease production of sexual offspring, especially males Hakkarainen 2005, 2007b), decreasing worker body size and body energy levels (Sorvari and Hakkarainen 2009), disturbing immunity of workers and gynes (Sorvari et al. 2008), increasing overwintering mortality (Sorvari et al. 2011), disturbing nest temperature and moisture regulation (Sorvari and Hakkarainen 2009;Sorvari et al. 2016) and finally causing abandonment of nest mounds (Sorvari and Hakkarainen 2007a).
To conclude, forest logging as a habitat degradation for a forest-dwelling ant is associated with disproportionate growth of body parts and the allometry in wings was different in logged habitat than in the less disturbed habitats. To our knowledge, this is the first observation of habitat difference-related allometric scaling in insects.
Author contributions JS collected the field data and made the statistical analyses. MKH took photographs and morphological measurements. JS wrote the first draft of the manuscript, MKH provided comments and further editions to the manuscript.
Funding Open Access funding provided by University of Turku (UTU) including Turku University Central Hospital. Emil Aaltonen foundation (JS), the Kone Foundation (JS), the Finnish Cultural Foundation (MKH), the Entomological Society of Finland (MKH) and the Turku University Foundation (MKH).

Availability of data
The data is available from the corresponding author upon request.
Code availability Not applicable.

Conflict of interest
The authors declare no conflict of interest.
Ethical approval Not applicable.

Consent to publication Not applicable.
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.