The impact of mating systems and dispersal on fine-scale genetic structure at maternally, paternally and biparentally inherited markers

For decades, studies have focused on how dispersal and mating systems influence genetic structure across populations or social groups. However, we still lack a thorough understanding of how these processes and their interaction shape spatial genetic patterns over a finer scale (tens-hundreds of metres). Using uniparentally inherited markers may help answer these questions, yet their potential has not been fully explored. Here, we use individual-level simulations to investigate the effects of dispersal and mating system on fine-scale genetic structure at autosomal, mitochondrial and Y chromosome markers. Using genetic spatial autocorrelation analysis, we found that dispersal was the major driver of fine-scale genetic structure across maternally, paternally and biparentally inherited markers. However, when dispersal was restricted (mean distance = 100 m), variation in mating behaviour created strong differences in the comparative level of structure detected at maternally and paternally inherited markers. Promiscuity reduced spatial genetic structure at Y chromosome loci (relative to monogamy), whereas structure increased under polygyny. In contrast, mitochondrial and autosomal markers were robust to differences in the specific mating system, although genetic structure increased across all markers when reproductive success was skewed towards fewer individuals. Comparing males and females at Y chromosome vs. mitochondrial markers, respectively, revealed that some mating systems can generate similar patterns to those expected under sex-biased dispersal. This demonstrates the need for caution when inferring ecological and behavioural processes from genetic results. Comparing patterns between the sexes, across a range of marker types, may help us tease apart the processes shaping fine-scale genetic structure.

Converting quadratic entropy to diversity: Both animals and alleles are diverse, but some are more diverse than others

The use of diversity metrics has a long history in population ecology, while population genetic work has been dominated by variance-derived metrics instead, a technical gap that has slowed cross-communication between the fields. Interestingly, Rao’s Quadratic Entropy (RQE), comparing elements for ‘degrees of divergence’, was originally developed for population ecology, but has recently been deployed for evolutionary studies. We here translate RQE into a continuous diversity analogue, and then construct a multiply nested diversity partition for alleles, individuals, populations, and species, each component of which exhibits the behavior of proper diversity metrics, and then translate these components into [0,1]—scaled form. We also deploy non-parametric statistical tests of the among-stratum components and novel tests of the homogeneity of within-stratum diversity components at any hierarchical level. We then illustrate this new analysis with eight nSSR loci and a pair of close Australian marsupial (Antechinus) congeners, using both ‘different is different’ and ‘degree of difference’ distance metrics. The total diversity in the collection is larger than that within either species, but most of the within-species diversity is resident within single populations. The combined A. agilis collection exhibits more diversity than does the combined A. stuartii collection, possibly attributable to localized differences in either local ecological disturbance regimes or differential levels of population isolation. Beyond exhibiting different allelic compositions, the two congeners are becoming more divergent for the arrays of allele sizes they possess.

Non-linear growth in tree ferns, Dicksonia antarctica and Cyathea australis

Tree ferns are an important structural component of forests in many countries. However, because their regeneration is often unrelated to major disturbances, their age is often difficult to determine. In addition, rates of growth may not be uniform, which further complicates attempts to determine their age. In this study, we measured 5 years of growth of Cyathea australis and Dicksonia antarctica after a large wildfire in 2009 in south-eastern Australia. We found growth rates of these two species were unaffected by aspect and elevation but slope had a minor effect with D. antarctica growing 0.3mm faster for each additional degree of slope. Geographic location influenced growth in both species by up to 12 - 14mm/yr. The most consistent factor influencing growth rate, however, was initial height at the time of the 2009 fire; a finding consistent in both species and all geographic locations. For both tree fern species, individuals that were taller at the commencement of the study had greater overall growth for the duration of the study. This effect did not decrease even among the tallest tree ferns in our study (up to 6 metres tall). Overall, Cyathea australis averaged 73 (± 22)mm/year of growth (± 1SD), with the rate increasing 5mm/yr per metre of additional height. Dicksonia antarctica averaged 33 (± 13)mm/year, increasing by 6mm/yr/m. Growth rates dependent on initial height were unexpected and we discuss possible reasons for this finding. Variable growth rates also suggest that common age estimation methods of dividing height by average growth rate are likely to underestimate the age of short tree ferns, while overestimating the age of tall tree ferns, particularly if they have been subject to a fire.

Short- and long-term effects of habitat fragmentation differ but are predicted by response to the matrix

Habitat loss and fragmentation are major threats to biodiversity and ecosystem processes. Our current understanding of the impacts of habitat loss and fragmentation is based largely on studies that focus on either short-term or long-term responses. Short-term responses are often used to predict long-term responses and make management decisions. The lack of studies comparing short- and long-term responses to fragmentation means we do not adequately understand when and how well short-term responses can be extrapolated to predict long-term responses, and when or why they cannot. To address this gap, we used data from one of the world's longest-running fragmentation experiments, The Wog Wog Habitat Fragmentation Experiment. Using data for carabid beetles, we found that responses in the long term (more than 22 yr post-fragmentation ≈22 generations) often contrasted markedly with those in the short term (5 yr post-fragmentation). The total abundance of all carabids, species richness and the occurrence of six species declined in the short term in the fragments but increased over the long term. The occurrence of three species declined initially and continued to decline, whilst another species was positively affected initially but decreased in the long term. Species' responses to the matrix that surrounds the fragments strongly predicted both the direction (increase/decline in occurrence) and magnitude of their responses to fragmentation. Additionally, species' responses to the matrix were somewhat predicted by their preferences for different types of native habitat (open vs. shaded). Our study highlights the degree of the matrix's influence in fragmented landscapes, and how this influence can change over time. We urge caution in using short-term responses to forecast long-term responses in cases where the matrix (1) impacts species' responses to fragmentation (by isolating them, creating new habitat or altering fragment habitat) and (2) is likely to change through time.

A morphological and phylogenetic investigation into divergence among sympatric Australian southern bull kelps (Durvillaea potatorum and D. amatheiae sp. nov.)

Genetic analyses can reveal a wealth of hitherto undiscovered cryptic biodiversity. For co-occurring and morphologically similar species, the combination of molecular, ecological and morphological analyses provides an excellent opportunity for understanding some of the processes that can lead to divergence and speciation. The Australian endemic brown macroalga Durvillaea potatorum (Phaeophyceae) was examined with a combination of genetic and morphological approaches to confirm the presence of two separate species and to infer the processes that led to their divergence. A total of 331 individuals from 11 sites around coastal Tasmania were collected and measured in situ for a range of morphological and ecological characteristics. Tissue samples were also collected for each individual to allow genetic analyses using mitochondrial (COI) and nuclear (28S) markers. Genetic analyses confirmed the presence of two deeply divergent clades. The significant morphological differentiation, despite high levels of intra-lineage variability, further supported their recognition as distinct species. We describe a new species, D. amatheiae sp. nov., which is characterised by a narrower and proportionately shorter stipe, shorter total length, and higher number of stipitate lateral blades and branches than D. potatorum (sensu stricto). The occurrence of both species in sympatry along Tasmania's eastern and western coasts, as well as their contrasting patterns of haplotype diversity, supports a hypothesis of geographical isolation, allopatric speciation and subsequent secondary contact in response to sea level and ocean current change throughout the Pleistocene glaciation cycles. This research contributes to resolving the phylogenetic relationships, taxonomy and evolution of the ecologically keystone kelp genus Durvillaea.

Dispersal responses override density effects on genetic diversity during post-disturbance succession

Dispersal fundamentally influences spatial population dynamics but little is known about dispersal variation in landscapes where spatial heterogeneity is generated predominantly by disturbance and succession. We tested the hypothesis that habitat succession following fire inhibits dispersal, leading to declines over time in genetic diversity in the early successional gecko Nephrurus stellatus We combined a landscape genetics field study with a spatially explicit simulation experiment to determine whether successional patterns in genetic diversity were driven by habitat-mediated dispersal or demographic effects (declines in population density leading to genetic drift). Initial increases in genetic structure following fire were likely driven by direct mortality and rapid population expansion. Subsequent habitat succession increased resistance to gene flow and decreased dispersal and genetic diversity inN. stellatus Simulated changes in population density alone did not reproduce these results. Habitat-mediated reductions in dispersal, combined with changes in population density, were essential to drive the field-observed patterns. Our study provides a framework for combining demographic, movement and genetic data with simulations to discover the relative influence of demography and dispersal on patterns of landscape genetic structure. Our results suggest that succession can inhibit connectivity among individuals, opening new avenues for understanding how disturbance regimes influence spatial population dynamics.

Do temporal changes in vegetation structure additional to time since fire predict changes in bird occurrence?

Fire is a major ecological process in ecosystems globally. Its impacts on fauna can be both direct (e.g., mortality) and indirect (e.g., altered habitat), resulting in population recovery being driven by several possible mechanisms. Separating direct from indirect impacts of fire on faunal population recovery can be valuable in guiding management of biodiversity in fire-prone environments. However, resolving the influence of direct and indirect processes remains a key challenge because many processes affecting fauna can change concomitantly with time since fire. We explore the mechanisms influencing bird response to fire by posing the question, can temporal changes in vegetation structure predict changes in bird occurrence on sites, and can these be separated from other temporal changes using the surrogate of time since fire? We conducted a 12-yr study of bird and vegetation responses to fire at 124 sites across six vegetation classes in Booderee National Park, Australia. Approximately half of these sites, established in 2002, were burned by a large (>3000 ha) wildfire in 2003. To disentangle collinear effects of temporal changes in vegetation and direct demographic effects on population recovery that are subsumed by time since fire, we incorporated both longitudinal and cross-sectional vegetation effects in addition to time since fire within logistic structural equation models. We identified temporal changes in vegetation structure and richness of plant and bird species that characterized burned and unburned sites in all vegetation classes. For nine bird species, a significant component of the year trend was driven by temporal trends in one of three vegetation variables (number of understory or midstory plant species, or midstory cover). By contrast, we could not separate temporal effects between time since fire and vegetation attributes for bird species richness, reporting rate, and the occurrence of 11 other bird species. Our findings help identify species for which indirect effects of vegetation dominate recovery and thus may benefit from vegetation management where conservation actions are required and, conversely, those species for which direct effects of time since fire drive recovery, where simply leaving a system to recover following the last disturbance will be sufficient.

Disturbance gradient shows logging affects plant functional groups more than fire

Understanding the impacts of natural and human disturbances on forest biota is critical for improving forest management. Many studies have examined the separate impacts on fauna and flora of wildfire, conventional logging, and salvage logging, but empirical comparisons across a broad gradient of simultaneous disturbances are lacking. We quantified species richness and frequency of occurrence of vascular plants, and functional group responses, across a gradient of disturbances that occurred concurrently in 2009 in the mountain ash forests of southeastern Australia. Our study encompassed replicated sites in undisturbed forest (~70 yr post fire), forest burned at low severity, forest burned at high severity, unburned forest that was clearcut logged, and forest burned at high severity that was clearcut salvage logged post-fire. All sites were sampled 2 and 3 yr post fire. Mean species richness decreased across the disturbance gradient from 30.1 species/site on low-severity burned sites and 28.9 species/site on high-severity burned sites, to 25.1 species/site on clearcut sites and 21.7 species/site on salvage logged sites. Low-severity burned sites were significantly more species-rich than clearcut sites and salvage logged sites; high-severity burned sites supported greater species richness than salvage logged sites. Specific traits influenced species' sensitivity to disturbance. Resprouting species dominated undisturbed mountain ash forests, but declined significantly across the gradient. Fern and midstory trees decreased significantly in frequency of occurrence across the gradient. Ferns (excluding bracken) decreased from 34% of plants in undisturbed forest to 3% on salvage logged sites. High-severity burned sites supported a greater frequency of occurrence and species richness of midstory trees compared to clearcut and salvage logged sites. Salvage logging supported fewer midstory trees than any other disturbance category, and were distinctly different from clearcut sites. Plant life form groups, including midstory trees, shrubs, and ferns, were dominated by very few species on logged sites. The differences in biotic response across the gradient of natural and human disturbances have significant management implications, particularly the need to reduce mechanical disturbance overall and to leave specific areas with no mechanical disturbance across the cut area during logging operations, to ensure the persistence of resprouting taxa.

Temporal trends in mammal responses to fire reveals the complex effects of fire regime attributes

Fire is a major ecological process in many ecosystems worldwide. We sought to identify which attributes of fire regimes affect temporal change in the presence and abundance of Australian native mammals. Our detailed study was underpinned by time series data on 11 mammal species at 97 long-term sites in southeastern Australia between 2003 and 2013. We explored how temporal aspects of fire regimes influenced the presence and conditional abundance of species. The key fire regime components examined were: (1) severity of a major fire in 2003, (2) interval between the last major fire (2003) and the fire prior to that, and (3) number of past fires. Our long-term data set enabled quantification of the interactions between survey year and each fire regime variable: an ecological relationship missing from temporally restricted studies. We found no evidence of any appreciable departures from the assumption of independence of the sites. Multiple aspects of fire regimes influenced temporal variation in the presence and abundance of mammals. The best models indicated that six of the 11 species responded to two or more fire regime variables, with two species influenced by all three fire regime attributes. Almost all species responded to time since fire, either as an interaction with survey year or as a main effect. Fire severity or its interaction with survey year was important for most terrestrial rodents. The number of fires at a site was significant for terrestrial rodents and several other species. Our findings contain evidence of the effects on native mammals of heterogeneity in fire regimes. Temporal response patterns of mammal species were influenced by multiple fire regime attributes, often in conjunction with survey year. This underscores the critical importance of long-term studies of biota that are coupled with data sets characterized by carefully documented fire history, severity, and frequency. Long-term studies are essential to predict animal responses to fires and guide management of when and where (prescribed) fire or, conversely, long-unburned vegetation is needed. The complexity of observed responses highlights the need for large reserves in which patterns of heterogeneity in fire regimes can be sustained in space and over time.

Implications of recurrent disturbance for genetic diversity

Exploring interactions between ecological disturbance, species’ abundances and community composition provides critical insights for ecological dynamics. While disturbance is also potentially an important driver of landscape genetic patterns, the mechanisms by which these patterns may arise by selective and neutral processes are not well‐understood. We used simulation to evaluate the relative importance of disturbance regime components, and their interaction with demographic and dispersal processes, on the distribution of genetic diversity across landscapes. We investigated genetic impacts of variation in key components of disturbance regimes and spatial patterns that are likely to respond to climate change and land management, including disturbance size, frequency, and severity. The influence of disturbance was mediated by dispersal distance and, to a limited extent, by birth rate. Nevertheless, all three disturbance regime components strongly influenced spatial and temporal patterns of genetic diversity within subpopulations, and were associated with changes in genetic structure. Furthermore, disturbance‐induced changes in temporal population dynamics and the spatial distribution of populations across the landscape resulted in disrupted isolation by distance patterns among populations. Our results show that forecast changes in disturbance regimes have the potential to cause major changes to the distribution of genetic diversity within and among populations. We highlight likely scenarios under which future changes to disturbance size, severity, or frequency will have the strongest impacts on population genetic patterns. In addition, our results have implications for the inference of biological processes from genetic data, because the effects of dispersal on genetic patterns were strongly mediated by disturbance regimes.

Identifying the location of fire refuges in wet forest ecosystems

The increasing frequency of large, high-severity fires threatens the survival of old-growth specialist fauna in fire-prone forests. Within topographically diverse montane forests, areas that experience less severe or fewer fires compared with those prevailing in the landscape may present unique resource opportunities enabling old-growth specialist fauna to survive. Statistical landscape models that identify the extent and distribution of potential fire refuges may assist land managers to incorporate these areas into relevant biodiversity conservation strategies. We used a case study in an Australian wet montane forest to establish how predictive fire simulation models can be interpreted as management tools to identify potential fire refuges. We examined the relationship between the probability of fire refuge occurrence as predicted by an existing fire refuge model and fire severity experienced during a large wildfire. We also examined the extent to which local fire severity was influenced by fire severity in the surrounding landscape. We used a combination of statistical approaches, including generalized linear modeling, variogram analysis, and receiver operating characteristics and area under the curve analysis (ROC AUC). We found that the amount of unburned habitat and the factors influencing the retention and location of fire refuges varied with fire conditions. Under extreme fire conditions, the distribution of fire refuges was limited to only extremely sheltered, fire-resistant regions of the landscape. During extreme fire conditions, fire severity patterns were largely determined by stochastic factors that could not be predicted by the model. When fire conditions were moderate, physical landscape properties appeared to mediate fire severity distribution. Our study demonstrates that land managers can employ predictive landscape fire models to identify the broader climatic and spatial domain within which fire refuges are likely to be present. It is essential that within these envelopes, forest is protected from logging, roads, and other developments so that the ecological processes related to the establishment and subsequent use of fire refuges are maintained.

Fine-scale refuges can buffer demographic and genetic processes against short-term climatic variation and disturbance: a 22-year case study of an arboreal marsupial

Ecological disturbance and climate are key drivers of temporal dynamics in the demography and genetic diversity of natural populations. Microscale refuges are known to buffer species' persistence against environmental change, but the effects of such refuges on demographic and genetic patterns in response to short-term environmental variation are poorly understood. We quantified demographic and genetic responses of mountain brushtail possums (Trichosurus cunninghami) to rainfall variability (1992-2013) and to a major wildfire. We hypothesized that there would be underlying differences in demographic and genetic processes between an unburnt mesic refuge and a topographically exposed zone that was burnt in 2009. Fire caused a 2-year decrease in survival in the burnt zone, but the population grew after the fire due to immigration, leading to increased expected heterozygosity. We documented a fire-related behavioural shift, where the rate of movement by individuals in the unburnt refuge to the burnt zone decreased after fire. Irrespective of the fire, there were long-term differences in demographic and genetic parameters between the mesic/unburnt refuge and the nonmesic/burnt zone. Survival was high and unaffected by rainfall in the refuge, but lower and rainfall-dependent in the nonmesic zone. Net movement of individuals was directional, from the mesic refuge to the nonmesic zone, suggesting fine-scale source-sink dynamics. There were higher expected heterozygosity (HE ) and temporal genetic stability in the refuge, but lower HE and marked temporal genetic structure in the exposed habitat, consistent with reduced generational overlap caused by elevated mortality and immigration. Thus, fine-scale refuges can mediate the short-term demographic and genetic effects of climate and ecological disturbance.

Phylogeographic Structure in Penguin Ticks across an Ocean Basin Indicates Allopatric Divergence and Rare Trans-Oceanic Dispersal

The association of ticks (Acarina) and seabirds provides an intriguing system for assessing the influence of long-distance dispersal on the evolution of parasitic species. Recent research has focused on host-parasite evolutionary relationships and dispersal capacity of ticks parasitising flighted seabirds. Evolutionary research on the ticks of non-flighted seabirds is, in contrast, scarce. We conducted the first phylogeographic investigation of a hard tick species (Ixodes eudyptidis) that parasitises the Little Blue Penguin (Eudyptula minor). Using one nuclear (28S) and two mitochondrial (COI and 16S) markers, we assessed genetic diversity among several populations in Australia and a single population on the South Island of New Zealand. Our results reveal two deeply divergent lineages, possibly representing different species: one comprising all New Zealand samples and some from Australia, and the other representing all other samples from Australian sites. No significant population differentiation was observed among any Australian sites from within each major clade, even those separated by hundreds of kilometres of coastline. In contrast, the New Zealand population was significantly different to all samples from Australia. Our phylogenetic results suggest that the New Zealand and Australian populations are effectively isolated from each other; although rare long-distance dispersal events must occur, these are insufficient to maintain trans-Tasman gene flow. Despite the evidence for limited dispersal of penguin ticks between Australia and New Zealand, we found no evidence to suggest that ticks are unable to disperse shorter distances at sea with their hosts, with no pattern of population differentiation found among Australian sites. Our results suggest that terrestrial seabird parasites may be quite capable of short-distance movements, but only sporadic longer-distance (trans-oceanic) dispersal.

The effect of sex-biased dispersal on opposite-sexed spatial genetic structure and inbreeding risk

Natal sex-biased dispersal has long been thought to reduce the risk of inbreeding by spatially separating opposite-sexed kin. Yet, comprehensive and quantitative evaluations of this hypothesis are lacking. In this study, we quantified the effectiveness of sex-biased dispersal as an inbreeding avoidance strategy by combining spatially explicit simulations and empirical data. We quantified the extent of kin clustering by measuring the degree of spatial autocorrelation among opposite-sexed individuals (FM structure). This allowed us to systematically explore how the extent of sex-biased dispersal, generational overlap, and mate searching distance, influenced both kin clustering, and the resulting inbreeding in the absence of complementary inbreeding avoidance strategies. Simulations revealed that when sex-biased dispersal was limited, positive FM genetic structure developed quickly and increased as the mate searching distance decreased or as generational overlap increased. Interestingly, complete long-range sex-biased dispersal did not prevent the development of FM genetic structure when generations overlapped. We found a very strong correlation between FM genetic structure and both FIS under random mating, and pedigree-based measures of inbreeding. Thus, we show that the detection of FM genetic structure can be a strong indicator of inbreeding risk. Empirical data for two species with different life history strategies yielded patterns congruent with our simulations. Our study illustrates a new application of spatial genetic autocorrelation analysis that offers a framework for quantifying the risk of inbreeding that is easily extendable to other species. Furthermore, our findings provide other researchers with a context for interpreting observed patterns of opposite-sexed spatial genetic structure.

The need for a comprehensive reassessment of the Regional Forest Agreements in Australia

Regional Forest Agreements (RFAs) are State–Federal agreements underpinning the management of the majority of Australia’s commercially productive native forests. Introduced between 1997 and 2001, they were designed to deliver certainty to forest industries while, simultaneously, guaranteeing environmental protection, including the conservation of biodiversity. Using examples, we argue that RFAs in some jurisdictions have failed to do either. We strongly recommend a comprehensive reassessment of RFAs. This is needed to: (1) take into account significant new knowledge on forest ecology and management that has been gathered in the past 20 years, including updated prognoses for some critically endangered species; (2) better evaluate the full range of wood and non-wood products and services provided by forests; (3) accommodate new methods of forest inventory and more environmentally sensitive silvicultural systems; and (4) better account for the impacts of natural disturbances, such as fires, on the area available for logging, sustained yield, and forest ecosystem integrity per se. Without a substantial overhaul of the RFAs, there is a significant risk of undervaluing the full range of native forest values, exacerbating species declines, and permanently damaging forest ecosystems.

Geographical variation in body size and sexual size dimorphism in an Australian lizard, Boulenger's Skink (Morethia boulengeri)

Ecogeographical rules help explain spatial and temporal patterns in intraspecific body size. However, many of these rules, when applied to ectothermic organisms such as reptiles, are controversial and require further investigation. To explore factors that influence body size in reptiles, we performed a heuristic study to examine body size variation in an Australian lizard, Boulenger's Skink Morethia boulengeri from agricultural landscapes in southern New South Wales, south-eastern Australia. We collected tissue and morphological data on 337 adult lizards across a broad elevation and climate gradient. We used a model-selection procedure to determine if environmental or ecological variables best explained body size variation. We explored the relationship between morphology and phylogenetic structure before modeling candidate variables from four broad domains: (1) geography (latitude, longitude and elevation), (2) climate (temperature and rainfall), (3) habitat (vegetation type, number of logs and ground cover attributes), and (4) management (land use and grazing history). Broad phylogenetic structure was evident, but on a scale larger than our study area. Lizards were sexually dimorphic, whereby females had longer snout-vent length than males, providing support for the fecundity selection hypothesis. Body size variation in M. boulengeri was correlated with temperature and rainfall, a pattern consistent with larger individuals occupying cooler and more productive parts of the landscape. Climate change forecasts, which predict warmer temperature and increased aridity, may result in reduced lizard biomass and decoupling of trophic interactions with potential implications for community organization and ecosystem function.