A molecular framework for the phylogeny of the ant subfamily dolichoderinae

Complete Nucleotide Sequence of the Mitogenome of Tapinoma ibericum (Hymenoptera: Formicidae: Dolichoderinae), Gene Organization and Phylogenetics Implications for the Dolichoderinae Subfamily

The ant Tapinoma ibericum Santschi, 1925 is native to the Iberian Peninsula. This species, as well as other species from the Tapinoma nigerrimum complex, could form supercolonies that make these species potentially invasive and could give rise to pests. Recently a mature colony from this species has been found in the Isle of Wight (United Kingdom). Mitogenomes have been used to study the taxonomy, biogeography and genetics of species, improving the development of strategies against pest invasion. However, the number of available mitogenomes from the subfamily Dolichoderinae is still scarce and only two of these mitogenomes belong to Tapinoma species. Herein, the complete mitogenome of T. ibericum is presented in order to increase the molecular information of the genus. The T. ibericum mitogenome, retrieved by Next-Generation Sequencing data, is 15,715 bp in length. It contains the typical set of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNAs and the A + T-rich control region. Comparisons of the T. ibericum mitogenome with other dolichoderine mitogenomes revealed the existence of four gene rearrangements in relation with the ancestral insect mitogenome. One of these rearrangements, involving the tRNA-Ile, tRNA-Gln and tRNA-Met genes, was found in most of the analyzed ant mitogenomes. Probably this rearrangement was an ancestral or plesiomorphic character in Formicidae. Interestingly, another rearrangement that affects to tRNA-Trp, tRNA-Cys and tRNA-Tyr genes was found only in Tapinoma species. This change could be a synapomorphic character for the genus Tapinoma, and could be used as a phylogenetic marker. Additionally, a phylogenetic analysis was performed using the protein-coding gene sequences from available Dolichoderinae mitogenomes, as well as mitogenomes from representative species from other Formicidae subfamilies. Results support the monophyletic nature of the genus Tapinoma placing it within the same clade as the rest of Dolichoderinae species.

Social structure of perennial Vespula squamosa wasp colonies

Many social species show variation in their social structure in response to different environmental conditions. For example, colonies of the yellowjacket wasp Vespula squamosa are typically headed by a single reproductive queen and survive for only a single season. However, in warmer climates, V. squamosa colonies sometimes persist for multiple years and can grow to extremely large size. We used genetic markers to understand patterns of reproduction and recruitment within these perennial colonies. We genotyped V. squamosa workers, pre-reproductive queens, and males from perennial colonies in the southeastern United States at 10 polymorphic microsatellite loci and one mitochondrial DNA locus. We found that V. squamosa from perennial nests were produced by multiple reproductives, in contrast to typical annual colonies. Relatedness of nestmates from perennial colonies was significantly lower than relatedness of nestmates from annual colonies. Our analyses of mitochondrial DNA indicated that most V. squamosa perennial colonies represented semiclosed systems whereby all individuals belonged to a single matriline despite the presence of multiple reproductive females. However, new queens recruited into perennial colonies apparently mated with non-nestmate males. Notably, perennial and annual colonies did not show significant genetic differences, supporting the hypothesis that perennial colony formation represents an instance of social plasticity. Overall, our results indicate that perennial V. squamosa colonies show substantial changes to their social biology compared to typical annual colonies and demonstrate variation in social behaviors in highly social species.

A new phylogenetic protocol: dealing with model misspecification and confirmation bias in molecular phylogenetics

Molecular phylogenetics plays a key role in comparative genomics and has increasingly significant impacts on science, industry, government, public health and society. In this paper, we posit that the current phylogenetic protocol is missing two critical steps, and that their absence allows model misspecification and confirmation bias to unduly influence phylogenetic estimates. Based on the potential offered by well-established but under-used procedures, such as assessment of phylogenetic assumptions and tests of goodness of fit, we introduce a new phylogenetic protocol that will reduce confirmation bias and increase the accuracy of phylogenetic estimates.

Evaluating the tradeoffs of a generalist parasitoid fungus, Ophiocordyceps unilateralis, on different sympatric ant hosts

It is essential for the survival and reproduction of parasitoids to adapt to the fluctuating host resources. Phenotypic plasticity may enable a parasitoid species to successfully achieve its control over a range of host species to maximize fitness in different hosts that may each require dissimilar, possibly conflicting, specific adaptations. However, there is limited information on how the fitness effects of host switching partition into costs due to the novelty of host species, where unfamiliarity with host physiological and morphological changes and its anti-parasite defenses reduces parasitoid growth, survivorship and/or reproductive success. In this study, the parasitoid fungus Ophiocordyceps unilateralis sensu lato was found to sympatrically infect a principal host ant species and other alternative sympatric hosts in the forest of central Taiwan. We herein report that the occurrence of ant infections by O. unilateralis s.l. shows spatial and temporal variation patterns on different host species. Results showed that the height from the ground to the leaf where the infected ants grip on, perithecia-forming ability, and growth rate of the stroma of the parasitoid fungus were dissimilar on different host species. These host range expansions not only related the fitness of O. unilateralis s.l. but also influenced the expression of extended phenotypic traits. Our findings revealed that a generalist parasitoid fungus suffered an evolutionary tradeoff between host breadth expansion and host-use efficiency.

DNA barcoding of aphid-associated ants (Hymenoptera, Formicidae) in a subtropical area of southern China

As one of the most abundant and complex groups of terrestrial insects, ants have associations with many other organismal groups, such as hemipteran insects producing honeydew. With the aim of expanding the knowledge base of ant species associated with aphids, this study analyzed mitochondrial COI barcodes of 301 ant samples for 37 aphid-associated ant species in a subtropical area of southern China. Sequence analyses revealed that the intraspecific and interspecific distances ranged from zero to 7.7%% and 0.2 to 31.7%, respectively. Three barcoding approaches - Automatic Barcode Gap Discovery, Bayesian Poisson Tree Processes and Generalized Mixed Yule-coalescent - were used to help delimit ant species based on COI sequences, and their results corresponded well with most of the morphospecies. All three approaches indicate cryptic diversity may exist within Tetramorium bicarinatum and Technomyrmex albipes, with intraspecific genetic distances of 7.7% and 6.24%, respectively. Our analyses also reported five species for the first time from Fujian Province of China, and the COI sequences of nine species are newly added into the GenBank. This study provides information about species diversity of aphid-associated ants in subtropical China and compiles a DNA barcode reference library for future ant barcoding work.

Genetic diversity and colony structure of Tapinoma melanocephalum on the islands and mainland of South China

AIM

Tapinoma melanocephalum is listed as one of the most important invasive pest species in China. Information regarding the patterns of invasion and effects of geographic isolation on the population genetics of this species is largely lacking.

LOCATION

South China.

METHODS

To address this problem, we genotyped 39 colonies (two colonies were collapsed due to genetic similarity) using microsatellite markers and mitochondrial DNA sequencing to compare colony genetic structure of T. melanocephalum on the mainland and islands of South China.

RESULTS

An analysis of the colony genotypes showed that the genetic diversity of the mainland population was slightly higher than that of the island populations but not significantly so. However, the observed heterozygosity on Shangchuan Island (SCD) was significantly lower than that of the other colonies. We also found six haplotypes in 111 mitochondrial DNA COI sequences. The relatedness (r) value between colonies of SCD was 0.410, higher than that of the other populations. The genetic clusters among colonies were not related to geographic locations and exhibited admixture likely due to frequent human-mediated dispersal associated with trade between the mainland population and the islands. Pairwise FSTs between populations showed differentiation among mainland populations, while SCD displayed high levels of divergence (FST > 0.15) from most mainland populations. There was no significant isolation by distance among colonies. Most populations showed signs of a bottleneck effect.

MAIN CONCLUSIONS

Our study suggests that there was no significant difference in the genetic diversity among the islands and the mainland; however, the lower genetic diversity, the higher degree of genetic divergence from other colonies, and the higher relatedness among nestmates made the SCD population stand out from all the others.

ModelFinder: fast model selection for accurate phylogenetic estimates

Model-based molecular phylogenetics plays an important role in comparisons of genomic data, and model selection is a key step in all such analyses. We present ModelFinder, a fast model-selection method that greatly improves the accuracy of phylogenetic estimates by incorporating a model of rate heterogeneity across sites not previously considered in this context and by allowing concurrent searches of model space and tree space.

Molecular Profile of the Brazilian Weaver Ant Camponotus textor Forel (Hymenoptera, Formicidae)

Camponotus textor Forel is, to date, the only weaver ant recorded from Brazil, and all existing studies on the species are restricted to describing its weaving and nesting behaviors. The aim of this work is to establish the molecular profile of the species. We sampled eight different colonies by sequencing mitochondrial genes (COI, transfer DNA (tRNA), and an intergenic spacer) and the nuclear gene 28S ribosomal DNA (rDNA). We then assessed haplotype diversity and also analyzed distribution patterns of this species based on the correlation between genetic and geographic distances. Our results provide an additional tool for species identification by identifying new regions that can be used as molecular markers for barcoding (such as the intergenic spacer (IGS) and tRNA-Leu). In addition, the phylogenetic analysis revealed that C. textor has features that could be associated with deep population divergences. We identified a wide range of mitotypes and three distinct groups, suggesting a possible reduction of gene flow between colonies.

The mitochondrial genome of the leaf-cutter ant Atta laevigata: a mitogenome with a large number of intergenic spacers

In this paper we describe the nearly complete mitochondrial genome of the leaf-cutter ant Atta laevigata, assembled using transcriptomic libraries from Sanger and Illumina next generation sequencing (NGS), and PCR products. This mitogenome was found to be very large (18,729 bp), given the presence of 30 non-coding intergenic spacers (IGS) spanning 3,808 bp. A portion of the putative control region remained unsequenced. The gene content and organization correspond to that inferred for the ancestral pancrustacea, except for two tRNA gene rearrangements that have been described previously in other ants. The IGS were highly variable in length and dispersed through the mitogenome. This pattern was also found for the other hymenopterans in particular for the monophyletic Apocrita. These spacers with unknown function may be valuable for characterizing genome evolution and distinguishing closely related species and individuals. NGS provided better coverage than Sanger sequencing, especially for tRNA and ribosomal subunit genes, thus facilitating efforts to fill in sequence gaps. The results obtained showed that data from transcriptomic libraries contain valuable information for assembling mitogenomes. The present data also provide a source of molecular markers that will be very important for improving our understanding of genomic evolutionary processes and phylogenetic relationships among hymenopterans.

Occurrence of three haplotypes of Linepithema micans (Forel) (Hymenoptera: Formicidae) in Southern Brazil

Linepithema micans (Forel) (Hymenoptera: Formicidae) is reported to occur from eastern Brazil to central Argentina in pasture or grassland, forest and second growth riparian forest, nesting under stones, rotting wood, and sandy soil. However, information on this species is poor and its ecological interactions and role as pests are unknown. Linepithema humile (Mayr), a closely related species to L. micans, known as the Argentine ant, is native to South America, and was accidentally introduced to several regions of the world. Recent studies have shown that other related species, such as L. micans, could become as pestiferous as L. humile because of its phylogenetic proximity. Samples of L. micans from different habitats in Southern Brazil were characterized by sequencing of the mitochondrial DNA. Sequences were compared to previously obtained sequences from samples of L. humile and the genetic distance and differences in the tRNALeu structure were investigated. Our data identified three haplotypes of L. micans, two of which were observed in ant populations closely associated with the Brazilian ground pearl Eurhizococcus brasiliensis (Hempel) (Hemiptera: Margarodidae), a soil scale that is a serious pest of vineyards. The third haplotype was identified in ants from populations invading residences in urban habitats.

Molecular phylogeny and biogeography of the spider ants, genus Leptomyrmex Mayr (Hymenoptera: Formicidae)

This study provides the first phylogenetic reconstruction of the ant genus Leptomyrmex Mayr, a prominent endemic component of rain forest and wet sclerophyll forest in Australia, New Guinea and New Caledonia. Five genes are used to reconstruct phylogeny and estimate of ages of diversification in order to test congruence of the history of nuclear and mitochondrial genes: three protein-coding nuclear genes: arginine kinase (argK, 897 bp), long wavelength rhodopsin (LW Rh, 546 bp) and wingless (Wg, 409 bp), as well as the large subunit ribosomal gene 28S (482 bp) and the mitochondrial gene cytochrome oxidase I (COI, 658 bp). Four different partitioning schemes were tested for optimal resolving power; results show that partitioning by gene, translational pattern and codon position were uniformly favoured over less complex partitions. Nuclear markers showed relatively minor sequence divergence and provided strongly supported topology; phylogeny based solely on mtDNA produced somewhat conflicting topology but offered little power to resolve species complexes. Monophyly of the genus Leptomyrmex was recovered, as was the sister-group relationship of 'micro-' and 'macro-'Leptomyrmex species. Divergence dating analyses estimate that Leptomyrmex arose in the Eocene (stem age ∼ 44 million years ago (ma)), and that the 'macro-' species diverged from the 'micro-' species in the early Oligocene (∼ 31 ma). Diversification of the crown group 'macro-' and 'micro-'Leptomyrmex occurred in the Miocene (∼ 15 ma and 7.9 ma, respectively). New Guinean and New Caledonian lineages appear to have diverged from Australian lineages only recently (∼ 4.7 ma and 10.3 ma, respectively), and the latter clade is inferred to have reached New Caledonia from Australia via long distance dispersal. These results challenge previous hypotheses of Leptomyrmex classification and assumptions about their historical dispersal, but are in agreement with the current knowledge of the geological history of Melanesia.

Phylogeny and biogeography of dolichoderine ants: effects of data partitioning and relict taxa on historical inference

Ants (Hymenoptera: Formicidae) are conspicuous organisms in most terrestrial ecosystems, often attaining high levels of abundance and diversity. In this study, we investigate the evolutionary history of a major clade of ants, the subfamily Dolichoderinae, whose species frequently achieve ecological dominance in ant communities. This group has also produced some of the world's most successful invasive ants. We use an extensive molecular data set ( approximately 9 kb of sequence data from 10 nuclear genes, covering 48 dolichoderine species and 6 outgroup taxa) to infer the phylogenetic relationships, divergence dates, and biogeographic history of these ants. We evaluate the effects of data partitioning and outgroup composition on phylogenetic inference by estimating relationships under a series of increasingly partitioned data sets and by running analyses both with and without Aneuretus simoni, a rare and localized species that is the nearest living relative of Dolichoderinae. We also examine the effects of excluding 2 data partitions with significant base composition heterogeneity. Our results reveal 4 well-supported and mutually exclusive clades of dolichoderines, corresponding to 4 newly defined tribes: Bothriomyrmecini (B), Dolichoderini (D), Leptomyrmecini (L), and Tapinomini (T). All Bayesian and likelihood analyses yield the same unrooted (ingroup-only) topology, ((D,L),(B,T)), with the outgroups attaching either on the Dolichoderini branch or on the Tapinomini branch. Placement of the root is highly sensitive to choice of model partition and to inclusion/exclusion of Aneuretus. Bayes' factors strongly favor the more partitioned models, and in these Tapinomini is recovered as sister to the remaining dolichoderines, but only if Aneuretus is included. Exclusion of Aneuretus precludes recovery of this topology in all but the most highly partitioned Bayesian analyses and then only with nonsignificant support, underscoring the importance of relict, taxonomically isolated taxa for phylogenetic inference. Removal of 2 partitions with heterogeneous base composition also markedly increases support for placement of the root on the Tapinomini branch. Our divergence date estimates and biogeographic analyses indicate that crown-group dolichoderines arose about 65 million years ago (Ma), although this was preceded by a substantial period (30 million years) of stem group evolution. The 4 extant tribes are estimated to have crown-group origins in the late Paleocene or Eocene (40-60 Ma). Tapinomini and Bothriomyrmecini originated in the Paleotropics and subsequently dispersed to other biogeographic regions. Crown-group Leptomyrmecini arose and diversified in the Neotropics, but they also gave rise to one clade that colonized Australia about 30 Ma and subsequently experienced a massive radiation on that continent. This event occurred later than the diversification of dolichoderines in the northern hemisphere, so that by the time dolichoderines came to dominate the Australian fauna they had already declined in abundance in the Holarctic region.

Molecular phylogeny of the benthic shallow-water octopuses (Cephalopoda: Octopodinae)

Octopus has been regarded as a "catch all" genus, yet its monophyly is questionable and has been untested. We inferred a broad-scale phylogeny of the benthic shallow-water octopuses (subfamily Octopodinae) using amino acid sequences of two mitochondrial DNA genes: Cytochrome oxidase subunit III and Cytochrome b apoenzyme, and the nuclear DNA gene Elongation Factor-1alpha. Sequence data were obtained from 26 Octopus species and from four related genera. Maximum likelihood and Bayesian approaches were implemented to estimate the phylogeny, and non-parametric bootstrapping was used to verify confidence for Bayesian topologies. Phylogenetic relationships between closely related species were generally well resolved, and groups delineated, but the genes did not resolve deep divergences well. The phylogenies indicated strongly that Octopus is not monophyletic, but several monophyletic groups were identified within the genus. It is therefore clear that octopodid systematics requires major revision.

The evolution of worker caste diversity in social insects

Morphological diversification of workers is predicted to improve the division of labor within social insect colonies, yet many species have monomorphic workers. Individual-level selection on the reproductive capacities of workers may counter colony-level selection for diversification, and life-history differences between species (timing of caste determination, colony size, genetic variation available) may mediate the strength of this selection. We tested this through phylogenetically independent contrast analyses on a new data set for 35 ant species. Evidence was found that early divergence of queen-worker developmental pathways may facilitate the evolution of worker diversity because queen-worker dimorphism was strongly positively associated with diversity. By contrast, risks for colonies that invest in specialized workers and colony size effects on costs of worker reproduction seem unlikely to strongly affect the evolution of worker diversity because there was no significant association between colony size and diversity when controlling statistically for queen-worker dimorphism. Finally, worker diversity was greater in species with multiple lineages per colony, and it was negatively associated with relatedness in monogynous species. This could be due to high intracolonial genetic variance favoring the expression and evolution of great worker diversity or to diversity evolving more easily when there is selection for repression of worker reproduction (worker policing).

Ecological stoichiometry of ants in a New World rain forest

C:N stoichiometry was investigated in relation to diet (delta(15)N), N-deprivation, and worker body size for a diverse assemblage of tropical Amazonian ants. Relative nitrogen (N) deprivation was assayed for 54 species as an exchange ratio (ER), defined as SUCmin/AAmin, or the minimum sucrose concentration, divided by the minimum amino acid concentration, accepted as food by >/=50% of tested workers. On average, N-deprivation (ER) was almost fivefold greater for N-omnivorous and N-herbivorous (N-OH) taxa than for N-carnivores. In two-way ANOVAs at three taxonomic levels (species and species groups, genera, and tribes), higher ER was associated with small body size and (marginally) with less carnivorous diets. ERs did not differ systematically between trophobiont-tending and "leaf-foraging" functional groups, but specialized wound-feeders and coccid-tenders were prominent among high ER taxa. Paradoxically, some high ER taxa were among the most predatory members of their genera or subfamilies. Biomass % N was lower in N-OH taxa than in carnivores and varied inversely with N-deprivation (log ER) in the former taxa only. In an expanded data set, N-content increased allometrically in N-OHs, N-carnivores, and all ants combined, and with carnivory in large-bodied ants only. Exceptional taxa included small-bodied and predaceous Wasmannia, with high % N despite high ER, and Linepithema, with the lowest % N despite high delta(15)N. Patterns in C:N stoichiometry are explained largely at the genus level and above by elemental composition of alarm/defensive/offensive chemical weaponry and, perhaps in some cases, by reduced N investment in cuticle in taxa with high surface:volume ratios. Several consequences of C:N stoichiometry identify Azteca, and possibly Crematogaster, as taxa preadapted for their roles as prominent associates of myrmecophytes. C:N stoichiometry of ants should be incorporated into models of strategic colony design and examined in a phylogenetic context as opportunities permit.

Phylogeny of ants (Formicidae) based on morphology and DNA sequence data

In order to reconstruct ants' phylogeny, we analysed DNA sequences for two nuclear genes, abdominal-A and Ultrabithorax, from 49 species of ants and two outgroups. As these genes control the development of the first segments of the abdomen in insects, which are very variable in ants (petiole, postpetiole, and gaster constriction), we hypothesized that the morphological variations between the subfamilies may be correlated with mutations of some abd-A or Ubx regions. Contrarily to our hypothesis, these sequences are highly conserved. The differences observed concern mainly third codon positions and present some saturation. Phylogenetic reconstructions were carried out using the genetic raw sequence data and by combining them with a set of morphological data (Total Evidence). Relations among subfamilies of ants remains poorly resolved with molecular data only, but adding these data to morphological characters confirms and reinforce the topology of : a Poneroid complex [Ponerinae, Cerapachyinae, Leptanillinae and army ants], a Formicoid complex [Dolichoderinae, Formicinae] and a Myrmecoid complex [Myrmicinae, Myrmeciinae, Pseudomyrmecinae, Nothomyrmeciinae]. Our molecular results allow resolution near the branch tips and three subfamilies (Dolichoderinae, Formicinae and Pseudomyrmecinae) always appear as monophyletic. The Formicinae and the Dolichoderinae have close relationships. The Camponotini appear as a strong clade inside the Formicinae. The Ponerinae are separated in two parts: the Ectatommini and all other tribes. The Cerapachyinae, Dorylinae, and Ecitoninae belong to the same clade, the Cerapachyinae being confirmed in their subfamily status. The Myrmicinae appears to be very heterogeneous, with the Attini forming a very stable and well-separated group.

Liquid-feeding performances of ants (Formicidae): ecological and evolutionary implications

Disparities in liquid-feeding performances of major ant taxa have likely been important to resource partitioning among ants, to the nature and composition of ant partnerships with plants and sap-feeding trophobionts, and to ecological and evolutionary diversification of ant taxa. We measured performance volumetrically for individual workers of 77 ant species from lowland rain forests of Amazonia and Borneo and three key North American taxa. In trials with 9% sucrose solution, performances were strongly related to body size (and alitrunk length) and to proventricular structure at generic to subfamilial levels. Highly modified proventriculi were associated with disproportionately large load sizes in Formicinae and certain small-bodied Dolichoderinae. These same taxa also ingested liquids more rapidly during foraging than did similar-sized species with plesiomorphic proventriculi. Secondarily reduced foraging performances of several formicines likely reflect ecological or evolutionary trade-offs related to dietary specialization or anti-predator defenses. Across formicines and dolichoderines, performances differed by functional group. Relatively small loads and slow uptake characterized species tending trophobionts (mainly Hemiptera) day and night in large worker aggregations. Large loads and rapid uptake typified solitary, diurnal "leaf-foragers." Intermediate feeding performances characterized a few species that both tended trophobionts in small aggregates and frequently foraged alone.

A tree island approach to inferring phylogeny in the ant subfamily Formicinae, with especial reference to the evolution of weaving

The ant subfamily Formicinae is a large assemblage (2458 species (J. Nat. Hist. 29 (1995) 1037), including species that weave leaf nests together with larval silk and in which the metapleural gland-the ancestrally defining ant character-has been secondarily lost. We used sequences from two mitochondrial genes (cytochrome b and cytochrome oxidase 2) from 18 formicine and 4 outgroup taxa to derive a robust phylogeny, employing a search for tree islands using 10000 randomly constructed trees as starting points and deriving a maximum likelihood consensus tree from the ML tree and those not significantly different from it. Non-parametric bootstrapping showed that the ML consensus tree fit the data significantly better than three scenarios based on morphology, with that of Bolton (Identification Guide to the Ant Genera of the World, Harvard University Press, Cambridge, MA) being the best among these alternative trees. Trait mapping showed that weaving had arisen at least four times and possibly been lost once. A maximum likelihood analysis showed that loss of the metapleural gland is significantly associated with the weaver life-pattern. The graph of the frequencies with which trees were discovered versus their likelihood indicates that trees with high likelihoods have much larger basins of attraction than those with lower likelihoods. While this result indicates that single searches are more likely to find high- than low-likelihood tree islands, it also indicates that searching only for the single best tree may lose important information.

Relationships among native and introduced populations of the Argentine ant (Linepithema humile) and the source of introduced populations

The Argentine ant (Linepithema humile) is a damaging invasive species that has become established in many Mediterranean-type ecosystems worldwide. To identify likely sources of introduced populations we examined the relationships among native Linepithema populations from Argentina and Brazil and introduced populations of L. humile using mitochondrial cytochrome b sequence data and nuclear microsatellite allele frequencies. The mitochondrial phylogeny revealed that the populations in Brazil were only distantly related to both the introduced populations and the native populations in Argentina, and confirmed that populations in Brazil, previously identified as L. humile, are likely a different species. The microsatellite-based analysis provided resolution among native and introduced populations of L. humile that could not be resolved using the mitochondrial sequences. In the native range, colonies that were geographically close to one another tended to be genetically similar, whereas more distant colonies were genetically different. Most samples from the introduced range were genetically similar, although some exceptions were noted. Most introduced populations were similar to native populations from the southern Rio Parana and were particularly similar to a population from Rosario, Argentina. These findings implicate populations from the southern Rio Parana as the most likely source of introduced populations. Moreover, these data suggest that current efforts to identify natural enemies of the Argentine ant for biological control should focus on native populations in the southern Rio Parana watershed.