Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Population Genetic Diversity in the Australian ‘Seascape’: A Bioregion Approach

  • Lisa C. Pope ,

    l.pope@uq.edu.au

    Affiliation School of Biological Sciences, The University of Queensland, Brisbane, Australia

  • Cynthia Riginos,

    Affiliation School of Biological Sciences, The University of Queensland, Brisbane, Australia

  • Jennifer Ovenden,

    Affiliation Molecular Fisheries Laboratory, School of Biomedical Sciences, The University of Queensland, Brisbane, Australia

  • Jude Keyse,

    Affiliation School of Biological Sciences, The University of Queensland, Brisbane, Australia

  • Simon P. Blomberg

    Affiliation School of Biological Sciences, The University of Queensland, Brisbane, Australia

Abstract

Genetic diversity within species may promote resilience to environmental change, yet little is known about how such variation is distributed at broad geographic scales. Here we develop a novel Bayesian methodology to analyse multi-species genetic diversity data in order to identify regions of high or low genetic diversity. We apply this method to co-distributed taxa from Australian marine waters. We extracted published summary statistics of population genetic diversity from 118 studies of 101 species and > 1000 populations from the Australian marine economic zone. We analysed these data using two approaches: a linear mixed model for standardised data, and a mixed beta-regression for unstandardised data, within a Bayesian framework. Our beta-regression approach performed better than models using standardised data, based on posterior predictive tests. The best model included region (Integrated Marine and Coastal Regionalisation of Australia (IMCRA) bioregions), latitude and latitude squared. Removing region as an explanatory variable greatly reduced model performance (delta DIC 23.4). Several bioregions were identified as possessing notably high genetic diversity. Genetic diversity increased towards the equator with a ‘hump’ in diversity across the range studied (−9.4 to −43.7°S). Our results suggest that factors correlated with both region and latitude play a role in shaping intra-specific genetic diversity, and that bioregion can be a useful management unit for intra-specific as well as species biodiversity. Our novel statistical model should prove useful for future analyses of within species genetic diversity at broad taxonomic and geographic scales.

Introduction

Genetic diversity can be considered the most basic level of biodiversity [1,2]. The need to conserve genetic diversity is being increasingly recognised by management agencies globally (e.g. [3]), with growing evidence that populations with higher levels of genetic diversity have greater resilience to changing and unpredictable environments [4,5]. Unlike global patterns of species diversity, which show strong latitudinal clines across a broad range of taxa, (e.g. [6,7]), it is yet to be established whether such patterns are common, or consistent, for within species genetic diversity across taxa (but see [814]).

Many factors will influence genetic diversity within populations. Neutral population genetic diversity is expected to be proportional to long-term effective population size and mutation [15], which is influenced by a wide range of factors such as: dispersal [16], breeding system [17], historic population size [18], and present-day abundance (e.g. [19] (nuclear only); [20]). It has been suggested that some of the same factors hypothesised to influence species richness could also influence genetic diversity. Such factors are often correlated with latitude, and include: increased mutation rates or relaxed metabolism constraints caused by higher temperatures (e.g. [21,22]); more stable, older populations [23]; and higher density due to greater food availability/niche diversity (e.g. [24]). As for many other species, marine species richness is strongly correlated with latitude, declining away from the equator in coastal species [7,25]. If there is a correlation between species richness and population genetic diversity, termed the species genetic diversity correlation (SGDC; [12,26]), we expect a decline in genetic diversity away from the equator.

Conversely, recent work has demonstrated that marine community evenness (number of individuals of each species) decreased towards the equator due to the presence of a greater number of rare species [27]. We might therefore predict a weaker relationship, or even a reduction in average intra-specific genetic diversity towards the equator due to greater variance in population size. Alternatively, differences among species in dispersal characteristics (e.g. [28]), life history strategy [29], or differences among populations due to factors such as range position (central vs edge populations, e.g. [30]), may not have strong geographic patterns, making concordant geographic patterns in population genetic diversity uncommon (e.g. [11]).

Bioregions represent an attempt to delineate areas containing distinct species assemblages, common across a broad range of taxa. Such regions are likely to possess similarities in both evolutionary history and environmental characteristics, making them a useful unit for testing biogeographic hypotheses [31]. Common evolutionary processes may result in broad-scale patterns of genetic diversity [32], and regions may therefore be more informative than latitude for describing intra-specific diversity. Extensive effort has been placed into designating bioregions within the Australian marine environment, based on the distributions of fish and other species, and bathymetric and environmental variables [33]. This regional framework, the Integrated Marine and Coastal Regionalisation of Australia (IMCRA), has been used to determine representative regions for conservation and management around Australia. However, the relevance of these regions to the distribution of intra-specific genetic diversity is unknown.

In order to determine within species genetic diversity patterns at a broad geographic scale, large amounts of data are required. A synthetic analysis that makes use of existing data from previously published studies has many advantages. Most notably, information can be obtained for a large number of species across a broad geographic extent, for very little cost and low effort. Without such a comprehensive data set, the description of ecosystem scale patterns is compromised. Combining data across multiple studies involves substantial challenges that need to be recognised. These include: combining data from different genetic markers and measures of genetic diversity, uneven geographic sampling, and uneven taxonomic sampling. However, the benefit of being able to generalise across a broad range of species and geographic extent should outweigh these disadvantages, where care is taken in the analysis (see [3437]).

Australia’s marine exclusive economic zone (EEZ) is one of the most biodiverse regions on earth [38]. There has been significant progress towards cataloguing Australian marine species diversity in recent years (e.g. fishmap, Atlas of Living Australia, http://www.ala.org.au; [39]), however there remain large numbers of species and ecosystems about which little is known [39]. Identification of poorly studied regions is considered a research priority [40]. Here, we perform an analysis of a random sample of published population genetic statistics from marine species around Australia. The majority of these studies provide estimates of neutral population genetic diversity, as few studies to date have estimated functional genetic diversity. We develop a methodology to analyse multi-species genetic diversity data, and use this to identify regions of high/low genetic diversity. We use this information to suggest priorities for future population genetic research in Australian waters. In addition, we test whether there are common patterns of within species genetic diversity that correlate with bioregion or latitude

Materials and Methods

Database construction

We performed a literature search on 11th April 2012. We combined searches from five databases: Web of Science, Zoological Record, BIOSIS (Thomson Reuters), Scopus (Elsevier), and ASFA1 (Proquest). Details of these searches are included in S1 Appendix.

The following selection criteria were then applied:

  1. Marine organism with at least three populations sampled from the Australian EEZ
  2. Minimum of five samples per population
  3. Study extent greater than 20 km
  4. Information on population HE (expected heterozygosity, [41]), or h (haplotype diversity, [42])
  5. Some information on the population geographic locations, either through coordinates, place names or clear placement in a figure allowing approximate locations to be derived using Google Earth v7.1.2.2041 (http://www.google.com/earth/).
  6. Studies that potentially included multiple cryptic species or hybridization, as indicated by the primary literature, were discarded.

Studies that met these criteria were entered into a database, with each species/marker combination per paper given a unique ‘dataID’. Species were divided into ‘broad taxonomic units’ (BTUs), generally based on Class but sometimes Phylum or Order where there were insufficient species to separate Classes (see Table 1). Genetic markers were assigned to six categories: mitochondrial (mtDNA) sequence, mtDNA restriction enzyme, allozyme, microsatellite, nuclear sequence, other (which included nuclear and mitochondrial single-strand conformation polymorphism (SSCP), inter-simple sequence repeat (ISSR), and denaturing gradient gel electrophoresis (DGGE) studies).

thumbnail
Table 1. The number of populations and species included in this regional analysis of genetic diversity, relative to the number of known Australian species.

BTU represents a broad taxonomic unit, # pops, # species and % sp. represent the data used in present study (%). Known Aus. species represents the number of species in public databases as reported by [39], presented both as total number (#) and as a % of the total recorded marine species for the Australian economic zone (%). Studied (%) is the percentage of recorded Australian species used in the present study. The table does not include all major marine groups in Australia, but lists the three largest groups not studied.

https://doi.org/10.1371/journal.pone.0136275.t001

We attempted to extract all population genetic studies performed in Australian marine waters by our literature search, though inevitably some studies were missed. While our sample is not exhaustive, it should be unbiased both taxonomically and geographically, and so should reflect genuine taxonomic and geographic sampling biases. To gain an idea as to the proportion of marine population genetic studies undertaken prior to 2012 that were included in our analysis, we performed ad hoc searches of authors that occurred at high frequency in our database, and selected 100 papers. We then determined the proportion of papers identified searching by author that were missing from our original search. We also performed ad hoc searches of under represented taxonomic groups (Platyhelminthes and Annelida). To examine if there has been a change in sampling bias in recent times we implemented the literature search used in mid 2012 on the 20th of January 2015. For the papers that matched our criteria we determined taxonomic group and rough geographic region studied, in order to determine if there had been a major taxonomic or geographic shift in the more recent literature.

Geographic and taxonomic analyses

The geographic distribution of population samples we obtained was visualised as the number of populations within a radius of 100 km, and density of points per IMCRA region. Maps were produced using ArcMap10.2.2 (ESRI 2011). IMCRA provincial bioregion boundaries (version 4) were obtained from the Atlas of Living Australia website, provided by the Environmental Resources Information Network (c) Commonwealth of Australia, Australian Government Department of the Environment and Heritage 2006).

To assess species similarity among regions we performed a community-style cluster analysis using sampling region as the group of interest and sampled species as presence/absence information, as performed by Keyse et al. [43] (‘region co-sampling cluster analysis’). We calculated Euclidean distances among sites using the vegan package for R ([44,45]), and clustered them into groups using Ward’s Minimum Variance criterion, with input distances squared and output branch lengths square rooted, as recommended by Murtagh and Legendre [46]. Only regions with at least five species sampled were considered.

Estimating regional mean genetic diversity

A common approach to comparing within species diversity is to standardise data in some way (e.g. [11,47]). While this does allow direct comparison, is has several potential disadvantages. Standardised data will contain much greater ‘noise’ due to the difficulty in obtaining accurate species means, and may be less influenced by environmental gradients due to the restricted geographic distribution of many species and the limitation that all species are assigned the same mean genetic diversity. Given these restrictions, we developed and tested an alternative method to analyse unstandardised population genetic diversity data, and compared this to the standard method. We did this using a Bayesian approach, which is particularly useful when applied to complex problems with no exact analytical solution, and which can be difficult to solve using standard ‘frequentist’ techniques [48]. Due to low sample sizes in most offshore regions, which caused instability in our models, we performed this analysis for coastal regions only.

Model construction.

To allow a more direct comparison, both analysis methods 1) ‘standardised’ (Z) and 2) ‘unstandardised’ (H), were executed within a Bayesian framework using the software WinBUGS version 1.4.3 [49], called from R using the library R2WinBUGS [50]. The response variable, population genetic diversity, was measured as either haplotype diversity (h; ie. mtDNA) or expected heterozygosity (HE; i.e. nuclear DNA; [41]). ‘Standardised’ models used ‘z-scored’ population genetic diversity as the response variable. These were calculated for each dataID (study:genetic marker combination per paper) by subtracting the mean and dividing by the standard deviation for that dataID group. In the ‘unstandardised’ models, the response variable was raw population genetic diversity. A logit link function was used to link population genetic diversity to the explanatory variables. We fitted linear mixed models to the standardised data, and mixed beta regression models to the unstandardised data. Beta regression is appropriate where the response variable is bounded by 0−1, is continuous and is highly skewed, as is the case with unstandardised population genetic diversity data [51].

We tested the utility of all combinations of five explanatory variables for each analysis type. These variables were: sample size (n), marker (M), species (S), latitude (L), latitude squared (L2) and region (R). Sample size was included as an explanatory term for the dispersion of our measure of genetic diversity (following [20,52]). Genetic marker type (M) was included as a factor with six levels. To attempt to account for similarities in genetic diversity among closely related species, and to weight each species equally, we constructed a correlation matrix from the branch lengths of a taxonomic ‘tree’, based on membership to Class, Order, Family, Genus and species, constructed using Grafen’s [53] branch lengths, with the R package ape and assuming a Brownian motion model of evolution [54]. The matrix was used as a multivariate normal prior for the species effect (S). We assessed both a linear and squared relationship between population genetic diversity and latitude, as covariates. Regional effects can be spatially autocorrelated if the factors causing high or low genetic diversity operate at a broader scale than the selected regions [55]. Potential spatial autocorrelation among IMCRA regions was taken into account using a conditional autoregression (CAR) prior [56]. Adjoining regions were identified and a distance matrix based on the extent of sharing of regional boundaries was constructed using WinBUGS [49].

The model for standardized data (Z, linear mixed effects model) can be formally written as:

The model for unstandardized data (H, beta–regression model) can be written as: where Xi contains the marker matrix, and latitude variables where appropriate, R is the region effects (with Ri CAR prior) and S is the species effect (with Si multivariate Normal prior, with zero mean and correlation matrix derived from the ‘taxonomic tree’). For further details please see S2 Appendix.

Model selection and checking.

All models were assessed based on: 1) Deviance Information Criteria score (DIC; [57]) and 2) p-values from posterior predictive tests (PPP). The DIC was calculated within WinBUGS. DIC weights were determined as calculated for AIC weights [58]. We used posterior predictive p values as a form of model checking [59]. P values from posterior predictive tests were determined by using our models to generate expected population genetic diversity values, and comparing these to the observed data [60]. Extreme p values (e.g p < 0.05, p > 0.95) indicate a contradiction of the model by the observed data [61]. We also performed sensitivity tests to determine the influence of various priors, and chose those that gave the lowest DIC scores. Model test runs were performed using three chains for 100 000 iterations, thinning every 100 iterations and with a burnin of 5000. All models showed high levels of convergence based on low autocorrelation, Brooks-Gelman-Rubin scores [62] and mixing profiles of chain traces after this number of iterations.

Once the best model was determined (based on DIC and PPP results), a final run of 500 000 iterations, with thinning every 100 iterations and a burnin of 50 000 was performed. We report the mean as a measure of central tendency and 95% central posterior density intervals for regional mean genetic diversity and latitude effects where appropriate. Models were run for data from coastal regions only, ensuring all data sets used had a minimum of 3 populations once outer regions were excluded. To assess variation in geographic patterns between different taxonomic groups, in addition to the coastal region data set containing all species (959 populations, 84 species), models were also run for a dataset comprised of Actinopterygii and Chondrichthyes (ray-finned and cartilaginous fish) (299 populations, 33 species), and Mollusca (148 populations, 9 species), separately.

Results

From 118 articles in the primary literature that met our criteria, we obtained 153 data records (each given a unique data ID; multiple dataIDs from a paper resulted from multiple markers and/or multiple species per article). From these we obtained 1154 population data points containing either HE or h from 101 species (S1 Table). Coastal regions contained data from 84 species, comprising 959 populations. The percentage of information gained from the six different marker categories was as follows: allozyme 32.8%; microsatellite 31.4%; mtDNA sequence 23.5%; mtDNA restriction enzyme 8.0%; nuclear sequence 1.2%; other 2.9%. The average genetic diversity per marker category is shown in S2 Fig.

From 100 papers identified through an ad hoc search by author, 76 were already captured in our database, suggesting we may have missed 24% of relevant literature in the time period we surveyed. Additional ad hoc searches of taxonomic groups that were under represented in our database (Platyhelminthes and Annelida) did not yield any new articles. Our search of articles from mid 2012 to the 20th of January 2015 detected 978 new publications, of which 57 matched our criteria. These included major taxonomic groups in similar proportions to those identified previously, with Molluscs (7/57) and Crustaceans (6/57) still under represented based on the number of known species around Australia. However, 25% of studies had sites located in either the Australian Bight or Joseph Bonaparte Gulf, suggesting that progress is being made towards a better understanding of these regions.

Geographic and taxonomic distribution of population genetic studies

Sample sites for population genetic studies were not evenly distributed around Australia (Fig 1A and 1B). The majority of sites occurred in coastal areas, particularly the eastern coast. Few of the non-coastal IMCRA regions were well sampled, although four out of 24 non-coastal regions (2, 3, 10, 18) had at least 15 sampling points. The four least studied coastal regions in terms of density were: 26, 32, 25, 27 (Fig 1A and 1B). These included the most northern regions apart from the Great Barrier Reef (GBR): the Gulf of Carpentaria and the ‘Top End’ (25), the Joseph Bonaparte Gulf (26) and northern Western Australia (WA) (27), as well a southern region in the Great Australian Bight (32).

thumbnail
Fig 1. The spatial distribution of populations used to estimate intraspecific genetic diversity for this study.

A) The sampling density of populations for estimates of genetic diversity a) over a 100km circle, determined using spatial analyst in ArcMap (ESRI). b) IMCRA region. IMCRA regions are identified by number: regions 1–24 are oceanic or outer shelf, regions 25–41 are coastal and are numbered anti-clockwise starting in the Gulf of Carpentaria. c) A dendrogram of the co-sampling of species among regions, constructed using Euclidean distance. Only regions with at least five species are shown. All maps were produced using GDA94 Australian Albers projection: EPSG 3577.

https://doi.org/10.1371/journal.pone.0136275.g001

The region co-sampling cluster analysis indicated that neighbouring regions commonly shared surveyed species (Fig 1C). An exception to this was the Great Barrier Reef, regions 40 and 41, which were dissimilar in species inclusion. Region 40 (mid-South GBR) appeared to be rarely co-sampled with any other region, and thus formed a distinct group. There was little co-sampling across the tropical/temperate regions. Regions in the south formed a separate sampling cluster (31−39), with regions along the south-east coast forming a distinct cluster within this group (37−39). There was high co-sampling across the top of Australia (25−27), but less co-sampling across the three ‘tropical’ regions: the ‘Top End’; the GBR; and mid WA (including Ashmore and Scott Reefs, region 2; and Rowley Shoals, region 3).

The taxonomic spread of species targeted for population genetic studies was wide, however most ‘broad operational taxonomic units’ contained small sample sizes (Table 1). The top four taxonomic units based on number of species were: Actinopterygii (41), Anthozoa (10), Chondricthyes (8), Mammalia (6). We compared the species studied in our database to those recorded for the extended Australian marine EEZ ([39]; see Table 1). Groups with a difference of > 5%, which could potentially be considered as well represented in the database included: Pisces, and Cnidirians. Groups with a difference of < 5%, which could be considered to be poorly represented were: Crustacea, and Mollusca (Table 1).

Patterns of regional marine population genetic diversity

No model for standardised data passed the posterior predictive tests for fit (Table 2, S2 Table). This indicated that none of these models were able to adequately reproduce the observed data. The best model for standardised data that was closest to passing our PPP test (in that the observed value fell closest to the expected distribution) contained sample size and regional effects (SnR). No region had a mean value that was notably different from zero for all standardised models that contained sample size and region, based on 95% central posterior density intervals of regional means, (including the best model). However, standardised models that did not contain sample size (e.g. SR), consistently found that region 25 (Gulf of Carpentaria) had notably high genetic diversity, and region 27 (Nth West coastal) had notably low genetic diversity (6 models).

thumbnail
Table 2. Results from models containing increasing numbers of explanatory variables, using both unstandardised population genetic diversity measures (Beta) and standardized genetic diversity (Z) as the response.

Explanatory variables are as follows: S = Species, M = genetic marker, L = latitude, L2 = latitude squared, R = IMCRA region, n = sample per population. DIC is the deviance information criterion, and delta represents the difference of a model’s DIC from the best model. PPP is the p value from our posterior predictive tests ([60]) derived from a comparison of the observed data with data generated from the model. Extreme probability values (e.g. p < 0.05, p > 0.95) indicate a contradiction of the model by the observed data [61]. DICw represents DIC weight calculated according to [58]. The best model from each analysis is indicated with an asterisk.

https://doi.org/10.1371/journal.pone.0136275.t002

The best overall model, based on both DIC and PPP scores, was the ‘unstandardised’ model that contained species, marker, latitude, latitude squared and regional effects (S + M + n + L + L2 + R; Table 2). This model had a high DIC weight (0.86), indicating it was a great improvement over all other models. The model containing region without latitude (S+M+n+R delta DIC = 18.3), was a better fit than the model with latitude and without region (S+M+n+L+L2 DIC = 23.4), however DIC weight and PPp results indicated that the model containing both region and latitude variables was superior.

For all unstandardised models where region was included as an explanatory variable, regions 26 (Joseph Bonaparte Gulf) and 29 (Central West coastal) had notably high genetic diversity and region 39 (Southern Qld / N NSW) was notably low (Table 3, Fig 2). In some regional models regions 36 (Tasmania coastal), 31 and 32 (Great Australian Bight) were also notably low. Region 27 (NW coastal) was notably high in the best model only. Where both latitude and latitude squared were included together as explanatory variables, diversity increased towards the equator. Latitude squared had a negative relationship indicating a hump in diversity across the latitudes studied (~ 9.4 to 43.7°S).

thumbnail
Fig 2. Regional genetic diversity determined using a beta-regression approach with unstandardized population genetic diversity data.

The larger map represents regional means, the smaller map the standard deviation within regions. Regions were colour coded by dividing values for each map into five equal intervals; values falling in the upper 20% interval are coloured red/pink, the lowest 20% dark green/light blue etc.

https://doi.org/10.1371/journal.pone.0136275.g002

thumbnail
Table 3. The mean intra population genetic diversity for IMCRA regions (see Fig 2), estimated from models using either standardised (Z) or un-standardised data (H).

The best model for standardised genetic data contained a sample size effect, and regional effect in addition to species (SnR). The best model for unstandardized data contained in addition genetic marker, latitude, and latitude squared effects (SnML2R). Results from beta models for fish and molluscs only are also shown. Significant means (those whose MC 95% CI do not span zero) are indicated with an asterisk. Means for regions with less than 10 populations sampled are not shown.

https://doi.org/10.1371/journal.pone.0136275.t003

There was little correspondence between regional means estimated using the standardised approach, compared to the unstandardised approach based on correlations between regional means (R2 = 0.02, p = 0.26)(Table 3, S1 Fig). Focusing on the unstandardised approach, regional means for ray-finned and cartilaginous fish were similar to those found for the entire data set (for regions where at least 10 ‘fish’ populations were sampled (n = 12, R2 = 0.51, p = 0.005). However, when comparing regions where Molluscs were sampled with at least 10 populations there was no relationship between the entire dataset and Molluscs (n = 7, R2 = 0.14, p = 0.22; S1 Fig).

Discussion

This study represents the largest synthesis to date of marine population genetic studies at an Australia wide scale. We have developed statistical methods to allow the regional analysis of population genetic data from multiple sources. This well parameterised Bayesian method appears to be an improvement over the more commonly used ‘standardisation’ approach, based on posterior predictive testing. Region was found to be a useful explanatory variable in describing the broad scale distribution of genetic diversity, in addition to latitude. Genetic diversity increased towards the equator, however the relationship between diversity and latitude was best described as ‘humped’, rather than linear. Our study has identified several regions as having either higher or lower genetic diversity than those around them (particularly high diversity in the Joseph Bonaparte Gulf (IMCRA 26) and low diversity in southern Queensland, northern NSW (IMCRA 39). Describing this regional variation represents a first step towards understanding the broad scale processes that might be driving them.

Geographic and taxonomic representation

The distribution of information on population genetic diversity was uneven around Australia. Sampling was primarily coastal and highest on the east coast of Australia, particularly the Great Barrier Reef. While well studied, the southern GBR has not often been co-sampled with other regions,. This lack of co-sampling of the GBR with other regions has previously been identified based on a dataset of tropical species spanning the Indo-Pacific [43]. The ‘Top End’, Gulf of Carpentaria, Joseph Bonaparte Gulf and the Great Australian Bight were relatively poorly studied.

There are over 32000 species recorded in major databases as occurring in the Australian marine economic zone, though, based on the rate of discovery of undescribed species, the actual number is far higher [39]. While we attempt to describe regional genetic diversity across a wide range of species, this study captures an extremely small part of recorded species diversity (101 species, 0.3%). The majority of species included in our analyses were ray-finned and cartilaginous fish (48%), which were over-represented, in that only ~16% of recorded Australian marine species are Pisces. Yet despite this <1% of recorded ‘fish’ species were included in our database. Cnidirians (10.9% vs 5.3%) and Mammals (~10% of known species) were also relatively well represented according to the proportions in our database compared to that of known Australian species. On the same basis, Crustacea (5.9% vs 19.3%) and Mollusca (8.9% vs 25.9%) were under represented.

Genetic information, in combination with other methods, is extremely useful for defining biological stocks [63,64], an understanding of which is necessary to achieve sustainable harvesting [65]. Australia is moving towards the management of its commercial marine resources as biological stocks rather than by political boundaries [66], however for many species information is lacking. Of the 49 most important commercial taxa in Australia identified in a National Stock Structure of Australia report [66], only 20% were found to have existing population genetic data that fitted our search categories (as defined in the methods).

Latitudinal patterns of genetic diversity

Genetic diversity followed the expected negative relationship with latitude, increasing towards the equator. Species richness has a strong correlation with latitude (reviewed in [67,68]), which is in turn correlated with many environmental variables, such as sea surface temperature, and coastal complexity [7]. While not directly tested here, as species richness is correlated with latitude, our results provide some support for the ‘species genetic diversity correlation’ in the marine environment [12,26]. However, models with a non-linear relationship between latitude and genetic diversity were a better fit than a simple linear relationship. A ‘humped’ relationship between genetic diversity and latitude has previously been reported in a study of tropical birds [10] where it was suggested this pattern might be caused by reduced diversity at range edges occurring over both ecological (see [69]) and evolutionary (i.e. geometric contraints; [70]) time frames. This argument requires congruence in the range edge of many species, and/or that range edge effects are extensive, and thus may only apply to certain regions and species groups.

Regional genetic diversity

IMCRA region was a useful concept to explain the distribution of intra-specific genetic variation, both for standardised and unstandardised data. Regional patterns persisted after latitude had been accounted for, and so additional explanations for these differences are needed. Three regions were found to possess ‘high’ levels of genetic diversity (The Joseph Bonaparte Gulf (26), NW coast (27) and the Central West coast (29)), and one possessed ‘low’ diversity (Queensland/northern NSW (39), based on the best fitting unstandardised model.

Areas of high diversity might represent regions of admixture, where distinct evolutionarily significant units overlap (e.g. ‘centre of overlap’, [71]). High within species genetic divergence has been reported either side of Torres Strait (e.g. [72]), which was impassable at times of low sea level, and disparate species groups have been reported for some species along the western coast [7375]. However, neither high diversity region is exactly congruent with these potential ‘admixture zones’. A related explanation for the low diversity observed in Sth Queensland/Nth NSW (IMCRA 39) may be the presence of large numbers of range edge populations, potentially lowering average intra-specific genetic diversity. However, meta-analyses suggest that range edge populations do not always contain lower diversity [30], weakening this argument. Now that patterns of regional intra-specific genetic diversity have been described, further testing of their generality can take place, and hypotheses to explain their origin can be generated and tested.

Limitations of our approach

It was surprisingly difficult to limit a search to the marine environment in a specific geographic area, despite attempting to do so across five different databases. While we may have extracted data from approximately 75% of relevant publications, it is unlikely that there was a strong taxonomic bias in this missing data due to our search methodology. It is possible that there was some geographic bias, as some databases required the use of place names to restrict to location. However, as we included the terms ‘Joseph Bonaparte Gulf’ and ‘Australian Bight’ we believe our finding that these regions were under sampled remains valid. Our survey of recent literature suggests that while taxonomic biases remain, some progress is being made towards a better understanding of the Great Australian Bight and Joseph Bonaparte Gulf. We hope that this first attempt to summarise population genetic research in the Australian marine environment will provide a useful starting point for future researchers.

The majority of population genetic analyses performed on wild species to date have utilised neutral genetic markers, due to their greater tractability for analyses. It may be that for conservation purposes, functional genetic diversity is of greater relevance, as this diversity will result in a change in an organism’s phenotype, however there generally appears to be weak correlation between functional and neutral genetic diversity [76,77]. Furthermore, while our focus here in on average population genetic diversity, with the idea that high diversity populations are of higher priority for conservation, it has also been argued that a suite of populations should be selected based on the overall proportion of the species’ genetic diversity that they represent, potentially prioritising populations that have low, but unique, genetic diversity [78,79]. It was not possible for us to consider population ‘unique’ diversity, or functional genetic diversity with the summary data we obtained. Future research may combine both these approaches across multiple species to assist in identification of priority regions for conservation.

The majority of studies considering broad patterns of population genetic diversity, by necessity, focus on common, widely distributed species (e.g. [9]). Rare species may possess lower levels of genetic diversity due to lower effective population size. Our likely bias towards common species may have inflated regional means in tropical areas, where a greater proportion of rare species occur [27]. Thus, the possibility remains that a different relationship between latitude and ‘average’ intraspecific genetic diversity could be observed if the sampling of species across regions was a better representation of both the common and rare species present.

Taxonomic groups were not evenly represented across regions and, despite the inclusion of species as an explanatory variable, this taxonomic bias will have influenced our results. This is evidenced by the low regional mean obtained for molluscs from the highest diversity region in the entire and ‘fish’ only data sets (Table 1, S1 Fig). Differences in regional diversity among taxa may result from many factors such as: life history differences [29], dispersal abilities [16], body size [20], and mutation rates [80], all of which often have a taxonomic component. Increasing species numbers, and improving the taxonomic spread would improve results. Depending on the management objectives, selecting particular species subsets may be of greater relevance.

Pelagic marine species richness was shown to peak in the mid latitudes (20–40 degrees), whereas coastal species peaked in the western Pacific, close to the equator [7]. Our study focused on coastal species and regions, as data for outer regions were limited. Pelagic species may show different regional trends to those we describe for coastal species, and more information on these species might provide greater insight into evolutionary processes.

We focus on regional means, and do find some significant results, yet the large standard deviations we obtained indicate extensive variation exists within regions (Table 3, Fig 2). While our results support the occurrence of general patterns [12], given the preceding caveats, we acknowledge that the extent to which these patterns apply across all species remains unsure [11].

A new statistical model for genetic diversity indices

Posterior predictive testing provided greater support for the unstandardised Bayesian models, indicating that these models were better able to reproduce results comparable to the observed data. It is perhaps not surprising that there was little correspondence between estimates of regional means using the two methods. Using the standardised approach, we removed the effects of marker and species by giving each study the same mean. While this does allow studies to be compared directly, we lose information as to whether some species have higher ‘absolute’ levels of genetic diversity. It is possible that the Bayesian analysis has captured geographic variation in ‘absolute’ genetic diversity to some extent, allowing a broader exploration of geographic variation in intra-specific diversity.

Additionally, to accurately standardise data requires a good understanding of the underlying mean and distribution for the species. For studies with few populations sampled, z-scored data will not be accurate, and this inaccuracy may detract from the underlying signal obtained from better sampled species. This may be the reason our standardised models failed all posterior predictive tests. If sufficient populations were available to calculate a reliable mean and standard deviation for a study, z-score transformed data may provide useful insight into whether areas possess unusually high or low genetic diversity relative to the species as a whole. However, if only smaller sample sizes are available, this may generate so much noise in the data that the results may be of limited use.

While the unstandardised Bayesian approach represents an advance in statistical methods for combining summary statistics calculated from diverse genetic marker types and species, further improvements are possible. If more data were available, a larger number of genetic marker categories might improve model performance. McCusker and Bentzen [20] found substantial differences between genetic diversity estimates derived from mtDNA control region sequences, as compared to other sections of mtDNA, which they attributed to a faster mutation rate. A different measure of population genetic diversity, such as nucleotide diversity, might provide more insight into evolutionary processes. Nucleotide diversity provides information on both the number of haplotypes and how divergent those haplotypes are within a population. Areas of admixture would be highlighted using this type of measure, which may provide greater insight into evolutionary processes. Finally, a phylogenetic correction for the species effect, the inclusion of a larger number of life history traits [29], and factors such as range position (edge/middle) might help explain a greater proportion of the variation observed.

Management implications and future research priorities

  1. Data collected for this analysis suggests that there are gaps in available information for some marine regions around Australia, particularly for most offshore regions, the Great Australian Bight and Gulf of Carpentaria (though this may be improving), and for taxonomic groups such as Crustacea and Mollusca. Analysis of co-sampling of species across regions indicated that species from the southern GBR (region 40) were often not sampled in other regions. Only 20% of Australia’s top 49 commercial marine species were included in this analysis. Expansion of this database would allow these gaps to be confirmed and prioritised for future research.
  2. In addition to information on intraspecific genetic diversity, population genetic data on connectivity is often included in population genetic studies, in the form of summary statistics such as FST. While careful statistical analysis would be required to take into account factors such as distance between populations, and population genetic diversity and its influence on estimates of FST, this wealth of information could allow a broad scale insight into connectivity among multiple species around Australia. We advocate public data archiving for marine genetic studies that make population datasets more visible and enable deeper connections between data, regions and study species–including analysis of connectivity [81]. We encourage the use and expansion of the database associated with this paper.
  3. Information on average regional genetic diversity could be incorporated into the decision framework for the designation of marine park areas, with the potential for increased protection afforded to areas of consistently high diversity. For example, the area around Joseph Bonaparte Gulf, was identified as possessing high genetic variation, yet this region is under-represented in the existing Commonwealth reserve system [82].
  4. Genetic marker systems that assay a larger proportion of the genome are becoming increasingly cost-effective for non-model organisms. Increased use of such datasets, particularly if there was also geographic consistency in sampling locations across species, may provide a more direct avenue towards describing broad scale patterns of population genetic diversity and connectivity in the future by using estimates that are more directly comparable and accurate, and include both functional and neutral estimates of population genetic diversity e.g. [83].

Supporting Information

S2 Appendix. Winbug code for the best unstandardised and standardised models.

https://doi.org/10.1371/journal.pone.0136275.s002

(ZIP)

S3 Appendix. Raw data file containing article ID, marker category, individual population heterozygosity values, and xy coordinates.

See associated readme.txt for a full explanation of column headings.

https://doi.org/10.1371/journal.pone.0136275.s003

(ZIP)

S1 Fig. Regional means for standardised data, and for fish and molluscs only, using unstandardised genetic diversity in a Bayesian model.

https://doi.org/10.1371/journal.pone.0136275.s004

(EPS)

S2 Fig. Population heterozygosity or haplotype diversity for each genetic marker category.

https://doi.org/10.1371/journal.pone.0136275.s005

(TIFF)

S1 Table. A list of the species and references from which population genetic diversity data were obtained.

https://doi.org/10.1371/journal.pone.0136275.s006

(DOCX)

S2 Table. Results from all models for both unstandardised (H) and standardized (Z) population genetic diversity measures.

https://doi.org/10.1371/journal.pone.0136275.s007

(DOCX)

Acknowledgments

We would like to acknowledge UQ for providing a part-time women’s fellowship and funding for the main author, following a career break. We thank Greg George from the UQ Biological Sciences library, for help in designing our literature search. We thank Jake Gratten for R-coding advice and Shane Lavery, Dean Blower and Penny Mills for comments on earlier drafts of the manuscript.

Author Contributions

Conceived and designed the experiments: CR JO LCP SPB. Performed the experiments: LCP. Analyzed the data: SPB LCP. Contributed reagents/materials/analysis tools: LCP SPB JK. Wrote the paper: LCP CR JO JK SPB.

References

  1. 1. Wilcox BA. In situ conservation of genetic resources: determinants of minimum area requirements. In: McNeely JA, Miller KR, editors. National Parks, Conservation and Development, Proceedings of the World Congress on National Parks. Washington, DC: Smithsonian Institution Press; 1984. p. 18–30.
  2. 2. Campbell AK. Save those molecules! Molecular biodiversity and life. Journal of Applied Ecology. 2003;40(2): 193–203. pmid:ISI:000182118700001
  3. 3. UNCSD. The future we want. Rio+20; Rio de Janeiro, Brazil; 2012.
  4. 4. Hoffmann AA, Sgro CM. Climate change and evolutionary adaptation. Nature. 2011;470(7335): 479–85. Epub 2011/02/26. pmid:21350480
  5. 5. Sgro CM, Lowe AJ, Hoffmann AA. Building evolutionary resilience for conserving biodiversity under climate change. Evol Appl. 2011;4(2): 326–37. pmid:25567976; PubMed Central PMCID: PMC3352557
  6. 6. Hillebrand H. On the generality of the latitudinal diversity gradient. American Naturalist. 2004;163(2): 192–211. pmid:14970922
  7. 7. Tittensor DP, Mora C, Jetz W, Lotze HK, Ricard D, Berghe EV, et al. Global patterns and predictors of marine biodiversity across taxa. Nature. 2010;466(7310): 1098–101. Epub 2010/07/30. pmid:20668450
  8. 8. Gillman LN, Keeling DJ, Ross HA, Wright SD. Latitude, elevation and the tempo of molecular evolution in mammals. Proceedings of the Royal Society B: Biological Sciences. 2009;276(1671): 3353–9. pmid:WOS:000268798600017; PubMed Central PMCID: PMC2817169
  9. 9. Adams RI, Hadly EA. Genetic diversity within vertebrate species is greater at lower latitudes. Evolutionary Ecology. 2013;27(1): 133–43. pmid:WOS:000312876500010
  10. 10. Miller MJ, Bermingham E, Klicka J, Escalante P, Winker K. Neotropical birds show a humped distribution of within-population genetic diversity along a latitudinal transect. Ecology Letters. 2010;13(5): 576–86. pmid:WOS:000276790200006
  11. 11. Taberlet P, Zimmermann NE, Englisch T, Tribsch A, Holderegger R, Alvarez N, et al. Genetic diversity in widespread species is not congruent with species richness in alpine plant communities. Ecology Letters. 2012;15r(12): 1439–48. Epub 2012/09/26. pmid:23006492
  12. 12. Vellend M, Lajoie G, Bourret A, Murria C, Kembel SW, Garant D. Drawing ecological inferences from coincident patterns of population- and community-level biodiversity. Molecular Ecology. 2014;23(12): 2890–901. Epub 2014/04/23. pmid:24750409
  13. 13. April J, Hanner RH, Mayden RL, Bernatchez L. Metabolic rate and climatic fluctuations shape continental wide pattern of genetic divergence and biodiversity in fishes. PLOS ONE. 2013;8(7): e70296. Epub 2013/08/08. pmid:23922969; PubMed Central PMCID: PMC3726496
  14. 14. Martin PR, McKay JK. Latitudinal variation in genetic divergence of populations and the potential for future speciation. Evolution. 2004;58(5): 938–45. pmid:WOS:000221632700004
  15. 15. Kimura MT. The neutral theory of molecular evolution. Cambridge: Cambridge University Press; 1983.
  16. 16. Ryman N, Allendorf FW, Jorde PE, Laikre L, Hossjer O. Samples from subdivided populations yield biased estimates of effective size that overestimate the rate of loss of genetic variation. Molecular Ecology Resoures. 2014;14(1): 87–99. pmid:WOS:000328372500010
  17. 17. Nunney L. The influence of mating system and overlapping generations on effective population size. Evolution. 1993;47(5): 1329–41. pmid:ISI:A1993NM70700004
  18. 18. Frankham R. Effective population size/adult population size ratios in wildlife: A review. Genetical Research. 1995;66(2): 95–107.
  19. 19. Bazin E, Glemin S, Galtier N. Population size does not influence mitochondrial genetic diversity in animals. Science. 2006;312(5773): 570–2. Epub 2006/04/29. pmid:16645093
  20. 20. McCusker MR, Bentzen P. Positive relationships between genetic diversity and abundance in fishes. Molecular Ecology. 2010;19(22): 4852–62. pmid:WOS:000283986700004
  21. 21. Evans KL, Warren PH, Gaston KJ. Species-energy relationships at the macroecological scale: a review of the mechanisms. Biological Reviews. 2005;80(1): 1–25. pmid:WOS:000227418900001
  22. 22. Gillman LN, Wright SD. Species richness and evolutionary speed: the influence of temperature, water and area. Journal of Biogeography. 2014;41(1): 39–51. pmid:WOS:000328271500005
  23. 23. Jansson R, Dynesius M. The fate of clades in a world of recurrent climatic change: Milankovitch oscillations and evolution. Annual Review of Ecology and Systematics. 2002;33: 741–77. pmid:WOS:000180007000025
  24. 24. Gaston KJ. Global patterns in biodiversity. Nature. 2000;405(6783): 220–7. Epub 2000/05/23. pmid:10821282
  25. 25. Last PR, White WT, Gledhill DC, Pogonoski JJ, Lyne V, Bax NJ. Biogeographical structure and affinities of the marine demersal ichthyofauna of Australia. Journal of Biogeography. 2011;38(8): 1484–96. pmid:WOS:000292696500005
  26. 26. Vellend M. Island biogeography of genes and species. The American Naturalist. 2003;162(3): 358–65. Epub 2003/09/13. pmid:12970843
  27. 27. Stuart-Smith RD, Bates AE, Lefcheck JS, Duffy JE, Baker SC, Thomson RJ, et al. Integrating abundance and functional traits reveals new global hotspots of fish diversity. Nature. 2013;501(7468): 539–42. Epub 2013/09/27. pmid:24067714
  28. 28. Marshall JL, Camp CD. Environmental correlates of species and genetic richness in lungless salamanders (family Plethodontidae). Acta Oecol. 2006;29(1): 33–44. pmid:WOS:000235545700005
  29. 29. Romiguier J, Gayral P, Ballenghien M, Bernard A, Cahais V, Chenuil A, et al. Comparative population genomics in animals uncovers the determinants of genetic diversity. Nature. 2014;515(7526): 261–3. Epub 2014/08/21. pmid:25141177
  30. 30. Eckert CG, Samis KE, Lougheed SC. Genetic variation across species' geographical ranges: the central-marginal hypothesis and beyond. Molecular Ecology. 2008;17(5): 1170–88. Epub 2008/02/28. pmid:18302683
  31. 31. Spalding MD, Fox HE, Allen GR, Davidson N, Ferdana ZA, Finlayson M, et al. Marine ecoregions of the world: a bioregionalization of coastal and shelf areas. BioScience. 2007;57(7): 573–83.
  32. 32. Avise JC. Phylogeography: The History and Formation of Species. Cambridge, MA: Harvard University Press; 2000.
  33. 33. Last PR, Lyne VD, Williams A, Davies CR, Butler AJ, Yearsley GK. A hierarchical framework for classifying seabed biodiversity with application to planning and managing Australia’s marine biological resources. Biological Conservation. 2010;143(7): 1675–86.
  34. 34. O'Corry-Crowe G. Climate change and the molecular ecology of arctic marine mammals. Ecological Applications. 2008;18(2): S56–S76. pmid:WOS:000255381500005
  35. 35. Weersing K, Toonen RJ. Population genetics, larval dispersal, and connectivity in marine systems. Marine Ecology Progress Series. 2009;393: 1–12. pmid:WOS:000272187800001
  36. 36. Selkoe KA, Henzler CM, Gaines SD. Seascape genetics and the spatial ecology of marine populations. Fish and Fisheries. 2008;9(4): 363–77. pmid:WOS:000261060300003
  37. 37. Riginos C, Douglas KE, Jin Y, Shanahan DF, Treml EA. Effects of geography and life history traits on genetic differentiation in benthic marine fishes. Ecography. 2011;34(4): 566–75. pmid:WOS:000293031700004
  38. 38. Costello MJ, Coll M, Danovaro R, Halpin P, Ojaveer H, Miloslavich P. A census of marine biodiversity knowledge, resources, and future challenges. PLOS ONE. 2010;5(8): e12110. Epub 2010/08/07. pmid:20689850; PubMed Central PMCID: PMC2914025
  39. 39. Butler AJ, Rees T, Beesley P, Bax NJ. Marine biodiversity in the Australian region. PLOS ONE. 2010;5(8): e11831. Epub 2010/08/07. pmid:20689847; PubMed Central PMCID: PMC2914019
  40. 40. Mapstone B, Appleford P, Broderick K, Connolly R, Higgins J, Hobday A, et al. National Climate Change Adaptation Research Plan for Marine Biodiversity and Resources. National Climate Change Adaptation Research Facility, editor. Gold Coast2010. 68 p.
  41. 41. Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89(3): 583–90. pmid:17248844; PubMed Central PMCID: PMC1213855
  42. 42. Nei M, Tajima F. DNA polymorphism detectable by restriction endonucleases. Genetics. 1981;97: 145–63. pmid:6266912
  43. 43. Keyse J, Crandall ED, Toonen RJ, Meyer CP, Treml EA, Riginos C. The scope of published population genetic data for Indo-Pacific marine fauna and future research opportunities in the region. Bulletin of Marine Science. 2014;90(1): 47–78.
  44. 44. Oksanen J, Guillaume Blanchet F, Kindt R, Legendre P, Minchin PR, O'Hara RB, et al. Vegan: Community Ecology Package. R package version 2.0–10. Available: http://CRAN.R-project.org/package=vegan2013.
  45. 45. R Core Team. R: A language and environment for statistical computing. Vienna, Austria. R Foundation for Statistical Computing; 2014. Available: http://www.R-project.org/.
  46. 46. Murtagh F, Legendre P. Ward's hierarchical clustering method: clustering criterion and agglomerative algorithm. Available: http://arxivorg/abs/11116285. 2011.
  47. 47. Vandergast AG, Perry WM, Lugo RV, Hathaway SA. Genetic Landscapes GIS Toolbox: tools to map patterns of genetic divergence and diversity. Molecular Ecology Resoures. 2011;11(1): 158–61. Epub 2011/03/25. pmid:21429115
  48. 48. Lawson AB, Browne WJ, Vidal Rodeiro CL. WinBUGS Basics. Disease Mapping with WinBUGS and MLwiN: John Wiley & Sons, Ltd; 2003. pp. 45–73.
  49. 49. Lunn DJ, Thomas A, Best N, Spiegelhalter D. WinBUGS—a Bayesian modelling framework: concepts, structure, and extensibility. Statistics and Computing. 2000;10(4): 325–37. pmid:ISI:000089242200005
  50. 50. Sturtz S, Ligges U, Gelman A. R2WinBUGS: A Package for Running WinBUGS from R. Journal of Statistical Software. 2005;12(3): 1–16.
  51. 51. Ferrari SLP, Cribari-Neto F. Beta regression for modelling rates and proportions. Journal of Applied Statistics. 2004;31(7): 799–815.
  52. 52. Pinsky ML, Palumbi SR. Meta-analysis reveals lower genetic diversity in overfished populations. Molecular Ecology. 2014;23(1): 29–39. pmid:WOS:000330950900006
  53. 53. Grafen A. The phylogenetic regression. Philosophical Transactions of the Royal Society of London Series B Biological Sciences. 1989;326(1233): 119–57 pmid:2575770
  54. 54. Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20(2): 289–90. Epub 2004/01/22. pmid:14734327
  55. 55. Figueroa-Zuniga JI, Arellano-Valle RB, Ferrari SLP. Mixed beta regression: a Bayesian perspective. Computational Statistics and Data Analysis. 2012; arXiv:1201.2375v4.
  56. 56. Besag J, York J, Mollie A. Bayesian image restoration, with two applications in spatial statistics Annals of the Institute of Statistical Mathematics. 1991;43(1): 1–20. pmid:WOS:A1991FH21800001
  57. 57. Spiegelhalter DJ, Best NG, Carlin BP, Van der Linde A. Bayesian measures of model complexity and fit (with discussion). Journal of the Royal Statistical Society, Series B. 2002;64(4): 583–616.
  58. 58. Burnham KP, Anderson DR. Model selection and mul-timodel inference: A practical information-theoretic approach. New York: Springer-Verlag; 2002.
  59. 59. Gelman A, Carlin JB, Stern HS, Dunson DB, Vehtari A, Rubin DB. Bayesian data analysis. 3rd ed. Dominici F, Faraway JJ, Tanner M, Zidek J, editors. Boca Raton, FL: Chapman & Hall/CRC; 2014.
  60. 60. Gelman A, Carlin JB, Stern HS, Rubin DB. Bayesian data analysis. 2nd ed. Chatfield C, Tanner M, Zidek J, editors. Boca Raton, FL: Chapman and Hall/CRC Press; 2004. 157–96 p.
  61. 61. Rubin DB. Bayesianly justifiable and relecant frequency calculations for the applied statistician. The Annals of Statistics. 1984;12(4): 1151–72.
  62. 62. Brooks SP, Gelman A. General methods for monitoring convergence of iterative simulations. Journal of Computational and Graphical Statistics. 1998;7(4): 434–55.
  63. 63. Coyle T. Stock identification and fisheries management: the importance of using several methods in a stock identification study. Australian Society for Fish Biology. 1997: 173–82.
  64. 64. Begg GA, Waldman JR. An holistic approach to fish stock identification. Fish Res. 1999;43(1–3): 35–44. pmid:ISI:000083505000004
  65. 65. Ovenden JR, Berry O, Welch DJ, Buckworth RC, Dichmont CM. Ocean’s eleven: a critical evaluation of the role of population, evolutionary and molecular genetics in the management of wild fisheries. Fish and Fisheries. 2013;16(1): 125–59.
  66. 66. Flood M, Stobutzki I, Andrews J, Begg G, Fletcher W, Gardner C, et al., editors. Status of key Australian fish stocks reports 2012. Canberra: Fisheries Research and Development Corporation; 2012.
  67. 67. Rohde K. Latitudinal gradients in species diversity—the search for the primary cause. Oikos. 1992;65(3): 514–27. pmid:WOS:A1992KB34700020
  68. 68. Hawkins BA. Ecology's oldest pattern? Trends in Ecology & Evolution. 2001;16(8): 470–. pmid:WOS:000170031300014
  69. 69. Vucetich JA, Waite TA. Spatial patterns of demography and genetic processes across the species' range: Null hypotheses for landscape conservation genetics. Conservation Genetics. 2003;4(5): 639–45. pmid:ISI:000185261600009
  70. 70. Colwell RK, Lees DC. The mid-domain effect: geometric constraints on the geography of species richness. Trends in Ecology & Evolution. 2000;15(2): 70–6. pmid:10652559
  71. 71. Gaither MR, Rocha LA. Origins of species richness in the Indo-Malay-Philippine biodiversity hotspot: evidence for the centre of overlap hypothesis. Journal of Biogeography. 2013;40(9): 1638–48. pmid:ISI:000322981900002
  72. 72. Mirams AGK, Treml EA, Shields JL, Liggins L, Riginos C. Vicariance and dispersal across an intermittent barrier: population genetic structure of marine animals across the Torres Strait land bridge. Coral Reefs. 2011;30(4): 937–49. pmid:WOS:000296085300009
  73. 73. McCallum AW, Poore GCB, Williams A, Althaus F, O'Hara T. Environmental predictors of decapod species richness and turnover along an extensive Australian continental margin (13–35°S). Marine Ecology. 2013;34(3): 298–312.
  74. 74. O'Hara TD, Rowden AA, Bax NJ. A Southern Hemisphere bathyal fauna is distributed in latitudinal bands. Current Biology. 2011;21(3): 226–30. Epub 2011/01/25. pmid:21256017
  75. 75. Woolley SNC, McCallum AW, Wilson R, O'Hara TD, Dunstan PK, Rouget M. Fathom out: biogeographical subdivision across the Western Australian continental margin—a multispecies modelling approach. Diversity and Distributions. 2013;19(12): 1506–17.
  76. 76. Reed DH, Frankham R. How closely correlated are molecular and quantitative measures of genetic variation? A meta-analysis. Evolution. 2001;55(6): 1095. pmid:11475045
  77. 77. Whitlock R. Relationships between adaptive and neutral genetic diversity and ecological structure and functioning: a meta-analysis. J Ecol. 2014;102(4): 857–72. pmid:WOS:000338027500004; PubMed Central PMCID: PMC4142011
  78. 78. Faith DP. Conservation evaluation and phylogenetic diversity. Biological Conservation. 1992;61(1):1–10. pmid:WOS:A1992HX12300001.
  79. 79. Volkmann L, Martyn I, Moulton V, Spillner A, Mooers AO. Prioritizing populations for conservation using phylogenetic networks. PLOS ONE. 2014;9(2): e88945. pmid:WOS:000332396200030; PubMed Central PMCID: PMC3938429
  80. 80. Nabholz B, Glemin S, Galtier N. The erratic mitochondrial clock: variations of mutation rate, not population size, affect mtDNA diversity across birds and mammals. BMC Evolutionary Biology. 2009;9: 54. pmid:WOS:000264931800001; PubMed Central PMCID: PMC2660308
  81. 81. Pope LC, Liggins L, Keyes J, Carvalho SB, Riginos C. Not the time or the place: the missing spatio-temporal link in publicly available genetic data. Molecular Ecology. 2015;24(15): 3802–9. pmid:26033415
  82. 82. Barr LM, Possingham HP. Are outcomes matching policy commitments in Australian marine conservation planning? Marine Policy. 2013;42: 39–48.
  83. 83. Harrisson KA, Pavlova A, Telonis-Scott M, Sunnucks P. Using genomics to characterize evolutionary potential for conservation of wild populations. Evol Appl. 2014;7(9): 1008–25. pmid:WOS:000344595700005; PubMed Central PMCID: PMC4231592