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

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

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

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

The genome assembly of Island Oak (Quercus tomentella), a relictual island tree species.

Island oak (Quercus tomentella) is a rare relictual island tree species that exists only on six islands off the coast of California and Mexico, but was once widespread throughout mainland California. Currently, this species is endangered by threats such as non-native plants, grazing animals, and human removal. Efforts for conservation and restoration of island oak currently underway could benefit from information about its range-wide genetic structure and evolutionary history. Here we present a high-quality genome assembly for Q. tomentella, assembled using PacBio HiFi and Omni-C sequencing, developed as part of the California Conservation Genomics Project (CCGP). The resulting assembly has a length of 781 Mb, with a contig N50 of 22.0 Mb and a scaffold N50 of 63.4 Mb. This genome assembly will provide a resource for genomics-informed conservation of this rare oak species. Additionally, this reference genome will be the first one available for a species in Quercus section Protobalanus, a unique oak clade present only in western North America.

A highly contiguous genome assembly for the pocket mouse Perognathus longimembris longimembris.

The little pocket mouse, Perognathus longimembris, and its nine congeners are small heteromyid rodents found in arid and seasonally arid regions of Western North America. The genus is characterized by behavioral and physiological adaptations to dry and often harsh environments, including nocturnality, seasonal torpor, food caching, enhanced osmoregulation, and a well-developed sense of hearing. Here we present a genome assembly of Perognathus longimembris longimembris generated from PacBio HiFi long read and Omni-C chromatin-proximity sequencing as part of the California Conservation Genomics Project. The assembly has a length of 2.35 Gb, contig N50 of 11.6 Mb, scaffold N50 of 73.2 Mb, and includes 93.8% of the BUSCO Glires genes. Interspersed repetitive elements constitute 41.2% of the genome. A comparison with the highly endangered Pacific pocket mouse, P. l. pacificus, reveals broad synteny. These new resources will enable studies of local adaptation, genetic diversity, and conservation of threatened taxa.

The reference genome assembly of the bright cobblestone lichen, Acarospora socialis.

Acarospora socialis, the bright cobblestone lichen, is commonly found in southwestern North America. This charismatic yellow lichen is a species of key ecological significance as it is often a pioneer species in new environments. Despite their ecological importance virtually no research has been conducted on the genomics of A. socialis. To address this, we used long-read sequencing to generate the first high-quality draft genome of A. socialis. Lichen thallus tissue was collected from Pinkham Canyon in Joshua Tree National Park, California and deposited in the UC Riverside herbarium under accession #295874. The de novo assembly of the mycobiont partner of the lichen was generated from Pacific Biosciences HiFi long reads and Dovetail Omni-C chromatin capture data. After removing algal and bacterial contigs, the fungal genome was approximately 31.2 Mb consisting of 38 scaffolds with contig and scaffold N50 of 2.4 Mb. The BUSCO completeness score of the assembled genome was 97.5% using the Ascomycota gene set. Information on the genome of A. socialis is important for California conservation purposes given that this lichen is threatened in some places locally by wildfires due to climate change. This reference genome will be used for understanding the genetic diversity, population genomics, and comparative genomics of A. socialis species. Genomic resources for this species will support population and landscape genomics investigations, exploring the use of A. socialis as a bioindicator species for climate change, and in studies of adaptation by comparing populations that occur across aridity gradients in California.

A reference genome assembly for the continentally distributed ring-necked snake, Diadophis punctatus.

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.

De novo assembly of a chromosome-level reference genome for the California Scrub-Jay, Aphelocoma californica.

We announce the assembly of the first de novo reference genome for the California Scrub-Jay (Aphelocoma californica). The genus Aphelocoma comprises four currently recognized species including many locally adapted populations across Mesoamerica and North America. Intensive study of Aphelocoma has revealed novel insights into the evolutionary mechanisms driving diversification in natural systems. Additional insights into the evolutionary history of this group will require continued development of high-quality, publicly available genomic resources. We extracted high molecular weight genomic DNA from a female California Scrub-Jay from northern California and generated PacBio HiFi long-read data and Omni-C chromatin conformation capture data. We used these data to generate a de novo partially phased diploid genome assembly, consisting of two pseudo-haplotypes, and scaffolded them using inferred physical proximity information from the Omni-C data. The more complete pseudo-haplotype assembly (arbitrarily designated "Haplotype 1") is 1.35 Gb in total length, highly contiguous (contig N50 = 11.53 Mb), and highly complete (BUSCO completeness score = 97%), with comparable scaffold sizes to chromosome-level avian reference genomes (scaffold N50 = 66.14 Mb). Our California Scrub-Jay assembly is highly syntenic with the New Caledonian Crow reference genome despite ~10 million years of divergence, highlighting the temporal stability of the avian genome. This high-quality reference genome represents a leap forward in publicly available genomic resources for Aphelocoma, and the family Corvidae more broadly. Future work using Aphelocoma as a model for understanding the evolutionary forces generating and maintaining biodiversity across phylogenetic scales can now benefit from a highly contiguous, in-group reference genome.

A highly contiguous reference genome for the Steller's jay (Cyanocitta stelleri).

The Steller's jay is a familiar bird of western forests from Alaska south to Nicaragua. Here, we report a draft reference assembly for the species generated from PacBio HiFi long-read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 352 scaffolds totaling 1.16 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 7.8 Mb, scaffold N50 of 25.8 Mb, and BUSCO completeness score of 97.2%. Repetitive elements span 16.6% of the genome including nearly 90% of the W chromosome. Compared with high-quality assemblies from other members of the family Corvidae, the Steller's jay genome contains a larger proportion of repetitive elements than 4 crow species (Corvus), but a lower proportion of repetitive elements than the California scrub-jay (Aphelocoma californica). This reference genome will serve as an essential resource for future studies on speciation, local adaptation, phylogeography, and conservation genetics in this species of significant biological interest.

A genome assembly for the southern Pacific rattlesnake, Crotalus oreganus helleri, in the western rattlesnake species complex.

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.

A draft reference genome assembly of the Pipevine Swallowtail butterfly, Battus philenor hirsuta.

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.

Chromosome-level reference genome of stinkwort, Dittrichia graveolens (L.) Greuter: A resource for studies on invasion, range expansion, and evolutionary adaptation under global change.

Dittrichia graveolens (L.) Greuter, or stinkwort, is a weedy annual plant within the family Asteraceae. The species is recognized for the rapid expansion of both its native and introduced ranges: in Europe, it has expanded its native distribution northward from the Mediterranean basin by nearly 7 °C latitude since the mid-20th century, while in California and Australia the plant is an invasive weed of concern. Here, we present the first de novo D. graveolens genome assembly (1N = 9 chromosomes), including complete chloroplast (151,013 bp) and partial mitochondrial genomes (22,084 bp), created using Pacific Biosciences HiFi reads and Dovetail Omni-C data. The final primary assembly is 835 Mbp in length, of which 98.1% are represented by 9 scaffolds ranging from 66 to 119 Mbp. The contig N50 is 74.9 Mbp and the scaffold N50 is 96.9 Mbp, which, together with a 98.8% completeness based on the BUSCO embryophyta10 database containing 1,614 orthologs, underscores the high quality of this assembly. This pseudo-molecule-scale genome assembly is a valuable resource for our fundamental understanding of the genomic consequences of range expansion under global change, as well as comparative genomic studies in the Asteraceae.

Reference genome of California walnut, Juglans californica, and resemblance with other genomes in the order Fagales.

Juglans californica, California walnut, is a vulnerable small tree that is locally abundant but restricted to woodland and chaparral habitats of Southern California threatened by urbanization and land use change. This species is the dominant species in a unique woodland ecosystem in California. It is one of 2 endemic California walnut species (family Juglandaceae). The other species, Northern California black walnut (J. hindsii), has been suggested controversially to be a variety of J. californica. Here, we report a new, chromosome-level assembly of J. californica as part of the California Conservation Genomics Project (CCGP). Consistent with the CCGP common methodology across ~150 genomes, we used Pacific Biosciences HiFi long reads and Omni-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 137 scaffolds spanning 551,065,703 bp, has a contig N50 of 30 Mb, a scaffold N50 of 37 Mb, and BUSCO complete score of 98.9%. Additionally, the mitochondrial genome has 701,569 bp. In addition, we compare this genome with other existing high-quality Juglans and Quercus genomes, which are in the same order (Fagales) and show relatively high synteny within the Juglans genomes. Future work will utilize the J. californica genome to determine its relationship with the Northern California walnut and assess the extent to which these 2 endemic trees might be at risk from fragmentation and/or climate warming.

Assembly of the largest squamate reference genome to date: The western fence lizard, Sceloporus occidentalis.

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.

Reference genome of an iconic lizard in western North America, Blainville's horned lizard Phrynosoma blainvillii.

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.

A highly contiguous genome assembly for the California quail (Callipepla californica).

The California quail (Callipepla californica) is an iconic native bird of scrub and oak woodlands in California and the Baja Peninsula of Mexico. Here, we report a draft reference assembly for the species generated from PacBio HiFi long read and Omni-C chromatin-proximity sequencing data as part of the California Conservation Genomics Project (CCGP). Sequenced reads were assembled into 321 scaffolds totaling 1.08 Gb in length. Assembly metrics indicate a highly contiguous and complete assembly with a contig N50 of 5.5 Mb, scaffold N50 of 19.4 Mb, and BUSCO completeness score of 96.5%. Transposable elements (TEs) occupy 16.5% of the genome, more than previous Odontophoridae quail assemblies but in line with estimates of TE content for recent long-read assemblies of chicken and Peking duck. Together these metrics indicate that the present assembly is more complete than prior reference assemblies generated for Odontophoridae quail. This reference will serve as an essential resource for studies on local adaptation, phylogeography, and conservation genetics in this species of significant biological and recreational interest.

A draft reference genome of the Vernal Pool Fairy Shrimp, Branchinecta lynchi.

We present the reference genome of the Vernal Pool Fairy Shrimp Branchinecta lynchi. This branchiopod crustacean is endemic to California's freshwater ephemeral ponds. It faces enormous habitat loss and fragmentation as urbanization and agriculture have fundamentally changed the vernal pool landscape over the past 3 centuries. The assembled genome consists of 22 chromosome-length scaffolds that account for 96.85% of the total sequence. One hundred and ninety-five unscaffolded contigs comprise the rest of the genome's 575.6 Mb length. The genome is substantially complete with a BUSCO score of 90.0%. There is no immediately identifiable sex chromosome, typical for this class of organism. This new resource will permit researchers to better understand the adaptive capacity of this imperiled species, as well as answer lingering questions about anostracan physiology, sex determination, and development.

A chromosome-level reference genome for the Versatile Fairy Shrimp, Branchinecta lindahli.

We present the novel reference genome of the Versatile Fairy Shrimp, Branchinecta lindahli. The Versatile Fairy Shrimp is a freshwater anostracan crustacean found across the western United States from Iowa to Oregon and from Alberta to Baja California. It is an ephemeral pool specialist, living in prairie potholes, irrigation ditches, tire treads, vernal pools, and other temporary freshwater wetlands. Anostracan fairy shrimp are facing global declines with 3 species in California on the Endangered Species list. This species was included in the California Conservation Genomics Project to provide an easily accessible reference genome, and to provide whole-genome resources for a generalist species, which may lead to new insights into Anostracan resiliency in the face of climate change. The final gapped genome comprises 15 chromosome-length scaffolds covering 98.63% of the 384.8 Mb sequence length, and an additional 55 unscaffolded contigs.

Reference genome of the Monkeyface Prickleback, Cebidichthys violaceus.

Pricklebacks (Family Stichaeidae) are generally cold-temperate fishes most commonly found in the north Pacific. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Monkeyface Prickleback, Cebidichthys violaceus, to establish a genomic model for understanding phylogeographic patterns of marine organisms in California. These patterns, in turn, may inform the design of marine protected areas using dispersal models based on forthcoming population genomic data. The genome of C. violaceus is typical of many marine fishes at less than 1 Gb (genome size = 575.6 Mb), and our assembly is near-chromosome level (contig N50 = 1 Mb, scaffold N50 = 16.4 Mb, BUSCO completeness = 93.2%). Within the context of the CCGP, the genome will be used as a reference for future whole genome resequencing projects, enhancing our knowledge of the population structure of the species and more generally, the efficacy of marine protected areas as a primary conservation tool across California's marine ecosystems.

Reference genome of the Woolly Sculpin, Clinocottus analis.

Sculpins (Family Cottidae) are generally cold-temperate intertidal reef fishes most commonly found in the North Pacific. As part of the California Conservation Genomics Project (CCGP), we sequenced the genome of the Woolly Sculpin, Clinocottus analis, to establish a genomic model for understanding phylogeographic structure of inshore marine taxa along the California coast. These patterns, in turn, should further inform the design of marine protected areas using dispersal models based on genomic data. The small genome of C. analis is typical of marine fishes at less than 1 Gb (genome size = 538 Mb), and our assembly is near-chromosome level (contig N50 = 9.1 Mb, scaffold N50 = 21 Mb, BUSCO completeness = 97.9%). Within the context of the CCGP, the Woolly Sculpin genome will be used as a reference for future whole-genome resequencing projects aimed at enhancing our knowledge of the population structure of the species, and efficacy of marine protected areas across the state.

A draft reference genome of the red abalone, Haliotis rufescens, for conservation genomics.

Red abalone, Haliotis rufescens, are herbivorous marine gastropods that primarily feed on kelp. They are the largest and longest-lived of abalone species with a range distribution in North America from central Oregon, United States, to Baja California, MEX. Recently, red abalone have been in decline as a consequence of overharvesting, disease, and climate change, resulting in the closure of the commercial fishery in the 1990s and the recreational fishery in 2018. Protecting this ecologically and economically important species requires an understanding of their current population dynamics and connectivity. Here, we present a new red abalone reference genome as part of the California Conservation Genomics Project (CCGP). Following the CCGP genome strategy, we used Pacific Biosciences HiFi long reads and Dovetail Omni-C data to generate a scaffold-level assembly. The assembly comprises 616 scaffolds for a total size of 1.3 Gb, a scaffold N50 of 45.7 Mb, and a BUSCO complete score of 97.3%. This genome represents a significant improvement over a previous assembly and will serve as a powerful tool for investigating seascape genomic diversity, local adaptation to temperature and ocean acidification, and informing management strategies.