Genomics can play important roles in biodiversity conservation, especially for Extinct-in-the-Wild species where genetic factors greatly influence risk of total extinction and probability of successful reintroductions. The Christmas Island blue-tailed skink (Cryptoblepharus egeriae) and Lister's gecko (Lepidodactylus listeri) are two endemic reptile species that went extinct in the wild shortly after the introduction of a predatory snake. After a decade of management, captive populations have expanded from 66 skinks and 43 geckos to several thousand individuals; however, little is known about patterns of genetic variation in these species. Here, we use PacBio HiFi long-read and Hi-C sequencing to generate highly contiguous reference genomes for both reptiles, including the XY chromosome pair in the skink. We then analyse patterns of genetic diversity to infer ancient demography and more recent histories of inbreeding. We observe high genome-wide heterozygosity in the skink (0.007 heterozygous sites per base-pair) and gecko (0.005), consistent with large historical population sizes. However, nearly 10% of the blue-tailed skink reference genome falls within long (>1 Mb) runs of homozygosity (ROH), resulting in homozygosity at all major histocompatibility complex (MHC) loci. In contrast, we detect a single ROH in Lister's gecko. We infer from the ROH lengths that related skinks may have established the captive populations. Despite a shared recent extinction in the wild, our results suggest important differences in these species' histories and implications for management. We show how reference genomes can contribute evolutionary and conservation insights, and we provide resources for future population-level and comparative genomic studies in reptiles.

Understanding the evolutionary consequences of anthropogenic change is imperative for estimating long-term species resilience. While contemporary genomic data can provide us with important insights into recent demographic histories, investigating past change using present genomic data alone has limitations. In comparison, temporal genomics studies, defined herein as those that incorporate time series genomic data, utilize museum collections and repeated field sampling to directly examine evolutionary change. As temporal genomics is applied to more systems, species and questions, best practices can be helpful guides to make the most efficient use of limited resources. Here, we conduct a systematic literature review to synthesize the effects of temporal genomics methodology on our ability to detect evolutionary changes. We focus on studies investigating recent change within the past 200 years, highlighting evolutionary processes that have occurred during the past two centuries of accelerated anthropogenic pressure. We first identify the most frequently studied taxa, systems, questions and drivers, before highlighting overlooked areas where further temporal genomic studies may be particularly enlightening. Then, we provide guidelines for future study and sample designs while identifying key considerations that may influence statistical and analytical power. Our aim is to provide recommendations to a broad array of researchers interested in using temporal genomics in their work.

There is growing interest in the role of structural variants (SVs) as drivers of local adaptation and speciation. From a biodiversity genomics perspective, the characterization of genome-wide SVs provides an exciting opportunity to complement single nucleotide polymorphisms (SNPs). However, little is known about the impacts of SV discovery and genotyping strategies on the characterization of genome-wide SV diversity within and among populations. Here, we explore a near whole-species resequence data set, and long-read sequence data for a subset of highly represented individuals in the critically endangered kākāpō (Strigops habroptilus). We demonstrate that even when using a highly contiguous reference genome, different discovery and genotyping strategies can significantly impact the type, size and location of SVs characterized genome-wide. Further, we found that the mean number of SVs in each of two kākāpō lineages differed both within and across generations. These combined results suggest that genome-wide characterization of SVs remains challenging at the population-scale. We are optimistic that increased accessibility to long-read sequencing and advancements in bioinformatic approaches including multireference approaches like genome graphs will alleviate at least some of the challenges associated with resolving SV characteristics below the species level. In the meantime, we address caveats, highlight considerations, and provide recommendations for the characterization of genome-wide SVs in biodiversity genomic research.

For both undivided and subdivided populations, the consensus method to maintain genetic diversity is the Optimal Contribution (OC) method. For subdivided populations, this method determines the optimal contribution of each candidate to each subpopulation to maximize global genetic diversity (which implicitly optimizes migration between subpopulations) while balancing the relative levels of coancestry between and within subpopulations. Inbreeding can be controlled by increasing the weight given to within-subpopulation coancestry (λ). Here we extend the original OC method for subdivided populations that used pedigree-based coancestry matrices, to the use of more accurate genomic matrices. Global levels of genetic diversity, measured as expected heterozygosity and allelic diversity, their distributions within and between subpopulations, and the migration pattern between subpopulations, were evaluated via stochastic simulations. The temporal trajectory of allele frequencies was also investigated. The genomic matrices investigated were (i) the matrix based on deviations of the observed number of alleles shared by two individuals from the expected number under Hardy-Weinberg equilibrium; and (ii) a matrix based on a genomic relationship matrix. The matrix based on deviations led to higher global and within-subpopulation expected heterozygosities, lower inbreeding and similar allelic diversity than the second genomic and pedigree-based matrices when a relatively high weight was given to the within-subpopulation coancestries (λ ≥ 5). Under this scenario, allele frequencies moved only slightly away from the initial frequencies. Therefore, the recommended strategy is to use the former matrix in the OC methodology giving a high weight to the within-subpopulation coancestry.