Molecular ecology meets remote sensing: environmental drivers to population structure of humpback dolphins in the Western Indian Ocean

Genetic analyses of population structure can be placed in explicit environmental contexts if appropriate environmental data are available. Here, we use high-coverage and high-resolution oceanographic and genetic sequence data to assess population structure patterns and their potential environmental influences for humpback dolphins in the Western Indian Ocean. We analyzed mitochondrial DNA data from 94 dolphins from the coasts of South Africa, Mozambique, Tanzania and Oman, employing frequency-based and maximum-likelihood algorithms to assess population structure and migration patterns. The genetic data were combined with 13 years of remote sensing oceanographic data of variables known to influence cetacean dispersal and population structure. Our analyses show strong and highly significant genetic structure between all putative populations, except for those in South Africa and Mozambique. Interestingly, the oceanographic data display marked environmental heterogeneity between all sampling areas and a degree of overlap between South Africa and Mozambique. Our combined analyses therefore suggest the occurrence of genetically isolated populations of humpback dolphins in areas that are environmentally distinct. This study highlights the utility of molecular tools in combination with high-resolution and high-coverage environmental data to address questions not only pertaining to genetic population structure, but also to relevant ecological processes in marine species.     

Population Genetic Studies Revealed Local Adaptation in a High Gene-Flow Marine Fish,the Small Yellow Croaker (Larimichthys polyactis)

The genetic differentiation of many marine fish species is low. Yet local adaptation may be common in marine fish species as the vast and changing marine environment provides more chances for natural selection. Here, we used anonymous as well as known protein gene linked microsatellites and mitochondrial DNA to detect the population structure of the small yellow croaker (Larimichthys polyactis) in the Northwest Pacific marginal seas. Among these loci, we detected at least two microsatellites, anonymous H16 and HSP27 to be clearly under diversifying selection in outlier tests. Sequence cloning and analysis revealed that H16 was located in the intron of BAHCC1 gene. Landscape genetic analysis showed that H16 mutations were significantly associated with temperature, which further supported the diversifying selection at this locus. These marker types presented different patterns of population structure: (i) mitochondrial DNA phylogeny showed no evidence of genetic divergence and demonstrated only one glacial linage; (ii) population differentiation using putatively neutral microsatellites presented a pattern of high gene flow in the L. polyactis. In addition, several genetic barriers were identified; (iii) the population differentiation pattern revealed by loci under diversifying selection was rather different from that revealed by putatively neutral loci. The results above suggest local adaptation in the small yellow croaker. In summary, population genetic studies based on different marker types disentangle the effects of demographic history, migration, genetic drift and local adaptation on population structure and also provide valuable new insights for the design of management strategies in L. polyactis.     

Isolation by environmental distance in mobile marine species: molecular ecology of franciscana dolphins at their southern range

The assessment of population structure is a valuable tool for studying the ecology of endangered species and drafting conservation strategies. As we enhance our understanding about the structuring of natural populations, it becomes important that we also understand the processes behind these patterns. However, there are few rigorous assessments of the influence of environmental factors on genetic patterns in mobile marine species. Given their dispersal capabilities and localized habitat preferences, coastal cetaceans are adequate study species for evaluating environmental effects on marine population structure. The franciscana dolphin, a rare coastal cetacean endemic to the Western South Atlantic, was studied to examine these issues. We analysed genetic data from the mitochondrial DNA and 12 microsatellite markers for 275 franciscana samples utilizing frequency‐based, maximum‐likelihood and Bayesian algorithms to assess population structure and migration patterns. This information was combined with 10 years of remote sensing environmental data (chlorophyll concentration, water turbidity and surface temperature). Our analyses show the occurrence of genetically isolated populations within Argentina, in areas that are environmentally distinct. Combined evidence of genetic and environmental structure suggests that isolation by distance and a process here termed isolation by environmental distance can explain the observed correlations. Our approach elucidated important ecological and conservation aspects of franciscana dolphins, and has the potential to increase our understanding of ecological processes influencing genetic patterns in other marine species.     

Population genetic analysis identifies source-sink dynamics for two sympatric garter snake species (Thamnophis elegans and Thamnophis sirtalis)

Population genetic structure can be shaped by multiple ecological and evolutionary factors, but the genetic consequences of these factors for multiple species inhabiting the same environment remain unexplored. We used microsatellite markers to examine the population structures of two coexisting species of garter snake, Thamnophis elegans and Thamnophis sirtalis, to determine if shared landscape and biology imposed similar population genetic structures. These snakes inhabit a series of ponds, lakes and flooded meadows in northern California and tend to converge on prey type wherever they coexist. Both garter snakes had comparable effective population sizes and bidirectional migration rates (estimated using a maximum-likelihood method based on the coalescent) with low but significant levels of genetic differentiation (F(ST) = 0.024 for T. elegans and 0.035 for T. sirtalis). Asymmetrical gene flow revealed large source populations for both species as well as potential sinks, suggesting frequent extinction-recolonization and metapopulation dynamics. In addition, we found a significant correlation between their genetic structures based on both pairwise F(ST)s for shared populations (P = 0.009) and for bidirectional migration rates (P = 0.024). Possible ecological and evolutionary factors influencing similarities and differences in genetic structure for the two species are discussed. Genetic measures of effective population size and migration rates obtained in this study are also compared with estimates obtained from mark-recapture data.     

Conservation genetics of the franciscana dolphin in Northern Argentina: population structure,by-catch impacts,and management implications

Evaluating population structure in the marine environment is a challenging task when the species of interest is continuously distributed, and yet the use of population or stock structure is a crucial component of management and conservation strategies. The franciscana dolphin (Pontoporia blainvillei), a rare endangered coastal cetacean, suffers high levels of by-catch all along its distribution range in the Western South Atlantic, and questions have been raised about boundaries or divisions for population management. Here we apply genetic tools to better understand population structure and migration, sex-biased dispersal, and to assess potential genetic and demographic impacts of by-catch. Our analyses, based on mtDNA control region sequences, reveal significant genetic division at the regional level and fine-scale structure within our study area. These results suggest that the population in northern Buenos Aires is the most isolated population in Argentina. We found no significant departure from an equal sex ratio among the by-caught animals. A few cases of multiple entanglements appeared to be mother–calf pairs based on field observations and individuals sharing the same mtDNA control region lineage. The distribution of haplotype frequencies observed could imply that some maternal lineages are more prone to be subject to higher rates of by-catch, although biopsy sampling is necessary to fully evaluate whether maternal lineage distributions are the same for biopsy sampled and by-caught animals. A genetic indication of population size disequilibrium was detected for all populations in Argentina, which is consistent with available rates of by-catch and abundance estimates. Collectively, our findings support the current scheme of larger recognized Franciscana Management Areas (FMA), but argue for a finer-scale subdivision within Northern Buenos Aires region (FMA IV). Finally, an integrated approach to promote conservation of this endangered small cetacean has to involve identification of genetic and demographic threats, a more sustainable fishery strategy to reduce by-catch, and designation of protected areas that are supported by underlying population structure for franciscana dolphins.     

Detecting populations in the 'ambiguous' zone: kinship-based estimation of population structure at low genetic divergence

Identifying population structure is one of the most common and important objectives of spatial analyses using population genetic data. Population structure is detected either by rejecting the null hypothesis of a homogenous distribution of genetic variation, or by estimating low migration rates. Issues arise with most current population genetic inference methods when the genetic divergence is low among putative populations. Low levels of genetic divergence may be as a result of either high ongoing migration or historic high migration but no current, ongoing migration. We direct attention to recent developments in the use of the tempo-spatial distribution of closely related individuals to detect population structure or estimate current migration rates. These 'kinship-based' approaches complement more traditional population-based genetic inference methods by providing a means to detect population structure and estimate current migration rates when genetic divergence is low. However, for kinship-based methods to become widely adopted, formal estimation procedures applicable to a range of species life histories are needed.     

Biogeography of the smooth snake (Coronella austriaca): origin and conservation of the northernmost population

Understanding historical range expansions and population demography can be crucial for the conservation and management of endangered species. In doing so, valuable information can be obtained regarding, for example, the identification of isolated populations, associations to particular habitats and distribution range shifts. As poikilotherms, snakes are vulnerable to environmental changes that can greatly shape their distribution ranges. Here we used mitochondrial data to elucidate the origin of the smooth snake population in Åland island, which is the northernmost location where the species is found. In Åland, we used mitochondrial and microsatellite data to fine‐map its spatial genetic structure, infer its demographic dynamics and determine its effective population size. We found three independent lineages, which expanded north from Iberian, the Balkans and Caucasus regions. The central lineage originating in the Balkans was the only one that reached Scandinavia. The Åland population belongs to this lineage and potentially colonized the island from the west via Sweden. This population appeared to be critically small and fragmented into two genetically isolated subpopulations. We discuss our results in light of previous findings regarding colonization routes in Europe and Scandinavia. Moreover, we discuss the origin and current genetic status of the Åland population relative to other co‐occurring snakes and suggest conservation measures based on our findings. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 426–435.     

Relative information content of polymorphic microsatellites and mitochondrial DNA for inferring dispersal and population genetic structure in the olive sea snake, Aipysurus laevis

Polymorphic microsatellites are widely considered more powerful for resolving population structure than mitochondrial DNA (mtDNA) markers, particularly for recently diverged lineages or geographically proximate populations. Weaker population subdivision for biparentally inherited nuclear markers than maternally inherited mtDNA may signal male-biased dispersal but can also be attributed to marker-specific evolutionary characteristics and sampling properties. We discriminated between these competing explanations with a population genetic study on olive sea snakes, Aipysurus laevis. A previous mtDNA study revealed strong regional population structure for A. laevis around northern Australia, where Pleistocene sea-level fluctuations have influenced the genetic signatures of shallow-water marine species. Divergences among phylogroups dated to the Late Pleistocene, suggesting recent range expansions by previously isolated matrilines. Fine-scale population structure within regions was, however, poorly resolved for mtDNA. In order to improve estimates of fine-scale genetic divergence and to compare population structure between nuclear and mtDNA, 354 olive sea snakes (previously sequenced for mtDNA) were genotyped for five microsatellite loci. F statistics and Bayesian multilocus genotype clustering analyses found similar regional population structure as mtDNA and, after standardizing microsatellite F statistics for high heterozygosities, regional divergence estimates were quantitatively congruent between marker classes. Over small spatial scales, however, microsatellites recovered almost no genetic structure and standardized F statistics were orders of magnitude smaller than for mtDNA. Three tests for male-biased dispersal were not significant, suggesting that recent demographic expansions to the typically large population sizes of A. laevis have prevented microsatellites from reaching mutation-drift equilibrium and local populations may still be diverging.     

Genetic structure reveals management units for the yellow cardinal (<Emphasis Type="Italic">Gubernatrix cristata</Emphasis>), endangered by habitat loss and illegal trapping

The yellow cardinal, Gubernatrix cristata, is an endangered passerine from southern South America. Populations are declining due to the loss of their natural habitat, which has caused a fragmented distribution, and the continuous extraction of individuals from the wild, mainly males, to sell them as cage birds. In this study, we assess the genetic variability of remaining yellow cardinal’s populations and determine whether these populations represent independent management units. We found that the degree of geographic isolation of the remaining populations parallels the genetic differentiation of these populations for both mitochondrial and nuclear markers, and supports the delimitation of four management units for the yellow cardinal (three in Argentina and one in Uruguay). Assignment tests showed that geographic genetic differentiation can be used to assign seized individuals from illegal pet trade to their original populations and thus manage their release.     

The high Andes, gene flow and a stable hybrid zone shape the genetic structure of a wide-ranging South American parrot

Background

While the gene flow in some organisms is strongly affected by physical barriers and geographical distance, other highly mobile species are able to overcome such constraints. In southern South America, the Andes (here up to 6,900 m) may constitute a formidable barrier to dispersal. In addition, this region was affected by cycles of intercalating arid/moist periods during the Upper/Late Pleistocene and Holocene. These factors may have been crucial in driving the phylogeographic structure of the vertebrate fauna of the region. Here we test these hypotheses in the burrowing parrot Cyanoliseus patagonus (Aves, Psittaciformes) across its wide distributional range in Chile and Argentina.

Results

Our data show a Chilean origin for this species, with a single migration event across the Andes during the Upper/Late Pleistocene, which gave rise to all extant Argentinean mitochondrial lineages. Analyses suggest a complex population structure for burrowing parrots in Argentina, which includes a hybrid zone that has remained stable for several thousand years. Within this zone, introgression by expanding haplotypes has resulted in the evolution of an intermediate phenotype. Multivariate regressions show that present day climatic variables have a strong influence on the distribution of genetic heterogeneity, accounting for almost half of the variation in the data.

Conclusions

Here we show how huge barriers like the Andes and the regional environmental conditions imposed constraints on the ability of a parrot species to colonise new habitats, affecting the way in which populations diverged and thus, genetic structure. When contact between divergent populations was re-established, a stable hybrid zone was formed, functioning as a channel for genetic exchange between populations.     

The phylogeography of dusky dolphins (Lagenorhynchus obscurus): a critical examination of network methods and rooting procedures

We investigated the phylogeography and evolutionary history of dusky dolphins (Lagenorhynchus obscurus) using DNA sequences of the full mitochondrial cytochrome b gene in 124 individuals from the putative stocks off Peru, Argentina and Southwest Africa. While genetic differentiation within oceans is surprisingly low, there is no evidence for recent female gene flow between Atlantic and Pacific waters. Highest genetic variability in terms of sequence divergence and number of haplotypes is found in the Atlantic. Our analyses also indicate that the eastern South Pacific dusky dolphins stock should be considered a separate management unit. Given the high level of mortality experienced by the Peruvian dusky dolphin in local fishery activities, these findings have important implications for an objective management of the species. Furthermore, we analysed our mitochondrial sequence data with several widely used network estimation and rooting methods. The resulting intraspecific gene genealogies and rooting inferences exhibited substantial differences, underlying the limitations of some algorithms. Given that scientific hypotheses and management decisions depend strongly on inferred tree or network topologies, there is a clear need for a systematic comparative analysis of available methods. Finally, the present study indicates that (i) the dusky and the Pacific white-sided dolphins are sister species and (ii) not only the Westwind Drift hypothesis but also other models of dispersion are compatible with the current geographical distribution of dusky dolphins.     

Ecological correlates of population genetic structure: a comparative approach using a vertebrate metacommunity

Identifying ecological factors associated with population genetic differentiation is important for understanding microevolutionary processes and guiding the management of threatened populations. We identified ecological correlates of several population genetic parameters for three interacting species (two garter snakes and an anuran) that occupy a common landscape. Using multiple regression analysis, we found that species interactions were more important in explaining variation in population genetic parameters than habitat and nearest-neighbour characteristics. Effective population size was best explained by census size, while migration was associated with differences in species abundance. In contrast, genetic distance was poorly explained by the ecological correlates that we tested, but geographical distance was prominent in models for all species. We found substantially different population dynamics for the prey species relative to the two predators, characterized by larger effective sizes, lower gene flow and a state of migration-drift equilibrium. We also identified an escarpment formed by a series of block faults that serves as a barrier to dispersal for the predators. Our results suggest that successful landscape-level management should incorporate genetic and ecological data for all relevant species, because even closely associated species can exhibit very different population genetic dynamics on the same landscape.     

Genetic diversity in the US ex situ populations of the endangered Puerto Rican Boa,Chilabothrus inornatus

The Puerto Rican Boa (Chilabothrus inornatus) was placed on the US Endangered Species List in 1970. Progress has been made since to clarify the recovery status of this species, though the design of a new recovery plan must include information regarding genetic variation within and among populations of this species. While measures of genetic diversity in wild populations of this species are finally becoming available, relative genetic diversity represented in ex situ populations is unknown, which hampers efforts to develop an ex situ species management plan. Here, we provide an analysis of genetic diversity in US public and private collections (zoos and breeders) using mitochondrial sequence data and five highly polymorphic nuclear microsatellite loci. We analyzed 50 boas from the US ex situ population and determined overall genetic diversity and relatedness among these individuals. We then compared these data to mitochondrial and microsatellite data obtained from 176 individuals from wild populations across the native range of the species. We found little inbreeding and a large amount of retained genetic diversity in the US ex situ population of C. inornatus relative to wild populations. Genetic diversity in the ex situ population is similar to that found in wild populations. Ours is only the second explicit attempt to characterize genetic diversity at the molecular level in ex situ populations of boid snakes. We anticipate that these results will inform current breeding strategies as well as offer additional information that will facilitate the continuation of ex situ conservation breeding or management in boas.     

Historical range contraction,and not taxonomy,explains the contemporary genetic structure of the Australian tree <Emphasis Type="Italic">Acacia dealbata</Emphasis> Link

Irrespective of its causes, strong population genetic structure indicates a lack of gene flow. Understanding the processes that underlie such structure, and the spatial patterns it causes, is valuable for conservation efforts such as restoration. On the other hand, when a species is invasive outside its native range, such information can aid management in the non-native range. Here we explored the genetic characteristics of the Australian tree Acacia dealbata in its native range. Two subspecies of A. dealbata have previously been described based on morphology and environmental requirements, but recent phylogeographic data raised questions regarding the validity of this taxonomic subdivision. The species has been widely planted within and outside its native Australian range and is also a highly successful invasive species in many parts of the world. We employed microsatellite markers to investigate the population genetic diversity and structure among 42 A. dealbata populations from across the species’ native range. We also tested whether environmental variables purportedly relevant for the putative separation of subspecies are linked with population genetic differentiation. We found no relationship between population genetic structure of A. dealbata in Australia and these environmental features. Rather, we identified two geographically distinct genetic clusters that corresponded with populations in the northeastern part of mainland Australia, and the southern mainland and Tasmanian range of the species. Our results do not support the taxonomic subdivision of the species into two distinct subspecies based on environmental features. We therefore assume that the observed morphological differences between the putative subspecies are plastic phenotypic responses. This study provides population genetic information that will be useful for the conservation of the species within Australia as well as to better understand the invasion dynamics of A. dealbata.     

Historical diversification of migration patterns in a passerine bird

Migratory strategies of birds require complex orientation mechanisms, morphological adaptations, and life-history adjustments. From an evolutionary perspective, it is important to know how fast this complex combination of traits can evolve. We analyzed mitochondrial control-region DNA sequences in 241 blackcaps (Sylvia atricapilla) from 12 populations with different migratory behaviors. The sample included sedentary populations in Europe and Atlantic archipelagos and migratory populations with different distances of migration, from regional to intercontinental migrations, and different heading directions (due to a migratory divide in central Europe). There was no genetic structure between migratory and sedentary populations, or among populations from different biogeographic areas (Atlantic islands, the Iberian Peninsula, or the continent), however we found evidence of a genetic structure when comparing populations located on either side of the migratory divide. These findings support an independent evolution of highly divergent migratory strategies in blackcaps, occurring after a postglacial colonization of the continent along western and eastern routes. Accordingly, mismatch-distribution analyses suggested an expansion of blackcaps from a very small population size, and time estimates dated such an expansion during the last postglacial period. However, the populations in Gibraltar, located in a putative Mediterranean refuge, appeared to be independent of these processes, showing evidence of restricted gene flow with other populations and demonstrating insignificant historical changes in effective population size. Our results show that the interruption of gene flow between migratory and sedentary populations is not necessary for the maintenance of such a polymorphism, and that even the most divergent migratory strategies of a bird species are susceptible to evolution in response to historical environmental changes.     

High interpopulation homogeneity in Central Argentina as assessed by Ancestry Informative Markers (AIMs)

The population of Argentina has already been studied with regard to several genetic markers, but much more data are needed for the appropriate definition of its genetic profile. This study aimed at investigating the admixture patterns and genetic structure in Central Argentina, using biparental markers and comparing the results with those previously obtained by us with mitochondrial DNA (mtDNA) in the same samples. A total of 521 healthy unrelated individuals living in 13 villages of the Córdoba and San Luis provinces were tested. The individuals were genotyped for ten autosomal ancestry informative markers (AIMs). Allele frequencies were compared with those of African, European and Native American populations, chosen to represent parental contributions. The AIM estimates indicated a greater influence of the Native American ancestry as compared to previous studies in the same or other Argentinean regions, but smaller than that observed with the mtDNA tests. These differences can be explained, respectively, by different genetic contributions between rural and urban areas, and asymmetric gene flow occurred in the past. But a most unexpected finding was the marked interpopulation genetic homogeneity found in villages located in diverse geographic environments across a wide territory, suggesting considerable gene flow.     

High gene flow across large geographic scales reduces extinction risk for a highly specialised coral feeding butterflyfish

The vulnerability of ecologically specialised species to environmental fluctuations has been well documented. However, population genetic structure can influence vulnerability to environmental change and recent studies have indicated that specialised species may have lower genetic diversity and greater population structuring compared to their generalist counterparts. To examine whether there were differences in population genetic structure between a dietary specialist (Chaetodon trifascialis) and a dietary generalist (Chaetodon lunulatus) we compared the demographic history and levels of gene flow of two related coral-feeding butterflyfishes. Using allele frequencies of ≥11 microsatellite loci and >350 bases of mitochondrial control region sequence our analyses of C. trifascialis and C. lunulatus from five locations across the Pacific Ocean revealed contrasting demographic histories and levels of genetic structure. Heterozygosity excess tests, neutrality tests and mismatch distributions were all highly significant in the dietary specialist C. trifascialis (all P < 0.01), suggesting genetic bottlenecks have occurred in all locations. In contrast, we found little evidence of genetic bottlenecks for the dietary generalist C. lunulatus. High gene flow and low genetic structuring was detected among locations for C. trifascialis (amova: R(ST) = 0.0027, P = 0.371; Φ(ST) = 0.068, P < 0.0001). Contrary to our expectations, a greater level of genetic structuring between locations was detected for C. lunulatus (amova: R(ST) = 0.0277, Φ(ST) = 0.166, both P < 0.0001). These results suggest that dietary specialisation may affect demographic history through reductions in population size following resource declines, without affecting population structure through reductions in gene flow in the same way that habitat specialisation appears to. Although C. trifascialis is highly vulnerable to coral loss, the high gene flow detected here suggests populations will be able to recover from local declines through the migration of individuals.     

Discordant patterns of population structure for two co‐distributed snake species across a fragmented Ontario landscape

Aim The goal of our study was to investigate the effects of a fragmented landscape on the genetic population structure of two sympatric snake species that differ in habitat preference. The eastern garter snake (Thamnophis sirtalis sirtalis) is a common, habitat generalist, whereas the endangered eastern foxsnake (Mintonius [Elaphe] gloydi) is rarer, geographically restricted, and a marsh‐specialist. We were most interested in comparing the genetic population structure of both species and identifying any natural and human‐created features of the landscape that overlap with genetic disjunctions. Location Southwestern Ontario, Canada, surveying over half of the remaining range of the eastern foxsnake. Methods We utilized DNA microsatellite markers to examine genetic population structure of both species. The number of genetically distinct clusters for each species was determined using both Bayesian spatial assignment and spatial principal component analyses (sPCA). Genetic clusters were overlaid onto a habitat map to deduce possible physiognomic barriers to gene flow. Results Spatial assignment revealed three genetic clusters for garter snakes and five for foxsnakes. Each individual garter snake had a near equal probability of membership to two or more clusters with no cluster mapping onto a discrete geographic region, indicating that garter snakes comprise a single genetic population. The identified foxsnake clusters correspond to geographically circumscribed locations on the landscape, roughly coincident with isolated patches of suitable habitat. sPCAs revealed significant global allelic structure for foxsnakes, but not for garter snakes. No significant local structure was found for either species. Main Conclusions Our results imply that foxsnakes and garter snakes are differentially impacted by the same landscape or have dramatically different effective population sizes. Unsuitable intervening habitat such as agricultural tracts and roads between existing populations of foxsnakes appears to act as barriers to gene flow, while garter snake movement appears unrestricted by these features. Our findings have important implications for the management of eastern foxsnakes.     

Genetic structure and diversity of the koala population in South Gippsland,Victoria: a remnant population of high conservation significance

In the Australian state of Victoria, the history of koalas and their management has resulted in the homogenisation and reduction of genetic diversity in many contemporary populations. Decreased genetic diversity may reduce a species’ ability to adapt to future environmental pressures such as climate change or disease. The South Gippsland koala population is considered to be unique in Victoria, as it is believed to be a remnant population, not originating from managed populations that have low genetic variation. This study investigated genetic structure and diversity of koalas in South Gippsland, with comparison to other populations in Victoria (French Island/Cape Otway, FI and Raymond Island, RI), New South Wales and south east Queensland. Population analyses were undertaken using both microsatellite genotype and mitochondrial DNA sequence data. Non-invasive sampling of koala scats was used to source koala DNA, allowing 222 South Gippsland koalas to be genotyped. Using nuclear data the South Gippsland koala population was found to be significantly differentiated (D_jost 95% CI SG–RI?=?0.03–0.06 and SG–FI?=?0.08–012) and more diverse (A_R 95% CI SG?=?4.7–5.6, RI?=?3.1–3.3, FI?=?3.0–3.3; p?=?0.001) than other Victorian koala populations, supporting the premise that koalas in the South Gippsland region are part of a remnant population, not derived from translocated island stock. These results were also supported by mitochondrial data where eight haplotypes (Pc4, Pc17, Pc26, Pc27, and Pc56–Pc59) were identified in South Gippsland while a single haplotype (Pc27) was found in all island koalas tested. Compared to other Victorian koala populations, greater genetic diversity found in South Gippsland koalas, may provide this population with a greater chance of survival in the face of future environmental pressures. The South Gippsland koala population is, therefore, of high conservation significance, warranting the implementation of strategies to conserve this population and its diversity into the future.