<title>Background</title> <p>Ants of the Northern-hemispheric, temperate genus <italic>Lasius </italic>(Formicinae) are scientifically significant, in terms of relative abundance and ecological impact [##UREF##0##1##,##UREF##1##2##]. Because of the diversity of their signal and defense chemistry, <italic>Lasius </italic>ants are organisms widely used in chemical ecology [##UREF##1##2##, ####UREF##2##3##, ##UREF##3##4##, ##UREF##4##5##, ##UREF##5##6##, ##UREF##6##7####6##7##] and the wide range of colony organisations makes the genus an ideal system for exploring social evolution [##REF##17964165##8##, ####UREF##7##9##, ##UREF##8##10##, ##REF##14635899##11##, ##UREF##9##12####9##12##]. Two further complex traits found in <italic>Lasius </italic>are yet to be adequately understood: social parasitism and fungiculture.</p> <p>Social parasitism implies that one eusocial species depends on the labour force of another [##REF##21227804##13##, ####UREF##10##14##, ##UREF##11##15####11##15##]. The social parasitism exhibited in <italic>Lasius </italic>is temporary in that it is confined to the early stages of the parasite's colony: the parasitic queen founds her colony through entering a host colony where she kills the resident queen and takes over the worker force [##UREF##0##1##,##UREF##1##2##,##REF##21227804##13##]. The study of social parasitism has become virtually a little discipline of entomology in itself [##UREF##1##2##], but the conditions for social parasitism to arise remain poorly understood [##REF##21227804##13##, ####UREF##10##14##, ##UREF##11##15##, ##UREF##12##16##, ##UREF##13##17##, ##REF##16674574##18##, ##UREF##14##19##, ##UREF##15##20####15##20##]. Social parasitism has evolved many times independently in ants [##UREF##1##2##,##REF##21227804##13##,##UREF##11##15##,##UREF##16##21##], but the evolutionary trajectories at finer systematic scale, e.g., whether it evolved once or multiply within genera, have only recently received detailed attention [##UREF##11##15##,##REF##16674574##18##,##UREF##16##21##, ####REF##12740442##22##, ##UREF##17##23##, ##UREF##18##24##, ##UREF##19##25##, ##UREF##20##26##, ##UREF##21##27####21##27##]. Fungiculture by ants, termites and beetles, on the other hand, provides a powerful study system for studying the origin and maintenance of mutualism [##UREF##22##28##]. In ants, fungiculture has evolved independently at least twice: in attines (members of the Myrmicinae), which culture the fungi for food, and in <italic>Lasius </italic>ants, which use fungi to build composite nest walls [##UREF##23##29##, ####UREF##24##30##, ##UREF##25##31##, ##UREF##26##32##, ##UREF##27##33####27##33##]. The patterns of diversification in the intensely studied attine fungiculture are only gradually starting to be understood, as brought out by recent papers on leaf-cutter fungiculture which reverse earlier impressions of certainty for some important issues [##REF##16815974##34##,##REF##18362345##35##]. For inferences on the evolution of the outstanding ecological and social traits including social parasitism and fungiculture in <italic>Lasius</italic>, a well-founded phylogeny of the genus is needed. There have been three previous studies to resolve the phylogeny of <italic>Lasius </italic>[##UREF##28##36##, ####UREF##29##37##, ##REF##15522790##38####15522790##38##], but these have disagreed with each other in significant respects (Fig. ##FIG##0##1##).</p> <p>In the present study we attempt to establish a robust phylogenetic framework for the relationships of the <italic>Lasius </italic>subgenera. We apply Bayesian analysis, a powerful tool in phylogenetic reconstruction of combined data [##REF##14965900##39##] and not previously applied to <italic>Lasius</italic>. In addition, we also apply Maximum Parsimony analysis (MP); MP represents a completely different computational technique for phylogenetic reconstruction [##REF##16701310##40##] and agreement of the reconstructions by the two independent methods would increase confidence in the tree. Our approach comprised five steps. (<italic>i</italic>) We combined evidence from different data sets, which for many organisms, including ants, often yields a stronger phylogenetic signal than using the data sets singly [##REF##16674574##18##,##REF##14965900##39##,##UREF##30##41##, ####REF##17079492##42##, ##REF##9479692##43##, ##UREF##31##44##, ##UREF##32##45##, ##REF##16351964##46##, ##REF##16631389##47##, ##UREF##33##48##, ##REF##16601190##49##, ##REF##17408491##50####17408491##50##]; specifically, we combined mitochondrial DNA sequence and morphological data. (<italic>ii</italic>) We explored potential causes of distortion of the molecular phylogenetic signal, namely substitution saturation [##REF##10486008##51##], positive selection [##REF##16716603##52##,##REF##17257111##53##], and compositional heterogeneity [##UREF##34##54##]. We also assessed which morphological characters may be functionally coupled with social parasitism [##UREF##11##15##,##UREF##34##54##,##UREF##35##55##], because similarities in those characters could reflect convergent adaptations to social parasitic life rather than reflect common ancestry [##UREF##36##56##,##REF##16690327##57##]. We then excluded any suspected cause of distortion in the phylogenetic reconstruction. To address potential issues of character exclusion [##UREF##37##58##], we explored the effect of excluding those morphological characters suspected to be functionally coupled with social parasitism by repeating the reconstruction when including them. (<italic>iii</italic>) We explored whether the topologies as inferred from the different data sets are in statistically significant conflict with the phylogenetic framework inferred from the concatenated data which would tend to reduce confidence in the latter [##REF##18432550##59##,##REF##17366134##60##]. Our approach was to handle any conflict arising by collapsing the affected node. (<italic>iv</italic>) We explored whether any of the previous <italic>Lasius </italic>phylogenies fit our data set as well as our resulting topology. (<italic>v</italic>) Finally, we used posterior mapping [##UREF##38##61##,##REF##12746144##62##] to define hypotheses on the evolution of social parasitism and fungiculture in <italic>Lasius</italic>. In all, the resulting topology provides a solid basis for studying the evolution of the various ecologically and sociobiologically relevant traits of <italic>Lasius </italic>across subgenera.</p>
<title>Methods</title> <title>The study system</title> <p>The species of the ant genus <italic>Lasius </italic>are currently placed in six subgenera, <italic>Acanthomyops</italic>, <italic>Austrolasius</italic>, <italic>Cautolasius</italic>, <italic>Chthonolasius</italic>, <italic>Dendrolasius</italic>, and <italic>Lasius sensu stricto </italic>[##UREF##44##70##]. The most recent taxonomic revision at the genus level dates back to 1955 [##UREF##28##36##]. It recognised four subgenera: <italic>Cautolasius</italic>, <italic>Chthonolasius</italic>, <italic>Dendrolasius</italic>, and <italic>Lasius sensu stricto</italic>. The fifth subgenus, <italic>Austrolasius</italic>, was established in 1967 [##UREF##44##70##] and includes one species which had up to then been placed in <italic>Chthonolasius</italic>. The sixth subgenus, <italic>Acanthomyops</italic>, has an unstable history in that its status changed several times between that of a separate genus and being a subgenus, mostly of <italic>Lasius </italic>[##UREF##44##70##]. Evidence for its inclusion into <italic>Lasius </italic>has accumulated [##REF##15522790##38##,##UREF##48##76##] and it was therefore formally returned to <italic>Lasius </italic>[##UREF##42##66##]. Monophyly of the subgenera is supported by various pieces of evidence [##UREF##49##77##], except for <italic>Lasius sensu stricto</italic>. <italic>Lasius sensu stricto </italic>harbours one taxon, <italic>L. pallitarsis</italic>, which on morphological and molecular grounds has been hypothesised to best constitute a separate subgenus [##REF##15522790##38##]. We here validate these earlier findings [##REF##15522790##38##].</p> <p>Social parasitism is confined to four subgenera, <italic>Acanthomyops</italic>, <italic>Austrolasius</italic>, <italic>Chthonolasius</italic>, and <italic>Dendrolasius</italic>, with all species of these subgenera obligatorily displaying this lifestyle [##UREF##0##1##,##UREF##42##66##]. One subgenus, <italic>Dendrolasius</italic>, is hyperparasitic in that it parasitises parasitic <italic>Chthonolasius </italic>species [##UREF##0##1##,##UREF##28##36##]. Fungiculture is known from two subgenera, <italic>Chthonolasius</italic>, and <italic>Dendrolasius </italic>[##REF##16815974##34##] and, as far as known, all members of the two subgenera use fungi (B.C. Schlick-Steiner, unpublished data; M. Maruyama, unpublished data; [##UREF##0##1##]), although New World species still await study.</p> <title>Taxon sampling</title> <p>The material of this study is listed in Table ##TAB##0##1##. It comprises 27 species of <italic>Lasius </italic>(including one undescribed) representing all 6 of the subgenera currently recognised as well as <italic>Lasius pallitarsis</italic>. We thus present data on more than a quarter of the currently 100 valid extant <italic>Lasius </italic>species [##UREF##44##70##], covering both the Palaearctic and the Nearctic regions. The number of species per subgenus ranged from one to 11. Only two and three species each of <italic>Chthonolasius </italic>and <italic>Acanthomyops</italic>, respectively, were analysed, with each sample carefully chosen to be unambiguously identified to species (because the species of these subgenera are known for their habitual multidirectional hybridisation [##REF##16995614##78##,##UREF##50##79##] which is likely to compromise the resolving power of gene sequences [##REF##16503279##80##] and distorts morphological characters [##UREF##50##79##]). To account for intraspecific diversity, up to four colonies per species were included wherever possible. Two species of <italic>Myrmecocystus</italic>, which genus is sister to <italic>Lasius </italic>[[##REF##17079492##42##,##REF##16601190##49##]; confirmed by a personal communication by P.S. Ward, of February 2008], as well as <italic>Formica japonica </italic>were used as outgroup. Species were identified according to [##UREF##0##1##,##UREF##28##36##,##UREF##49##77##,##UREF##51##81##, ####UREF##52##82##, ##UREF##53##83####53##83##]. Voucher specimens are deposited at the National Science Museum, Tokyo under the voucher numbers listed in Table ##TAB##0##1##.</p> <title>Molecular protocols</title> <p>Some samples were irreplaceable dried museum specimens. Therefore, genomic DNA was extracted from whole body using a DNAeasy tissue kit (Qiagen, Hilden, Germany) using established protocols [##UREF##54##84##] without any damage to the voucher specimens. We used 1 μl of DNA (25 – 50 ng/μl) as template for PCR amplification. A 490 – 550 bp region of <italic>16S rRNA </italic>was amplified and sequenced using primers "16Sar-L" 5'-CGCCTGTTTATCAAAAACAT-3' and "16Sar-L2" 5'-CCGGTCTGAACTCAGATCATG-3' originally taken from [##UREF##55##85##] but the latter one slightly altered and thus renamed. A ca. 900 bp region of <italic>cox1 </italic>was amplified and sequenced using the new, degenerated primers "Lasius-L" 5'-TAYCCGCCATTAGCTTCAAA-3' and "Lasius-R" 5'-TGAAATTAAGGATCCAATWGA-3'. Reactions were carried out at 10 μl volumes in a PCR Thermal Cycler MP (TaKaRa Bio Inc.) under the following conditions: a first cycle of 94°C for 3 min, followed by 35 cycles of 94°C for 30 s, annealing at 50°C for 50 s, and finally 72°C for 1 min for the <italic>16S rRNA</italic>; for <italic>cox1 </italic>all settings were identical except for annealing which was set to 42°C for 1 min 15 s. PCR products were purified with 0.5 μl of ExoSap-IT (GE Healthcare Life Sciences). All products were sequenced in both directions using BigDye Terminator v3.1 (Applied Biosystems) on an ABI 3100 Avant DNA Sequencer (Applied Biosystems) at the National Science Museum, Tokyo. The sequence data were deposited at DNA Data Base of Japan, DDBJ (see Table ##TAB##0##1## for accession numbers).</p> <title>Exploration of molecular data concerning potential causes for phylogenetic distortion</title> <p>The <italic>16S rRNA </italic>and <italic>cox1 </italic>sequences were aligned with default settings of the program Clustal X v1.83 [##REF##9396791##86##]; ambiguously aligned sites in the <italic>16S rRNA </italic>alignment were excluded. We partitioned the <italic>cox1 </italic>sequences into the first, second and third codon positions using the program DAMBE v4.2.13 [##UREF##56##87##]. This program was also used to perform tests for the saturation of substitutions [##REF##12470932##88##] on the <italic>cox1 </italic>and <italic>16S rRNA </italic>data. For <italic>cox1 </italic>all codon positions were tested simultaneously, as well as separately. We found no indication of substitution saturation (see Results for detail), but we additionally performed the <italic>cox1 </italic>MCMC analysis without the third codon position.</p> <p>To detect potential positive selection we used the program HYPHY [##REF##15713735##89##] accessed through the Datamonkey interface <ext-link ext-link-type="uri" xlink:href="http://www.datamonkey.org">http://www.datamonkey.org</ext-link>. Mean numbers of nonsynonymous substitutions (dN) and synonymous substitutions (dS) per site (ratio dN/dS) were estimated in the <italic>cox1 </italic>data using the fixed effect (two-rate FEL) method and basing estimates on a Neighbour-Joining tree under the HKY substitution model. We used a nominal alpha level of 0.1.</p> <p>To avoid any effect from compositional heterogeneity of sequences [##REF##15371251##90##] on the phylogenetic reconstruction we separately tested each codon position of <italic>cox1 </italic>as well as <italic>16S rRNA </italic>using the program TREE-PUZZLE 5.2 [##REF##11934758##91##]. When we found indications for compositional heterogeneity we recoded the sequences into purines and pyrimidines and repeated the test.</p> <title>Morphological analysis</title> <p>The morphological characters considered are presented in Additional File ##SUPPL##0##1##. Our rationale in composing the character set was aiming at, first, a comprehensive capturing of variation at the species, subgenus and genus level, and, second, exclusion of any characters potentially distorting the phylogenetic reconstructions. We started from the complete set of character definitions of Janda et al. [##REF##15522790##38##]. Pursuing our first aim, we applied 67 characters and their states exactly as described by [##REF##15522790##38##], adapted – because analysis of our material revealed the necessity to do so – the definitions of character states, partly also concerning the number of character states, for another 16 characters, out of which 7 characters also were adapted in the character definitions themselves, and added 20 entirely new characters. Pursuing our second aim, we excluded from the set presented by [##REF##15522790##38##] the three behavioural/ecological characters (among others the occurrence of social parasitism), and on the basis of information from [##UREF##0##1##,##UREF##1##2##,##REF##21227804##13##,##UREF##34##54##,##UREF##35##55##,##UREF##39##63##,##UREF##40##64##] we further excluded 37 morphological characters, including four of the new ones. The excluded characters were characters which we suspected to be functionally coupled with temporary social parasitism, either (<italic>i</italic>) because of direct functional reasoning, or (<italic>ii</italic>) because they are known to be correlated with social parasitism in other ant genera which contain both parasitic and non-parasitic species. The characters excluded under (<italic>i</italic>) concerned the mandible, which is needed by parasitic <italic>Lasius </italic>queens to dismember host workers and strangle the host queen; the maxillary palp and the scape, which may be under selective pressure to be short and robust so as to escape damage by aggressive hosts in the initial stage of colony take-over; and the mesosoma size, because parasitic, i.e., dependently founding, queens need less tissue storage than independently founding queens. The characters excluded under (<italic>ii</italic>) concerned the length of hairs, which can be either extremely long or extremely short in parasites; the pubescence, which is absent in some parasites; the size of the head and the overall body, which both are frequently small in parasitic queens; and the shape of the petiole, which is frequently aberrant in parasites. We finally also excluded three characters for which we observed no variation in our material (their general lack of variation for the species analysed confirmed by a personal communication by B. Seifert, of April 2008). For details of how we composed our morphological data set including information on the provenance of characters and whether we excluded them, see Additional File ##SUPPL##0##1##. Overall, we included 64 morphological characters in our reconstructions; 35 of these concern adult workers, 18 adult queens, 23 adult males, and 2 worker larvae; 47 are binary and 17 are multi-state. To explore the effect of excluding characters, as we had done, on phylogenetic reconstructions, we also subjected the complete data set to reconstructions (but excluding invariant characters); this data set included 99 characters. We assessed any morphological data by analysis of voucher material housed in the National Science Museum, Tokyo, except the three characters concerning larval morphology. A total of 155 specimens were analysed. All multistate characters were treated as unordered.</p> <title>Phylogenetic reconstructions</title> <p>To select the best-fitting nucleotide substitution models for <italic>cox1 </italic>and <italic>16S rRNA </italic>we used the hierarchical likelihood ratio test (hLRT) and the Akaike Information Criterion (AIC) implemented in the program MrModeltest 2.2 [##UREF##57##92##]. When sequence partitions had been recoded into purines and pyrimidines, models were adjusted to account for the two state character of the data. AIC and hLRT in one instance differed in the selection of models for the single DNA sequence data partitions (Table ##TAB##1##2##). AIC selected the more parameter rich model and we opted for this solution as erring on the side of overparameterisation is preferable over the opposite [##UREF##58##93##,##REF##15764559##94##].</p> <p>For the morphological data the Markov <italic>k </italic>(Mk) model [##REF##12116640##95##] was applied both with (+Γ) and without gamma-distributed rates of character change in separate MCMC runs. We used Bayes factors for model selection as they are established to provide good orientation tools in this [##REF##14965900##39##] and calculated them as follows, using the outcomes of the single MCMC runs: 2LnB<sub>10 </sub>= 2 × (Harmonic Mean Ln likelihood for Mk – Harmonic Mean Ln likelihood for Mk+Γ). In the interpretation of the yielded absolute value of 2.5 in favour of the Mk+Γ model we followed published recommendations [##UREF##59##96##], on page 777, and consequently used the Mk+Γ model for all reconstructions.</p> <p>Bayesian analysis using MCMC was performed with MrBayes 3.1.2 [##REF##12912839##97##] on the individual data sets (<italic>cox1</italic>, <italic>16S rRNA</italic>, morphology) and the combined, concatenated data set (<italic>cox1 </italic>plus <italic>16S rRNA </italic>plus morphology). We also analysed the combined, concatenated data set including those morphological characters suspected to be coupled with social parasitism. In addition, we analysed the concatenated molecular data (<italic>cox1 </italic>plus <italic>16S rRNA</italic>) without any morphological data. Data partitions were established to allow model parameters to be separately estimated for all partitions and additionally for the single codon positions of <italic>cox1</italic>. 10,000,000 generations with a sample frequency set to 100 were run. As after 9,000,000 generations stationarity was achieved with average standard deviation of split frequencies in all cases constantly below 0.002 except the reconstruction with the W55 constraint for which it was 0.063 (Table ##TAB##2##3##), we always used the last 10,000 trees of each run to compute a majority rule consensus tree assigning posterior probabilities of tree topology. We also confirmed that true convergence had been reached and that the MCMC was sampling from the posterior distribution by repeating all runs three times and checking for congruence across the runs. All runs were performed using parallel versions of MrBayes, implemented on a SGI Origin 3800 under IRIX version 6.21m, of HPC, James Cook University. All MCMC runs achieved stationarity and detailed statistics on the runs are presented in Table ##TAB##2##3##. In the interpretation of the MCMC trees we followed previous authors [##REF##15764559##94##,##REF##12414316##98##,##UREF##60##99##] to regard only nodes with node support of p &gt; 0.95 as significantly supported in Bayesian analysis. We also applied this cutoff when comparing the MCMC trees based on the individual data sets and that based on the concatenated data.</p> <p>We also performed MP analysis of the combined, concatenated data, as well as of the combined concatenated data adding those morphological characters suspected to be coupled with social parasitism. All MP analyses were unweighted and performed with PAUP* 4.0b10 [##UREF##61##100##] using the heuristic search algorithm with tree bisection reconnection branch swapping and 10 random stepwise additions. All characters were treated as unordered and polymorphic states were taken into account. Node support was calculated by 1,000 bootstrap replicates. In the interpretation of the MP trees, we applied the widely accepted node support threshold of &gt; 70 [##UREF##43##67##].</p> <p>We deposited the aligned, concatenated data matrix with TreeBase (Study accession number S2136).</p> <title>Comparison with previous <italic>Lasius </italic>phylogenies</title> <p>To compare the subgenus relationships of W55, H98 and J04 directly with the new, Bayesian framework, we enforced the various topologies as constraints (Fig. ##FIG##0##1##) on our concatenated data. We then performed additional MCMC runs of our concatenated data under these constraints and compared the outcomes of the single MCMC runs using Bayes factors as given under morphological analyses. In extracting the topologies from the literature we proceeded as follows. For W55 we used the tree of "Fig. ##FIG##1##2##" of [##UREF##28##36##] slightly modified. We adopted the position of <italic>Acanthomyops</italic>, explicitly treated as ingroup of <italic>Lasius </italic>in the tree, although Wilson did not formally treat <italic>Acanthomyop</italic>s as a subgenus of <italic>Lasius </italic>in the taxonomic revision itself. We allocated <italic>Lasius sitkaensis</italic>, now treated as a junior synonym of <italic>L. pallitarsis</italic>, to <italic>Lasius sensu stricto</italic>. A similar situation pertains to <italic>Austrolasius</italic>: the subgenus had not yet been established in 1955, but the then only known species which today is treated under <italic>Austrolasius</italic>, <italic>L. carniolicus</italic>, was allocated to <italic>Chthonolasius </italic>and we accounted for this in our treatment of W55. For H98 we applied the node support threshold of &gt; 70 [##UREF##43##67##] to the Maximum Parsimony reconstruction presented in "Fig. ##FIG##0##1##" of [##UREF##29##37##], and we applied the same threshold for J04, to the Maximum Parsimony reconstruction presented in "Fig. 6" of [##REF##15522790##38##].</p> <title>Posterior mapping analysis</title> <p>To estimate the probabilities of the possible ancestral states at each well supported node of the concatenated Bayesian topology we chose the Bayesian approach of posterior mapping [##UREF##38##61##,##REF##12746144##62##], using the program SIMMAP 1.0 [##REF##16504105##101##] freely available online <ext-link ext-link-type="uri" xlink:href="http://www.SIMMAP.com">http://www.SIMMAP.com</ext-link>. In contrast to parsimony approaches to character mapping this is a probabilistic approach, which (<italic>i</italic>) does not assume that only a single change has occurred along any branch, and (<italic>ii</italic>) is not prone to underestimation of the variance in ancestral state assignments [##REF##16504105##101##]. Further, the SIMMAP approach allows uncertainty in phylogenetic reconstruction. A single stochastic mapping was done per tree using the last 2,000 post-burnin trees of the 20,000 trees used to derive the consensus tree of Fig. ##FIG##1##2## and the ancestral states were inferred for the consensus tree from the MrBayes analysis.</p>
<title>Results</title> <p>We found no evidence for saturation of substitutions, neither for the <italic>16S ribosomal RNA </italic>(<italic>16S rRNA</italic>) nor for any of the codon positions, single or in combination, of the <italic>cytochrome c oxidase subunit I </italic>(<italic>cox1</italic>) data. Tests for positive selection within the <italic>cox1 </italic>data indicated that none of the sites was subject to positive selection. But we did detect compositional heterogeneity of sequences for the third codon position of <italic>cox1</italic>. Recoding the nucleotides of these sites to purines and pyrimidines successfully eliminated this effect and we used the recoded sequences for all phylogenetic reconstructions. Scrutinising the morphological data set for characters potentially coupled functionally with social parasitism [##UREF##0##1##,##UREF##1##2##,##REF##21227804##13##,##UREF##34##54##,##UREF##35##55##,##UREF##39##63##,##UREF##40##64##] yielded 37 characters which we hence excluded from our morphological data set. The final data set consisted of 48 samples of 30 species including three outgroup species (Table ##TAB##0##1##) for which a total of 1,265 base pairs (bp), and 64 morphological characters were used for phylogenetic reconstruction (Table ##TAB##1##2##).</p> <p>Topology, branch lengths, and Bayesian posterior probabilities of the Markov Chain Monte Carlo (MCMC) analyses of the concatenated data (<italic>cox1 </italic>plus <italic>16S rRNA </italic>plus morphology) are given in Fig. ##FIG##1##2##. Monophyly of all subgenera was strongly supported (0.98 – 1.00), as were all nodes defining subgenus relationships (0.99 – 1.00), with exception of one (0.83), connecting <italic>Cautolasius </italic>to (<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>). Strictly applying a cutoff for node support of &gt; 0.95, we retrieved two major lineages, (((<italic>Austrolasius </italic>+ <italic>Acanthomyops</italic>) + <italic>Chthonolasius</italic>) + <italic>Lasius pallitarsis</italic>) and ((<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>), but the position of <italic>Cautolasius </italic>relative to these two lineages remains unresolved.</p> <p>The three individual data set phylogenies (<italic>cox1</italic>, <italic>16S rRNA</italic>, morphology) differed considerably in phylogenetic resolution at subgenus level and above and none of them achieved resolution of all nodes scored by the concatenated data reconstruction. Repeating reconstruction of <italic>cox1 </italic>without the third codon position did not yield any well supported nodes that contradicted well supported nodes in the reconstruction using all three codon positions, which confirms the reliability of the result from the applied test for substitution saturation. Comparison of the individual data set topologies with the concatenated data topology at subgenus level and above revealed not a single significant (for MCMC, posterior probability &gt; 0.95) disagreement but rather agreement on well supported nodes. The following nodes of the concatenated topology were supported by the individual data set topologies: <italic>cox1 </italic>- <italic>Cautolasius </italic>monophyly, <italic>Chthonolasius </italic>monophyly, and monophyly of the two <italic>Myrmecocystus </italic>outgroup species; <italic>16S rRNA </italic>- <italic>Cautolasius </italic>monophyly, (<italic>Austrolasius </italic>+ <italic>Acanthomyops</italic>), and (<italic>Chthonolasius </italic>+ <italic>Lasius pallitarsis </italic>+ (<italic>Austrolasius </italic>+ <italic>Acanthomyops</italic>)); morphology - <italic>Dendrolasius </italic>monophyly and (<italic>Austrolasius </italic>+ <italic>Acanthomyops</italic>). Including those morphological characters suspected to be functionally coupled with social parasitism, which we therefore had excluded before, resulted in an identical topology, with very similar node support values (Fig. ##FIG##1##2##), but with a decrease of the value for the node connecting <italic>Cautolasius </italic>to (<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>) from 0.83 to 0.74. MP reconstructions based on the concatenated data set significantly (node support threshold of &gt; 70 for MP) confirmed all significant nodes of the Bayesian reconstructions, including retrieving as significant the node connecting <italic>Cautolasius </italic>to (<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>). This pattern persisted in the MP analysis when using the concatenated data including those morphological characters suspected to be functionally coupled with social parasitism, with one exception: the node connecting <italic>Dendrolasius </italic>to <italic>Lasius sensu stricto </italic>was then no longer supported (Fig. ##FIG##1##2##). Considering the lack of significant disagreement between the signals in the individual data set phylogenies and the confirmation of the Bayesian topology by Maximum Parsimony reconstructions, we regard the Bayesian concatenated data topology as a phylogenetic framework sufficiently robust for use in subsequent evolutionary hypotheses ("new" topology in Fig. ##FIG##0##1##).</p> <p>Subsequently, we also made a Bayesian reconstruction of the concatenated molecular data sets (<italic>cox1 </italic>plus <italic>16S rRNA</italic>), without any morphological data, to allow for DNA based estimates of branch length. This tree (inset in Fig. ##FIG##1##2##) was in large agreement with the other Bayesian trees, in topology and branch lengths, but in addition to the node connecting <italic>Cautolasius </italic>to (<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>), three previously well supported nodes now lacked significant support (<italic>Acanthomyops </italic>monophyly; node connecting <italic>Chthonolasius </italic>to (<italic>Austrolasius </italic>+ <italic>Acanthomyops</italic>); <italic>Lasius sensu stricto </italic>monophyly).</p> <p>In additional rounds of reconstruction, we enforced the subgenus relationships recovered in three previous phylogenetic analyses [##UREF##28##36##, ####UREF##29##37##, ##REF##15522790##38####15522790##38##] as constraints (Fig. ##FIG##0##1##) on our concatenated data. All reconstructions with constraints had lower likelihood values than the reconstruction without constraint (Table ##TAB##2##3##) and, on the basis of our concatenated data, Bayes factor analysis revealed very strong evidence against all previously suggested subgenus relationships (Table ##TAB##3##4##).</p> <p>The Bayesian posterior probabilities for the occurrence of social parasitism and fungiculture at the nodes of the phylogenetic framework above subgenus level ("new" topology in Fig. ##FIG##0##1##) are shown in Table ##TAB##4##5##. For eight of the ten possible inferences a state was significantly inferred (p &gt; 0.95), for the remaining two the probability values were 0.70 and 0.82.</p>
<title>Discussion and Conclusion</title> <title>The new phylogenetic framework</title> <p>The new phylogenetic framework we present clarifies the relationships of all but one subgenus, namely the position of <italic>Cautolasius </italic>relative to the other subgenera: MP reconstructions significantly support the topology, <italic>Cautolasius </italic>+ (<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>), and Bayesian reconstructions do not contradict it. Given the lack of significant support in the Bayesian tree, we nevertheless subsequently refrain from considering the node as resolved ("new" topology in Fig. ##FIG##0##1##). This lack of significant resolution does not affect considerations on the evolution of social parasitism and fungiculture, though, because the ancestral state of <italic>Lasius </italic>was absence of both traits (Table ##TAB##4##5##), and because also <italic>Cautolasius </italic>displays neither trait. A significantly supported sister-status to (<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>) thus would not alter any conclusions (see discussion below). The previous phylogenies were recovered by various methods and we discuss below how the methods applied may have influenced the respective results. Taken together, the reasons for the increased robustness that we postulate for the new phylogenetic framework, which robustness is highlighted also by the confirmation of our Bayesian topology by our MP reconstructions, may be that (1) we excluded potential causes of phylogenetic distortion and (2) used additional analysis methods, that (3) the phylogenetic signal of combined data is potentially stronger than that of individual data sets [##REF##16674574##18##,##REF##14965900##39##,##UREF##30##41##, ####REF##17079492##42##, ##REF##9479692##43##, ##UREF##31##44##, ##UREF##32##45##, ##REF##16351964##46##, ##REF##16631389##47##, ##UREF##33##48##, ##REF##16601190##49##, ##REF##17408491##50####17408491##50##], especially when applying Bayesian inference [##REF##14965900##39##], that (4) all nodes but one at subgenus level or above were significantly supported in our concatenated analysis (Fig. ##FIG##1##2##), and that (5) we did not recover any disagreement between the individual data set phylogenies and the concatenated data topologies. Nevertheless, because the molecular data of this paper all derive from mitochondrial DNA, we cannot absolutely exclude that the tree is influenced by introgression at very shallow levels. Nuclear pseudogenes are also a concern, but we consider it unlikely that we amplified these because neither reading frame shifts nor sequence ambiguities were apparent. The possibilities of incomplete lineage sorting and of selection driven by symbionts which are in disequilibrium with mtDNA [##UREF##41##65##] cannot be ruled out, but such effects are unlikely to confound clade history at a deeper level, such as that of subgenera. It remains true that future studies using multiple nuclear genes are desirable to confirm our findings.</p> <p>The new phylogenetic framework presented here confirms the monophyly of the six subgenera. The new framework also affirms the taxonomic placement of <italic>Acanthomyops </italic>as a subgenus of <italic>Lasius </italic>[##UREF##42##66##]. Moreover, there is additional evidence for treating <italic>Lasius pallitarsis </italic>as a separate, monotypic subgenus, as suggested earlier [##REF##15522790##38##]. Detailed morphological characterisation of the new subgenus and taxonomic implications will be followed up elsewhere in the frame of a formal taxonomic revision.</p> <p>The evolution of the <italic>Lasius </italic>subgenera occurred in two major lineages, the first lineage comprising <italic>Acanthomyops</italic>, <italic>Austrolasius</italic>, <italic>Chthonolasius</italic>, and <italic>Lasius pallitarsis</italic>, and the second lineage comprising <italic>Dendrolasius</italic>, and <italic>Lasius sensu stricto</italic>, with <italic>Cautolasius </italic>probably belonging to the second lineage (Fig. ##FIG##1##2##). Within the first lineage <italic>Acanthomyops </italic>and <italic>Austrolasius </italic>form a crown-group which is sister to <italic>Chthonolasius</italic>.</p> <title>Comparison with previous phylogenies</title> <p>The new framework disagrees with previous topologies [##UREF##28##36##, ####UREF##29##37##, ##REF##15522790##38####15522790##38##]. Given that we actively sought to exclude potential factors causing phylogenetic distortion from our data, and that the previous phylogenies disagreed with each other (Fig. ##FIG##0##1##), we suggest that where our topology disagrees with previous ones the new one is preferable (Table ##TAB##3##4##). Examination of the reasons for disagreements between the new topology and previous schemes is desirable and we now address this.</p> <p>The phylogeny by Wilson published in 1955 [##UREF##28##36##] (W55) is the only one of the previous phylogenies that hypothesised the existence of two major lineages as recovered in the new framework. In terms of subgenera, the two W55 lineages agree with those of the new framework with the exception of <italic>Lasius pallitarsis </italic>(treated as <italic>L. sitkaensis </italic>and believed to belong to <italic>Lasius sensu stricto </italic>at the time [##UREF##28##36##]). However, there are disagreements within the two lineages. Within the first lineage, the situation is ambiguous because <italic>Austrolasius</italic>, found in the new framework to be sister to <italic>Acanthomyops</italic>, was treated as part of <italic>Chthonolasius</italic>. Within the second lineage, the relations differ in that <italic>Lasius sensu stricto </italic>is considered sister to <italic>Cautolasius </italic>in W55, whereas it is sister to <italic>Dendrolasius </italic>in the new framework. Reasons for the disagreements could include that W55 is based on morphological information only, and that it lacked a formal reconstruction algorithm.</p> <p>In the phylogeny by Hasegawa of 1998 [##UREF##29##37##] (H98), there was only one significantly supported (i.e., for MP, node support &gt; 70 [##UREF##43##67##]) subgenus relationship, i.e., that <italic>Cautolasius</italic>, <italic>Chthonolasius</italic>, and <italic>Lasius sensu stricto </italic>form a crown-group sister to <italic>Dendrolasius</italic>, and this is not supported by the new framework. Reasons for the disagreement could include that H98 had very limited taxon sampling, and that no check for compositional heterogeneity was undertaken: The reconstruction methods applied (Neighbour Joining and Maximum Parsimony) could also contribute as the genetic distance based Neighbour Joining reconstruction is known to have the limitation that rate variation among sites cannot be accurately accounted for [##REF##11675604##68##], and as MP ignores the possible existence of a range of alternative topologies that are not significantly less likely than the most parsimonious one, even though they would require more evolutionary changes [##REF##16701310##40##].</p> <p>The subgenus relationships recovered by Janda and coworkers in 2004 [##REF##15522790##38##] (J04), placing <italic>Lasius sensu stricto </italic>in a sister clade to the rest, and within this clade <italic>Lasius pallitarsis </italic>sister to the rest, are not supported by the new framework. However, our data confirm two important aspects of J04, the allocation of <italic>Acanthomyops </italic>to the genus <italic>Lasius</italic>, and the discovery that <italic>Lasius pallitarsis </italic>does not belong to <italic>Lasius sensu stricto </italic>or any other of the established subgenera. Reasons for the disagreements with the new framework could include that the DNA data of J04 add up to a total of only 568 bp and that no check for compositional heterogeneity was undertaken. Moreover, while the combination of molecular with other data by Janda et al. marks an advance in the history of <italic>Lasius </italic>phylogeny, the morphological data included 33 characters (Additional File ##SUPPL##0##1##) that, as indicated by other studies [##UREF##0##1##,##UREF##1##2##,##UREF##9##12##,##REF##10486008##51##,##REF##16716603##52##,##REF##17366134##60##,##UREF##38##61##], may possibly be coupled functionally with social parasitism. Scrutiny in morphological character selection was recently confirmed as crucial in phylogenetic reconstruction [##REF##17366142##69##]. Also, the occurrence of social parasitism itself was included as character. These characters might contribute to grouping all parasitic subgenera together in an unresolved crown group as indeed is the case in J04 (additionally including <italic>Cautolasius</italic>), whereas they were distributed over the two major lineages in the new framework. The phylogenetic reconstruction methods applied (Maximum Parsimony) could possibly have contributed to suboptimal reconstruction as compared to Bayesian inference, especially for combined data sets, demonstrated after the publication of Janda et al. [##REF##14965900##39##]. Possibly none of these various causes had a strong effect <italic>per se</italic>, as suggested by rather minor differences between the Bayesian trees based on the concatenated data excluding <italic>versus </italic>including those morphological characters suspected to be coupled with social parasitism as well as between the Bayesian and the MP trees when excluding those morphological characters (Fig. ##FIG##1##2##). The various causes may, however, have added up to a significant effect, as is suggested by the lack of significant support for the sister-group relationship in the MP tree when including those morphological characters (0.68), which node scored a value of 1.0 in the Bayesian tree independently of excluding or including the characters.</p> <title>Hypotheses on the evolution of social parasitism and fungiculture</title> <p>The new phylogenetic framework ("new" topology in Fig. ##FIG##0##1##) and the reconstructed ancestral states (Table ##TAB##4##5##) suggest that both social parasitism and fungiculture evolved two times independently within <italic>Lasius</italic>, once in each of the major lineages (Fig. ##FIG##2##3##). For social parasitism such parallel evolution within a monophyletic group is not surprising in general. It has long been established that the combination of certain traits resulted in a predisposition for social parasitism in two (Myrmicinae and Formicinae) of the 21 [##UREF##44##70##] ant subfamilies [##UREF##1##2##,##UREF##11##15##,##UREF##45##71##]. Evidence that social parasitism also evolved multiple times within tribes [##UREF##12##16##,##UREF##18##24##,##REF##8529270##72##] or within genera [##UREF##11##15##,##UREF##16##21##,##REF##12740442##22##,##UREF##19##25##,##UREF##21##27##] confirms the principle for lower taxonomic levels. Having a reliable phylogenetic framework for <italic>Lasius </italic>facilitates examination of temporary social parasitism, considered less derived than other types of social parasitism [##REF##21227804##13##,##UREF##11##15##,##UREF##18##24##]. One factor that may have favoured the rise of parasitism in <italic>Lasius </italic>may derive from their colony organisation: In several species of the non-parasitic subgenera, <italic>Lasius sensu stricto </italic>and <italic>Cautolasius</italic>, the lack of aggression between different single-queened colonies has been reported [##REF##17964165##8##,##UREF##13##17##,##UREF##46##73##,##UREF##47##74##]. Whereas this behaviour concerns intraspecific interactions, the predisposition to reduce aggression may have been important for social parasitism to arise. The paramount significance of reduced aggression is illustrated by the chemical disguise of founding queens of <italic>Chthonolasius </italic>to appease host workers [##UREF##0##1##,##UREF##1##2##] and the exceptional winter activity of <italic>L</italic>. (<italic>Ch</italic>.) <italic>mixtus </italic>to ease entry to the then less aggressive host colonies [##UREF##0##1##]. Characterisation of such potentially preadaptive traits and search for them in extant non-parasitic <italic>Lasius </italic>species might help in finding potential early stages of incipient social parasitism. Such discoveries would then contribute to resolving mechanisms in the evolution of social parasitism in general [##UREF##10##14##, ####UREF##11##15##, ##UREF##12##16##, ##UREF##13##17####13##17##]. Fungiculture on the other hand only is known from one other group of ants, the Attini (Myrmicinae). It is established that fungiculture arose only once in attines [##REF##18362345##35##,##REF##11409051##75##] and the <italic>Lasius </italic>situation might thus indicate a stronger predisposition to evolve fungiculture for ants generally.</p> <p>For the evolution of both fungiculture and social parasites two pairs of alternative hypotheses remain (Fig. ##FIG##2##3##), both concerning the lineage of (((<italic>Lasius pallitarsis </italic>+ ((<italic>Chthonolasius </italic>+ (<italic>Austrolasius </italic>+ <italic>Acanthomyops</italic>))). This is because for both traits the probability value for the ancestral state for one node was below 0.95 (Table ##TAB##4##5##). For each trait there is one scenario considering an earlier origin with subsequent loss in the respective lineage, and the other considering an origin at a later stage in the respective lineage. Concerning both traits the hypothesis of a reversal is less likely from the point of view of biological realism because it is less parsimonious, and also because the probability values for the respective nodes (2 and 3) indicate, though not significantly, that social parasitism and fungiculture had not yet evolved at those stages (0.70 and 0.82). Still, we cannot formally refute the reversal scenarios here, given the lack of statistical support, and briefly discuss them in the below.</p> <p>For social parasitism, a reversal has been considered generally unlikely [##UREF##1##2##,##UREF##13##17##]. On the other hand, reversal from temporary social parasitism may be more likely than reversal from any other, more derived type of social parasitism. For example, temporarily socially parasitic ant species of <italic>Formica sensu stricto </italic>that are capable of alternative nest foundation through colony budding or adoption into another colony show a certain flexibility in the nest foundation mode in this type of social parasitism [##UREF##0##1##]. On the other hand, the fact that <italic>Acanthomyops</italic>, <italic>Austrolasius</italic>, and <italic>Chthonolasius </italic>form an entirely socially parasitic clade shows that temporary social parasitism may in fact be a stable trait and makes reversal appear less likely.</p> <p>Once arisen, fungiculture is not known to have been lost in insects and loss has been regarded as improbable [##UREF##22##28##], which provides an additional argument against the scenario involving a loss in <italic>Lasius</italic>. It may, however, be more probable in the <italic>Lasius </italic>case: fungiculture for nest building is probably less likely to entail a high dependence of the ants as compared to fungiculture for nutrition, which involves physiological adaptations [##REF##11409051##75##].</p> <p>It is not possible to definitely decide in favour of any of the alternative hypotheses at present, neither concerning social parasitism nor fungiculture, but it is worth considering the question further because a plausible loss of either would be of considerable evolutionary significance. For both traits the fossil record may offer answers in combination with a sound molecular clock for <italic>Lasius</italic>, for social parasitism by revealing the queen morphology of the ancestor of the entire major lineage, for fungiculture by a preserved fungal nest structure. Laboratory experiments may also be helpful for the analysis of both traits, for social parasitism by checking whether any of the temporary social parasites of the lineage are capable of independent colony foundation, and for fungiculture whether <italic>Chthonolasius </italic>colonies can be maintained when deprived of their fungus.</p> <p>This study provides a basis for studying the evolution of the various ecologically and sociobiologically relevant traits across the ca. 100 [##UREF##44##70##] species of <italic>Lasius</italic>, by establishing a phylogenetic framework and resolving the position of six of the seven lineages at subgenus level. We have used the framework to define hypotheses of the evolution of two outstanding traits, social parasitism and fungiculture, the evolution of which continues to pose riddles to evolutionary biology. Our results suggest that both traits arose twice in <italic>Lasius </italic>which opens up new opportunities for comparative analyses in a close phylogenetic relationship. We present competitive hypotheses that either do or do not involve reversal from the traits.</p>
<title>Discussion and Conclusion</title> <title>The new phylogenetic framework</title> <p>The new phylogenetic framework we present clarifies the relationships of all but one subgenus, namely the position of <italic>Cautolasius </italic>relative to the other subgenera: MP reconstructions significantly support the topology, <italic>Cautolasius </italic>+ (<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>), and Bayesian reconstructions do not contradict it. Given the lack of significant support in the Bayesian tree, we nevertheless subsequently refrain from considering the node as resolved ("new" topology in Fig. ##FIG##0##1##). This lack of significant resolution does not affect considerations on the evolution of social parasitism and fungiculture, though, because the ancestral state of <italic>Lasius </italic>was absence of both traits (Table ##TAB##4##5##), and because also <italic>Cautolasius </italic>displays neither trait. A significantly supported sister-status to (<italic>Dendrolasius </italic>+ <italic>Lasius sensu stricto</italic>) thus would not alter any conclusions (see discussion below). The previous phylogenies were recovered by various methods and we discuss below how the methods applied may have influenced the respective results. Taken together, the reasons for the increased robustness that we postulate for the new phylogenetic framework, which robustness is highlighted also by the confirmation of our Bayesian topology by our MP reconstructions, may be that (1) we excluded potential causes of phylogenetic distortion and (2) used additional analysis methods, that (3) the phylogenetic signal of combined data is potentially stronger than that of individual data sets [##REF##16674574##18##,##REF##14965900##39##,##UREF##30##41##, ####REF##17079492##42##, ##REF##9479692##43##, ##UREF##31##44##, ##UREF##32##45##, ##REF##16351964##46##, ##REF##16631389##47##, ##UREF##33##48##, ##REF##16601190##49##, ##REF##17408491##50####17408491##50##], especially when applying Bayesian inference [##REF##14965900##39##], that (4) all nodes but one at subgenus level or above were significantly supported in our concatenated analysis (Fig. ##FIG##1##2##), and that (5) we did not recover any disagreement between the individual data set phylogenies and the concatenated data topologies. Nevertheless, because the molecular data of this paper all derive from mitochondrial DNA, we cannot absolutely exclude that the tree is influenced by introgression at very shallow levels. Nuclear pseudogenes are also a concern, but we consider it unlikely that we amplified these because neither reading frame shifts nor sequence ambiguities were apparent. The possibilities of incomplete lineage sorting and of selection driven by symbionts which are in disequilibrium with mtDNA [##UREF##41##65##] cannot be ruled out, but such effects are unlikely to confound clade history at a deeper level, such as that of subgenera. It remains true that future studies using multiple nuclear genes are desirable to confirm our findings.</p> <p>The new phylogenetic framework presented here confirms the monophyly of the six subgenera. The new framework also affirms the taxonomic placement of <italic>Acanthomyops </italic>as a subgenus of <italic>Lasius </italic>[##UREF##42##66##]. Moreover, there is additional evidence for treating <italic>Lasius pallitarsis </italic>as a separate, monotypic subgenus, as suggested earlier [##REF##15522790##38##]. Detailed morphological characterisation of the new subgenus and taxonomic implications will be followed up elsewhere in the frame of a formal taxonomic revision.</p> <p>The evolution of the <italic>Lasius </italic>subgenera occurred in two major lineages, the first lineage comprising <italic>Acanthomyops</italic>, <italic>Austrolasius</italic>, <italic>Chthonolasius</italic>, and <italic>Lasius pallitarsis</italic>, and the second lineage comprising <italic>Dendrolasius</italic>, and <italic>Lasius sensu stricto</italic>, with <italic>Cautolasius </italic>probably belonging to the second lineage (Fig. ##FIG##1##2##). Within the first lineage <italic>Acanthomyops </italic>and <italic>Austrolasius </italic>form a crown-group which is sister to <italic>Chthonolasius</italic>.</p> <title>Comparison with previous phylogenies</title> <p>The new framework disagrees with previous topologies [##UREF##28##36##, ####UREF##29##37##, ##REF##15522790##38####15522790##38##]. Given that we actively sought to exclude potential factors causing phylogenetic distortion from our data, and that the previous phylogenies disagreed with each other (Fig. ##FIG##0##1##), we suggest that where our topology disagrees with previous ones the new one is preferable (Table ##TAB##3##4##). Examination of the reasons for disagreements between the new topology and previous schemes is desirable and we now address this.</p> <p>The phylogeny by Wilson published in 1955 [##UREF##28##36##] (W55) is the only one of the previous phylogenies that hypothesised the existence of two major lineages as recovered in the new framework. In terms of subgenera, the two W55 lineages agree with those of the new framework with the exception of <italic>Lasius pallitarsis </italic>(treated as <italic>L. sitkaensis </italic>and believed to belong to <italic>Lasius sensu stricto </italic>at the time [##UREF##28##36##]). However, there are disagreements within the two lineages. Within the first lineage, the situation is ambiguous because <italic>Austrolasius</italic>, found in the new framework to be sister to <italic>Acanthomyops</italic>, was treated as part of <italic>Chthonolasius</italic>. Within the second lineage, the relations differ in that <italic>Lasius sensu stricto </italic>is considered sister to <italic>Cautolasius </italic>in W55, whereas it is sister to <italic>Dendrolasius </italic>in the new framework. Reasons for the disagreements could include that W55 is based on morphological information only, and that it lacked a formal reconstruction algorithm.</p> <p>In the phylogeny by Hasegawa of 1998 [##UREF##29##37##] (H98), there was only one significantly supported (i.e., for MP, node support &gt; 70 [##UREF##43##67##]) subgenus relationship, i.e., that <italic>Cautolasius</italic>, <italic>Chthonolasius</italic>, and <italic>Lasius sensu stricto </italic>form a crown-group sister to <italic>Dendrolasius</italic>, and this is not supported by the new framework. Reasons for the disagreement could include that H98 had very limited taxon sampling, and that no check for compositional heterogeneity was undertaken: The reconstruction methods applied (Neighbour Joining and Maximum Parsimony) could also contribute as the genetic distance based Neighbour Joining reconstruction is known to have the limitation that rate variation among sites cannot be accurately accounted for [##REF##11675604##68##], and as MP ignores the possible existence of a range of alternative topologies that are not significantly less likely than the most parsimonious one, even though they would require more evolutionary changes [##REF##16701310##40##].</p> <p>The subgenus relationships recovered by Janda and coworkers in 2004 [##REF##15522790##38##] (J04), placing <italic>Lasius sensu stricto </italic>in a sister clade to the rest, and within this clade <italic>Lasius pallitarsis </italic>sister to the rest, are not supported by the new framework. However, our data confirm two important aspects of J04, the allocation of <italic>Acanthomyops </italic>to the genus <italic>Lasius</italic>, and the discovery that <italic>Lasius pallitarsis </italic>does not belong to <italic>Lasius sensu stricto </italic>or any other of the established subgenera. Reasons for the disagreements with the new framework could include that the DNA data of J04 add up to a total of only 568 bp and that no check for compositional heterogeneity was undertaken. Moreover, while the combination of molecular with other data by Janda et al. marks an advance in the history of <italic>Lasius </italic>phylogeny, the morphological data included 33 characters (Additional File ##SUPPL##0##1##) that, as indicated by other studies [##UREF##0##1##,##UREF##1##2##,##UREF##9##12##,##REF##10486008##51##,##REF##16716603##52##,##REF##17366134##60##,##UREF##38##61##], may possibly be coupled functionally with social parasitism. Scrutiny in morphological character selection was recently confirmed as crucial in phylogenetic reconstruction [##REF##17366142##69##]. Also, the occurrence of social parasitism itself was included as character. These characters might contribute to grouping all parasitic subgenera together in an unresolved crown group as indeed is the case in J04 (additionally including <italic>Cautolasius</italic>), whereas they were distributed over the two major lineages in the new framework. The phylogenetic reconstruction methods applied (Maximum Parsimony) could possibly have contributed to suboptimal reconstruction as compared to Bayesian inference, especially for combined data sets, demonstrated after the publication of Janda et al. [##REF##14965900##39##]. Possibly none of these various causes had a strong effect <italic>per se</italic>, as suggested by rather minor differences between the Bayesian trees based on the concatenated data excluding <italic>versus </italic>including those morphological characters suspected to be coupled with social parasitism as well as between the Bayesian and the MP trees when excluding those morphological characters (Fig. ##FIG##1##2##). The various causes may, however, have added up to a significant effect, as is suggested by the lack of significant support for the sister-group relationship in the MP tree when including those morphological characters (0.68), which node scored a value of 1.0 in the Bayesian tree independently of excluding or including the characters.</p> <title>Hypotheses on the evolution of social parasitism and fungiculture</title> <p>The new phylogenetic framework ("new" topology in Fig. ##FIG##0##1##) and the reconstructed ancestral states (Table ##TAB##4##5##) suggest that both social parasitism and fungiculture evolved two times independently within <italic>Lasius</italic>, once in each of the major lineages (Fig. ##FIG##2##3##). For social parasitism such parallel evolution within a monophyletic group is not surprising in general. It has long been established that the combination of certain traits resulted in a predisposition for social parasitism in two (Myrmicinae and Formicinae) of the 21 [##UREF##44##70##] ant subfamilies [##UREF##1##2##,##UREF##11##15##,##UREF##45##71##]. Evidence that social parasitism also evolved multiple times within tribes [##UREF##12##16##,##UREF##18##24##,##REF##8529270##72##] or within genera [##UREF##11##15##,##UREF##16##21##,##REF##12740442##22##,##UREF##19##25##,##UREF##21##27##] confirms the principle for lower taxonomic levels. Having a reliable phylogenetic framework for <italic>Lasius </italic>facilitates examination of temporary social parasitism, considered less derived than other types of social parasitism [##REF##21227804##13##,##UREF##11##15##,##UREF##18##24##]. One factor that may have favoured the rise of parasitism in <italic>Lasius </italic>may derive from their colony organisation: In several species of the non-parasitic subgenera, <italic>Lasius sensu stricto </italic>and <italic>Cautolasius</italic>, the lack of aggression between different single-queened colonies has been reported [##REF##17964165##8##,##UREF##13##17##,##UREF##46##73##,##UREF##47##74##]. Whereas this behaviour concerns intraspecific interactions, the predisposition to reduce aggression may have been important for social parasitism to arise. The paramount significance of reduced aggression is illustrated by the chemical disguise of founding queens of <italic>Chthonolasius </italic>to appease host workers [##UREF##0##1##,##UREF##1##2##] and the exceptional winter activity of <italic>L</italic>. (<italic>Ch</italic>.) <italic>mixtus </italic>to ease entry to the then less aggressive host colonies [##UREF##0##1##]. Characterisation of such potentially preadaptive traits and search for them in extant non-parasitic <italic>Lasius </italic>species might help in finding potential early stages of incipient social parasitism. Such discoveries would then contribute to resolving mechanisms in the evolution of social parasitism in general [##UREF##10##14##, ####UREF##11##15##, ##UREF##12##16##, ##UREF##13##17####13##17##]. Fungiculture on the other hand only is known from one other group of ants, the Attini (Myrmicinae). It is established that fungiculture arose only once in attines [##REF##18362345##35##,##REF##11409051##75##] and the <italic>Lasius </italic>situation might thus indicate a stronger predisposition to evolve fungiculture for ants generally.</p> <p>For the evolution of both fungiculture and social parasites two pairs of alternative hypotheses remain (Fig. ##FIG##2##3##), both concerning the lineage of (((<italic>Lasius pallitarsis </italic>+ ((<italic>Chthonolasius </italic>+ (<italic>Austrolasius </italic>+ <italic>Acanthomyops</italic>))). This is because for both traits the probability value for the ancestral state for one node was below 0.95 (Table ##TAB##4##5##). For each trait there is one scenario considering an earlier origin with subsequent loss in the respective lineage, and the other considering an origin at a later stage in the respective lineage. Concerning both traits the hypothesis of a reversal is less likely from the point of view of biological realism because it is less parsimonious, and also because the probability values for the respective nodes (2 and 3) indicate, though not significantly, that social parasitism and fungiculture had not yet evolved at those stages (0.70 and 0.82). Still, we cannot formally refute the reversal scenarios here, given the lack of statistical support, and briefly discuss them in the below.</p> <p>For social parasitism, a reversal has been considered generally unlikely [##UREF##1##2##,##UREF##13##17##]. On the other hand, reversal from temporary social parasitism may be more likely than reversal from any other, more derived type of social parasitism. For example, temporarily socially parasitic ant species of <italic>Formica sensu stricto </italic>that are capable of alternative nest foundation through colony budding or adoption into another colony show a certain flexibility in the nest foundation mode in this type of social parasitism [##UREF##0##1##]. On the other hand, the fact that <italic>Acanthomyops</italic>, <italic>Austrolasius</italic>, and <italic>Chthonolasius </italic>form an entirely socially parasitic clade shows that temporary social parasitism may in fact be a stable trait and makes reversal appear less likely.</p> <p>Once arisen, fungiculture is not known to have been lost in insects and loss has been regarded as improbable [##UREF##22##28##], which provides an additional argument against the scenario involving a loss in <italic>Lasius</italic>. It may, however, be more probable in the <italic>Lasius </italic>case: fungiculture for nest building is probably less likely to entail a high dependence of the ants as compared to fungiculture for nutrition, which involves physiological adaptations [##REF##11409051##75##].</p> <p>It is not possible to definitely decide in favour of any of the alternative hypotheses at present, neither concerning social parasitism nor fungiculture, but it is worth considering the question further because a plausible loss of either would be of considerable evolutionary significance. For both traits the fossil record may offer answers in combination with a sound molecular clock for <italic>Lasius</italic>, for social parasitism by revealing the queen morphology of the ancestor of the entire major lineage, for fungiculture by a preserved fungal nest structure. Laboratory experiments may also be helpful for the analysis of both traits, for social parasitism by checking whether any of the temporary social parasites of the lineage are capable of independent colony foundation, and for fungiculture whether <italic>Chthonolasius </italic>colonies can be maintained when deprived of their fungus.</p> <p>This study provides a basis for studying the evolution of the various ecologically and sociobiologically relevant traits across the ca. 100 [##UREF##44##70##] species of <italic>Lasius</italic>, by establishing a phylogenetic framework and resolving the position of six of the seven lineages at subgenus level. We have used the framework to define hypotheses of the evolution of two outstanding traits, social parasitism and fungiculture, the evolution of which continues to pose riddles to evolutionary biology. Our results suggest that both traits arose twice in <italic>Lasius </italic>which opens up new opportunities for comparative analyses in a close phylogenetic relationship. We present competitive hypotheses that either do or do not involve reversal from the traits.</p>
<title>Background</title> <p>Ants of the genus <italic>Lasius </italic>are ecologically important and an important system for evolutionary research. Progress in evolutionary research has been hindered by the lack of a well-founded phylogeny of the subgenera, with three previous attempts disagreeing. Here we employed two mitochondrial genes (<italic>cytochrome c oxidase subunit I, 16S ribosomal RNA</italic>), comprising 1,265 bp, together with 64 morphological characters, to recover the phylogeny of <italic>Lasius </italic>by Bayesian and Maximum Parsimony inference after exploration of potential causes of phylogenetic distortion. We use the resulting framework to infer evolutionary pathways for social parasitism and fungiculture.</p> <title>Results</title> <p>We recovered two well supported major lineages. One includes <italic>Acanthomyops</italic>, <italic>Austrolasius</italic>, <italic>Chthonolasius</italic>, and <italic>Lasius pallitarsis</italic>, which we confirm to represent a seventh subgenus, the other clade contains <italic>Dendrolasius</italic>, and <italic>Lasius sensu stricto</italic>. The subgenus <italic>Cautolasius</italic>, displaying neither social parasitism nor fungiculture, probably belongs to the second clade, but its phylogenetic position is not resolved at the cutoff values of node support we apply. Possible causes for previous problems with reconstructing the <italic>Lasius </italic>phylogeny include use of other reconstruction techniques, possibly more prone to instabilities in some instances, and the inclusion of phylogenetically distorting characters.</p> <title>Conclusion</title> <p>By establishing an updated phylogenetic framework, our study provides the basis for a later formal taxonomic revision of subgenera and for studying the evolution of various ecologically and sociobiologically relevant traits of <italic>Lasius</italic>, although there is need for future studies to include nuclear genes and additional samples from the Nearctic. Both social parasitism and fungiculture evolved twice in <italic>Lasius</italic>, once in each major lineage, which opens up new opportunities for comparative analyses. The repeated evolution of social parasitism has been established for other groups of ants, though not for temporary social parasitism as found in <italic>Lasius</italic>. For fungiculture, the independent emergence twice in a monophyletic group marks a novel scenario in ants. We present alternative hypotheses for the evolution of both traits, with one of each involving loss of the trait. Though less likely for both traits than later evolution without reversal, we consider reversal as sufficiently plausible to merit independent testing.</p>
<title>Authors' contributions</title> <p>MM initiated the study, participated in its design and coordination, carried out the DNA sequencing work, performed the sequence alignment, and sampled the morphological characters. FMS participated in the design and coordination of the study, in the data analysis, in drafting a first version of the manuscript, and in revising it. TA helped with the design of the study. CS helped with the design of the study. RHC helped to design the data analysis and the format of the manuscript. BCS–S participated in the design and coordination of the study, in the data analysis, in drafting a first version of the manuscript, and in revising it. All authors read and approved the final version.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank A. Françoeur, N. Fujiwara, K. Hamaguti, Y. Ikeshita, T. Itino, F. Ito, Y. Kida, K. Kinomura, T. Kishimoto, T. Komatsu, S. Nagashima, K. Ogata, M. Ôhara, J.M. Raczkowski, C.A. Schmidt, T. Shimada, T. Toida for providing material. This work was performed at the molecular laboratory of Division of Mammals and Birds, National Science Museum, Tokyo, under supervision of I. Nishiumi, S. Hoshino and other staffs. For valuable information, we thank J. Felsenstein, K. Ogata, B. Seifert, and P.S. Ward. This project was partially supported by Grant-in-Aid of JSPS Postdoctoral Fellowship to MM. BCS and FMS were supported by the Austrian Science Fund [P-17219-B06, J2639-B17, J2642-B17]. RHC's work on evolutionary genetics is supported by the Australian Research Council [DP0665890]. We are also grateful to five anonymous referees and the Assistant Editor, E. Alexandersson, for constructive and inspiring suggestions.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Subgenus relationships in the previous and new phylogenetic reconstructions of the genus <italic>Lasius</italic></bold>. Subgenera abbreviations: <italic>Acanthomyops </italic>(<italic>Ac</italic>), <italic>Austrolasius </italic>(<italic>Au</italic>), <italic>Cautolasius </italic>(<italic>Ca</italic>), <italic>Chthonolasius </italic>(<italic>Ch</italic>), <italic>Dendrolasius </italic>(<italic>D</italic>), <italic>Lasius sensu stricto </italic>(<italic>L</italic>). The topologies were extracted from papers by Wilson [##UREF##28##36##], "W55", Hasegawa [##UREF##29##37##], "H98", and Janda and coworkers [##REF##15522790##38##], "J04", as well as from the Bayesian tree of combined, concatenated data in Fig. 2 of this paper, "new"; see Methods section for details of the procedure used for inferring the topologies W55, H98, and J04. A dotted line indicates that node support for monophyly of the subgenus was not significant. White squares indicate constraints enforced in constraint analyses using our concatenated data set in order to test the subgenus relationships of W55, H98, and J04.</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Bayesian topology from the analysis of the combined, concatenated data</bold>. Subgenera are abbreviated as in Fig. ##FIG##0##1##. The tree is a consensus tree resulting from a Bayesian analysis of our concatenated data set based on <italic>cox1 </italic>plus 16S rRNA plus morphology. The credibility values are posterior probabilities derived from 20,000 trees representing 2 million generations after burnin (upper left), bootstrap values from the 50% majority-rule consensus MP tree of the same data (lower right, in quotation marks); values for nodes following the basal divergence within subgenera are omitted. The node numbers refer to the inferred ancestral character states given in Table 5. The inset tree on grey background is a Bayesian tree based on <italic>cox1 </italic>plus 16S rRNA and the credibility values are posterior probabilities derived from 20,000 trees representing 2 million generations after burnin.</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Hypotheses on the evolution of social parasitism and fungiculture in <italic>Lasius </italic>ants</bold>. Subgenera are abbreviated as in Fig. 1. Asterisks indicate the emergence of the trait, crossed circles its loss. Subgenera currently displaying social parasitism or fungiculture are indicated by frames filled black. Alternative hypotheses are offered for the evolution of social parasitism and fungiculture due to two insignificant results of the ancestral character state reconstruction in Table 5, with probabilites for the competing scenarios given, based on the values in Table 5.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>List of samples used for DNA sequencing and morphological analysis</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Species</th><th align="left">Subgenus</th><th align="left">Collection locality; collector</th><th align="center" colspan="2">DDBJ accession numbers</th><th align="left">Museum voucher no</th></tr><tr><th/><th/><th/><th colspan="2"><hr/></th><th/></tr><tr><th/><th/><th/><th align="left"><italic>cox1</italic></th><th align="left"><italic>16S rRNA</italic></th><th/></tr></thead><tbody><tr><td align="left"><italic>Lasius arizonicus</italic></td><td align="left"><italic>Ac.</italic></td><td align="left">USA: Arizona, Madera Canyon; C.A. Schmidt</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370982">AB370982</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371028">AB371028</ext-link></td><td align="left">MMANT12</td></tr><tr><td align="left"><italic>L. interjectus</italic></td><td align="left"><italic>Ac.</italic></td><td align="left">USA: Arizona, West Turkey Creek; C.A. Schmidt</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370981">AB370981</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371027">AB371027</ext-link></td><td align="left">MMANT13</td></tr><tr><td align="left"><italic>L. latipes</italic></td><td align="left"><italic>Ac.</italic></td><td align="left">USA: Wisconsin, Milwaukee; J.M. Raczkowski</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB433922">AB433922</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB433927">AB433927</ext-link></td><td align="left">NMANT120</td></tr><tr><td align="left"><italic>L. reginae</italic></td><td align="left"><italic>Au.</italic></td><td align="left">Austria: Trandorf; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370983">AB370983</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371029">AB371029</ext-link></td><td align="left">MMANT23</td></tr><tr><td align="left"><italic>L. flavus</italic></td><td align="left"><italic>Ca.</italic></td><td align="left">Austria: Leiser Berge; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370984">AB370984</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371030">AB371030</ext-link></td><td align="left">MMANT22</td></tr><tr><td align="left"><italic>L. flavus</italic></td><td align="left"><italic>Ca.</italic></td><td align="left">Russia: Ussurisky, Kaimanovka; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370985">AB370985</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371031">AB371031</ext-link></td><td align="left">MMANT38</td></tr><tr><td align="left"><italic>L. flavus</italic></td><td align="left"><italic>Ca.</italic></td><td align="left">Japan: Gifu-ken, Takayama-shi, M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370986">AB370986</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371032">AB371032</ext-link></td><td align="left">MMANT45</td></tr><tr><td align="left"><italic>L. nearcticus</italic></td><td align="left"><italic>Ca.</italic></td><td align="left">USA: Arizona, Rustler Park; C.A. Schmidt</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370987">AB370987</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371033">AB371033</ext-link></td><td align="left">MMANT14</td></tr><tr><td align="left"><italic>L. mixtus</italic></td><td align="left"><italic>Ch.</italic></td><td align="left">Austria: Göpfritz; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370988">AB370988</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371034">AB371034</ext-link></td><td align="left">MMANT30</td></tr><tr><td align="left"><italic>L. umbratus</italic></td><td align="left"><italic>Ch.</italic></td><td align="left">Japan: Tôkyô-to, Koganei-shi; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370989">AB370989</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371035">AB371035</ext-link></td><td align="left">MMANT6</td></tr><tr><td align="left"><italic>L. capitatus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Nagano-ken, Matsumoto-shi; T. Komatsu</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370990">AB370990</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371036">AB371036</ext-link></td><td align="left">MMANT44</td></tr><tr><td align="left"><italic>L. capitatus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Gifu-ken, Shôkawa-mura; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370993">AB370993</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371039">AB371039</ext-link></td><td align="left">MMANT47</td></tr><tr><td align="left"><italic>L. capitatus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Yamanashi-ken, Kitakoma-gun; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370991">AB370991</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371037">AB371037</ext-link></td><td align="left">MMANT58</td></tr><tr><td align="left"><italic>L. capitatus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Tochigi-ken, Haga-gun; S. Nagashima</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370992">AB370992</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371038">AB371038</ext-link></td><td align="left">MMANT62</td></tr><tr><td align="left"><italic>L. fuji</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Hokkaidô, Maruseppu-chô; Y. Kida</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370994">AB370994</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371040">AB371040</ext-link></td><td align="left">MMANT1</td></tr><tr><td align="left"><italic>L. fuji</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Russia: Ussurisky, Kaimanovka; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370995">AB370995</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371041">AB371041</ext-link></td><td align="left">MMANT34</td></tr><tr><td align="left"><italic>L. fuliginosus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Austria: Urschendorf; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370996">AB370996</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371042">AB371042</ext-link></td><td align="left">MMANT24</td></tr><tr><td align="left"><italic>L. fuliginosus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Austria: Vienna; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370997">AB370997</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371043">AB371043</ext-link></td><td align="left">MMANT70</td></tr><tr><td align="left"><italic>L. nipponensis</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Russia: Ussursky, Vityas; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371001">AB371001</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371047">AB371047</ext-link></td><td align="left">MMANT33</td></tr><tr><td align="left"><italic>L. nipponensis</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Hokkaidô, Sapporo-shi; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370998">AB370998</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371044">AB371044</ext-link></td><td align="left">MMANT63</td></tr><tr><td align="left"><italic>L. nipponensis</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Nagano-ken, Fujimi-chô; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB370999">AB370999</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371045">AB371045</ext-link></td><td align="left">MMANT64</td></tr><tr><td align="left"><italic>L. nipponensis</italic></td><td align="left"><italic>D</italic>.</td><td align="left">China: Hubei, Xianfeng; T. Kishimoto</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371000">AB371000</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371046">AB371046</ext-link></td><td align="left">MMANT67</td></tr><tr><td align="left"><italic>L. orientalis</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Hokkaidô, Shari-chô; Y. Kida</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371002">AB371002</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371048">AB371048</ext-link></td><td align="left">MMANT4</td></tr><tr><td align="left"><italic>L. orientalis</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Gifu-ken, Kamitakara-mura; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371003">AB371003</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371049">AB371049</ext-link></td><td align="left">MMANT60</td></tr><tr><td align="left"><italic>L. spathepus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Shimane-ken, Oki-shotô; T. Shimada</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371006">AB371006</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371052">AB371052</ext-link></td><td align="left">MMANT32</td></tr><tr><td align="left"><italic>L. spathepus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Yamanashi-ken, Nagasaka-chô; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371005">AB371005</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371051">AB371051</ext-link></td><td align="left">MMANT74</td></tr><tr><td align="left"><italic>L. spathepus</italic></td><td align="left"><italic>D</italic>.</td><td align="left">Japan: Kyôto-fu, Kyôto-shi; N. Fujiwara</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371007">AB371007</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371053">AB371053</ext-link></td><td align="left">MMANT77</td></tr><tr><td align="left"><italic>L. alienus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Austria: Braunsberg; B.C. Schlick-Steiner &amp; F.M.Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371008">AB371008</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371054">AB371054</ext-link></td><td align="left">MMANT21</td></tr><tr><td align="left"><italic>L. austriacus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Austria: Feldberg; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371009">AB371009</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371055">AB371055</ext-link></td><td align="left">MMANT27</td></tr><tr><td align="left"><italic>L. brunneus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Austria: Rassing; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371010">AB371010</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371056">AB371056</ext-link></td><td align="left">MMANT25</td></tr><tr><td align="left"><italic>L. emarginatus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Austria: Vienna; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371011">AB371011</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371057">AB371057</ext-link></td><td align="left">MMANT41</td></tr><tr><td align="left"><italic>L. hayashi</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Japan: Gifu-ken, Kamitakara-mura; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371013">AB371013</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371059">AB371059</ext-link></td><td align="left">MMANT46</td></tr><tr><td align="left"><italic>L. hayashi</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Japan: Chiba-ken, Kimitsu-shi; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371012">AB371012</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371058">AB371058</ext-link></td><td align="left">MMANT54</td></tr><tr><td align="left"><italic>L. japonicus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Japan: Kagawa-ken, Takamatsu-shi; F. Ito &amp; Y. Ikeshita</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371015">AB371015</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371061">AB371061</ext-link></td><td align="left">MMANT19</td></tr><tr><td align="left"><italic>L. japonicus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Russia: Ussurisky, Kaimanovka; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371017">AB371017</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371063">AB371063</ext-link></td><td align="left">MMANT37</td></tr><tr><td align="left"><italic>L. japonicus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Japan: Chiba-ken, Kimitsu-shi; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371014">AB371014</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371060">AB371060</ext-link></td><td align="left">MMANT55</td></tr><tr><td align="left"><italic>L. japonicus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Japan: Hokkaidô, Sapporo-shi; T. Toida</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371016">AB371016</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371062">AB371062</ext-link></td><td align="left">MMANT76</td></tr><tr><td align="left"><italic>L. neglectus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Hungary: Budapest; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371018">AB371018</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371064">AB371064</ext-link></td><td align="left">MMANT20</td></tr><tr><td align="left"><italic>L. niger</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Austria: Vienna; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371019">AB371019</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371065">AB371065</ext-link></td><td align="left">MMANT26</td></tr><tr><td align="left"><italic>L. platythorax</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Austria: Moosbrunn; B.C. Schlick-Steiner &amp; F.M. Steiner</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371020">AB371020</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371066">AB371066</ext-link></td><td align="left">MMANT28</td></tr><tr><td align="left"><italic>L. productus</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Japan: Kagawa-ken, Takamatsu-shi; F. Ito &amp; Y. Ikeshita</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371021">AB371021</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371067">AB371067</ext-link></td><td align="left">MMANT18</td></tr><tr><td align="left"><italic>L. sakagamii</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Japan: Gifu-ken, Gifu-shi; J. Heinze</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371022">AB371022</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371068">AB371068</ext-link></td><td align="left">MMANT29</td></tr><tr><td align="left"><italic>L. sakagamii</italic></td><td align="left"><italic>L</italic>.</td><td align="left">Japan: Tôkyô-to, Edogawa-ku; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371023">AB371023</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371069">AB371069</ext-link></td><td align="left">MMANT56</td></tr><tr><td align="left"><italic>L</italic>. sp.3</td><td align="left"><italic>L</italic>.</td><td align="left">Russia: Ussurisky, Kaimanovka; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371024">AB371024</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371070">AB371070</ext-link></td><td align="left">MMANT40</td></tr><tr><td align="left"><italic>L. pallitarsis</italic></td><td align="left"><italic>L pallitarsis</italic></td><td align="left">USA: Arizona, Apache Ntl Forest; C.A. Schmidt</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371025">AB371025</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371071">AB371071</ext-link></td><td align="left">MMANT15</td></tr><tr><td align="left"><italic>Myrmecocystus mimicus</italic></td><td align="left">n.a.</td><td align="left">USA: California, Carrizo Plain; P.S. Ward</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB433923">AB433923</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB433928">AB433928</ext-link></td><td align="left">MMANT117</td></tr><tr><td align="left"><italic>Myrmecocystus mendux</italic></td><td align="left">n.a.</td><td align="left">USA: Arizona, Pima Canyon; C.A. Schmidt</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB433920">AB433920</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB433925">AB433925</ext-link></td><td align="left">MMANT66</td></tr><tr><td align="left"><italic>Formica japonica</italic></td><td align="left">n.a.</td><td align="left">Japan: Tôkyô-to, Shinjuku-ku; M. Maruyama</td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371026">AB371026</ext-link></td><td align="left"><ext-link ext-link-type="gen" xlink:href="AB371072">AB371072</ext-link></td><td align="left">MMANT7</td></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>Character counts and substitution models for partitions</p></caption><table frame="hsides" rules="groups"><thead><tr><th/><th align="left">Characters total</th><th align="left">Characters variable but parsimony uninformative</th><th align="left">Characters parsimony informative</th><th align="left">AIC model selection</th><th align="left">hLRT model selection</th><th align="left">Model used in final MCMC runs</th></tr></thead><tbody><tr><td align="left"><italic>cox1 </italic>position1</td><td align="left">281</td><td align="left">17</td><td align="left">22</td><td align="left">GTR+I+Γ</td><td align="left">GTR+Γ</td><td align="left">GTR+I+Γ</td></tr><tr><td align="left"><italic>cox1 </italic>position2</td><td align="left">281</td><td align="left">8</td><td align="left">2</td><td align="left">F81</td><td align="left">F81</td><td align="left">F81</td></tr><tr><td align="left"><italic>cox1 </italic>position3 RY</td><td align="left">281</td><td align="left">28</td><td align="left">87</td><td align="left">-</td><td align="left">-</td><td align="left">F81</td></tr><tr><td align="left"><italic>16s rRNA</italic></td><td align="left">422</td><td align="left">45</td><td align="left">83</td><td align="left">GTR+I+Γ</td><td align="left">GTR+Γ</td><td align="left">GTR+I+Γ</td></tr><tr><td align="left">morphology</td><td align="left">64</td><td align="left">10</td><td align="left">54</td><td align="left">-</td><td align="left">-</td><td align="left">Mk+Γ</td></tr></tbody></table></table-wrap> <table-wrap id="T3" position="float"><label>Table 3</label><caption><p>Summary of results from Bayesian analyses</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Data</th><th align="left">ngens</th><th align="left">ln(Ar)</th><th align="left">ln(Hr)</th><th align="left">asdsf</th><th align="left">burnin</th><th align="left">99%</th></tr></thead><tbody><tr><td align="left">concatenated data: <italic>cox1</italic>_position 3 RY + <italic>16S rRNA </italic>+ morphology (Mk+Γ)</td><td align="left">10.0</td><td align="left">-5474.9</td><td align="left">-5533.3</td><td align="left">0.002</td><td align="left">9.0</td><td align="left">13369</td></tr><tr><td align="left">W55 constraint on concatenated data</td><td align="left">10.0</td><td align="left">-5593.4</td><td align="left">-5671.6</td><td align="left">0.002</td><td align="left">9.0</td><td align="left">11893</td></tr><tr><td align="left">H98 constraint on concatenated data</td><td align="left">10.0</td><td align="left">-5623.8</td><td align="left">-5717.0</td><td align="left">0.063</td><td align="left">9.0</td><td align="left">15960</td></tr><tr><td align="left">J04 constraint on concatenated data</td><td align="left">10.0</td><td align="left">-5576.0</td><td align="left">-5630.7</td><td align="left">0.002</td><td align="left">9.0</td><td align="left">13706</td></tr></tbody></table></table-wrap> <table-wrap id="T4" position="float"><label>Table 4</label><caption><p>Comparing previous Lasius phylogenies with the new phylogenetic framework</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Data</th><th align="left">Bayes factor</th><th align="left">Interpretation</th></tr></thead><tbody><tr><td align="left">concatenated vs. W55 constraint on concatenated</td><td align="left">276.7</td><td align="left">very strong evidence against W55 constraint</td></tr><tr><td align="left">concatenated vs. H98 constraint on concatenated</td><td align="left">367.4</td><td align="left">very strong evidence against H98 constraint</td></tr><tr><td align="left">concatenated vs. J04 constraint on concatenated</td><td align="left">194.9</td><td align="left">very strong evidence against J04 constraint</td></tr></tbody></table></table-wrap> <table-wrap id="T5" position="float"><label>Table 5</label><caption><p>Bayesian posterior probabilities for the occurrence of social parasitism and fungiculture as ancestral character states</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Node</th><th align="center" colspan="2">Social parasitism</th><th align="center" colspan="2">Fungiculture</th></tr><tr><th/><th colspan="2"><hr/></th><th colspan="2"><hr/></th></tr><tr><th/><th align="center">no</th><th align="center">yes</th><th align="center">no</th><th align="center">yes</th></tr></thead><tbody><tr><td align="center">1</td><td align="center"><bold>0.98</bold></td><td align="center">0.02</td><td align="center"><bold>1.00</bold></td><td align="center">0.00</td></tr><tr><td align="center">2</td><td align="center">0.70</td><td align="center">0.30</td><td align="center"><bold>0.98</bold></td><td align="center">0.02</td></tr><tr><td align="center">3</td><td align="center">0.00</td><td align="center"><bold>1.00</bold></td><td align="center">0.82</td><td align="center">0.18</td></tr><tr><td align="center">4</td><td align="center">0.00</td><td align="center"><bold>1.00</bold></td><td align="center"><bold>1.00</bold></td><td align="center">0.00</td></tr><tr><td align="center">5</td><td align="center"><bold>0.97</bold></td><td align="center">0.03</td><td align="center"><bold>0.98</bold></td><td align="center">0.02</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional File 1</title><p><bold>Definition and data matrix of morphological characters</bold>. L worker larva, M adult male, Q adult queen, W adult worker. The characters are consecutively numbered in the definitions which equal the numbers used in the matrix. The numeric codes for the character state definitions are given in parentheses. If a species is polymorphic in character states, both states are given in parantheses. In the penultimate row, the corresponding character name of Janda et al. [##REF##15522790##38##] is given and the following signs describe the relation to [##REF##15522790##38##] concerning the definitions, "*" indicates that character and states are identical, "**" indicates that the characters are identical but the states have been adapted, "***" that characters and states have been adapted, and "****" indicates new characters, not used in [##REF##15522790##38##]. The letters in ultimate row indicate: "a" character used for reconstruction of new phylogenetic framework, "b" character suspected to be functionally coupled with social parasitism or a behavioural/ecological character, "c" invariant character. The samples of the matrix are identical with those given in Table ##TAB##0##1##. Subgenera are abbreviated as in Table ##TAB##0##1##.</p></caption></supplementary-material>
<table-wrap-foot><p>Voucher specimens have been deposited at the National Science Museum, Tokyo, with the numbers indicated and a reference to this publication (Maruyama et al. 2008/voucher no). Abbreviations of <italic>Lasius </italic>subgenera are <italic>Acanthomyops </italic>(<italic>Ac</italic>), <italic>Austrolasius </italic>(<italic>Au</italic>), <italic>Cautolasius </italic>(<italic>Ca</italic>), <italic>Chthonolasius </italic>(<italic>Ch</italic>), <italic>Dendrolasius </italic>(<italic>D</italic>), <italic>Lasius sensu stricto </italic>(<italic>L</italic>). Subgenus was not applicable (n.a.) for the outgroup taxa, <italic>Myrmecocystus mimicus </italic>and <italic>M. mendux</italic>, and <italic>Formica japonica</italic>.</p></table-wrap-foot> <table-wrap-foot><p>RY purine + pyrimidine coding, otherwise the original nucleotide sequence was used. "-" under model selection indicates that AIC and hLRT model selection were not applicable for the partition.</p></table-wrap-foot> <table-wrap-foot><p>ngens (number of generations) and burnin are given in units of a million; Ar and Hr refer to the arithmetic and harmonic means, averaged over the simultaneous runs; asdsf = average standard deviation of split frequencies; 99% refers to the number of trees sampled from the 99% credible set.</p></table-wrap-foot> <table-wrap-foot><p>Summary of Bayes factor comparisons and interpretation after [##UREF##59##96##].</p></table-wrap-foot> <table-wrap-foot><p>"no" indicates absence, "yes" indicates presence; the posterior probabilities were estimated using SIMMAP and the last 2000 post-burnin trees of the 20,000 used to derive the consensus tree of Fig. 2. Significantly positive values (p &gt; 0.95) are given in bold. Node numbers refer to those shown in Fig. 2.</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-237-1"/> <graphic xlink:href="1471-2148-8-237-2"/> <graphic xlink:href="1471-2148-8-237-3"/>
<media xlink:href="1471-2148-8-237-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Seifert"], "given-names": ["B"], "source": ["Die Ameisen Mittel- und Nordeuropas"], "year": ["2007"], "publisher-name": ["Tauer: lutra Verlags- und Vertriebsgesellschaft"]}, {"surname": ["H\u00f6lldobler", "Wilson"], "given-names": ["B", "EO"], "source": ["The Ants"], "year": ["1990"], "publisher-name": ["Cambridge, MA: The Belknap Press of Harvard University Press"]}, {"surname": ["Akino", "Yamaoka"], "given-names": ["T", "R"], "article-title": ["Trail discrimination signal of "], "italic": ["Lasius japonicus "], "source": ["Chemoecology"], "year": ["2005"], "volume": ["15"], "fpage": ["21"], "lpage": ["30"], "pub-id": ["10.1007/s00049-005-0288-6"]}, {"surname": ["Kern", "Klein", "Janssen", "Bestmann", "Attygalle", "Sch\u00e4fer", "Maschwitz"], "given-names": ["F", "RW", "E", "H-J", "AB", "D", "U"], "article-title": ["Mellein, a trail pheromone component of the ant "], "italic": ["Lasius fuliginosus"], "source": ["J Chem Ecol"], "year": ["1997"], "volume": ["23"], "fpage": ["779"], "lpage": ["792"], "pub-id": ["10.1023/B:JOEC.0000006410.35938.49"]}, {"surname": ["Wu", "Mori"], "given-names": ["J", "K"], "article-title": ["Synthesis of the enantiomers of 2,6-dimethyl-6-helpten-1-ol, a mandibular gland secretion of the male ant "], "italic": ["Lasius niger"], "source": ["Agric Biol Chem"], "year": ["1991"], "volume": ["55"], "fpage": ["2667"], "lpage": ["2668"]}, {"surname": ["Akino", "Yamaoka"], "given-names": ["T", "R"], "article-title": ["Evidence for volatile and contact signals of nestmate recognition in the black shining ant "], "italic": ["Lasius fuliginosus "], "source": ["Entomol Sci"], "year": ["2000"], "volume": ["3"], "fpage": ["1"], "lpage": ["8"]}, {"surname": ["Morgan"], "given-names": ["ED"], "article-title": ["Chemical sorcery for sociality: Exocrine secretions of ants (Hymenoptera: Formicidae)"], "source": ["Myrmecol News"], "year": ["2008"], "volume": ["11"], "fpage": ["79"], "lpage": ["90"]}, {"surname": ["Bourke", "Franks"], "given-names": ["AFG", "NR"], "source": ["Social Evolution in Ants"], "year": ["1995"], "publisher-name": ["Princeton, NJ: Princeton University Press"]}, {"surname": ["Crozier", "Pamilo"], "given-names": ["RH", "P"], "source": ["Evolution of Social Insect Colonies"], "year": ["1996"], "publisher-name": ["New York: Oxford University Press"]}, {"surname": ["Boomsma", "Have"], "given-names": ["JJ", "TM van der"], "article-title": ["Queen mating and paternity variation in the ant "], "italic": ["Lasius niger"], "source": ["Mol Ecol"], "year": ["1998"], "volume": ["7"], "fpage": ["1709"], "lpage": ["1718"], "pub-id": ["10.1046/j.1365-294x.1998.00504.x"]}, {"surname": ["Bourke", "Franks"], "given-names": ["AFG", "NR"], "article-title": ["Alternative adaptations, sympatric speciation and the evolution of parasitic, inquiline ants"], "source": ["Biol J Linn Soc"], "year": ["1991"], "volume": ["43"], "fpage": ["157"], "lpage": ["178"], "pub-id": ["10.1111/j.1095-8312.1991.tb00591.x"]}, {"surname": ["Ward"], "given-names": ["PS"], "article-title": ["A new workerless social parasite in the ant genus "], "italic": ["Pseudomyrmex "], "source": ["Syst Entomol"], "year": ["1996"], "volume": ["21"], "fpage": ["253"], "lpage": ["263"], "pub-id": ["10.1046/j.1365-3113.1996.d01-12.x"]}, {"surname": ["Buschinger"], "given-names": ["A"], "article-title": ["Sympatric speciation and radiative evolution of socially parasitic ants \u2013 heretic hypotheses and their factual background"], "source": ["Z zool Syst Evolut-forsch"], "year": ["1990"], "volume": ["28"], "fpage": ["241"], "lpage": ["260"]}, {"surname": ["Wilson"], "given-names": ["EO"], "source": ["The Insect Societies"], "year": ["1971"], "publisher-name": ["Cambridge, MA: Harvard University Press"]}, {"surname": ["Emery"], "given-names": ["C"], "article-title": ["\u00dcber den Ursprung der dulotischen, parasitischen und myrmekophilen Ameisen"], "source": ["Biol Centralbl"], "year": ["1909"], "volume": ["29"], "fpage": ["352"], "lpage": ["362"]}, {"surname": ["Anderson", "Linksvayer", "Smith"], "given-names": ["KE", "TA", "CR"], "article-title": ["The causes and consequences of genetic caste determination in ants (Hymenoptera: Formicidae)"], "source": ["Myrmecol News"], "year": ["2008"], "volume": ["11"], "fpage": ["119"], "lpage": ["132"]}, {"surname": ["Agosti"], "given-names": ["D"], "article-title": ["A new inquiline ant (Hymenoptera: Formicidae) in "], "italic": ["Cataglyphis "], "source": ["J Nat Hist"], "year": ["1994"], "volume": ["28"], "fpage": ["913"], "lpage": ["919"], "pub-id": ["10.1080/00222939400770481"]}, {"surname": ["Beibl", "Buschinger", "Foitzik", "Heinze"], "given-names": ["J", "A", "S", "J"], "article-title": ["Phylogeny and phylogeography of the Mediterranean species of the parasitic ant genus "], "italic": ["Chalepoxenus ", "Temnothorax "], "source": ["Insect Soc"], "year": ["2007"], "volume": ["54"], "fpage": ["189"], "lpage": ["199"], "pub-id": ["10.1007/s00040-007-0932-y"]}, {"surname": ["Beibl", "Stuart", "Heinze", "Foitzik"], "given-names": ["J", "RJ", "J", "S"], "article-title": ["Six origins of slavery in formicoxenine ants"], "source": ["Insect Soc"], "year": ["2005"], "volume": ["52"], "fpage": ["291"], "lpage": ["297"], "pub-id": ["10.1007/s00040-005-0808-y"]}, {"surname": ["Sumner", "Aanen", "Delabie", "Boomsma"], "given-names": ["S", "DK", "J", "JJ"], "article-title": ["The evolution of social parasitism in "], "italic": ["Acromyrmex "], "source": ["Insect Soc"], "year": ["2004"], "volume": ["51"], "fpage": ["37"], "lpage": ["42"], "pub-id": ["10.1007/s00040-003-0723-z"]}, {"surname": ["Sanetra", "Buschinger"], "given-names": ["M", "A"], "article-title": ["Phylogenetic relationships among social parasites and their hosts in the ant tribe Tetramoriini (Hymenoptera: Formicidae)"], "source": ["Eur J Entomol"], "year": ["2000"], "volume": ["97"], "fpage": ["95"], "lpage": ["117"]}, {"surname": ["Baur", "Sanetra", "Chalwatzis", "Buschinger", "Zimmermann"], "given-names": ["A", "M", "N", "A", "FK"], "article-title": ["Sequence comparisons of the internal transcribed spacer region of ribosomal genes support close relationships between parasitic ants and their respective host species (Hymenoptera: Formicidae)"], "source": ["Insect Soc"], "year": ["1996"], "volume": ["43"], "fpage": ["53"], "lpage": ["67"], "pub-id": ["10.1007/BF01253956"]}, {"surname": ["Mueller", "Gerardo", "Aanen", "Six", "Schultz"], "given-names": ["UG", "NM", "DK", "DL", "TR"], "article-title": ["The evolution of agriculture in insects"], "source": ["Annu Rev Ecol Evol Syst"], "year": ["2005"], "volume": ["36"], "fpage": ["563"], "lpage": ["595"], "pub-id": ["10.1146/annurev.ecolsys.36.102003.152626"]}, {"surname": ["Maschwitz", "H\u00f6lldobler"], "given-names": ["U", "B"], "article-title": ["Der Kartonnestbau bei "], "italic": ["Lasius fuliginosus "], "source": ["Z vergl Physiol"], "year": ["1970"], "volume": ["66"], "fpage": ["176"], "lpage": ["189"], "pub-id": ["10.1007/BF00297777"]}, {"surname": ["Elliott"], "given-names": ["JSB"], "article-title": ["Fungi in the nests of ants"], "source": ["Trans Br Mycol Soc"], "year": ["1915"], "volume": ["5"], "fpage": ["138"], "lpage": ["142"]}, {"surname": ["Fresenius"], "given-names": ["JBGW"], "source": ["Beitr\u00e4ge zur Mykologie 2"], "year": ["1852"], "publisher-name": ["Frankfurt: HL Br\u00f6nner Verlag"]}, {"surname": ["Lagerheim"], "given-names": ["G"], "article-title": ["Ueber "], "italic": ["Lasius fuliginosus "], "source": ["Entomol Tidskr Arg"], "year": ["1900"], "volume": ["21"], "fpage": ["2"], "lpage": ["29"]}, {"surname": ["Schlick-Steiner", "Steiner", "Konrad", "Seifert", "Christian", "Moder", "Stauffer", "Crozier"], "given-names": ["BC", "FM", "H", "B", "E", "K", "C", "RH"], "article-title": ["Specificity and transmission mosaic of ant nest wall fungi"], "source": ["Proc Natl Acad Sci USA"], "year": ["2008"], "volume": ["105"], "fpage": ["941"], "lpage": ["944"], "pub-id": ["10.1073/pnas.0708320105"]}, {"surname": ["Wilson"], "given-names": ["EO"], "article-title": ["A monographic revision of the ant genus "], "italic": ["Lasius"], "source": ["Bull Mus Comp Zool"], "year": ["1955"], "volume": ["113"], "fpage": ["3"], "lpage": ["205"]}, {"surname": ["Hasegawa"], "given-names": ["E"], "article-title": ["Phylogeny and host-parasite relationships in social parasitism in "], "italic": ["Lasius "], "source": ["Entomol Sci"], "year": ["1998"], "volume": ["1"], "fpage": ["133"], "lpage": ["135"]}, {"surname": ["Ward", "Downie"], "given-names": ["PS", "DA"], "article-title": ["The ant subfamily Pseudomyrmecinae (Hymenoptera: Formicidae): phylogeny and evolution of big-eyed arboreal ants"], "source": ["Syst Entomol"], "year": ["2005"], "volume": ["30"], "fpage": ["310"], "lpage": ["335"], "pub-id": ["10.1111/j.1365-3113.2004.00281.x"]}, {"surname": ["Lucas", "Fresneau", "Kolmer", "Heinze", "Delabie", "Pho"], "given-names": ["C", "D", "K", "J", "JHC", "DB"], "article-title": ["A multidisciplinary approach to discriminating different taxa in the species complex "], "italic": ["Pachycondyla villosa "], "source": ["Biol J Linn Soc"], "year": ["2002"], "volume": ["75"], "fpage": ["249"], "lpage": ["259"], "pub-id": ["10.1111/j.1095-8312.2002.tb01425.x"]}, {"surname": ["Ward", "Brady"], "given-names": ["PS", "SG"], "article-title": ["Phylogeny and biogeography of the ant subfamily Myrmeciinae (Hymenoptera: Formicidae)"], "source": ["Invertebr Syst"], "year": ["2003"], "volume": ["17"], "fpage": ["361"], "lpage": ["386"], "pub-id": ["10.1071/IS02046"]}, {"surname": ["Brady"], "given-names": ["SG"], "article-title": ["Evolution of the army ant syndrome: The origin and long-term evolutionary stasis of a complex of behavioral and reproductive adaptations"], "source": ["Proc Natl Acad Sci USA"], "year": ["2003"], "volume": ["100"], "fpage": ["6576"], "lpage": ["6579"], "pub-id": ["10.1073/pnas.1137809100"]}, {"surname": ["Seifert", "Buschinger"], "given-names": ["B", "A"], "article-title": ["Pleometrotische Koloniegr\u00fcndung von "], "italic": ["Lasius meridionalis ", "Lasius paralienus "], "source": ["Myrmecologische Nachrichten"], "year": ["2001"], "volume": ["4"], "fpage": ["11"], "lpage": ["25"]}, {"surname": ["Wilson"], "given-names": ["EO"], "italic": ["Leptothorax duloticus "], "source": ["Evolution"], "year": ["1975"], "volume": ["29"], "fpage": ["108"], "lpage": ["119"], "pub-id": ["10.2307/2407145"]}, {"surname": ["Jenner"], "given-names": ["R"], "article-title": ["Presumed primitive until proving derived: appreciating evolutionary loss"], "source": ["The Systematist"], "year": ["2004"], "volume": ["23"], "fpage": ["12"], "lpage": ["16"]}, {"surname": ["Grandcolas", "Deleporte", "Desutter-Grandcolas", "Daugeron"], "given-names": ["P", "P", "L", "C"], "article-title": ["As many characters as possible should be included in the cladistic analysis"], "source": ["Cladistics"], "year": ["2001"], "volume": ["17"], "fpage": ["104"], "lpage": ["110"]}, {"surname": ["Schultz", "Churchill"], "given-names": ["TR", "GA"], "article-title": ["The role of subjectivity in reconstructing ancestral character states: A Bayesian approach to unknown rates, states, and transformation asymmetries"], "source": ["Syst Biol"], "year": ["1999"], "volume": ["48"], "fpage": ["651"], "lpage": ["664"], "pub-id": ["10.1080/106351599260229"]}, {"surname": ["Stille"], "given-names": ["M"], "article-title": ["Queen/worker thorax volume ratios and nest-founding strategies in ants"], "source": ["Oecologia"], "year": ["1996"], "volume": ["105"], "fpage": ["87"], "lpage": ["93"], "pub-id": ["10.1007/BF00328795"]}, {"surname": ["Keller", "Passera"], "given-names": ["L", "L"], "article-title": ["Size and fat content of gynes in relation to the mode of colony founding in ants (Hymenoptera; Formicidae)"], "source": ["Oecologia"], "year": ["1989"], "volume": ["80"], "fpage": ["236"], "lpage": ["240"]}, {"surname": ["Hurst", "Jiggins"], "given-names": ["GDD", "FM"], "article-title": ["Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts"], "source": ["Proc R Soc B"], "year": ["2005"], "volume": ["272"], "fpage": ["1525"], "lpage": ["1534"], "pub-id": ["10.1098/rspb.2005.3056"]}, {"surname": ["Ward"], "given-names": ["PS"], "article-title": ["A synoptic review of the ants of California (Hymenoptera: Formicidae)"], "source": ["Zootaxa"], "year": ["2005"], "volume": ["936"], "fpage": ["1"], "lpage": ["68"]}, {"surname": ["Hillis", "Bull"], "given-names": ["DM", "JJ"], "article-title": ["An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis"], "source": ["Syst Biol"], "year": ["1993"], "volume": ["42"], "fpage": ["182"], "lpage": ["192"], "pub-id": ["10.2307/2992540"]}, {"surname": ["Bolton", "Alpert", "Ward", "Naskrecki"], "given-names": ["B", "G", "PS", "P"], "source": ["Bolton's Catalogue of Ants of the World: 1758\u20132005"], "year": ["2007"], "publisher-name": ["Cambridge, MA: Harvard University Press"]}, {"surname": ["Wilson"], "given-names": ["EO"], "article-title": ["Tropical social parasites in the ant genus "], "italic": ["Pheidole"], "source": ["Insect Soc"], "year": ["1984"], "volume": ["31"], "fpage": ["316"], "lpage": ["334"], "pub-id": ["10.1007/BF02223615"]}, {"surname": ["Steiner", "Schlick-Steiner", "Sch\u00f6dl", "Espadaler", "Seifert", "Christian", "Stauffer"], "given-names": ["FM", "BC", "S", "X", "B", "E", "C"], "article-title": ["Phylogeny and bionomics of "], "italic": ["Lasius austriacus "], "source": ["Insect Soc"], "year": ["2004"], "volume": ["51"], "fpage": ["24"], "lpage": ["29"], "pub-id": ["10.1007/s00040-003-0699-8"]}, {"surname": ["Pontin"], "given-names": ["AJ"], "article-title": ["Population stabilization and competition between the ants "], "italic": ["Lasius flavus ", "L. niger "], "source": ["J Anim Ecol"], "year": ["1961"], "volume": ["30"], "fpage": ["47"], "lpage": ["54"], "pub-id": ["10.2307/2112"]}, {"surname": ["Savolainen"], "given-names": ["R"], "article-title": ["Phylogenetic relationships of "], "italic": ["Lasius"], "source": ["XIV International Congress of IUSSI. The Golden Jubilee Proceedings. 27 July \u2013 3 August 2002, Hokkaido University, Sapporo, Japan"], "year": ["2002"], "publisher-name": ["IUSSI, Sapporo"], "fpage": ["159"]}, {"surname": ["Seifert"], "given-names": ["B"], "article-title": ["A taxonomic revision of the Palaearctic members of the ant subgenus "], "italic": ["Lasius "], "source": ["Abh Ber Naturkundemus G\u00f6rlitz"], "year": ["1992"], "volume": ["66"], "fpage": ["1"], "lpage": ["67"]}, {"surname": ["Seifert"], "given-names": ["B"], "article-title": ["Social cleptogamy in the ant subgenus "], "italic": ["Chthonolasius "], "source": ["Abh Ber Naturkundemus G\u00f6rlitz"], "year": ["2006"], "volume": ["77"], "fpage": ["251"], "lpage": ["276"]}, {"surname": ["Wing"], "given-names": ["MW"], "article-title": ["Taxonomic revision of the Nearctic ant genus "], "italic": ["Acanthomyops "], "source": ["Mem Cornell Univ Agr Exp Stn"], "year": ["1968"], "volume": ["405"], "fpage": ["1"], "lpage": ["173"]}, {"surname": ["Yamauchi", "Hayashida"], "given-names": ["K", "K"], "article-title": ["Taxonomic studies on the genus "], "italic": ["Lasius ", "Lasius"], "source": ["J Fac Sci Hokkaido Univ Ser 6, Zool"], "year": ["1970"], "volume": ["17"], "fpage": ["501"], "lpage": ["519"]}, {"surname": ["Seifert"], "given-names": ["B"], "article-title": ["A revision of the European species of the ant subgenus "], "italic": ["Chthonolasius "], "source": ["Entomol Abh Staatl Mus Tierkd Dres"], "year": ["1988"], "volume": ["51"], "fpage": ["143"], "lpage": ["180"]}, {"surname": ["Johnson", "Yoshizawa", "Smith"], "given-names": ["KP", "K", "VS"], "article-title": ["Multiple origins of parasitism in lice"], "source": ["Proc R Soc London (B)"], "year": ["2004"], "volume": ["271"], "fpage": ["1771"], "lpage": ["1776"], "pub-id": ["10.1098/rspb.2004.2798"]}, {"surname": ["Palumbi", "Martin", "Romano", "McMillan", "Stice", "Grabowski"], "given-names": ["SR", "A", "S", "WO", "L", "G"], "source": ["The Simple Fool's Guide to PCR, Version 2.0. Privately published"], "year": ["1991"], "publisher-name": ["University of Hawaii"]}, {"surname": ["Xia", "Xie"], "given-names": ["X", "Z"], "article-title": ["DAMBE: Data analysis in molecular biology and evolution"], "source": ["J Heredity"], "year": ["2001"], "volume": ["92"], "fpage": ["371"], "lpage": ["373"], "pub-id": ["10.1093/jhered/92.4.371"]}, {"surname": ["Lemmon", "Moriarty"], "given-names": ["AR", "EC"], "article-title": ["The importance of proper model assumption in Bayesian phylogenetics"], "source": ["Syst Biol"], "year": ["2004"], "volume": ["53"], "fpage": ["267"], "lpage": ["277"], "pub-id": ["10.1080/10635150490423520"]}, {"surname": ["Nylander"], "given-names": ["JAA"], "source": ["MrModeltest v2. Program distributed by the author"], "year": ["2004"], "publisher-name": ["Evolutionary Biology Centre, Uppsala University"]}, {"surname": ["Kass", "Raftery"], "given-names": ["RE", "AE"], "article-title": ["Bayes factors"], "source": ["J Am Stat Assoc"], "year": ["1995"], "volume": ["90"], "fpage": ["773"], "lpage": ["795"], "pub-id": ["10.2307/2291091"]}, {"surname": ["Hillis", "Bull"], "given-names": ["DM", "JJ"], "article-title": ["An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis"], "source": ["Syst Biol"], "year": ["1993"], "volume": ["42"], "fpage": ["182"], "lpage": ["192"], "pub-id": ["10.2307/2992540"]}, {"surname": ["Swofford"], "given-names": ["DL"], "source": ["PAUP*: Phylogenetic analysis using parsimony (*and other methods), ver. 4.0b4a"], "year": ["2003"], "publisher-name": ["Sunderland, MA: Sinauer"]}]
BMC Evol Biol. 2008 Aug 19; 8:237