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

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


Background & Summary
The CaliPopGen database consists of four datasets that contain estimates of population genetic diversity, differentiation, and life history traits for 448 eukaryotic species sampled across California, USA. The state has exceptionally high plant and animal biodiversity, and a correspondingly large number of endangered taxa 1 . It is often divided into 19 terrestrial and three marine ecoregions, reflecting its tremendous geologic and ecological diversity 2,3 , including the highest and lowest elevations in the contiguous U.S., extreme deserts and temperate rainforests, and mean annual precipitation ranging from 150 mm-1200 mm 4 . California is the most populous state in the USA, accommodating roughly 12% of the nation's human population; the third largest state geographically, encompassing 5% of the country's continental land area; and is a major agricultural producer. This combination of high species richness and human-mediated pressures constitute a persistent threat to the short-and long-term persistence of biodiversity, and has led to California's inclusion as one of only two global biodiversity hotspots in the USA 5,6 . Perhaps unsurprisingly, California has the greatest number of documented and possibly extinct species of vascular plants 7 , and more than twice as many federally protected species (total of 287) as any other state 8 . It has also been the focus of more population genetic studies, including states with similarly high numbers of threatened/endangered species like Florida and Hawaii 1 . However, this wealth of genetic information has never been adequately summarized or made publicly available. The few broadly comparative analyses for California have largely been based on inferences derived from fewer than 10 species [9][10][11][12] , with the exceptions of one review 13 , and one empirical study 14 , both of which were restricted to marine taxa. California is a perennial leader in biodiversity management, and our compilation of genetic data for the state aligns with the administrative level at which environmental legislation and biodiversity management is implemented, increasing the likelihood that the CaliPopGen database will inform conservation actions.
In compiling CaliPopGen, we examined 4,942 published studies identified by our search criteria in the Web of Science, of which 450 met our final inclusion criteria and are included in the database. The majority of genetic samples represented in this database were collected from 1995-2015 (ranging from 1888-2019), and all studies were published between 1983-2020. CaliPopGen contains information on more than 187,394 individuals from 5,453 unique populations, of which 5,276 are spatially georeferenced. These populations include terrestrial (46.6%), marine (21.9%), freshwater (14.1%), amphibious (9.7%), and diadromous (7.7%) populations of fungi (<2% of unique species), chromists (<2%), plants (23%), and animals (73%; Fig. 1). CaliPopGen includes population level data with broader taxonomic coverage than recent, more global compilations, which have focused on freshwater and marine fishes 15 , mammals 16 , mammals and amphibians 17 , vertebrates 18 , and birds, fishes, insects and mammals 19 . Its focus at the regional (state) level is unique. The CaliPopGen database also includes a wider range of molecular markers (Fig. 2), populations and species than these previous population genetic compilations. Molecular markers in our database include RFLPs, AFLPs, allozymes and isozymes, microsatellites, mitochondrial, and other nuclear markers, whereas previously published datasets frequently focussed on one or a few loci (e.g. [15][16][17]19 ) or single marker types (e.g. 18 ,). Our inclusion of all available marker types both reflects the change in methodological approaches through time (for example, the temporal replacement of allozymes with microsatellites in the early 2000's), and presents opportunities for quantitative comparisons among different marker types.
Expressed as a proportion of the study area, CaliPopGen contains at least an order of magnitude more species (0.83/1,000 km²), populations (9.59/1,000 km²), and individuals (284.04/1,000 km 2 ) than the aforementioned studies and databases. This high spatial density of samples across the full ecological scope of California (Fig. 3) should facilitate future analyses of ecological trends at the population level where biological processes actually occur, and is well suited to help identify important mechanisms shaping genetic diversity, connectivity and fragmentation. CaliPopGen should also serve as a point of departure for future studies, providing a genetic baseline against which researchers can contrast and quantify future population genetic impacts resulting from changes in climate or land use. As such, CaliPopGen is an historical antecedent to ongoing genomic initiatives to study www.nature.com/scientificdata www.nature.com/scientificdata/ the diversity and distribution of California's flora and fauna, including the California Conservation Genomics Project 20 , and other projects using landscape genomic approaches.
To supplement the genetic data in CaliPopGen, we also compiled datasets containing life history information for all plant and animal species in the database, including adult body size, lifespan, reproductive and dispersal traits, and conservation status.

Methods
Population genetic data collection from primary data sources. Figure 4 describes the overall data collection workflow for the four datasets that comprise CaliPopGen. We first identified literature potentially containing population genetic data for California by querying the Web of Science Core Collection (https:// webofknowledge.com/) for relevant literature from 1900 to 2020 with the terms: topic = (California*) AND topic = (genetic* OR genomic*) AND topic = (species OR taxa* OR population*). We included only empirical peer-reviewed literature and excluded unreviewed preprints. In using these search terms, our goal was to broadly identify genetic papers focused on California with population or species-level analyses, while avoiding purely phylogenetic studies or those focused on agricultural or model species. This resulted in 4,942 unique records.
We next screened titles and abstracts to retain articles that: (1) provided data on populations of species which are self-sustaining without anthropogenic involvement; (2) included at least some eukaryote species; (3) www.nature.com/scientificdata www.nature.com/scientificdata/ included population(s) sampled within California; (4) mentioned measures of genetic diversity or differentiation; and (5) were not reviews (thus restricting our search to only primary literature). We retained 1869 studies after this first pass of literature screening (see Technical Validation for estimate of inter-and intra-screener bias).
Our second, more in-depth screening pass involved reading the full text of these 1869 studies. We had two goals. First, we confirmed that retained papers fully met all five of our inclusion criteria (the first screen was very liberal with respect to these criteria, and many papers failed to meet at least one criterion after close reading). Second, we eliminated papers where the data were not presented in a way that allowed us to extract population-level information. For example, many of the more systematics-focused studies pooled samples from large, somewhat ill-defined regions ("Sierra Nevada" or "Southern California"); if such regions were larger than 50 km in a linear dimension, we deemed them unusable for making geographically-informative inferences. Other studies presented summaries of population data, often in the form of phylogenetic networks or trees, but did not include information on actual population genetic parameters and therefore were not relevant to our database. We retained 528 publications after this second pass.
From this set of papers, we extracted species, locality, and genetic data for each California population or sampling locality described in each study (Fig. 3A). This included Latin binomial/trinomial, English common name, population identifiers, and geographic coordinates of sampling sites. We also noted population/sampling www.nature.com/scientificdata www.nature.com/scientificdata/ localities that were interpreted as comprised of interspecific hybrids, and listed both parental species. We collected population genetic diversity and differentiation statistics for each unique genetic marker for each population/sampling locality; as a result, a sampling locality may have multiple entry rows, one for each locus or marker type. Parameters extracted for each population/marker combination include sample size, genetic marker type, gene targets, number of loci, years of sampling, and reported values for effective population size (N e ), expected (H E ) and observed (H O ,) heterozygosity, nucleotide diversity (π, pi), alleles-per-locus (APL), allelic richness (A R ), percent polymorphic loci (PPL), haplotype diversity (HDIV), inbreeding coefficient (e.g. F IS , F IT , G IS ), and pairwise population genetic comparison parameters (F ST , G ST , D ST , Nei's D, Jost's D, or phi). We note that while there are technical differences between allelic richness and alleles-per-locus, source literature often used the terms interchangeably, and we include the parameters and their values as named in the source. We define marker type as the general category of genetic marker used (e.g., "microsatellite" or "nuclear"), while gene targets are the specific locus/loci (e.g., "COI"). We present these data in two separate datasets, one containing all population-level genetic summary statistics (Dataset 1 21 , see Fig. 3C and detailed description in Table 1) and a second for estimates of pairwise genetic differentiation (Dataset 2 21 , see Fig. 3D and detailed description in Table 2).
All genetic data were extracted directly from the source literature. However, we also updated or added to the metadata for these population genetic values in several ways. We included kingdom, phylum, and a lower-level taxonomic grouping for each species (usually class), and updated scientific and common names based on the currently accepted taxonomy of the Global Biodiversity Information Facility 22 . When geographic coordinates were not provided for a sampling locality, as was frequently the case in the older literature, we used Google Maps (https://www.google.com/maps) to georeference localities based on either in-text descriptions or embedded figure maps guided by permanent landmarks like a bend in a river or administrative boundaries. Because this can only yield approximate coordinates, we recorded estimated accuracy as the radius of our best estimate of possible error in kilometers. If coordinates were provided in degree/minute/seconds, we used Google Maps to translate them to decimal degrees. In cases where coordinates were not provided and locality descriptions were too vague to determine coordinates with less than 50 km estimated coordinate error, we did not attempt to extract coordinates but still provide the genetic data. All coordinates are provided in the web Mercator projection (EPSG:3857). We excluded studies that reported genetic parameter values only for samples aggregated regionally ("Southern California" or "Sierra Nevada"). If marker type was not explicitly included, we classified marker type based on the gene targets reported, if provided.

Life history trait data collection.
To increase the utility of CaliPopGen, we also assembled data on life history traits for all animal (Dataset 3 21 ) and plant (Dataset 4 21 ) species contained in Datasets 1 21 and 2 21 . We assembled trait data that have previously been shown to correlate with genetic diversity, including those related to reproduction, life cycle, and body size, as well as conservation status (e.g. 23-26 ,). Life history data were compiled by first referencing large online repositories, often specific to taxonomic groups, like the TRY plant trait database 27 , and the Royal Botanic Gardens Kew Seed Information Database 28 . If trait data for species of interest were unavailable from these compilations, we conducted keyword literature searches for each combination of species and life history trait, and extracted data from the primary literature. When data were not available for the subspecies or species for which we had genetic data, we report values for the next closest taxonomic level, up to and including family, as available in the literature.
For both animals and plants, we defined habitat types as marine, freshwater, diadromous, amphibious, or terrestrial. Marine species include those that are found in brackish or wetland-marine habitats, as well as bird species that primarily reside in marine habitats. Freshwater species include those that are found in wetland-freshwater habitats, as well as species that primarily reside in freshwater. The diadromous category www.nature.com/scientificdata www.nature.com/scientificdata/ includes fish species that are catadromous or anadromous. We considered species to be amphibious if they have an obligatory aquatic stage in their life cycle, but also spend a significant portion of their life cycle on land. Terrestrial species were defined as those that spend most of their life cycle on land and are not aquatic for any portion of their life cycle. In a few cases (e.g., waterbirds that are both freshwater and marine, semi-aquatic reptiles), a species could reasonably be placed in more than one category, and we did our best to identify the primary life history category for such taxa. If the taxonomic identity of an entry was hybrid between species or subspecies, this was noted in the speciesID column and no life history data were reported.
The CaliPopGen Animal Life History Traits Dataset 3 21 (description of dataset in Table 3) includes habitat type, lifespan, fecundity, lifetime reproductive success, age at sexual maturity, number of breeding events per year, mode of reproduction, adult length and mass, California native status, listing status under the US Endangered Species Act (ESA), listing status under the California Endangered Species Act (CESA), and status as a California Species of Special Concern (SSC). For some traits, value ranges were recorded-for example, minimum to maximum lifespan. In other cases, we recorded single values and, when available, a definition of this single value, (for example, minimum, average, or maximum lifespan). We report either the range of the age of sexual maturity (minimum to maximum), or a single value, depending on the available literature. For sexually dimorphic species, we report female adult length and weight when available, because female body size often correlates with fecundity. Across animal taxonomic groups, different measures of body size and length measurements are often used, reflecting community consensus on how to measure size. Given this variation, we report the type of length measurement, if available, as Standard Length (SL), Fork Length (FL), Total Length (TL), Snout-to-Vent Length (SVL), Straight-Line Carapace (SLC), or Wingspan (WS).
The CaliPopGen Plant Life History Traits Dataset 4 21 (description of dataset in Table 4) includes habitat type, lifespan, life cycle, adult height, self-compatibility, monoecious or dioecious, mode of reproduction, pollination and seed dispersal modes, mass per seed, California native status, NatureServe 29 element ranks (global and state ranks, see Table 5 Table 3) Columns 71-101 Plant life history data (see Table 4) Table 1. Description of the population genetic data in Dataset 1 21 . www.nature.com/scientificdata www.nature.com/scientificdata/ status under the California Endangered Species Act (CESA). In contrast to most animal species, plant lifespan was typically reported as a single value. We define life cycles as the following: Annual: completes full life cycle in one year; Biennial: completes full life cycle in two years; Perennial: completes full life cycle in more than two years; Perennial-Evergreen: perennial and retains functional leaves throughout the year; Perennial-Deciduous: perennial and loses all leaves synchronously for part of the year. Some species are variable (for example, have annual and biennial individuals), and in those cases we attempted to characterize the most common modality.
Because of the paucity of data available for chromists and fungi, we did not extract life history trait data for the relatively few species in these taxonomic groups.
Data visualization and summary. We used the R-package raster (v3.1-5) to visualize the spatial extent of the data in CaliPopGen in Fig. 3. Panel (A) shows a summary plot of all unique populations of both the Population Genetic Diversity in Dataset 1 21 and the Pairwise Population Differentiation in Dataset 2 21 . Panel (B) shows the total number of unique populations in each California terrestrial ecoregion. Panel (C) depicts all data entries of Population Genetic Diversity Dataset 1 21 , summed for each 20x20 km grid cell. Panel (D) shows the density of pairwise straight lines drawn between pairs of localities in the Pairwise Population Differentiation Dataset 2 21 , depicted as the total number of lines per 20x20 km grid cell. The number of populations and species of both Datasets 1 21 & 2 21 are summarized for each marine and terrestrial ecoregion in Table 6.

Column Name Description
CitationID Unique ID assigned to each source article.

EntryID
Unique ID assigned to each unique entry in the entire CaliPopGen database When multiple analyses were conducted, the level by which data were split is noted here (e.g. species or sampling year)

Pop1ComparisonCharacteristic
Characteristic of special comparison

Pop2ComparisonCharacteristic
Characteristic of special comparison Pop1YearStart First year of sample collection Pop1YearEnd Last year of sample collection Pop2YearStart First year of sample collection Pop2YearEnd Last year of sample collection  Table 3) Columns 75-101 Plant life history data (see Table 4) Table 2. Description of the pairwise genetic distance data in Dataset 2 21 . www.nature.com/scientificdata www.nature.com/scientificdata/

Data Records
The CaliPopGen database comprises four datasets, which are hosted at Figshare and can be downloaded as XLSX, TSV and CSV files. For convenience, the life history trait data for both animals (Dataset 3 21 ) and plants (Dataset 4 21 ) have also been included in Dataset 1 21 and Dataset 2 21 . We combined the genetic and life history data under the assumption that potential users may want to examine correlations between these two classes of data. Dataset 1 21 : The Population Genetic Diversity dataset consists of 101 columns, described in Table 1, and is comprised of data from 401 studies on 446 (sub-)species and 4,697 unique species-population-marker type combinations, with the latter equaling the number of rows in the dataset. The first 31 columns summarize taxonomic, population, marker type, and genetic data, while the remaining 70 columns contain data on animal and plant life history (Dataset 3 21 and Dataset 4 21 , respectively, see below). Dataset 2 21 : The Pairwise Population Differentiation dataset consists of 106 columns, described in Table 2, and is comprised of data from 199 studies on 197 (sub-)species and 14,703 pairwise population comparisons, with the latter equaling the number of rows. The first 36 columns summarize taxonomic, population, marker type, and pairwise population comparison data, while the remaining 70 columns contain data on animal and plant life history (Dataset 3 21 and Dataset 4 21 , respectively, see below). Dataset 3 21 : The Animal Life History Traits dataset consists of 45 columns, containing data for 432 species and subspecies, and is described in Table 3. The first 29 columns describe the life history of species and subspecies, and give details on their conservation status, while columns 30-45 provide information on the sources of these data. Dataset 4 21 : The Plant Life History Traits dataset consists of 37 columns containing data for 177 species and is described in Table 4. The first 23 columns describe the life-history of species and subspecies, and give details on their conservation status, while columns 24-37 provide sources of data. Total species numbers of Dataset 3 21 & 4 21 are higher than the number of species of Dataset 1 21 & 2 21 because we left species in these datasets even though their genetic entries may have been excluded based on the criteria set out in the Methods.  Table 3. Description of the animal life-history data in Dataset 3 21 . www.nature.com/scientificdata www.nature.com/scientificdata/ Technical Validation Article classification. During the first step in our screening protocol based on titles and abstracts (see Fig. 4), we examined the repeatability (intra-individual variation), and reproducibility (inter-individual variation) of article classification. Given that multiple individuals were doing the article screening, we recognize that understanding variation at both of these levels is important. During this first screening, six screeners assigned a non-overlapping set of articles into three broad categories ("reject", "include", or "possibly include"), based on our five screening criteria (see Methods); we used "possibly include" if it was unclear from the title and abstract if a paper contained appropriate data. Each screener independently evaluated 777-782 articles (total screened = 4,942). To quantify the repeatability of our screeners, all of whom were UCLA undergraduates, each individual re-screened a subset of their original set of articles. 54 randomly selected papers were re-screened  G1/S1 Critically imperiled; at very high risk of extirpation in the jurisdiction due to very restricted range, very few populations or occurrences, very steep declines, severe threats, or other factors. G2/S2 Imperiled; at high risk of extirpation in the jurisdiction due to restricted range, few populations or occurrences, steep declines, severe threats, or other factors.

G3/S3
Vulnerable; at moderate risk of extirpation in the jurisdiction due to a fairly restricted range, relatively few populations or occurrences, recent and widespread declines, threats, or other factors.

G4/S4
Apparently secure; at a fairly low risk of extirpation in the jurisdiction due to an extensive range and/or many populations or occurrences, but with possible vii cause for some concern as a result of local recent declines, threats, or other factors.

G5/S5
Secure; at very low or no risk of extirpation in the jurisdiction due to a very extensive range, abundant populations or occurrences, with little to no concern from declines or threats. www.nature.com/scientificdata www.nature.com/scientificdata/ by the same person (6 screeners, range 6-13 papers per person, mean = 10.8 papers re-screened/screener). We allowed 10 weeks between the initial and re-screening procedures, which all screeners felt was a sufficiently long time that they would not remember their initial classification, and papers were randomly chosen by the senior authors. To quantify the reproducibility of the screening process across individuals, 421 papers were re-screened by a different individual than the original screener (8 re-screeners, range 46-60 per person, mean = 50.33 re-screened/screener). Each of the 421 papers was re-screened by exactly one new person. This procedure included JB and EMT in addition to the original six undergraduates.
As might be expected, intra-individual repeatability (agreement between the initial and re-screened classification of a paper screened by the same person) was higher than inter-individual reproducibility (agreement between the initial and re-screened classification of a paper screened by two different people): 92.6% (50/54) of papers re-screened by the same individual received an identical score whereas 74.8% (315/421) of papers re-screened by a different individual received an identical score. Across both of these exercises, 17.5% of articles that were re-screened by either the same or different individual (total = 475) were assigned to different categories between the first and second screening. For the inter-individual analyses, 27.7% of "possibly include" articles changed status when screened by different individuals, while only 16.0% of "reject" and 16.3% of "include" decisions changed. However, when we subsequently attempted to extract data during the Data Collection Phase, we did so from both "include" and "possibly include" papers, so the relatively low change of "reject" status makes us comfortable that screener variability and its potential bias had at most a very limited impact in our decision pipeline.

Data validation.
To identify and correct potential recording errors in the datasets after the initial round of data extraction, we flagged numerical outliers and values outside of theoretical expectations for all genetic parameters and life history traits. Both outliers and values outside theoretical bounds may represent values as reported in the original publication, or they may be transcription errors as we compiled datasets. To increase the likelihood of identifying errors via outlier analysis, we examined each genetic parameter distribution separately for each marker type and taxonomic group (for example, H O of microsatellite markers in Aves was examined separately from H O of mitochondrial markers in Reptilia), and we examined life history trait distributions separately for each taxonomic group. In all cases, we identified outliers as values greater or less than the upper or lower quartiles +/− 1.5 * IQR (IQR = inter-quartile range), using the function boxplot.stats in the R-package grDevices. For all identified outliers we returned to the original source publication to confirm that values were as reported, or corrected them if they were a data-entry error. Correctly transcribed values falling outside of their theoretical bounds (H E , H O , π, F ST , G ST , D ST are constrained between zero and 1, F IS is bounded by −1 and 1, N e must be greater than zero) were left unaltered, which users of the CaliPopGen databases should consider carefully in using these results.  Table 6. Summary of total numbers of populations and species per California ecoregion, separately for population genetic and pairwise datasets. The first three are marine, followed by the 19 USDA-defined ecoregions.