Search

When invasive and endangered native taxa hybridize, the resulting admixture introduces novel conservation challenges. Across a large region of central California, a hybrid swarm consisting of admixed endangered California tiger salamanders (CTS) (Ambystoma californiense) and introduced barred tiger salamanders (BTS) (Ambystoma mavortium) has replaced native populations, threatening the genetic integrity of CTS and the vernal pool systems they inhabit. We employed a large‐scale, genomically informed field experiment to test whether shortening breeding pond hydroperiod would favor native CTS genotypes. We constructed 14 large, seminatural ponds to evaluate the effect of hydroperiod duration on larval survival and mass at metamorphosis. We tracked changes in non‐native allele frequencies with a 5237‐gene exon capture array and employed a combination of custom Bayesian and generalized linear models to quantify the effect of pond duration on salamander fitness. Earlier work on this system showed hybrid superiority under many conditions and suggested that hybrids are favored in human‐modified ponds with artificially long hydroperiods. Consistent with these earlier studies, we found overwhelming evidence for hybrid superiority. Very short hydroperiods substantially reduced the mass (1.1–1.5 fold) and survival probability (10–13 fold) of both native and hybrid larvae, confirming that hydroperiod likely exerts a strong selective pressure in the wild. We identified 86 genes, representing 1.8% of 4723 screened loci, that significantly responded to this hydroperiod‐driven selection. In contrast to earlier work, under our more natural experimental conditions, native CTS survival and size at metamorphosis were always less than hybrids, suggesting that hydroperiod management alone will not shift selection to favor native larval genotypes. However, shortening pond hydroperiod may limit productivity of hybrid ponds, complementing other strategies to remove hybrids while maintaining vernal pool ecosystems. This study confirms and expands on previous work that highlights the importance of hydroperiod management to control invasive aquatic species. [ABSTRACT FROM AUTHOR]

The California Conservation Genomics Project (CCGP) is a unique, critically important step forward in the use of comprehensive landscape genetic data to modernize natural resource management at a regional scale. We describe the CCGP, including all aspects of project administration, data collection, current progress, and future challenges. The CCGP will generate, analyze, and curate a single high-quality reference genome and 100–150 resequenced genomes for each of 153 species projects (representing 235 individual species) that span the ecological and phylogenetic breadth of California's marine, freshwater, and terrestrial ecosystems. The resulting portfolio of roughly 20 000 resequenced genomes will be analyzed with identical informatic and landscape genomic pipelines, providing a comprehensive overview of hotspots of within-species genomic diversity, potential and realized corridors connecting these hotspots, regions of reduced diversity requiring genetic rescue, and the distribution of variation critical for rapid climate adaptation. After 2 years of concerted effort, full funding ($12M USD) has been secured, species identified, and funds distributed to 68 laboratories and 114 investigators drawn from all 10 University of California campuses. The remaining phases of the CCGP include completion of data collection and analyses, and delivery of the resulting genomic data and inferences to state and federal regulatory agencies to help stabilize species declines. The aspirational goals of the CCGP are to identify geographic regions that are critical to long-term preservation of California biodiversity, prioritize those regions based on defensible genomic criteria, and provide foundational knowledge that informs management strategies at both the individual species and ecosystem levels. [ABSTRACT FROM AUTHOR]

Hybridization between native and non‐native species is an ongoing global conservation threat. Hybrids that exhibit traits and tolerances that surpass parental values are of particular concern, given their potential to outperform native species. Effective management of hybrid populations requires an understanding of both physiological performance and the underlying mechanisms that drive transgressive hybrid traits. Here, we explore several aspects of the hybridization between the endangered California tiger salamander (Ambystoma californiense; CTS) and the introduced barred tiger salamander (Ambystoma mavortium; BTS). We assayed critical thermal maximum (CTMax) to compare the ability of CTS, BTS and F1 hybrids to tolerate acute thermal stress, and found that hybrids exhibit a wide range of CTMax values, with 33% (4/12) able to tolerate temperatures greater than either parent. We then quantified the genomic response, measured at the RNA transcript level, of each salamander, to explore the mechanisms underlying thermal tolerance strategies. We found that CTS and BTS have strikingly different values and tissue‐specific patterns of overall gene expression, with hybrids expressing intermediate values. F1 hybrids display abundant and variable degrees of allele‐specific expression (ASE), likely arising from extensive compensatory evolution in gene regulatory mechanisms between CTS and BTS. We found evidence that the proportion of genes with allelic imbalance in individual hybrids correlates with their CTMax, suggesting a link between ASE and expanded thermal tolerance that may contribute to the success of hybrid salamanders in California. Future climate change may further complicate management of CTS if hybrid salamanders are better equipped to deal with rising temperatures. [ABSTRACT FROM AUTHOR]

Introductions of non‐native predators often reduce biodiversity and affect natural predator–prey relationships and may increase the abundance of potential disease vectors (e.g., mosquitoes) indirectly through competition or predation cascades. The Santa Monica Mountains (California, U.S.A.), situated in a global biodiversity hotspot, is an area of conservation concern due to climate change, urbanization, and the introduction of non‐native species. We examined the effect of non‐native crayfish (Procambarus clarkii) on an existing native predator, dragonfly nymphs (Aeshna sp.), and their mosquito larvae (Anopheles sp.) prey. We used laboratory experiments to compare the predation efficiency of both predators, separately and together, and field data on counts of dragonfly nymphs and mosquito larvae sampled from 13 local streams. We predicted a lower predation efficiency of crayfish compared with native dragonfly nymphs and a reduced predation efficiency of dragonfly nymphs in the presence of crayfish. Dragonfly nymphs were an order of magnitude more efficient predators than crayfish, and dragonfly nymph predation efficiency was reduced in the presence of crayfish. Field count data showed that populations of dragonfly nymphs and mosquito larvae were strongly correlated with crayfish presence in streams, such that sites with crayfish tended to have fewer dragonfly nymphs and more mosquito larvae. Under natural conditions, it is likely that crayfish reduce the abundance of dragonfly nymphs and their predation efficiency and thereby, directly and indirectly, lead to higher mosquito populations and a loss of ecosystem services related to disease vector control. Article impact statement: The efficiency of native predators to control potential disease vectors is diminished in the presence of non‐native predators. [ABSTRACT FROM AUTHOR]

Accurate time-calibrated phylogenies are the centerpiece of many macroevolutionary studies, and the relationship between the size and scale of molecular data sets and the density and accuracy of fossil calibrations is a key element of time tree studies. Here, we develop a target capture array specifically for living turtles, compare its efficiency to an ultraconserved element (UCE) dataset, and present a time-calibrated molecular phylogeny based on 539 nuclear loci sequenced from 26 species representing the breadth of living turtle diversity plus outgroups. Our gene array, based on three fully sequenced turtle genomes, is 2.4 times more variable across turtles than a recently published UCE data set for an identical subset of 13 species, confirming that taxon-specific arrays return more informative data per sequencing effort than UCEs. We used our genomic data to estimate the ages of living turtle clades including a mid-late Triassic origin for crown turtles and a mid-Carboniferous split of turtles from their sister group, Archosauria. By specifically excluding several of the earliest potential crown turtle fossils and limiting the age of fossil calibration points to the unambiguous crown lineage Caribemys oxfordiensis from the Late Jurassic (Oxfordian, 163.5–157.3 Ma) we corroborate a relatively ancient age for living turtles. We also provide novel age estimates for five of the ten testudine families containing more than a single species, as well as several intrafamilial clades. Most of the diversity of crown turtles appears to date to the Paleogene, well after the Cretaceous-Paleogene mass extinction 66 mya. [ABSTRACT FROM AUTHOR]

Wide‐ranging species are seldom considered conservation priorities, yet they have the potential to harbour genetically deeply differentiated units across environments or ecological barriers, including some that warrant taxonomic recognition. Documenting such cryptic genetic diversity is especially important for wide‐ranging species that are in decline, as they may comprise a set of even more endangered lineages or species with small distributions. However, studies of wide‐ranging species, particularly when they cross political borders, are extremely challenging. One approach to overcoming these challenges is to conduct detailed local analyses in combination with less detailed, range‐wide studies. We used this approach with the red‐footed tortoise (Chelonoidis carbonarius), a threatened species likely to contain cryptic diversity given its vast range and the distinctive ecoregions that it inhabits. Previous single‐gene molecular studies indicated the presence of at least five lineages, two of which occur in different ecoregions separated by the Andes within Colombia. We used a comprehensive genomic analysis to test the hypothesis of cryptic diversity within the single jurisdiction of Colombia. We used a combination of restriction‐site‐associated DNA sequencing and environmental niche modelling to provide three independent lines of evidence that support the presence of important cryptic diversity that may deserve taxonomic recognition: allopatric reproductive isolation, local adaptation and ecological divergence. We also provide a fine‐scale genetic map with the distribution of conservation units in Colombia. As we complete ongoing range‐wide analyses and make taxonomic adjustments, we recommend that the two lineages in Colombia be treated as separate units for conservation purposes. [ABSTRACT FROM AUTHOR]

The explosive growth of empirical population genetics has seen a proliferation of analytical methods leading to a steady increase in our ability to accurately measure key population parameters, including genetic isolation, effective population size, and gene flow, in natural systems. Assuming they yield similar results, population genetic methods offer an attractive complement to, or replacement of, traditional field-ecological studies. However, empirical assessments of the concordance between direct field-ecological and indirect population genetic studies of the same populations are uncommon in the literature. In this study, we investigate genetic isolation, rates of dispersal, and population sizes for the endangered California tiger salamander, Ambystoma californiense, across multiple breeding seasons in an intact vernal pool network. We then compare our molecular results to a previously published study based on multiyear, mark-recapture data from the same breeding sites. We found that field and genetic estimates of population size were only weakly correlated, but dispersal rates were remarkably congruent across studies and methods. In fact, dispersal probability functions derived from genetic data and traditional field-ecological data were a significant match, suggesting that either method can be used effectively to assess population connectivity. These results provide one of the first explicit tests of the correspondence between landscape genetic and field-ecological approaches to measuring functional population connectivity and suggest that even single-year genetic samples can return biologically meaningful estimates of natural dispersal and gene flow. [ABSTRACT FROM AUTHOR]

After decades of discussion, there is little consensus on the extent to which hybrids between endangered and nonendangered species should be protected by US law. As increasingly larger, genome-scale data sets are developed, we can identify individuals and populations with even trace levels of genetic admixture, making the 'hybrid problem' all the more difficult. We developed a decision-tree framework for evaluating hybrid protection, including both the processes that produced hybrids (human-mediated or natural) and the ecological impact of hybrids on natural ecosystems. We then evaluated our decision tree for four case studies drawn from our own work and briefly discuss several other cases from the literature. Throughout, we highlight the management outcomes that our approach provides and the nuances of hybridization as a conservation problem. [ABSTRACT FROM AUTHOR]

Defining species' niches is central to understanding their distributions and is thus fundamental to basic ecology and climate change projections. Ecological niche models (ENMs) are a key component of making accurate projections and include descriptions of the niche in terms of both response curves and rankings of variable importance. In this study, we evaluate Maxent's ranking of environmental variables based on their importance in delimiting species' range boundaries by asking whether these same variables also govern annual recruitment based on long-term demographic studies. We found that Maxent-based assessments of variable importance in setting range boundaries in the California tiger salamander (Ambystoma californiense; CTS) correlate very well with how important those variables are in governing ongoing recruitment of CTS at the population level. This strong correlation suggests that Maxent's ranking of variable importance captures biologically realistic assessments of factors governing population persistence. However, this result holds only when Maxent models are built using best-practice procedures and variables are ranked based on permutation importance. Our study highlights the need for building high-quality niche models and provides encouraging evidence that when such models are built, they can reflect important aspects of a species' ecology. [ABSTRACT FROM AUTHOR]

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems. [ABSTRACT FROM AUTHOR]

Many studies employ ecological niche models (ENMs) to predict species' occurrences in undersampled regions, generally without field confirmation. Here, we use field surveys to test the relative utility of four potential refinements to the standard ENM approach: 1) altering model complexity based on AICc, 2) selecting background points from a biologically informed region, 3) using target-group background to account for sampling bias in existing localities, and 4) using many rangewide localities (global model) versus fewer proximal localities (local model) to construct geographically restricted range predictions. We used Maxent to predict new localities for the California tiger salamander Ambystoma californiense, an endangered species that often goes undocumented due to its cryptic lifestyle. We followed this with a field survey of 260 previously unsampled potential breeding sites in Solano County, CA and used the resulting presence/absence data to compare all factorial combinations of the four model refinements using a new application of the Kruskal-Wallis test for ENM outputs. Our field surveys led to the discovery of 81 previously undocumented breeding localities for the California tiger salamander and demonstrated that ENMs could be significantly improved by utilizing target-group background to account for spatial sampling bias and local models to focus model output on the subregion of the range being surveyed. Our results clearly demonstrate the potential for local models to outperform global models, and we recommend supplementing traditional Maxent global models that utilize all known localities with local models, particularly when species occupy geographically structured, heterogeneous habitat types. We also recommend using target-group background since the improvement we observed when including it in our models was significant and very similar to that documented by previous studies. Most importantly, we emphasize the importance of field verification to enable rigorous statistical comparisons among models. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Background: Among the greatest challenges for biology in the 21st century is inference of the tree of life. Interest in, and progress toward, this goal has increased dramatically with the growing availability of molecular sequence data. However, we have very little sense, for any major clade, of how much progress has been made in resolving a full tree of life and the scope of work that remains. A series of challenges stand in the way of completing this task but, at the most basic level, progress is limited by data: a limited fraction of the world's biodiversity has been incorporated into a phylogenetic analysis. More troubling is our poor understanding of what fraction of the tree of life is understood and how quickly research is adding to this knowledge. Here we measure the rate of progress on the tree of life for one clade of particular research interest, the vertebrates. Results: Using an automated phylogenetic approach, we analyse all available molecular data for a large sample of vertebrate diversity, comprising nearly 12,000 species and 210,000 sequences. Our results indicate that progress has been rapid, increasing polynomially during the age of molecular systematics. It is also skewed, with birds and mammals receiving the most attention and marine organisms accumulating far fewer data and a slower rate of increase in phylogenetic resolution than terrestrial taxa. We analyse the contributors to this phylogenetic progress and make recommendations for future work. Conclusions: Our analyses suggest that a large majority of the vertebrate tree of life will: (1) be resolved within the next few decades; (2) identify specific data collection strategies that may help to spur future progress; and (3) identify branches of the vertebrate tree of life in need of increased research effort. [ABSTRACT FROM AUTHOR]

As phylogenetic data sets grow in size and number, objective methods to summarize this information are becoming increasingly important. Supermatrices can combine existing data directly and in principle provide effective syntheses of phylogenetic information that may reveal new relationships. However, several serious difficulties exist in the construction of large supermatrices that must be overcome before these approaches will enjoy broad utility. We present analyses that examine the performance of sparse supermatrices constructed from large sequence databases for the reconstruction of species-level phylogenies. We develop a largely automated informatics pipeline that allows for the construction of sparse supermatrices from GenBank data. In doing so, we develop strategies for alleviating some of the outstanding impediments to accurate phylogenetic inference using these approaches. These include taxonomic standardization, automated alignment, and the identification of rogue taxa. We use turtles as an exemplar clade and present a well-supported species-level phylogeny for two-thirds of all turtle species based on a ∼50 kb supermatrix consisting of 93% missing data. Finally, we discuss some of the remaining pitfalls and concerns associated with supermatrix analyses, provide comparisons to supertree approaches, and suggest areas for future research. [ABSTRACT FROM PUBLISHER]

Aposematism is one of the great mysteries of evolutionary biology. The evolution of aposematic coloration is poorly understood, but even less understood is the evolution of polymorphism in aposematic signals. Here, we use a phylogeographic approach to investigate the evolution of color polymorphism in Dendrobates pumilio, a well-known poison-dart frog (family Dendrobatidae), which displays perhaps the most striking color variation of any aposematic species. With over a dozen color morphs, ranging from bright red to dull green, D. pumilio provides an ideal opportunity to examine the evolution of color polymorphism and evolutionary shifts to cryptic coloration in an otherwise aposematic species. We constructed a phylogenetic tree for all D. pumilio color morphs from 3051bp of mtDNA sequence data, reconstructed ancestral states using parsimony and Bayesian methods, and tested the recovered tree against constraint trees using parametric bootstrapping to determine the number of changes to each color type. We find strong evidence for nearly maximal numbers of changes in all color traits, including five independent shifts to dull dorsal coloration. Our results indicate that shifts in coloration in aposematic species may occur more regularly than predicted and that convergence in coloration may indicate that similar forces are repeatedly driving these shifts. [ABSTRACT FROM AUTHOR]

Current conservation mitigation plans often fail to ensure full in-kind habitat replacement for endangered species, which suggests the need for improved methods for calculating mitigation credits. A simple, yet biologically meaningful method for calculating mitigation credits would be to let the number of mitigation credits assigned to a parcel of land scale with the reproductive value of the individuals occupying that parcel. This can be accomplished by dividing the population into 2 or more subdivisions with different reproductive values, calculating the densities of these subdivisions as a function of one or more habitat parameters, and then forming a weighted sum of these densities such that each density distribution is weighted by the reproductive value of its respective subdivision of the population. This weighted sum is the density distribution of reproductive value, and by integrating it over a particular parcel, one can determine the mitigation value of that parcel. We carried out this procedure for a population of California tiger salamanders (Ambystoma californiense ), with distance from breeding site as our habitat parameter and the 3 visually identifiable age classes (adults, juveniles, and metamorphs) as our population subdivisions. This led to a density distribution of reproductive value that decreased exponentially with increasing distance from a breeding site. Mitigation strategies derived from this function will be more likely to ensure the persistence of California tiger salamander populations than current approaches, which assign all land within 1.6 km of a breeding site the same mitigation value. Use of the density distribution of reproductive value as a basis for mitigation plans is a procedure that can be applied to all endangered species, and it should improve the quality of mitigation decisions. [ABSTRACT FROM AUTHOR]

Despite considerable effort from the systematics community, delimiting species boundaries in recent radiations remains a daunting challenge. We argue that genealogical approaches, although sometimes useful, may not solve this important problem, because recently derived species often have not had sufficient time to achieve monophyly. Instead, we suggest that population genetic approaches that rely on large sets of informative markers like single nucleotide polymorphisms (SNPs) provide an alternative framework for delimiting very recently derived species. We address two major challenges in applying such markers to species delimitation: discovering markers in nonmodel systems and using them to delimit recently derived species. Using turtles as a test case, we explore the utility of a single, relatively low-coverage genomic resource as an aid in gene and marker discovery. We exploit an end-sequenced bacterial artificial chromosome (BAC) library from an individual painted turtle (Chrysemys picta) and outline a novel protocol that efficiently identifies primer pairs that amplify homologous sequences across the tree of living turtles. Preliminary data using this library to discover SNPs in Emydura macquarii, a species that diverged from C. picta ∼ 210 million years ago, indicate that sequences identified from the Chrysemys BAC library provide useful SNPs even in this very distantly related taxon. Several recent methods in wide use in the population genetics literature allow one to discover potential species, or test existing species hypotheses, with SNP data and may be particularly informative for very recently derived species. As BAC and other genomic resources become increasingly available for scattered taxa across the tree of life, we are optimistic that these resources will provide abundant, inexpensive markers that will help delimit boundaries in problematic, recent species radiations. [ABSTRACT FROM AUTHOR]

Hybridization between differentiated lineages can have many different consequences depending on fitness variation among hybrid offspring. When introduced organisms hybridize with natives, the ensuing evolutionary dynamics may substantially complicate conservation decisions. Understanding the fitness consequences of hybridization is an important first step in predicting its evolutionary outcome and conservation impact. Here, we measured natural selection caused by differential viability of hybrid larvae in wild populations where native California Tiger Salamanders (Ambystoma californiense) and introduced Barred Tiger Salamanders (Ambystoma tigrinum mavortium) have been hybridizing for 50–60 years. We found strong evidence of hybrid vigor; mixed-ancestry genotypes had higher survival rates than genotypes containing mostly native or mostly introduced alleles. Hybrid vigor may be caused by heterozygote advantage (over-dominance) or recombinant hybrid vigor (due to epistasis or complementation). These genetic mechanisms are not mutually exclusive, and we find statistical support for both overdominant and recombinant contributions to hybrid vigor in larval tiger salamanders. Because recombinant homozygous genotypes can breed true, a single highly fit genotype with a mosaic of native and introduced alleles may eventually replace the historically pure California Tiger Salamander (listed as Threatened under the U.S. Endangered Species Act). The management implications of this outcome are complex: Genetically pure populations may not persist into the future, but average fitness and population viability of admixed California Tiger Salamanders may be enhanced. The ecological consequences for other native species are unknown. [ABSTRACT FROM AUTHOR]

Delimiting the boundaries of species involved in radiations is critical to understanding the tempo and mode of lineage formation. Single locus gene trees may or may not reflect the underlying pattern of population divergence and lineage formation, yet they constitute the vast majority of the empirical data in species radiations. In this study we make use of an expressed sequence tag (EST) database to perform nuclear (nDNA) and mitochondrial (mtDNA) genealogical tests of species boundaries in Ambystoma ordinarium, a member of an adaptive radiation of metamorphic and paedomorphic salamanders (the Ambystoma tigrinum complex) that have diversified across terrestrial and aquatic environments. Gene tree comparisons demonstrate extensive nonmonophyly in the mtDNA genealogy of A. ordinarium, while seven of eight independent nuclear loci resolve the species as monophyletic or nearly so, and diagnose it as a well-resolved genealogical species. A differential introgression hypothesis is supported by the observation that western A. ordinarium localities contain mtDNA haplotypes that are identical or minimally diverged from haplotypes sampled from a nearby paedomorphic species, Ambystoma dumerilii, while most nDNA trees place these species in distant phylogenetic positions. These results provide a strong example of how historical introgression can lead to radical differences between gene trees and species histories, even among currently allopatric species with divergent life history adaptations and morphologies. They also demonstrate how EST-based nuclear resources can be used to more fully resolve the phylogenetic history of species radiations. [ABSTRACT FROM AUTHOR]

We analysed phylogeography and population genetic variation across the range of the western pond turtle (Emys marmorata) using rapidly evolving mitochondrial and nuclear DNA sequence data. Nuclear DNA sequences from two unlinked introns displayed extremely low levels of variation, but phylogenetic analyses based on mtDNA recovered four well-supported and geographically coherent clades. These included a large Northern clade composed of populations from Washington south to San Luis Obispo County, California, west of the Coast Ranges; a San Joaquin Valley clade from the southern Great Central Valley; a geographically restricted Santa Barbara clade from a limited region in Santa Barbara and Ventura counties; and a Southern clade that occurs south of the Tehachapi Mountains and west of the Transverse Range south to Baja California, Mexico. An analysis of molecular variance (amova) based on regional hydrographic units revealed that populations from the Sacramento Valley north to Washington were virtually invariant, with no evidence of population substructure among northern river drainage basins. In other areas,E. marmoratacontains considerable unrecognized variation, particularly in central and southern California and in northern Baja California, Mexico. Our northern clade is congruent with the distribution of the subspeciesEmys marmorata marmorata(Washington–central California). However, no clade is congruent with the distribution of the southern subspeciesEmys marmorata pallidafrom central California–Baja. Thus, recognition of the current subspecies split is not warranted, based on the available genetic evidence. Ouramovaand phylogenetic results, in conjunction with a growing comparative database for other codistributed aquatic taxa, confirm the occurrence of genetic breaks across the Tehachapi Mountains and Transverse Range bounding the southern end of the Great Central Valley, and point to southern California as a rich source of cryptic genetic variation. [ABSTRACT FROM AUTHOR]

The red-legged frog, Rana aurora, has been recognized as both a single, polytypic species and as two distinct species since its original description 150 years ago. It is currently recognized as one species with two geographically contiguous subspecies, aurora and draytonii; the latter is protected under the US Endangered Species Act. We present the results of a survey of 50 populations of red-legged frogs from across their range plus four outgroup species for variation in a phylogenetically informative, ∼400 base pairs (bp) fragment of the mitochondrial cytochrome b gene. Our mtDNA analysis points to several major results. (1) In accord with several other lines of independent evidence, aurora and draytonii are each diagnosably distinct, evolutionary lineages; the mtDNA data indicate that they do not constitute a monophyletic group, but rather that aurora and R. cascadae from the Pacific northwest are sister taxa; (2) the range of the draytonii mtDNA clade extends about 100 km further north in coastal California than was previously suspected, and corresponds closely with the range limits or phylogeographical breaks of several codistributed taxa; (3) a narrow zone of overlap exists in southern Mendocino County between aurora and draytonii haplotypes, rather than a broad intergradation zone; and (4) the critically endangered population of draytonii in Riverside County, CA forms a distinct clade with frogs from Baja California, Mexico. The currently available evidence favours recognition of aurora and draytonii as separate species with a narrow zone of overlap in northern California. [ABSTRACT FROM AUTHOR]

After an estimated five million years of independent evolution, the barred tiger salamander (Ambystoma tigrinum mavortium) was introduced by bait dealers into the native range of the California tiger salamander (A. californiense). Hybridization and backcrossing have been occurring in central California for 50-60 years, or an estimated 15-30 generations. We studied genetic and ecological factors influencing admixture of these two divergent gene pools by analyzing frequencies of hybrid genotypes in three kinds of breeding habitats: natural vernal pools, ephemeral man-made cattle ponds, and perennial man-made ponds. Perennial ponds tended to have higher frequencies of nonnative alleles than either type of seasonal pond, even in cases where perennial and seasonal ponds are within a few hundred meters. Thus, the hybrid zone has a mosaic structure that depends on pond hydrology or ecology. The presence of some broadly acting constraints on admixture is suggested by linkage disequilibria between physically unlinked molecular markers within ponds. In addition, we found several marker-specific deviations from Hardy-Weinberg equilibrium. One marker showed a consistent deficit of heterozygotes across pond types. Another showed heterozygote deficits only in vernal pools. A third was more likely to have heterozygote excess in ephemeral cattle ponds. These patterns indicate that admixture is influenced by complex genotype-by-environment interactions. [ABSTRACT FROM AUTHOR]

Provides information on a study that analyzed variation in the rapidly evolving mitochondrial control region in 241 turtles from 117 localities across North America from southern Canada to extreme northern Mexico to examine whether the painted turtle, Chrysemys picta represents a continentally distributed species based on molecular analysis. Methodology of the study; Results and discussion on the study.

The number, size, and scope of phylogenetic analyses using molecular data has increased steadily in recent years. This has simultaneously led to a dramatic improvement in our understanding of phylogenetic relationships and a better appreciation for an array of methodological problems that continue to hinder progress in phylogenetic studies of certain data sets and/or particular parts of the tree of life. This review focuses on several persistent problems, including rooting, conflict among data sets, weak support in trees, strong but evidently incorrect support, and the computational issues arising when methods are applied to the large data sets that are becoming increasingly commonplace. We frame each of these issues as a specific problem to be overcome, review the relevant theoretical and empirical literature, and suggest solutions, or at least strategies, for further investigation of the issues involved. [ABSTRACT FROM AUTHOR]

Wind-borne pesticides have long been suggested as a cause of amphibian declines in areas without obvious habitat destruction. In California, the transport and deposition of pesticides from the agriculturally intensive Central Valley to the adjacent Sierra Nevada is well documented, and pesticides have been found in the bodies of Sierra frogs. Pesticides are therefore a plausible cause of declines, but to date no direct links have been found between pesticides and actual amphibian population declines. Using a geographic information system, we constructed maps of the spatial pattern of declines for eight declining California amphibian taxa, and compared the observed patterns of decline to those predicted by hypotheses of wind-borne pesticides, habitat destruction, ultraviolet radiation, and climate change. In four species, we found a strong positive association between declines and the amount of upwind agricultural land use, suggesting that wind-borne pesticides may be an important factor in declines. For two other species, declines were strongly associated with local urban and agricultural land use, consistent with the habitat-destruction hypothesis. The patterns of decline were not consistent with either the ultraviolet radiation or climate-change hypotheses for any of the species we examined. [ABSTRACT FROM AUTHOR]

The intense exploitation of turtles in Asian markets has contributed to declines in turtle populations across the continent. Three-quarters of Asia’s turtles are threatened and half are endangered. A recent workshop on the Asian turtle crisis identified taxonomic studies of widespread species as a priority for research because these low risk species may include unrecognized, narrowly distributed taxa of much higher concern. Chitra indica is a widely exploited softshell turtle (family Trionychidae) found across southern Asia. Individuals from Thailand have been described as a separate species, Chitra chitra, but this has not been universally accepted, and many sources consider Chitra monotypic. Phylogenetic analysis of sequence data from the mitochondrial ND4 gene revealed three deeply divergent, monophyletic lineages within Chitra: C. indica, C. chitra, and a third unnamed form from Myanmar. This new form is probably Critically Endangered, which highlights the importance of systematic studies in determining conservation priorities. [Copyright &y& Elsevier]

We used two different experimental approaches to test the hypothesis that thyroid hormone receptor (TR) variation is associated with alternate life cycles modes in ambystomatid salamanders. In the first experiment, the inheritance of TRα and TRβ genotypes was determined for metamorphic and non metamorphic offspring from backcrosses between Ambystoma mexicanum (an obligate metamorphic-failure species) and metamorphic F1 hybrids (A. mexicanum × A. tigrinum tigrinum). The segregation of TR genotype was independent of the expression of life cycle mode phenotype, and neither TR locus was linked to DNA markers that flank a major-effect locus for life cycle mode. In the second experiment, a portion of the ligand-binding domain of TRα and TRβ was cloned and sequenced for DNA samples from 14 different ambystomatid salamander populations, including obligate metamorphic, facultative metamorphic, and obligate metamorphic-failure taxa. Nucleotide sequence variation was found for both TRα and TRβ, with several nonsynonomous substitutions that presumably code for nonconservative amino acid replacements. However, no general relationship was found between TR allelic variation and life cycle mode among populations or species. These data do not implicate TRs as candidate loci involved in the current maintenance or past evolution of alternate life cycle modes in members of the tiger salamander complex. [ABSTRACT FROM AUTHOR]

Abstract We present a comprehensive survey of genetic variation across the range of the narrowly distributed endemic Yosemite toad Bufo canorus, a declining amphibian restricted to the Sierra Nevada of California. Based on 322 bp of mitochondrial cytochrome b sequence data, we found limited support for the monophyly of B. canorus and its closely related congener B. exsul to the exclusion of the widespread western toad B. boreas. However, B. exsul was always phylogenetically nested within B. canorus, suggesting that the latter may not be monophyletic. SSCP (single-strand conformation polymorphism) analysis of 372 individual B. canorus from 28 localities in Yosemite and Kings Canyon National Parks revealed no shared haplotypes among these two regions and lead us to interpret these two parks as distinct management units for B. canorus. Within Yosemite, we found significant genetic substructure both at the level of major drainages and among breeding ponds. Kings Canyon samples show a different pattern, with substantial variation among breeding sites, but no substructure among drainages. Across the range of B. canorus as well as among Yosemite ponds, we found an isolation-by-distance pattern suggestive of a stepping stone model of migration. However, in Kings Canyon we found no hint of such a pattern, suggesting that movement patterns of toads may be quite different in these nearby parklands. Our data imply that management for B. canorus should focus at the individual pond level, and effective management may necessitate reintroductions if local extirpations occur. A brief review of other pond-breeding anurans suggests that highly structured populations are often the case, and thus that our results for B. canorus may be general for other species of frogs and toads. [ABSTRACT FROM AUTHOR]

Presents a study which used paedomorphosis in the Mexican axolotl to analyse adaptive evolution. Description of paedomorphosis; Information on the gene hypothesis for the evolution of paedomorphosis in the Mexican axolotl; Methods used to conduct the study; Discussion on results of study.

Presents a study on the genetic differentiation on the tiger salamander found in the United States and Mexico. Materials and methods; Results; Discussion.

The relationship between genomics and the core disciplines of ecology, evolution, and systematics is rapidly developing, computationally intensive, and changing the ways that our science asks and answers its most fundamental questions. As the six papers under this theme in this volume of AREES demonstrate, comparative and population genomics stand to fundamentally change how we study the nature of adaptation, the tempo and mode of evolution, and the history of life on Earth. Genomics is enabling our understanding of questions ranging from human historical demography to conservation biology of poorly known plants and animals, enhancing the role of organismal biology across basic and applied science. [ABSTRACT FROM AUTHOR]

As genomic-scale data sets become economically feasible for most organisms, a key question for conservation biology is whether the increased resolution offered by new genomic approaches justifies repeating earlier studies based on traditional markers, rather than investing those same time and monetary resources in less-known species. Genomic studies offer clear advantages when the objective is to identify adaptive loci that may be critical to conservation policy-makers. However, the answer is far less certain for the population and landscape studies based on neutral loci that dominate the conservation genetics research agenda. We used Restriction-site Associated DNA sequencing (RADseq) to revisit earlier molecular studies of the IUCN Critically Endangered Magdalena River turtle (Podocnemis lewyana), documenting the conservation insights gained by increasing the number of neutral markers by several orders of magnitude. Earlier research indicated that P. lewyana has the lowest genetic diversity known for any chelonian, and little or no population differentiation among independent rivers. In contrast, the RADseq data revealed discrete population structure with isolation-by-distance within river segments and identified precise population breaks clearly delineating management units. It also confirmed that the species does not have extremely low heterozygosity and that effective population sizes are probably sufficient to maintain long-term evolutionary potential. Contrary to earlier inferences from more limited population genetic markers, our genomic data suggest that management strategies should shift from active genetic rescue to more passive protection without extreme interventions. We conclude with a list of examples of conservation studies in other vertebrates indicating that for many systems a genomic update is worth the investment. [ABSTRACT FROM AUTHOR]

The phylogeny of California Slender Salamanders (Batrachoseps attenuatus), which inhabits Northern California's Coast Ranges and Sierra Nevada foothills, was previously investigated by Martínez-Solano et al. (2007), who recovered a monophyletic species that contained five geographically cohesive major clades. However, a population isolate of B. attenuatus found in the Sutter Buttes, a geographically isolated volcanic mountain range within the Sacramento Valley, was not available in their phylogenetic analysis. We sequenced the ∼780-base pair segment of the mitochondrial DNA cytochrome b (cytb) gene used in Martínez-Solano et al. (2007) to examine the phylogenetic placement of this isolated population within the B. attenuatus complex. Our results, based on a single mitochondrial marker, indicate that the Sutter Buttes population of B. attenuatus is genetically distinct and falls within the Martínez-Solano et al. (2007) "East Bay/Sierra subclade" that, in turn, is nested within their "Eastern" major clade. Although we cannot accurately estimate how long the Sutter Buttes population has been isolated, B. attenuatus was likely able to colonize the Sutter Buttes during pluvial Pleistocene interglacials when more contiguous woodlands may have opened dispersal corridors across the Great Central Valley. [ABSTRACT FROM AUTHOR]

Living turtles are characterized by extraordinarily low species diversity given their age. The clade's extensive fossil record indicates that climate and biogeography may have played important roles in determining their diversity. We investigated this hypothesis by collecting a molecular dataset for 591 individual turtles that, together, represent 80% of all turtle species, including representatives of all families and 98% of genera, and used it to jointly estimate phylogeny and divergence times. We found that the turtle tree is characterized by relatively constant diversification (speciation minus extinction) punctuated by a single threefold increase. We also found that this shift is temporally and geographically associated with newly emerged continental margins that appeared during the Eocene-Oligocene transition about 30 million years before present. In apparent contrast, the fossil record from this time period contains evidence for a major, but regional, extinction event. These seemingly discordant findings appear to be driven by a common global process: global cooling and drying at the time of the Eocene-Oligocene transition. This climatic shift led to aridification that drove extinctions in important fossilbearing areas, while simultaneously exposing new continental margin habitat that subsequently allowed for a burst of speciation associated with these newly exploitable ecological opportunities. [ABSTRACT FROM AUTHOR]

Carpenter ants in the genus Camponotus are large, conspicuous ants that are abundant and ecologically influential in many terrestrial ecosystems. The bicolored carpenter ant, Camponotus vicinus Mayr, is distributed across a wide range of elevations and latitudes in western North America, where it is a prominent scavenger and predator. Here, we present a high-quality genome assembly of C. vicinus from a sample collected in Sonoma County, California, near the type locality of the species. This genome assembly consists of 38 scaffolds spanning 302.74 Mb, with contig N50 of 15.9 Mb, scaffold N50 of 19.9 Mb, and BUSCO completeness of 99.2%. This genome sequence will be a valuable resource for exploring the evolutionary ecology of C. vicinus and carpenter ants generally. It also provides an important tool for clarifying cryptic diversity within the C. vicinus species complex, a genetically diverse set of populations, some of which are quite localized and of conservation interest. [ABSTRACT FROM AUTHOR]

The Yuma myotis bat (Myotis yumanensis) is a small vespertilionid bat and one of 52 species of new world Myotis bats in the subgenus Pizonyx. While M. yumanensis populations currently appear relatively stable, it is one of 12 bat species known or suspected to be susceptible to white-nose syndrome, the fungal disease causing declines in bat populations across North America. Only two of these 12 species have genome resources available, which limits the ability of resource managers to use genomic techniques to track the responses of bat populations to white-nose syndrome generally. Here we present the first de novo genome assembly for Yuma myotis, generated as a part of the California Conservation Genomics Project. The M. yumanensis genome was generated using a combination of PacBio HiFi long reads and Omni-C chromatin-proximity sequencing technology. This high-quality genome is one of the most complete bat assemblies available, with a contig N50 of 28.03 Mb, scaffold N50 of 99.14 Mb, and BUSCO completeness score of 93.7%. The Yuma myotis genome provides a high-quality resource that will aid in comparative genomic and evolutionary studies, as well as inform conservation management related to white-nose syndrome. [ABSTRACT FROM AUTHOR]

Rattlesnakes play important roles in their ecosystems by regulating prey populations, are involved in complex coevolutionary dynamics with their prey, and exhibit a variety of unusual adaptations, including maternal care, heat-sensing pit organs, hinged fangs, and medically-significant venoms. The western rattlesnake (Crotalus oreganus) is one of the widest ranging rattlesnake species, with a distribution from British Columbia, where it is listed as threatened, to Baja California and east across the Great Basin to western Wyoming, Colorado and New Mexico. Here, we report a new reference genome assembly for one of six currently recognized subspecies, C. oreganus helleri , as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomic sequencing 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. The assembly comprises a total of 698 scaffolds spanning 1,564,812,557 base pairs, has a contig N50 of 64.7 Mb, a scaffold N50 of 110.8 Mb, and BUSCO complete score of 90.5%. This reference genome will be valuable for studies on the genomic basis of venom evolution and variation within Crotalus , in resolving the taxonomy of C. oreganus and its relatives, and for the conservation and management of rattlesnakes in general. [ABSTRACT FROM AUTHOR]

Snakes in the family Colubridae include more than 2,000 currently recognized species, and comprise roughly 75% of the global snake species diversity on Earth. For such a spectacular radiation, colubrid snakes remain poorly understood ecologically and genetically. Two subfamilies, Colubrinae (788 species) and Dipsadinae (833 species), comprise the bulk of colubrid species richness. Dipsadines are a speciose and diverse group of snakes that largely inhabit Central and South America, with a handful of small-body-size genera that have invaded North America. Among them, the ring-necked snake, Diadophis punctatus , has an incredibly broad distribution with 14 subspecies. Given its continental distribution and high degree of variation in coloration, diet, feeding ecology, and behavior, the ring-necked snake is an excellent species for the study of genetic diversity and trait evolution. Within California, six subspecies form a continuously distributed "ring species" around the Central Valley, while a seventh, the regal ring-necked snake, Diadophis punctatus regalis is a disjunct outlier and Species of Special Concern in the state. Here, we report a new reference genome assembly for the San Diego ring-necked snake, D. p. similis , as part of the California Conservation Genomics Project. This assembly comprises a total of 444 scaffolds spanning 1,783 Mb and has a contig N50 of 8.0 Mb, scaffold N50 of 83 Mb, and BUSCO completeness score of 94.5%. This reference genome will be a valuable resource for studies of the taxonomy, conservation, and evolution of the ring-necked snake across its broad, continental distribution. [ABSTRACT FROM AUTHOR]

The California Pipevine Swallowtail Butterfly, Battus philenor hirsuta , and its host plant, the California Pipevine or Dutchman's Pipe, Aristolochia californica Torr., are an important California endemic species pair. While this species pair is an ideal system to study co-evolution, genomic resources for both are lacking. Here, we report a new, chromosome-level assembly of B. philenor hirsuta 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 109 scaffolds spanning 443 mega base (Mb) pairs, with a contig N50 of 14.6 Mb, a scaffold N50 of 15.2 Mb, and BUSCO complete score of 98.9%. In combination with the forthcoming A. californica reference genome, the B. philenor hirsuta genome will be a powerful tool for documenting landscape genomic diversity and plant–insect co-evolution in a rapidly changing California landscape. [ABSTRACT FROM AUTHOR]

Genomic data are increasingly used for high resolution population genetic studies including those at the forefront of biological conservation. A key methodological challenge is determining sequence similarity clustering thresholds for RADseq data when no reference genome is available. These thresholds define the maximum permitted divergence among allelic variants and the minimum divergence among putative paralogues and are central to downstream population genomic analyses. Here we develop a novel set of metrics to determine sequence similarity thresholds that maximize the correct separation of paralogous regions and minimize oversplitting naturally occurring allelic variation within loci. These metrics empirically identify the threshold value at which true alleles at opposite ends of several major axes of genetic variation begin to incorrectly separate into distinct clusters, allowing researchers to choose thresholds just below this value. We test our approach on a recently published data set for the protected foothill yellow‐legged frog (Rana boylii). The metrics recover a consistent pattern of roughly 96% similarity as a threshold above which genetic divergence and data missingness become increasingly correlated. We provide scripts for assessing different clustering thresholds and discuss how this approach can be applied across a wide range of empirical data sets. [ABSTRACT FROM AUTHOR]

We developed a species-specific, quantitative polymerase chain reaction assay multiplex for the detection of 2 threatened and endangered anuran species, Arroyo Toad (Anaxyrus californicus Camp, 1915) and Western Spadefoot (Spea hammondii Baird, 1859), from environmental DNA (eDNA). Both species have experienced range wide declines over the last century, mostly because of habitat loss. More accurate methods to detect and monitor both species are needed given their declines and cryptic life histories. We tested species-specificity and sensitivity of eDNA assays in the laboratory and validated the multiplexed assays with field-collected eDNA samples. Both assays were species-specific and sensitive and effectively detected the target species from streams, rivers, and ponds, although detection was imperfect. These assays provide an efficient and economical tool to detect, monitor, manage, and conserve these at-risk species. [ABSTRACT FROM AUTHOR]

Spiny lizards (genus Sceloporus) have long served as important systems for studies of behavior, thermal physiology, dietary ecology, vector biology, speciation, and biogeography. The western fence lizard, Sceloporus occidentalis , is found across most of the major biogeographical regions in the western United States and northern Baja California, Mexico, inhabiting a wide range of habitats, from grassland to chaparral to open woodlands. As small ectotherms, Sceloporus lizards are particularly vulnerable to climate change, and S. occidentalis has also become an important system for studying the impacts of land use change and urbanization on small vertebrates. Here, we report a new reference genome assembly for S. occidentalis , as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomics 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. The assembly comprises a total of 608 scaffolds spanning 2,856 Mb, has a contig N50 of 18.9 Mb, a scaffold N50 of 98.4 Mb, and BUSCO completeness score of 98.1% based on the tetrapod gene set. This reference genome will be valuable for understanding ecological and evolutionary dynamics in S. occidentalis , the species status of the California endemic island fence lizard (S. becki), and the spectacular radiation of Sceloporus lizards. [ABSTRACT FROM AUTHOR]

A critical consideration when using molecular ecological methods to detect trends and parameterize models at very fine spatial and temporal scales has always been the technical limits of resolution. Key landscape features, including most anthropogenic modifications, can cause biologically important, but very recent changes in gene flow that require substantial statistical power to detect. The problem is one of temporal scale: Human change is rapid and recent, while genetic changes accumulate slowly. We generated SNPs from thousands of nuclear loci to characterize the population structure of New York‐endangered eastern tiger salamanders (Ambystoma tigrinum) on Long Island and quantify the impacts of roads on population fragmentation. In stark contrast to a recent microsatellite study, we uncovered highly structured populations over an extremely small spatial scale (approximately 40 km2) in an increasingly human‐modified landscape. Geographic distance and the presence of roads between ponds were both strong predictors of genetic divergence, suggesting that both natural and anthropogenic factors contribute to the observed patterns of genetic variation. All ponds supported small to modest effective breeding populations, and pond surface area showed a strong positive correlation with population size. None of these patterns emerged in an earlier study of the same system using microsatellite loci, and we determined that at least 300–400 SNPs were needed to recover the fine‐scale population structure present in this system. Conservation assessments using earlier genetic techniques in other species may similarly lack the statistical power for small‐scale inferences and benefit from reassessments using genomic tools. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Genome assemblies are increasingly being used to identify adaptive genetic variation that can help prioritize the population management of protected species. This approach may be particularly relevant to species like Blainville's horned lizard, Phrynosoma blainvillii , due to its specialized diet on noxious harvester ants, numerous adaptative traits for avoiding predation (e.g. cranial horns, dorsoventrally compressed body, cryptic coloration, and blood squirting from the orbital sinuses), and status as Species of Special Concern in California. Rangewide decline since the early 20th century, the basis of its conservation status, has been driven mainly by habitat conversion, over-collecting, and invasion of a non-native ant that displaces its native ant prey base. Here, we report on a scaffold-level genome assembly for P. blainvillii as part of the California Conservation Genomics Project (CCGP), produced using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology. The de novo assembly has 78 scaffolds, a total length of ~2.21 Gb, a scaffold N50 length of ~352 Mb, and BUSCO score of 97.4%. This is the second species of Phrynosoma for which a reference genome has been assembled and represents a considerable improvement in terms of contiguity and completeness. Combined with the landscape genomics data being compiled by the CCGP, this assembly will help strategize efforts to maintain and/or restore local genetic diversity, where interventions like genetic rescue, translocation, and strategic land preservation may be the only means by which P. blainvillii and other low-vagility species can survive in the fragmented habitats of California. [ABSTRACT FROM AUTHOR]

The Virginia rail, Rallus limicola , is a member of the family Rallidae, which also includes many other species of secretive and poorly studied wetland birds. It is recognized as a single species throughout its broad distribution in North America where it is exploited as a game bird, often with generous harvest limits, despite a lack of systematic population surveys and evidence of declines in many areas due to wetland loss and degradation. To help advance understanding of the phylogeography, biology, and ecology of this elusive species, we report the first reference genome assembly for the Virginia rail, produced as part of the California Conservation Genomics Project (CCGP). We produced a de novo genome assembly using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology with an estimated sequencing error rate of 0.191%. The assembly consists of 1,102 scaffolds spanning 1.39 Gb, with a contig N50 of 11.0 Mb, scaffold N50 of 25.3 Mb, largest contig of 45 Mb, and largest scaffold of 128.4 Mb. It has a high BUSCO completeness score of 96.9% and represents the first genome assembly available for the genus Rallus. This genome assembly will help resolve questions about the complex evolutionary history of rails and evaluate the potential of rails for adaptive evolution in the face of growing threats from climate change and habitat loss and fragmentation. It will also provide a valuable resource for rail conservation efforts by quantifying Virginia rail vagility, population connectivity, and effective population sizes. [ABSTRACT FROM AUTHOR]

Climate-driven changes in hydrological regimes are of global importance and are particularly significant in riparian ecosystems. Riparian ecosystems in California provide refuge to many native and vulnerable species within a xeric landscape. California Tetragnatha spiders play a key role in riparian ecosystems, serving as a link between terrestrial and aquatic elements. Their tight reliance on water paired with the widespread distributions of many species make them ideal candidates to better understand the relative role of waterways versus geographic distance in shaping the population structure of riparian species. To assist in better understanding population structure, we constructed a reference genome assembly for Tetragnatha versicolor using long-read sequencing, scaffolded with proximity ligation Omni-C data. The near-chromosome-level assembly is comprised of 174 scaffolds spanning 1.06 Gb pairs, with a scaffold N50 of 64.1 Mb pairs and BUSCO completeness of 97.6%. This reference genome will facilitate future study of T. versicolor population structure associated with the rapidly changing environment of California. [ABSTRACT FROM AUTHOR]

The black rail, Laterallus jamaicensis , is one of the most secretive and poorly understood birds in the Americas. Two of its five subspecies breed in North America: the Eastern black rail (L. j. jamaicensis), found primarily in the southern and mid-Atlantic states, and the California black rail (L. j. coturniculus), inhabiting California and Arizona, are recognized across the highly disjunct distribution. Population declines, due primarily to wetland loss and degradation, have resulted in conservation status listings for both subspecies. To help advance understanding of the phylogeography, biology, and ecology of this elusive species, we report the first reference genome assembly for the black rail, produced as part of the California Conservation Genomics Project (CCGP). We produced a de novo genome assembly using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology with an estimated sequencing error rate of 0.182%. The assembly consists of 964 scaffolds spanning 1.39 Gb, with a contig N50 of 7.4 Mb, scaffold N50 of 21.4 Mb, largest contig of 44.8 Mb, and largest scaffold of 101.2 Mb. The assembly has a high BUSCO completeness score of 96.8% and represents the first genome assembly available for the genus Laterallus. This genome assembly can help resolve questions about the complex evolutionary history of rails, assess black rail vagility and population connectivity, estimate effective population sizes, and evaluate the potential of rails for adaptive evolution in the face of growing threats from climate change, habitat loss and fragmentation, and disease. [ABSTRACT FROM AUTHOR]