<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Ecol. Evol.</journal-id>
<journal-title>Frontiers in Ecology and Evolution</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Ecol. Evol.</abbrev-journal-title>
<issn pub-type="epub">2296-701X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fevo.2022.882357</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Ecology and Evolution</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Body Size and Symmetry Properties of Termite Soldiers Under Two Intraspecific Competition Scenarios</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Aguilera-Olivares</surname> <given-names>Daniel</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/916633/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Torres-Guzm&#x00E1;n</surname> <given-names>B&#x00E1;rbara</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Arab</surname> <given-names>Alberto</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/910174/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Flores-Prado</surname> <given-names>Luis</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1490864/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Niemeyer</surname> <given-names>Hermann M.</given-names></name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
<xref ref-type="corresp" rid="c002"><sup>&#x002A;</sup></xref>
<xref ref-type="author-notes" rid="fn002"><sup>&#x2020;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1830769/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Departamento de Silvicultura, Universidad de Concepci&#x00F3;n</institution>, <addr-line>Concepci&#x00F3;n</addr-line>, <country>Chile</country></aff>
<aff id="aff2"><sup>2</sup><institution>Departamento de Ciencias Ecol&#x00F3;gicas, Universidad de Chile</institution>, <addr-line>Santiago</addr-line>, <country>Chile</country></aff>
<aff id="aff3"><sup>3</sup><institution>Centro de Ci&#x00EA;ncias Naturais e Humanas, Universidade Federal do ABC (UFABC)</institution>, <addr-line>Santo Andr&#x00E9;</addr-line>, <country>Brazil</country></aff>
<aff id="aff4"><sup>4</sup><institution>Instituto de Entomolog&#x00ED;a, Universidad Metropolitana de Ciencias de la Educaci&#x00F3;n</institution>, <addr-line>Santiago</addr-line>, <country>Chile</country></aff>
<aff id="aff5"><sup>5</sup><institution>Departamento de Qu&#x00ED;mica, Universidad de Chile</institution>, <addr-line>Santiago</addr-line>, <country>Chile</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: David Andrew Gray, California State University, Los Angeles, United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Thomas Chouvenc, University of Florida, United States; Jo&#x00EB;l Meunier, UMR 7261 Institut de Recherche sur la Biologie de l&#x2019;Insecte (IRBI), France</p></fn>
<corresp id="c001">&#x002A;Correspondence: Daniel Aguilera-Olivares, <email>daguilera@udec.cl</email></corresp>
<corresp id="c002">Hermann M. Niemeyer, <email>niemeyer@abulafia.ciencias.uchile.cl</email></corresp>
<fn fn-type="other" id="fn002"><p><sup>&#x2020;</sup>ORCID: Daniel Aguilera-Olivares, <ext-link ext-link-type="uri" xlink:href="http://orcid.org/0000-0003-0776-0275">orcid.org/0000-0003-0776-0275</ext-link>; Alberto Arab, <ext-link ext-link-type="uri" xlink:href="http://orcid.org/0000-0003-0009-6658">orcid.org/0000-0003-0009-6658</ext-link>; Luis Flores-Prado, <ext-link ext-link-type="uri" xlink:href="http://orcid.org/0000-0002-3373-8154">orcid.org/0000-0002-3373-8154</ext-link>; Hermann M. Niemeyer, <ext-link ext-link-type="uri" xlink:href="http://orcid.org/0000-0003-2273-1608">orcid.org/0000-0003-2273-1608</ext-link></p></fn>
<fn fn-type="other" id="fn004"><p>This article was submitted to Behavioral and Evolutionary Ecology, a section of the journal Frontiers in Ecology and Evolution</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>01</day>
<month>07</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>882357</elocation-id>
<history>
<date date-type="received">
<day>23</day>
<month>02</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>18</day>
<month>05</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2022 Aguilera-Olivares, Torres-Guzm&#x00E1;n, Arab, Flores-Prado and Niemeyer.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Aguilera-Olivares, Torres-Guzm&#x00E1;n, Arab, Flores-Prado and Niemeyer</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Single-piece nesting termites live and forage in the same piece of wood throughout their life, which limit their colony size. In certain species, more than one colony thrive in a given piece of wood (multicolonial substrate) and intraspecific competition become important in this limited resource, as has been reported in <italic>Zootermopsis nevadensis</italic> (Hagen, 1858) and <italic>Neotermes chilensis</italic> (Blattodea: Kalotermitidae) (Blanchard, 1851). The effects of such competition have been described mainly at population and colony levels rather than at the individual level. In eusocial insects such as termites, intraspecific competition constitutes a stress factor imposed to a colony as a whole and should also cause developmental instability in soldiers produced under such conditions. Investment in the production of soldiers involves a trade-off between colony maintenance costs and defense benefits. Hence, we hypothesize that body size and fluctuating asymmetry, two indicators of developmental instability, will increase when two or more colonies of <italic>N. chilensis</italic> share a piece of wood (high intraspecific competition scenario). Our results showed that soldiers developing in multicolonial substrates were indeed larger and more asymmetric than soldiers developing in unicolonial substrates. The large body size in a soldier could improve its chance to win a physical contest with a non-nestmate opponent; thus, despite the high cost to produce large soldiers in small colonies, larger soldier production could be an adaptative strategy to avoid being outcompeted. However, the effects of deviations from perfect symmetry on soldier performance are not clear.</p>
</abstract>
<kwd-group>
<kwd>developmental instability</kwd>
<kwd>drywood termite</kwd>
<kwd>fluctuating asymmetry</kwd>
<kwd>Kalotermitidae</kwd>
<kwd><italic>Neotermes chilensis</italic></kwd>
<kwd>stress</kwd>
<kwd>Isoptera</kwd>
<kwd>Blattodea</kwd>
</kwd-group>
<contract-num rid="cn001">1120210</contract-num>
<contract-num rid="cn001">3190381</contract-num>
<contract-sponsor id="cn001">Fondo Nacional de Desarrollo Cient&#x00ED;fico y Tecnol&#x00F3;gico <named-content content-type="fundref-id">10.13039/501100002850</named-content></contract-sponsor>
<counts>
<fig-count count="3"/>
<table-count count="3"/>
<equation-count count="2"/>
<ref-count count="83"/>
<page-count count="10"/>
<word-count count="8309"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="intro">
<title>Introduction</title>
<p>Intraspecific competition studies in social insects have focused mainly at the colony level (<xref ref-type="bibr" rid="B59">Passera et al., 1996</xref>; <xref ref-type="bibr" rid="B75">Thomas et al., 2005</xref>; <xref ref-type="bibr" rid="B67">Sanada-Morimura et al., 2006</xref>; <xref ref-type="bibr" rid="B20">Cronin et al., 2012</xref>; <xref ref-type="bibr" rid="B47">Lloyd and Poulin, 2014</xref>; <xref ref-type="bibr" rid="B10">Blight et al., 2016</xref>) and at the population level (<xref ref-type="bibr" rid="B38">Holway et al., 1998</xref>; <xref ref-type="bibr" rid="B2">Adams and Tschinkel, 2001</xref>; <xref ref-type="bibr" rid="B44">Korb and Linsenmair, 2001</xref>; <xref ref-type="bibr" rid="B12">Boulay et al., 2007</xref>; <xref ref-type="bibr" rid="B29">Grohmann et al., 2010</xref>; <xref ref-type="bibr" rid="B13">Bourguignon et al., 2011</xref>; <xref ref-type="bibr" rid="B61">Perdereau et al., 2011</xref>; <xref ref-type="bibr" rid="B78">Tsuji, 2013</xref>; <xref ref-type="bibr" rid="B63">Pringle and Tarnita, 2017</xref>); however, the effects of intraspecific competition at the individual level are poorly studied.</p>
<p>Termites are eusocial species that live in colonies organized into three castes (<xref ref-type="bibr" rid="B76">Thorne, 1996</xref>; <xref ref-type="bibr" rid="B43">Korb and Hartfelder, 2008</xref>; <xref ref-type="bibr" rid="B23">Eggleton, 2011</xref>): (i) primary reproductives involved in reproduction, (ii) soldiers involved in defense against intruders, and (iii) true workers (or pseudergates) involved in nest keeping and taking care of juveniles and other castes and also, together with soldiers, in defense. Mandibles of termite soldiers are morphologically adapted to defense (<xref ref-type="bibr" rid="B69">Scholtz et al., 2008</xref>); hence, soldiers are unable to feed themselves and they depend on pseudergates or workers for survival (<xref ref-type="bibr" rid="B32">Haverty, 1977</xref>; <xref ref-type="bibr" rid="B35">Henderson, 1998</xref>). Consequently, investment in the production of soldiers involves a trade-off between colony maintenance costs and defense benefits (<xref ref-type="bibr" rid="B55">Noirot, 1989</xref>; <xref ref-type="bibr" rid="B15">Chouvenc et al., 2015</xref>). Pseudergates are undifferentiated and totipotent individuals who may develop into other castes such as soldiers or primary reproductives, or they may spend their whole life as pseudergates, depending on environmental conditions (<xref ref-type="bibr" rid="B56">Ogino et al., 1993</xref>; <xref ref-type="bibr" rid="B66">Roisin and Korb, 2011</xref>). Thus, environmental factors such as colony conditions (colony size, reproductive status, caste ratio, resource availability, etc.), temperature and seasonality significantly affect hormone and gene expressions (<xref ref-type="bibr" rid="B68">Scharf et al., 2007</xref>; <xref ref-type="bibr" rid="B52">Miura and Scharf, 2011</xref>) which trigger the differentiation from pseudergate to presoldier instar. This process may also be enhanced by ecological factors such as inter and/or intraspecific competition as a form of induced defense at the colony level (<xref ref-type="bibr" rid="B59">Passera et al., 1996</xref>; <xref ref-type="bibr" rid="B6">Aguilera-Olivares et al., 2017</xref>).</p>
<p>Termite species may be classified according to their nesting behavior in two main categories: separate-piece nesting termites, where foraging substrates and nesting substrate are different; and single-piece nesting termites, who spent their whole life in the same piece of wood (<xref ref-type="bibr" rid="B1">Abe, 1991</xref>; <xref ref-type="bibr" rid="B70">Shellman-Reeve, 1997</xref>). The colonies of separate-piece nesting termites are composed commonly by true workers and their size ranges from thousands to millions of individuals who can forage and move for thousands of square meters. On the other hand, colonies of single-piece nesting termites are characterized by the presence of pseudergates and colony size (usually no larger than a few thousand individuals) is limited by the volume of the piece of wood where they thrive (<xref ref-type="bibr" rid="B1">Abe, 1991</xref>; <xref ref-type="bibr" rid="B70">Shellman-Reeve, 1997</xref>; <xref ref-type="bibr" rid="B43">Korb and Hartfelder, 2008</xref>; <xref ref-type="bibr" rid="B53">Mizumoto and Bourguignon, 2021</xref>). If several colonies of single-piece nesting termites share the same substrate, inter and/or intra-specific competition is expected to occur, particularly when nesting resources become limited (<xref ref-type="bibr" rid="B77">Thorne et al., 2003</xref>; <xref ref-type="bibr" rid="B65">Ripa and Luppichini, 2004</xref>; <xref ref-type="bibr" rid="B7">Amarillo-Su&#x00E1;rez et al., 2011</xref>).</p>
<p><italic>Neotermes chilensis</italic> (Blattodea: Kalortermitidae) (Blanchard, 1851) is a single-piece nesting termite (i.e., it completes its life cycle within a finite resource) which uses dry scapes (stems of inflorescences) of the bromeliads <italic>Puya alpestris</italic> ssp. <italic>zoellneri</italic> [ex <italic>P. berteroniana</italic>&#x2014;<xref ref-type="bibr" rid="B82">Zizka et al. (2013)</xref>] and <italic>P. chilensis</italic> (Molina, 1782) among its hosts (<xref ref-type="bibr" rid="B5">Aguilera-Olivares et al., 2015</xref>). In previous work at the same study site (<xref ref-type="bibr" rid="B6">Aguilera-Olivares et al., 2017</xref>), we reported that about half of the scapes were occupied by a single colony (unicolonial substrates) and the other half contained between 2 and 9 colonies (multicolonial substrates). Additionally, we observed galleries interconnecting two colonies within a scape which had been blocked with sawdust, probably following an agonistic interaction. Under this intraspecific competition scenario, the soldiers/non-soldiers ratio significantly increased. This implies that colonies exposed to competition have an inherently higher cost due to the increased production and maintenance of soldiers, and since resources are finite in this termite, an increased developmental stress is imposed on these soldiers. <italic>N. chilensis</italic> stands as an interesting species to study the effect of intraspecific competition in an eusocial insect because: 1) it lives on a limited resource; 2) it is not able to move to another resource, so it must compete when the nesting substrate is shared with other colonies; 3) due to low colony sizes (no more than 500 individuals) it is easier to study than other termite species with typical colony sizes in the thousands to millions of individuals; and 4) about 90% of individuals inside a colony are pseudergates, i.e. totipotential individuals which could develop as soldiers, alates or remain as pseudergates their whole life according environmental signals. Moreover, the soldier caste of <italic>N. chilensis</italic> appears as particularly relevant in terms of intraspecific competition at the individual level, especially the consequences of stress during development.</p>
<p>Developmental stability is the ability of an organism to produce its developmentally programed phenotype despite epigenetic perturbations; contrastingly, developmental instability occurs when an organism is unable to buffer those perturbations (<xref ref-type="bibr" rid="B49">Markow, 1995</xref>). In species with symmetric morphological traits, the developmentally programed phenotype is one with perfect symmetry of such traits; and any random deviations from this symmetry on individuals inside a population leads to fluctuating asymmetry (<xref ref-type="bibr" rid="B58">Palmer, 1994</xref>; <xref ref-type="bibr" rid="B28">Graham et al., 2010</xref>). These deviations may be associated with genetic stress such as hybridization (<xref ref-type="bibr" rid="B27">Graham and Felley, 1985</xref>; <xref ref-type="bibr" rid="B30">Handy et al., 2004</xref>), and with environmental stress such as food limitation (<xref ref-type="bibr" rid="B74">Swaddle and Witter, 1994</xref>), heat shock (<xref ref-type="bibr" rid="B39">Hosken et al., 2000</xref>), and intraspecific competition (<xref ref-type="bibr" rid="B81">Witter and Swaddle, 1994</xref>). In insects, there is evidence that high levels of larval density constitute an environmental stressor inducing fluctuating asymmetry (<xref ref-type="bibr" rid="B18">Clarke and McKenzie, 1992</xref>; <xref ref-type="bibr" rid="B26">Gibbs and Breuker, 2006</xref>; <xref ref-type="bibr" rid="B8">Beasley et al., 2013</xref>), which may be caused by the effects of intraspecific competition due to a reduction of food availability (<xref ref-type="bibr" rid="B40">Hunt and Allen, 2000</xref>). Nevertheless, the effect on fluctuating asymmetry of another scenario of intraspecific competition caused by a restriction in the use of a nesting substrate, as occurs when two or more insect colonies share the same nest, has not been assessed.</p>
<p>On the other hand, body size is a life-history trait which has been used in insects as an indicator of developmental stress (<xref ref-type="bibr" rid="B79">Warren et al., 2006</xref>; <xref ref-type="bibr" rid="B19">Couret and Benedict, 2014</xref>). Thus, the effect of factors such as altitude (<xref ref-type="bibr" rid="B21">Cushman et al., 1993</xref>; <xref ref-type="bibr" rid="B73">Smith et al., 2007</xref>; <xref ref-type="bibr" rid="B36">Hoiss et al., 2012</xref>), temperature (<xref ref-type="bibr" rid="B11">Bochdanovits and De Jong, 2003</xref>), and competition (<xref ref-type="bibr" rid="B33">Heinrich, 1993</xref>; <xref ref-type="bibr" rid="B79">Warren et al., 2006</xref>; <xref ref-type="bibr" rid="B7">Amarillo-Su&#x00E1;rez et al., 2011</xref>; <xref ref-type="bibr" rid="B80">Wills et al., 2014</xref>; <xref ref-type="bibr" rid="B42">Korallo-Vinarskaya et al., 2015</xref>) during insect development are important predictors of body size outcome. When interference competition occurs, body size generally determines who will dominate the resources; thus, larger individuals are frequently more successful than smaller ones (<xref ref-type="bibr" rid="B34">Heinrich and Bartholomew, 1979</xref>; <xref ref-type="bibr" rid="B57">Otronen, 1988</xref>; <xref ref-type="bibr" rid="B33">Heinrich, 1993</xref>; <xref ref-type="bibr" rid="B83">Zobel and Paxton, 2007</xref>; <xref ref-type="bibr" rid="B9">Bespalova and Helms, 2014</xref>). In termites, significantly higher fluctuating asymmetry has been detected in soldiers of <italic>Coptotermes formosanus</italic> Shiraki, which develop within incipient colonies (instable and stressful scenario) in comparison with soldiers which develop within mature colonies (more stable environment) (<xref ref-type="bibr" rid="B14">Chouvenc et al., 2014a</xref>). Furthermore, studies in three <italic>Reticulitermes</italic> (<italic>R. speratus</italic>, <italic>R. virginicus</italic>, and <italic>R. flavipes</italic>) and two <italic>Cryptotermes</italic> (<italic>C. secundus</italic> and <italic>C. domesticus</italic>) termite species have shown increased female-to-male ratio because female soldiers are larger than male soldiers and thus have a better chance to win in a physical encounter (<xref ref-type="bibr" rid="B50">Matsuura, 2006</xref>; <xref ref-type="bibr" rid="B54">Muller and Korb, 2008</xref>). In summary, intraspecific competition is a stressful scenario for termite colonies which could affect body size and symmetry properties of developing individuals and hence their performance.</p>
<p>The aim of this work was to study the effects of intraspecific competition at the individual level in an eusocial insect. If competition in the context of finite resources induces colonies to invest in larger and/or more numerous soldiers, it would directly impose a developmental stress on these individuals. We hypothesized that soldiers of <italic>N. chilensis</italic> developing in colonies that share a nesting resource with other colonies (high intraspecific competition scenario) would show significantly higher values of fluctuating asymmetry and they would be larger than soldiers whose colonies occur singly within a nesting resource (low intraspecific competition scenario).</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="S2.SS1">
<title>Species and Study Area</title>
<p><italic>N. chilensis</italic> is an endemic termite from Chile distributed between 26 and 33.5&#x00B0;S (<xref ref-type="bibr" rid="B65">Ripa and Luppichini, 2004</xref>). The study area (Las Chilcas: 32&#x00B0;52&#x2032;S; 70&#x00B0;52&#x2032;W) is located within the sclerophyllous shrub community of central Chile (<xref ref-type="bibr" rid="B25">Gajardo, 1994</xref>) whose predominant species are <italic>Adesmia arborea</italic> Bert. ex Savi (Fabaceae), <italic>Colliguaya odorifera</italic> Mol. (Euphorbiaceae), <italic>Echinopsis chiloensis</italic> (Colla) H. Friedrich and G.D. Rowley (Cactaceae), <italic>Puya chilensis</italic> and <italic>P. alpestris</italic> ssp. <italic>zoellneri</italic> (Bromeliaceae). In the study area, <italic>N. chilensis</italic> builds its nests inside the dry scapes of <italic>P. alpestris</italic> ssp. <italic>zoellneri</italic>, with colonies containing up to 500 individuals and it is the only termite species present. Scapes were severed from the rest of the plant, brought to the laboratory in Santiago after being enclosed within a mesh to avoid loss of individuals, and maintained in a breeding room at 16&#x00B1;1&#x00B0;C under darkness for a maximum of 2 weeks until they were dissected.</p>
</sec>
<sec id="S2.SS2">
<title>Collection of Samples</title>
<p>A total of sixteen scapes of <italic>P. alpestris</italic> ssp. <italic>zoellneri</italic> were dissected and the number of colonies and of soldiers and non-soldiers were recorded. Eighteen colonies and fifty-three individuals in total were used in this study. Individuals collected were stored in 70% v/v alcohol. <italic>Multicolonial soldiers</italic>: five scapes had 2 or 3 colonies (mean &#x00B1; sem: 2.75 &#x00B1; 0.25 colonies per scape) some of which had soldiers (6.14 &#x00B1; 1.62 soldiers per colony). A total of seven colonies from multicolonial substrates with 29&#x2013;178 individuals (104.71 &#x00B1; 23.3 individuals per colony) were used in the study. The analysis included 22 multicolonial soldiers (2&#x2013;5 soldiers were extracted from each colony, 3.14 &#x00B1; 0.46). <italic>Unicolonial soldiers</italic>: 11 scapes had a unique colony. In these scapes, colony size ranged from 48 to 492 individuals (136.1 &#x00B1; 38.4 individuals per colony) and all had soldiers (7.91 &#x00B1; 2.07 soldiers per colony). All the 11 colonies in unicolonial substrates were used in the study. The analysis included 31 unicolonial soldiers (2&#x2013;4 soldiers were extracted from each colony, 2.82 &#x00B1; 0.18).</p>
</sec>
<sec id="S2.SS3">
<title>Measurement of Traits</title>
<p>Since termites have a non-quitinized abdomen, body size changes dramatically with hydration level. Hence, only quitinized traits were used (17 in total) to estimate body size (<xref ref-type="bibr" rid="B54">Muller and Korb, 2008</xref>; <xref ref-type="bibr" rid="B41">Johnson et al., 2011</xref>): six mandibular, seven head and four thoracic traits (described in <xref ref-type="table" rid="T1">Table 1</xref> and <xref ref-type="fig" rid="F1">Figures 1A&#x2013;E</xref>). On the other hand, the left and the right side of four traits were used to assess developmental instability through fluctuating asymmetry (<xref ref-type="table" rid="T1">Table 1</xref> and <xref ref-type="fig" rid="F1">Figure 1C</xref>). One photograph was taken of each trait for body size estimation and one photograph was taken of each side of traits selected for the fluctuating asymmetry study. An Olympus<sup>&#x00AE;</sup> Trinocular Stereo Zoom Microscope Model SZ61 was used with an integrated MSHOT 30 camera; the zoom used for each trait is given in <xref ref-type="table" rid="T1">Table 1</xref>. Using the software M-shot Digital Imaging System (<xref ref-type="bibr" rid="B51">Micro-shot Technology Co, 2010</xref>), the traits were measured in triplicate in three different days to avoid bias and reduce errors.</p>
<table-wrap position="float" id="T1">
<label>TABLE 1</label>
<caption><p>Morphological traits used to evaluate body size and fluctuating asymmetry in soldiers of <italic>N. chilensis</italic> from unicolonial and multicolonial substrates.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<tbody>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #43acc7"><bold>Morphological traits</bold></td>
<td valign="top" align="center" style="background-color: #43acc7"><bold>Abbreviation</bold></td>
<td valign="top" align="center" style="background-color: #43acc7"><bold>Zoom</bold></td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc">Right mandibular area</td>
<td valign="top" align="center" style="background-color: #ecf1dc">RMA</td>
<td valign="top" align="center" style="background-color: #ecf1dc">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc">Left mandibular area</td>
<td valign="top" align="center" style="background-color: #ecf1dc">LMA</td>
<td valign="top" align="center" style="background-color: #ecf1dc">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc">Right mandibular perimeter</td>
<td valign="top" align="center" style="background-color: #ecf1dc">RMP</td>
<td valign="top" align="center" style="background-color: #ecf1dc">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc">Left mandibular perimeter</td>
<td valign="top" align="center" style="background-color: #ecf1dc">LMP</td>
<td valign="top" align="center" style="background-color: #ecf1dc">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc">Interdental distance 1&#x2013;2 (right)</td>
<td valign="top" align="center" style="background-color: #ecf1dc">IDR</td>
<td valign="top" align="center" style="background-color: #ecf1dc">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"><bold>Mandibles</bold></td>
<td valign="top" align="left" style="background-color: #ecf1dc">Interdental distance 1&#x2013;2 (left)</td>
<td valign="top" align="center" style="background-color: #ecf1dc">IDL</td>
<td valign="top" align="center" style="background-color: #ecf1dc">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd">Head area</td>
<td valign="top" align="center" style="background-color: #d7e4bd">HA</td>
<td valign="top" align="center" style="background-color: #d7e4bd">15x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd">Head perimeter</td>
<td valign="top" align="center" style="background-color: #d7e4bd">HP</td>
<td valign="top" align="center" style="background-color: #d7e4bd">15x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd">Head length</td>
<td valign="top" align="center" style="background-color: #d7e4bd">HL</td>
<td valign="top" align="center" style="background-color: #d7e4bd">10.5x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd">Head maximum width</td>
<td valign="top" align="center" style="background-color: #d7e4bd">HMxW</td>
<td valign="top" align="center" style="background-color: #d7e4bd">15x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd">Head minimum width</td>
<td valign="top" align="center" style="background-color: #d7e4bd">HMnW</td>
<td valign="top" align="center" style="background-color: #d7e4bd">15x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd">Postmentum length</td>
<td valign="top" align="center" style="background-color: #d7e4bd">PmL</td>
<td valign="top" align="center" style="background-color: #d7e4bd">22.5x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"><bold>Head</bold></td>
<td valign="top" align="left" style="background-color: #d7e4bd">Postmentum width</td>
<td valign="top" align="center" style="background-color: #d7e4bd">PmW</td>
<td valign="top" align="center" style="background-color: #d7e4bd">22.5x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc">Pronotum length</td>
<td valign="top" align="center" style="background-color: #ecf1dc">PnL</td>
<td valign="top" align="center" style="background-color: #ecf1dc">18x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc">Pronotum width</td>
<td valign="top" align="center" style="background-color: #ecf1dc">PnW</td>
<td valign="top" align="center" style="background-color: #ecf1dc">18x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"><bold>Thorax</bold></td>
<td valign="top" align="left" style="background-color: #ecf1dc">Mesonotum length</td>
<td valign="top" align="center" style="background-color: #ecf1dc">MsL</td>
<td valign="top" align="center" style="background-color: #ecf1dc">18x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"><bold>Body size analyses</bold></td>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc">Metanotum length</td>
<td valign="top" align="center" style="background-color: #ecf1dc">MtL</td>
<td valign="top" align="center" style="background-color: #ecf1dc">18x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce">Prothoracic femur length</td>
<td valign="top" align="center" style="background-color: #c5bbce">PFL</td>
<td valign="top" align="center" style="background-color: #c5bbce">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce">Prothoracic femur width</td>
<td valign="top" align="center" style="background-color: #c5bbce">PFW</td>
<td valign="top" align="center" style="background-color: #c5bbce">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #c5bbce">Thorax</td>
<td valign="top" align="left" style="background-color: #c5bbce">Prothoracic tibiae length</td>
<td valign="top" align="center" style="background-color: #c5bbce">PTL</td>
<td valign="top" align="center" style="background-color: #c5bbce">30x</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"><bold>Fluctuating asymmetry analyses</bold></td>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce">Prothoracic femur area</td>
<td valign="top" align="center" style="background-color: #c5bbce">PFA</td>
<td valign="top" align="center" style="background-color: #c5bbce">30x</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn><p><italic>Zoom used for taking photographs to measure traits are informed.</italic></p></fn>
</table-wrap-foot>
</table-wrap>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Morphological traits used to evaluate body size and fluctuating asymmetry in soldiers of <italic>N. chilensis</italic> from unicolonial and multicolonial substrates. (A) Ventral view of the head; (B) dorsal view of the head; (C) prothoracic legs; (D) mandibles and (E) dorsal view of the thorax. Abbreviations used are described in <xref ref-type="table" rid="T1">Table 1</xref>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-10-882357-g001.tif"/>
</fig>
</sec>
<sec id="S2.SS4">
<title>Statistical Analyses</title>
<p>All the analyses and graphics were performed and generated using R version 3.4.4 (<xref ref-type="bibr" rid="B64">R Core Team, 2019</xref>).</p>
<sec id="S2.SS4.SSS1">
<title>Body Size</title>
<p>A Principal Component Analyses (PCA) was performed with the 17 traits measured using the FactoMineR package (<xref ref-type="bibr" rid="B45">Le et al., 2008</xref>) to retrieve a set of principal component variables that explained major differences in soldier morphology. The total contribution of each trait was calculated using the function <italic>fviz.contrib</italic>, the most contributed traits were selected as explained in <xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 1</xref>. The principal components that most explained variation in soldier morphology were then used as predictor variables in Generalized Mixed Models (GLMM), in order to test significant differences in body size between soldiers from unicolonial and multicolonial substrates using the glmmTMB package (<xref ref-type="bibr" rid="B48">Magnusson et al., 2017</xref>). Coloniality (unicolonial and multicolonial substrates) was the fixed effect and colony size, i.e., the total number of individuals per colony, was the covariate. In order to control the indirect effect of the colony of origin, this parameter was used as the random factor in the models. The significance of the random effect was tested in reduced models. For each model, we assessed residual distribution using the DHARMa package in R (<xref ref-type="bibr" rid="B31">Hartig, 2021</xref>).</p>
</sec>
<sec id="S2.SS4.SSS2">
<title>Fluctuating Asymmetry</title>
<p>The difference between the mean right side and the mean left side of each trait (R-L) was calculated. Following <xref ref-type="bibr" rid="B58">Palmer (1994)</xref>, the following analyses were performed: (i) normality of the (R-L) distribution was assessed using the Kolmogorov-Smirnov test with Lilliefors correction to discard significant antisymmetry; (ii) significant deviation from zero was tested using a one-sample <italic>t</italic>-test to discard significant directional asymmetry; (iii) a regression was performed between R-L and trait size to discard size dependence; and (iv) a two-way ANOVA was performed with soldiers and side (triplicate data for each side) as factors to discard measurement error.</p>
<p>In order to assess deviations from perfect symmetry, two fluctuating asymmetry indexes were chosen (<xref ref-type="bibr" rid="B58">Palmer, 1994</xref>):</p>
<list list-type="simple">
<list-item>
<label>1)</label>
<p>FA1= mean | R-L |, calculated for each four of traits in <xref ref-type="table" rid="T1">Table 1</xref>. We performed a PCA with the four FA1 index values obtained. We performed a GLMM to test significant differences in fluctuating asymmetry between soldiers from unicolonial and multicolonial substrates with PC1 and PC2 as predictor variables, similar to body size analyses described above.</p>
</list-item>
<list-item>
<label>2)</label>
<p>FA10 = &#x03C3;<sup>2</sup>, which gives an estimate of the variance between right and left sides after removing the effects of measurement errors and possible directional asymmetry. Thus, from the result of the two-way ANOVA mentioned above, the mean square of the interaction (MS<sub><italic>int</italic></sub>), mean square of the error (MS<sub><italic>error</italic></sub>), number of replicates (N<sub><italic>r</italic></sub>), number of soldiers (N<sub><italic>i</italic></sub>), and number of sides (N<sub><italic>s</italic></sub>) were used to calculate the variance (<inline-formula><mml:math id="INEQ8"><mml:msubsup><mml:mi mathvariant="normal">&#x03C3;</mml:mi><mml:mi>i</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:math></inline-formula>) for each trait (Equation 1) and their approximate degrees of freedom (Equation 2):</p>
</list-item>
</list>
<disp-formula id="S2.E1">
<label>(1)</label>
<mml:math id="M1" display="block">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="normal">&#x03C3;</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo rspace="7.5pt">=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>o</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="S2.E2">
<label>(2)</label>
<mml:math id="M2" display="block">
<mml:mrow>
<mml:mtext>df</mml:mtext>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>o</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>o</mml:mi>
<mml:mo>&#x2062;</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
</sec>
</sec>
</sec>
<sec id="S3" sec-type="results">
<title>Results</title>
<sec id="S3.SS1">
<title>Body Size Analysis</title>
<p>Our preliminary PCA showed that the variation in morphology between soldiers from unicolonial and multicolonial substrates was mainly explained by 10 traits (RMA, RMP, IDR, PnW, PnL, LMA, HP, HA, HMxW, and HL) (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 1A</xref> and <xref ref-type="supplementary-material" rid="DS1">Supplementary Table 1</xref>). These traits were selected and a second PCA was run, in which the first two principal components explained 77% of total variance of soldier morphology (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 1B</xref> and <xref ref-type="supplementary-material" rid="DS1">Supplementary Table 2</xref>). The first principal component (PC1) was named &#x201C;head scores&#x201D; because head variables (head area, head perimeter, head length, and head maximum width) were the most contributory variables (51.6%) to this component (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 1B</xref> and <xref ref-type="supplementary-material" rid="DS1">Supplementary Table 2</xref>). These variables showed strong positive correlations with this component (correlation coefficients from 0.81 to 0.90&#x2014;<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 2</xref>). High values of PC1 represent soldiers with larger heads. The second axis (PC2) was the &#x201C;mandible and head scores&#x201D; because mandible traits (right mandible perimeter, right mandible area, interdental distance 1&#x2013;2 right) and a head trait (head length) were the most contributory variables (67.11%) to this component (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 1B</xref> and <xref ref-type="supplementary-material" rid="DS1">Supplementary Table 2</xref>). These variables showed positive correlation with mandibular traits (correlation coefficients from 0.45 to 0.70 for RMP&#x2014;<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 2</xref>) and a negative correlation with the head trait (correlation coefficient &#x2212;0.44&#x2014;<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 2</xref>). The GLMM on PC1 showed that soldiers from multicolonial substrates had larger heads than soldiers from unicolonial substrates (z = &#x2212;2.59; <italic>p</italic> = 0.009) (<xref ref-type="fig" rid="F2">Figure 2</xref> and <xref ref-type="table" rid="T2">Table 2</xref>). It was also found that colony size had a significant positive effect on head traits (z = 3.49; <italic>p</italic> &#x003C; 0.002) but did not have a significant effect on the second principal component (<xref ref-type="table" rid="T2">Table 2</xref>). Similar results were obtained when the random effect was excluded from the models (<xref ref-type="table" rid="T2">Table 2</xref>). No significant deviations from model assumptions were detected (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 3A</xref>).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p>Soldiers&#x2019; head <bold>(A)</bold> and mandible/head <bold>(B)</bold> components among soldiers from multicolonial and unicolonial substrates. In the plot, the modeled means are represented by large circles, bars are the modeled 95% confidence intervals, and the model-adjusted individual response values are represented by the small dots. Metric scores were obtained from PCA of body size metrics. <italic>P</italic>-values indicate significance at <italic>p</italic> &#x003C; 0.05.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-10-882357-g002.tif"/>
</fig>
<table-wrap position="float" id="T2">
<label>TABLE 2</label>
<caption><p>Fixed effects of the GLMM analyses for body size metrics (green shadow) and fluctuating asymmetric indexes 1 (FA1) (purple shadow) of soldiers of <italic>N. chilensis</italic> from multicolonial and unicolonial substrates.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<tbody>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #43acc7"><bold>Analyses</bold></td>
<td valign="top" align="left" style="background-color: #43acc7"><bold>Estimate</bold></td>
<td valign="top" align="left" style="background-color: #43acc7"><bold>SE</bold></td>
<td valign="top" align="left" style="background-color: #43acc7"><bold><italic>z</italic>-value</bold></td>
<td valign="top" align="left" style="background-color: #43acc7"><bold>Pr (&#x003E;| z|)</bold></td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"><bold><italic>PC1 scores of body size metrics</italic></bold></td>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x00A0;&#x00A0;Coloniality (fixed effect)</td>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x2212;1.70</td>
<td valign="top" align="left" style="background-color: #ecf1dc">0.62</td>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x2212;2.72</td>
<td valign="top" align="left" style="background-color: #ecf1dc"><bold>0.007</bold></td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x00A0;&#x00A0;Colony size (covariate)</td>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x00A0;&#x00A0;0.01</td>
<td valign="top" align="left" style="background-color: #ecf1dc">0.003</td>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x00A0;&#x00A0;3.29</td>
<td valign="top" align="left" style="background-color: #ecf1dc"><bold>0.001</bold></td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"><bold><italic>PC2 scores of body size metrics</italic></bold></td>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x00A0;&#x00A0;Coloniality (fixed effect)</td>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x2212;0.36</td>
<td valign="top" align="left" style="background-color: #d7e4bd">0.40</td>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x2212;0.90</td>
<td valign="top" align="left" style="background-color: #d7e4bd">0.366</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x00A0;&#x00A0;Colony size (covariate)</td>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x2212;0.002</td>
<td valign="top" align="left" style="background-color: #d7e4bd">0.002</td>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x2212;0.94</td>
<td valign="top" align="left" style="background-color: #d7e4bd">0.349</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #e8e1ec"><bold><italic>PC1 scores of FA1</italic></bold></td>
<td valign="top" align="left" style="background-color: #e8e1ec"/>
<td valign="top" align="left" style="background-color: #e8e1ec"/>
<td valign="top" align="left" style="background-color: #e8e1ec"/>
<td valign="top" align="left" style="background-color: #e8e1ec"/>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x00A0;&#x00A0;Coloniality (fixed effect)</td>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x2212;0.92</td>
<td valign="top" align="left" style="background-color: #e8e1ec">0.32</td>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x2212;2.94</td>
<td valign="top" align="left" style="background-color: #e8e1ec"><bold>0.003</bold></td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x00A0;&#x00A0;Colony size (covariate)</td>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x2212;0.0001</td>
<td valign="top" align="left" style="background-color: #e8e1ec">0.001</td>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x00A0;&#x00A0;3.49</td>
<td valign="top" align="left" style="background-color: #e8e1ec">0.794</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #c5bbce"><bold><italic>PC2 scores of FA1</italic></bold></td>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce"/>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"><bold>Without Random Factor</bold></td>
<td valign="top" align="left" style="background-color: #c5bbce">&#x00A0;&#x00A0;Coloniality (fixed effect)</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.45</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.29</td>
<td valign="top" align="left" style="background-color: #c5bbce">1.57</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.117</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #c5bbce">&#x00A0;&#x00A0;Colony size (covariate)</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.001</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.001</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.79</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.428</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc"><bold><italic>PC1 scores of body size metrics</italic></bold></td>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
<td valign="top" align="left" style="background-color: #ecf1dc"/>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x00A0;&#x00A0;Coloniality (fixed effect)</td>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x2212;1.70</td>
<td valign="top" align="left" style="background-color: #ecf1dc">0.65</td>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x2212;2.59</td>
<td valign="top" align="left" style="background-color: #ecf1dc"><bold>0.009</bold></td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x00A0;&#x00A0;Colony size (covariate)</td>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x00A0;&#x00A0;0.01</td>
<td valign="top" align="left" style="background-color: #ecf1dc">0.003</td>
<td valign="top" align="left" style="background-color: #ecf1dc">&#x00A0;&#x00A0;3.49</td>
<td valign="top" align="left" style="background-color: #ecf1dc"><bold>0.002</bold></td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd"><bold><italic>PC2 scores of body size metrics</italic></bold></td>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
<td valign="top" align="left" style="background-color: #d7e4bd"/>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x00A0;&#x00A0;Coloniality (fixed effect)</td>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x2212;0.86</td>
<td valign="top" align="left" style="background-color: #d7e4bd">0.39</td>
<td valign="top" align="left" style="background-color: #d7e4bd">1.16</td>
<td valign="top" align="left" style="background-color: #d7e4bd">0.248</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x00A0;&#x00A0;Colony size (covariate)</td>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x2212;0.38</td>
<td valign="top" align="left" style="background-color: #d7e4bd">0.42</td>
<td valign="top" align="left" style="background-color: #d7e4bd">&#x2212;0.89</td>
<td valign="top" align="left" style="background-color: #d7e4bd">0.374</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #e8e1ec"><bold><italic>PC1 scores of FA1</italic></bold></td>
<td valign="top" align="left" style="background-color: #e8e1ec"/>
<td valign="top" align="left" style="background-color: #e8e1ec"/>
<td valign="top" align="left" style="background-color: #e8e1ec"/>
<td valign="top" align="left" style="background-color: #e8e1ec"/>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x00A0;&#x00A0;Coloniality (fixed effect)</td>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x2212;0.86</td>
<td valign="top" align="left" style="background-color: #e8e1ec">0.39</td>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x2212;2.22</td>
<td valign="top" align="left" style="background-color: #e8e1ec"><bold>0.027</bold></td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x00A0;&#x00A0;Colony size (covariate)</td>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x2212;0.0001</td>
<td valign="top" align="left" style="background-color: #e8e1ec">0.001</td>
<td valign="top" align="left" style="background-color: #e8e1ec">&#x2212;0.22</td>
<td valign="top" align="left" style="background-color: #e8e1ec">0.830</td>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #c5bbce"><bold><italic>PC2 scores of FA1</italic></bold></td>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce"/>
<td valign="top" align="left" style="background-color: #c5bbce"/>
</tr>
<tr>
<td valign="top" align="left" style="background-color: #43acc7"><bold>Random factor: Soldier nested in Colony</bold></td>
<td valign="top" align="left" style="background-color: #c5bbce">&#x00A0;&#x00A0;Coloniality (fixed effect)</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.45</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.29</td>
<td valign="top" align="left" style="background-color: #c5bbce">1.57</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.117</td>
</tr>
<tr>
<td valign="top" align="right" style="background-color: #43acc7"/>
<td valign="top" align="left" style="background-color: #c5bbce">&#x00A0;&#x00A0;Colony size (covariate)</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.001</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.001</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.79</td>
<td valign="top" align="left" style="background-color: #c5bbce">0.428</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn><p><italic>Colony size was used as covariate. All analyses were run without any random effect or with soldier nested in colony origin as random effect. Probabilities under 0.05 are highlighted in bold.</italic></p></fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="S3.SS2">
<title>Fluctuating Asymmetry Analysis</title>
<p>Preliminary tests showed non-significant antisymmetry (<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 3</xref>), non-significant directional asymmetry (<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 4</xref>), non-significant dependence of the trait on size (<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 4</xref>); and non-significant measurement errors (<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 5</xref>) in the four traits analyzed.</p>
<p>The first principal component (PC1) explained the 36.2% of variation among soldiers. Since prothoracic femur area, prothoracic femur length and prothoracic femur width FA1 index values had the highest factor scores for this principal component (<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 6</xref>), we named it &#x201C;prothoracic femur&#x201D; component (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 2</xref>). The second principal component (PC2) explained the 27.7% of variation among soldiers and prothoracic tibia length FA1 index values had the highest factor score (<xref ref-type="supplementary-material" rid="DS1">Supplementary Table 6</xref>); hence, this principal component was named &#x201C;prothoracic tibia&#x201D; (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 2</xref>). GLMM showed that soldiers from multicolonial substrates had significantly higher fluctuating asymmetry indexes than soldiers from unicolonial substrates in prothoracic femur component (z = &#x2212;2.22; <italic>p</italic> = 0.027); however, non-significant differences (z = 1.57; <italic>p</italic> = 0.117) were found in the prothoracic tibia component (<xref ref-type="table" rid="T2">Table 2</xref> and <xref ref-type="fig" rid="F3">Figure 3</xref>). Colony size had non-significant effect on both components (<xref ref-type="table" rid="T2">Table 2</xref>). Similar results were obtained when no random factors were included in the model (<xref ref-type="table" rid="T2">Table 2</xref>). No significant deviations from model assumptions were detected (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 3B</xref>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p>Soldiers&#x2019; prothoracic femur <bold>(A)</bold> and prothoracic tibia <bold>(B)</bold> principal components from FA1 indexes among soldiers from multicolonial and unicolonial substrates. In the plot, the modeled means are represented by large circles, bars are the modeled 95% confidence intervals of each mean, and the model-adjusted individual response values are represented by the small dots. Scores were obtained from PCA of fluctuating asymmetry indexes (FA1). <italic>P</italic>-values indicate significance at <italic>p</italic> &#x003C; 0.05.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-10-882357-g003.tif"/>
</fig>
<p>Finally, the FA10 index showed a significantly more fluctuating asymmetry in soldiers from multicolonial substrates than in soldiers from unicolonial substrates in all femur traits (<xref ref-type="table" rid="T3">Table 3</xref>).</p>
<table-wrap position="float" id="T3">
<label>TABLE 3</label>
<caption><p>Fluctuating asymmetry index 10 (FA10) of four morphological traits of soldiers of <italic>N. chilensis</italic> from multicolonial and unicolonial substrates.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="center" style="background-color: #43acc7">Morphological trait</td>
<td valign="top" align="center" style="background-color: #43acc7">Soldier</td>
<td valign="top" align="center" style="background-color: #43acc7">FA10 (&#x03C3;<sup>2</sup>)</td>
<td valign="top" align="center" style="background-color: #43acc7">Statistics</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="center" style="background-color: #e8e1ec">Prothoracic;</td>
<td valign="top" align="center" style="background-color: #e8e1ec">Multicolonial</td>
<td valign="top" align="center" style="background-color: #e8e1ec">0.00231</td>
<td valign="top" align="center" style="background-color: #e8e1ec"><italic>F</italic><sub>(20</sub>, <sub>27)</sub> = 2.20</td>
</tr>
<tr>
<td valign="top" align="center" style="background-color: #e8e1ec">femur length</td>
<td valign="top" align="center" style="background-color: #e8e1ec">Unicolonial</td>
<td valign="top" align="center" style="background-color: #e8e1ec">0.00105</td>
<td valign="top" align="center" style="background-color: #e8e1ec"><bold><italic>P</italic> = 0.029</bold></td>
</tr>
<tr>
<td valign="top" align="center" style="background-color: #c5bbce">Prothoracic</td>
<td valign="top" align="center" style="background-color: #c5bbce">Multicolonial</td>
<td valign="top" align="center" style="background-color: #c5bbce">0.00291</td>
<td valign="top" align="center" style="background-color: #c5bbce"><italic>F</italic><sub>(20</sub>, <sub>19)</sub> = 5.65</td>
</tr>
<tr>
<td valign="top" align="center" style="background-color: #c5bbce">femur width</td>
<td valign="top" align="center" style="background-color: #c5bbce">Unicolonial</td>
<td valign="top" align="center" style="background-color: #c5bbce">0.00052</td>
<td valign="top" align="center" style="background-color: #c5bbce"><bold><italic>P</italic> &#x003C; 0.001</bold></td>
</tr>
<tr>
<td valign="top" align="center" style="background-color: #e8e1ec">Prothoracic</td>
<td valign="top" align="center" style="background-color: #e8e1ec">Multicolonial</td>
<td valign="top" align="center" style="background-color: #e8e1ec">0.00487</td>
<td valign="top" align="center" style="background-color: #e8e1ec"><italic>F</italic><sub>(30</sub>, <sub>20)</sub> = 1.50</td>
</tr>
<tr>
<td valign="top" align="center" style="background-color: #e8e1ec">tibia length</td>
<td valign="top" align="center" style="background-color: #e8e1ec">Unicolonial</td>
<td valign="top" align="center" style="background-color: #e8e1ec">0.00731</td>
<td valign="top" align="center" style="background-color: #e8e1ec"><italic>P</italic> = 0.172</td>
</tr>
<tr>
<td valign="top" align="center" style="background-color: #c5bbce">Prothoracic</td>
<td valign="top" align="center" style="background-color: #c5bbce">Multicolonial</td>
<td valign="top" align="center" style="background-color: #c5bbce">0.01635</td>
<td valign="top" align="center" style="background-color: #c5bbce"><italic>F</italic><sub>(21</sub>, <sub>19)</sub> = 2.12</td>
</tr>
<tr>
<td valign="top" align="center" style="background-color: #c5bbce">femur area</td>
<td valign="top" align="center" style="background-color: #c5bbce">Unicolonial</td>
<td valign="top" align="center" style="background-color: #c5bbce">0.00771</td>
<td valign="top" align="center" style="background-color: #c5bbce"><bold><italic>P</italic> = 0.049</bold></td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn><p><italic>P-values under 0.05 are highlighted in bold.</italic></p></fn>
</table-wrap-foot>
</table-wrap>
</sec>
</sec>
<sec id="S4" sec-type="discussion">
<title>Discussion</title>
<p>Many abiotic and biotic factors may affect the normal ontogenetic development of insects (<xref ref-type="bibr" rid="B58">Palmer, 1994</xref>; <xref ref-type="bibr" rid="B49">Markow, 1995</xref>; <xref ref-type="bibr" rid="B28">Graham et al., 2010</xref>; <xref ref-type="bibr" rid="B19">Couret and Benedict, 2014</xref>). In this work, termite colonies were collected from a common site; hence, those present in unicolonial and multicolonial substrates were exposed to the same abiotic factors. Additionally, <italic>N. chilensis</italic> was the only termite species present in the study site and <italic>N. chilensis</italic> termites were the only insects found inside the scapes. Hence, intraspecific competition operating on a termite colony as a whole is expected to be the most important stress factor that may differentially affect unicolonial and multicolonial soldier individual development. All the indexes and analyses used in this work showed significantly higher levels of fluctuating asymmetry in soldiers developed in multicolonial substrates than soldiers developed in unicolonial substrates. Thus, our results showed that the presence of neighboring colonies generated developmental instability in soldier differentiation from pseudergates.</p>
<p>When two or more colonies share a nesting substrate, several scenarios could follow, e.g., their galleries could meet and colonies could fuse, they could fight against each other with the possible consequence of the elimination of one of them, or they could avoid each other. In exotic invasive species, nestmate recognition mechanisms generally fail because genetic diversity in the population is low and relatedness among nestmates is high. This may lead to colony fusion, as has been shown in termites (<xref ref-type="bibr" rid="B60">Perdereau et al., 2010</xref>, <xref ref-type="bibr" rid="B61">2011</xref>; <xref ref-type="bibr" rid="B46">Lee et al., 2019</xref>) and ants (<xref ref-type="bibr" rid="B38">Holway et al., 1998</xref>; <xref ref-type="bibr" rid="B24">Eyer et al., 2018</xref>). On the other hand, in native species, where genetic diversity in the population is high and relatedness between colonies is low, colony fusion is rare (<xref ref-type="bibr" rid="B71">Shelton and Grace, 1996</xref>; <xref ref-type="bibr" rid="B22">Deheer and Vargo, 2004</xref>; <xref ref-type="bibr" rid="B72">Smith et al., 2012</xref>). In the population of the native species <italic>N. chilensis</italic> under study, the mean relatedness of nestmate pairs was shown to be 0.465 &#x00B1; 0.0085 (mean &#x00B1; sem), corresponding closely to that expected for full siblings (<italic>r</italic> = 0.50) and the relatedness between non-nestmate pairs was close to zero (<xref ref-type="bibr" rid="B5">Aguilera-Olivares et al., 2015</xref>). This suggests that if two or more colonies meet within the nesting substrate, colony fusion would be rare or inexistent. In fact, nestmate recognition was reported in <italic>N. chilensis</italic>, soldiers being more aggressive toward non-nestmates than toward nestmates from different castes (<xref ref-type="bibr" rid="B3">Aguilera-Olivares et al., 2016a</xref>,<xref ref-type="bibr" rid="B4">b</xref>). These observations suggest that individuals from two colonies of <italic>N. chilensis</italic> within a scape either fight and eliminate opponent colonies or avoid them. Additionally, the colony of origin had a very slight if any effect when its inclusion or exclusion as random factor was considered in the models; thus, despite the fact that individuals inside a colony of <italic>N. chilensis</italic> are closely related, there is enough variation among soldiers to be considered independent samples, suggesting that the response to intraspecific competition in our analyses is independent of their degree of relatedness.</p>
<p>Body size is a life-history trait that has been used in insects as an indicator of developmental stress (<xref ref-type="bibr" rid="B79">Warren et al., 2006</xref>; <xref ref-type="bibr" rid="B19">Couret and Benedict, 2014</xref>). Intraspecific competition as a stress factor has been described in literature (<xref ref-type="bibr" rid="B33">Heinrich, 1993</xref>; <xref ref-type="bibr" rid="B79">Warren et al., 2006</xref>; <xref ref-type="bibr" rid="B7">Amarillo-Su&#x00E1;rez et al., 2011</xref>; <xref ref-type="bibr" rid="B80">Wills et al., 2014</xref>; <xref ref-type="bibr" rid="B42">Korallo-Vinarskaya et al., 2015</xref>). In territorial species where interference competition occurs with physical encounters displayed, larger individuals have better chances to win a fight and monopolize resources than smaller individuals (<xref ref-type="bibr" rid="B62">Price et al., 2011</xref>; <xref ref-type="bibr" rid="B9">Bespalova and Helms, 2014</xref>; <xref ref-type="bibr" rid="B37">Holland et al., 2021</xref>). In our study, we demonstrated that soldiers from multicolonial substrates had significantly larger heads than those from unicolonial substrates, strongly suggesting that in a competitive scenario, colonies produce larger soldiers which increase the chance to win fights and monopolize resources. Furthermore, <xref ref-type="bibr" rid="B54">Muller and Korb (2008)</xref> showed that females of <italic>Zootermopsis nevadensis</italic> (Hagen, 1858), a single-piece nesting termite, develop more often into soldiers than males because they are significantly larger, which is an advantage when blocking holes/galleries or fighting with other individuals (<xref ref-type="bibr" rid="B83">Zobel and Paxton, 2007</xref>). In our work, quitinized traits of <italic>N. chilensis</italic> that showed differences in size between uni- and multicolonial substrates are closely related to defensive tasks. Thus, head traits (head area, head perimeter, head length, and head maximum width) could be related to the blocking of galleries with the head and to more developed muscles to produce stronger bites (<xref ref-type="bibr" rid="B9">Bespalova and Helms, 2014</xref>). When the scape used by two or more colonies is depleted and the tunnels dug by members of different colonies finally meet, all of these traits could improve the chance to win a fight in a hypothetical scenario of interference intraspecific competition (<xref ref-type="bibr" rid="B34">Heinrich and Bartholomew, 1979</xref>; <xref ref-type="bibr" rid="B57">Otronen, 1988</xref>; <xref ref-type="bibr" rid="B33">Heinrich, 1993</xref>; <xref ref-type="bibr" rid="B83">Zobel and Paxton, 2007</xref>; <xref ref-type="bibr" rid="B9">Bespalova and Helms, 2014</xref>).</p>
<p>Soldiers totally depend on workers or pseudergates for survival (<xref ref-type="bibr" rid="B32">Haverty, 1977</xref>; <xref ref-type="bibr" rid="B35">Henderson, 1998</xref>); thus, the soldier production represent a trade-off at colony level between benefits (colony defense) and costs (investment in soldier nurturing) (<xref ref-type="bibr" rid="B55">Noirot, 1989</xref>; <xref ref-type="bibr" rid="B15">Chouvenc et al., 2015</xref>). Additionally, in incipient colonies some larvae must develops into soldiers instead of workers in order to maintain a relatively stable soldier proportion, which could have consequences in colony growth as a result of a reduction of worker number (<xref ref-type="bibr" rid="B15">Chouvenc et al., 2015</xref>). In <italic>Coptotermes</italic> species, a subterranean and separate-piece nesting termite, incipient colonies with low numbers of workers and hence a reduced nurturing capacity produce smaller and more asymmetric soldiers (i.e., low quality soldiers produced at low cost) as an strategy to keep a stable caste ratio which produces a delay on the growth of colonies (<xref ref-type="bibr" rid="B17">Chouvenc and Su, 2014b</xref>; <xref ref-type="bibr" rid="B15">Chouvenc et al., 2015</xref>, <xref ref-type="bibr" rid="B16">2017</xref>); while soldiers from mature colonies which contain high number of workers and hence exhibit high nurturing capacity, are larger and symmetric (high quality soldiers produced at high cost) (<xref ref-type="bibr" rid="B17">Chouvenc and Su, 2014b</xref>; <xref ref-type="bibr" rid="B16">Chouvenc et al., 2017</xref>). In single-piece nesting termites, such as <italic>N. chilensis</italic>, colony survival is directly related with its capacity to get and defend resources; thus, despite the high cost to produce large soldiers in small colonies, larger soldier production could be an adaptative strategy to avoid being outcompeted.</p>
</sec>
<sec id="S5" sec-type="conclusion">
<title>Conclusion</title>
<p>This study shows for first time the effect of intraspecific competition at the individual level in an eusocial insect. Soldiers of <italic>N. chilensis</italic> which developed from pseudergates in a high intraspecific competition scenario (multicolonial substrates) were larger and exhibited more asymmetric traits than soldiers which developed in a low intraspecific competition scenario (unicolonial substrates). While the large body size in a soldier could improve its chance to win a physical contest with a non-nestmate opponent, the effects of deviations from perfect symmetry on soldier performance are not clear. In single-piece nesting termites, such as <italic>N. chilensis</italic>, colony survival is directly related with its capacity to get and defend resources; thus, despite the high cost to produce large soldiers in small colonies, larger soldier production could be an adaptative strategy to avoid being outcompeted. Having shown the effects of intraspecific competition on body size and asymmetry of soldiers, it would be interesting to test the fitness consequences of larger body size and fluctuating asymmetry in soldiers of <italic>N. chilensis</italic> and its consequences for colony life. Additionally, the effects could be tested in other single-piece nesting termites such as the worldwide distributed (except in Asia) <italic>Cryptotermes brevis</italic> (Walker, 1853) and in <italic>Z. nevadensis</italic>, where multicolonial substrates have been reported.</p>
</sec>
<sec id="S6" sec-type="data-availability">
<title>Data Availability Statement</title>
<p>The original contributions presented in this study are included in the article/<xref ref-type="supplementary-material" rid="DS1">Supplementary Material</xref>, further inquiries can be directed to the corresponding author/s.</p>
</sec>
<sec id="S7">
<title>Author Contributions</title>
<p>DA-O, LF-P, BT-G, and HN: study conception and design. BT-G and DA-O: material preparation and data collection. AA, BT-G, and DA-O: statistical analyses. DA-O: writing first draft of the manuscript. All authors commented on previous versions of the manuscript and read and approved the final manuscript.</p>
</sec>
<sec id="conf1" sec-type="COI-statement">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="pudiscl1" sec-type="disclaimer">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<sec id="S8" sec-type="funding-information">
<title>Funding</title>
<p>This work was financially supported by FONDECYT grant 1120210 to HN. DA-O was grateful for FONDECYT postdoctoral project 3190381. AA was supported by the S&#x00E3;o Paulo Research Foundation (FAPESP), grant #2018/22839-6.</p>
</sec>
<ack>
<p>We thank the two reviewers for helpful comments and suggestions to improve the previous version of the manuscript.</p>
</ack>
<sec id="S10" sec-type="supplementary-material">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fevo.2022.882357/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fevo.2022.882357/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="DS1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_1.XLSX" id="TS1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Abe</surname> <given-names>T.</given-names></name></person-group> (<year>1991</year>). <article-title>Ecological factors associated with the evolution of worker and soldier castes in termites.</article-title> <source><italic>Ann. Entomol.</italic></source> <volume>9</volume> <fpage>101</fpage>&#x2013;<lpage>107</lpage>.</citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Adams</surname> <given-names>E. S.</given-names></name> <name><surname>Tschinkel</surname> <given-names>W. R.</given-names></name></person-group> (<year>2001</year>). <article-title>Mechanisms of population regulation in the fire ant Solenopsis invicta: an experimental study.</article-title> <source><italic>J. Anim. Ecol.</italic></source> <volume>70</volume> <fpage>355</fpage>&#x2013;<lpage>369</lpage>. <pub-id pub-id-type="doi">10.1046/j.1365-2656.2001.00501.x</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aguilera-Olivares</surname> <given-names>D.</given-names></name> <name><surname>Burgos-Lefimil</surname> <given-names>C.</given-names></name> <name><surname>Melendez</surname> <given-names>W.</given-names></name> <name><surname>Flores-Prado</surname> <given-names>L.</given-names></name> <name><surname>Niemeyer</surname> <given-names>H. M.</given-names></name></person-group> (<year>2016a</year>). <article-title>Chemical basis of nestmate recognition in a defense context in a one-piece nesting termite.</article-title> <source><italic>Chemoecology</italic></source> <volume>26</volume> <fpage>163</fpage>&#x2013;<lpage>172</lpage>. <pub-id pub-id-type="doi">10.1007/s00049-016-0217-x</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aguilera-Olivares</surname> <given-names>D.</given-names></name> <name><surname>Rizo</surname> <given-names>J. F.</given-names></name> <name><surname>Burgos-Lefimil</surname> <given-names>C.</given-names></name> <name><surname>Flores-Prado</surname> <given-names>L.</given-names></name> <name><surname>Niemeyer</surname> <given-names>H. M.</given-names></name></person-group> (<year>2016b</year>). <article-title>Nestmate recognition in defense against nest invasion by conspecifics during swarming in a one-piece nesting termite.</article-title> <source><italic>Rev. Chil. Hist. Nat.</italic></source> <volume>89</volume>:<issue>11</issue>. <pub-id pub-id-type="doi">10.1186/s40693-016-0063-9</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aguilera-Olivares</surname> <given-names>D.</given-names></name> <name><surname>Flores-Prado</surname> <given-names>L.</given-names></name> <name><surname>V&#x00E9;liz</surname> <given-names>D.</given-names></name> <name><surname>Niemeyer</surname> <given-names>H. M.</given-names></name></person-group> (<year>2015</year>). <article-title>Mechanisms of inbreeding avoidance in the one-piece drywood termite Neotermes chilensis.</article-title> <source><italic>Insectes Soc.</italic></source> <volume>62</volume> <fpage>237</fpage>&#x2013;<lpage>245</lpage>. <pub-id pub-id-type="doi">10.1007/s00040-015-0399-1</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aguilera-Olivares</surname> <given-names>D.</given-names></name> <name><surname>Palma-Onetto</surname> <given-names>V.</given-names></name> <name><surname>Flores-Prado</surname> <given-names>L.</given-names></name> <name><surname>Niemeyer</surname> <given-names>H. M.</given-names></name></person-group> (<year>2017</year>). <article-title>X-ray computed tomography reveals that intraspecific competition promotes soldier differentiation in a one-piece nesting termite.</article-title> <source><italic>Entomol. Exp. Appl.</italic></source> <volume>163</volume> <fpage>26</fpage>&#x2013;<lpage>34</lpage>. <pub-id pub-id-type="doi">10.1111/eea.12557</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Amarillo-Su&#x00E1;rez</surname> <given-names>A. R.</given-names></name> <name><surname>Stillwell</surname> <given-names>R. C.</given-names></name> <name><surname>Fox</surname> <given-names>C. W.</given-names></name></person-group> (<year>2011</year>). <article-title>Natural selection on body size is mediated by multiple interacting factors: a comparison of beetle populations varying naturally and experimentally in body size.</article-title> <source><italic>Ecol. Evol.</italic></source> <volume>1</volume> <fpage>1</fpage>&#x2013;<lpage>14</lpage>. <pub-id pub-id-type="doi">10.1002/ece3.1</pub-id> <pub-id pub-id-type="pmid">22393478</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Beasley</surname> <given-names>D. A. E.</given-names></name> <name><surname>Bonisoli-Alquati</surname> <given-names>A.</given-names></name> <name><surname>Mousseau</surname> <given-names>T. A.</given-names></name></person-group> (<year>2013</year>). <article-title>The use of fluctuating asymmetry as a measure of environmentally induced developmental instability: a meta-analysis.</article-title> <source><italic>Ecol. Indic.</italic></source> <volume>30</volume> <fpage>218</fpage>&#x2013;<lpage>226</lpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2013.02.024</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bespalova</surname> <given-names>I.</given-names></name> <name><surname>Helms</surname> <given-names>K. R.</given-names></name></person-group> (<year>2014</year>). <article-title>Social selection and the evolution of a female weapon in queens of the ant Messor pergandei (Hymenoptera: Formicidae).</article-title> <source><italic>Biol. J. Linn. Soc.</italic></source> <volume>113</volume> <fpage>1011</fpage>&#x2013;<lpage>1020</lpage>. <pub-id pub-id-type="doi">10.1111/bij.12376</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blight</surname> <given-names>O.</given-names></name> <name><surname>D&#x00ED;az-Mariblanca</surname> <given-names>A. G.</given-names></name> <name><surname>Cerd&#x00E1;</surname> <given-names>X.</given-names></name> <name><surname>Boulay</surname> <given-names>R.</given-names></name></person-group> (<year>2016</year>). <article-title>A proactive-reactive syndrome affects group success in an ant species.</article-title> <source><italic>Behav. Ecol.</italic></source> <volume>27</volume> <fpage>118</fpage>&#x2013;<lpage>125</lpage>. <pub-id pub-id-type="doi">10.1093/beheco/arv127</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bochdanovits</surname> <given-names>Z.</given-names></name> <name><surname>De Jong</surname> <given-names>G.</given-names></name></person-group> (<year>2003</year>). <article-title>Temperature dependent larval resource allocation shaping adult body size in Drosophila melanogaster.</article-title> <source><italic>J. Evol. Biol.</italic></source> <volume>16</volume> <fpage>1159</fpage>&#x2013;<lpage>1167</lpage>. <pub-id pub-id-type="doi">10.1046/j.1420-9101.2003.00621.x</pub-id> <pub-id pub-id-type="pmid">14640407</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Boulay</surname> <given-names>R.</given-names></name> <name><surname>Cerd&#x00E1;</surname> <given-names>X.</given-names></name> <name><surname>Simon</surname> <given-names>T.</given-names></name> <name><surname>Roldan</surname> <given-names>M.</given-names></name> <name><surname>Hefetz</surname> <given-names>A.</given-names></name></person-group> (<year>2007</year>). <article-title>Intraspecific competition in the ant Camponotus cruentatus: should we expect the &#x201C;dear enemy&#x201D; effect?</article-title> <source><italic>Anim. Behav.</italic></source> <volume>74</volume> <fpage>985</fpage>&#x2013;<lpage>993</lpage>. <pub-id pub-id-type="doi">10.1016/j.anbehav.2007.02.013</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bourguignon</surname> <given-names>T.</given-names></name> <name><surname>Leponce</surname> <given-names>M.</given-names></name> <name><surname>Roisin</surname> <given-names>Y.</given-names></name></person-group> (<year>2011</year>). <article-title>Are the spatio-temporal dynamics of soil-feeding termite colonies shaped by intra-specific competition?</article-title> <source><italic>Ecol. Entomol.</italic></source> <volume>36</volume> <fpage>776</fpage>&#x2013;<lpage>785</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-2311.2011.01328.x</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chouvenc</surname> <given-names>T.</given-names></name> <name><surname>Basille</surname> <given-names>M.</given-names></name> <name><surname>Li</surname> <given-names>H.</given-names></name> <name><surname>Su</surname> <given-names>N.-Y.</given-names></name></person-group> (<year>2014a</year>). <article-title>Developmental instability in incipient colonies of social insects.</article-title> <source><italic>PLoS One</italic></source> <volume>9</volume>:<issue>e113949</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0113949</pub-id> <pub-id pub-id-type="pmid">25423502</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chouvenc</surname> <given-names>T.</given-names></name> <name><surname>Basille</surname> <given-names>M.</given-names></name> <name><surname>Su</surname> <given-names>N.-Y.</given-names></name></person-group> (<year>2015</year>). <article-title>The production of soldiers and the maintenance of caste proportions delay the growth of termite incipient colonies.</article-title> <source><italic>Insect. Soc.</italic></source> <volume>62</volume> <fpage>23</fpage>&#x2013;<lpage>29</lpage>. <pub-id pub-id-type="doi">10.1007/s00040-014-0369-z</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chouvenc</surname> <given-names>T.</given-names></name> <name><surname>Basille</surname> <given-names>M.</given-names></name> <name><surname>Su</surname> <given-names>N.-Y.</given-names></name></person-group> (<year>2017</year>). <article-title>Role of accelerated developmental pathway and limited nurturing capacity on soldier developmental instability in subterranean termite incipient colonies.</article-title> <source><italic>Insect. Soc.</italic></source> <volume>64</volume> <fpage>477</fpage>&#x2013;<lpage>483</lpage>. <pub-id pub-id-type="doi">10.1007/s00040-017-0566-7</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chouvenc</surname> <given-names>T.</given-names></name> <name><surname>Su</surname> <given-names>N.-Y.</given-names></name></person-group> (<year>2014b</year>). <article-title>Colony age-dependent pathway in caste development of Coptotermes formosanus Shiraki.</article-title> <source><italic>Insect. Soc.</italic></source> <volume>61</volume> <fpage>171</fpage>&#x2013;<lpage>182</lpage>. <pub-id pub-id-type="doi">10.1007/s00040-014-0343-9</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Clarke</surname> <given-names>G. M.</given-names></name> <name><surname>McKenzie</surname> <given-names>L. J.</given-names></name></person-group> (<year>1992</year>). <article-title>Fluctuating asymmetry as a quality control indicator for insect mass rearing processes.</article-title> <source><italic>J. Econ. Entomol.</italic></source> <volume>85</volume> <fpage>2045</fpage>&#x2013;<lpage>2050</lpage>. <pub-id pub-id-type="doi">10.1093/jee/85.6.2045</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Couret</surname> <given-names>J.</given-names></name> <name><surname>Benedict</surname> <given-names>M. Q.</given-names></name></person-group> (<year>2014</year>). <article-title>A meta-analysis of the factors influencing development rate variation in Aedes aegypti (Diptera: Culicidae).</article-title> <source><italic>BMC Ecol.</italic></source> <volume>14</volume>:<issue>3</issue>. <pub-id pub-id-type="doi">10.1186/1472-6785-14-3</pub-id> <pub-id pub-id-type="pmid">24495345</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cronin</surname> <given-names>A. L.</given-names></name> <name><surname>F&#x00E9;d&#x00E9;rici</surname> <given-names>P.</given-names></name> <name><surname>Doums</surname> <given-names>C.</given-names></name> <name><surname>Monnin</surname> <given-names>T.</given-names></name></person-group> (<year>2012</year>). <article-title>The influence of intraspecific competition on resource allocation during dependent colony foundation in a social insect.</article-title> <source><italic>Oecologia</italic></source> <volume>168</volume> <fpage>361</fpage>&#x2013;<lpage>369</lpage>. <pub-id pub-id-type="doi">10.1007/s00442-011-2098-6</pub-id> <pub-id pub-id-type="pmid">21833638</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cushman</surname> <given-names>J. H.</given-names></name> <name><surname>Lawton</surname> <given-names>J. H.</given-names></name> <name><surname>Manly</surname> <given-names>B. F. J.</given-names></name></person-group> (<year>1993</year>). <article-title>Latitudinal patterns in European ant assemblages: variation in species richness and body size.</article-title> <source><italic>Oecologia</italic></source> <volume>95</volume> <fpage>30</fpage>&#x2013;<lpage>37</lpage>. <pub-id pub-id-type="doi">10.1007/BF00649503</pub-id> <pub-id pub-id-type="pmid">28313308</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Deheer</surname> <given-names>C. J.</given-names></name> <name><surname>Vargo</surname> <given-names>E. L.</given-names></name></person-group> (<year>2004</year>). <article-title>Colony genetic organization and colony fusion in the termite Reticulitermes flavipes as revealed by foraging patterns over time and space.</article-title> <source><italic>Mol. Ecol.</italic></source> <volume>13</volume> <fpage>431</fpage>&#x2013;<lpage>441</lpage>. <pub-id pub-id-type="doi">10.1046/j.1365-294X.2003.2065.x</pub-id> <pub-id pub-id-type="pmid">14717897</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Eggleton</surname> <given-names>P.</given-names></name></person-group> (<year>2011</year>). &#x201C;<article-title>An introduction to termites: biology, taxonomy and functional morphology</article-title>,&#x201D; in <source><italic>Biology of Termites: a Modern Synthesis</italic></source>, <role>eds</role> <person-group person-group-type="editor"><name><surname>Bignell</surname> <given-names>D. E.</given-names></name> <name><surname>Roisin</surname> <given-names>Y.</given-names></name> <name><surname>Lo</surname> <given-names>N.</given-names></name></person-group> (<publisher-loc>Netherlands</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>1</fpage>&#x2013;<lpage>26</lpage>.</citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Eyer</surname> <given-names>P. A.</given-names></name> <name><surname>McDowell</surname> <given-names>B.</given-names></name> <name><surname>Johnson</surname> <given-names>L. N. L.</given-names></name> <name><surname>Calcaterra</surname> <given-names>L. A.</given-names></name> <name><surname>Fernandez</surname> <given-names>M. B.</given-names></name> <name><surname>Shoemaker</surname> <given-names>D.</given-names></name><etal/></person-group> (<year>2018</year>). <article-title>Supercolonial structure of invasive populations of the tawny crazy ant Nylanderia fulva in the US.</article-title> <source><italic>BMC Evol. Biol.</italic></source> <volume>18</volume>:<issue>209</issue>. <pub-id pub-id-type="doi">10.1186/s12862-018-1336-5</pub-id> <pub-id pub-id-type="pmid">30594137</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gajardo</surname> <given-names>R.</given-names></name></person-group> (<year>1994</year>). <source><italic>La Vegetaci&#x00F3;n Natural de Chile. Clasificaci&#x00F3;n y Distribuci&#x00F3;n.</italic></source> <publisher-loc>Santiago</publisher-loc>: <publisher-name>Editorial Universitaria</publisher-name>.</citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gibbs</surname> <given-names>M.</given-names></name> <name><surname>Breuker</surname> <given-names>C.</given-names></name></person-group> (<year>2006</year>). <article-title>Effect of larval-rearing density on adult life-history traits and developmental stability of the dorsal eyespot pattern in the speckled wood butterfly, Pararge aegeria.</article-title> <source><italic>Entomol. Exp. Appl.</italic></source> <volume>118</volume> <fpage>41</fpage>&#x2013;<lpage>47</lpage>. <pub-id pub-id-type="doi">10.1111/j.1570-7458.2006.00361.x</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Graham</surname> <given-names>J. H.</given-names></name> <name><surname>Felley</surname> <given-names>J. D.</given-names></name></person-group> (<year>1985</year>). <article-title>Genomic coadaptation and developmental stability within introgressed populations of Enneacanthus gloriosus and E. obesus (Pisces. Centrarchidae).</article-title> <source><italic>Evolution</italic></source> <volume>39</volume>:<issue>104</issue>. <pub-id pub-id-type="doi">10.2307/2408520</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Graham</surname> <given-names>J. H.</given-names></name> <name><surname>Raz</surname> <given-names>S.</given-names></name> <name><surname>Hel-Or</surname> <given-names>H.</given-names></name> <name><surname>Nevo</surname> <given-names>E.</given-names></name></person-group> (<year>2010</year>). <article-title>Fluctuating asymmetry: methods, theory, and applications.</article-title> <source><italic>Symmetry</italic></source> <volume>2</volume> <fpage>466</fpage>&#x2013;<lpage>540</lpage>. <pub-id pub-id-type="doi">10.3390/sym2020466</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Grohmann</surname> <given-names>C.</given-names></name> <name><surname>Oldeland</surname> <given-names>J.</given-names></name> <name><surname>Stoyan</surname> <given-names>D.</given-names></name> <name><surname>Linsenmair</surname> <given-names>K. E.</given-names></name></person-group> (<year>2010</year>). <article-title>Multi-scale pattern analysis of a mound-building termite species.</article-title> <source><italic>Insect. Soc.</italic></source> <volume>57</volume> <fpage>477</fpage>&#x2013;<lpage>486</lpage>. <pub-id pub-id-type="doi">10.1007/s00040010-0107-0</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Handy</surname> <given-names>S. M.</given-names></name> <name><surname>McBreen</surname> <given-names>K.</given-names></name> <name><surname>Cruzan</surname> <given-names>M. B.</given-names></name></person-group> (<year>2004</year>). <article-title>Patterns of fitness and fluctuating asymmetry across a broad hybrid zone.</article-title> <source><italic>Int. J. Plant Sci.</italic></source> <volume>165</volume> <fpage>973</fpage>&#x2013;<lpage>981</lpage>. <pub-id pub-id-type="doi">10.1086/423872</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hartig</surname> <given-names>F.</given-names></name></person-group> (<year>2021</year>). <article-title>DHARMa | Diagnostics for HierArchical Regession Models. R Package Version 0.4.4</article-title>. <ext-link ext-link-type="uri" xlink:href="http://florianhartig.github.io/DHARMa/">http://florianhartig.github.io/DHARMa/</ext-link></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Haverty</surname> <given-names>M. I.</given-names></name></person-group> (<year>1977</year>). <article-title>The proportion of soldiers in termite colonies: a list and a bibliography (Isoptera).</article-title> <source><italic>Sociobiology</italic></source> <volume>2</volume> <fpage>199</fpage>&#x2013;<lpage>216</lpage>.</citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Heinrich</surname> <given-names>B.</given-names></name></person-group> (<year>1993</year>). &#x201C;<article-title>Beetles large and small</article-title>&#x201D;, in <source><italic>The Hot-Blooded Insects</italic></source>, (<publisher-loc>Heidelberg</publisher-loc>: <publisher-name>Springer Berlin Heidelberg</publisher-name>), <fpage>191</fpage>&#x2013;<lpage>226</lpage>.</citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Heinrich</surname> <given-names>B.</given-names></name> <name><surname>Bartholomew</surname> <given-names>G. A.</given-names></name></person-group> (<year>1979</year>). <article-title>Roles of endothermy and size in inter- and intraspecific competition for elephant dung in an african dung beetle, Scarabaeus laevistriatus.</article-title> <source><italic>Physiol. Zool.</italic></source> <volume>52</volume> <fpage>484</fpage>&#x2013;<lpage>496</lpage>. <pub-id pub-id-type="doi">10.1086/physzool.52.4.30155939</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Henderson</surname> <given-names>G.</given-names></name></person-group> (<year>1998</year>). &#x201C;<article-title>Primer pheromones and possible soldier caste influence on the evolution of sociality in lower termites</article-title>,&#x201D; in <source><italic>Pheromone Communication in Social Insects: Ants, Wasps, Bees, and Termites</italic></source>, <role>eds</role> <person-group person-group-type="editor"><name><surname>Vander Meer</surname> <given-names>R. K.</given-names></name> <name><surname>Breed</surname> <given-names>M. D.</given-names></name> <name><surname>Winston</surname> <given-names>M. L.</given-names></name> <name><surname>Espelie</surname> <given-names>K. E.</given-names></name></person-group> (<publisher-loc>Boulder, CO</publisher-loc>: <publisher-name>Westview Press</publisher-name>), <fpage>314</fpage>&#x2013;<lpage>330</lpage>.</citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hoiss</surname> <given-names>B.</given-names></name> <name><surname>Krauss</surname> <given-names>J.</given-names></name> <name><surname>Potts</surname> <given-names>S. G.</given-names></name> <name><surname>Roberts</surname> <given-names>S.</given-names></name> <name><surname>Steffan-Dewenter</surname> <given-names>I.</given-names></name></person-group> (<year>2012</year>). <article-title>Altitude acts as an environmental filter on phylogenetic composition, traits and diversity in bee communities.</article-title> <source><italic>Proc. R Soc. B Biol. Sci.</italic></source> <volume>279</volume> <fpage>4447</fpage>&#x2013;<lpage>4456</lpage>. <pub-id pub-id-type="doi">10.1098/rspb.2012.1581</pub-id> <pub-id pub-id-type="pmid">22933374</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Holland</surname> <given-names>J. B.</given-names></name> <name><surname>Nakayama</surname> <given-names>S.</given-names></name> <name><surname>Porfiri</surname> <given-names>M.</given-names></name> <name><surname>Nov</surname> <given-names>O.</given-names></name> <name><surname>Bloch</surname> <given-names>G.</given-names></name></person-group> (<year>2021</year>). <article-title>Body size and behavioural plasticity interact to influence the performance of free-foraging bumble bee colonies.</article-title> <source><italic>Insects</italic></source> <volume>12</volume>:<issue>236</issue>. <pub-id pub-id-type="doi">10.3390/insects12030236</pub-id> <pub-id pub-id-type="pmid">33802199</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Holway</surname> <given-names>D. A.</given-names></name> <name><surname>Suarez</surname> <given-names>A. V.</given-names></name> <name><surname>Case</surname> <given-names>T. J.</given-names></name></person-group> (<year>1998</year>). <article-title>Loss of intraspecific aggression in the success of a widespread invasive social insect.</article-title> <source><italic>Science</italic></source> <volume>282</volume> <fpage>949</fpage>&#x2013;<lpage>952</lpage>. <pub-id pub-id-type="doi">10.1126/science.282.5390.949</pub-id> <pub-id pub-id-type="pmid">9794767</pub-id></citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hosken</surname> <given-names>D. J.</given-names></name> <name><surname>Blanckenhorn</surname> <given-names>W. U.</given-names></name> <name><surname>Ward</surname> <given-names>P. I.</given-names></name></person-group> (<year>2000</year>). <article-title>Developmental stability in yellow dung flies (Scathophaga stercoraria): fluctuating asymmetry, heterozygosity and environmental stress.</article-title> <source><italic>J. Evol. Biol.</italic></source> <volume>13</volume> <fpage>919</fpage>&#x2013;<lpage>926</lpage>. <pub-id pub-id-type="doi">10.1046/j.1420-9101.2000.00239.x</pub-id></citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hunt</surname> <given-names>J.</given-names></name> <name><surname>Allen</surname> <given-names>G. R.</given-names></name></person-group> (<year>2000</year>). <article-title>Larval density and developmental instability in the acoustically orienting parasitoid Homotrixa alleni.</article-title> <source><italic>Acta Ethol.</italic></source> <volume>2</volume> <fpage>129</fpage>&#x2013;<lpage>132</lpage>. <pub-id pub-id-type="doi">10.1007/s102110000018</pub-id></citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Johnson</surname> <given-names>S. E.</given-names></name> <name><surname>Breisch</surname> <given-names>N. L.</given-names></name> <name><surname>Momen</surname> <given-names>B.</given-names></name> <name><surname>Thorne</surname> <given-names>B. L.</given-names></name></person-group> (<year>2011</year>). <article-title>Morphology and gonad development of normal soldiers and reproductive soldiers of the termite Zootermopsis nevadensis nevadensis (Isoptera, Archotermopsidae).</article-title> <source><italic>Zookeys</italic></source> <volume>148</volume> <fpage>15</fpage>&#x2013;<lpage>30</lpage>. <pub-id pub-id-type="doi">10.3897/zookeys.148.1672</pub-id> <pub-id pub-id-type="pmid">22287888</pub-id></citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Korallo-Vinarskaya</surname> <given-names>N. P.</given-names></name> <name><surname>Vinarski</surname> <given-names>M. V.</given-names></name> <name><surname>Khokhlova</surname> <given-names>I. S.</given-names></name> <name><surname>Shenbrot</surname> <given-names>G. I.</given-names></name> <name><surname>Krasnov</surname> <given-names>B. R.</given-names></name></person-group> (<year>2015</year>). <article-title>Intraspecific variation of body size in a gamasid mite Laelaps clethrionomydis: environment, geography and host dependence.</article-title> <source><italic>Parasitol. Res.</italic></source> <volume>114</volume> <fpage>3767</fpage>&#x2013;<lpage>3774</lpage>. <pub-id pub-id-type="doi">10.1007/s00436-015-4606-9</pub-id> <pub-id pub-id-type="pmid">26141434</pub-id></citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Korb</surname> <given-names>J.</given-names></name> <name><surname>Hartfelder</surname> <given-names>K.</given-names></name></person-group> (<year>2008</year>). <article-title>Life history and development - A framework for understanding developmental plasticity in lower termites.</article-title> <source><italic>Biol. Rev.</italic></source> <volume>83</volume> <fpage>295</fpage>&#x2013;<lpage>313</lpage>. <pub-id pub-id-type="doi">10.1111/j.1469-185X.2008.00044.x</pub-id> <pub-id pub-id-type="pmid">18979593</pub-id></citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Korb</surname> <given-names>J.</given-names></name> <name><surname>Linsenmair</surname> <given-names>K. E.</given-names></name></person-group> (<year>2001</year>). <article-title>The causes of spatial patterning of mounds of a fungus-cultivating termite: results from nearest-neighbour analysis and ecological studies.</article-title> <source><italic>Oecologia</italic></source> <volume>127</volume> <fpage>324</fpage>&#x2013;<lpage>333</lpage>. <pub-id pub-id-type="doi">10.1007/s004420000597</pub-id> <pub-id pub-id-type="pmid">28547103</pub-id></citation></ref>
<ref id="B45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Le</surname> <given-names>S.</given-names></name> <name><surname>Josse</surname> <given-names>J.</given-names></name> <name><surname>Husson</surname> <given-names>F.</given-names></name></person-group> (<year>2008</year>). <article-title>FactoMineR: An R Package for Multivariate Analysis.</article-title> <source><italic>J. Stat. Softw.</italic></source> <volume>25</volume> <fpage>1</fpage>&#x2013;<lpage>18</lpage>. <pub-id pub-id-type="doi">10.18637/jss.v025.i01</pub-id></citation></ref>
<ref id="B46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lee</surname> <given-names>S.-B.</given-names></name> <name><surname>Mullins</surname> <given-names>A.</given-names></name> <name><surname>Aguilera-Olivares</surname> <given-names>D.</given-names></name> <name><surname>Chouvenc</surname> <given-names>T.</given-names></name> <name><surname>Su</surname> <given-names>N.-Y.</given-names></name></person-group> (<year>2019</year>). <article-title>Fused colonies of the Formosan Subterranean Termite (Blattodea: Rhinotermitidae) for laboratory experiments.</article-title> <source><italic>J. Econ. Entomol.</italic></source> <volume>112</volume> <fpage>2311</fpage>&#x2013;<lpage>2315</lpage>. <pub-id pub-id-type="doi">10.1093/jee/toz154</pub-id> <pub-id pub-id-type="pmid">31165146</pub-id></citation></ref>
<ref id="B47"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lloyd</surname> <given-names>M. M.</given-names></name> <name><surname>Poulin</surname> <given-names>R.</given-names></name></person-group> (<year>2014</year>). <article-title>Geographic variation in caste ratio of trematode colonies with a division of labour reflect local adaptation.</article-title> <source><italic>Parasitol. Res.</italic></source> <volume>113</volume> <fpage>2593</fpage>&#x2013;<lpage>2602</lpage>. <pub-id pub-id-type="doi">10.1007/s00436-014-3913-x</pub-id> <pub-id pub-id-type="pmid">24770673</pub-id></citation></ref>
<ref id="B48"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Magnusson</surname> <given-names>A.</given-names></name> <name><surname>Skaug</surname> <given-names>H.</given-names></name> <name><surname>Nielsen</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2017</year>). <source><italic>Package &#x2018;glmmTMB.&#x2019; R Package Version 0.2. 0.</italic></source></citation></ref>
<ref id="B49"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Markow</surname> <given-names>T. A.</given-names></name></person-group> (<year>1995</year>). <article-title>Evolutionary ecology and developmental instability.</article-title> <source><italic>Annu. Rev. Entomol.</italic></source> <volume>40</volume> <fpage>105</fpage>&#x2013;<lpage>120</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.en.40.010195.000541</pub-id></citation></ref>
<ref id="B50"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Matsuura</surname> <given-names>K.</given-names></name></person-group> (<year>2006</year>). <article-title>A novel hypothesis for the origin of the sexual division of labor in termites: which sex should be soldiers?</article-title> <source><italic>Evol. Ecol.</italic></source> <volume>20</volume> <fpage>565</fpage>&#x2013;<lpage>574</lpage>. <pub-id pub-id-type="doi">10.1007/s10682-006-9117-9</pub-id></citation></ref>
<ref id="B51"><citation citation-type="journal"><collab>Micro-shot Technology Co</collab> (<year>2010</year>). <source><italic>M-shot Digital Imaging System.</italic></source> <publisher-loc>Guangzhou</publisher-loc>: <publisher-name>Micro-shot Technology Co., Ltd</publisher-name>.</citation></ref>
<ref id="B52"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Miura</surname> <given-names>T.</given-names></name> <name><surname>Scharf</surname> <given-names>M. E.</given-names></name></person-group> (<year>2011</year>). &#x201C;<article-title>Molecular basis underlying caste differentiation in termites</article-title>,&#x201D; in <source><italic>Biology of termites: a modern synthesis</italic></source>, <role>eds</role> <person-group person-group-type="editor"><name><surname>Bignell</surname> <given-names>D. E.</given-names></name> <name><surname>Roisin</surname> <given-names>Y.</given-names></name> <name><surname>Lo</surname> <given-names>N.</given-names></name></person-group> (<publisher-loc>Netherlands</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>211</fpage>&#x2013;<lpage>253</lpage>.</citation></ref>
<ref id="B53"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mizumoto</surname> <given-names>N.</given-names></name> <name><surname>Bourguignon</surname> <given-names>T.</given-names></name></person-group> (<year>2021</year>). <article-title>The evolution of body size in termites.</article-title> <source><italic>Proc. R Soc. B</italic></source> <volume>288</volume>:<issue>20211458</issue>. <pub-id pub-id-type="doi">10.1098/rspb.2021.1458</pub-id> <pub-id pub-id-type="pmid">34784763</pub-id></citation></ref>
<ref id="B54"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Muller</surname> <given-names>H.</given-names></name> <name><surname>Korb</surname> <given-names>J.</given-names></name></person-group> (<year>2008</year>). <article-title>Male or female soldiers? An evaluation of several factors which may influence soldier sex ratio in lower termites.</article-title> <source><italic>Insectes Soc.</italic></source> <volume>55</volume> <fpage>213</fpage>&#x2013;<lpage>219</lpage>. <pub-id pub-id-type="doi">10.1007/s00040-008-0996-3</pub-id></citation></ref>
<ref id="B55"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Noirot</surname> <given-names>C. H.</given-names></name></person-group> (<year>1989</year>). <article-title>Social structure in termite societies.</article-title> <source><italic>Ethol. Ecol. Evol.</italic></source> <volume>1</volume> <fpage>1</fpage>&#x2013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1080/08927014.1989.9525528</pub-id></citation></ref>
<ref id="B56"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ogino</surname> <given-names>K.</given-names></name> <name><surname>Hirono</surname> <given-names>Y.</given-names></name> <name><surname>Matsumoto</surname> <given-names>T.</given-names></name> <name><surname>Ishikawa</surname> <given-names>H.</given-names></name></person-group> (<year>1993</year>). <article-title>Juvenile hormone analogue, S-31183, causes a high level induction of presoldier differentiation in the japanese damp-wood termite.</article-title> <source><italic>Zoolog. Sci.</italic></source> <volume>10</volume> <fpage>361</fpage>&#x2013;<lpage>366</lpage>.</citation></ref>
<ref id="B57"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Otronen</surname> <given-names>M.</given-names></name></person-group> (<year>1988</year>). <article-title>The effect of body size on the outcome of fights in burying beetles (Nicrophorus).</article-title> <source><italic>Ann. Zool. Fennici</italic></source> <volume>25</volume> <fpage>191</fpage>&#x2013;<lpage>201</lpage>. <pub-id pub-id-type="doi">10.1086/687392</pub-id> <pub-id pub-id-type="pmid">27501089</pub-id></citation></ref>
<ref id="B58"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Palmer</surname> <given-names>A. R.</given-names></name></person-group> (<year>1994</year>). &#x201C;<article-title>Fluctuating asymmetry analyses: a primer</article-title>,&#x201D; in <source><italic>Developmental Instability: Its Origins and Evolutionary Implications</italic></source>, <role>ed.</role> <person-group person-group-type="editor"><name><surname>Markow</surname> <given-names>T. A.</given-names></name></person-group> (<publisher-loc>Dordrecht</publisher-loc>: <publisher-name>Kluwer Academic Publishers</publisher-name>), <fpage>335</fpage>&#x2013;<lpage>364</lpage>.</citation></ref>
<ref id="B59"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Passera</surname> <given-names>L.</given-names></name> <name><surname>Roncin</surname> <given-names>E.</given-names></name> <name><surname>Kaufmann</surname> <given-names>B.</given-names></name> <name><surname>Keller</surname> <given-names>L.</given-names></name></person-group> (<year>1996</year>). <article-title>Increased soldier production in ant colonies exposed to intraspecific competition.</article-title> <source><italic>Nature</italic></source> <volume>379</volume> <fpage>630</fpage>&#x2013;<lpage>631</lpage>. <pub-id pub-id-type="doi">10.1038/379630a0</pub-id></citation></ref>
<ref id="B60"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Perdereau</surname> <given-names>E.</given-names></name> <name><surname>Bagn&#x00E8;res</surname> <given-names>A.-G.</given-names></name> <name><surname>Dupont</surname> <given-names>S.</given-names></name> <name><surname>Dedeine</surname> <given-names>F.</given-names></name></person-group> (<year>2010</year>). <article-title>High occurrence of colony fusion in a European population of the American termite Reticulitermes flavipes.</article-title> <source><italic>Insectes Soc.</italic></source> <volume>57</volume> <fpage>393</fpage>&#x2013;<lpage>402</lpage>. <pub-id pub-id-type="doi">10.1007/s00040-010-0096-z</pub-id></citation></ref>
<ref id="B61"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Perdereau</surname> <given-names>E.</given-names></name> <name><surname>Dedeine</surname> <given-names>F.</given-names></name> <name><surname>Christid&#x00E8;s</surname> <given-names>J. P.</given-names></name> <name><surname>Dupont</surname> <given-names>S.</given-names></name> <name><surname>Bagn&#x00E8;res</surname> <given-names>A.-G.</given-names></name></person-group> (<year>2011</year>). <article-title>Competition between invasive and indigenous species: an insular case study of subterranean termites.</article-title> <source><italic>Biol. Invasions</italic></source> <volume>13</volume> <fpage>1457</fpage>&#x2013;<lpage>1470</lpage>. <pub-id pub-id-type="doi">10.1007/s10530-010-9906-5</pub-id></citation></ref>
<ref id="B62"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Price</surname> <given-names>P. W.</given-names></name> <name><surname>Denno</surname> <given-names>R. F.</given-names></name> <name><surname>Eubanks</surname> <given-names>M. D.</given-names></name> <name><surname>Finke</surname> <given-names>D. L.</given-names></name> <name><surname>Kaplan</surname> <given-names>I.</given-names></name></person-group> (<year>2011</year>). <source><italic>Insect Ecology: Behavior, Populations and Communities.</italic></source> <publisher-loc>Cambridge</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>.</citation></ref>
<ref id="B63"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pringle</surname> <given-names>R. M.</given-names></name> <name><surname>Tarnita</surname> <given-names>C. E.</given-names></name></person-group> (<year>2017</year>). <article-title>Spatial self-organization of ecosystems: integrating multiple mechanisms of regular-pattern formation.</article-title> <source><italic>Annu. Rev. Entomol.</italic></source> <volume>62</volume> <fpage>359</fpage>&#x2013;<lpage>377</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-ento-031616-035413</pub-id> <pub-id pub-id-type="pmid">28141964</pub-id></citation></ref>
<ref id="B64"><citation citation-type="journal"><collab>R Core Team</collab> (<year>2019</year>). <source><italic>R: A Language and Environment for Statistical Computing.</italic></source> <publisher-loc>Vienna</publisher-loc>: <publisher-name>R Foundation for Statistical Computing</publisher-name>.</citation></ref>
<ref id="B65"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ripa</surname> <given-names>R.</given-names></name> <name><surname>Luppichini</surname> <given-names>P.</given-names></name></person-group> (<year>2004</year>). &#x201C;<article-title>Especies de termitas de importancia econ&#x00F3;mica presentes en Chile</article-title>,&#x201D; in <source><italic>Termitas y otros insectos xil&#x00F3;fagos en Chile: especies, biolog&#x00ED;a y manejo</italic></source>, <role>eds</role> <person-group person-group-type="editor"><name><surname>Ripa</surname> <given-names>R.</given-names></name> <name><surname>Luppichini</surname> <given-names>P.</given-names></name></person-group> (<publisher-loc>Quillota</publisher-loc>: <publisher-name>Instituto de Investigaciones Agropecuarias</publisher-name>), <fpage>73</fpage>&#x2013;<lpage>105</lpage>.</citation></ref>
<ref id="B66"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Roisin</surname> <given-names>Y.</given-names></name> <name><surname>Korb</surname> <given-names>J.</given-names></name></person-group> (<year>2011</year>). &#x201C;<article-title>Social organisation and the status of workers in termites</article-title>,&#x201D; in <source><italic>Biology of Termites: a Modern Synthesis</italic></source>, <role>eds</role> <person-group person-group-type="editor"><name><surname>Bignell</surname> <given-names>D. E.</given-names></name> <name><surname>Roisin</surname> <given-names>Y.</given-names></name> <name><surname>Lo</surname> <given-names>N.</given-names></name></person-group> (<publisher-loc>Dordrecht</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>133</fpage>&#x2013;<lpage>164</lpage>. <pub-id pub-id-type="doi">10.1111/brv.12330</pub-id> <pub-id pub-id-type="pmid">28508537</pub-id></citation></ref>
<ref id="B67"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sanada-Morimura</surname> <given-names>S.</given-names></name> <name><surname>Satoh</surname> <given-names>T.</given-names></name> <name><surname>Obara</surname> <given-names>Y.</given-names></name></person-group> (<year>2006</year>). <article-title>Territorial behavior and temperature preference for nesting sites in a pavement ant Tetramorium tsushimae.</article-title> <source><italic>Insectes Soc.</italic></source> <volume>53</volume> <fpage>141</fpage>&#x2013;<lpage>148</lpage>. <pub-id pub-id-type="doi">10.1007/s00040-005-0849-2</pub-id></citation></ref>
<ref id="B68"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Scharf</surname> <given-names>M. E.</given-names></name> <name><surname>Buckspan</surname> <given-names>C. E.</given-names></name> <name><surname>Grzymala</surname> <given-names>T. L.</given-names></name> <name><surname>Zhou</surname> <given-names>X.</given-names></name></person-group> (<year>2007</year>). <article-title>Regulation of polyphenic caste differentiation in the termite Reticulitermes flavipes by interaction of intrinsic and extrinsic factors.</article-title> <source><italic>J. Exp. Biol.</italic></source> <volume>210</volume> <fpage>4390</fpage>&#x2013;<lpage>4398</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.010876</pub-id> <pub-id pub-id-type="pmid">18232118</pub-id></citation></ref>
<ref id="B69"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Scholtz</surname> <given-names>O. I.</given-names></name> <name><surname>Macleod</surname> <given-names>N.</given-names></name> <name><surname>Eggleton</surname> <given-names>P.</given-names></name></person-group> (<year>2008</year>). <article-title>Termite soldier defence strategies: a reassessment of Prestwich&#x2019;s classification and an examination of the evolution of defence morphology using extended eigenshape analyses of head morphology.</article-title> <source><italic>Zool. J. Linn. Soc.</italic></source> <volume>153</volume> <fpage>631</fpage>&#x2013;<lpage>650</lpage>. <pub-id pub-id-type="doi">10.1111/j.1096-3642.2008.00396.x</pub-id></citation></ref>
<ref id="B70"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shellman-Reeve</surname> <given-names>J. S.</given-names></name></person-group> (<year>1997</year>). &#x201C;<article-title>The spectrum of eusociality in termites</article-title>,&#x201D; in <source><italic>The Evolution of Social Behavior in Insects and Arachnids</italic></source>, <role>eds</role> <person-group person-group-type="editor"><name><surname>Choe</surname> <given-names>J. C.</given-names></name> <name><surname>Crespi</surname> <given-names>B. J.</given-names></name></person-group> (<publisher-loc>Cambridge</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>), <fpage>52</fpage>&#x2013;<lpage>93</lpage>.</citation></ref>
<ref id="B71"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shelton</surname> <given-names>T. G.</given-names></name> <name><surname>Grace</surname> <given-names>J. K.</given-names></name></person-group> (<year>1996</year>). <article-title>Review of agonistic behaviors in the Isoptera.</article-title> <source><italic>Sociobiology</italic></source> <volume>28</volume> <fpage>155</fpage>&#x2013;<lpage>176</lpage>.</citation></ref>
<ref id="B72"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>K. F.</given-names></name> <name><surname>Stefaniak</surname> <given-names>L.</given-names></name> <name><surname>Saito</surname> <given-names>Y.</given-names></name> <name><surname>Gemmill</surname> <given-names>C. E. C.</given-names></name> <name><surname>Cary</surname> <given-names>S. C.</given-names></name> <name><surname>Fidler</surname> <given-names>A. E.</given-names></name></person-group> (<year>2012</year>). <article-title>Increased inter-colony fusion rates are associated with reduced COI haplotype diversity in an invasive colonial ascidian Didemnum vexillum.</article-title> <source><italic>PLoS One</italic></source> <volume>7</volume>:<issue>e30473</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0030473</pub-id> <pub-id pub-id-type="pmid">22303442</pub-id></citation></ref>
<ref id="B73"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>R. J.</given-names></name> <name><surname>Hines</surname> <given-names>A.</given-names></name> <name><surname>Richmond</surname> <given-names>S.</given-names></name> <name><surname>Merrick</surname> <given-names>M.</given-names></name> <name><surname>Drew</surname> <given-names>A.</given-names></name> <name><surname>Fargo</surname> <given-names>R.</given-names></name></person-group> (<year>2007</year>). <article-title>Altitudinal variation in body size and population density of Nicrophorus investigator (Coleoptera: Silphidae).</article-title> <source><italic>Environ. Entomol.</italic></source> <volume>29</volume> <fpage>290</fpage>&#x2013;<lpage>298</lpage>. <pub-id pub-id-type="doi">10.1093/ee/29.2.290</pub-id></citation></ref>
<ref id="B74"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Swaddle</surname> <given-names>J. P.</given-names></name> <name><surname>Witter</surname> <given-names>M. S.</given-names></name></person-group> (<year>1994</year>). <article-title>Food, feathers and fluctuating asymmetries.</article-title> <source><italic>Proc. R Soc. Lond. B</italic></source> <volume>255</volume> <fpage>147</fpage>&#x2013;<lpage>152</lpage>. <pub-id pub-id-type="doi">10.1098/rspb.1994.0021</pub-id></citation></ref>
<ref id="B75"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thomas</surname> <given-names>M. L.</given-names></name> <name><surname>Tsutsui</surname> <given-names>N. D.</given-names></name> <name><surname>Holway</surname> <given-names>D. A.</given-names></name></person-group> (<year>2005</year>). <article-title>Intraspecific competition influences the symmetry and intensity of aggression in the Argentine ant.</article-title> <source><italic>Behav. Ecol.</italic></source> <volume>16</volume> <fpage>472</fpage>&#x2013;<lpage>481</lpage>. <pub-id pub-id-type="doi">10.1093/beheco/ari014</pub-id></citation></ref>
<ref id="B76"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thorne</surname> <given-names>B. L.</given-names></name></person-group> (<year>1996</year>). <article-title>Termite Terminology.</article-title> <source><italic>Sociobiology</italic></source> <volume>28</volume> <fpage>253</fpage>&#x2013;<lpage>263</lpage>. <pub-id pub-id-type="doi">10.1111/j.1469-185x.2008.00044.x</pub-id> <pub-id pub-id-type="pmid">18979593</pub-id></citation></ref>
<ref id="B77"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thorne</surname> <given-names>B. L.</given-names></name> <name><surname>Breisch</surname> <given-names>N. L.</given-names></name> <name><surname>Muscedere</surname> <given-names>M. L.</given-names></name></person-group> (<year>2003</year>). <article-title>Evolution of eusociality and the soldier caste in termites: influence of intraspecific competition and accelerated inheritance.</article-title> <source><italic>Proc. Natl. Acad. Sci. U. S. A.</italic></source> <volume>100</volume> <fpage>12808</fpage>&#x2013;<lpage>12813</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.2133530100</pub-id> <pub-id pub-id-type="pmid">14555764</pub-id></citation></ref>
<ref id="B78"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tsuji</surname> <given-names>K.</given-names></name></person-group> (<year>2013</year>). <article-title>Kin selection, species richness and community.</article-title> <source><italic>Biol. Lett.</italic></source> <volume>9</volume> <fpage>11</fpage>&#x2013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.1098/rsbl.2013.0491</pub-id> <pub-id pub-id-type="pmid">24132096</pub-id></citation></ref>
<ref id="B79"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Warren</surname> <given-names>M.</given-names></name> <name><surname>McGeoch</surname> <given-names>M. A.</given-names></name> <name><surname>Nicolson</surname> <given-names>S. W.</given-names></name> <name><surname>Chown</surname> <given-names>S. L.</given-names></name></person-group> (<year>2006</year>). <article-title>Body size patterns in Drosophila inhabiting a mesocosm: interactive effects of spatial variation in temperature and abundance.</article-title> <source><italic>Oecologia</italic></source> <volume>149</volume> <fpage>245</fpage>&#x2013;<lpage>255</lpage>. <pub-id pub-id-type="doi">10.1007/s00442-006-0434-z</pub-id> <pub-id pub-id-type="pmid">16773332</pub-id></citation></ref>
<ref id="B80"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wills</surname> <given-names>B. D.</given-names></name> <name><surname>Moreau</surname> <given-names>C. S.</given-names></name> <name><surname>Wray</surname> <given-names>B. D.</given-names></name> <name><surname>Hoffmann</surname> <given-names>B. D.</given-names></name> <name><surname>Suarez</surname> <given-names>A. V.</given-names></name></person-group> (<year>2014</year>). <article-title>Body size variation and caste ratios in geographically distinct populations of the invasive big-headed ant, Pheidole megacephala (Hymenoptera: Formicidae).</article-title> <source><italic>Biol. J. Linn. Soc.</italic></source> <volume>113</volume> <fpage>423</fpage>&#x2013;<lpage>438</lpage>. <pub-id pub-id-type="doi">10.1111/bij.12386</pub-id></citation></ref>
<ref id="B81"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Witter</surname> <given-names>M. S.</given-names></name> <name><surname>Swaddle</surname> <given-names>J. P.</given-names></name></person-group> (<year>1994</year>). <article-title>Fluctuating asymmetries, competition and dominance.</article-title> <source><italic>Proc. R Soc. Lond. B</italic></source> <volume>256</volume> <fpage>299</fpage>&#x2013;<lpage>303</lpage>. <pub-id pub-id-type="doi">10.1098/rspb.1994.0084</pub-id></citation></ref>
<ref id="B82"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zizka</surname> <given-names>G.</given-names></name> <name><surname>Schneider</surname> <given-names>J. V.</given-names></name> <name><surname>Schulte</surname> <given-names>K.</given-names></name> <name><surname>Novoa</surname> <given-names>P.</given-names></name></person-group> (<year>2013</year>). <article-title>Taxonomic revision of the Chilean Puya species (Puyoideae, Bromeliaceae), with special notes on the Puya alpestris-Puya berteroniana species complex.</article-title> <source><italic>Brittonia</italic></source> <volume>65</volume> <fpage>387</fpage>&#x2013;<lpage>407</lpage>. <pub-id pub-id-type="doi">10.1007/s12228-012-9290-9</pub-id></citation></ref>
<ref id="B83"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zobel</surname> <given-names>M. U.</given-names></name> <name><surname>Paxton</surname> <given-names>R. J.</given-names></name></person-group> (<year>2007</year>). <article-title>Is big the best? Queen size, usurpation and nest closure in a primitively eusocial sweat bee (Lasioglossum malachurum).</article-title> <source><italic>Behav. Ecol. Sociobiol.</italic></source> <volume>6&#x2018;</volume> <fpage>35</fpage>&#x2013;<lpage>447</lpage>. <pub-id pub-id-type="doi">10.1007/s00265-006-0271-4</pub-id></citation></ref>
</ref-list>
</back>
</article>
