Resurrection of genus Phocanema Myers, 1959, as a genus independent from Pseudoterranova Mozgovoĭ, 1953, for nematode species (Anisakidae) parasitic in pinnipeds and cetaceans, respectively

Species of the genus Pseudoterranova, infect kogiid cetaceans and pinnipeds. However, there is mounting molecular evidence that those from cetaceans and pinnipeds are not congeneric. Here, we provide further evidence of the non-monophyly of members of Pseudoterranova from phylogenetic analyses of the conserved nuclear LSU rDNA gene, entire ITS rDNA region and mtDNA cox2 gene, and identify morphological characters that may be used to distinguish the members of the two clades. We propose the resurrection of the genus Phocanema, with Ph. decipiens (sensu stricto) as the type species, to encompass Ph. decipiens, Ph. azarasi, Ph. bulbosa, Ph. cattani and Ph. krabbei, all parasites of pinnipeds. We propose to restrict the conception of genus Pseudoterranova, which now harbours two species infecting kogiid whales; Ps. kogiae (type species) and Ps. ceticola. Members of the genera Phocanema and Pseudoterranova differ by the shape and orientation of the lips, relative tail lengths, adult size, type of final host (pinniped vs. cetacean) and phylogenetic placement based on nuclear rDNA and mtDNA cox2 sequences.

Genetic diversity and demographic history of the leopard seal: A Southern Ocean top predator

Leopard seals (Hydrurga leptonyx) are top predators that can exert substantial top-down control of their Antarctic prey species. However, population trends and genetic diversity of leopard seals remain understudied, limiting our understanding of their ecological role. We investigated the genetic diversity, effective population size and demographic history of leopard seals to provide fundamental data that contextualizes their predatory influence on Antarctic ecosystems. Ninety leopard seals were sampled from the northern Antarctic Peninsula during the austral summers of 2008-2019 and a 405bp segment of the mitochondrial control region was sequenced for each individual. We uncovered moderate levels of nucleotide (π = 0.013) and haplotype (Hd = 0.96) diversity, and the effective population size was estimated at around 24,000 individuals (NE = 24,376; 95% CI: 16,876-33,126). Consistent with findings from other ice-breeding pinnipeds, Bayesian skyline analysis also revealed evidence for population expansion during the last glacial maximum, suggesting that historical population growth may have been boosted by an increase in the abundance of sea ice. Although leopard seals can be found in warmer, sub-Antarctic locations, the species' core habitat is centered on the Antarctic, making it inherently vulnerable to the loss of sea ice habitat due to climate change. Therefore, detailed assessments of past and present leopard seal population trends are needed to inform policies for Antarctic ecosystems.

Spectral tuning and deactivation kinetics of marine mammal melanopsins

In mammals, the photopigment melanopsin (Opn4) is found in a subset of retinal ganglion cells that serve light detection for circadian photoentrainment and pupil constriction (i.e., mydriasis). For a given species, the efficiency of photoentrainment and length of time that mydriasis occurs is determined by the spectral sensitivity and deactivation kinetics of melanopsin, respectively, and to date, neither of these properties have been described in marine mammals. Previous work has indicated that the absorbance maxima (λmax) of marine mammal rhodopsins (Rh1) have diversified to match the available light spectra at foraging depths. However, similar to the melanopsin λmax of terrestrial mammals (~480 nm), the melanopsins of marine mammals may be conserved, with λmax values tuned to the spectrum of solar irradiance at the water's surface. Here, we investigated the Opn4 pigments of 17 marine mammal species inhabiting diverse photic environments including the Infraorder Cetacea, as well as the Orders Sirenia and Carnivora. Both genomic and cDNA sequences were used to deduce amino acid sequences to identify substitutions most likely involved in spectral tuning and deactivation kinetics of the Opn4 pigments. Our results show that there appears to be no amino acid substitutions in marine mammal Opn4 opsins that would result in any significant change in λmax values relative to their terrestrial counterparts. We also found some marine mammal species to lack several phosphorylation sites in the carboxyl terminal domain of their Opn4 pigments that result in significantly slower deactivation kinetics, and thus longer mydriasis, compared to terrestrial controls. This finding was restricted to cetacean species previously found to lack cone photoreceptor opsins, a condition known as rod monochromacy. These results suggest that the rod monochromat whales rely on extended pupillary constriction to prevent photobleaching of the highly photosensitive all-rod retina when moving between photopic and scotopic conditions.

Predicting sample success for large-scale ancient DNA studies on marine mammals

In recent years, nonhuman ancient DNA studies have begun to focus on larger sample sizes and whole genomes, offering the potential to reveal exciting and hitherto unknown answers to ongoing biological and archaeological questions. However, one major limitation to such studies is the substantial financial and time investments still required during sample screening, due to uncertainty regarding successful sample selection. This study investigates the effect of a wide range of sample properties including latitude, sample age, skeletal element, collagen preservation, and context on endogenous content and DNA damage profiles for 317 ancient and historic pinniped samples collected from across the North Atlantic and surrounding regions. Using generalised linear and mixed-effect models, we found that a range of factors affected DNA preservation within each of the species under consideration. The most important findings were that endogenous content varied significantly within species according to context, the type of skeletal element, the collagen content and collection year. There also appears to be an effect of the sample's geographic origin, with samples from the Arctic generally showing higher endogenous content and lower damage rates. Both latitude and sample age were found to have significant relationships with damage levels, but only for walrus samples. Sex, ontogenetic age and extraction material preparation were not found to have any significant relationship with DNA preservation. Overall, skeletal element and sample context were found to be the most influential factors and should therefore be considered when selecting samples for large-scale ancient genome studies.

Low guanylyl cyclase activity in Weddell seals: implications for peripheral vasoconstriction and perfusion of the brain during diving

Nitric oxide (NO) is a potent vasodilator, which improves perfusion and oxygen delivery during tissue hypoxia in terrestrial animals. The vertebrate dive response involves vasoconstriction in select tissues, which persists despite profound hypoxia. Using tissues collected from Weddell seals at necropsy, we investigated whether vasoconstriction is aided by downregulation of local hypoxia signaling mechanisms. We focused on NO-soluble guanylyl cyclase (GC)-cGMP signaling, a well-known vasodilatory transduction pathway. Seals have a lower GC protein abundance, activity, and capacity to respond to NO stimulation than do terrestrial mammals. In seal lung homogenates, GC produced less cGMP (20.1 ± 3.7 pmol·mg protein-1·min-1) than the lungs of dogs (-80 ± 144 pmol·mg protein-1·min-1 less than seals), sheep (-472 ± 96), rats (-664 ± 104) or mice (-1,160 ± 104, P < 0.0001). Amino acid sequences of the GC enzyme α-subunits differed between seals and terrestrial mammals, potentially affecting their structure and function. Vasoconstriction in diving Weddell seals is not consistent across tissues; perfusion is maintained in the brain and heart but decreased in other organs such as the kidney. A NO donor increased median GC activity 49.5-fold in the seal brain but only 27.4-fold in the kidney, consistent with the priority of cerebral perfusion during diving. Nos3 expression was high in the seal brain, which could improve NO production and vasodilatory potential. Conversely, Pde5a expression was high in the seal renal artery, which may increase cGMP breakdown and vasoconstriction in the kidney. Taken together, the results of this study suggest that alterations in the NO-cGMP pathway facilitate the diving response.

Automated masking of AFLP markers improves reliability of phylogenetic analyses

The amplified fragment length polymorphisms (AFLP) method has become an attractive tool in phylogenetics due to the ease with which large numbers of characters can be generated. In contrast to sequence-based phylogenetic approaches, AFLP data consist of anonymous multilocus markers. However, potential artificial amplifications or amplification failures of fragments contained in the AFLP data set will reduce AFLP reliability especially in phylogenetic inferences. In the present study, we introduce a new automated scoring approach, called “AMARE” (AFLP MAtrix REduction). The approach is based on replicates and makes marker selection dependent on marker reproducibility to control for scoring errors. To demonstrate the effectiveness of our approach we record error rate estimations, resolution scores, PCoA and stemminess calculations. As in general the true tree (i.e. the species phylogeny) is not known, we tested AMARE with empirical, already published AFLP data sets, and compared tree topologies of different AMARE generated character matrices to existing phylogenetic trees and/or other independent sources such as morphological and geographical data. It turns out that the selection of masked character matrices with highest resolution scores gave similar or even better phylogenetic results than the original AFLP data sets.

Transactivation potencies of Baikal seal constitutive active/androstane receptor by persistent organic pollutants and brominated flame retardants

To characterize ligand-dependent transcriptional activation of constitutive active/androstane receptor (CAR) in aquatic mammals, transactivation potentials of the Baikal seal (Pusa sibirica) CAR (bsCAR) by environmental pollutants, including persistent organic pollutants (POPs) and brominated flame retardants (BFRs), were investigated using an in vitro reporter gene assay, and compared with those of the mouse CAR (mCAR). Measurement of luciferase reporter gene activities demonstrated that the seal CAR was activated by POPs, including a technical mixture of PCBs (Kanechlor-500), certain individual PCB congeners, DDT compounds, and trans-nonachlor. No or slight bsCAR-dependent activity was detected in experiments with PBDE congeners and HBCDs. The interspecies comparison of lowest observed effect concentration (LOEC) for CAR transactivation byeach compound revealed that bsCAR responds more sensitively to PCBs than mCAR. In addition, bsCAR was weakly deactivated by PBDE99, whereas mCAR transcriptional activity decreased weakly by PBDE100, PBDE154, and PBDE187. Comparison of reporter gene activities by the congeners with the same IUPAC numbers among PCBs and PBDEs revealed that both bsCAR and mCAR were not activated by PBDE99 and PBDE153, but were activated by PCB99 and PCB153. The small ligand-binding pocket in CAR may contribute to difference in response between PCBs and PBDEs. Given that ethical rationale prevents dosing studies with such organohalogens in aquatic mammals, our in vitro assay system constructed with CAR cDNA from a species of interest provides a useful and realistic alternative approach in ecotoxicology.

Evolution of Pinnipedia lactation strategies: a potential role for alpha-lactalbumin?

Despite the considerable variation in milk composition found among mammals, a constituent common across all groups is lactose, the main sugar and osmole in most eutherians milk. Exceptions to this are the families Otariidae (fur seals and sea lions) and Odobenidae (walruses), where lactose has not been detected. We investigated the molecular basis for this by cloning alpha-lactalbumin, the modifier protein of the lactose synthase complex. A mutation was observed which, in addition to preventing lactose production, may enable otariids to maintain lactation despite the extremely long inter-suckling intervals during the mother's time at sea foraging (more than 23 days in some species).

Sequence variation in the Tbx4 gene in marine mammals

The amino-acid sequences of the T-domain region of the Tbx4 gene, which is required for hindlimb development, are 100% identical in humans and mice. Cetaceans have lost most of their hindlimb structure, although hindlimb buds are present in very early cetacean embryos. To examine whether the Tbx4 gene has the same function in cetaceans as in other mammals, we analyzed Tbx4 sequences from cetaceans, dugong, artiodactyls and marine carnivores. A total of 39 primers were designed using human and dog Tbx4 nucleotide sequences. Exons 3, 4, 5, 6, 7, and 8 of the Tbx4 genes from cetaceans, artiodactyls, and marine carnivores were sequenced. Non-synonymous substitution sites were detected in the T-domain regions from some cetacean species, but were not detected in those from artiodactyls, the dugong, or the carnivores. The C-terminal regions contained a number of non-synonymous substitutions. Although some indels were present, they were in groups of three nucleotides and therefore did not cause frame shifts. The dN/dS values for the T-domain and C-terminal regions of the cetacean and artiodactylous Tbx4 genes were much lower than 1, indicating that the Tbx4 gene maintains it function in cetaceans, although full expression leading to hindlimb development is suppressed.

Positive selection in the N-terminal extramembrane domain of lung surfactant protein C (SP-C) in marine mammals

Maximum-likelihood models of codon and amino acid substitution were used to analyze the lung-specific surfactant protein C (SP-C) from terrestrial, semi-aquatic, and diving mammals to identify lineages and amino acid sites under positive selection. Site models used the nonsynonymous/synonymous rate ratio (omega) as an indicator of selection pressure. Mechanistic models used physicochemical distances between amino acid substitutions to specify nonsynonymous substitution rates. Site models strongly identified positive selection at different sites in the polar N-terminal extramembrane domain of SP-C in the three diving lineages: site 2 in the cetaceans (whales and dolphins), sites 7, 9, and 10 in the pinnipeds (seals and sea lions), and sites 2, 9, and 10 in the sirenians (dugongs and manatees). The only semi-aquatic contrast to indicate positive selection at site 10 was that including the polar bear, which had the largest body mass of the semi-aquatic species. Analysis of the biophysical properties that were influential in determining the amino acid substitutions showed that isoelectric point, chemical composition of the side chain, polarity, and hydrophobicity were the crucial determinants. Amino acid substitutions at these sites may lead to stronger binding of the N-terminal domain to the surfactant phospholipid film and to increased adsorption of the protein to the air-liquid interface. Both properties are advantageous for the repeated collapse and reinflation of the lung upon diving and resurfacing and may reflect adaptations to the high hydrostatic pressures experienced during diving.

Sexing pinnipeds with ZFX and ZFY loci

We developed and tested a protocol for determining the sex of individual pinnipeds using the sex-chromosome-specific genes ZFX and ZFY. We screened a total of 368 seals (168 crabeater, Lobodon carcinophaga; 159 Weddell, Leptonychotes weddellii; and 41 Ross, Ommatophoca rossii) of known or unknown sex and compared the molecular sex to the sex assigned at the time of biopsy sample collection in the Ross and Amundsen seas, Antarctica. We also screened 6 captive northern elephant seals (Mirounga angustirostris) and 2 captive California sea lions (Zalophus californianus) of known sex. The assigned sex and genetic sex agreed for virtually all seals. Indeed, discrepancies ranged from 0.0% to 6.7% among species. It is not clear, however, if the few mis-assignments of sex occurred in situ or in the laboratory. The assigned morphological and molecular sex might both be correct with the discrepancies owing perhaps to developmental effects of environmental pollution. A subset of individuals sequenced at both loci revealed no intraspecific sequence variation. There was, however, sequence variation among species at both loci, which allowed them to be uniquely identified with as few as 2 and as many as 31 nucleotides.

Pinniped phylogeny and a new hypothesis for their origin and dispersal

The relationships and the zoogeography of the three extant pinniped families, Otariidae (sea lions and fur seals), Odobenidae (one extant species, the walrus), and Phocidae (true seals), have been contentious. Here, we address these topics in a molecular study that includes all extant species of true seals and sea lions, four fur seals and the walrus. Contrary to prevailing morphological views the analyses conclusively showed monophyletic Pinnipedia with a basal split between Otarioidea (Otariidae+Odobenidae) and Phocidae. The northern fur seal was the sister to all remaining otariids and neither sea lions nor arctocephaline fur seals were recognized as monophyletic entities. The basal Phocidae split between Monachinae (monk seals and southern true seals) and Phocinae (northern true seals) was strongly supported. The phylogeny of the Phocinae suggests that the ancestors of Cystophora (hooded seal) and the Phocini (e.g. harp seal, ringed seal) adapted to Arctic conditions and ice-breeding before 12 MYA (million years ago) as supported by the white natal coat of these lineages. The origin of the endemic Caspian and Baikal seals was dated well before the onset of major Pleistocene glaciations. The current findings, together with recent advances in pinniped paleontology, allow the proposal of a new hypothesis for pinniped origin and early dispersal. The hypothesis posits that pinnipeds originated on the North American continent with early otarioid and otariid divergences taking place in the northeast Pacific and those of the phocids in coastal areas of southeast N America for later dispersal to colder environments in the N Atlantic and the Arctic Basin, and in Antarctic waters.

Evidence from nuclear DNA sequences sheds light on the phylogenetic relationships of Pinnipedia: single origin with affinity to Musteloidea

Considerable long-standing controversy and confusion surround the phylogenetic affinities of pinnipeds, the largely marine group of "fin-footed" members of the placental mammalian order Carnivora. Until most recently, the two major competing hypotheses were that the pinnipeds have a single (monophyletic) origin from a bear-like ancestor, or that they have a dual (diphyletic) origin, with sea lions (Otariidae) derived from a bear-like ancestor, and seals (Phocidae) derived from an otter-, mustelid-, or musteloid-like ancestor. We examined phylogenetic relationships among 29 species of arctoid carnivorans using a concatenated sequence of 3228 bp from three nuclear loci (apolipoprotein B, APOB; interphotoreceptor retinoid-binding protein, IRBP; recombination-activating gene 1, RAG1). The species represented Pinnipedia (Otariidae: Callorhinus, Eumetopias; Phocidae: Phoca), bears (Ursidae: Ursus, Melursus), and Musteloidea (Mustelidae: Mustela, Enhydra, Melogale, Martes, Gulo, Meles; Procyonidae: Procyon; Ailuridae: Ailurus; Mephitidae: Mephitis). Maximum parsimony, maximum likelihood, and Bayesian inference phylogenetic analyses of separate and combined datasets produced trees with largely congruent topologies. The analyses of the combined dataset resulted in well-resolved and well-supported phylogeny reconstructions. Evidence from nuclear DNA evolution presented here contradicts the two major hypotheses of pinniped relationships and strongly suggests a single origin of the pinnipeds from an arctoid ancestor shared with Musteloidea to the exclusion of Ursidae.

Molecular genetic diversity and evolution at the MHC DQB locus in four species of pinnipeds

Variation was investigated at exon 2 (including part of the putative peptide-binding region) of the class II major histocompatibility complex (MHC) DQB locus for two congeneric phocid seal species and two congeneric otariid seal species. Polymorphism in one phocid species, the southern elephant seal (Mirounga leonina), was comparable to that seen in human populations, while the other phocid, the northern elephant seal (Mirounga angustirostris), has been through a severe population bottleneck and exhibited much less variation at this locus. A phylogenetic comparison of the four species was consistent with the trans-specific pattern of evolution described for other taxa at this locus, and relative nonsynonymous and synonymous substitution rates suggest the maintenance of polymorphisms by natural selection. A comparison of sequence patterns also suggested that some variation could have been generated through recombinational events, primarily within genera. These results suggest a pattern of evolution of the immune response in pinnipeds similar to that in terrestrial mammal species.

Interspecific microsatellite markers for the study of pinniped populations

Microsatellites have rapidly become the marker of choice for a wide variety of population genetic studies. Here we describe 20 pinniped microsatellite markers which have been tested across 18 pinniped species. The majority of these markers have broad utility in all pinnipeds and provide a strong base for detailed population genetic studies in the Pinnipedia.

Use of spectral analysis to test hypotheses on the origin of pinnipeds

The evolutionary origin of the pinnipeds (seals, sea lions, and walruses) is still uncertain. Most authors support a hypothesis of a monophyletic origin of the pinnipeds from a caniform carnivore. A minority view suggests a diphyletic origin with true seals being related to the mustelids (otters and ferrets). The phylogenetic relationships of the walrus to other pinniped and carnivore families are also still particularly problematic. Here we examined the relative support for mono- and diphyletic hypotheses using DNA sequence data from the mitochondrial small subunit (12S) rRNA and cytochrome b genes. We first analyzed a small group of taxa representing the three pinniped families (Phocidae, Otariidae, and Odobenidae) and caniform carnivore families thought to be related to them. We inferred phylogenetic reconstructions from DNA sequence data using standard parsimony and neighbor-joining algorithms for phylogenetic inference as well as a new method called spectral analysis (Hendy and Penny) in which phylogenetic information is displayed independently of any selected tree. We identified and compensated for potential sources of error known to lead to selection of incorrect phylogenetic trees. These include sampling error, unequal evolutionary rates on lineages, unequal nucleotide composition among lineages, unequal rates of change at different sites, and inappropriate tree selection criteria. To correct for these errors, we performed additional transformations of the observed substitution patterns in the sequence data, applied more stringent structural constraints to the analyses, and included several additional taxa to help resolve long, unbranched lineages in the tree. We find that there is strong support for a monophyletic origin of the pinnipeds from within the caniform carnivores, close to the bear/raccoon/panda radiation. Evidence for a diphyletic origin was very weak and can be partially attributed to unequal nucleotide compositions among the taxa analyzed. Subsequently, there is slightly more evidence for grouping the walrus with the eared seals versus the true seals. A more conservative interpretation, however, is that the walrus is an early, but not the first, independent divergence from the common pinniped ancestor.

Multiple nuclear-gene phylogenies: application to pinnipeds and comparison with a mitochondrial DNA gene phylogeny

Phylogenetic analyses of closely related species should use information from multiple, independent genes with relatively high rates of sequence evolution. To investigate species for which there are few prior sequence data for single-copy nuclear (scnDNA) genes, primers for gene amplification can be designed to highly conserved regions of exons in order to amplify both coding (exons) and noncoding (introns) sequences. We have explored this approach in a phylogenetic analysis of six species of pinnipeds that, together with terrestrial carnivore outgroups, encompass divergence times < or = 40-50 Mya. We sequenced one intron from each of the aldolase A (ALD-A), aldolase C (ALD-C), and histone H2AF genes; one exon from the major-histocompatibility-complex DQA gene; a H2AF processed pseudogene (psi H2AF); and, for comparison with the nuclear genes, the 5' portion of the mitochondrial DNA (mtDNA) control region. The pinniped psi H2AF genes were found to be of limited use because they were paralogous with the gene in the outgroup. The rate of silent substitution in scnDNA (primarily introns) was 5-10-fold lower than that for mtDNA control region I, and scnDNA sequence divergence increased linearly with time < or = 40-50 Mya. Alleles at three polymorphic scnDNA loci (ALD-A, H2AF, and DQA) in the southern elephant seal were paraphyletic with respect to the allele from the closely related northern elephant seal, while the more numerous mtDNA alleles were monophyletic. This we attribute to the consequences of a higher mutation rate rather than to a lower effective population size of mtDNA compared with scnDNA. Within the short (i.e., < 500-bp) sequences of individual scnDNA sequences, phylogenetically informative variation was insufficient to obtain robust phylogenies. However, the combined scnDNA sequences produced a well-supported phylogeny congruent with that derived from mtDNA. This analysis illustrates the high resolution of mtDNA sequences compared with a similar length of scnDNA sequence, but it also demonstrates the utility of combining information from multiple short scnDNA sequences obtained using broadly applicable primers.

The seal myoglobin gene: an unusually long globin gene

The tetrameric haemoglobins of vertebrates are encoded by alpha- and beta-globin gene families which arose during evolution by a succession of gene duplications, commencing with an alpha beta globin gene duplication which occurred about 500 Myr ago, early in the evolution of the vertebrates. All functional alpha- and beta- globin genes analysed so far share common features, including coding sequence homologies and the presence of two introns within the coding sequence at locations which correlate with interdomain boundaries within globin polypeptides. Here we describe the isolation and characterization of an additional diverged member of the globin gene family, the seal myoglobin gene. We show that monomeric myoglobin, which diverged from haemoglobin about 600-800 Myr ago before the appearance of tetrameric haemoglobins, is also specified by a gene containing two introns at positions precisely homologous to haemoglobin introns. Unlike vertebrate haemoglobin genes, however, the noncoding regions of the seal myoglobin gene are remarkably long.

Molecular cloning of seal myoglobin mRNA

Grey seal skeletal muscle containing high levels of myoglobin was used to prepare poly(A)+ RNA. In vitro translation of this RNA produced a range of polypeptides including myoglobin. cDNA was prepared by reverse transcription of muscle poly(A)+ RNA and cloned into the plasmid pAT 153. 4% of cDNA recombinants were shown to contain myoglobin cDNA inserts. DNA sequence analysis of one clone (pSM 178) which contained a relatively large myoglobin cDNA insert showed an incomplete cDNA comprising the terminal 293 nucleotides of 3' non-translated mRNA sequences. Hybridization experiments using this myoglobin cDNA indicated that seal myoglobin is coded by a single gene which is transcribed to give a 1400 nucleotide mRNA considerably longer than related haemoglobin mRNAs.

Karyotype stability in marine mammals

A considerable amount of information is now available on the karyology of the marine mammals, Pinnipedia and Cetacea. This information emanates both from banded and unbanded karyotypes. A majority of the pinnipeds has now been investigated karyologically, and representatives of all extant cetacean families have also been studied. The karyotype stability among the marine mammals has been stressed earlier, and the present discussion emphasizes this point further.

The chromosomes of the southern elephant seal, Mirounga leonina (Phocidae: Mammalia)

Chromosome analyses were carried out on bone marrow cells from three elephant seal pups (Mirounga leonina), a species characterized by a 2n = 34 chromosome number. Arm ratios for the different chromosomes were calculated to facilitate the classification of chromosomes into groups. The karyotype of the southern elephant seal differs from other reported 34-chromosome phocids in the morphology of the Y chromosome and in that of one of the autosomal chromosome pairs.