Quantifying uncertainty in inferences of landscape genetic resistance due to choice of individual-based genetic distance metric

Estimates of gene flow resulting from landscape resistance inferences frequently inform conservation management decision-making processes. Therefore, results must be robust across approaches and reflect real-world gene flow instead of methodological artefacts. Here, we tested the impact of 32 individual-based genetic distance metrics on the robustness and accuracy of landscape resistance modelling results. We analysed three empirical microsatellite datasets and 36 simulated datasets that varied in landscape resistance and genetic spatial autocorrelation. We used ResistanceGA to generate optimised multi-feature resistance surfaces for each of these datasets using 32 different genetic distance metrics. Results of the empirical dataset demonstrated that the choice of genetic distance metric can have strong impacts on inferred optimised resistance surfaces. Simulations showed accurate parametrisation of resistance surfaces across most genetic distance metrics only when a small number of environmental features was impacting gene flow. Landscape scenarios with many features impacting gene flow led to a generally poor recovery of true resistance surfaces. Simulation results also emphasise that choosing a genetic distance metric should not be based on marginal R2 -based model fit. Until more robust methods are available, resistance surfaces can be optimised with different genetic distance metrics and the convergence of results needs to be assessed via pairwise matrix correlations. Based on the results presented here, high correlation coefficients across different genetic distance categories likely indicate accurate inference of true landscape resistance. Most importantly, empirical results should be interpreted with great caution, especially when they appear counter-intuitive in light of the ecology of a species.

Reference genome for the American rubyspot damselfly, Hetaerina americana

Damselflies and dragonflies (Order: Odonata) play important roles in both aquatic and terrestrial food webs and can serve as sentinels of ecosystem health and predictors of population trends in other taxa. The habitat requirements and limited dispersal of lotic damselflies make them especially sensitive to habitat loss and fragmentation. As such, landscape genomic studies of these taxa can help focus conservation efforts on watersheds with high levels of genetic diversity, local adaptation, and even cryptic endemism. Here, as part of the California Conservation Genomics Project (CCGP), we report the first reference genome for the American rubyspot damselfly, Hetaerina americana, a species associated with springs, streams and rivers throughout California. Following the CCGP assembly pipeline, we produced two de novo genome assemblies. The primary assembly includes 1,630,044,487 base pairs, with a contig N50 of 5.4 Mb, a scaffold N50 of 86.2 Mb, and a BUSCO completeness score of 97.6%. This is the seventh Odonata genome to be made publicly available and the first for the subfamily Hetaerininae. This reference genome fills an important phylogenetic gap in our understanding of Odonata genome evolution, and provides a genomic resource for a host of interesting ecological, evolutionary, and conservation questions for which the rubyspot damselfly genus Hetaerina is an important model system.

A brief history of population genetic research in California and an evaluation of its utility for conservation decision-making

A recently published macrogenetic dataset of California's flora and fauna, CaliPopGen, comprehensively summarizes population genetic research published between 1985-2020. Integrating these genetic data into the requisite "best available science" upon which conservation professionals rely should facilitate the prioritization of populations based on genetic health. We evaluate the extent to which the CaliPopGen Dataset provides genetic diversity estimates that are 1) unbiased, 2) sufficient in quantity, 3) cover entire species' ranges, and 4) include potentially adaptive loci. We identified genetic diversity estimates for 4,462 spatially-referenced populations of 432 species, confirming California's rich published history of population genetics research. Most recent studies used microsatellites markers, which have uniquely high levels of variation, and estimates of all genetic metrics varied significantly across marker types. Most studies used less than 10 loci for inferences, rendering parameter estimates potentially unreliable, and covered small spatial extents that include only a fraction of the studied species' California distribution (median 16.3%). In contrast, the ongoing California Conservation Genomics Project (CCGP) aims to cover the full geographical and environmental breadth of each species' occupied habitats, and uses a consistent approach based on whole-genome data. However, the CCGP will sequence only 12% of the number of individuals, and covers only about half the evolutionary diversity, of the CaliPopGen Database. There is clearly a place in the evaluation of the genetic health of California for both approaches going forward, especially if differences among studies can be minimized, and overlap emphasized. A complementary use of both datasets is warranted to inform optimal conservation decision-making.

The phylogeny of California, and how it informs setting multi-species conservation priorities

Incorporating measures of taxonomic diversity into research and management plans has long been a tenet of conservation science. Increasingly, active conservation programs are turning towards multi-species landscape and regional conservation actions, and away from single species approaches. This is both a reflection of changing trends in conservation science and advances in foundational technologies, including genomics and geospatial science. Multi-species approaches may provide more fundamental insights into evolutionary processes and equip managers with a more holistic understanding of the landscapes under their jurisdiction. Central to this approach are data generation and analyses which embrace and reflect a broad range of taxonomic diversity. Here we examine the family-level phylogenetic breadth of the California Conservation Genomics Project (CCGP) based on family-level phylogenetic diversity, family-level phylogenetic distinctness, and family richness. We place this in the context of the diversity present in California and compare it to the 35-plus years of genetic research compiled in the CaliPopGen Database. We found that the family-level phylogenetic diversity in the CCGP reflected that of California very well, slightly over-representing chordates and under-representing arthropods, and that 42% of CCGP phylogenetic diversity represented new contributions to genetic data for the state. In one focused effort, the CCGP was able to achieve roughly half the family-level phylogenetic diversity studied over the last several decades. To maximize studied phylogenetic diversity, future work should focus on arthropods, a conclusion that likely reflects the overall lack of attention to this hyper diverse clade.

Beyond the landscape: Resistance modelling infers physical and behavioural gene flow barriers to a mobile carnivore across a metropolitan area

Urbanization affects key aspects of wildlife ecology. Dispersal in urban wildlife species may be impacted by geographical barriers but also by a species' inherent behavioural variability. There are no functional connectivity analyses using continuous individual-based sampling across an urban-rural continuum that would allow a thorough assessment of the relative importance of physical and behavioural dispersal barriers. We used 16 microsatellite loci to genotype 374 red foxes (Vulpes vulpes) from the city of Berlin and surrounding rural regions in Brandenburg in order to study genetic structure and dispersal behaviour of a mobile carnivore across the urban-rural landscape. We assessed functional connectivity by applying an individual-based landscape genetic optimization procedure. Three commonly used genetic distance measures yielded different model selection results, with only the results of an eigenvector-based multivariate analysis reasonably explaining genetic differentiation patterns. Genetic clustering methods and landscape resistance modelling supported the presence of an urban population with reduced dispersal across the city border. Artificial structures (railways, motorways) served as main dispersal corridors within the cityscape, yet urban foxes avoided densely built-up areas. We show that despite their ubiquitous presence in urban areas, their mobility and behavioural plasticity, foxes were affected in their dispersal by anthropogenic presence. Distinguishing between man-made structures and sites of human activity, rather than between natural and artificial structures, is thus essential for better understanding urban fox dispersal. This differentiation may also help to understand dispersal of other urban wildlife and to predict how behaviour can shape population genetic structure beyond physical barriers.

Admixture of hybrid swarms of native and introduced lizards in cities is determined by the cityscape structure and invasion history

Introductions of non-native lineages increase opportunities for hybridization. Non-native lineages of the common wall lizard, Podarcis muralis, are frequently introduced in cities where they hybridize with native populations. We aimed at unravelling the invasion history and admixture of native and non-native wall lizards in four German cities using citywide, comprehensive sampling. We barcoded and genotyped 826 lizards and tested if gene flow in populations composed of admixed native and introduced lineages is facilitated by similar environmental factors to those in native populations by comparing fine-scale landscape genetic patterns. In cities with non-native lineages, lizards commonly occurred in numerous clusters of hybrid swarms, which showed variable lineage composition, consisting of up to four distinct evolutionary lineages. Hybrid swarms held vast genetic diversity and showed recent admixture with other hybrid swarms. Landscape genetic analyses showed differential effects of cityscape structures across cities, but identified water bodies as strong barriers to gene flow in both native and admixed populations. By contrast, railway tracks facilitated gene flow of admixed populations only. Our study shows that cities represent unique settings for hybridization, caused by multiple introductions of non-native taxa. Cityscape structure and invasion histories of cities will determine future evolutionary pathways at these novel hybrid zones.

Description of Apistogramma paulmuelleri sp. n., a new geophagine cichlid species (Teleostei: Perciformes) from the Amazon river basin in Loreto, Peru

A new species of Apistogramma is described from Peru, based on a total of 28 specimens collected in a small forest stream in the catchment of a nameless tributary of the Rio Amazonas about 80 kilometres south of Iquitos, Departamento Loreto (approximately 73° 34′ W / 04°24′ S). At first sight Apistogramma paulmuelleri sp. n. resembles A. regani, but is differentiated from the latter and all other Apistogramma species by the combination of a large band-like spot on the caudal-fin base, four distinct abdominal stripes, a roundish, banded caudal fin, a low dorsal fin without any striking features, in adult males yellow on the chin between the gill covers and on parts of the branchiostegal membranes, ivory sides to the head below the cheeks, and reversal of the band pattern during aggression and courtship display. Apistogramma paulmuelleri sp. n. is currently thought to be a representative of the Apistogramma eunotus complex within the Apistogramma regani lineage.

Female fitness optimum at intermediate mating rates under traumatic mating

Traumatic mating behaviors often bear signatures of sexual conflict and are then typically considered a male strategy to circumvent female choice mechanisms. In an extravagant mating ritual, the hermaphroditic sea slug Siphopteron quadrispinosum pierces the integument of their mating partners with a syringe-like penile stylet that injects prostate fluids. Traumatic injection is followed by the insertion of a spiny penis into the partner’s gonopore to transfer sperm. Despite traumatic mating, field mating rates exceed those required for female fertilization insurance, possibly because costs imposed on females are balanced by direct or indirect benefits of multiple sperm receipt. To test this idea, we exposed animals to a relevant range of mating opportunity regimes and assessed the effects on mating behavior and proxies of female fitness. We find penis intromission duration to decrease with mating rates, and a female fecundity maximum at intermediate mating rates. The latter finding indicates that benefits beyond fertilization insurance can make higher mating rates also beneficial from a female perspective in this traumatically mating species.

Apistogramma playayacu sp. n.: Description of a new cichlid species (Teleostei: Perciformes: Geophaginae) from the Rio Napo system, Ecuador

Apistogramma playayacu sp. n. is described based on a total of 41 specimens from tributaries of the Rio Napo from the border area between Ecuador and Peru (about 76°36′ W // 00°22′ S). The species is distinguished from all other Apistogramma species by the combination of 7 to 8 horizontal rows of mostly triangular black spots along the sides, horizontally devided spot on caudal base, rounded caudal fin, and straight narrow cheek stripe. The species is a moderately large (males up to 49 mm, females to 38 mm SL), high-backed, laterally compressed, and only slightly elongated Apistogramma species. Apistogramma playayacu sp. n. probably belongs to the Apistogramma macmasteri complex within the Apistogramma regani lineage (sensu Römer, 2006), and exhibits pronounced sexual dimorphism.