Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization.

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

Whale shark rhodopsin adapted to deep-sea lifestyle by a substitution associated with human disease.

Spectral tuning of visual pigments often facilitates adaptation to new environments, and it is intriguing to study the visual ecology of pelagic sharks with secondarily expanded habitats. The whale shark, which dives into the deep sea of nearly 2,000 meters besides near-surface filter feeding, was previously shown to possess the 'blue-shifted' rhodopsin (RHO), which is a signature of deep-sea adaptation. In this study, our spectroscopy of recombinant whale shark RHO mutants revealed that this blue shift is caused dominantly by an unprecedented spectral tuning site 94. In humans, the mutation at the site causes congenital stationary night blindness (CSNB) by reducing the thermal stability of RHO. Similarly, the RHO of deep-diving whale shark has reduced thermal stability, which was experimentally shown to be achieved by site 178 and 94. RHOs having the natural substitution at site 94 are also found in some Antarctic fishes, suggesting that the blue shift by the substitution at the CSNB site associated with the reduction in thermal stability might be allowed in cold-water deep-sea habitats.

A shark-inspired general model of tooth morphogenesis unveils developmental asymmetries in phenotype transitions.

Developmental complexity stemming from the dynamic interplay between genetic and biomechanic factors canalizes the ways genotypes and phenotypes can change in evolution. As a paradigmatic system, we explore how changes in developmental factors generate typical tooth shape transitions. Since tooth development has mainly been researched in mammals, we contribute to a more general understanding by studying the development of tooth diversity in sharks. To this end, we build a general, but realistic, mathematical model of odontogenesis. We show that it reproduces key shark-specific features of tooth development as well as real tooth shape variation in small-spotted catsharks Scyliorhinus canicula. We validate our model by comparison with experiments in vivo. Strikingly, we observe that developmental transitions between tooth shapes tend to be highly degenerate, even for complex phenotypes. We also discover that the sets of developmental parameters involved in tooth shape transitions tend to depend asymmetrically on the direction of that transition. Together, our findings provide a valuable base for furthering our understanding of how developmental changes can lead to both adaptive phenotypic change and trait convergence in complex, phenotypically highly diverse, structures.

A singular shark bitter taste receptor provides insights into the evolution of bitter taste perception.

Many animal and plant species synthesize toxic compounds as deterrent; thus, detection of these compounds is of vital importance to avoid their ingestion. Often, such compounds are recognized by taste 2 receptors that mediate bitter taste in humans. Until now, bitter taste receptors have only been found in bony vertebrates, where they occur as a large family already in coelacanth, a "living fossil" and the earliest-diverging extant lobe-finned fish. Here, we have revisited the evolutionary origin of taste 2 receptors (T2Rs) making use of a multitude of recently available cartilaginous fish genomes. We have identified a singular T2R in 12 cartilaginous fish species (9 sharks, 1 sawfish, and 2 skates), which represents a sister clade to all bony fish T2Rs. We have examined its ligands for two shark species, a catshark and a bamboo shark. The ligand repertoire of bamboo shark represents a subset of that of the catshark, with roughly similar thresholds. Amarogentin, one of the most bitter natural substances for humans, also elicited the highest signal amplitudes with both shark receptors. Other subsets of ligands are shared with basal bony fish T2Rs indicating an astonishing degree of functional conservation over nearly 500 mya of separate evolution. Both shark receptors respond to endogenous steroids as well as xenobiotic compounds, whereas separate receptors exist for xenobiotics both in early- and late-derived bony vertebrates (coelacanth, zebrafish, and human), consistent with the shark T2R reflecting the original ligand repertoire of the ancestral bitter taste receptor at the evolutionary origin of this family.

An Ancient MHC-Linked Gene Encodes a Nonrearranging Shark Antibody, UrIg, Convergent with IgG.

Gnathostome adaptive immunity is defined by the Ag receptors, Igs and TCRs, and the MHC. Cartilaginous fish are the oldest vertebrates with these adaptive hallmarks. We and others have unearthed nonrearranging Ag receptor-like genes in several vertebrates, some of which are encoded in the MHC or in MHC paralogous regions. One of these genes, named UrIg, was detected in the class III region of the shark MHC that encodes a protein with typical V and C domains such as those found in conventional Igs and TCRs. As no transmembrane region was detected in gene models or cDNAs, the protein does not appear to act as a receptor. Unlike some other shark Ig genes, the UrIg V region shows no evidence of RAG-mediated rearrangement, and thus it is likely related to other V genes that predated the invasion of the RAG transposon. The UrIg gene is present in all elasmobranchs and evolves conservatively, unlike Igs and TCRs. Also, unlike Ig/TCR, the gene is not expressed in secondary lymphoid tissues, but mainly in the liver. Recombinant forms of the molecule form disulfide-linked homodimers, which is the form also detected in many shark tissues by Western blotting. mAbs specific for UrIg identify the protein in the extracellular matrix of several shark tissues by immunohistochemistry. We propose that UrIg is related to the V gene invaded by the RAG transposon, consistent with the speculation of emergence of Ig/TCR within the MHC or proto-MHC.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae.

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Ancient and Nonuniform Loss of Olfactory Receptor Expression Renders the Shark Nose a De Facto Vomeronasal Organ.

Cartilaginous fishes are renowned for a keen sense of smell, a reputation based on behavioral observations and supported by the presence of large and morphologically complex olfactory organs. At the molecular level, genes belonging to the four families coding for most olfactory chemosensory receptors in other vertebrates have been identified in a chimera and a shark, but it was unknown whether they actually code for olfactory receptors in these species. Here, we describe the evolutionary dynamics of these gene families in cartilaginous fishes using genomes of a chimera, a skate, a sawfish, and eight sharks. The number of putative OR, TAAR, and V1R/ORA receptors is very low and stable, whereas the number of putative V2R/OlfC receptors is higher and much more dynamic. In the catshark Scyliorhinus canicula, we show that many V2R/OlfC receptors are expressed in the olfactory epithelium in the sparsely distributed pattern characteristic for olfactory receptors. In contrast, the other three vertebrate olfactory receptor families are either not expressed (OR) or only represented with a single receptor (V1R/ORA and TAAR). The complete overlap of markers of microvillous olfactory sensory neurons with pan-neuronal marker HuC in the olfactory organ suggests the same cell-type specificity of V2R/OlfC expression as for bony fishes, that is, in microvillous neurons. The relatively low number of olfactory receptors in cartilaginous fishes compared with bony fishes could be the result of an ancient and constant selection in favor of a high olfactory sensitivity at the expense of a high discrimination capability.

Antigen specific VNAR screening in whitespotted bamboo shark (Chiloscyllium plagiosum) with next generation sequencing.

IgNAR exhibits significant promise in the fields of cancer and anti-virus biotherapies. Notably, the variable regions of IgNAR (VNAR) possess comparable antigen binding affinity with much smaller molecular weight (∼12 kDa) compared to IgNAR. Antigen specific VNAR screening is a changeling work, which limits its application in medicine and therapy fields. Though phage display is a powerful tool for VNAR screening, it has a lot of drawbacks, such as small library coverage, low expression levels, unstable target protein, complicating and time-consuming procedures. Here we report VANR screening with next generation sequencing (NGS) could effectively overcome the limitations of phage display, and we successfully identified approximately 3000 BAFF-specific VNARs in Chiloscyllium plagiosum vaccinated with the BAFF antigen. The results of modelling and molecular dynamics simulation and ELISA assay demonstrated that one out of the top five abundant specific VNARs exhibited higher binding affinity to the BAFF antigen than those obtained through phage display screening. Our data indicates NGS would be an alternative way for VNAR screening with plenty of advantages.

Sharks Provide Evidence for a Highly Complex TNFSF Repertoire in the Jawed Vertebrate Ancestor.

Cytokines of the TNF superfamily (TNFSF) control many immunological processes and are implicated in the etiology of many immune disorders and diseases. Despite their obvious biological importance, the TNFSF repertoires of many species remain poorly characterized. In this study, we perform detailed bioinformatic, phylogenetic, and syntenic analyses of five cartilaginous fish genomes to identify their TNFSF repertoires. Strikingly, we find that shark genomes harbor ∼30 TNFSF genes, more than any other vertebrate examined to date and substantially more than humans. This is due to better retention of the ancestral jawed vertebrate TNFSF repertoire than any other jawed vertebrate lineage, combined with lineage-specific gene family expansions. All human TNFSFs appear in shark genomes, except for lymphotoxin-α (LTA; TNFSF1) and TNF (TNFSF2), and CD70 (TNFSF7) and 4-1BBL (TNFSF9), which diverged by tandem duplications early in tetrapod and mammalian evolution, respectively. Although lacking one-to-one LTA and TNF orthologs, sharks have evolved lineage-specific clusters of LTA/TNF co-orthologs. Other key findings include the presence of two BAFF (TNFSF13B) genes along with orthologs of APRIL (TNFSF13) and BALM (TNFSF13C) in sharks, and that all cartilaginous fish genomes harbor an ∼400-million-year-old cluster of multiple FASLG (TNFSF6) orthologs. Finally, sharks have retained seven ancestral jawed vertebrate TNFSF genes lost in humans. Taken together, our data indicate that the jawed vertebrate ancestor possessed a much larger and diverse TNFSF repertoire than previously hypothesized and oppose the idea that the cartilaginous fish immune system is "primitive" compared with that of mammals.

Detailed characterization of the complete mitochondrial genome of the oceanic whitetip shark Carcharhinus longimanus (Poey, 1861).

BACKGROUND: The oceanic whitetip shark Carcharhinus longimanus (family Carcharhinidae) is one of the largest sharks inhabiting all tropical and subtropical oceanic regions. Due to their life history traits and mortality attributed to pelagic longline fishing practices, this species is experiencing substantial population decline. Currently, C. longimanus is considered by the IUCN Red List of Threatened Species as "vulnerable" throughout its range and "critically endangered" in the western north Atlantic. This study sequences and describes the complete mitochondrial genome of C. longimanus in detail. METHODS AND RESULTS: The mitochondrial genome of C. longimanus was assembled through next-generation sequencing and then analyzed using specialized bioinformatics tools. The circular, double-stranded AT-rich mitogenome of C. longimanus is 16,704 bp long and contains 22 tRNA genes, 2 rRNA genes, 13 protein coding genes and a 1,065 bp long control region (CR). Out of the 22 tRNA genes, only one (tRNA-Ser1) lacked a typical 'cloverleaf' secondary structure. The prevalence of TTA (Leu), ATT (Ile) and CTA (Leu) codons in the PCGs likely contributes to the AT-rich nature of this mitogenome. In the CR, ten microsatellites were detected but no tandem repeats were found. Stem-and-loop secondary structures were common along the entire length of the CR. Ka/Ks values estimated for all PCGs were < 1, indicating that all the PCGs experience purifying selection. A phylomitogenomic analysis based on translated PCGs confirms the sister relationship between C. longimanus and C. obscurus. The analysis did not support the monophyly of the genus Carcharhinus. CONCLUSIONS: The assembled mitochondrial genome of this pelagic shark can provide insight into the phylogenetic relationships in the genus Carcharhinus and aid conservation and management efforts in the Central Pacific Ocean.

First molecular identification and phylogenetic illustration of Sarcocystis species infection in Red Sea shortfin mako shark (Isurus oxyrinchus Rafinesque, 1810).

BACKGROUND: members of the genus Sarcocystis are intracellular obligate protozoan parasites classified within the phylum Apicomplexa and have an obligate heteroxenous life cycle involving two hosts. A more comprehensive understanding of the prevalence and geographic range of different Sarcocystis species in marine ecosystems is needed globally and nationally. Hence, the objective of this study was to document the incidence of Sarcocystis infection in sharks within the aquarium ecosystem of Egypt and to identify the species through the characterization of the SSU rDNA gene. METHODS: All organs of the mako shark specimen underwent macroscopic screening to detect the existence of a Sarcocystis cyst. Ten cysts were collected from the intestine and processed separately to extract the genomic DNA. The polymerase chain reaction (PCR) was accomplished by amplifying a specific 18S ribosomal RNA (rRNA) gene fragment. Subsequently, the resulting amplicons were subjected to purification and sequencing processes. RESULTS: Macroscopic examination of the mako shark intestinal wall sample revealed the presence of Sarcocystis cysts of various sizes and shapes, and sequencing of the amplicons from Sarcocystis DNA revealed a 100% nucleotide identity with the sequence of Sarcocystis tenella recorded from sheep in Iran; The mako shark sequence has been deposited in the GeneBank with the accession number OQ721979. This study presents the first scientific evidence demonstrating the presence of the Sarcocystis parasite in sharks, thereby documenting this specific marine species as a novel intermediate host in the Sarcocystis life cycle. CONCLUSIONS: This is the first identification of Sarcocystis infection in sharks, and we anticipate it will be an essential study for future screenings and establishing effective management measures for this disease in aquatic ecosystems.

Parallel Evolution of Ameloblastic scpp Genes in Bony and Cartilaginous Vertebrates.

In bony vertebrates, skeletal mineralization relies on the secretory calcium-binding phosphoproteins (Scpp) family whose members are acidic extracellular proteins posttranslationally regulated by the Fam20°C kinase. As scpp genes are absent from the elephant shark genome, they are currently thought to be specific to bony fishes (osteichthyans). Here, we report a scpp gene present in elasmobranchs (sharks and rays) that evolved from local tandem duplication of sparc-L 5' exons and show that both genes experienced recent gene conversion in sharks. The elasmobranch scpp is remarkably similar to the osteichthyan scpp members as they share syntenic and gene structure features, code for a conserved signal peptide, tyrosine-rich and aspartate/glutamate-rich regions, and harbor putative Fam20°C phosphorylation sites. In addition, the catshark scpp is coexpressed with sparc-L and fam20°C in tooth and scale ameloblasts, similarly to some osteichthyan scpp genes. Despite these strong similarities, molecular clock and phylogenetic data demonstrate that the elasmobranch scpp gene originated independently from the osteichthyan scpp gene family. Our study reveals convergent events at the sparc-L locus in the two sister clades of jawed vertebrates, leading to parallel diversification of the skeletal biomineralization toolkit. The molecular evolution of sparc-L and its coexpression with fam20°C in catshark ameloblasts provides a unifying genetic basis that suggests that all convergent scpp duplicates inherited similar features from their sparc-L precursor. This conclusion supports a single origin for the hypermineralized outer odontode layer as produced by an ancestral developmental process performed by Sparc-L, implying the homology of the enamel and enameloid tissues in all vertebrates.

Philopatry influences the genetic population structure of the blacktip shark (Carcharhinus limbatus) at multiple spatial scales.

Understanding how interactions among microevolutionary forces generate genetic population structure of exploited species is vital to the implementation of management policies that facilitate persistence. Philopatry displayed by many coastal shark species can impact gene flow and facilitate selection, and has direct implications for the spatial scales of management. Here, genetic structure of the blacktip shark (Carcharhinus limbatus) was examined using a mixed-marker approach employing mitochondrial control region sequences and 4339 SNP-containing loci generated using ddRAD-Seq. Genetic variation was assessed among young-of-the-year sampled in 11 sites in waters of the United States in the western North Atlantic Ocean, including the Gulf of Mexico. Spatial and environmental analyses detected 68 nuclear loci putatively under selection, enabling separate assessments of neutral and adaptive genetic structure. Both mitochondrial and neutral SNP data indicated three genetically distinct units-the Atlantic, eastern Gulf, and western Gulf-that align with regional stocks and suggest regional philopatry by males and females. Heterogeneity at loci putatively under selection, associated with temperature and salinity, was observed among sites within Gulf units, suggesting local adaptation. Furthermore, five pairs of siblings were identified in the same site across timescales corresponding with female reproductive cycles. This indicates that females re-used a site for parturition, which has the potential to facilitate the sorting of adaptive variation among neighbouring sites. The results demonstrate differential impacts of microevolutionary forces at varying spatial scales and highlight the importance of conserving essential habitats to maintain sources of adaptive variation that may buffer species against environmental change.

Stepping up to genome scan allows stock differentiation in the worldwide distributed blue shark Prionace glauca.

The blue shark Prionace glauca is a top predator with one of the widest geographical distributions of any shark species. It is classified as Critically Endangered in the Mediterranean Sea, and Near Threatened globally. Previous genetic studies did not reject the null hypothesis of a single global population. The blue shark was proposed as a possible archetype of the "grey zone of population differentiation," coined to designate cases where population structure may be too recent or too faint to be detected using a limited set of markers. Here, blue shark samples collected throughout its global range were sequenced using a specific RAD method (DArTseq), which recovered 37,655 genome-wide single nucleotide polymorphisms (SNPs). Two main groups emerged, with Mediterranean Sea and northern Atlantic samples (Northern population) differentiated significantly from the Indo-west Pacific samples (Southern population). Significant pairwise FST values indicated further genetic differentiation within the Atlantic Ocean, and between the Atlantic Ocean and the Mediterranean Sea. Reconstruction of recent demographic history suggested divergence between Northern and Southern populations occurred about 500 generations ago and revealed a drastic reduction in effective population size from a large ancestral population. Our results illustrate the power of genome scans to detect population structure and reconstruct demographic history in highly migratory marine species. Given that the management plans of the blue shark (targeted or bycatch) fisheries currently assume panmictic regional stocks, we strongly recommend that the results presented here be considered in future stock assessments and conservation strategies.

Use of vivo-morpholinos for gene knockdown in the postnatal shark retina.

Work in the catshark Scyliorhinus canicula has shown that the evolutionary origin of postnatal neurogenesis in vertebrates is earlier than previously thought. Thus, the catshark can serve as a model of interest to understand postnatal neurogenic processes and their evolution in vertebrates. One of the best characterized neurogenic niches of the catshark CNS is found in the peripheral region of the retina. Unfortunately, the lack of genetic tools in sharks limits the possibilities to deepen in the study of genes involved in the neurogenic process. Here, we report a method for gene knockdown in the juvenile catshark retina based on the use of Vivo-Morpholinos. To establish the method, we designed Vivo-Morpholinos against the proliferation marker PCNA. We first evaluated the possible toxicity of 3 different intraocular administration regimes. After this optimization step, we show that a single intraocular injection of the PCNA Vivo-Morpholino decreases the expression of PCNA in the peripheral retina, which leads to reduced mitotic activity in this region. This method will help in deciphering the role of other genes potentially involved in postnatal neurogenesis in this animal model.

Evolution of reproductive modes in sharks and rays.

The ecological and life history drivers of the diversification of reproductive modes in early vertebrates are not fully understood. Sharks, rays and chimaeras (group Chondrichthyes) have an unusually diverse variety of reproductive modes and are thus an ideal group to test the factors driving the evolution of reproductive complexity. Here, using 960 species representing all major Chondrichthyes taxa, we reconstruct the evolution of their reproduction modes and investigate the ecological and life history predictors of reproduction. We show that the ancestral Chondrichthyes state was egg-laying and find multiple independent transitions between egg-laying and live-bearing via an intermediate state of yolk-only live-bearing. Using phylogenetically informed analysis, we also show that live-bearing species have larger body size and larger offspring than egg-laying species. In addition, live-bearing species are distributed over shallow to intermediate depths, while egg-layers are typically found in deeper waters. This suggests that live-bearing is more closely associated with pelagic, rather than demersal habitats. Taken together, using a basal vertebrate group as a model, we demonstrat how reproductive mode co-evolves with environmental conditions and life-history traits.

IgNAR characterization and gene loci identification in whitespotted bamboo shark (Chiloscyllium plagiosum) genome.

Single domain antibodies (sdAb) are promising candidates in cancer and anti-virus biotherapies for their unique structure characters. Though VHH and IgNAR have been discovered in camelidae and nurse shark (Ginlymostoma cirratum) respectively serval decades ago, expense of these large animals still limits the studies and applications of sdAb. Recently, IgNAR has been found in whitespotted bamboo shark (Chiloscyllium plagiosum), a small-sized sharks, while how to characterize and achieved the IgNAR of whitespotted bamboo shark is still unclear. In our research, we identified four IgNAR coding gene loci in whitespotted bamboo shark chromosome 44 (NC_057753.1), and primers were designed for single domain variable regions of IgNAR (VNAR) libraries preparation. Following sequencing results revealed that all plasmids constructed with our predicted VNAR libraries contained VNAR coding sequences, which confirmed the specificities of our primers in VNAR amplification. To our surprise, ≥90% VNAR sequences were encoded by IgNAR1, which suggests IgNAR1 is the most active IgNAR transcription locus in whitespotted bamboo shark. Interestingly, we found IgNAR(ΔC2-C3) were encoded by IgNAR3. Our findings gave a new sight of whitespotted bamboo shark IgNAR, which would broad the way of IgNAR studies and applications in biotherapies.

Shark VNAR phage display libraries: An alternative source for therapeutic and diagnostic recombinant antibody fragments.

The development of recombinant antibody fragments as promising alternatives to full-length immunoglobulins offers vast opportunities for biomedicine. Antibody fragments have important advantages compared with conventional monoclonal antibodies that make them attractive for the biotech industry: superior stability and solubility, reduced immunogenicity, higher specificity and affinity, capacity to target the hidden epitope and cross the blood-brain barrier, the ability to refold after heat denaturation and inexpensive and easy large-scale production. Different antibody formats such as antigen-binding fragments (Fab), single-chain fragment variable (scFv) consisting of the antigen-binding domains of Ig heavy (VH) and light (VL) chain regions connected by a flexible peptide linker, single-domain antibody fragments (sdAbs) like camelid heavy-chain variable domains (VHHs) and shark variable new antigen receptor (VNARs), and bispecific antibodies (bsAbs) are currently being evaluated as diagnostics or therapeutics in preclinical studies and clinical trials. In the present review, we summarize and discuss studies on VNARs, the smallest recombinant antibody fragment, obtained after the screening of different types of phage display antibody libraries. Results published until March 2023 are discussed.

A Highly Complex, MHC-Linked, 350 Million-Year-Old Shark Nonclassical Class I Lineage.

Cartilaginous fish, or Chondrichthyes, are the oldest extant vertebrates to possess the MHC and the Ig superfamily-based Ag receptors, the defining genes of the gnathostome adaptive immune system. In this work, we have identified a novel MHC lineage, UEA, a complex multigene nonclassical class I family found in sharks (division Selachii) but not detected in chimaeras (subclass Holocephali) or rays (division Batoidea). This new lineage is distantly related to the previously reported nonclassical class I lineage UCA, which appears to be present only in dogfish sharks (order Squaliformes). UEA lacks conservation of the nine invariant residues in the peptide (ligand)-binding regions (PBR) that bind to the N and C termini of bound peptide in most vertebrate classical class I proteins, which are replaced by relatively hydrophobic residues compared with the classical UAA. In fact, UEA and UCA proteins have the most hydrophobic-predicted PBR of all identified chondrichthyan class I molecules. UEA genes detected in the whale shark and bamboo shark genome projects are MHC linked. Consistent with UEA comprising a very large gene family, we detected weak expression in different tissues of the nurse shark via Northern blotting and RNA sequencing. UEA genes fall into three sublineages with unique characteristics in the PBR. UEA shares structural and genetic features with certain nonclassical class I genes in other vertebrates, such as the highly complex XNC nonclassical class I genes in Xenopus, and we anticipate that each shark gene, or at least each sublineage, will have a unique function, perhaps in bacterial defense.

Characterization of 35 new microsatellite markers for the blacktip reef shark (Carcharhinus melanopterus) and cross-species amplification in eight other shark species.

BACKGROUND: Shark species are overfished at a global scale, as they are poached for the finning industry or are caught as bycatch. Efficient conservation measures require fine-scale spatial and temporal studies to characterize shark habitat use, infer migratory habits, analyze relatedness, and detect population genetic differentiation. Gathering these types of data is costly and time-consuming, especially when it requires collection of shark tissue samples. METHODS AND RESULTS: Genetic tools, such as microsatellite markers, are the most economical sampling method for collecting genetic data, as they enable the estimation of genetic diversity, population structure and parentage relationships and are thus an efficient way to inform conservation strategies. Here, a set of 45 microsatellite loci was tested on three blacktip reef shark (Carcharhinus melanopterus) populations from three Polynesian islands: Moorea, Morane and Tenararo. The set was composed of 10 previously published microsatellite markers and 35 microsatellite markers that were developed specifically for C. melanopterus as part of the present study. The 35 novel and 10 existing loci were cross-amplified on eight additional shark species (Carcharhinus amblyrhynchos, C. longimanus, C. sorrah, Galeocerdo cuvier, Negaprion acutidens, Prionacea glauca, Rhincodon typus and Sphyrna lewini). These species had an average of 69% of successful amplification, considered if at least 50% of the individual samples being successfully amplified per species and per locus. CONCLUSIONS: This novel microsatellite marker set will help address numerous knowledge gaps that remain, concerning genetic stock identification, shark behavior and reproduction via parentage analysis.