Using a Handful of Transcriptomes to Detect Sex-Linked Markers and Develop Molecular Sexing Assays in a Species with Homomorphic Sex Chromosomes

To understand the biology of a species, it is often crucial to be able to differentiate males and females. However, many species lack easily identifiable sexually dimorphic traits. In those that possess sex chromosomes, molecular sexing offers a good alternative, and molecular sexing assays can be developed through the comparison of male and female genomic sequences. However, in many nonmodel species, sex chromosomes are poorly differentiated, and identifying sex-linked sequences and developing sexing assays can be challenging. In this study, we highlight a simple transcriptome-based procedure for the detection of sex-linked markers suitable for the development of sexing assays that circumvents limitations of more commonly used approaches. We apply it to the spotted snow skink Carinascincus ocellatus, a viviparous lizard with homomorphic XY chromosomes that has environmentally induced sex reversal. With transcriptomes from three males and three females alone, we identify thousands of putative Y-linked sequences. We confirm linkage through alignment of assembled transcripts to a distantly related lizard genome and readily design multiple single locus polymerase chain reaction primers to sex C. ocellatus and related species. Our approach also facilitates valuable comparisons of sex determining systems on a broad taxonomic scale.

Disease-driven top predator decline affects mesopredator population genomic structure

Top predator declines are pervasive and often have dramatic effects on ecological communities via changes in food web dynamics, but their evolutionary consequences are virtually unknown. Tasmania's top terrestrial predator, the Tasmanian devil, is declining due to a lethal transmissible cancer. Spotted-tailed quolls benefit via mesopredator release, and they alter their behaviour and resource use concomitant with devil declines and increased disease duration. Here, using a landscape community genomics framework to identify environmental drivers of population genomic structure and signatures of selection, we show that these biotic factors are consistently among the top variables explaining genomic structure of the quoll. Landscape resistance negatively correlates with devil density, suggesting that devil declines will increase quoll genetic subdivision over time, despite no change in quoll densities detected by camera trap studies. Devil density also contributes to signatures of selection in the quoll genome, including genes associated with muscle development and locomotion. Our results provide some of the first evidence of the evolutionary impacts of competition between a top predator and a mesopredator species in the context of a trophic cascade. As top predator declines are increasing globally, our framework can serve as a model for future studies of evolutionary impacts of altered ecological interactions.

From rivers to ocean basins: The role of ocean barriers and philopatry in the genetic structuring of a cosmopolitan coastal predator

The Bull Shark (Carcharhinus leucas) faces varying levels of exploitation around the world due to its coastal distribution. Information regarding population connectivity is crucial to evaluate its conservation status and local fishing impacts. In this study, we sampled 922 putative Bull Sharks from 19 locations in the first global assessment of population structure of this cosmopolitan species. Using a recently developed DNA-capture approach (DArTcap), samples were genotyped for 3400 nuclear markers. Additionally, full mitochondrial genomes of 384 Indo-Pacific samples were sequenced. Reproductive isolation was found between and across ocean basins (eastern Pacific, western Atlantic, eastern Atlantic, Indo-West Pacific) with distinct island populations in Japan and Fiji. Bull Sharks appear to maintain gene flow using shallow coastal waters as dispersal corridors, whereas large oceanic distances and historical land-bridges act as barriers. Females tend to return to the same area for reproduction, making them more susceptible to local threats and an important focus for management actions. Given these behaviors, the exploitation of Bull Sharks from insular populations, such as Japan and Fiji, may instigate local decline that cannot readily be replenished by immigration, which can in turn affect ecosystem dynamics and functions. These data also supported the development of a genetic panel to ascertain the population of origin, which will be useful in monitoring the trade of fisheries products and assessing population-level impacts of this harvest.

Habitat connectivity and host relatedness influence virus spread across an urbanising landscape in a fragmentation-sensitive carnivore

Spatially heterogeneous landscape factors such as urbanisation can have substantial effects on the severity and spread of wildlife diseases. However, research linking patterns of pathogen transmission to landscape features remains rare. Using a combination of phylogeographic and machine learning approaches, we tested the influence of landscape and host factors on feline immunodeficiency virus (FIV_Lru) genetic variation and spread among bobcats (Lynx rufus) sampled from coastal southern California. We found evidence for increased rates of FIV_Lru lineage spread through areas of higher vegetation density. Furthermore, single-nucleotide polymorphism (SNP) variation among FIV_Lru sequences was associated with host genetic distances and geographic location, with FIV_Lru genetic discontinuities precisely correlating with known urban barriers to host dispersal. An effect of forest land cover on FIV_Lru SNP variation was likely attributable to host population structure and differences in forest land cover between different populations. Taken together, these results suggest that the spread of FIV_Lru is constrained by large-scale urban barriers to host movement. Although urbanisation at fine spatial scales did not appear to directly influence virus transmission or spread, we found evidence that viruses transmit and spread more quickly through areas containing higher proportions of natural habitat. These multiple lines of evidence demonstrate how urbanisation can change patterns of contact-dependent pathogen transmission and provide insights into how continued urban development may influence the incidence and management of wildlife disease.

Population genomics of a predatory mammal reveals patterns of decline and impacts of exposure to toxic toads

Mammal declines across northern Australia are one of the major biodiversity loss events occurring globally. There has been no regional assessment of the implications of these species declines for genomic diversity. To address this, we conducted a species-wide assessment of genomic diversity in the northern quoll (Dasyurus hallucatus), an Endangered marsupial carnivore. We used next generation sequencing methods to genotype 10,191 single nucleotide polymorphisms (SNPs) in 352 individuals from across a 3220-km length of the continent, investigating patterns of population genomic structure and diversity, and identifying loci showing signals of putative selection. We found strong heterogeneity in the distribution of genomic diversity across the continent, characterized by (i) biogeographical barriers driving hierarchical population structure through long-term isolation, and (ii) severe reductions in diversity resulting from population declines, exacerbated by the spread of introduced toxic cane toads (Rhinella marina). These results warn of a large ongoing loss of genomic diversity and associated adaptive capacity as mammals decline across northern Australia. Encouragingly, populations of the northern quoll established on toad-free islands by translocations appear to have maintained most of the initial genomic diversity after 16 years. By mapping patterns of genomic diversity within and among populations, and investigating these patterns in the context of population declines, we can provide conservation managers with data critical to informed decision-making. This includes the identification of populations that are candidates for genetic management, the importance of remnant island and insurance/translocated populations for the conservation of genetic diversity, and the characterization of putative evolutionarily significant units.

Sarcoptic mange changes bacterial and fungal microbiota of bare-nosed wombats (Vombatus ursinus)

Background

Sarcoptes scabiei is globally distributed and one of the most impactful mammalian ectoparasites. Sarcoptic mange, caused by infection with S. scabiei, causes disruption of the epidermis and its bacterial microbiota, but its effects on host fungal microbiota and on the microbiota of marsupials in general have not been studied. Here, we (i) examine bacterial and fungal microbiota changes associated with mange in wild bare-nosed wombats (BNWs) and (ii) evaluate whether opportunistic pathogens are potentiated by S. scabiei infection in this species.

Methods

Using Amplicon Sequencing of the 16S rRNA and ITS2 rDNA genes, we detected skin microbiota changes of the bare-nosed wombat (Vombatus ursinus). We compared the alpha and beta diversity among healthy, moderate, and severe disease states using ANOVA and PERMANOVA with nesting. Lastly, we identified taxa that differed between disease states using analysis of composition of microbes (ANCOM) testing.

Results

We detected significant changes in the microbial communities and diversity with mange in BNWs. Severely affected BNWs had lower amplicon sequence variant (ASV) richness compared to that of healthy individuals, and the microbial communities were significantly different between disease states with higher relative abundance of potentially pathogenic microbial taxa in mange-affected BNWs including Staphylococcus sciuri, Corynebacterium spp., Brevibacterium spp., Brachybacterium spp., and Pseudogymnascus spp. and Debaryomyces spp.

Conclusion

This study represents the first investigation of microbial changes in association with sarcoptic mange in a marsupial host, as well as the first investigation of fungal microbial changes on the skin of any host suffering from sarcoptic mange. Our results are broadly consistent with bacterial microbiota changes observed in humans, pigs, canids, and Iberian ibex, suggesting the epidermal microbial impacts of mange may be generalisable across host species. We recommend that future studies investigating skin microbiota changes include both bacterial and fungal data to gain a more complete picture of the effects of sarcoptic mange.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05452-y.

Sex reversal explains some, but not all, climate-mediated sex ratio variation within a viviparous reptile

Evolutionary transitions in sex-determining systems have occurred frequently yet understanding how they occur remains a major challenge. In reptiles, transitions from genetic to temperature-dependent sex determination can occur if the gene products that determine sex evolve thermal sensitivity, resulting in sex-reversed individuals. However, evidence of sex reversal is limited to oviparous reptiles. Here we used thermal experiments to test whether sex reversal is responsible for differences in sex determination in a viviparous reptile, Carinascincus ocellatus, a species with XY sex chromosomes and population-specific sex ratio response to temperature. We show that sex reversal is occurring and that its frequency is related to temperature. Sex reversal was unidirectional (phenotypic males with XX genotype) and observed in both high- and low-elevation populations. We propose that XX-biased genotypic sex ratios could produce either male- or female-biased phenotypic sex ratios as observed in low-elevation C. ocellatus under variable rates of XX sex reversal. We discuss reasons why sex reversal may not influence sex ratios at high elevation. Our results suggest that the mechanism responsible for evolutionary transitions from genotypic to temperature-dependent sex determination is more complex than can be explained by a single process such as sex reversal.

Climate Shapes the Geographic Distribution and Introgressive Spread of Color Ornamentation in Common Wall Lizards

AbstractClimate can exert an effect on the strength of sexual selection, but empirical evidence is limited. Here, we tested whether climate predicts the geographic distribution and introgressive spread of sexually selected male color ornamentation across 114 populations of the common wall lizard, Podarcis muralis. Coloration was highly structured across the landscape and did not reflect genetic differentiation. Instead, color ornamentation was consistently exaggerated in hot and dry environments, suggesting that climate-driven selection maintains geographic variation in spite of gene flow. Introgression of color ornamentation into a distantly related lineage appears to be ongoing and was particularly pronounced in warm climates with wet winters and dry summers. Combined, these results suggest that sexual ornamentation is consistently favored in climates that allow a prolonged reproductive season and high and reliable opportunities for lizard activity. This pattern corroborates theoretical predictions that such climatic conditions reduce the temporal clustering of receptive females and increase male-male competition, resulting in strong sexual selection. In summary, we provide compelling evidence for the importance of climate for the evolution of color ornamentation, and we demonstrate that geographic variation in the strength of sexual selection influences introgression of this phenotype.

Lack of detectable genetic isolation in the cyclic rodent Microtus arvalis despite large landscape fragmentation owing to transportation infrastructures

Although roads are widely seen as dispersal barriers, their genetic consequences for animals that experience large fluctuations in population density are poorly documented. We developed a spatially paired experimental design to assess the genetic impacts of roads on cyclic voles (Microtus arvalis) during a high-density phase in North-Western Spain. We compared genetic patterns from 15 paired plots bisected by three different barrier types, using linear mixed models and computing effect sizes to assess the importance of each type, and the influence of road features like width or the age of the infrastructure. Evidence of effects by roads on genetic diversity and differentiation were lacking. We speculate that the recurrent (each 3-5 generations) episodes of massive dispersal associated with population density peaks can homogenize populations and mitigate the possible genetic impact of landscape fragmentation by roads. This study highlights the importance of developing spatially replicated experimental designs that allow us to consider the large natural spatial variation in genetic parameters. More generally, these results contribute to our understanding of the not well explored effects of habitat fragmentation on dispersal in species showing "boom-bust" dynamics.

Contrasting population manipulations reveal resource competition between two large marsupials: bare-nosed wombats and eastern grey kangaroos

Resource competition is an important interaction that can structure ecological communities, but is difficult to demonstrate in nature, and rarely demonstrated for large mammals including marsupials. We analysed 10 years of population survey data to investigate resource competition between bare-nosed wombats (Vombatus ursinus) and eastern grey kangaroos (Macropus giganteus) at two sites to assess whether resource competition is occurring. At one site, wombat abundance was reduced by increased mortality from mange disease, whereas at the other site, kangaroo abundance was reduced primarily by culling. We used the modified Lotka-Volterra competition (LVC) models to describe the mechanism of resource competition and fitted those models to the empirical data by maximum likelihood estimation. We found strong negative relationships between the abundance of wombats and kangaroos at each site, and resource competition was also mechanistically supported by the modified LVC models. The estimated competition coefficients indicate that bare-nosed wombats are a slightly superior competitor of eastern grey kangaroos than vice versa, and that intraspecific competition is almost twice as strong as interspecific competition. In addition, this study facilitated the calculation of the transmission rate associated with mange disease at one site (0.011), and the removal rate owing to culling, the introduction of a predator species, and drought at the other site (0.0006). Collectively, this research represents a rare empirical demonstration of resource competition between large mammals and contributes new insight into the ecology of two of Australia's largest grazing marsupials.

Pleistocene divergence in the absence of gene flow among populations of a viviparous reptile with intraspecific variation in sex determination

Polymorphisms can lead to genetic isolation if there is differential mating success among conspecifics divergent for a trait. Polymorphism for sex-determining system may fall into this category, given strong selection for the production of viable males and females and the low success of heterogametic hybrids when sex chromosomes differ (Haldane's rule). Here we investigated whether populations exhibiting polymorphism for sex determination are genetically isolated, using the viviparous snow skink Carinascincus ocellatus. While a comparatively high elevation population has genotypic sex determination, in a lower elevation population there is an additional temperature component to sex determination. Based on 11,107 SNP markers, these populations appear genetically isolated. "Isolation with Migration" analysis also suggests these populations diverged in the absence of gene flow, across a period encompassing multiple Pleistocene glaciations and likely greater geographic proximity of populations. However, further experiments are required to establish whether genetic isolation may be a cause or consequence of differences in sex determination. Given the influence of temperature on sex in one lineage, we also discuss the implications for the persistence of this polymorphism under climate change.

The effects of weather variability on patterns of genetic diversity in Tasmanian bettongs

While the effects of climate (long-term, prevailing weather) on species abundance, range and genetic diversity have been widely studied, short-term, localized variations in atmospheric conditions (i.e., weather) can also rapidly alter species' geographical ranges and population sizes, but little is known about how they affect genetic diversity. We investigated the relationship between weather and range-wide genetic diversity in a marsupial, Bettongia gaimardi, using dynamic species distribution models (SDMs). Genetic diversity was lower in parts of the range where the weather-based SDM predicted high variability in probability of B. gaimardi occurrence during 1950-2009. This is probably an effect of lower population sizes and extinction-recolonization cycles in places with highly variable weather. Spatial variation in genetic diversity was also better predicted by mean probabilities of B. gaimardi occurrence from weather- than climate-based SDMs. Our results illustrate the importance of weather in driving population dynamics and species distributions on decadal timescales and thereby in affecting genetic diversity. Modelling the links between changing weather patterns, species distributions and genetic diversity will allow researchers to better forecast biological impacts of climate change.

Morphometric and population genomic evidence for species divergence in the Chimarrichthys fish complex of the Tibetan Plateau

The uplift of the Tibetan Plateau altered the environmental conditions of the local area substantially. Here, we conducted a comprehensive investigation based on morphometrics, population genomics, and climatic factors to evaluate phenotypic and genome-level variations in a radiation of Chimarrichthys catfish endemic to the Plateau. Discriminant function analysis showed phenotypic differences of Chimarrichthys between rivers with respect to elevation. Genetic structure analysis based on 6606 single nucleotide polymorphisms (SNPs) deduced genetic differences between rivers, and species delimitation indicated that the Chimarrichthys fish complex could be divided into three species. Restriction site-associated DNA tags were mapped to the gene sets of Glyptosternon maculatum, and matches were searched against databases for Gene Ontology annotation. Genomic regions exhibiting marked differences among localities represented a range of biological functions, including growth (gdf11), bone development (bmp8a), cellular response to light stimulus (opn3), regulation of the rhodopsin-mediated signalling pathway (grk1), immune response (rag1 and ung), reproductive process (antxr2), and regulation of intracellular iron levels (ireb2). The tag44126, where gene gdf11 is located, was identified as an outlier exhibiting divergence between rivers with altitude differences, and the SNP is thymine (T) in Dadu and Yalong River (~2700 m), but guanine (G) in Jinsha and Qingyi rivers (~2200 and ~ 684 m), suggesting a possible effect of altitude on its differentiation.

Differences in Homomorphic Sex Chromosomes Are Associated with Population Divergence in Sex Determination in Carinascincus ocellatus (Scincidae: Lygosominae)

Sex determination directs development as male or female in sexually reproducing organisms. Evolutionary transitions in sex determination have occurred frequently, suggesting simple mechanisms behind the transitions, yet their detail remains elusive. Here we explore the links between mechanisms of transitions in sex determination and sex chromosome evolution at both recent and deeper temporal scales (<1 Myr; ~79 Myr). We studied a rare example of a species with intraspecific variation in sex determination, Carinascincus ocellatus, and a relative, Liopholis whitii, using c-banding and mapping of repeat motifs and a custom Y chromosome probe set to identify the sex chromosomes. We identified both unique and conserved regions of the Y chromosome among C. ocellatus populations differing in sex determination. There was no evidence for homology of sex chromosomes between C. ocellatus and L. whitii, suggesting independent evolutionary origins. We discuss sex chromosome homology between members of the subfamily Lygosominae and propose links between sex chromosome evolution, sex determination transitions, and karyotype evolution.

Does the virus cross the road? Viral phylogeographic patterns among bobcat populations reflect a history of urban development

Urban development has major impacts on connectivity among wildlife populations and is thus likely an important factor shaping pathogen transmission in wildlife. However, most investigations of wildlife diseases in urban areas focus on prevalence and infection risk rather than potential effects of urbanization on transmission itself. Feline immunodeficiency virus (FIV) is a directly transmitted retrovirus that infects many felid species and can be used as a model for studying pathogen transmission at landscape scales. We investigated phylogenetic relationships among FIV isolates sampled from five bobcat (Lynx rufus) populations in coastal southern California that appear isolated due to major highways and dense urban development. Divergence dates among FIV phylogenetic lineages in several cases reflected historical urban growth and construction of major highways. We found strong FIV phylogeographic structure among three host populations north-west of Los Angeles, largely coincident with host genetic structure. In contrast, relatively little FIV phylogeographic structure existed among two genetically distinct host populations south-east of Los Angeles. Rates of FIV transfer among host populations did not vary significantly, with the lack of phylogenetic structure south-east of Los Angeles unlikely to reflect frequent contemporary transmission among populations. Our results indicate that major barriers to host gene flow can also act as barriers to pathogen spread, suggesting potentially reduced susceptibility of fragmented populations to novel directly transmitted pathogens. Infrequent exchange of FIV among host populations suggests that populations would best be managed as distinct units in the event of a severe disease outbreak. Phylogeographic inference of pathogen transmission is useful for estimating the ability of geographic barriers to constrain disease spread and can provide insights into contemporary and historical drivers of host population connectivity.

Identifying mechanisms of genetic differentiation among populations in vagile species: historical factors dominate genetic differentiation in seabirds

Elucidating the factors underlying the origin and maintenance of genetic variation among populations is crucial for our understanding of their ecology and evolution, and also to help identify conservation priorities. While intrinsic movement has been hypothesized as the major determinant of population genetic structuring in abundant vagile species, growing evidence indicates that vagility does not always predict genetic differentiation. However, identifying the determinants of genetic structuring can be challenging, and these are largely unknown for most vagile species. Although, in principle, levels of gene flow can be inferred from neutral allele frequency divergence among populations, underlying assumptions may be unrealistic. Moreover, molecular studies have suggested that contemporary gene flow has often not overridden historical influences on population genetic structure, which indicates potential inadequacies of any interpretations that fail to consider the influence of history in shaping that structure. This exhaustive review of the theoretical and empirical literature investigates the determinants of population genetic differentiation using seabirds as a model system for vagile taxa. Seabirds provide a tractable group within which to identify the determinants of genetic differentiation, given their widespread distribution in marine habitats and an abundance of ecological and genetic studies conducted on this group. Herein we evaluate mitochondrial DNA (mtDNA) variation in 73 seabird species. Lack of mutation-drift equilibrium observed in 19% of species coincided with lower estimates of genetic differentiation, suggesting that dynamic demographic histories can often lead to erroneous interpretations of contemporary gene flow, even in vagile species. Presence of land across the species sampling range, or sampling of breeding colonies representing ice-free Pleistocene refuge zones, appear to be associated with genetic differentiation in Tropical and Southern Temperate species, respectively, indicating that long-term barriers and persistence of populations are important for their genetic structuring. Conversely, biotic factors commonly considered to influence population genetic structure, such as spatial segregation during foraging, were inconsistently associated with population genetic differentiation. In light of these results, we recommend that genetic studies should consider potential historical events when identifying determinants of genetic differentiation among populations to avoid overestimating the role of contemporary factors, even for highly vagile taxa.

Australian lizards are outstanding models for reproductive biology research

Australian lizards are a diverse group distributed across the continent and inhabiting a wide range of environments. Together, they exhibit a remarkable diversity of reproductive morphologies, physiologies, and behaviours that is broadly representative of vertebrates in general. Many reproductive traits exhibited by Australian lizards have evolved independently in multiple lizard lineages, including sociality, complex signalling and mating systems, viviparity, and temperature-dependent sex determination. Australian lizards are thus outstanding model organisms for testing hypotheses about how reproductive traits function and evolve, and they provide an important basis of comparison with other animals that exhibit similar traits. We review how research on Australian lizard reproduction has contributed to answering broader evolutionary and ecological questions that apply to animals in general. We focus on reproductive traits, processes, and strategies that are important areas of current research, including behaviours and signalling involved in courtship; mechanisms involved in mating, egg production, and sperm competition; nesting and gestation; sex determination; and finally, birth in viviparous species. We use our review to identify important questions that emerge from an understanding of this body of research when considered holistically. Finally, we identify additional research questions within each topic that Australian lizards are well suited for reproductive biologists to address.

Urbanization reduces genetic connectivity in bobcats (Lynx rufus) at both intra- and interpopulation spatial scales

Urbanization is a major factor driving habitat fragmentation and connectivity loss in wildlife. However, the impacts of urbanization on connectivity can vary among species and even populations due to differences in local landscape characteristics, and our ability to detect these relationships may depend on the spatial scale at which they are measured. Bobcats (Lynx rufus) are relatively sensitive to urbanization and the status of bobcat populations is an important indicator of connectivity in urban coastal southern California. We genotyped 271 bobcats at 13,520 SNP loci to conduct a replicated landscape resistance analysis in five genetically distinct populations. We tested urban and natural factors potentially influencing individual connectivity in each population separately, as well as study-wide. Overall, landscape genomic effects were most frequently detected at the study-wide spatial scale, with urban land cover (measured as impervious surface) having negative effects and topographic roughness having positive effects on gene flow. The negative effect of urban land cover on connectivity was also evident when populations were analyzed separately despite varying substantially in spatial area and the proportion of urban development, confirming a pervasive impact of urbanization largely independent of spatial scale. The effect of urban development was strongest in one population where stream habitat had been lost to development, suggesting that riparian corridors may help mitigate reduced connectivity in urbanizing areas. Our results demonstrate the importance of replicating landscape genetic analyses across populations and considering how landscape genetic effects may vary with spatial scale and local landscape structure.

Isolation, marine transgression and translocation of the bare-nosed wombat (Vombatus ursinus)

Island populations can represent genetically distinct and evolutionarily important lineages relative to mainland conspecifics. However, phenotypic divergence of island populations does not necessarily reflect genetic divergence, particularly for lineages inhabiting islands periodically connected during Pleistocene low sea stands. Marine barriers may also not be solely responsible for any divergence that is observed. Here, we investigated genetic divergence among and within the three phenotypically distinct subspecies of bare-nosed wombats (Vombatus ursinus) in south-east Australia that are presently-but were not historically-isolated by marine barriers. Using genome-wide single nucleotide polymorphisms, we identified three genetically distinct groups (mainland Australia, Bass Strait island, and Tasmania) corresponding to the recognized subspecies. However, isolation by distance was observed in the Tasmanian population, indicating additional constraints on gene flow can contribute to divergence in the absence of marine barriers, and may also explain genetic structuring among fragmented mainland populations. We additionally confirm origins and quantify the genetic divergence of an island population 46 years after the introduction of 21 individuals from the Vulnerable Bass Strait subspecies. In the light of our findings, we make recommendations for the maintenance of genetic variation and fitness across the species range.

Pathogens in space: Advancing understanding of pathogen dynamics and disease ecology through landscape genetics

Landscape genetics has provided many insights into how heterogeneous landscape features drive processes influencing spatial genetic variation in free-living organisms. This rapidly developing field has focused heavily on vertebrates, and expansion of this scope to the study of infectious diseases holds great potential for landscape geneticists and disease ecologists alike. The potential application of landscape genetics to infectious agents has garnered attention at formative stages in the development of landscape genetics, but systematic examination is lacking. We comprehensively review how landscape genetics is being used to better understand pathogen dynamics. We characterize the field and evaluate the types of questions addressed, approaches used and systems studied. We also review the now established landscape genetic methods and their realized and potential applications to disease ecology. Lastly, we identify emerging frontiers in the landscape genetic study of infectious agents, including recent phylogeographic approaches and frameworks for studying complex multihost and host-vector systems. Our review emphasizes the expanding utility of landscape genetic methods available for elucidating key pathogen dynamics (particularly transmission and spread) and also how landscape genetic studies of pathogens can provide insight into host population dynamics. Through this review, we convey how increasing awareness of the complementarity of landscape genetics and disease ecology among practitioners of each field promises to drive important cross-disciplinary advances.

Conservation of Sex-Linked Markers among Conspecific Populations of a Viviparous Skink, Niveoscincus ocellatus, Exhibiting Genetic and Temperature-Dependent Sex Determination

Sex determination systems are exceptionally diverse and have undergone multiple and independent evolutionary transitions among species, particularly reptiles. However, the mechanisms underlying these transitions have not been established. Here, we tested for differences in sex-linked markers in the only known reptile that is polymorphic for sex determination system, the spotted snow skink, Niveoscincus ocellatus, to quantify the genomic differences that have accompanied this transition. In a highland population, sex is determined genetically, whereas in a lowland population, offspring sex ratio is influenced by temperature. We found a similar number of sex-linked loci in each population, including shared loci, with genotypes consistent with male heterogamety (XY). However, population-specific linkage disequilibrium suggests greater differentiation of sex chromosomes in the highland population. Our results suggest that transitions between sex determination systems can be facilitated by subtle genetic differences.

The cascading pathogenic consequences of Sarcoptes scabiei infection that manifest in host disease

Sarcoptic mange, caused by the parasitic mite Sarcoptes scabiei, causes a substantive burden of disease to humans, domestic animals and wildlife, globally. There are many effects of S. scabiei infection, culminating in the disease which hosts suffer. However, major knowledge gaps remain on the pathogenic impacts of this infection. Here, we focus on the bare-nosed wombat host (Vombatus ursinus) to investigate the effects of mange on: (i) host heat loss and thermoregulation, (ii) field metabolic rates, (iii) foraging and resting behaviour across full circadian cycles, and (iv) fatty acid composition in host adipose, bone marrow, brain and muscle tissues. Our findings indicate that mange-infected V. ursinus lose more heat to the environment from alopecia-affected body regions than healthy individuals. Additionally, mange-infected individuals have higher metabolic rates in the wild. However, these metabolic demands are difficult to meet, because infected individuals spend less time foraging and more time inactive relative to their healthy counterparts, despite being outside of the burrow for longer. Lastly, mange infection results in altered fatty acid composition in adipose tissue, with increased amounts of omega-6 acids, and decreased amounts of omega-3 acids, a consequence of chronic cutaneous inflammation and inhibition of anti-inflammatory responses. These findings highlight the interactions of mange-induced physiological and behavioural changes, and have implications for the treatment and rehabilitation of infected individuals.

Detecting Selection on Temporal and Spatial Scales: A Genomic Time-Series Assessment of Selective Responses to Devil Facial Tumor Disease

Detecting loci under selection is an important task in evolutionary biology. In conservation genetics detecting selection is key to investigating adaptation to the spread of infectious disease. Loci under selection can be detected on a spatial scale, accounting for differences in demographic history among populations, or on a temporal scale, tracing changes in allele frequencies over time. Here we use these two approaches to investigate selective responses to the spread of an infectious cancer—devil facial tumor disease (DFTD)—that since 1996 has ravaged the Tasmanian devil (Sarcophilus harrisii). Using time-series ‘restriction site associated DNA’ (RAD) markers from populations pre- and post DFTD arrival, and DFTD free populations, we infer loci under selection due to DFTD and investigate signatures of selection that are incongruent among methods, populations, and times. The lack of congruence among populations influenced by DFTD with respect to inferred loci under selection, and the direction of that selection, fail to implicate a consistent selective role for DFTD. Instead genetic drift is more likely driving the observed allele frequency changes over time. Our study illustrates the importance of applying methods with different performance optima e.g. accounting for population structure and background selection, and assessing congruence of the results.

Coalescent Modelling Suggests Recent Secondary-Contact of Cryptic Penguin Species

Molecular genetic analyses present powerful tools for elucidating demographic and biogeographic histories of taxa. Here we present genetic evidence showing a dynamic history for two cryptic lineages within Eudyptula, the world's smallest penguin. Specifically, we use a suite of genetic markers to reveal that two congeneric taxa ('Australia' and 'New Zealand') co-occur in southern New Zealand, with only low levels of hybridization. Coalescent modelling suggests that the Australian little penguin only recently expanded into southern New Zealand. Analyses conducted under time-dependent molecular evolutionary rates lend support to the hypothesis of recent anthropogenic turnover, consistent with shifts detected in several other New Zealand coastal vertebrate taxa. This apparent turnover event highlights the dynamic nature of the region’s coastal ecosystem.

Extensive population decline in the Tasmanian devil predates European settlement and devil facial tumour disease

The Tasmanian devil (Sarcophilus harrisii) was widespread in Australia during the Late Pleistocene but is now endemic to the island of Tasmania. Low genetic diversity combined with the spread of devil facial tumour disease have raised concerns for the species' long-term survival. Here, we investigate the origin of low genetic diversity by inferring the species' demographic history using temporal sampling with summary statistics, full-likelihood and approximate Bayesian computation methods. Our results show extensive population declines across Tasmania correlating with environmental changes around the last glacial maximum and following unstable climate related to increased 'El Niño-Southern Oscillation' activity.

Isolation via next-generation sequencing of microsatellites from the Tasmanian macroalgae Lessonia corrugata (Lessoniaceae)

PREMISE OF THE STUDY

Microsatellite markers for the macroalgae Lessonia corrugata (Lessoniaceae) were developed, for the first time, to enable population genetic assessment of this important foundation species.

METHODS AND RESULTS

Ion Torrent sequencing identified 16,622 loci, 29 of which were trialed in L. corrugata. Seven loci were found to be polymorphic and screened for variation in 76 individuals from two populations in Tasmania, Australia. Observed heterozygosity ranged from 0.086 to 0.686 (mean 0.386) and the number of alleles per locus ranged from two to five (mean 2.57). Heterozygosity was not significantly different from expected values.

CONCLUSIONS

These loci can be used to study the population genetics of L. corrugata, a key habitat-forming species in the Tasmanian marine ecosystem, and will help to assess gene flow among spatially discrete populations such as those in marine protected areas.

Persistence and dispersal in a Southern Hemisphere glaciated landscape: the phylogeography of the spotted snow skink (Niveoscincus ocellatus) in Tasmania

Background

The aim of this research was to identify the effects of Pleistocene climate change on the distribution of fauna in Tasmania, and contrast this with biotic responses in other temperate regions in the Northern and Southern Hemisphere that experienced glacial activity during this epoch. This was achieved by examining the phylogeographic patterns in a widely distributed Tasmanian endemic reptile, Niveoscincus ocellatus. 204 individuals from 29 populations across the distributional range of N. ocellatus were surveyed for variation at two mitochondrial genes (ND2, ND4), and two nuclear genes (β-globin, RPS8). Phylogenetic relationships were reconstructed using a range of methods (maximum parsimony, Bayesian inference and haplotype networks), and the demographic histories of populations were assessed (AMOVA, Tajima’s D, Fu’s Fs, mismatch distributions, extended Bayesian skyline plots, and relaxed random walk analyses).

Results

There was a high degree of mitochondrial haplotype diversity (96 unique haplotypes) and phylogeographic structure, where spatially distinct groups were associated with Tasmania’s Northeast and a large area covering Southeast and Central Tasmania. Phylogeographic structure was also present within each major group, but the degree varied regionally, being highest in the Northeast. Only the Southeastern group had a signature of demographic expansion, occurring during the Pleistocene but post-dating the Last Glacial Maximum. In contrast, nuclear DNA had low levels of variation and a lack of phylogeographic structure, and further loci should be surveyed to corroborate the mitochondrial inferences.

Conclusions

The phylogeographic patterns of N. ocellatus indicate Pleistocene range and demographic expansion in N. ocellatus, particularly in the Southeast and Central areas of Tasmania. Expansion in Central and Southeastern areas appears to have been more recent in both demographic and spatial contexts, than in Northeast Tasmania, which is consistent with inferences for other taxa of greater stability and persistence in Northeast Tasmania during the Last Glacial Maximum. These phylogeographic patterns indicate contrasting demographic histories of populations in close proximity to areas directly affected by glaciers in the Southern Hemisphere during the LGM.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0397-y) contains supplementary material, which is available to authorized users.

Nonequilibrium Conditions Explain Spatial Variability in Genetic Structuring of Little Penguin (Eudyptula minor)

Factors responsible for spatial structuring of population genetic variation are varied, and in many instances there may be no obvious explanations for genetic structuring observed, or those invoked may reflect spurious correlations. A study of little penguins (Eudyptula minor) in southeast Australia documented low spatial structuring of genetic variation with the exception of colonies at the western limit of sampling, and this distinction was attributed to an intervening oceanographic feature (Bonney Upwelling), differences in breeding phenology, or sea level change. Here, we conducted sampling across the entire Australian range, employing additional markers (12 microsatellites and mitochondrial DNA, 697 individuals, 17 colonies). The zone of elevated genetic structuring previously observed actually represents the eastern half of a genetic cline, within which structuring exists over much shorter spatial scales than elsewhere. Colonies separated by as little as 27 km in the zone are genetically distinguishable, while outside the zone, homogeneity cannot be rejected at scales of up to 1400 km. Given a lack of additional physical or environmental barriers to gene flow, the zone of elevated genetic structuring may reflect secondary contact of lineages (with or without selection against interbreeding), or recent colonization and expansion from this region. This study highlights the importance of sampling scale to reveal the cause of genetic structuring.

Implications of macroalgal isolation by distance for networks of marine protected areas

The global extent of macroalgal forests is declining, greatly affecting marine biodiversity at broad scales through the effects macroalgae have on ecosystem processes, habitat provision, and food web support. Networks of marine protected areas comprise one potential tool that may safeguard gene flow among macroalgal populations in the face of increasing population fragmentation caused by pollution, habitat modification, climate change, algal harvesting, trophic cascades, and other anthropogenic stressors. Optimal design of protected area networks requires knowledge of effective dispersal distances for a range of macroalgae. We conducted a global meta-analysis based on data in the published literature to determine the generality of relation between genetic differentiation and geographic distance among macroalgal populations. We also examined whether spatial genetic variation differed significantly with respect to higher taxon, life history, and habitat characteristics. We found clear evidence of population isolation by distance across a multitude of macroalgal species. Genetic and geographic distance were positively correlated across 49 studies; a modal distance of 50-100 km maintained F(ST) < 0.2. This relation was consistent for all algal divisions, life cycles, habitats, and molecular marker classes investigated. Incorporating knowledge of the spatial scales of gene flow into the design of marine protected area networks will help moderate anthropogenic increases in population isolation and inbreeding and contribute to the resilience of macroalgal forests.

Did postglacial sea-level changes initiate the evolutionary divergence of a Tasmanian endemic raptor from its mainland relative?

Populations on continental islands are often distinguishable from mainland conspecifics with respect to body size, appearance, behaviour or life history, and this is often congruent with genetic patterns. It is commonly assumed that such differences developed following the complete isolation of populations by sea-level rise following the Last Glacial Maximum (LGM). However, population divergence may predate the LGM, or marine dispersal and colonization of islands may have occurred more recently; in both cases, populations may have also diverged despite ongoing gene flow. Here, we test these alternative hypotheses for the divergence between wedge-tailed eagles from mainland Australia (Aquila audax audax) and the threatened Tasmanian subspecies (Aquila audax fleayi), based on variation at 20 microsatellite loci and mtDNA. Coalescent analyses indicate that population divergence appreciably postdates the severance of terrestrial habitat continuity and occurred without any subsequent gene flow. We infer a recent colonization of Tasmania by marine dispersal and cannot discount founder effects as the cause of differences in body size and life history. We call into question the general assumption of post-LGM marine transgression as the initiator of divergence of terrestrial lineages on continental islands and adjacent mainland, and highlight the range of alternative scenarios that should be considered.

Optimizing the use of shed feathers for genetic analysis

Shed feathers obtained by noninvasive genetic sampling (NGS) are a valuable source of DNA for genetic studies of birds. They can be collected across a large geographical range and facilitate research on species that would otherwise be extremely difficult to study. A limitation of this approach is uncertainty concerning the quality of the extracted DNA. Here we investigate the relationship between feather type, feather condition and DNA quality (amplification success) in order to provide a simple, cost-effective method for screening samples prior to genetic analysis. We obtained 637 shed feathers of the powerful owl (Ninox strenua) from across its range in southeastern Australia. The extracted DNA was amplified using polymerase chain reaction for a range of markers including mitochondrial DNA, ND3 and nuclear DNA, a simple sequence repeat (Nst02) and a portion of the CHD-1 gene (P2/P8). We found that feather condition significantly influenced the amplification success of all three loci, with feathers characterized as 'good' having greater success. Feather type was found to be of lower importance, with good quality feathers of all types consistently producing high success for all three loci. We also found that the successful amplification of multilocus genotypes was dependant on the condition of the starting material and was highly correlated with successful amplification of the sex-linked CHD-1 locus. Samples with low DNA quality have a higher probability of amplification failure and are more likely to produce incorrect genotypes; therefore, identifying samples with high DNA quality can save substantial time and cost associated with the genetic analysis of NGS. As a result, we propose a method for screening shed feathers in order to provide a subset of samples which will have a greater probability of containing high quality DNA suitable for the amplification of multilocus genotypes.

Understanding age-specific dispersal in fishes through hydrodynamic modelling, genetic simulations and microsatellite DNA analysis

Many marine species have vastly different capacities for dispersal during larval, juvenile and adult life stages, and this has the potential to complicate the identification of population boundaries and the implementation of effective management strategies such as marine protected areas. Genetic studies of population structure and dispersal rarely disentangle these differences and usually provide only lifetime-averaged information that can be considered by managers. We address this limitation by combining age-specific autocorrelation analysis of microsatellite genotypes, hydrodynamic modelling and genetic simulations to reveal changes in the extent of dispersal during the lifetime of a marine fish. We focus on an exploited coral reef species, Lethrinus nebulosus, which has a circum-tropical distribution and is a key component of a multispecies fishery in northwestern Australia. Conventional population genetic analyses revealed extensive gene flow in this species over vast distances (up to 1,500 km). Yet, when realistic adult dispersal behaviours were modelled, they could not account for these observations, implying adult dispersal does not dominate gene flow. Instead, hydrodynamic modelling showed that larval L. nebulosus are likely to be transported hundreds of kilometres, easily accounting for the observed gene flow. Despite the vast scale of larval transport, juvenile L. nebulosus exhibited fine-scale genetic autocorrelation, which declined with age. This implies both larval cohesion and extremely limited juvenile dispersal prior to maturity. The multidisciplinary approach adopted in this study provides a uniquely comprehensive insight into spatial processes in this marine fish.

Resolution of the Acanthopagrus black seabream complex based on mitochondrial and amplified fragment-length polymorphism analyses

In this study, DNA analyses were employed to verify the identity of six morphologically similar species that occur in the coastal waters of Taiwan: the black seabream complex (Acanthopagrus latus, Acanthopagrus schlegelii, Acanthopagrus sivicolus, Acanthopagrus taiwanensis, Acanthopagrus chinshira and Acanthopagrus pacificus). Amplified fragment-length polymorphism (AFLP) analyses clearly distinguished the same six species that are morphologically diagnosable based on subtle differences in scale counts and anal-fin colouration. In contrast, mitochondrial DNA analyses based on cytochrome b gene sequences did not distinguish individuals of A. schlegelii and A. sivicolus, reflecting either historical introgression or recent speciation and incomplete sorting of their mitochondrial lineages. Phylogenetic relationships among these six north-west Pacific Ocean species of Acanthopagrus analysed using AFLP data were consistent with scale rows above the lateral line (TRac), sperm ultrastructure and geographical distribution. The study provides molecular tools for future research relevant to improved management of these resources, and an increased understanding of the evolutionary history of this radiation.

A hybrid zone and bidirectional introgression between two catadromous species: Australian bass Macquaria novemaculeata and estuary perch Macquaria colonorum

The presence and distribution of hybrid individuals and the existence of a hybrid zone between the catadromous Australian bass Macquaria novemaculeata and estuary perch Macquaria colonorum were investigated throughout the range of both species in Australia. Bayesian analyses and genotypic simulations identified 140 putative hybrids (11·5% of the total sample) with varying levels of introgression. Most hybrids were observed in an area extending from the Snowy River to the Albert River suggesting a hybrid zone in the eastern Bass Strait region. Sixteen hybrids, however, were found outside this zone, possibly reflecting the movement of hybrid offspring between estuaries or their inadvertent release during fish stocking programmes. Biparental backcrossing was found to occur suggesting that hybrids were fertile. These results have implications for the management of the extensive stocking programme in M. novemaculeata and for understanding the potential role of habitat degradation and reduced water flow in facilitating hybridization in species with migratory life histories.

Two spinefoot colour morphs: mottled spinefoot Siganus fuscescens and white-spotted spinefoot Siganus canaliculatus are synonyms

Mottled spinefoot Siganus fuscescens and white-spotted spinefoot Siganus canaliculatus are two similar species that differ subtly in colouration and morphology. Three major mtDNA clades were identified for these species, but individuals were clustered by amplified fragment length polymorphism (AFLP) according to geography rather than morphology, suggesting that the colour morphs are interbreeding.