**Diversity and Distribution of the Dominant Ant Genus** *Anonychomyrma* **(Hymenoptera: Formicidae) in the Australian Wet Tropics**

#### **Lily Leahy 1,\*, Brett R. Sche**ff**ers 2, Stephen E. Williams <sup>1</sup> and Alan N. Andersen <sup>3</sup>**


Received: 18 November 2020; Accepted: 7 December 2020; Published: 14 December 2020

**Abstract:** *Anonychomyrma* is a dolichoderine ant genus of cool-temperate Gondwanan origin with a current distribution that extends from the north of southern Australia into the Australasian tropics. Despite its abundance and ecological dominance, little is known of its species diversity and distribution throughout its range. Here, we describe the diversity and distribution of *Anonychomyrma* in the Australian Wet Tropics bioregion, where only two of the many putative species are described. We hypothesise that the genus in tropical Australia retains a preference for cool wet rainforests reminiscent of the Gondwanan forests that once dominated Australia, but now only exist in upland habitats of the Wet Tropics. Our study was based on extensive recent surveys across five subregions and along elevation and vertical (arboreal) gradients. We integrated genetic (CO1) data with morphology to recognise 22 species among our samples, 20 of which appeared to be undescribed. As predicted, diversity and endemism were concentrated in uplands above 900 m a.s.l. Distribution modelling of the nine commonest species identified maximum temperature of the warmest month, rainfall seasonality, and rainfall of the wettest month as correlates of distributional patterns across subregions. Our study supported the notion that *Anonychomyrma* radiated from a southern temperate origin into the tropical zone, with a preference for areas of montane rainforest that were stably cool and wet over the late quaternary.

**Keywords:** ant diversity; biogeography; species delimitation; Dolichoderinae; species distribution models; climatic gradients; wet tropics; climate change

#### **1. Introduction**

The ubiquity, high species richness, and ecological dominance of ants make them model organisms for understanding spatial patterns of diversity and community assembly [1]. Ant distribution is tightly coupled to climate [2], so that an understanding of ant species distributions across climatic gradients provides clues as to how ectothermic species coped with climatic conditions in the past and provides important insights into how they are likely to respond to a future climate.

*Anonychomyrma* Donisthorpe 1947 is an ecologically dominant ant genus in relatively cool and mesic habitats of Australia and New Guinea, extending northwest into Southeast Asia and east to the Solomon Islands [3]. The genus contains 27 described species and 5 subspecies, with most of these from New Guinea [4,5]. The genus is a member of an Australian dolichoderine subgroup that is thought to be of Gondwanan origin and that most likely originated in Australia approximately 30 million years ago from an ancestor shared with southern South America through a land connection via Antarctica [6]. *Anonychomyrma* is particularly common in heathlands, woodlands, and open forests of temperate southern Australia ([7,8]; in both referenced articles, the genus is referred to as the *nitidiceps* group of *Iridomyrmex*) and in rainforest of subtropical and tropical mountains of Queensland and New Guinea [5,9,10]. Most species in drier habitats nest in the ground, whereas rainforest species appear to be primarily arboreal.

The Australian Wet Tropics (AWT) bioregion of northeastern Australia represents the largest remaining remnant of the Gondwanan mesothermal rainforests that once dominated Australia [11]. Two species of *Anonychomyrma* have been described from the region. One of these is *Anonychomyrma gilberti* (Forel) (Figure 1), a widely distributed species in the region where it is a dominant ant in the rainforest canopy [10,12]. The second species is *Anonychomyrma malandana* (Forel), which was described from Malanda on the Atherton Tablelands in 1915. Since its original description, we are aware of only one other publication documenting it, where it was commonly recorded on Mt. Windsor (Figure 2) and referred to as *Anonychomyrma* sp. D [9]. Despite just two described species, the AWT *Anonychomyrma* fauna is highly diverse, with [9] 15 species reported during ant surveys of just four of the thirteen upland subregions.

**Figure 1.** *Anonychomyrma gilberti*, showing head (**A**) and lateral (**B**) views.

**Figure 2.** Map of the Australian Wet Tropics (AWT) showing survey site locations for each subregion (naming of subregions follows [13]). Rainforest distribution is shown in green. BMC is the Black Mountain Corridor and represents a well-known biogeographic barrier for numerous taxa within the AWT.

Both contemporary and historical factors drive the distribution and diversity of species and thereby, the geographic structure of biological communities [14,15]. The biogeography of the AWT is strongly driven by historic processes that occurred during climate fluctuations of the late quaternary (last ~18,000 years) when rainforests contracted and expanded several times [16–19]. As rainforests contracted to wetter uplands, many lowland species became locally extinct, leading to a substantial loss of species richness in the lowlands [20]. As a result, many vertebrate and invertebrate species that are endemic to the AWT are restricted to elevations above 300 m a.s.l. that have had historically stable climates and consistently supported rainforests throughout this period of high climatic fluctuations [21]. In addition, there is an older biogeographic barrier along the Black Mountain Corridor (BMC), just south of the Carbine uplands (Figure 2), that has possibly existed since the Pliocene more than 2.5 Ma

ago [22] and did not support rainforest until ~8 ka ago [11,23]. There is evidence in multiple taxa of a phylogeographic break across the Black Mountain Corridor that structures phylogenetic diversity at interspecific as well as intraspecific levels [19,24,25], such that many species have either a northern or southern restricted distribution [24].

Biogeographic patterns in the region are also strongly influenced by contemporary climatic conditions at both regional and local scales. At a regional scale, seasonality of rainfall and temperature vary markedly among subregions [26]. Rainfall seasonality is negatively correlated with the abundance of birds in the AWT, with dry season severity creating bottlenecks in critical food resources such as insects, fruits, and nectar [26]. Similarly, there are marked seasonal fluctuations in the abundance of homopteran insects in the AWT [27], which, through their production of honeydew, provide a particularly important food resource for ants, including *A. gilberti* [28]. As such, honeydew-reliant ants, such as species of *Anonychomyrma*, might be expected to show strong biogeographic patterns in response to rainfall seasonality.

At a finer geographic scale, temperature and rainfall both vary markedly with elevation. Mean temperature in the AWT declines by approximately 1 ◦C per 200 m increase in elevation [29]. Elevations above 600–800 m a.s.l. receive year-round moisture via the orographic cloud layer that creates a misty and cool environment [30]. The orographic cloud layer has been linked to a significant change in ant species composition [9]. Finally, at a local scale, there are strong microclimatic gradients across vertical space from the ground to the canopy, a general feature of rainforest environments globally [31]. A recent study of ant distributions (including nine species of *Anonychomyrma*) in the AWT found a strong positive correlation between vertical niche breadth, and therefore exposure to microclimatic variation, and elevation range size [32], indicating that vertical microhabitat associations of ants influence climatic niches and distributional patterns along elevation gradients.

Here, we used extensive collections from five mountain ranges to document the AWT *Anonychomyrma* fauna and to investigate patterns of species diversity and distributions across latitude, elevation, and vertical habitat space. Like many other dolichoderine genera [3], *Anonychomyrma* is morphologically conservative and therefore taxonomically challenging, so we used an integrated taxonomic approach that included CO1 gene barcoding [33] to inform species boundaries. We identified the key climatic drivers of species distributions and employed species distribution modelling based on climatic niches to compare actual with potential distributions and thus identified factors that may have shaped contemporary subregional patterns of distribution. Specifically, we address the following questions:


We predicted that given the Gondwanan origin of *Anonychomyrma* and its prevalence in cool temperate Australia, most species will have montane distributions in the AWT, with high levels of diversity and endemism in the cool upland regions that have retained stable rainforest vegetation and climate conditions more reminiscent of the Gondwanan rainforest that once dominated the east coast of Australia [21]. This is the pattern for other faunal lineages of Gondwanan origin, such as *Terrisswalkerius* earthworms [24], rainforest possums [34], and myobatrachid frogs [35]. We also expect species distributions to be strongly influenced by the paleogeography of the AWT, with most species having restricted northern or southern distribution across the Black Mountain Corridor [23]. Finally, in addition to high temperature, we expected rainfall seasonality to be a key driver of species distributions.

#### **2. Materials and Methods**

#### *2.1. Study Region*

The Australian Wet Tropics bioregion is a World Heritage Area in far northeastern Australia covering approximately 36,000 km<sup>2</sup> (20–15◦ S and 147–145◦ E) consisting of approximately 12,000 km2 of rainforest. Rainfall is highly seasonal, with 75–90% of the annual 2000–8000 mm rainfall occurring in the wet season between November and April. Elevations above 1000 m a.s.l. can also receive up to 66% of monthly water input from cloud stripping [30]. Rainforests of the AWT harbour a distinctive ant fauna that is highly disjunct from that of surrounding savanna [36,37]. The fauna has strong affinities with that of Indo-Malayan rainforests, containing many genera that are rainforest specialists and whose Australian distributions are restricted to North Queensland. The fauna notably lacks arid-adapted taxa, such as *Iridomyrmex* and *Melophorus*, that dominate ant communities of Australia's open sclerophyll habitats [38]. Tropical rain forests are generally regarded as supporting the world's richest ant faunas [39,40], but Australia's rainforest ant fauna is relatively depauperate [36], in striking contrast to its exceptionally rich savanna fauna [41].

#### *2.2. Ant Sampling*

Collections of *Anonychomyrma* were assembled from three sources. The first was the study of [9], which was based on surveys at 26 of the long-term biodiversity monitoring sites established by Stephen E. Williams at James Cook University [35]. These sites were distributed across six subregions, covering the full latitudinal and elevational range of the AWT and representing approximately 94% of the available environmental space in the region [34]. The subregions spanned approximately 500 km from north to south and sites were placed approximately every 200 m along the elevation gradient as follows: Finnegan (200, 500, 600, 800 m a.s.l.), Windsor (900, 1100, 1300 m a.s.l.), Carbine (100, 400, 600, 800, 1000, 1200 m a.s.l.), Lamb (700, 900, 1100 m a.s.l.), Atherton (200, 400, 600, 800, 1000 m a.s.l.), and Spec (350, 600, 800, 1000 m a.s.l.) (Figure 2). Sampling occurred within six plots at each site separated by a distance of 200 m; only three plots were located at each of the 350 m a.s.l. site at Spec and 100 m a.s.l. site at Atherton due to limited rainforest cover. Each site was sampled using a combination of pitfall trapping, litter extractions, and baiting (both on the ground and on trees at a height of 1.5 m). At Windsor, Carbine, Atherton, and Spec, sampling was conducted on three occasions from 2011 to 2013, covering two wet seasons (November–January) and one dry season (June–September). At Finnegan and Lamb, sampling occurred only during one wet season. The study also included ants collected in pitfall traps during previous beetle surveys at the Windsor, Carbine, Atherton, and Spec sites [42].

The second source was from [32], a study involving additional surveys at 15 of the above sites: Finnegan (200, 500, 700 m a.s.l.), Windsor (900, 1100, 1300 m a.s.l.), Carbine (100, 600, 1000, 1200 m a.s.l.), and Atherton (200, 400, 600, 800, 1000 m a.s.l.). At each of these sites, either two (1 site), three (4 sites), four (4 sites), or five (6 sites) trees were sampled, totalling 60 trees. Trees were at least 50 m apart and were chosen for surveying based on size and climbing accessibility. Trees were sampled using tuna bait traps and accessed using the single-rope climbing technique. At each tree, five bait traps were set on the ground and at every 3 m above ground to the maximum accessible height of the tree, which ranged from 15 to 27 m. Traps were set in the morning and collected 2–3 h later. Finnegan, Windsor, and Carbine were surveyed from October to December 2012 and Atherton was surveyed from December 2017 to February 2018.

The third source was data from additional surveys conducted in 2018–2019. Carbine sites were resampled in February 2018 using the same methodology as [32], and from June to October 2019, but with surveys during both day and night. At each site, five trees (at least 50 m apart) were surveyed once in the daytime and once at night. Day and night surveys of the same trees were not conducted sequentially and were at least a day apart to allow resumption of normal ant activity in the tree in case of disruption. Daytime surveys started at 10:00 h, night-time surveys started at 21:00 h, and baits were

collected three hours later. Ants were also opportunistically hand collected immediately adjacent to baited vials. In addition, a collection was made over five days during the same period at the Daintree Rainforest Observatory (100 m a.s.l.) placing baited vials at 25–30 m height at ten trees (accessed by the canopy crane) for one hour and conducting one-hour hand searches on the ground.

#### *2.3. Gene Barcoding*

DNA was extracted from foreleg tissue and sequences were obtained for 97 specimens of *Anonychomyrma* collected during the study (Table S1). DNA extraction and CO1 sequencing were conducted through the Barcode of Life Data (BOLD) System (for extraction details, see http://ccdb.ca/ resources). Each sequenced specimen was assigned a unique identification code that combined the batch within which it was processed and its number within the batch (e.g., ANONC006-20), and all specimens were labelled with their respective BOLD identification numbers in the ant collection held at the CSIRO Tropical Ecosystems Research Centre in Darwin.

DNA sequences were checked and edited in MEGA 7 [43]. Sequences were aligned using the UPGMB clustering method in MUSCLE [44], and translated into (invertebrate) proteins to check for stop codons and nuclear paralogues. The aligned sequences were trimmed accordingly, resulting in 657 base pairs. A sequence from the *anceps* complex of *Iridomyrmex* (sp. A, IRIDX092-18 from [45]) was used as an outgroup for rooting the gene tree. Tree inference by maximum likelihood was conducted through the IQTREE web server (http://iqtree.cibiv.univie.ac.at/; [46]) using ultrafast bootstrap approximation [47]. IQTREE has been shown to be a robust algorithm for tree inference that compares favourably with other methods [48]. Model selection was inferred using a 3-codon partition file and linked branch lengths with the AutoMRE 'ModelFinder' function to find the best-fit model for tree inference [49]. Trees were viewed and edited in FigTree v1.4.3 [50] and annotated using Photoshop CS5.1.

#### *2.4. Species Delimitation*

There is no specific level of CO1 divergence that can be used to define a species, but the level of CO1 variation within ant species is typically 1–3% [51]. However, some ant species can show substantially higher variation (for example, [51,52]), and in other cases two clear species can show no CO1 differentiation (for example, [53]). Some ant species from other genera are known to have workers that are virtually identical morphologically, and they can only be separated by detailed morphometric analysis or through reproductive castes [54]. When delimiting species, we focused on morphological differentiation between sister (i.e., most closely related) CO1 clades, considering all available samples from the same collections as sequenced specimens. A full set of voucher specimens of recognised species was held at the CSIRO Tropical Ecosystems Research Centre in Darwin and a duplicate set in the ant collection at James Cook University.

#### *2.5. Patterns of Abundance, Diversity, and Endemism*

We used the total number of records in survey plots including repeat surveys (*n* = 525 surveys) across all sites (*n* = 23) as a measure of abundance. Species were then ranked by abundance to document the species abundance pattern. We excluded the Daintree subregion from all analyses of diversity as there was only one site sampled. *Anonychomyrma* was not recorded from the Lamb subregion and these sites were also excluded from the analysis. To assess latitudinal diversity patterns, we assessed variation in species richness among subregions by plotting the mean number of species across sites per subregion along with the total number of surveys per subregion. We likewise plotted the number of subregion endemics per subregion. To assess elevation diversity patterns, we investigated differences in species richness among elevations by pooling across the five subregions and plotting the mean number of species across plots per elevation along with the total number of surveys per elevation. We tested for a correlation in the number of species with survey effort for both subregion diversity and elevation diversity using Pearson's product-moment correlation test. To assess the vertical distributions (foraging activity from ground to canopy) of species, we considered only the vertical (arboreal) surveys

conducted by [32]. We selected species for which there were two or more survey records and ten or more sample records (as there were 5 samples per vertical height band). For each of the resultant eight species, we plotted the number of sample records in a 3-m band divided by the total number of sample records for that species, to provide a relative proportion of occurrence in each 3-m band from ground to canopy.

#### *2.6. Species Distribution Modelling*

We used species distribution models as an exploratory tool to three purposes: first, to investigate the potential climatic drivers of species distributions; second, to identify potential additional areas of suitable habitat that were not sampled in this study; and third, to explore how species-predicted distributions correlate with historical patterns of rainforest refugia. We included the Daintree subregion in our species records and conducted distribution modelling over only the rainforest regions of the AWT (~12,000 km2) as we were interested in predicting the distribution of *Anonychomyrma* species within their primary habitat of rainforest vegetation. Species distribution models based on species occurrences and climatic data were derived using MaxEnt, a maximum entropy algorithm, in the program MaxEnt using default settings (version 3.4.1, [55]). We acknowledge that the use of default settings in MaxEnt may not lead to the optimal model in all cases, but given the exploratory purpose of modelling here, these settings were sufficient for our stated goals [56]. MaxEnt performs well in modelling species with small sample sizes if the ecological niche (e.g., environmental tolerance) of those species is sufficiently covered by sampling and they have small geographic ranges [57,58]. However, MaxEnt can perform poorly when species with small number of records are geographically widespread or have a wide environmental tolerance that is not sufficiently covered by sampling [59,60]. For example, [57] reported increasing model performance for species with increasing sample size, but that useful models are still produced for rare and narrow range species with 5–10 records, while [58] strongly emphasised that sample size is relevant to the total area over which modelling is being conducted. They found 14 records to be a minimum for species prevalence (fraction of raster cells occupied) of 0.1 over an area that encompassed most of tropical Africa, which is a far larger area than we model here [58]. To explore the issue of sample size, we first filtered species for ≥ 5 location points (not counting repeated surveys and counting a record in a survey plot as one occurrence/location point), which resulted in a subset of nine species. The number of location points ranged from 96 (*A. gilberti*) to 6 (sp. H). To investigate the relationship between sample size and sample coverage of ecological niche, we looked at how much of the available environmental space of the rainforest region of the AWT was covered by each species' location points. We regressed mean annual temperature against mean annual precipitation to create the environmental space of all survey locations and plotted each species' location points within that plot. Species with low (< 15) numbers of points were highly concentrated in environmental space, indicating that the number of location points sufficiently captured each species' climatic niche volume (Figure S1). As such, our species with low samples sizes were very likely to have narrow environmental niches and therefore could be modelled with the number of location records available. However, we acknowledge that caution should still be applied in interpreting the models with <15 location points [61].

We used the accuCLIM climate variables derived by [62] that were based off a distributed network of microclimate loggers in the AWT and produced using a boosted regression tree approach that statistically downscaled existing coarse weather layers to fine-scale weather layers at 250 m<sup>2</sup> resolution for the AWT region comparable with current best-practice climate layers (e.g., ANUCLIM, [63]). These layers are highly accurate in relation to regional topography and vegetation [62]. We clipped all climate variables to rainforest extent in the AWT. We selected nine of the seventeen climate variables that were most ecologically relevant [42,64–66]. We looked for collinearity in variables by looking at all pairwise interactions between continuous covariates using Pearson's correlation coefficient. Variables with an *R*<sup>2</sup> value of >0.7 were considered for removal [67]. Variance inflation factors were then calculated and any variables that had a value >10 were excluded from further analysis using the package *usdm* in R [68]. A final set of five variables was selected as follows: maximum temperature

of the warmest month, temperature seasonality, rainfall of the wettest month, rainfall seasonality, and isothermality (an indicator of temperature variability: mean diurnal temperature range divided by annual temperature range; [62]). In MaxEnt, 10 replicates were used per species model with cross validation and 1000 iterations [69]. All models presented had values for the area under the receiver operating characteristic curve (AUC) greater than 0.9 and therefore performed adequately [70]. We then calculated summed habitat suitability for the rainforest extent by summing all nine species habitat suitability scores across all pixels and standardizing each pixel's value to a summed habitat suitability score between 0 and 1.

#### **3. Results**

#### *3.1. The Anonychomyrma Fauna*

We recognised 22 species of *Anonychomyrma* (sp. A–C and E–U, along with *A. gilberti* and *A. malandana*) among our samples, nine of which were recorded at a single subregion (Table 1). All but one of our recognised species were successfully sequenced. The CO1 tree (Figure 3) indicated that species C (Figure 4A,B) was the most phylogenetically divergent, with no close relatives (13–20% divergence from other species). It was a relatively large, gracile, and somewhat polymorphic species with an angular propodeum, and was common at high elevation at Carbine and Windsor. *Anonychomyrma malandana* was also indicated as highly distinctive phylogenetically (> 12% divergence from all other species) as well as being highly distinctive morphologically—it was an extremely shiny species with a globose head and very long antennal scapes (Figure 4C,D). It was a high-elevation species that was common at Windsor. Species H and sp. G were also phylogenetically distinct (Figure 3), both with > 11% divergence from all other species. Species H was a small, nondescript species with a biconvex mesosoma and relatively short scapes (Figure 5A,B) and was morphologically very similar to several unrelated species (see below); it was recorded at Carbine and Windsor at the highest elevations (1200 m and 1300 m, respectively). Species G had a short, prominently rounded propodeum, and its scape and first gastric tergite atypically lacked erect hairs (Figure 5C,D). These characters were shared by the smaller sp. J (Figure 5E,F), but despite their morphological similarity (they were considered conspecific in [9]), these species were widely separated in the CO1 tree (Figure 3). They had overlapping geographic distributions at Atherton: sp. G was recorded from all five subregions, whereas sp. J was recorded exclusively from Atherton (Table 1). Species J was shown as being most closely related to sp. I (Figure 6A,B) and sp. S (Figure 6C,D; Figure 3). Compared with sp. I, sp. S had a substantially larger head with a markedly concave occipital margin, more angular occipital corners, and shorter scapes. Both species were recorded only from Atherton and only from high elevation (800 m and 1000 m; Table 1).

**Figure 3.** *Cont.*

**Figure 3.** *Cont.*

**Figure 3.** (**A**–**C**) Maximum likelihood CO1 tree of 97 specimens of *Anonychomyrma* from the Australian Wet Tropics, with a specimen of the *anceps* complex of *Iridomyrmex* as an outgroup.

**Table 1.** The subregion locations and elevation range (maximum–minimum elevation) and number of survey points of 22 species of *Anonychomyrma* in the Australian Wet Tropics, species are ordered from their north to south subregion distribution, and then by number of survey points, *n* indicates number of survey points. BMC position indicates the positions of each subregion north or south of the Black Mountain Corridor, a well-known biogeographic barrier. The full elevation range of rainforest habitat is shown in brackets for each subregion.


**Figure 4.** Head and lateral views of *Anonychomyrma* sp. C (**A**,**B**) and *Anonychomyrma malandana* (**C**,**D**).

**Figure 5.** Head and lateral views of *Anonychomyrma* sp. H (**A**,**B**), sp. G (**C**,**D**), and sp. J (**E**,**F**).

**Figure 6.** Head and lateral views of *Anonychomyrma* sp. I (**A**,**B**) and sp. S (**C**,**D**).

Thirteen of the remaining fourteen species formed two clades. The first contained *A. gilberti* and species A, L, M, and R, which had 6–8% mean CO1 divergence among them (Figure 3). *Anonychomyrma gilberti* occurred in all subregions, but the other four species were each recorded in a single subregion (Table 1). Species A, L, and M (Figure 7A–F) had the general appearance of *A. gilberti* (Figure 1), but without such strongly golden pubescence or deeply V-shaped occipital margin; indeed, sp. L had relatively weak pubescence and rounded occipital corners (Figure 7C,D). Species R was morphologically very different to these taxa—it was much smaller, lacked any golden pubescence, had short scapes, and the occipital margin was only feebly concave. Despite its location on the CO1 tree, it seemed to be more closely allied to members of the second clade spanning sp. P to sp. K (Figure 3). Indeed, it seemed morphologically indistinguishable from sp. N from Atherton (Figure 8A,B). All species within this second clade were similar morphologically, and CO1 divergence among them was often < 2%. Compared with sp. N, sp. T (known only from Windsor) had a more conspicuously concave occipital margin and more angular occipital corners, and in sp. P the occipital margin was even more deeply concave. Species U had a narrow head, short scapes, and a prominently rounded propodeum. The propodeum was also prominently rounded in sp. F (Figure 8C,D), but the head was markedly broader and scapes were longer. Species Q had a flattened propodeum and short scapes. Species K (Figure 8E,F) had a very similar appearance to that of sp. N (Figure 7A,B), but the two taxa were very distinct genetically (11% mean CO1 divergence).

**Figure 7.** Head and lateral views of *Anonychomyrma* sp. A (**A**,**B**), sp. L (**C**,**D**), and sp. M (**E**,**F**), which are all closely related to *Anonychomyrma gilberti* (Figure 1).

**Figure 8.** Head and lateral views of *Anonychomyrma* sp. N (**A**,**B**), sp. F (**C**,**D**) and sp. K (**E**,**F**).

The final sequenced species was sp. E, which had a biconvex mesosoma, broad head, and long scapes (Figure 9A,B). It occurred exclusively at Windsor (Table 1). The species that was not sequenced, sp. O, was highly distinctive—it was very small, with a narrowly rectangular head and very short scapes (Figure 10). It appeared to be closely allied to *Anonychomyrma minuta* (Donisthorpe) from New Guinea. It was recorded only at Atherton, at mid elevation (400, 600, and 800 m sites).

**Figure 9.** Head and lateral views of *Anonychomyrma* sp. E (**A**,**B**).

**Figure 10.** Head and lateral views of *Anonychomyrma* sp. O (**A**,**B**).

#### *3.2. Patterns of Abundance, Diversity, and Endemism*

*Anonychomyrma gilberti* (185 records) and sp. G (73 records) dominated the samples, collectively representing 81.3% of all species records (Figure 11). Nine species were recorded in only one survey (Figure 11). The number of species within a subregion ranged from 0 at Lamb Range to 11 at Atherton. There was no relationship between survey effort and mean number of species per either subregion (Pearson's correlation = −0.02, *p* = 0.9; Figure 12A,B) or elevation (Pearson's correlation = 0.06, p = 0.8; Figure 13). Windsor (1100–1300 m a.s.l.) had the highest mean species richness per elevational site (5.7 ± 0.7 SE; Figure 12A), and the second highest number of subregion endemics with 4 out of its 10 species (Figure 12B). Atherton (100–1000 m a.s.l.) had the highest number of subregion endemics with 7 endemics out of its total of 11 species (Figure 12B). Carbine, despite having the greatest elevation range (100–1200 m a.s.l.) and the highest survey effort, had similar mean species richness to other low diversity subregions and had only 2 endemic species out of the 8 total species recorded there (Figure 12A,B). Overall, there was high subregion endemism with 16 (72%) species recorded in only one subregion (Table 1). However, all four species (*A. gilberti*, *A. malandana*, sp. G, and sp. K) that occurred in more than two subregions had distributions spanning the BMC (Table 1). The mean number of species increased with elevation peaking at 1100 m a.s.l. and slightly declining again at 1200 and 1300 m a.s.l. (Figure 13).

**Figure 11.** Ranked abundance as total number of samples for each species from pooled data for all survey sites across five subregions in the Australian Wet Tropics.

**Figure 12.** Variation in diversity and endemism among subregions. Subregions labelled from left to right along the north to south latitudinal gradient. (**A**) Mean (±SE) number of *Anonychomyrma* species per elevational site as grey bars and number of surveys per subregion as orange symbols. (**B**) Number of subregion endemics per subregion, and in brackets, the total species richness per subregion. Orange symbols are number of surveys per subregion.

**Figure 13.** Diversity across elevation, pooling across five subregions. Mean number of *Anonychomyrma* species (±SE) per elevation as gray bars and number of surveys per elevation as orange symbols.

All of the eight species recorded during vertical (tree) surveys nested arboreally (L. Leahy, pers. obs.), with most species foraging high into the canopy. However, there was variation in the way different species used the vertical gradient (Figure 14). Three species, *A. gilberti*, *A. malandana*, and sp. G foraged on the ground as well as arboreally. The two most commonly sampled species, *A. gilberti* and sp. G, were most evenly distributed in their foraging along the vertical gradient, with relatively equal foraging from ground to canopy. These two species were also the most geographically and elevationally widespread species (Table 1). *Anonychomyrma malandana* and sp. C also showed relatively even foraging across the vertical gradient. The less commonly sampled species showed increasing foraging concentration in the upper parts of the tree from the subcanopy (9–12 m) up to the high canopy (18–27 m; Figure 14).

#### *3.3. Species Distribution Modelling*

Maximum temperature of the warmest month was the most important predictor for seven species and the second most important predictor for the other two species modelled (Table 2). For all species, there was a negative relationship between maximum temperature of the warmest month and predicted habitat suitability, but the temperature at which habitat suitability declined (ranging from 26 ◦C for the high-elevation-restricted sp. H to 34 ◦C for the widespread *A. gilberti* and sp. G) and the rate of decline differed substantially (Figure S2). Rainfall of the wettest month was the second most important predictor variable for *A. gilberti,* sp. E, *A. malandana*, and sp. K, being strongly negatively related to habitat suitability in all cases (Table 2, Figure S3). Rainfall seasonality was also important in models, but generally ranked as the second or third most important variable, except for two Spec endemics, sp. A and sp. M, for which rainfall seasonality was the most important predictor. These two Spec endemic species had a positive relationship between rainfall seasonality and habitat suitability, whereas all other species had higher habitat suitability at intermediate levels of rainfall seasonality or had a negative relationship with rainfall seasonality (Figure S4).

**Figure 14.** Vertical distribution of *Anonychomyrma* species sampled in 105 surveys (out of 525 total surveys) from ground to canopy at 15 elevation sites at four subregions in the AWT. Pooled occurrence data from four subregions across the Australian Wet Tropics, with species ordered by number of survey occurrences from left to right. Each tile represents a 3-metre vertical band, coloured by proportion of total sample occurrence for each species. Numbers beneath species labels are total number of sample occurrences for that species (total samples = 3745). The maximum potential tree survey height for each species (based on occurrence in each surveyed tree) is: *gilberti* = 27 m, sp. G = 27 m, sp. C = 24 m, *malandana* = 24 m, sp. K = 24 m, sp. H = 27 m, sp. J = 21 m, sp. I = 24 m. Red represents high proportion of samples and yellow represents low proportion of samples.

Distribution modelling predicted that six out of nine species modelled had potential distributions limited to high elevations in either the north (Windsor and Carbine) or far south (Spec) (Figure 15). Summed rainforest habitat suitability (summing each grid cell for suitability for the nine modelled species) was highest at high elevation (Figure 15A), which was consistent with the occurrence of highest species richness (Figure 13). Species-specific habitat suitability generally matched surveyed occurrences, suggesting that actual distributions were close to the potential distributions based on climatic niches (Figure 15). However, our models indicated that highly suitable (≥0.75 probability of occurrence) environments occurred in the southern region for sp. C, sp. E, and sp. H, despite these species being recorded only in the northern subregions in our study. We note that sp. E and sp. H had <10 location points and so these models should be treated with some caution. None of the species modelled had a potential distribution restricted to the centre (Atherton or Lamb) or to the northern lowlands (Figure 15B–J). However, many of the species that were not sufficiently common for modelling were recorded from Atherton, and one species (sp. B) with one record was recorded only at Finnegan at 200 m a.s.l. (Table 1). Models predicted *A. gilberti* and sp. G to have the greatest potential geographic range, followed by sp. K (Figure 15). Distribution maps showed there was a high probability that these three species occur in other rainforest subregions in the AWT not surveyed here, including Bartle Frere, Bellenden Ker, Lamb, and Kirrama (Figure 15; for map of all AWT subregions, see [18]). Although sp. K had only seven location points for modelling, these were all high elevation (Table 1) and our exploration of the environmental niche space indicated a narrow environmental tolerance; we are therefore confident that sp. K is elevationally restricted despite being geographically widespread.

**Figure 15.** (**A**) Summed rainforest habitat suitability across all species in a grid cell based on the MaxEnt species distribution model for each species. (**B**–**J**) Modelled rainforest distributions for nine species of *Anonychomyrma* showing habitat suitability in the AWT bioregion. In all figures, blue is low habitat suitability, red is high habitat suitability, and grey background is the outline of the AWT region.

**Table 2.** Species distribution model outputs for nine species showing the three highest ranking environmental variables and their permutation importance in the model. AUC = area under the curve, a value between 0 and 1 indicates model performance, and higher AUC scores indicate better model fit. *n* = number of unique survey location points for each species. Max T hottest month = Maximum temperature of the hottest month, Rain seasonality = Rainfall seasonality, Rain wettest month = Rainfall of the wettest month, Temp seasonality = Temperature seasonality.


#### **4. Discussion**

The Australian Wet Tropics hosts a diverse *Anonychomyrma* fauna, whose species show strong geographic patterning. Our genetic and morphological analysis delimited 22 species from our samples, far more than previously recorded in the region. Only two of the species appeared to be described. One of these, *A. gilberti*, is an ecologically dominant ant in rainforest canopies throughout the region [9,10]. The other is *A. malandana*, which was originally described from Malanda on the Atherton Uplands, but in our samples, it was recorded in Windsor and Carbine as well (Figure 2). Our samples were from a restricted range of sites in the region and many of our species are known from just one or a few records, which suggests that many additional species likely occur in the AWT.

The distribution and diversity of *Anonychomyrma* within the AWT suggest that this Gondwanan genus has retained a preference for cooler climates. This is evidenced by the concentration of species at elevations above 900 m a.s.l., with only three of the 22 species being recorded at elevation below 300 m a.s.l. The uplands of Windsor and Atherton are centres of particularly high diversity and endemism. Species distribution models also supported this pattern, and summed habitat suitability showed a general trend that species had high habitat suitability in upland areas that are also the main refugial areas of long-term rainforest stability (for AWT map of paleo-geological habitat stability, see [21]). Such a preference for high elevation is in contrast to the distributional patterns of most other ant taxa in the AWT, and more broadly across the tropics, where ant diversity peaks at low to mid elevations and drops dramatically at high elevations [9,71–73]. A similar distributional pattern for

*Anonychomyrma* occurs in Papua New Guinea, where six of eight species surveyed were restricted to high elevation (≥900 m a.s.l.) [5]. Similarly, in neotropical sites of Central America, several lineages of ants that have a north American temperate climate origin, such as *Temnothorax* [74], are restricted to montane rainforests [71].

All eight species surveyed for vertical distribution were strongly arboreal, often foraging from the understorey up into the high canopy. Only a few species were recorded also foraging on the ground. This suggests that arboreality is a strong trait in the tropical clades of *Anonychomyrma* in Australia. Many of the species were noted to nest arboreally in trunk cavities (L. Leahy, pers. obs.). The two most abundant species (*A. gilberti* and sp. G) had the broadest vertical ranges and also had the broadest geographic and elevation ranges, occurring in all subregions and at all elevations. Species with foraging restricted to the canopy tended to have more restricted subregional and elevational distributions. These findings are consistent with the general pattern for AWT ants that vertical niche breadth is positively related to elevation range size [32]. This pattern also follows the general macroecological rule that locally common species tend to be widely distributed [75].

Species distribution modelling of nine species strongly supported maximum temperature during the warmest month as an important predictor of distribution, indicating that high temperature is likely to be an important factor in limiting species distributions. The two most abundant species (*A. gilberti* and sp. G) were more tolerant of higher maximum temperatures, occurring in all subregions and across the full elevation range. In our modelling, the estimated maximum temperature at which predicted habitat suitability declines was based on a downscaled climatic layer specific to the AWT and was therefore more accurate than a climate layer based on the limited coverage of weather stations in the region [62,76]. However, it is important to note that climate circulation models do not account for vertical climatic gradients, whereby the canopy can experience temperatures several degrees warmer than the forest floor [31,77]. This is an important consideration because the diminutive stature of ants places them close to surface temperatures, which can be far hotter under direct solar radiation, particularly in the canopy [78]. Therefore, given the arboreal habits of *Anonychomyrma*, maximum temperature exposure is likely to be much higher than the predicted estimates from our species distribution models.

For several species, there was a negative relationship between habitat suitability and rainfall of the wettest month and rainfall seasonality. This could explain the relative depauperate species representation in the uplands of Carbine, where elevations above 1000 m a.s.l. have high rainfall seasonality [26] and very high rainfall during the wet season. Windsor, in comparison, which had the highest species richness, sits in the rain shadow of Carbine and therefore has less rainfall with less seasonality [26]. Similarly, Atherton, also with high species richness and endemism, has relatively low seasonality of both rainfall and temperature. Spec in the extreme south of the AWT had low species richness and may be too dry for most species as its dry seasons are substantially more severe than in the other subregions [26]. Rainfall seasonality is a strong driver of insect abundance throughout the tropics [27,79,80], including in the AWT [26]. High rainfall seasonality is likely to limit ant distributions through its effect on food availability, particularly honeydew from homopteran insects.

A number of our distribution models predicted species to occur in subregions that were not sampled in this study. Of particular note, although we did not record any *Anonychomyrma* species in the Lamb subregion, several species were predicted to occur there based on the availability of a suitable habitat. The absence of *Anonychomyrma* at Lamb was therefore likely to be at least partly an artifact of low sampling intensity at that subregion. Other subregions that had high predicted suitability, but were not sampled here, include Bartle Frere and Bellenden-Ker. These are the two highest mountain ranges in the Australian Wet Tropics and have not been systematically sampled for ants. They are centres of diversity and subregion endemism for several Wet Tropics species, including the Gondwanan plant species *Eucryphia wilkei,* which is endemic to elevations above 1500 m a.s.l. on Bartle Frere, and whose congeners occur only in Tasmania and Chile [35,81]. It is highly likely that the tops of these two mountains harbour additional species of *Anonychomyrma*.

There was limited evidence of a biogeographic barrier across the Black Mountain Corridor. Four out of six species that occurred in more than one subregion had distributions spanning north and south of the BMC. Invertebrate taxa, such as dung beetles [82], earthworms [24], flightless insects [83], and schizophoran flies [84], that show a biogeographic divide between northern and southern distributions, all have limited dispersal ability compared to ants with winged queens (such as that occurs in *Anonychomyrma*), which can be dispersed by wind over long distances, including across biogeographic barriers (e.g., [85]).

#### **5. Conclusions**

The biogeographic pattern of *Anonychomyrma* in the AWT supports the hypothesis that this genus has a Gondwanan origin and has radiated from a Australian southern temperate distribution into the tropics [6]. *Anonychomyrma* has persisted and diversified in rainforest areas of cool climate with only a few species reaching into the hotter lowland environments and spreading across multiple subregions [18,86,87]. This mirrors the distributional pattern of the genus in Papua New Guinea [5] and we would expect a similar pattern to occur in other parts of the genus' tropical range.

Tropical species restricted to upland habitats with a limited geographic range are considered highly vulnerable to anthropogenic climate change [88,89]. The arboreal foraging habits of *Anonychomyrma* may provide some options to behaviourally regulate climate exposure, given the strong thermal gradient along trees [90–92]. However, the increasingly warm and more seasonal climate that is forecasted for the Australian Wet Tropics [34], along with significant predicted shifts of habitat from rainforest to drier vegetation types [65], places these Gondwanan mountain-top relics under increasing threat.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/1424-2818/12/12/474/s1, Table S1: Collection locations of specimens from the Anonychomyrma group that were CO1-barcoded for this study, Figure S1: Environmental space occupied by nine species used in species distribution modelling based on our sampling across the AWT, Figure S2: Species distribution modelling response curves showing Anonychomyrma species responses to accuCLIM 1996–2015 Variable 05: Maximum temperature of the warmest month, Figure S3: Species distribution modelling response curves showing Anonychomyrma species responses to accuCLIM 1996–2015 Variable 13: Rainfall of the wettest month, Figure S4: Species distribution modelling response curves showing Anonychomyrma species responses to accuCLIM 1996–2015 Variable 15: Rainfall seasonality.

**Author Contributions:** A.N.A. conceived of the study, managed the TERC collection, undertook the analysis of the CO1 and morphological data, and contributed substantially to the writing of the paper, L.L. collected specimens in the field, conducted all other analyses, took specimen photos, and wrote the first draft of the manuscript. B.R.S. collected specimens in the field and contributed to the writing of the paper, S.E.W. provided funding for field work and contributed to the writing of the paper. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was supported by multiple grants from L.L.: The Explorer's Club, Wet Tropics Management Authority, Skyrail Rainforest Foundation, Holsworth Wildlife Research Endowment—Equity Trustees Charitable Foundation, and Ecological Society of Australia. Funding from S.E.W.: National Environmental Research Program—Tropical Ecosystem Hub from the Australian Government. L.L. is also supported by a PhD scholarship from the Australian Government. Field research was conducted under permit WITK15811415 issued by the Queensland Government of Australia.

**Acknowledgments:** We thank the many different Aboriginal traditional owner groups on whose lands the study was conducted. We thank Somayeh Nowrouzi for collecting specimens analysed in this study and all volunteers involved in collections. We thank Stefanie K. Oberprieler for help with analysis of the CO1 data.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
