"The Reference Genome Of The Kidnapper Ant, Polyergus Mexicanus".

Polyergus kidnapper ants are widely distributed, but relatively uncommon, throughout the Holarctic, spanning an elevational range from sea level to over 3000 m. These species are well known for their obligate social parasitism with various Formica ant species, which they kidnap in dramatic, highly coordinated raids. Kidnapped Formica larvae and pupae become integrated into the Polyergus colony where they develop into adults and perform nearly all of the necessary colony tasks for the benefit of their captors. In California, Polyergus mexicanus is the most widely distributed Polyergus, but recent evidence has identified substantial genetic polymorphism within this species, including genetically divergent lineages associated with the use of different Formica host species. Given its unique behavior and genetic diversity, Polyergus mexicanus plays a critical role in maintaining ecosystem balance by influencing the population dynamics and genetic diversity of its host ant species, Formica, highlighting its conservation value and importance in the context of biodiversity preservation. Here, we present a high-quality genome assembly of P. mexicanus from a sample collected in Plumas County, CA, USA, in the foothills of the central Sierra Nevada. This genome assembly consists of 364 scaffolds spanning 252.31 Mb, with contig N50 of 481,250 kb, scaffold N50 of 10.36 Mb, and BUSCO completeness of 95.4%. We also assembled the genome of the Wolbachia endosymbiont of P. mexicanus - a single, circular contig spanning 1.23 Mb. These genome sequences provide essential resources for future studies of conservation genetics, population genetics, speciation, and behavioral ecology in this charismatic social insect.

Reversing the decline of threatened koala (Phascolarctos cinereus) populations in New South Wales: Using genomics to enhance conservation outcomes.

Genetic management is a critical component of threatened species conservation. Understanding spatial patterns of genetic diversity is essential for evaluating the resilience of fragmented populations to accelerating anthropogenic threats. Nowhere is this more relevant than on the Australian continent, which is experiencing an ongoing loss of biodiversity that exceeds any other developed nation. Using a proprietary genome complexity reduction-based method (DArTSeq), we generated a data set of 3239 high quality Single Nucleotide Polymorphisms (SNPs) to investigate spatial patterns and indices of genetic diversity in the koala (Phascolarctos cinereus), a highly specialised folivorous marsupial that is experiencing rapid and widespread population declines across much of its former range. Our findings demonstrate that current management divisions across the state of New South Wales (NSW) do not fully represent the distribution of genetic diversity among extant koala populations, and that care must be taken to ensure that translocation paradigms based on these frameworks do not inadvertently restrict gene flow between populations and regions that were historically interconnected. We also recommend that koala populations should be prioritised for conservation action based on the scale and severity of the threatening processes that they are currently faced with, rather than placing too much emphasis on their perceived value (e.g., as reservoirs of potentially adaptive alleles), as our data indicate that existing genetic variation in koalas is primarily partitioned among individual animals. As such, the extirpation of koalas from any part of their range represents a potentially critical reduction of genetic diversity for this iconic Australian species.

Evolutionary drivers of reproductive fitness in two endangered forest trees.

Population genetics theory predicts a relationship between fitness, genetic diversity (H0) and effective population size (Ne), which is often tested through heterozygosity-fitness correlations (HFCs). We tested whether population and individual fertility and heterozygosity are correlated in two endangered Mexican spruces (Picea martinezii and Picea mexicana) by combining genomic, demographic and reproductive data (seed development and germination traits). For both species, there was a positive correlation between population size and seed development traits, but not germination rate. Individual genome-wide heterozygosity and seed traits were only correlated in P. martinezii (general-effects HFC), and none of the candidate single nucleotide polymorphisms (SNPs) associated with individual fertility showed heterozygote advantage in any species (no local-effects HFC). We observed a single and recent (c. 30 thousand years ago (ka)) population decline for P. martinezii; the collapse of P. mexicana occurred in two phases separated by a long period of stability (c. 800 ka). Recruitment always contributed more to total population census than adult trees in P. mexicana, while this was only the case in the largest populations of P. martinezii. Equating fitness to either H0 or Ne, as traditionally proposed in conservation biology, might not always be adequate, as species-specific evolutionary factors can decouple the expected correlation between these parameters.

De novo genome assembly of a Geomyid rodent, Botta's pocket gopher (Thomomys bottae bottae).

Botta's pocket gopher (Thomomys bottae) is a common and widespread subterranean rodent of the North American west. The species has been of long interest to evolutionary biologists due to the phenotypic diversity across its range and unusual levels of variation in chromosome number and composition. Here, we present a high-quality reference genome from a male T. b. bottae individual captured in the San Francisco Bay Area. The assembly is comprised of 2,792 scaffolds, with a scaffold N50 value of 23.6 Mb and a BUSCO score of 91.0%. This genome helps fill a significant taxonomic sampling gap in rodent genome resources. With this reference genome, we envision new opportunities to investigate questions regarding the genomics of adaptation to the belowground niche. Further, we can begin to explore the impact of associated life history traits, such as limited dispersal and low population connectivity, on intraspecific genetic and phenotypic variation, genome evolution, speciation, and phylogenetic relationships across the Geomyoidea.

Beyond population size: whole-genome data reveal bottleneck legacies in the peninsular Italian wolf.

Preserving genetic diversity and adaptive potential while avoiding inbreeding depression is crucial for the long-term conservation of natural populations. Despite demographic increases, traces of past bottleneck events at the genomic level should be carefully considered for population management. From this perspective, the peninsular Italian wolf is a paradigmatic case. After being on the brink of extinction in the late 1960s, peninsular Italian wolves rebounded and recolonized most of the peninsula aided by conservation measures, including habitat and legal protection. Notwithstanding their demographic recovery, a comprehensive understanding of the genomic consequences of the historical bottleneck in Italian wolves is still lacking. To fill this gap, we sequenced whole genomes of thirteen individuals sampled in the core historical range of the species in Central Italy to conduct population genomic analyses, including a comparison with wolves from two highly-inbred wolf populations (i.e., Scandinavia and Isle Royale). We found that peninsular Italian wolves, despite their recent recovery, still exhibit relatively low genetic diversity, a small effective population size, signatures of inbreeding, and a non-negligible genetic load. Our findings indicate that the peninsular Italian wolf population is still susceptible to bottleneck legacies, which could lead to local inbreeding depression in case of population reduction or fragmentations. This study emphasizes the importance of considering key genetic parameters to design appropriate long-term conservation management plans.

Influence of dams on sauger population structure and hybridization with introduced walleye.

Dams have negatively affected freshwater biodiversity throughout the world. These negative effects tend to be exacerbated for aquatic taxa with migratory life histories, and for taxa whose habitat is fundamentally altered by the formation of large reservoirs. Sauger (Sander candadensis; Percidae), large-bodied migratory fishes native to North America, have seen population declines over much of the species' range, and dams are often implicated for their role in blocking access to spawning habitat and otherwise negatively affecting river habitat. Furthermore, hybridization appears to be more frequent between sauger and walleye in the reservoirs formed by large dams. In this study, we examine the role of dams in altering sauger population connectivity and facilitating hybridization with introduced walleye in Wyoming's Wind River and Bighorn River systems. We collected genomic data from individuals sampled over a large spatial scale and replicated sampling throughout the spawning season, with the intent to capture potential variation in hybridization prevalence or genomic divergence between sauger with different life histories. The timing of sampling was not related to hybridization prevalence or population divergence, suggesting limited genetic differences between sauger spawning in different time and places. Overall, there was limited hybridization detected, however, hybridization was most prevalent in Boysen Reservoir (a large impounded section of the Wind River). Dams in the lower Wind River and upper Bighorn River were associated with population divergence between sauger upstream and downstream of the dams, and demographic models suggest that this divergence has occurred in concordance with the construction of the dam. Sauger upstream of the dams exhibited substantially lower estimates of genetic diversity, which implies that disrupted connectivity between Wind River and Bighorn River sauger populations may already be causing negative demographic effects. This research points towards the importance of considering the evolutionary consequences of dams on fish populations in addition to the threats they pose to population persistence.

Chromosome-Level Reference Genome of the Ponza Grayling (Hipparchia sbordonii), an Italian Endemic and Endangered Butterfly.

Islands are crucial evolutionary hotspots, providing unique opportunities for differentiation of novel biodiversity and long-term segregation of endemic species. Islands are also fragile ecosystems, where biodiversity is more exposed to environmental and anthropogenic pressures than on continents. The Ponza grayling, Hipparchia sbordonii, is an endemic butterfly species that is currently found only in two tiny islands of the Pontine archipelago, off the coast of Italy, occupying an area smaller than 10 km2. It has been classified as Endangered (IUCN) because of the extremely limited area of occurrence, population fragmentation, and the recent demographic decline. Thanks to a combination of different assemblers of long and short genomic reads, bulk transcriptome RNAseq, and synteny analysis with phylogenetically close butterflies, we produced a highly contiguous, chromosome-scale annotated reference genome for the Ponza grayling, including 28 autosomes and the Z sexual chromosomes. The final assembly spanned 388.61 Gb with a contig N50 of 14.5 Mb and a BUSCO completeness score of 98.5%. Synteny analysis using four other butterfly species revealed high collinearity with Hipparchia semele and highlighted 10 intrachromosomal inversions longer than 10 kb, of which two appeared on the lineage leading to H. sbordonii. Our results show that a chromosome-scale reference genome is attainable also when chromatin conformation data may be impractical or present specific technical challenges. The high-quality genomic resource for H. sbordonii opens up new opportunities for the accurate assessment of genetic diversity and genetic load and for the investigations of the genomic novelties characterizing the evolutionary path of this endemic island species.

Clinging on the brink: Whole genomes reveal human-induced population declines and severe inbreeding in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus).

Human-mediated habitat destruction has had a profound impact on increased species extinction rates and population declines worldwide. The coastal development in the United Arab Emirates (UAE) over the last two decades, serves as an example of how habitat transformation can alter the landscape of a country in just a few years. Here, we study the genomic implications of habitat transformation in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus), the only endemic vertebrate of the UAE. We generate a high-quality reference genome for this gecko, representing the first reference genome for the family Phyllodactylidae, and produce whole-genome resequencing data for 23 specimens from 10 different species of leaf-toed geckos. Our results show that A. caudivolvulus has consistently lower genetic diversity than any other Arabian species of Asaccus, suggesting a history of ancient population declines. However, high levels of recent inbreeding are recorded among populations in heavily developed areas, with a more than 50% increase in long runs of homozygosity within a 9-year period. Moreover, results suggest that this species does not effectively purge deleterious mutations, hence making it more vulnerable to future stochastic threats. Overall, results show that A. caudivolvulus is in urgent need of protection, and habitat preservation must be warranted to ensure the species' survival.

Conservation genomics of the critically endangered Chinese pangolin.

The Chinese pangolin (Manis pentadactyla, MP) has been extensively exploited and is now on the brink of extinction, but its population structure, evolutionary history, and adaptive potential are unclear. Here, we analyzed 94 genomes from three subspecies of the Chinese pangolin and identified three distinct genetic clusters (MPA, MPB, and MPC), with MPB further divided into MPB1 and MPB2 subpopulations. The divergence of these populations was driven by past climate change. For MPB2 and MPC, recent human activities have caused dramatic population decline and small population size as well as increased inbreeding, but not decrease in genomic variation and increase in genetic load probably due to strong gene flow; therefore, it is crucial to strengthen in situ habitat management for these two populations. By contrast, although human activities have a milder impact on MPA, it is at high risk of extinction due to long-term contraction and isolation, and genetic rescue is urgently needed. MPB1 exhibited a relatively healthy population status and can potentially serve as a source population. Overall, our findings provide novel insights into the conservation of the Chinese pangolin and biogeography of the mammals of eastern Asia.

A genome assembly of the American black bear, Ursus americanus, from California.

The American black bear, Ursus americanus, is a widespread and ecologically important species in North America. In California, the black bear plays an important role in a variety of ecosystems and serves as an important species for recreational hunting. While research suggests that the populations in California are currently healthy, continued monitoring is critical, with genomic analyses providing an important surveillance tool. Here we report a high-quality, near chromosome-level genome assembly from a U. americanus sample from California. The primary assembly has a total length of 2.5 Gb contained in 316 scaffolds, a contig N50 of 58.9 Mb, a scaffold N50 of 67.6 Mb, and a BUSCO completeness score of 96%. This U. americanus genome assembly will provide an important resource for the targeted management of black bear populations in California, with the goal of achieving an appropriate balance between the recreational value of black bears and the maintenance of viable populations. The high quality of this genome assembly will also make it a valuable resource for comparative genomic analyses among black bear populations and among bear species.

Exploring phylogeny and genomic vulnerability of Melastoma (Melastomataceae) endemic to a World Natural Heritage site, the Bonin Islands.

Melastoma (Melastomataceae) includes ca. 100 species across tropical to subtropical regions of Asia and Oceania. The Bonin Islands harbor three endemic taxa: M. tetramerum, M. tetramerum var. pentapetalum, and M. candidum var. alessandrense. Of these, M. tetramerum is critically endangered and faces near extinction in the wild. This study investigates the phylogenetic relationships among these endemic Melastoma species in the Bonin Islands based on the whole chloroplast genome and nuclear SNPs. The results revealed that M. candidum var. alessandrense was placed in the clade of the widespread East Asian M. candidum and has a distinct evolutionary origin from the other two taxa. The population genomics analyses (heterozygosity, rates of deleterious mutations, and numbers and lengths of runs of homozygosity) indicated lower genetic diversity and more vulnerable genomes of endemic Melastoma, especially M. tetramerum var. pentapetalum. M. tetramerum var. pentapetalum is not a target of any protection programs, however, conservation plans might be required for this variety because M. tetramerum var. pentapetalum would have a more vulnerable genome than M. tetramerum, which faces near extinction in the wild. This information can facilitate the development of effective conservation strategies in a precautionary way that anticipates imminent threats to the survival of the species.

Bringing genomics to the field: An integrative approach to seed sourcing for forest restoration.

Premise: Global anthropogenic change threatens the health and productivity of forest ecosystems. Assisted migration and reforestation are tools to help mitigate these impacts. However, questions remain about how to approach sourcing seeds to ensure high establishment and future adaptability. Methods: Using exome-capture sequencing, we demonstrate a computational approach to finding the best n-sets from a candidate list of seed sources that collectively achieve high genetic diversity (GD) and minimal genetic load (GL), while also increasing evolvability in quantitative traits. The benefits of this three-part strategy (diversity-load-evolvability) are to increase near-term establishment success while also boosting evolutionary potential to respond to future stressors. Members of The Nature Conservancy and the Central Appalachian Spruce Restoration Initiative planted 58,000 seedlings across 255 acres. A subset of seedlings was monitored for establishment success and variation in growth. Results: The results show gains in GD relative to GL and increases in quantitative genetic variation in seedling growth for pooled vs. single-source restoration. No single "super source" was observed across planting sites; rather, monitoring results demonstrate that pooling of multiple sources helps achieve higher GD:GL and evolvability. Discussion: Our study shows the potential for integrating genomics into local-scale restoration and the importance of building partnerships between academic researchers and applied conservation managers.

Local adaptation of Pinus leiophylla under climate and land use change models in the Avocado Belt of Michoacán.

Climate change and land use change are two main drivers of global biodiversity decline, decreasing the genetic diversity that populations harbour and altering patterns of local adaptation. Landscape genomics allows measuring the effect of these anthropogenic disturbances on the adaptation of populations. However, both factors have rarely been considered simultaneously. Based on a set of 3660 SNPs from which 130 were identified as outliers by a genome-environment association analysis (LFMM), we modelled the spatial turnover of allele frequencies in 19 localities of Pinus leiophylla across the Avocado Belt in Michoacán state, Mexico. Then, we evaluated the effect of climate change and land use change scenarios, in addition to evaluating assisted gene flow strategies and connectivity metrics across the landscape to identify priority conservation areas for the species. We found that localities in the centre-east of the Avocado Belt would be more vulnerable to climate change, while localities in the western area are more threatened by land conversion to avocado orchards. Assisted gene flow actions could aid in mitigating both threats. Connectivity patterns among forest patches will also be modified by future habitat loss, with central and eastern parts of the Avocado Belt maintaining the highest connectivity. These results suggest that areas with the highest priority for conservation are in the eastern part of the Avocado Belt, including the Monarch Butterfly Biosphere Reserve. This work is useful as a framework that incorporates distinct layers of information to provide a more robust representation of the response of tree populations to anthropogenic disturbances.

Conservation genomics provides insights into genetic resilience and adaptation of the endangered Chinese hazelnut, Corylus chinensis.

Global climate change has increased concerns regarding biodiversity loss. However, many key conservation issues still required further research, including demographic history, deleterious mutation load, adaptive evolution, and putative introgression. Here we generated the first chromosome-level genome of the endangered Chinese hazelnut, Corylus chinensis, and compared the genomic signatures with its sympatric widespread C. kwechowensis-C. yunnanensis complex. We found large genome rearrangements across all Corylus species and identified species-specific expanded gene families that may be involved in adaptation. Population genomics revealed that both C. chinensis and the C. kwechowensis-C. yunnanensis complex had diverged into two genetic lineages, forming a consistent pattern of southwestern-northern differentiation. Population size of the narrow southwestern lineages of both species have decreased continuously since the late Miocene, whereas the widespread northern lineages have remained stable (C. chinensis) or have even recovered from population bottlenecks (C. kwechowensis-C. yunnanensis complex) during the Quaternary. Compared with C. kwechowensis-C. yunnanensis complex, C. chinensis showed significantly lower genomic diversity and higher inbreeding level. However, C. chinensis carried significantly fewer deleterious mutations than C. kwechowensis-C. yunnanensis complex, as more effective purging selection reduced the accumulation of homozygous variants. We also detected signals of positive selection and adaptive introgression in different lineages, which facilitated the accumulation of favorable variants and formation of local adaptation. Hence, both types of selection and exogenous introgression could have mitigated inbreeding and facilitated survival and persistence of C. chinensis. Overall, our study provides critical insights into lineage differentiation, local adaptation, and the potential for future recovery of endangered trees.

Estimation of contemporary effective population size in plant populations: Limitations of genomic datasets.

Effective population size (N e) is a pivotal evolutionary parameter with crucial implications in conservation practice and policy. Genetic methods to estimate N e have been preferred over demographic methods because they rely on genetic data rather than time-consuming ecological monitoring. Methods based on linkage disequilibrium (LD), in particular, have become popular in conservation as they require a single sampling and provide estimates that refer to recent generations. A software program based on the LD method, GONE, looks particularly promising to estimate contemporary and recent-historical N e (up to 200 generations in the past). Genomic datasets from non-model species, especially plants, may present some constraints to the use of GONE, as linkage maps and reference genomes are seldom available, and SNP genotyping is usually based on reduced-representation methods. In this study, we use empirical datasets from four plant species to explore the limitations of plant genomic datasets when estimating N e using the algorithm implemented in GONE, in addition to exploring some typical biological limitations that may affect N e estimation using the LD method, such as the occurrence of population structure. We show how accuracy and precision of N e estimates potentially change with the following factors: occurrence of missing data, limited number of SNPs/individuals sampled, and lack of information about the location of SNPs on chromosomes, with the latter producing a significant bias, previously unexplored with empirical data. We finally compare the N e estimates obtained with GONE for the last generations with the contemporary N e estimates obtained with the programs currentNe and NeEstimator.

Development and validation of a GT-seq panel for genetic monitoring in a threatened species using minimally invasive sampling.

Minimally invasive samples are often the best option for collecting genetic material from species of conservation concern, but they perform poorly in many genomic sequencing methods due to their tendency to yield low DNA quality and quantity. Genotyping-in-thousands by sequencing (GT-seq) is a powerful amplicon sequencing method that can genotype large numbers of variable-quality samples at a standardized set of single nucleotide polymorphism (SNP) loci. Here, we develop, optimize, and validate a GT-seq panel for the federally threatened northern Idaho ground squirrel (Urocitellus brunneus) to provide a standardized approach for future genetic monitoring and assessment of recovery goals using minimally invasive samples. The optimized panel consists of 224 neutral and 81 putatively adaptive SNPs. DNA collected from buccal swabs from 2016 to 2020 had 73% genotyping success, while samples collected from hair from 2002 to 2006 had little to no DNA remaining and did not genotype successfully. We evaluated our GT-seq panel by measuring genotype discordance rates compared to RADseq and whole-genome sequencing. GT-seq and other sequencing methods had similar population diversity and F ST estimates, but GT-seq consistently called more heterozygotes than expected, resulting in negative F IS values at the population level. Genetic ancestry assignment was consistent when estimated with different sequencing methods and numbers of loci. Our GT-seq panel is an effective and efficient genotyping tool that will aid in the monitoring and recovery of this threatened species, and our results provide insights for applying GT-seq for minimally invasive DNA sampling techniques in other rare animals.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery.

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

Museum genomics approach to study the taxonomy and evolution of Woolly-necked storks using historic specimens.

The accessibility of genomic tools in evolutionary biology has allowed for a thorough exploration of various evolutionary processes associated with adaptation and speciation. However, genomic studies in natural systems present numerous challenges, reflecting the inherent complexities of studying organisms in their native habitats. The utilization of museum specimens for genomics research has received increased attention in recent times, facilitated by advancements in ancient DNA techniques. In this study, we have utilized a museum genomics approach to analyze historic specimens of Woolly-necked storks (Ciconia spp.) and examine their genetic composition and taxonomic status and explore the evolutionary and adaptive trajectories of populations over the years. The Woolly-necked storks are distributed in Asia and Africa with a taxonomic classification that has been a matter of ambiguity. Asian and African Woollynecks were recently recognized as different species based on their morphological differences; however, their genomic validation was lacking. In this study, we have used ∼70-year-old museum samples for whole-genome population-scale sequencing. Our study has revealed that Asian and African Woollynecks are genetically distinct, consistent with the current taxonomic classification based on morphological features. However, we also found a high genetic divergence between the Asian subspecies Ciconia episcopus neglecta and Ciconia episcopus episcopus, suggesting this classification requires a detailed examination to explore processes of ongoing speciation. Because taxonomic classification directly impacts conservation efforts, and there is evidence of declining populations of Asian Woollynecks in Southeast Asia, our results highlight that population-scale studies are urgent to determine the genetic, ecological, and phylogenetic diversity of these birds.

A draft reference genome assembly of California Pipevine, Aristolochia californica Torr.

The California Pipevine, Aristolochia californica Torr., is the only endemic California species within the cosmopolitan birthwort family Aristolochiaceae. It occurs as an understory vine in riparian and chaparral areas and in forest edges and windrows. The geographic range of this plant species almost entirely overlaps with that of its major specialized herbivore, the California Pipevine Swallowtail Butterfly Battus philenor hirsuta. While this species pair is a useful, ecologically well-understood system to study co-evolution, until recently, genomic resources for both have been lacking. Here, we report a new, chromosome-level assembly of A. californica as part of the California Conservation Genomics Project (CCGP). Following the sequencing and assembly strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin proximity sequencing technology to produce a de novo assembled genome. Our genome assembly, the first for any species in the genus, contains 531 scaffolds spanning 661 megabase (Mb) pairs, with a contig N50 of 6.53 Mb, a scaffold N50 of 42.2 Mb, and BUSCO complete score of 98%. In combination with the recently published B. philenor hirsuta reference genome assembly, the A. californica reference genome assembly will be a powerful tool for studying co-evolution in a rapidly changing California landscape.

Variation in genomic vulnerability to climate change across temperate populations of eelgrass (Zostera marina).

A global decline in seagrass populations has led to renewed calls for their conservation as important providers of biogenic and foraging habitat, shoreline stabilization and carbon storage. Eelgrass (Zostera marina) occupies the largest geographic range among seagrass species spanning a commensurately broad spectrum of environmental conditions. In Canada, eelgrass is managed as a single phylogroup despite occurring across three oceans and a range of ocean temperatures and salinity gradients. Previous research has focused on applying relatively few markers to reveal population structure of eelgrass, whereas a whole-genome approach is warranted to investigate cryptic structure among populations inhabiting different ocean basins and localized environmental conditions. We used a pooled whole-genome re-sequencing approach to characterize population structure, gene flow and environmental associations of 23 eelgrass populations ranging from the Northeast United States to Atlantic, subarctic and Pacific Canada. We identified over 500,000 SNPs, which when mapped to a chromosome-level genome assembly revealed six broad clades of eelgrass across the study area, with pairwise F ST ranging from 0 among neighbouring populations to 0.54 between Pacific and Atlantic coasts. Genetic diversity was highest in the Pacific and lowest in the subarctic, consistent with colonization of the Arctic and Atlantic oceans from the Pacific less than 300 kya. Using redundancy analyses and two climate change projection scenarios, we found that subarctic populations are predicted to be potentially more vulnerable to climate change through genomic offset predictions. Conservation planning in Canada should thus ensure that representative populations from each identified clade are included within a national network so that latent genetic diversity is protected, and gene flow is maintained. Northern populations, in particular, may require additional mitigation measures given their potential susceptibility to a rapidly changing climate.