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 and 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. Finally, a synopsis of the available population genetic data for California, all other US states and 241 other countries, allows us to identify states and countries for which meaningful data summaries, such as CaliPopGen, could be collated and others, which have limited published data available and are prime targets for future, empirical work.

The phylogeny of California, and how it informs setting multispecies 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 toward multispecies 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. Multispecies 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 (PD), 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 PD in the CCGP reflected that of California very well, slightly overrepresenting chordates and underrepresenting arthropods, and that 42% of CCGP PD represented new contributions to genetic data for the state. In one focused effort, the CCGP was able to achieve roughly half the family-level PD studied over the last several decades. To maximize studied PD, future work should focus on arthropods, a conclusion that likely reflects the overall lack of attention to this hyperdiverse 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.

Cityscape genetics: structural vs. functional connectivity of an urban lizard population.

Functional connectivity is essential for the long-term persistence of populations. However, many studies assess connectivity with a focus on structural connectivity only. Cityscapes, namely urban landscapes, are particularly dynamic and include numerous potential anthropogenic barriers to animal movements, such as roads, traffic or buildings. To assess and compare structural connectivity of habitats and functional connectivity of gene flow of an urban lizard, we here combined species distribution models (SDMs) with an individual-based landscape genetic optimization procedure. The most important environmental factors of the SDMs are structural diversity and substrate type, with high and medium levels of structural diversity as well as open and rocky/gravel substrates contributing most to structural connectivity. By contrast, water cover was the best model of all environmental factors following landscape genetic optimization. The river is thus a major barrier to gene flow, while of the typical anthropogenic factors only buildings showed an effect. Nonetheless, using SDMs as a basis for landscape genetic optimization provided the highest ranked model for functional connectivity. Optimizing SDMs in this way can provide a sound basis for models of gene flow of the cityscape, and elsewhere, while presence-only and presence-absence modelling approaches showed differences in performance. Additionally, interpretation of results based on SDM factor importance can be misleading, dictating more thorough analyses following optimization of SDMs. Such approaches can be adopted for management strategies, for example aiming to connect native common wall lizard populations or disconnect them from non-native introduced populations, which are currently spreading in many cities in Central Europe.

Biodiversity in cities needs space: a meta-analysis of factors determining intra-urban biodiversity variation.

Understanding varying levels of biodiversity within cities is pivotal to protect it in the face of global urbanisation. In the early stages of urban ecology studies on intra-urban biodiversity focused on the urban-rural gradient, representing a broad generalisation of features of the urban landscape. Increasingly, studies classify the urban landscape in more detail, quantifying separately the effects of individual urban features on biodiversity levels. However, while separate factors influencing biodiversity variation among cities worldwide have recently been analysed, a global analysis on the factors influencing biodiversity levels within cities is still lacking. We here present the first meta-analysis on intra-urban biodiversity variation across a large variety of taxonomic groups of 75 cities worldwide. Our results show that patch area and corridors have the strongest positive effects on biodiversity, complemented by vegetation structure. Local, biotic and management habitat variables were significantly more important than landscape, abiotic or design variables. Large sites greater than 50 ha are necessary to prevent a rapid loss of area-sensitive species. This indicates that, despite positive impacts of biodiversity-friendly management, increasing the area of habitat patches and creating a network of corridors is the most important strategy to maintain high levels of urban biodiversity.

Sexual selection drives asymmetric introgression in wall lizards.

Hybridisation is increasingly recognised as an important cause of diversification and adaptation. Here, we show how divergence in male secondary sexual characters between two lineages of the common wall lizard (Podarcis muralis) gives rise to strong asymmetries in male competitive ability and mating success, resulting in asymmetric hybridisation upon secondary contact. Combined with no negative effects of hybridisation on survival or reproductive characters in F1-hybrids, these results suggest that introgression should be asymmetric, resulting in the displacement of sexual characters of the sub-dominant lineage. This prediction was confirmed in two types of secondary contact, across a natural contact zone and in two introduced populations. Our study illustrates how divergence in sexually selected traits via male competition can determine the direction and extent of introgression, contributing to geographic patterns of genetic and phenotypic diversity.

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.

CaliPopGen: A genetic and life history database for the fauna and flora of California.

CaliPopGen is a database of population genetic data for native and naturalized eukaryotic species in California, USA. It summarizes the published literature (1985-2020) for 5,453 unique populations with genetic data from more than 187,394 individuals and 448 species (513 species plus subspecies) across molecular markers including allozymes, RFLPs, mtDNA, microsatellites, nDNA, and SNPs. Terrestrial habitats accounted for the majority (46.4%) of the genetic data. Taxonomic groups with the greatest representation were Magnoliophyta (20.31%), Insecta (13.4%), and Actinopterygii (12.85%). CaliPopGen also reports life-history data for most included species to enable analyses of the drivers of genetic diversity across the state. The large number of populations and wide taxonomic breadth will facilitate explorations of ecological patterns and processes across the varied geography of California. CaliPopGen covers all terrestrial and marine ecoregions of California and has a greater density of species and georeferenced populations than any previously published population genetic database. It is thus uniquely suited to inform conservation management at the regional and state levels across taxonomic groups.

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.