Phylogeography of the Koala, (Phascolarctos cinereus), and Harmonising Data to Inform Conservation

The Australian continent exhibits complex biogeographic patterns but studies of the impacts of Pleistocene climatic oscillation on the mesic environments of the Southern Hemisphere are limited. The koala (Phascolarctos cinereus), one of Australia’s most iconic species, was historically widely distributed throughout much of eastern Australia but currently represents a complex conservation challenge. To better understand the challenges to koala genetic health, we assessed the phylogeographic history of the koala. Variation in the maternally inherited mitochondrial DNA (mtDNA) Control Region (CR) was examined in 662 koalas sampled throughout their distribution. In addition, koala CR haplotypes accessioned to Genbank were evaluated and consolidated. A total of 53 unique CR haplotypes have been isolated from koalas to date (including 15 haplotypes novel to this study). The relationships among koala CR haplotypes were indicative of a single Evolutionary Significant Unit and do not support the recognition of subspecies, but were separated into four weakly differentiated lineages which correspond to three geographic clusters: a central lineage, a southern lineage and two northern lineages co-occurring north of Brisbane. The three geographic clusters were separated by known Pleistocene biogeographic barriers: the Brisbane River Valley and Clarence River Valley, although there was evidence of mixing amongst clusters. While there is evidence for historical connectivity, current koala populations exhibit greater structure, suggesting habitat fragmentation may have restricted female-mediated gene flow. Since mtDNA data informs conservation planning, we provide a summary of existing CR haplotypes, standardise nomenclature and make recommendations for future studies to harmonise existing datasets. This holistic approach is critical to ensuring management is effective and small scale local population studies can be integrated into a wider species context.


Introduction
The Australian continent exhibits complex biogeographic patterns as a consequence of the interplay between continental drift, topography and climatic change [1][2][3]. In contrast to the Northern Hemisphere, where continental ice sheets restricted species to relatively few major refugia during glacial cycles [4], the Southern Hemisphere is characterised by biogeographic barriers which periodically hampered or prevented gene flow/movement during periods of increased aridity. While multiple biogeographic barriers have been identified across Australia, these often appear idiosyncratic, with impacts on the phylogeographic patterns varying according to species-specific characters such as dispersal ability and habitat requirements [1,[5][6][7]. As a result, biogeographic barriers are frequently detected in some species but not others, and the limited number of studies restricts our understanding of the impacts of past climate change on the Australian flora and fauna.
The koala (Phascolarctos cinereus) is currently widely distributed throughout eastern Australia, from north Queensland to South Australia, including several insular populations (Fig 1; [8]). The mainland range encompasses several putative biogeographic barriers, including the St Lawrence Gap (SLG), the Brisbane Valley Barrier (BVB), the Clarence River Barrier (CRB) and the Hunter Valley Barrier (HVB). Koalas also have a complex management history, with numerous translocations (sometimes not comprehensively documented) having occurred throughout the southern parts of their range since the 1920s [8]. In Victoria (Vic), small numbers of koalas were introduced onto offshore islands, notably French and Phillip Islands. These populations were subsequently used as a source for reintroductions on mainland Vic, and potentially also in western New South Wales (NSW; [9]). In South Australia (SA), koalas were extinct by the 1930s [10,11] but populations in the Mt Lofty Ranges were established using koalas from Queensland (Qld), NSW and Vic, and this population was subsequently used to establish other populations on mainland SA [10,12,13]. Additionally, koalas were introduced to Kangaroo Island, SA from French Island, and from there subsequently used to establish populations on Eyre Peninsula, SA [10,12,13].
Koalas exhibit substantial morphological variation throughout their range, which appears to reflect clinal variation in response to environmental conditions (e.g. size and colouration; [12,[14][15][16]). Initially, three koala subspecies were recognised on the basis of variation in size and colouration: P. c. adustus from northern Qld, P. c. cinereus from NSW, and P. c. victor from Vic. However, the geographical boundaries of these subspecies were poorly defined and therefore became associated with state political boundaries [17]. Although populations at the extremes of the distribution appear morphologically and genetically differentiated, this variation appears clinal and is insufficient to support subspecies classification [12,14,16,18].
As an internationally recognised wildlife icon, the koala presents complex and challenging management issues. The conundrum arises due to the differing impacts of threatening processes and previous management across their range [19]. Throughout most of their range koala populations are subject to several threatening processes including historical and ongoing habitat loss and fragmentation (due to agriculture and urbanisation) [8], climate change and drought [20][21][22], infectious diseases, notably chlamydia and Koala Retrovirus [23][24][25][26], and the effects of increased conflicts with humans [12,27,28]. At European settlement, koalas were not abundant but increased in numbers with reduced indigenous hunting [8]. Subsequent hunting of koalas in the 1900s to support international fur trade and extensive habitat loss resulted in declines across the range and the extirpation of many populations in the southern parts of the range by the 1930s [8,11]. As a result of substantial regional declines in Qld, NSW and the Australian Capital Territory koalas are listed as 'vulnerable' under the Environmental Protection and Biodiversity Conservation Act 1999 [29]. In contrast, koalas occur at high densities in much of the southern parts of their range, leading to over-browsing [30] requiring population control measures to prevent starvation [8,[31][32][33]. Koala management and conservation is often further complicated by jurisdictional issues between Australian State and Commonwealth laws. Thus, in order to effectively manage koalas at a national level, an holistic approach which incorporates information from a range of geographic scales and considers each area in a wider distributional context is required.  A series of studies have previously assessed the population genetics of koalas at a range of scales [18,[34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49]. Together these studies have revealed information on dispersal, phylogeography, population history and genetic structure, but many have focused at local or regional scales, and few have attempted to place local populations into the broader national context. While the initial ground-breaking distribution-wide phylogeographic study by Houlden et al. [18] provided a strong basis for future studies, a subsequent lack of standard nomenclature and cross referencing for mitochondrial DNA (mtDNA) Control Region (CR) haplotypes resulted in difficulty placing subsequent smaller-scale studies into an overall national framework: a context which could have been used to effectively inform koala conservation and management.
This study revisits the question of koala phylogeography, examining the largest dataset gathered to date, including samples from a range of time periods (1870s to 2015) to provide information on the biogeographic patterns in eastern Australia and provide a robust framework for koala management and conservation. Specifically we aim to use mtDNA CR sequence data to investigate relationships among koala populations nationally, the impact of known biogeographic barriers, the extent of female mediated gene flow and signatures of demographic changes across the range. In addition, we aim to consolidate existing CR data to establish a coherent dataset which will enhance future studies and management decisions by providing a framework to enable local information to be considered in the context of the entire koala distribution.

Materials and Methods Sampling
Tissue samples were obtained from 662 wild koalas from throughout the distribution (Fig 1,  Table 1). Koala samples were obtained from researchers, consultants, as part of veterinary care at the Port Macquarie Koala Hospital and Australia Zoo Wildlife Hospital, and from the Australian Museum Tissue Collection. Sample collection was performed in accordance with methods approved by the Australian Museum Animal Ethics Committee (Permit Numbers: 11-03, 15-05). Australian Museum registration numbers are provided in S1 Table. Tissues samples were stored in 70-100% ethanol, RNA later or frozen at -20°C until DNA extraction.

DNA extraction and mtDNA amplification
Genomic DNA was extracted using Qiagen DNeasy Blood and Tissue kit (Qiagen GmbH, Hilden, Germany) following standard protocols or according to a high salt method [50]. An~850 base pair (bp) fragment of Domain 1 of the mtDNA CR (from the tRNA proline to the end of the central conserved region) was amplified using the marsupial-specific primers L15999M and H16498M [51]. PCRs were carried out in 25 μl reactions using 100-500 ng of genomic DNA, 1 x Reaction Buffer (Bioline My taq Red Reagent Buffer; Bioline, Australia), 2 pmol primers and Bioline My Taq Red DNA polymerase (0.5 unit). Negative controls were included in each PCR. Thermocycling was performed on an Eppendorf Mastercycler EpS (Eppendorf, Hamburg, Germany) under the following conditions; initial denaturation (94°C for 2 min); 38-45 cycles of denaturation (94°C for 20 s); annealing (60°C for 40 s) and extension (72°C for 50 s) followed by a final extension (5 min at 72°C). PCR products were cleaned using ExoSap-IT © (USB Corporation, Cleveland, Ohio, USA). Sequencing was resolved on an AB 3730xl Sequencer at AGRF Sydney.
Additionally, sequence data for the mtDNA CR was obtained from five museum specimens (GenBank KJ530551 to KJ530556; [38,52]). These individuals dated from 1870 to 1938 and the dried museum skin specimens were provided by the Bohusläns Museum, Göteborg Museum, Museum of Comparative Zoology at Harvard University, Queensland Museum and Stockholm Museum. All work was carried out in a laboratory dedicated to ancient DNA experiments (in the Leibniz Institute for Zoo and Wildlife Research) to avoid the risk of contamination with modern DNA. The DNA extraction, library and bait preparation and hybridisation procedures, and bioinformatics analyses for these samples are described in Tsangaras et al. [38] and Tsangaras et al. [52]. Hybridization capture PCR product baits were generated from modern koala genomic DNA (SN265; KJ530552.1) with CR primers (PCI-CR-NF:5 0 -CATCAACACCCAAA GCTGAT-3 0 and PCI-CR-NR: 5 0 -TTCTAGGTACGTCCGCAATCT-3 0 ). Subsequent library and bait preparation, and hybridisation capture are described in [52], with sequencing on an Illumina MiSeq platform at University of Copenhagen National High-throughput DNA Sequencing Center. Bioinformatics procedures are described in [52] with consensus sequences for all samples generated using Geneious v 6.1.8 [53].

Mitochondrial DNA analysis
Sequences were checked and edited with reference to chromatograms using Sequencher v 5.3.2 (Gene Codes Corporation, Ann Arbor, MI, USA). Any individuals with ambiguous sequence data or unique singleton haplotypes (i.e. it was the only individual possessing that sequence) were re-amplified and sequenced for verification. Unique haplotype sequences were lodged with GenBank under accession numbers KX618862 -KX618876).
In addition to the unique sequences generated in this study, we obtained existing published koala haplotypes available on GenBank (accession numbers AJ005846-AJ005863; KJ530551-KJ530556; KC505325; GQ851933-GQ851940; AJ012057-AJ012064; KF745869-KF745875). Since a number of koala mtDNA sequences of varying sizes exist, we compared all the existing haplotypes, identifying potential duplicates (i.e. identical sequences based on accessioned data). All haplotypes were then assigned standardised names as described in Table 2. All existing and novel sequences were then aligned using the CLUSTAL X algorithm implemented in MEGA 6 [54]. The most appropriate model of evolution was determined using jMODELTEST version 1.1 [55,56] using the Akaike Information Criterion (AIC). Phylogenetic analyses. Phylogenetic relationships amongst all 53 unique haplotypes (i.e. novel haplotypes reported here and those from previous studies) were estimated using Maximum Likelihood (ML) and Bayesian Inference (BI). ML analyses were implemented in GARLI 2.1 Web service [58]. The best ML tree was estimated using an adaptive search method with 1000 replicates to determine the optimal topology with 0.95 probability. The gamma distribution and proportion of invariant sites were estimated in GARLI. Support for the branching topology was evaluated with 1000 bootstrap replicates. BI analyses were calculated in BEAST v1.8.3. [59,60] using a Metropolis-coupled Markov chain Monte Carlo sampling approach run for 10 7 iterations, sampling every 1000 th . The most appropriate model of evolution available was selected, the Hasegawa-Kishino-Yano model (HKY; [61]), based on the results of jMO-DELTEST. Three independent replicates were conducted and inspected for consistency to check for local optima in TRACER. In addition a range of priors (including the default settings) were run and the robustness of the data assessed in TRACER (S1 File). Mitochondrial DNA CR sequence from the common wombat (Vombatus ursinus; NC_003322.1), the closest relative of the koala [62], was used as an outgroup in all analyses. The maximum credibility lineage tree and posterior probabilities were calculated in TREE ANNOTATOR, with the first 1000 trees discarded as burn-ins. A measure of within-group differences (D A ) was calculated between all identified lineages within the tree using MEGA.
A haplotype network was created including all 53 unique koala mtDNA CR haplotypes using the TCS procedure [63], which uses the probability of parsimony calculated for pairwise comparisons to create the network [64], implemented in PopART (http://popart.otago.ac.nz). Information on the frequency of haplotypes was only obtained from the sampling in the present study. Indels were coded as a single mutation. Two networks were generated: one including all 53 sequences with missing data masked and another including only sequences with >850 bp. To combine all the available information we manually compared and combined the two networks. We took a conservative approach to the number of mutations in short sequences, assuming they were identical to longer sequences where data was missing.
Mitochondrial DNA diversity and partitioning. For sites sampled in this study where at least 7 individuals were sampled, the levels of mtDNA CR diversity and differentiation were assessed. Haplotypic diversity (h), nucleotide diversity (π) and pairwise differentiation (F ST ) amongst populations were estimated using ARLEQUIN v3.5.2.2 [65]. Deviations from neutrality were examined for each of these 20 sampling localities, and overall using Fu's F S [66] and Tajima's D-statistic [67].
Partitioning of mtDNA diversity was assessed using an Analysis of Molecular Variance (AMOVA), implemented in ARLEQUIN. Partitioning was examined (i) within vs between populations across the range to determine general patterns of differentiation, and (ii) between populations from previously described subspecies (corresponding to Qld, NSW and Vic/SA) to determine if subspecies correspond to genetic partitioning and (iii) the same as (ii) but with all known reintroduced populations removed from the analysis, as translocated populations may disrupt natural levels of differentiation. A Spatial Analysis of Molecular Variance (SAMOVA), implemented in SAMOVA v2 was used to define homogenous groups [68]. This analysis partitions the populations into a specified number of clusters, maximising the variation between clusters (F CT ). The analysis was run with and without location data. To determine the most appropriate grouping for our dataset we estimated F CT for up to 20 populations (i.e. the number of sampled locations). The influence of geographic distance on the extent of differentiation (F ST ) was also assessed using Isolation By Distance, Web Service (IBDWS) version 3.23 [69]. The natural logarithm of F ST /(1-F ST ) [70] and geographic distance was used due to the level of variance in the dataset and the significance of this relationship was estimated using a Mantel's test, with 10000 permutations. This was tested using sites containing >7 individuals (see Table 1) and repeated using all sampled sites containing more than one individual.
Population history and demography. Mismatch analysis was performed in ARLEQUIN to test for historical demographic expansion events in our sampled populations [65]. We tested for historical expansion events in the species overall (i.e. all samples as a single population), and for each of the geographically defined lineage identified by phylogenetic analyses. Since koalas in the northern and southern regions of the range have different population histories [9,19], we tested for demographic events between populations where koalas have been extensively translocated (Vic/SA) and primarily natural populations (NSW/Qld). Concordance between these models and the empirical data was assessed by means of a least-squares approach [71]. To further assess demographic changes in koala populations Bayesian skyline plots (BSP) were generated in BEAST for each of the mtDNA CR lineages and the overall dataset. The HKY model of substitution (the most appropriate available based on jMODELTEST) was selected and a substitution rate of 15% per million years was selected, which has previously been applied to several macropod marsupials [5,6]. There are currently no suitable calibration points/fossils to calibrate a molecular clock for the koala. The analysis was run for 10 7 generations, sampling every 10000 th generation. The results were visualised in TRACER and calibrated for generation time (approximately 7 years; [72]).

Results
A total of 36 unique mtDNA CR haplotypes were identified in the 662 contemporary koalas sampled in this study. When compared to haplotypes previously accessioned to GenBank, 15 haplotypes were found to be novel and 21 had previously been reported ( Table 2). These 36 haplotypes were defined by 46 variable sites, of which 26 were parsimony informative. Searches of existing mtDNA CR data for koalas revealed 48 haplotypes have been published in the literature and accessioned into GenBank ( Table 2). These sequences vary in size and sections of the CR but fall into approximately 600 bp fragments and 850 bp fragments, where the larger fragment entirely overlaps the smaller fragments. We found that only 36 of these 48 haplotypes represent unique sequences. In order to align existing datasets a summary of these haplotypes, including matching haplotypes and our proposed standardised haplotype nomenclature is shown in Table 2. The standardised haplotype names outlined in Table 2 are referred to throughout the rest of this paper.

Phylogenetic analyses
The topology of the phylogenetic trees generated using BI and ML were similar. The BI maximum credibility lineage tree is shown in Fig 2. Four well supported but shallowly diverged lineages were present, with each mostly confined to discrete geographic areas. Lineages 1 and 2 were found north of Brisbane, Lineage 3 between Brisbane (i.e. Redlands; site 5) and Iluka (site 9), and Lineage 4 haplotypes occurred south of Iluka, although there is evidence of admixture across these boundaries (Figs 3 and 4). There was no consistent evidence of further substructure in the trees. The haplotype network (Fig 4) illustrates how shallow the divergences are amongst the four lineages, which ranged between 1-1.6%.

Mitochondrial DNA diversity and partitioning
Overall haplotype diversity of the koalas sampled in this study was 0.84 (±0.08) but varied substantially from 0.12 to 0.75 across the sites sampled. Nucleotide diversity was 0.73% across all koalas sampled, and again ranged widely among the sampling sites from no variation (i.e. a single haplotype) to 0.60% within-sites ( Table 1). The results from some populations should be treated with caution due to the limited number of samples. A comparison between the 850 bp fragments of the five historic samples and contemporary samples revealed that only one of these haplotypes was not represented in the contemporary samples. Significant deviations from neutrality were detected in each of the lineages and overall for Fu's F S , but Tajima's D values were not significant. This suggests the deviation from neutrality is relatively weak since Fu's F is more sensitive than Tajima's D to deviations from neutrality [73].
The results of the AMOVA revealed that genetic diversity was partitioned between sampling sites (88.84% of the variation) rather than within sites (11.16%). There was evidence of significant partitioning of genetic diversity between the previously described subspecies (both including and excluding introduced populations) but only moderate amounts of genetic diversity were partitioned between the 'subspecies' (35-43%), relative to between sites within (50-54%). The SAMOVA analysis identified three as the optimal number of groups. Regardless of incorporating geographic data the SAMOVA separated populations between Maryborough (site 4) and Redlands (site 5), and between Iluka (site 9) and Pine Creek (site 10).
Pairwise comparisons of population differentiation ranged from 0 to 1, and significant differences occurred between most sites ( Table 3). The exceptions were either geographically proximate (e.g. Coomera (site 6) and Ballina (site 8)) or occurred in the southern parts of the range where translocations have occurred (e.g. French Island (site 15) and Eyre Peninsula (site 19)). There were also instances of geographically proximate sites being genetically divergent (e.g. Pine Creek (site 10) and Iluka (site 9)). There was no significant relationship between genetic and geographic distances when sites containing seven or more individuals were  included (p = 0.98). There was a significant positive relationship when all the sampling sites (with n > 1) were included (p = 0.01), but this explained only 13% of the variation.

Population history and demography
The results of the mismatch analysis were consistent with a model of spatial expansion in all of the a priori groups tested, although the grouping of NSW and Qld haplotypes and the southern lineage failed to converge after 2000 steps. Models of demographic expansion could not be rejected as the non-linear least squares algorithm failed to converge after 2000 steps for all a priori groupings. BSP showed each lineage appeared relatively stable over time, except for the southern lineage, which exhibited an increase in population size (Fig 5). When assessed overall, BSP suggests koala numbers began increasing approximately 20000 years ago, and stabilised at the current size approximately 3000 years ago.

Phylogeography and genetic structure
Our study represents the largest and most extensive sample set yet used to investigate the phylogeography of the koala throughout its range, and suggests the koala has a more complex biogeographic history than previously considered. The four shallow genetic lineages identified here do not correspond to previously described sub-species, and it is now apparent that morphological variation within the koala represents a cline reflecting gradients in environmental conditions across the distribution [12,16]. In a previous phylogeographic study, Houlden et al. [18] found a similar pattern of genetic divergence, noting the presence of 'three clusters' in their phylogenetic analyses, which correspond to three of the four identified in the present study. However, further interpretation was hindered by their limited sampling (particularly from northern regions) and the presence of haplotypes from multiple lineages at a single site. Here, the more extensive sampling (and more powerful phylogenetic analyses) revealed additional fine scale details and the ability to delineate these lineages. The extent of mixing suggests these biogeographic features are no longer significant barriers to movement, and other factors such as habitat availability are likely to have a greater impact on contemporary gene flow as evidenced by several recent population genetics studies (e.g. [37,41,44]). For example, Dennison et al. [49] demonstrated dispersal throughout the region where the central (3) and southern (4) lineages adjoin. Thus, these results strengthen the conclusions of Houlden et al. [18] that the koala represents a single species (or Evolutionary Significant Unit (ESU); sensu Moritz [74]) with no sub-specific separation. Table 3. Levels of differentiation (ΦST) between samples localities (n > 7) of koalas.   Two biogeographic barriers, the BVB and CRB, previously identified in other eastern Australian species, have impacted the movement of koalas. The BVB appears to separate the two northern lineages (Lineages 1 and 2) from the central lineage (Lineage 3). This biogeographic barrier has previously been reported as impacting the phylogeography of a range of species including mammals [75], amphibians [76], reptiles [77] and invertebrates [78]. Detailed analysis of koalas in this region by Lee et al. [45] however, indicated haplotypes belonging to the central lineage (Lineage 3) occur along the coast north of the BVB (e.g. Brisbane area). This may be the result of movement following a reduction in the effectiveness of Pleistocene barriers, or suggest a role for the D'Aguilar Range, which runs roughly north-south and separates areas containing central lineage haplotypes (e.g. Brisbane area) from those with northern haplotypes (e.g. Brookfield (site d)) [45]. Further south, the central lineage (Lineage 3) and the southern lineage (Lineage 4) are separated by the CRB. Like the BVB this barrier no longer presents a substantial obstruction to movement as there is evidence of mixing of the two lineages over a large geographic area both in mtDNA (this study) and evidence of contemporary gene flow [49]. The CRB is apparent in several marsupials including long-nosed potoroos (Potorous tridactylus [6]), brush-tailed rock wallaby (Petrogale pencillata [79]) and Hastings River mouse (Pseudomys oralis [80]) as well as several reptile species [81]. Finally, our phylogenetic analyses indicated the presence of two lineages in the northern parts of the range, but there was extensive mixing across a broad geographic area and it was not possible to elucidate any potential barriers using our data. To date the genetics of the koalas in this region have been poorly investigated and further sampling in this area will be required to identify if any biogeographic barriers previously identified in other species have impacted on the koala in this region (e.g. SLG; [75,82,83]).

Whitsunday/ Mackay
There is no evidence of barriers further south, such as the HVB impacting on phylogeographic patterns in koalas, reflecting the idiosyncratic nature of many Pleistocene barriers [6,79]. It is possible that the widespread translocation of koalas throughout most of the southern parts of the range (i.e. Vic) has masked any potential phylogeographic patterns in this region. The spread of haplotypes throughout Vic and SA reflects koalas' more recent history of documented translocations, and suggests some potentially undocumented translocations. For example, koalas at Eyre Peninsula (site 19), Kangaroo Island (site 20) and Bessiebelle (site 17) are all similar to French Island (site 15) which, along with Phillip Island, was the reported source population [9]. In contrast, the Mt. Lofty koalas exhibit high levels of diversity, among the highest of the southern populations, consistent with their founding individuals being obtained from a wider area including NSW and Qld [9]. It has been proposed that the genetic diversity of koalas in the Strzelecki and south Gippsland areas may not have been as heavily impacted by translocations as other areas [18,39,46]. However, the levels of diversity present are still comparatively low. Interestingly, koalas from Cape Otway possessed a single haplotype not found in other Victorian sites, but that we documented in northern and western NSW locations. This is inconsistent with documented translocations from French and Phillip Islands to this general area [9]. It is likely this is indicative of undocumented translocations, but further work would be required to distinguish between remnant haplotypes and undocumented translocations.
Overall, the spread of haplotypes and lack of substructure within lineages suggests koalas have maintained effective (female) genetic connectivity historically (i.e. over evolutionary time). However, significant differentiation and partitioning of mtDNA CR diversity among sites is indicative of limited female movement on shorter timescales (i.e. ecological time). These results are consistent with previous regional and local population genetic studies of koalas and suggests that differentiation of koala populations is associated with contemporary (i.e. Post-European) barriers to dispersal [18,44,45,49]. Furthermore, ecological studies have suggested that koalas are capable of moving large distances [84,85]. In contrast to Qld and NSW, populations in Vic and SA exhibit low levels of divergence, which is consistent with previous genetic studies and the management history of koalas in these regions [9,18,39]. Together this information suggests that recent human-induced landscape changes and habitat loss coupled with founder and bottleneck effects are likely to have contributed to this pattern.

Historic genetic diversity and demography
The levels of mtDNA diversity detected in koalas in this study are typical of arboreal species and those with specific habitat requirements [6]. The combination of habitat loss, disease and intesense hunting pressure to support an international fur trade resulted in dramatic declines in koala populations and the exitrpation of several southern populations by the 1930s [8,9,11]. The impact of such declines would be expected to be apparent in the genetics of koalas [47,66] and several studies have found evidence of genetic bottlenecks in some populations [39,47]. Sites in SA and Vic exhibited lower levels of diversity compared with those in Qld and NSW, which may be associated with reintroduction from a very limited source following extirpation [9,18]. However, we found no evidence of recent genetic bottlenecks or declines using mtDNA CR data; rather our BSP indicated koala populations have remained relatively stable or increased during the past 1000 years. Tsangaras et al. [38] found little evidence of the loss of diversity in mtDNA CR haplotypes when comparing contemporary samples to those from the late 1800s and early 1900s. Furthermore, our extensive sampling of contemporary koalas revealed that only one of these haplotypes appears restricted to historic specimens.
The BSP does however suggest a post-LGM expansion, particularly into the southern regions of the range (Fig 5). This pattern of expansion is consistent with bioclimatic models, which suggest koalas were largely restricted to several fragmented areas in northeastern NSW and southeastern Qld during the LGM [86]. It is also possible that this earlier expansion event has led to low levels of diversity in this region which hamper the detection of this most recent bottleneck (particularly when combined with the lower effective size of mtDNA) [42,87]. Much more work, particularly studies including samples from Pleistocene deposits will be required to elucidate patterns of mtDNA CR diversity loss in koalas.

Aligning mtDNA datasets
Sequence variation in mtDNA CR has been widely employed to assess genetic diversity and phylogeography, and still forms the basis of management programmes for many species [88]. In koalas, mtDNA has been utilized at a range of scales, including local and distributional scales, using both contemporary and historic samples [18,38,41,45]. Sequence data is easily transferable between labs, and readily shared via online repositories (i.e. GenBank). However, to date different studies have employed different naming systems and there is some duplication of sequences (but with different names assigned by different authors), which can add to confusion about levels of diversity and structure. Table 2 matches new and existing haplotypes and integrates them into a single nomenclature. This will improve the integration of datasets, but it is essential that researchers adopt common nomenclature to ensure consistency and comparability in the future. We acknowledge that this is not straightforward where different length sequences are used. In the long term, the decreasing costs of Next Generation Sequencing (NGS) and use of whole mitogenome sequencing may reduce this problem but currently, targeted sequencing of mtDNA CR still forms the basis of many management strategies. We suggest researchers undertake simple measures such as ensuring only new unique sequences are named. For complex systems, it may be necessary for researchers to adopt strict nomenclature systems, such as those used in human genetics. One option is adopting the revised Cambridge Naming System (rCNS), where haplotype names are based on the SNPs which distinguish that haplotype from a reference sequence [89,90]. This minimises the length of names but also provides some information on the relationships among haplotypes. Adopting such conventions would require a shift in the approach of wildlife researchers to a more complex system, but one which may streamline phylogeographic studies, particularly where multiple local studies exist.

Management and Conservation
Koala conservation is not straightforward, due to the varying impacts of threatening processes and jurisdictional issues [19]. Effective management requires information at a range of scales but ensuring local or regional level studies can be placed into the broader distributional context. To date, this has not been achieved for koalas. To address this issue we have: (1) Synthesised existing data for mtDNA CR, providing a standardised nomenclature and framework for management and future work ( Table 2); and (2) Used the largest dataset to date to assess the phylogeography of koalas that can be used to inform management. Our results confirm that koalas should be managed as a single ESU, with no sub-specific separation. However, the presence of low levels of genetic divergence resulting from Pleistocene barriers and morphological variation, which most likely reflects environmental clines [14,16], should be taken into consideration in maximising the conservation of genetic diversity. These results also highlight the importance of nationally significant koala populations in southeastern Qld and northeastern NSW, as this area contains a high proportion of the mtDNA diversity present in koalas, but which is under threat from urbanisation. This is particularly important, given these results support previous inferences that movement and dispersal in koalas has been restricted by habitat fragmentation. This nationwide phylogeographic study provides a strong framework for fitting information from local studies into a national context.