Opportunities and challenges in assessing climate change vulnerability through genomics.

By investigating how past selection has affected allele frequencies across space, genomic tools are providing new insights into adaptive evolutionary processes. Now researchers are considering how this genomic information can be used to predict the future vulnerability of species under climate change. Genomic vulnerability assessments show promise, but challenges remain.

Conservation genetics as a management tool: The five best-supported paradigms to assist the management of threatened species.

About 50 y ago, Crow and Kimura [An Introduction to Population Genetics Theory (1970)] and Ohta and Kimura [Genet. Res. 22, 201-204 (1973)] laid the foundations of conservation genetics by predicting the relationship between population size and genetic marker diversity. This work sparked an enormous research effort investigating the importance of population dynamics, in particular small population size, for population mean performance, population viability, and evolutionary potential. In light of a recent perspective [J. C. Teixeira, C. D. Huber, Proc. Natl. Acad. Sci. U.S.A. 118, 10 (2021)] that challenges some fundamental assumptions in conservation genetics, it is timely to summarize what the field has achieved, what robust patterns have emerged, and worthwhile future research directions. We consider theory and methodological breakthroughs that have helped management, and we outline some fundamental and applied challenges for conservation genetics.

Heterogeneous patterns of heterozygosity loss in isolated populations of the threatened eastern barred bandicoot (Perameles gunnii).

Identifying and analysing isolated populations is critical for conservation. Isolation can make populations vulnerable to local extinction due to increased genetic drift and inbreeding, both of which should leave imprints of decreased genome-wide heterozygosity. While decreases in heterozygosity among populations are frequently investigated, fewer studies have analysed how heterozygosity varies among individuals, including whether heterozygosity varies geographically along lines of discrete population structure or with continuous patterns analogous to isolation by distance. Here we explore geographical patterns of differentiation and individual heterozygosity in the threatened eastern barred bandicoot (Perameles gunnii) in Tasmania, Australia, using genomic data from 85 samples collected between 2008 and 2011. Our analyses identified two isolated demes undergoing significant genetic drift, and several areas of fine-scale differentiation across Tasmania. We observed discrete genetic structures across geographical barriers and continuous patterns of isolation by distance, with little evidence of recent or historical migration. Using a recently developed analytical pipeline for estimating autosomal heterozygosity, we found individual heterozygosities varied within demes by up to a factor of two, and demes with low-heterozygosity individuals also still contained those with high heterozygosity. Spatial interpolation of heterozygosity scores clarified these patterns and identified the isolated Tasman Peninsula as a location where low-heterozygosity individuals were more common than elsewhere. Our results provide novel insights into the relationship between isolation-driven genetic structure and local heterozygosity patterns. These may help improve translocation efforts, by identifying populations in need of assistance, and by providing an individualised metric for identifying source animals for translocation.

A molecular method for biomonitoring of an exotic plant-pest: Leafmining for environmental DNA.

Understanding how invasive species respond to novel environments is limited by a lack of sensitivity and throughput in conventional biomonitoring methods. Arthropods in particular are often difficult to monitor due to their small size, rapid lifecycles, and/or visual similarities with co-occurring species. This is true for the agromyzid leafminer fly, Liriomyza sativae, a global pest of vegetable and nursery industries that has recently established in Australia. A robust method based on environmental DNA (eDNA) was developed exploiting traces of DNA left inside "empty" leaf mines, which are straightforward to collect and persist longer in the environment than the fly. This extends the window of possible diagnosis to at least 28 days after a leaf mine becomes empty. The test allowed for visually indistinguishable leafmining damage caused by L. sativae to be genetically differentiated from that of other flies. Field application resulted in the identification of new local plant hosts for L. sativae, including widely distributed weeds and common garden crops, which has important implications for the pest's ability to spread. Moreover, the test confirmed the presence of a previously unknown population of L. sativae on an island in the Torres Strait. The developed eDNA method is likely to become an important tool for L. sativae and other leafmining species of biosecurity significance, which, historically, have been difficult to detect, diagnose and monitor. More generally, eDNA is emerging as a highly sensitive and labour-efficient surveillance tool for difficult to survey species to improve outcomes for agricultural industries, global health, and the environment.

Multispecies models reveal that eDNA metabarcoding is more sensitive than backpack electrofishing for conducting fish surveys in freshwater streams.

Environmental DNA (eDNA) sampling can provide accurate, cost-effective, landscape-level data on species distributions. Previous studies have compared the sensitivity of eDNA sampling to traditional sampling methods for single species, but similar comparative studies on multispecies eDNA metabarcoding are rare. Using hierarchical site occupancy detection models, we examined whether key choices associated with eDNA metabarcoding (primer selection, low-abundance read filtering and the number of positive water samples used to classify a species as present at a site) affect the sensitivity of metabarcoding, relative to backpack electrofishing for fish in freshwater streams. Under all scenarios (teleostei and vertebrate primers; 0%, 0.1% and 1% read filtering thresholds; one or two positive samples required to classify species as present), we found that eDNA metabarcoding is, on average, more sensitive than electrofishing. Combining vertebrate and teleostei markers resulted in higher detection probabilities relative to the use of either marker in isolation. Increasing the threshold used to filter low-abundance reads decreased species detection probabilities but did not change our overall finding that eDNA metabarcoding was more sensitive than electrofishing. Using a threshold of two positive water samples (out of five) to classify a species as present typically had negligible effects on detection probabilities compared to using one positive water sample. Our findings demonstrate that eDNA metabarcoding is generally more sensitive than electrofishing for conducting fish surveys in freshwater streams, and that this outcome is not sensitive to methodological decisions associated with metabarcoding.

Heterogeneous genetic invasions of three insecticide resistance mutations in Indo-Pacific populations of Aedes aegypti (L.).

Nations throughout the Indo-Pacific region use pyrethroid insecticides to control Aedes aegypti, the mosquito vector of dengue, often without knowledge of pyrethroid resistance status of the pest or origin of resistance. Two mutations (V1016G + F1534C) in the sodium channel gene (Vssc) of Ae. aegypti modify ion channel function and cause target-site resistance to pyrethroid insecticides, with a third mutation (S989P) having a potential additive effect. Of 27 possible genotypes involving these mutations, some allelic combinations are never seen whereas others predominate. Here, five allelic combinations common in Ae. aegypti from the Indo-Pacific region are described and their geographical distributions investigated using genome-wide SNP markers. We tested the hypothesis that resistance allele combinations evolved de novo in populations versus the alternative that dispersal of Ae. aegypti between populations facilitated genetic invasions of allele combinations. We used latent factor mixed-models to detect SNPs throughout the genome that showed structuring in line with resistance allele combinations and compared variation at SNPs within the Vssc gene with genome-wide variation. Mixed-models detected an array of SNPs linked to resistance allele combinations, all located within or in close proximity to the Vssc gene. Variation at SNPs within the Vssc gene was structured by resistance profile, whereas genome-wide SNPs were structured by population. These results demonstrate that alleles near to resistance mutations have been transferred between populations via linked selection. This indicates that genetic invasions have contributed to the widespread occurrence of Vssc allele combinations in Ae. aegypti in the Indo-Pacific region, pointing to undocumented mosquito invasions between countries.

Shifting paradigms in restoration of the world's coral reefs.

Many ecosystems around the world are rapidly deteriorating due to both local and global pressures, and perhaps none so precipitously as coral reefs. Management of coral reefs through maintenance (e.g., marine-protected areas, catchment management to improve water quality), restoration, as well as global and national governmental agreements to reduce greenhouse gas emissions (e.g., the 2015 Paris Agreement) is critical for the persistence of coral reefs. Despite these initiatives, the health and abundance of corals reefs are rapidly declining and other solutions will soon be required. We have recently discussed options for using assisted evolution (i.e., selective breeding, assisted gene flow, conditioning or epigenetic programming, and the manipulation of the coral microbiome) as a means to enhance environmental stress tolerance of corals and the success of coral reef restoration efforts. The 2014-2016 global coral bleaching event has sharpened the focus on such interventionist approaches. We highlight the necessity for consideration of alternative (e.g., hybrid) ecosystem states, discuss traits of resilient corals and coral reef ecosystems, and propose a decision tree for incorporating assisted evolution into restoration initiatives to enhance climate resilience of coral reefs.

Conservation of genetic uniqueness of populations may increase extinction likelihood of endangered species: the case of Australian mammals.

BACKGROUND: As increasingly fragmented and isolated populations of threatened species become subjected to climate change, invasive species and other stressors, there is an urgent need to consider adaptive potential when making conservation decisions rather than focussing on past processes. In many cases, populations identified as unique and currently managed separately suffer increased risk of extinction through demographic and genetic processes. Other populations currently not at risk are likely to be on a trajectory where declines in population size and fitness soon appear inevitable. RESULTS: Using datasets from natural Australian mammal populations, we show that drift processes are likely to be driving uniqueness in populations of many threatened species as a result of small population size and fragmentation. Conserving and managing such remnant populations separately will therefore often decrease their adaptive potential and increase species extinction risk. CONCLUSIONS: These results highlight the need for a paradigm shift in conservation biology practise; strategies need to focus on the preservation of genetic diversity at the species level, rather than population, subspecies or evolutionary significant unit. The introduction of new genetic variants into populations through in situ translocation needs to be considered more broadly in conservation programs as a way of decreasing extinction risk by increasing neutral genetic diversity which may increase the adaptive potential of populations if adaptive variation is also increased.

Rapid sequential spread of two Wolbachia variants in Drosophila simulans.

The maternally inherited intracellular bacteria Wolbachia can manipulate host reproduction in various ways that foster frequency increases within and among host populations. Manipulations involving cytoplasmic incompatibility (CI), where matings between infected males and uninfected females produce non-viable embryos, are common in arthropods and produce a reproductive advantage for infected females. CI was associated with the spread of Wolbachia variant wRi in Californian populations of Drosophila simulans, which was interpreted as a bistable wave, in which local infection frequencies tend to increase only once the infection becomes sufficiently common to offset imperfect maternal transmission and infection costs. However, maternally inherited Wolbachia are expected to evolve towards mutualism, and they are known to increase host fitness by protecting against infectious microbes or increasing fecundity. We describe the sequential spread over approximately 20 years in natural populations of D. simulans on the east coast of Australia of two Wolbachia variants (wAu and wRi), only one of which causes significant CI, with wRi displacing wAu since 2004. Wolbachia and mtDNA frequency data and analyses suggest that these dynamics, as well as the earlier spread in California, are best understood as Fisherian waves of favourable variants, in which local spread tends to occur from arbitrarily low frequencies. We discuss implications for Wolbachia-host dynamics and coevolution and for applications of Wolbachia to disease control.

Detection of Low-Level Cardinium and Wolbachia Infections in Culicoides.

Bacterial endosymbionts have been identified as potentially useful biological control agents for a range of invertebrate vectors of disease. Previous studies of Culicoides (Diptera: Ceratopogonidae) species using conventional PCR assays have provided evidence of Wolbachia (1/33) and Cardinium (8/33) infections. Here, we screened 20 species of Culicoides for Wolbachia and Cardinium, utilizing a combination of conventional PCR and more sensitive quantitative PCR (qPCR) assays. Low levels of Cardinium DNA were detected in females of all but one of the Culicoides species screened, and low levels of Wolbachia were detected in females of 9 of the 20 Culicoides species. Sequence analysis based on partial 16S rRNA gene and gyrB sequences identified "Candidatus Cardinium hertigii" from group C, which has previously been identified in Culicoides from Japan, Israel, and the United Kingdom. Wolbachia strains detected in this study showed 98 to 99% sequence identity to Wolbachia previously detected from Culicoides based on the 16S rRNA gene, whereas a strain with a novel wsp sequence was identified in Culicoides narrabeenensis. Cardinium isolates grouped to geographical regions independent of the host Culicoides species, suggesting possible geographical barriers to Cardinium movement. Screening also identified Asaia bacteria in Culicoides. These findings point to a diversity of low-level endosymbiont infections in Culicoides, providing candidates for further characterization and highlighting the widespread occurrence of these endosymbionts in this insect group.

I Environmental DNA sampling is more sensitive than a traditional survey technique for detecting an aquatic invader.

Effective management of alien species requires detecting populations in the early stages of invasion. Environmental DNA (eDNA) sampling can detect aquatic species at relatively low densities, but few studies have directly compared detection probabilities of eDNA sampling with those of traditional sampling methods. We compare the ability of a traditional sampling technique (bottle trapping) and eDNA to detect a recently established invader, the smooth newt Lissotriton vulgaris vulgaris, at seven field sites in Melbourne, Australia. Over a four-month period, per-trap detection probabilities ranged from 0.01 to 0.26 among sites where L. v. vulgaris was detected, whereas per-sample eDNA estimates were much higher (0.29-1.0). Detection probabilities of both methods varied temporally (across days and months), but temporal variation appeared to be uncorrelated between methods. Only estimates of spatial variation were strongly correlated across the two sampling techniques. Environmental variables (water depth, rainfall, ambient temperature) were not clearly correlated with detection probabilities estimated via trapping, whereas eDNA detection probabilities were negatively correlated with water depth, possibly reflecting higher eDNA concentrations at lower water levels. Our findings demonstrate that eDNA sampling can be an order of magnitude more sensitive than traditional methods, and illustrate that traditional- and eDNA-based surveys can provide independent information on species distributions when occupancy surveys are conducted over short timescales.

Genome-wide SNPs lead to strong signals of geographic structure and relatedness patterns in the major arbovirus vector, Aedes aegypti.

BACKGROUND: Genetic markers are widely used to understand the biology and population dynamics of disease vectors, but often markers are limited in the resolution they provide. In particular, the delineation of population structure, fine scale movement and patterns of relatedness are often obscured unless numerous markers are available. To address this issue in the major arbovirus vector, the yellow fever mosquito (Aedes aegypti), we used double digest Restriction-site Associated DNA (ddRAD) sequencing for the discovery of genome-wide single nucleotide polymorphisms (SNPs). We aimed to characterize the new SNP set and to test the resolution against previously described microsatellite markers in detecting broad and fine-scale genetic patterns in Ae. aegypti. RESULTS: We developed bioinformatics tools that support the customization of restriction enzyme-based protocols for SNP discovery. We showed that our approach for RAD library construction achieves unbiased genome representation that reflects true evolutionary processes. In Ae. aegypti samples from three continents we identified more than 18,000 putative SNPs. They were widely distributed across the three Ae. aegypti chromosomes, with 47.9% found in intergenic regions and 17.8% in exons of over 2,300 genes. Pattern of their imputed effects in ORFs and UTRs were consistent with those found in a recent transcriptome study. We demonstrated that individual mosquitoes from Indonesia, Australia, Vietnam and Brazil can be assigned with a very high degree of confidence to their region of origin using a large SNP panel. We also showed that familial relatedness of samples from a 0.4 km2 area could be confidently established with a subset of SNPs. CONCLUSIONS: Using a cost-effective customized RAD sequencing approach supported by our bioinformatics tools, we characterized over 18,000 SNPs in field samples of the dengue fever mosquito Ae. aegypti. The variants were annotated and positioned onto the three Ae. aegypti chromosomes. The new SNP set provided much greater resolution in detecting population structure and estimating fine-scale relatedness than a set of polymorphic microsatellites. RAD-based markers demonstrate great potential to advance our understanding of mosquito population processes, critical for implementing new control measures against this major disease vector.

Microsatellite loci and the complete mitochondrial DNA sequence characterized through next generation sequencing and de novo genome assembly for the critically endangered orange-bellied parrot, Neophema chrysogaster.

A suite of polymorphic microsatellite markers and the complete mitochondrial genome sequence was developed by next generation sequencing (NGS) for the critically endangered orange-bellied parrot, Neophema chrysogaster. A total of 14 polymorphic loci were identified and characterized using DNA extractions representing 40 individuals from Melaleuca, Tasmania, sampled in 2002. We observed moderate genetic variation across most loci (mean number of alleles per locus = 2.79; mean expected heterozygosity = 0.53) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. De novo and reference-based genome assemblies performed using MIRA were used to assemble the N. chrysogaster mitochondrial genome sequence with mean coverage of 116-fold (range 89 to 142-fold). The mitochondrial genome consists of 18,034 base pairs, and a typical metazoan mitochondrial gene content consisting of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a single large non-coding region (control region). The arrangement of mitochondrial genes is also typical of Avian taxa. The annotation of the mitochondrial genome and the characterization of 14 microsatellite markers provide a valuable resource for future genetic monitoring of wild and captive N. chrysogaster populations. As found previously, NGS provides a rapid, low cost and reliable method for polymorphic nuclear genetic marker development and determining complete mitochondrial genome sequences when only a fraction of a genome is sequenced.

The development of 10 novel polymorphic microsatellite markers through next generation sequencing and a preliminary population genetic analysis for the endangered Glenelg spiny crayfish, Euastacus bispinosus.

The Glenelg spiny crayfish, Euastacus bispinosus, is an iconic freshwater invertebrate of south eastern Australia and listed as 'endangered' under the Environment Protection and Biodiversity Conservation Act 1999, and 'vulnerable' under the International Union for Conservation of Nature's Red List. The species has suffered major population declines as a result of over-fishing, low environmental flows, the introduction of invasive fish species and habitat degradation. In order to develop an effective conservation strategy, patterns of gene flow, genetic structure and genetic diversity across the species distribution need to be clearly understood. In this study we develop a suite of polymorphic microsatellite markers by next generation sequencing. A total of 15 polymorphic loci were identified and 10 characterized using 22 individuals from the lower Glenelg River. We observed low to moderate genetic variation across most loci (mean number of alleles per locus = 2.80; mean expected heterozygosity = 0.36) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. Individuals from two additional sites (Crawford River, Victoria; Ewens Ponds Conservation Park, South Australia) were genotyped at all 10 loci and a preliminary investigation of genetic diversity and population structure was undertaken. Analyses indicate high levels of genetic differentiation among sample locations (F ST = 0.49), while the Ewens Ponds population is genetically homogeneous, indicating a likely small founder group and ongoing inbreeding. Management actions will be needed to restore genetic diversity in this and possibly other at risk populations. These markers will provide a valuable resource for future population genetic assessments so that an effective framework can be developed for implementing conservation strategies for E. bispinosus.

High-throughput PCR assays to monitor Wolbachia infection in the dengue mosquito (Aedes aegypti) and Drosophila simulans.

We have developed and validated two new fluorescence-based PCR assays to detect the Wolbachia wMel strain in Aedes aegypti and the wRi and wAu strains in Drosophila simulans. The new assays are accurate, informative, and cost-efficient for large-scale Wolbachia screening.

Australian Bryobia mites (Trombidiformes: Tetranychidae) form a complex of cryptic taxa with unique climatic niches and insecticide responses.

BACKGROUND: Bryobia (Koch) mites belong to the economically important spider mite family, the Tetranychidae, with >130 species described worldwide. Due to taxonomic difficulties and most species being asexual, species identification relies heavily on genetic markers. Multiple putative Bryobia mite species have been identified attacking pastures and grain crops in Australia. In this study, we collected 79 field populations of Bryobia mites and combined these with 134 populations that were collected previously. We characterised taxonomic variation of mites using 28S rDNA amplicon-based DNA metabarcoding using next-generation sequencing approaches and direct Sanger sequencing. We then undertook species distribution modelling of the main genetic lineages and examined the chemical responses of multiple field populations. RESULTS: We identified 47 unique haplotypes across all mites sampled that grouped into four distinct genetic lineages. These lineages have different distributions, with three of the four putative lineages showing different climatic envelopes, as inferred from species distribution modelling. Bryobia mite populations also showed different responses to a widely used insecticide (the organophosphate, omethoate), but not to another chemical (the pyrethroid, bifenthrin) when examined using laboratory bioassays. CONCLUSION: Our findings indicate that cryptic diversity is likely to complicate the formulation of management strategies for Bryobia mites. Although focussed on Australia, this study demonstrates the challenges of studying Bryobia and highlights the importance of further research into this complex group of mites across the world. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Genetic structure and long-distance dispersal in populations of the wingless pest springtail, Sminthurus viridis (Collembola: Sminthuridae).

The lucerne flea, Sminthurus viridis (Collembola: Sminthuridae) (L.) is a major pest of broadacre agriculture across southern Australia. Few molecular studies have been conducted on S. viridis and none have examined its population genetics, despite the importance for developing effective control strategies. Here, we characterize the genetic structure of Australian populations using three allozyme and eight microsatellite loci, as well as sequencing a fragment of the mitochondrial DNA cytochrome oxidase I gene. We found that S. viridis in Australia are diploid, sexually reproducing and exhibit significant population structure as a result of limited gene flow. Despite significant differentiation between populations, there was very low cytochrome oxidase subunit I (COI) gene sequence variation, indicating the presence of a single species in Australia. The observed structure only marginally complied with an 'isolation by distance' model with human-mediated long-distance dispersal likely occurring. Allozymes and microsatellites gave very similar FST estimates, although differences found for novel alternative estimates of differentiation suggest that the allozymes did not capture the full extent of the population structure. These results highlight that control strategies may need to vary for locally adapted S. viridis populations and strategies aimed at limiting the spread of any future pesticide resistance will need to manage the effects of human-mediated dispersal.

Frequency-dependent selection maintains clonal diversity in an asexual organism.

Asexual organisms can be genetically variable and evolve through time, yet it is not known how genetic diversity is maintained in populations. In sexual organisms, negative frequency-dependent selection plays a role in maintaining diversity at some loci, but in asexual organisms, this mechanism could provide a general explanation for persistent genetic diversity because it acts on the whole genome and not just on some polymorphisms within a genome. Using field manipulations, we show that negative frequency-dependent selection maintains clonal diversity in an asexual mite species, and we link predicted equilibrium clonal frequencies to average frequencies in space and time. Intense frequency-dependent selection is likely to be a general mechanism for persistent genetic diversity in asexual organisms.

Population dynamics and diapause response of the springtail pest Sminthurus viridis (Collembola: Sminthuridae) in southeastern Australia.

The springtail, Sminthurus viridis (L.) (Collembola: Sminthuridae), is an important agricultural pest across southern Australia. We investigated the seasonal abundance patterns and summer diapause response of S. viridis in southeastern Australia by using field and shadehouse (a greenhouse that offers seedlings shade) experiments. Seasonal activity patterns of S. viridis were largely consistent with previous studies, with the pest active from autumn to spring. In addition, the timing and pattern of the summer-diapausing egg stage was established, with multiple generations probably producing diapause eggs. A strong relationship between soil moisture and temperature with autumn emergence also was observed. These results suggest that S. viridis autumn pest pressure can be predicted and indicate that late-season spraying strategies currently used for a sympatric agricultural pest are unlikely to be as effective against S. viridis.

Genetic mixing for population management: From genetic rescue to provenancing.

Animal and plant species around the world are being challenged by the deleterious effects of inbreeding, loss of genetic diversity, and maladaptation due to widespread habitat destruction and rapid climate change. In many cases, interventions will likely be needed to safeguard populations and species and to maintain functioning ecosystems. Strategies aimed at initiating, reinstating, or enhancing patterns of gene flow via the deliberate movement of genotypes around the environment are generating growing interest with broad applications in conservation and environmental management. These diverse strategies go by various names ranging from genetic or evolutionary rescue to provenancing and genetic resurrection. Our aim here is to provide some clarification around terminology and to how these strategies are connected and linked to underlying genetic processes. We draw on case studies from the literature and outline mechanisms that underlie how the various strategies aim to increase species fitness and impact the wider community. We argue that understanding mechanisms leading to species decline and community impact is a key to successful implementation of these strategies. We emphasize the need to consider the nature of source and recipient populations, as well as associated risks and trade-offs for the various strategies. This overview highlights where strategies are likely to have potential at population, species, and ecosystem scales, but also where they should probably not be attempted depending on the overall aims of the intervention. We advocate an approach where short- and long-term strategies are integrated into a decision framework that also considers nongenetic aspects of management.