EVOLUTIONARY BIOLOGY OF THE AUSTRALIAN CARNIVOROUS MARSUPIAL GENUS ANTECHINUS
Thomas Yeshe Mutton B. App. Sci. (Hons)
School of Earth, Environmental and Biological Sciences Science and Engineering Faculty
Queensland University of Technology
Submitted in fulfilment of the requirements for the degree of Doctor of Philosophy
i
Keywords
Antechinus, Australia, Australian mesic zone, biogeographical barriers, biogeography, breeding biology, congeneric competition, conservation, genetic structure, Dasyuridae, dasyurid, life-history, mammal, mammalian, Miocene, molecular systematics, Phascogalini, phylogenetics, phylogeography, Pleistocene, Pliocene, population genetics, Queensland, systematics (Davison )
ii
Abstract
Antechinus is an Australian genus of small carnivorous marsupials. Since , the number of described species in the genus has increased by % from ten to fifteen. The systematic relationships of these new species and others in the genus have not been well resolved and a broad phylogeographic study of the genus is lacking. Moreover, little ecological information is known about these new species. Therefore, the present study examined the evolutionary biology of Antechinus in two complimentary components. The first component aimed to resolve the systematics and phylogeography of the genus Antechinus. The second component, at a finer spatiotemporal scale, aimed to improve understanding of the autecology, habitat use and risk of extinction within the group, with a focus on the recently named buff-footed antechinus, A. mysticus and a partially sympatric congener, A. subtropicus.
To date, powerful (>two gene) molecular studies have only included eight Antechinus species. Here, the first comprehensive, multi-gene phylogenetic analysis of the genus using six genes was provided, incorporating all known species and subspecies of Antechinus. Four main lineages of Antechinus were reconstructed: ( ) The dusky antechinus (A. arktos, A. swainsonii, A. vandycki, A. mimetes) and A. minimus; ( ) A. godmani; ( ) A. agilis, A. stuartii and A. subtropicus; ( ) A. argentus, A. mysticus, A. adustus, A. flavipes, A. leo and A. bellus. The inclusion of A. adustus in lineage was surprising, as it was previously believed a member of lineage . One species, A. stuartii, was not monophyletic and may be more appropriately classified as two species. Timing of cladogenesis was estimated for all species of Antechinus, permitting an evolutionary scenario to be posited for the group.
The mesic zone is the ancestral biome of Australia and, by most measures, the zone of highest taxonomic diversity. However, the south-east region of the
iii
zone has been identified as lacking in mammalian phylogeographic studies. Sequencing of the mitochondrial Cytochrome-b (Cytb) gene was therefore undertaken broadly throughout the geographic range of all species of Antechinus, which collectively occur across the entire zone. High levels of intraspecific phylogeographic structure, largely corresponding to probable Pleistocene biogeographic barriers and comparable to the genetic divergence between named Antechinus subspecies, were revealed in a number of species. For example, the broadly distributed A. flavipes flavipes contained two divergent clades which appear to have been separated by the Great Dividing Range. Within the south-east Australian A. agilis, four phylogeographic clades were identified, each of which appears to have been divided by putative biogeographic barriers. Antechinus mimetes mimetes showed strong divergence between north-east NSW and southern populations. Victorian populations of this species also appear to have been divided by the East Gippsland Basin which may also have divided A. agilis clades and other taxa in the region. Considerable divergence was found between south-east Queensland populations of A. mysticus and the only known population in mid- east Queensland.
The buff-footed antechinus, A. mysticus, was only recently discovered and consequently it is not well understood. It has been suggested that this species utilises a broad range of habitats which over time may result in highly connected populations at a regional scale. This hypothesis was tested using a population genetic approach. Nine microsatellite loci were genotyped for six populations of A. mysticus, sampled from throughout its known range (Eungella in mid-east Qld [ME Qld] and south-east Qld [SE Qld]). Comparative data were sought from four populations of A. subtropicus, which can occur in sympatry across a similar area of SE Qld. Antechinus subtropicus is known to utilise altitudinal vine and rainforest communities which are more fragmented; consequently, it was expected to show greater population differentiation than A. mysticus. Two populations of A. mysticus were revealed to be deeply divergent from most SE Qld populations examined: Eungella (ME
iv
Qld), which is in geographic isolation (> km) from the other A. mysticus populations, and Cooloola (SE Qld). Cooloola is the northernmost known SE Qld population of A. mysticus and shares a mtDNA Cytb haplotype with other SE Qld populations. Yet this population had low microsatellite allelic diversity and appears to be in very low abundance, suggesting divergence is driven by population isolation; it may be at extinction risk. Genetic structure between the other SE Qld populations was also higher than expected, similar to A. subtropicus populations. This suggests A. mysticus is patchily distributed and may be limited to moderate altitude and wetter environments in SE Qld.
The first multi-year ecological study of breeding, growth and movement of A. mysticus was also undertaken. Over a two year period, monthly capture-mark- recapture data from two proximate (~ . km apart) sites were collected. At one site, the subtropical antechinus, A. subtropicus, occurred in sympatry. Antechinus mysticus followed the synchronous semelparous (suicidal breeding) breeding strategy seen in all Antechinus. Males were last caught in August at both sites and females appear to give birth in September. Average movement by A. mysticus was comparable to other similar-sized congeners. Antechinus subtropicus occurred at lower density although not lower body mass than proximate conspecifics, at this relatively low altitude study site. Competition from A. subtropicus may affect A. mysticus, as A. mysticus have a higher body mass and males moved further when not in sympatry with A. subtropicus.
Taken together, the present study permitted an investigation of spatio- temporal influences on Antechinus evolution, from the distant Miocene to a very recent inter-generational scale. Genetic diversity was found to be high in Antechinus at all levels investigated which suggests this mesic genus has experienced fragmentation throughout its existence. Several Antechinus species are at risk of extinction. In light of the observed population structuring and purported history of fragmentation, Antechinus must be closely monitored
v
in a future likely to include further habitat degradation, as temperatures and human abundance increases across Australia.
vi
Table of Contents
Keywords ...... ii Abstract ...... iii Table of Contents ...... vii List of Figures ...... ix List of Tables ...... xiii Statement of Original Authorship ...... xvi Acknowledgements ...... xvii Chapter 1: General introduction ...... 1 1.1 Taxonomy and systematics of Antechinus ...... 1 1.2 Biogeography...... 3 1.3 Life history and breeding biology ...... 4 1.4 Movement and home-range ...... 6 1.5 The present study ...... 7 Chapter 2: Molecular systematics and evolution of the genus Antechinus .. 11 2.1 Introduction ...... 12 2.2 Materials and methods ...... 16 2.3 Results ...... 22 2.4 Discussion ...... 31 Chapter 3: Phylogeography of the genus Antechinus ...... 39 3.1 Introduction ...... 40 3.2 Methods ...... 43 3.3 Results ...... 45 3.4 Discussion ...... 58 Chapter 4: Comparative population genetics of two small carnivorous marsupials in the genus Antechinus ...... 66 4.1 Introduction ...... 67 4.2 Methods ...... 71 4.3 Results ...... 77 4.4 Discussion ...... 88 Chapter 5: Life history and ecology of a new species of carnivorous marsupial, the buff-footed antechinus (A. mysticus) and a sympatric congener 97 5.1 Introduction ...... 98
vii
5.2 Methods ...... 102 5.1 Results ...... 106 5.2 Discussion ...... 120 Chapter 6: General discussion ...... 129 References ...... 143
viii
List of Figures
Figure . . BI phylogeny of the concatenated dataset with ML bootstrap values below and BI posterior probabilities above the line. Clades not reconstructed in the ML phylogeny are denoted with ‘*’. ML bootstrap support of and BI posterior probability values of . are not displayed. L - denote the main lineages of Antechinus...... 24 Figure . . BI phylogeny of the combined nuclear dataset with BI posterior probabilities shown above the node and ML bootstrap values shown below the node. Clades not reconstructed in the ML phylogeny are denoted with ‘*’. ML bootstrap support of and BI posterior probability values of . are not displayed...... 25 Figure . . Geographic distribution of Antechinus. (a) shows distribution of all Antechinus execpt the dusky antechinus and A. minimus lineage, (b) shows the distribution of the dusky antechinus and A. minimus lineage. A number of previously identified biogeographic barriers (see Bryant and Krosch ) are show on on the map...... 26 Figure . . BEAST phylogeny of the concatenated dataset. AS, cladogenesis dates and HPD values for lettered nodes are shown in Table . ...... 29 Figure . . BI phylogeny of Cytb dataset with BI posterior probabilities shown above the line and ML bootstrap values shown below the line. Clades not reconstructed in the ML phylogeny are denoted with ‘*’. Bootstrap support of , posterior probability values of . and all node support values for unmarked clades are not shown. For ease of reading node support values for southern A. f. flavipes subclades are also not shown. Tip labels and all support values are given in Fig. . - . L - denote the main Antechinus lineages identified in Chapter Two. See Supplementary Table . for more information on all haplotypes...... 50 Figure . . BI Cytb phylogeny of lineage extracted from Fig. . and map of corresponding sample locations. BI posterior probabilities are shown above the line and ML bootstrap values are shown below the line. Clades not reconstructed in the ML phylogeny are denoted with ‘*’. Colouring represents lineage membership. CB haplotype are displayed on the map with the letters ‘CB’ excluded. Shapes represent: triangles, A. swainsonii and A. vandycki; circles, A. minimus; rectangles, A. mimetes. The single haplotype of A. arktos is not shown on the map, see Supplementary Table . for more information on all haplotypes.
ix
Dotted lines delineate biogeographic barriers displayed in Fig. . ...... 50 Figure . . BI Cytb phylogeny of lineage extracted from Fig. . and map of corresponding sample locations. BI posterior probabilities are shown above the line and ML bootstrap values are shown below the line. Clades not reconstructed in the ML phylogeny are denoted with ‘*’. Colouring represents lineage membership. CB haplotype are displayed on the map with the letters ‘CB’ excluded. Shapes represent: triangles, A. agilis; circles, A. stuartii and A. subtropicus; red stars, sympatric sites of A. stuartii north and south. Dotted lines delineate biogeographic barriers displayed in Fig. . . NA and NA refer to A. subtropicus sites identified in Chapter Two which were not sequenced at Cytb. See Supplementary Table . for more information on all haplotypes...... 53 Figure . . BI Cytb phylogeny of lineage extracted from Fig. . and map of corresponding sample locations. BI posterior probabilities are shown above the line and ML bootstrap values are shown below the line. Clades not reconstructed in the ML phylogeny are denoted with ‘*’. Colouring represents lineage membership. CB haplotype are displayed on the map with the letters ‘CB’ excluded. Shapes represent: triangles, A. adustus and A. leo; circles, A. argentus and A. mysticus; rectangles, A. flavipes. Dotted lines delineate biogeographic barriers displayed in Fig. . . The Northern Territory A. bellus and Western Australian A. f. leucogaster are not shown on the map. See Supplementary Table . for more information on all haplotypes...... 55 Figure . . % connection limit parsimony network of A. stuartii south and A. stuartii north samples from the Cytb dataset. Red = New England NP, green = Werrikimbe NP, purple = Gosford, orange = Main Ranges NP, grey = Border Ranges NP. One haplotype (CB ) occurred at two sites, this is shown by the larger circle size of this haplotype. Only samples of a haplotype which occur > km apart are considered to be at separate sites. Blue does not represent a geographic location, see Supplementary Table. . for location information of all samples. Unfilled circles represent unsampled hypothetical haplotypes. Dotted lines around haplotypes indicate subclades shown Fig. . ...... 57 Figure . . Antechinus mysticus and Antechinus subtropicus populations sampled in the present study. Squares represent A. mysticus sites, circles represent A. subtropicus sites. The star represents a site at which both species were caught...... 72 Figure . . Population differentiation among the six A. mysticus populations sampled, visualised using Principal Coordinate Analysis...... 84
x
Figure 4.3. Graphical representation of membership coefficients of the Bayesian STRUCTURE analysis of microsatellite loci for A. mysticus obtained from sites across the known range of the species. Each plot represents different population assignments for K: (a) K = ; and (b) K = . Solid black lines delineate the different sites; each vertical line represents a single individual. Colours represent cluster assignments. A graph of the relationship of ∆K to K is shown in panel (c)...... 86 Figure 4.4. Population differentiation among the four A. subtropicus populations sampled, visualised using Principal Coordinate Analysis...... 87 Figure . . Graphical representation of membership coefficients of the Bayesian STRUCTURE analysis of microsatellite loci for A. subtropicus obtained from sites in SE Qld. Each plot represents different population assignments for K: (a) K = ; (b) K = . Solid black lines delineate the different sites; each vertical line represents a single individual. Colours represent cluster assignments. A graph of the relationship of ∆K to K is shown in panel (c)...... 87 Figure . . Trap success over a two year trapping period for (a) A. mysticus at AW, (b) A. mysticus at DW and (c) A. subtropicus at DW...... 111 Figure . . Minimum number of (a) A. mysticus and (b) A. subtropicus known to be alive (KTBA) at AW and DW. Antechinus subtropicus was only caught at DW...... 114 Figure 5.3. Change in monthly mean weight over the two year trapping period for (a) A. mysticus at AW, (b) A. mysticus at DW and (c) A. subtropicus at DW. Bars show standard error, values without standard error bars represents months in which only one individual was caught. Breaks in the connector lines indicate the change from the first to the second year of the study...... 116 Supplementary material Supplementary Figure . . ML phylogeny of the concatenated dataset with bootstrap values shown below the node. Bootstrap support is not shown for nodes with % support …………………….……… Supplementary Figure . . ML phylogeny of the nuclear dataset with bootstrap values shown below the node. Bootstrap support is not shown for nodes with % support…………………………….….…… Supplementary Figure . . ML phylogeny of the Cytb dataset with bootstrap values shown below the node. Bootstrap support is not shown for nodes with % support. Clades marked with blue lines represent recognised taxonomic units, purple lines mark putative taxonomic units clades revealed in the BI analyse shown in Fig. . ………………………………………………………………………….…….
xi
Supplementary Figure . . Graphical representation of membership coefficients of the Bayesian STRUCTURE analysis of microsatellite loci (Aa K is excluded) for A. mysticus obtained from sites across the known range of the species. Each plot represents different population assignments for K: (a) K = ; and (b) K = . Solid black lines delineate the different sites; each vertical line represents a single individual. Colours represent cluster assignments. Pairwise RST estimates of A. mysticus populations for these microsatellite loci are shown (c)………. Supplementary Figure . . Total rainfall in each month from - and mean rainfall per month since . Rainfall was recorded by the Australian Bureau of Meteorology at Samsonvale station ~ - km from the field sites ……………………………………………..……………
xii
List of Tables
Table . . Loci analysed for the concatenated and nuclear analyses. Loci sequenced are denoted with an ‘X’, detailed information on the samples is shown in Supplementary Table . . All samples denoted were analysed for the concatenated analyses. Only nuclear (IRBP, BRCA, vWF APOB) and mtDNA (Cytb, S) sequences were used in the respective combined analyses. Sequences for the sample marked with ‘*’ were taken from Genbank. Cytb and IRBP sequences from samples have previously been published in Baker et al. ( ; ; ; ), see Supplementary Table . ...... 17 Table . . List of primers used in this study. All primers were designed from primers used in Krajewski et al. ( ) (Cytb, S, IRBP) and Meredith et al. ( ) (APOB, BRCA , IRBP, vWF)...... 18 Table . . Nucleotide substitution models for BI (MrBayes and BEAST) analysis determined using JMODELTEST version . . ...... 20 Table 2.5. Results of BEAST analysis and AS reconstruction. Node lettering is shown on Fig. . . For the MP analysis = mesic- closed habitat; small/medium size; high altitude preferencing; = mesic-open habitat; large size; low altitude preferencing...... 30 Table . . Mean net proportion difference (%) between Antechinus species, subspecies and A. stuartii north and south...... 55 Table . . Geographic location of all sites. All sites are named after the National Park in which they occur, except Imbil, which is a State Forest and Crohamhurst, a Conservation Area. Abbreviations: ASL = above sea level; Am = A. mysticus; As = A. subtropicus...... 73 Table . . Primer sequence and references of the nine microsatellites loci genotyped in this study. The primers were divided into two multiplex groupings...... 75 Table . . Table of descriptive statistics and Hardy-Weinberg Equilibrium test for each microsatellite locus at each site for (a) A. mysticus and (b) A. subtropicus. N, number of samples; Na, number of alleles; Ho, observed heterozygosity; He, expected heterozygosity; HWE sig., Hardy-Weinberg Equilibrium significant at . level...... 79 Table . . Summary of genetic variation in sampled populations of (a) A. mysticus and (b) A. subtropicus based on eight and nine amplified microsatellite loci, respectively. N, sample size; A, total number of alleles; AR, allelic richness standardised for allele size;
xiii
uA, unique (private) alleles; rA, rare alleles (frequency ≤ %); He, expected heterozygosity; Ho, observed heterozygosity and FIS, inbreeding coefficient. FIS values significantly different (p = <