State of Shark and Ray Genomics in an Era of Extinction

Chondrichthyan species (sharks, rays, skates, and chimeras) are a class of high ecological, economic, and cultural significance, and yet they are the most threatened taxa in the marine environment. The creation of reference chromosome-length genome assemblies allows for conservation genomics methods, such as population and ecological genomics, to be utilized. Despite being greatly threatened and of great importance in maintaining ecosystem function, chondrichthyan species have been repeatedly absent from conservation-based genome sequencing projects. Less than 1% of these species have a genome sequence, despite their almost 50% either threatened or Data Deficient conservation status. Most notably, there are seven orders within this class without any genome representation. In this review, we identify gaps in chondrichthyan genomic resources and demonstrate how the lack of genomic resources for this major taxonomic class is limiting the conservation of these already difficult to conserve species. We highlight other applications for chondrichthyans genomics, such as evolutionary and developmental biology. Likely, the mismatching sampling protocols and limited computational skills and communication between fields have been preventing the integration of marine and molecular sciences. Here, we propose that this field is in dire need to move forward quickly to increase protection for marine species and ecosystems through improved collaboration between marine, molecular, and computer sciences.


INTRODUCTION
The marine environment is vital in providing a multitude of ecosystem services and is responsible for the function of our global climate (Barbier, 2017;Gattuso et al., 2018). Marine ecosystems display extensive signs of human interference that are reflected at multiple levels, ranging from primary producers that provide habitat (e.g., seagrasses and corals) to apex predators like sharks (Hoegh-Guldberg and Bruno, 2010;Smale et al., 2019). Cumulative pressures and impacts of multiple stressors (Sequeira et al., 2019) are likely contributing to species extinction rates which have been rapidly increasing far above historical averages, and exceeding threshold boundaries, resulting in the sixth mass extinction in the earth's history (Rockstrom, 2009;Ceballos et al., 2015). The conservation of species is of critical importance to maintain Earth's biodiversity and its ecological, economic, and social benefits (Hoffmann et al., 2015;Shafer et al., 2015). Novel techniques will be required to assist environmental stewardship in the face of increasing threats.
Emerging genomic technologies are powerful tools to assist conservation research (Corlett, 2017). For instance, ecological and environmental genomics uses genomic data to investigate relationships between organisms and their environment to predict species' responses and adaptions to environmental change (Matz, 2018;Ruegg et al., 2018). However, these methods require high-quality reference genomes to compare study-specific genetic and genomic data against to have accurate results (Shafer et al., 2015). Genomics technology has recently become widely accessible to a range of fields outside the traditional molecular sciences through increased power and affordability (Goodwin et al., 2016;Worley et al., 2017). Historically, genome sequencing involved generating many short reads (contigs) and piecing them together based on sequence overlap (contigging), where longer sequence reads allow more overlap (Schatz et al., 2010). These larger contiguous sequences can then be linked to create scaffolds and ordered according to the predicted genome structure (Dudchenko et al., 2017). However, due to regions of high repeat content in genomes or low sequencing coverage, consequently resulting in gaps between contigs, within and between scaffolds, most genome assemblies are highly fragmented (Kadota et al., 2020). Recent developments in long-read sequencing technology, such as Pacific Biosciences (PacBio) or the Oxford Nanopore, have made advancements in overcoming this issue by generating very long contigs initially (Lang et al., 2020). This creates higher sequence overlap and therefore less gaps with higher genome completeness.
Novel high-throughput chromosome conformation capture (Hi-C) technology, combined with contigging, can generate high-quality full chromosome-length genome assemblies in a relatively inexpensive and quick manner (Dudchenko et al., 2017). Hi-C data are produced based on the principle that the closer two DNA locations are to each other, the more frequently they will come into contact and so a proximity matrix of DNA locations can be created. Then, through a series of bioinformatics pipelines, chromosome-length scaffolds can be produced (Lieberman-Aiden et al., 2009;Rao et al., 2014;Dudchenko et al., 2018). The storage and analysis of biological data using computer sciences, termed "bioinformatics, " is a fastdeveloping field in molecular sciences with its ability to deal with big and complex data sets such as genetic/genomic data (Baxevanis et al., 2020). Bioinformatics pipelines can transfer data through multiple functions within and between software and programs. The combination of Hi-C data and bioinformatics pipelines has allowed chromosome-length genome assembly to become largely automated.
The addition of Hi-C data to traditional genome sequencing methods has allowed the production of high-quality genomes in a timely and cost-effective manner (Dudchenko et al., 2017). This development has led to multiple global projects aiming to sequence the genomes of Earth's biodiversity for a multitude of purposes. Ambitious consortiums like the Earth BioGenome Project (EBP) (Lewin et al., 2018), which aims to sequence the genomes of all eukaryotic species, work on a worldwide scale of networks and partnering organizations, such as the DNA Zoo (Hoencamp et al., 2021) and the Genome 10k/Vertebrate Genomes Project (VGP) (Koepfli et al., 2015). Despite almost 40 affiliated projects globally, there has been a heavy focus on terrestrial species, leading to a lack of representation of marine species. In particular, there is a severe lack of cartilaginous fish, such as sharks and rays (Class: Chondrichthyes), with less than 1% of these species sequenced, in spite of their high ecological, social, and economic importance (Figure 1).
Sharks and rays have vital roles in maintaining ecosystem health and function through top-down regulation (Yagnesh et al., 2020;Flowers et al., 2021). However, these chondrichthyan species have an extinction risk well above most other vertebrate classes on average with an estimated over one-third of species threatened by extinction (Dulvy et al., 2014. Despite their high-risk status and significance to ecosystem function, which should make these species high priority for conservation, an extensive gap in conservation-based genetic and genomic resources exists (Johri et al., 2019). Additionally, chondrichthyan species present opportunities for investigation into other potential applications, including biomedical research (Marra et al., 2019). The integration of genomics in the conservation and study of chondrichthyan species will provide more effective and efficient methods of protection through accurate population assessments (Attard et al., 2018;Hohenlohe et al., 2021). In this review, we discuss the state of chondrichthyan genomics research to date along with clear applications for conservation as well as other research fields. We indicate the obstacles currently hindering swift integration of molecular marine sciences in marine conservation and provide recommendations for a way forward through improved protocols and communication between fields.

THE STATE OF CHONDRICHTHYAN GENOMICS RESOURCES
The class of Chondrichthyes is split into two subclasses, one containing sharks, rays, and skates (Elasmobranchii) and the other chimeras (Holocephali). To date, only 11 chondrichthyan genomes have been sequenced out of ∼1,200 species, i.e., less than 1% (Figure 1). The first genome published for this class was of the Elephant "Shark" (Callorhinchus milii), which is in fact a chimera, member of the Holocephali subclass (Venkatesh et al., 2014). This was the first chondrichthyan genome chosen to be sequenced due to its small size [0.975 gigabase pairs (Gbp)] and potential to become a model organism for evolution and development studies (Venkatesh et al., 2005). Indeed, that study, which was published almost two decades ago, has prompted recent research into chondrichthyan genomes and their utilization in evolutionary and developmental biology as well as biomedical insights (Hara et al., 2018;Marra et al., 2019; Figure 1). Recently, a higher quality genome assembly for the elephant shark was made publicly available, however it is still only at scaffold level (NCBI, 2021a). The small-eyed rabbitfish (Hydrolagus affinis) is the only other species of chimera published (Fonseca et al., 2020).
Multiple genome assemblies have been published for a few chondrichthyan species, such as the whale shark and white shark, which demonstrate differences in sequencing technology and assembly techniques. A second de novo (from scratch) genome assembly for the whale shark was recently published (Weber et al., 2020) with a lower percentage of missing BUSCOs (Benchmarking Universal Single-Copy Orthologs), meaning more complete genes, despite the higher fragmentation compared to the older scaffold assembly for this species. The difference in these two assemblies is due to differences in the sequencing technologies used: Roche 454, Illumina, and Pacific Biosciences platforms were used originally (Read et al., 2017), compared with Illumina HiSeq 2500m platform (Weber et al., 2020) used in the recent assembly ( Table 1). The different quality of these genome assemblies highlights how both sequencing technology and assembly techniques can impact the final genome. A genome assembly with higher completeness can be utilized much more effectively and investigated deeper than a limited incomplete assembly (Hoffmann et al., 2015;Rhie and Project, 2021). For example, investigation into the white shark's ability to regulate cell growth to prevent tumors was hindered by the original incomplete genome, as a large section of the specific gene of interest was missing and therefore could not be investigated further (Marra et al., 2019). Therefore, the desired output quality should be considered during the sequencing process and technology choices made accordingly.
Frontiers in Marine Science | www.frontiersin.org genomes. Some of these methods may now be outdated and less accurate than current technology, or simply use different algorithms, and so produce varying results (Table 1). Notably, one of the major differences between these genomes is the level of assembly. While most recently released assemblies are at chromosome-length, earlier genomes, limited by the technology of the time, are only to a scaffold level and therefore more likely to be incomplete and fragmented. For example, the original 2019 scaffold genome assembly of the white shark is smaller, compared to the recent 2021 chromosome-length assembly, as it is less complete ( Table 1). This comparison between two assemblies of the same species demonstrates the difference quality can have for a reference genome. This highlights the issues associated with comparing genomes of different species of varying qualities, when lower quality assemblies are likely to limit downstream analysis due to missing fragments of genomic data and could end up impacting results.

GENOMICS RESOURCES AND CONSERVATION
Chondrichthyan species are an exceptionally difficult class for applying conservation strategies (Dulvy et al., 2017). This is likely due to their large geographical range, fishing pressure, and slow life-history traits such as reproduction rate. Conservation genomics has extensive applications to assist investigating and protecting challenging species like chondrichthyans, through studies such as population genomics which can accurately assess effective population sizes and structures (Johri et al., 2019;Oleksiak and Rajora, 2020). These methods, however, optimally require high-quality chromosome-length genome assemblies to function as important references from which to compare and build investigations (Matz, 2018;Cordier et al., 2020). For example, low coverage whole-genome resequencing (WGR) requires a species-specific reference assembly to align to individual's genomes, in order to assess relatedness and genetic structure within and between populations (Fuentes-Pardo and Ruzzante, 2017; Lou et al., 2020). Without a reference genome, other methods such as polymorphism analyses have lower resolution and less accurate results (Galla et al., 2019;Johri et al., 2019). Reference genomes are particularly useful for marine non-model species, such as chondrichthyans which are lacking prior genetic or genomic data, as it creates a baseline resource (Oleksiak and Rajora, 2020). Additionally, in the event of species extinction, a reference genome can assist with preserving genetic diversity and information for posterior studies. Hence, sequencing reference genomes for target conservation species provides the optimal and most comprehensive approach to conservation genomics (Oleksiak and Rajora, 2020). Significant advances in genomics have been made over the last decade, yet uptake in the conservation and management of species continues to be limited. Two branches of genomics are particularly relevant for conservation: ecological genomics and environmental genomics. Ecological aims to investigate the genome structure and function to understand the relationship between organisms and their environment. Environmental genomics aims to predict organism responses and adaptations to environmental changes (Martin, 2008;Ungerer et al., 2008;Van Straalen and Roelofs, 2011). Environmental and ecological genomics has the potential to detect species presence, local and potential future genetic adaptations, and responses to anthropogenic influences, such as climate change and ocean acidification (Ungerer et al., 2008;Savolainen et al., 2013;Ruegg et al., 2018;Cordier et al., 2020). Additionally, environmental genomics can investigate, monitor, and manage non-indigenous species in the marine environment through studying genetic traits that allow species to become highly adaptable and invasive (Bourne et al., 2018). Both offer high value to conservation efforts; environmental genomics has value for assessing species composition through presence detection, whereas ecological genomics is more valuable for investigating species responses and adaption through genome function. These applications for a single genome alone have vast potential for investigating adaptations or negative responses in chondrichthyan species, and other non-model organisms, to climate change and other threats (Matz, 2018). The inclusion of environmental and ecological genomics in applied conservation management, however, requires high-quality genomic resources (e.g., chromosome-length assemblies), again highlighting the need for their development.
With a chromosome-length genome assembly and the accompanying fine-scale genetic data it is also possible to track an individual's movements, population sizes (N e ), and kin relationships on a global scale (Oleksiak and Rajora, 2020;Hohenlohe et al., 2021). As many shark species are highly migratory and cosmopolitan, investigation at the global scale is essential to assess the effective protection of these highrisk species (Dulvy et al., 2017;Attard et al., 2018). Methods such as documenting genetic polymorphisms on a genome-scale, close-kin mark-recapture, and kinship analysis, have already been successful with marine species (Cipriano and Palumbi, 1999;Attard et al., 2018;Hillary et al., 2018). While genomewide polymorphism analyses do not necessarily require a reference genome, only with a high-quality genome can finescale resolution of genetic data be investigated, in particular adaptive genetic variation (Pazmiño et al., 2017;Galla et al., 2019;Oleksiak and Rajora, 2020). The ability to detect adaptive genetic variation, as opposed to only neutral variation, is one of the main advances of population genomics utilizing reference genomes over traditional genetic techniques (Funk et al., 2012;Attard et al., 2018;Hohenlohe et al., 2021). Additionally, the availability of a reference genome allows calculation of population N e based on linkage disequilibrium, which is more accurate and effective than other methods, and the reconstruction of demographic history (Hendricks et al., 2018;Hohenlohe et al., 2021). These techniques can be implemented in genomics-based ecosystem monitoring and informing management practices (Cordier et al., 2020;Hohenlohe et al., 2021) to assist the conservation of ecologically significant species and associated habitats.
Moreover, genomics' ability to monitor and assess populations leads to extensive applications in fisheries management and aquaculture (Bernatchez et al., 2017;Houston et al., 2020). Approximately one-third of chondrichthyan species are threatened due to overfishing as many are commercially exploited or caught accidently as by-catch (Dulvy et al., 2017. For instance, the elephant shark (chimera) (Finucci et al., 2021), brownbanded bamboo shark (Fahmi et al., 2020), whitespotted bamboo shark (Chen et al., 2007), catshark, smalltooth sawfish (Yan et al., 2021), little skate (Curtis and Sosebee, 2015), and thorny skate (Pennino et al., 2019) (all species with available genomes), are commercially fished species which require management to avoid overexploitation. Molecular methods are fast becoming a mainstream tool in fisheries management and have already shown underrepresented or unregulated catch of chondrichthyan species through DNA barcoding and analyzing genetic population structure (Vargas-Caro et al., 2017;Kuguru et al., 2018). Population genomics methods, such as utilizing reference genomes to estimate effective population size or genetic structure, could become critical in ensuring sustainable fisheries management. Therefore, the slow development and utilization of chondrichthyan genomics is currently limiting potential application in sustainable resource management.
There are 12 extant orders within the Elasmobranchii subclass (and one within Holocephali), yet only six have genome representation ( Table 2). If the purpose of mass genome sequencing projects is to comprehensively act as a seed bank of genomes for conservation purposes, then it is imperative that this gap be addressed. Chondrichthyans are a class of highly unique evolution, including many keystone species that exert top-down control on marine ecosystems, provide economic and cultural values, and are greatly threatened by a multitude of stressors on global, regional, and local scales (Dulvy et al., 2014(Dulvy et al., , 2017Stein et al., 2018). These factors alone should make them a target for biodiversity conservation regardless. Out of the ∼1,200 chondrichthyan species, 180 (14.7) are classified as Vulnerable, 124 (10.1%) are Endangered, and 89 (7.3%) are Critically Endangered on the International Union for Conservation of Nature (IUCN) Red List (IUCN Red List, 2021). Strikingly, 172 species (14%) within this major class are listed as Data Deficient (Table 2; IUCN Red List, 2021). This means that almost 50% of all chondrichthyan species are either threatened by extinction or are of unknown conservation status. In these cases, a reference genome would provide a baseline from which a multitude of analyses could be conducted and greatly assist in conservation status assessment utilizing methods previously discussed. Hence, conservation genomics has the potential to provide researchers and management with vast applications in the field of chondrichthyan conservation where there is limited baseline knowledge (Johri et al., 2019).
Despite the high percentage of threatened or Data Deficient species, there are many orders of elasmobranch species without genome representation ( Table 2). For instance, the order Squatiniformes (angel sharks -subclass Elasmobranchii) has a total of 22 species, of which 8 (36.4%) are Critically Endangered, 13 (59.1%) threatened overall, and 1 (4.5%) is Data Deficient ( Table 2), while none have genomic representation. Another order of Elasmobranchii, the Rhinopristiformes (sawfish and guitarfish), has only a single genome represented in the total of 64 species although 43.8% of these species are Critically Endangered, 71.9% threatened and 9.4% are Data Deficient ( Table 2). The contrast between the percentage of Endangered species and the proportion of species with genomes within this class is alarming, given the significant role genomic resources can have in conservation. For instance, in the event of species extinction, a reference genome can assist with preserving genetic diversity and information for posterior studies. Conservation genomics methods could assist in preventing these high-risk species from becoming extinct, through previously mentioned methods to assess effective population sizes, species distributions, adaptive capacity, and responses to stress, and incorporating these into management solutions such as genetic rescues (Attard et al., 2018;Chattopadhyay et al., 2019;Hohenlohe et al., 2021). For example, through genome-wide population analysis, the Critically Endangered Siamese crocodile was found to have extremely low genetic diversity and effective population size and allowed conservation efforts to pick individuals with the correct genetics for captive breeding programs and others that were reintroduced into the wild (Chattopadhyay et al., 2019).

CHONDRICHTHYAN EVOLUTIONARY AND DEVELOPMENTAL INVESTIGATIONS
Chondrichthyan species are distinctly unique in their evolution. They represent an ancient lineage with a high degree of evolutionary divergence between the two subclasses and different orders (refer to Figure 1 for taxonomic relationships). It is estimated the class diverged roughly 473 million years ago (MYA) from other Gnathostomes (jawed vertebrates), with the two subclasses, Holocephali (chimeras) and Elasmobranchii (sharks, rays, and skates), diverging roughly 399 MYA (Hedges and Kumar, 2009;Inoue et al., 2010;Chen et al., 2012;Licht et al., 2012;Betancur-R et al., 2015). Chondrichthyans are a very distinct class as they are the only vertebrate clade to possess a cartilaginous skeleton, rather than bone (Boisvert et al., 2019). They are also noted for their ability to lose and replaced teeth continuously throughout their lifespan, particularly in sharks and rays (Boisvert et al., 2019). The availability of multiple species' genome allows for the reconstruction of clade phylogenetic relationships and demographic history (Hohenlohe et al., 2021). Phylogenomics has risen as the optimal method for resolving the phylogeny of species, as simply using traditional gene analysis has led to inconsistency and ambiguity in the evolution of species (Johri et al., 2019). This has highlighted the need for genome-wide investigations to increase accuracy in these studies (Johri et al., 2019). The first chondrichthyan genome (elephant shark) was sequenced for evolutionary and developmental investigation (Venkatesh et al., 2005), but led to the surprising revelation that the elephant shark had a higher degree of sequence similarity to humans than teleost fish genomes. This is despite humans being more closely related to teleost (bony) fish (Venkatesh et al., 2007). The elephant shark genome also displayed the lowest evolutionary rate of any vertebrate species yet investigated, when compared to teleost and tetrapod (bony) species, which likely stems from their slow neutral nucleotide mutation rate and very little change to the organization of their genes (Venkatesh et al., 2007(Venkatesh et al., , 2014. The sequencing of the genome for this species alone provided insights into the applications for developing chondrichthyan genomic resources, as it has subsequently been used in a range of studies, including embryogenesis, evolutionary and developmental biology, and functional genetics (Venkatesh et al., 2014;Hara et al., 2018;Marra et al., 2019).
In these following recent investigations, chromosomelength chondrichthyan genome assemblies have allowed the reconstruction of ancestral vertebrate chromosomes, chromosomal rearrangements, and breakages as well as the evolution of immune responses that show potential for human disease treatments . Additionally, investigation of the small-eyed rabbitfish genome showed the evolution of the endocrine systems in early jawed vertebrates through a comparative analysis of nuclear receptors (NRs) in these chondrichthyan genomes (Fonseca et al., 2020). These investigations provide an example of how chondrichthyan genomes can be utilized to solve evolutionary mysteries into the evolution and development of vertebrates, jawed vertebrates, and chondrichthyan species . These studies show a glimpse of fundamental advances in the field of evolutionary biology that can be made from chondrichthyan genomes.
Chondrichthyan genomics also has applications in biomedical science and physiology, as demonstrated by the great white shark and whale shark genomes (Marra et al., 2019;Weber et al., 2020). Compared to other vertebrate species, both genomes displayed positive selection and enrichment for molecular adaptations relating to genome stability and wound healing (Marra et al., 2019). Additionally, they had markers for cancer resistance through enhanced DNA damage response and repair mechanisms, genome, and telomere maintenance (Marra et al., 2019;Weber et al., 2020). Whale sharks possess expanded neurodegeneration genes, among other neural connectivity genes, again illustrating chondrichthyan evolutionary resistance (Weber et al., 2020). Without the addition of chondrichthyan species to a wide range of other classes, it would not be possible to assess physiological trends across all animals. For example, comparative genomics revealed that body weight, lifespan, and metabolic rate are all positively correlated across all vertebrates (Weber et al., 2020). This clearly demonstrates the broader context of chondrichthyan species as an important component to address overarching questions across all vertebrate classes.
As demonstrated above, a high-quality chromosome-length scaffold assembly can have a multitude of applications spanning various disciplines (Lewin et al., 2018;Kadota et al., 2020). A single species genome can serve as a baseline for research into various aspects of shark biology, physiology, and therefore conservation, which can continue to be built upon and improved. Genomics studies of chondrichthyan species have contributed new knowledge to many fields of research and have provided a foundation for further in-depth study. An abundance of different research areas could utilize these reference genomes, including but not limited to conservation genetics, transcriptomics, proteomics, comparative genetics/genomics, functional genetics/genomics, and evolution and development studies (Worley et al., 2017;Lewin et al., 2018).

OVERCOMING OBSTACLES IN THE INTEGRATION OF MOLECULAR AND MARINE SCIENCES
To develop genomic resources efficiently and effectively for chondrichthyan species, there are necessary obstacles to overcome to achieve the best research and conservation outcomes. Genome assemblies at the chromosome length are typically of higher quality and completeness (Dudchenko et al., 2017;Kadota et al., 2020;Rhie and Project, 2021). Therefore, they are likely to be more informative in downstream analyses, such as conservation and population genomics. Chondrichthyan genomes, particularly larger ones such as that of the catsharks (>4 Gbp), require long-read sequencing (e.g., Pacific Biosciences HiFi reads) to produce an accurate chromosome-length assembly due to the high level of heterozygosity and repeat content (Kelley et al., 2016;Rhie and Project, 2021). High percentages of repetitive elements throughout a genome lead to a higher number of gaps, making chromosome-length scaffolds extremely difficult to assemble (Rhie and Project, 2021). Although utilizing long-read sequencing technology assists in overcoming this issue, it is not a perfect solution due to the higher degree of sequence inaccuracy. Scaffold gaps remain unavoidable, meaning that accurate genome assembly for these species requires manual curation, in addition to automated pipelines (Rhie and Project, 2021). This presents a challenge to developing genomic resources for this at-risk class of species as both long-read sequencing and manual curation are expensive and time intensive. Long-range sequencing remains the ideal approach for chondrichthyan species, as this sequencing technology will continue to increase in accuracy and decrease in cost. Developing this as a standard practice should be a priority to maximize effectiveness of genomic resources.
Another major challenge for developing molecular resources of any kind for this class is obtaining suitable species samples. Many chondrichthyan species are notoriously difficult to collect samples from, mostly due to their elusiveness and habitat choice making them hard to find (Kelley et al., 2016). For example, chimeras inhabit the deep sea, while many pelagic species of sharks spend their time in the open ocean. For cosmopolitan shark species often found near-shore, such as the great white shark, another point of difficulty is safety while collecting biopsy samples. Conversely, Hammerhead shark species exhibit high capture stress and post-release mortality, making sampling a delicate task (Drymon and Wells, 2017). These factors make collecting fresh tissue samples, necessary for genome sequencing, a complicated, resource-intensive, and time-consuming activity. We recommend that all future opportunities where sharks are already being handled, for research or government tagging purposes, collect samples suitable for genetic study, preferably in a standardized approach according to molecular practices.
Obtaining suitable samples for genomic DNA extraction, Hi-C preparation, and sequencing from previously collected samples can also be problematic. Extracting RNA, which is necessary for proper genome annotation, essentially requires completely undegraded and uncontaminated DNA samples. Historically in biological sciences, biopsy samples were preserved in ethanol. However, ethanol prevents Hi-C techniques from detecting DNA contacts and so these samples are often unusable for sequencing genomes. Ideally, for a good quality chromosome-length genome assembly, tissue samples must be flash frozen at −80 • C or liquid nitrogen. Although freezing samples has become more popular in recent years, storage in ethanol or at −20 • C is common, resulting in many samples with DNA that is still too degraded to utilize for genomics as high molecular weight DNA is essential (Matz, 2018). Biological sciences in general, particularly in marine sciences, must move toward biobanking and storing species samples in accordance with standardized molecular techniques to improve sample quality critical for developing genomic resources (Wong et al., 2012;Rhie and Project, 2021).
An emerging challenge lies not only in generating genomic data but also in processing it (Dupont et al., 2007;Schweizer et al., 2021). The advances in sequencing technology and lack of computational software for bioinformatics particularly has led to a "bottleneck, " where data are being produced much faster than can be analyzed (Corlett, 2017;Hohenlohe et al., 2021;Schweizer et al., 2021). Additionally, genomics realistically requires high-performance computing to process and analyze full genomes, often gigabytes of data. Many bioinformatics pipelines and programs require a reasonably high level of computer programming knowledge, therefore outsourcing bioinformatic and computational assistance is a common way to overcome this issue (Corlett, 2017). Whilst there are increasing resources, such as workshops, seminars, and online courses, now available for researchers to upskill on emerging techniques, these are currently underutilized (Hendricks et al., 2018;Schweizer et al., 2021). Increasing the user-friendliness of bioinformatic tools would accelerate the integration of bioinformatics and computer sciences in marine sciences, which will be essential for the efficient development of genomic resources for species conservation (Hohenlohe et al., 2021).

DISCUSSION
This review highlights the limited genomic resources for chondrichthyan species, with less than 1% having a sequenced genome, as well as the prominent gaps, most notably the seven orders in this class lacking genome representation. Many of the few available genomes are only at a fragmented scaffold level that is insufficient for providing full potential for conservation of these highly threatened taxa. Chromosome-length genome assemblies have higher gene completeness and are therefore more informative in analysis compared to scaffolds, proving to be a more valuable resource. The development of genomic resources, in particular the creation of chromosome-length assemblies, can aid in the conservation and management of chondrichthyan species, such as through population and ecological genomics methods. So far, however, the development of chondrichthyan genomes has largely been for other purposes, such as the study of laboratory species (e.g., bamboo and catsharks), resulting in the majority of the available genomes for this class being species with conservation status of Least Concern or not Endangered. This is important due to the vast applications for chondrichthyan genomes beyond conservation, including evolutionary and developmental biology and biomedical sciences.
Despite the progress of chondrichthyan genomics over the past decade and the potential to quickly generate new knowledge that could directly aid in the conservation of these threatened species, genomic resources for this class remain very limited. There is still a long way to go to match the resources available for other marine clades, and even further to match terrestrial taxa. In molecular science, chondrichthyan species have been underutilized, and likewise, genomics technologies remain limited in marine sciences. Previous investigations and reviews have highlighted the great applications for non-model organism genomic data (Ellegren, 2014;Matz, 2018;Rhie and Project, 2021). Others have clearly demonstrated the need for more effective conservation strategies for sharks and rays, which remain difficult to manage (Dulvy et al., 2017;Stein et al., 2018). However, there has been very limited discussion of utilizing genomic data in the context of chondrichthyan species conservation (Johri et al., 2019).
This review identifies the urgent needs and the prioritization of orders most lacking in genomic resources that also have high percentages of Endangered species, such as the sawfish and guitarfish (Rhinopristiformes) and angel sharks (Squatiniformes). At minimum, baseline data should first be established to include a representative chromosome-length genome for every order within the Chondrichthyes class, considering how different they are from other vertebrate classes, and then for Critically Endangered species. This is a critical first step; however, the full applications for these genomes could have much further potential, such as comprehensive genomic population assessments of chondrichthyan species. This will assist with targets for conservation and increase management effectiveness and efficiency, with resulting conservation assessments able to inform policy decisions.
Going forward, it will be imperative to work across disciplinary boundaries to address emerging global issues like climate change and biodiversity loss. In particular, developing suitable sample collection and storage for genetic study and long-read sequencing as standard practice, as well as taking advantage of all opportunities for sampling will be vital in integrating these biological and molecular fields. The increase in bioinformatic resources becoming available is a promising sign of progress toward multidisciplinary utilization. If these resources become more user-friendly for other disciplines, then these molecular techniques that require high levels of bioinformatic skillsets and computational power for analyses will have higher utility in marine science. Increasing accessibility of molecular and bioinformatic skills will assist in decreasing the resistance and hesitation to fully integrating these fields across different biological sciences. If marine, molecular and computer sciences come together to develop genomic resources for chondrichthyan species, there is hope to conserve their high ecological significance and unique evolutionary history before extinction takes place.

AUTHOR CONTRIBUTIONS
JP: writing. JP, MF, AS, and PK: review and editing. MF, AS, and PK: supervision. All authors have read and agreed to the published version of the manuscript.