Abstract
Jetsam ambergris, found on beaches worldwide, has always been assumed to originate as a natural product of sperm whales (Physeteroidea). However, only indirect evidence has ever been produced for this, such as the presence of whale prey remains in ambergris. Here, we extracted and analysed DNA sequences from jetsam ambergris from beaches in New Zealand and Sri Lanka, and sequences from ambergris of a sperm whale beached in The Netherlands. The lipid-rich composition of ambergris facilitated high preservation-quality of endogenous DNA, upon which we performed shotgun Illumina sequencing. Alignment of mitochondrial and nuclear genome sequences with open-access reference data for multiple whale species confirms that all three jetsam samples derived originally from sperm whales (Physeter macrocephalus). Shotgun sequencing here also provides implications for metagenomic insights into ambergris-preserved DNA. These results demonstrate significant implications for elucidating the origins of jetsam ambergris as a prized natural product, and also for the understanding of sperm whale metabolism and diet, and the ecological mechanisms underlying these coproliths.
1. Introduction
Ambergris, a known natural product of the sperm whale [1–4], is also found as jetsam on beaches worldwide [5], and has been highly prized for its utility in the perfume industry [6]. Although it has long been held that the jetsam ambergris collected on beaches originates from sperm whales [4], little or no evidence for this has ever been published, and distinctions exist between such samples and samples of ambergris directly taken from sperm whales. For example, jetsam ambergris samples generally contain much higher proportions of the triterpenoid alcohol ambrein and much lower proportions of sterols than do samples of ambergris from sperm whales [5,7]. Conversely though, jetsam ambergris sometimes does contain fragments of squid beaks [4], and since cephalopods, such as squid, constitute the major dietary component of sperm whales, this has been cited as evidence of an origin of the jetsam coproliths from sperm whales. It is even theorized ambergris may originate as a pathological secretion from the irritant of the hard squid beak chitin [8]. However, other marine mammal species (e.g. members of Globicephala and Ziphiidae) also predate on squid [9–11], and some (including dwarf and pygmy sperm whales) are also cited as potential sources of ambergris [4]. Therefore, to further elucidate the origin of jetsam ambergris, we analysed DNA from an ambergris sample collected from a sperm whale beached in The Netherlands and compared it with DNA sequences isolated from jetsam ambergris collected from beaches in New Zealand and Sri Lanka.
Ambergris is held to be predominantly composed of ambrein due to its production from squalene, a common metabolic product in many organisms [12]. This process may be induced by gut microbial influence, and precipitates in dense, solid masses within the whale colon [7]. The coprolitic accretions that result are compositionally well-suited to preserving DNA from the colon since ambrein is hydrophobic and apparently resistant to degradation within the acidic enteric environment. Evidence from radiocarbon dating certainly indicates resistance to microbial and photo-degradation in the marine environment for up to a millennium in some jetsam ambergris samples [13]. We hypothesized that such material might provide an opportune cache for preserving DNA, even after prolonged exposure to detrimental conditions at sea.
2. Material and methods
Jetsam ambergris specimens from the North Sea, the Indian Ocean and the Pacific were analysed, representing the material's global distribution [13]. Three jetsam ambergris specimens (one from Sri Lanka, two from Pitt Island, New Zealand) were subsampled for DNA extraction. A fourth specimen originated from dissection of a male sperm whale beached in December 2012, at Razende Bol near Texel, The Netherlands. The latter ‘fresh’ ambergris, from a confirmed sperm whale carcass, provided a known comparison to the jetsam specimens with unconfirmed biological history. Specimens of ambergris were obtained and analysed for ambrein and faecal sterol content in previous studies [5,7].
DNA extraction and sequencing were undertaken at the GLOBE Institute, University of Copenhagen, in a dedicated ancient DNA laboratory following strict procedures for minimization of contamination. Approximately 120 mg was subsampled (figure 1 and table 1) for DNA extraction. Samples were incubated in 400 µl proteinase K-containing buffer following Gilbert et al. [14] at 56°C for 10 h; supernatants were then treated using a phenol–chloroform step following Carøe et al. [15] and purified using Monarch DNA Cleanup Columns (5 µg) (New England Biolabs, Beverly, MA, USA) according to the manufacturer's guidelines. Double-stranded libraries were built from DNA extracts following the BEST protocol [15], designed and proven specifically for sequencing of ancient and degraded DNA. Libraries were amplified and indexed through PCR using PfuTurbo Cx Hotstart (Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer's guidelines. Products were pooled at equimolar concentration before sequencing on an Illumina HiSeq 4000 (Illumina, San Diego, CA, USA) platform, using 80 bp single end read chemistry at the Danish National High-throughput DNA Sequencing Centre, Copenhagen, Denmark.
sample | location | total mass (g) | analysed mass (mg) | % ambrein |
---|---|---|---|---|
S.01 | Pitt Island, New Zealand | 50 | 96 | 92 |
S.02 | 20 | 110 | 83 | |
S.03 | west Sri Lanka | 101 | 188 | 60 |
TEXEL151212 | Texel, Netherlands | 83000 | 92 | 93 |
sample | total retained reads | average retained read length (bp) | total aligned reads (mtDNA) | times of coverage (mtDNA) | total aligned reads (whole genome) | times of coverage (whole genome) |
---|---|---|---|---|---|---|
S.01 | 77 261 083 | 72.9 | 43 | 0.175 | 12 782 | 0.0000546 |
S.02 | 89 486 411 | 71.7 | 2440 | 1.648 | 26 169 | 0.000135 |
S.03 | 71 907 406 | 68.3 | 40 235 | 9.717 | 2 447 082 | 0.00426 |
TEXEL151212 | 92 385 587 | 62.6 | 71 190 | 19.654 | 3 099 642 | 0.00639 |
Sequence analysis was undertaken on the high-performance computing facility at the University of Copenhagen, with FASTQ files processed using the PALEOMIX pipeline (v. 1.2.13) [16]. FastQC v. 0.11.8 [17] was initially used for quality control of raw sequence data. Adapters were trimmed using AdapterRemoval v. 2.3.1 [18], with reads less than 25 bp also removed. Reads were then mapped to reference sequences using BWA [19], also applying mapDamage2.0 [20] for basic degradation quantification, producing alignments with references sequences. ANGSD [21] was then used to produce sequences in FASTA format.
Uncertainty around the origin and the biological mechanisms for the production of ambergris prompted us to consider multiple possible candidate cetacean and pinniped species in sequence analysis. Species identity was inferred by mapping success and phylogenetic relationship to 19 cetacean and pinniped candidate species in NCBI RefSeq (see the electronic supplementary material). These species were selected based on potential suitability as deep-diving marine mammals filling a similar ecological niche to sperm whales, to rule out such species being co-adapted to produce ambergris. Sample sequences were concatenated and aligned using MAFFT v. 7.392 [22]. Phylogenetic tree models were then produced in MEGA X [23] using the maximum-likelihood method with the Hasegawa–Kishino–Yano model [24], with distances estimated by the maximum composite likelihood approach (details of all reference sequences used are included in the electronic supplementary material).
3. Results
The phylogenetic analyses unequivocally supported the sperm whale origin of the four ambergris samples (figure 2; electronic supplementary material, figure S1). Likewise, alignment with the Physeter macrocephalus mitochondrial reference genome from NCBI (NC_002503.2) produced the highest coverage results for all samples of all the alignments made and provides a confident attribution, though with significant variations in success between samples. Sequencing of the sample from a stranded sperm whale (TEXEL151212) produced by far the highest coverage (approx. 20×) for sperm whale mitochondrion, while one of the jetsam samples from Pitt Island (S.01) only yielded approximately 0.2× coverage (see table 2). Alignments with Kogia sima (dwarf sperm whale) and Kogia breviceps (pygmy sperm whale) reference mitochondrial genomes (NC_041303.1, NC_005272.1) also yielded coverage (details in electronic supplementary material), though many highly conserved functional regions are shared between analysed species, resulting in high sequence similarity [25]. However, coverage for kogiid species was typically around a factor of 10 less than for Physeter. Alignment with the P. macrocephalus whole nuclear reference genome (ASM283717v2) was also successful, though this is more apparent in comparisons of total number of reads mapped to the genome. Alignment for Architeuthis dux (giant squid), a reputed common prey of sperm whales (e.g. [26]), was unsuccessful, but this is not their predominant prey [26].
Results from MapDamage show remarkably little inter-sample variation affecting C to T transitions at the 5′ strand ends, though a higher percentage of G to A transitions at the 3′ end exists for S.01 (details in electronic supplementary material, figure S2), indicating possible higher biomolecular degradation. Generally, however, very little chemical modification has occurred, and the distribution of alterations across strands remains uniform and flat.
4. Discussion
This study has demonstrated that three jetsam ambergris samples can confidently be attributed to sperm whale through DNA extraction. While confirmation of a sperm whale origin for jetsam ambergris is not a surprising result, the present study is the first in providing a significant proof-of-concept in retrieving endogenous DNA from ambergris and successfully using it for organism identification. Importantly, the origin of all three jetsam ambergris samples studied here can confidently be identified as sperm whale on the basis of not only genetic alignment success, but also modelling of mitochondrial genomes in phylogenetic relatedness trees, including for a large sample of outgroup marine mammal taxa. However, although all samples analysed here were identified as originating from sperm whale, it is still quite possible that other closely related deep-diving marine mammals (such as the dwarf and pygmy sperm whales) might produce ambergris [4] and have simply still not been recorded as doing so to date.
The predominant cause of the dramatic variation in genetic coverage seen between samples is unclear. Analysis of DNA degradation in mapDamage2.0 shows little correlation with alignment coverages, as might be expected, and there is also little variation between ambrein content in samples that might be expected to contribute to differential DNA preservation. The precise age of the present jetsam samples is unknown, although previous studies have successfully radiocarbon dated other ambergris samples [13]. However, radiocarbon dating of relatively recent samples is problematical owing to the impact of fossil fuel emissions [27], and radiocarbon dates since the increase of anthropic carbon release are unreliable. Producing a consistent degradation rate for G to A transitions in reliably dated older samples might, in future, aid a better understanding of differential DNA damage. Another option for future research might be studies of glutamine deamidation and aspartic acid racemization from analysis of organic peptides possibly also present in ambergris [28,29]. Alternatively, however, intra-sample variation in DNA and ambrein concentration might just as likely account for low coverage in sample S.01, while more recent exposure to sperm whale tissue undoubtedly accounts for the high coverage in beached whale sample TEXEL151212.
The preservational potential of ambrein precipitates for DNA extends not just to endogenous whale genetics, but also to metagenomic coverage of the whale gut microbiome, and potentially also the DNA of their prey. For example, DNA may also remain within partially or undigested squid beaks found in sperm whale faeces [30], and in ambergris [31], which are even theorized to be a pathological cause of ambergris secretion [8]. Understanding of the prokaryotic composition of the microbiome environment in ambergris could also further elucidate the origin of ambergris, particularly in the conversion of squalene to ambrein and the process by which ambergris appears to form in layers of accretion. Further analyses on endogenous DNA retrieval from jetsam ambergris, including also DNA from whale gut microbiota and prey, would yield significantly greater insights into sperm whale ecology, evolution and metabolism.
5. Conclusion
Jetsam ambergris has long been an enigmatic material, subject to discussion and analyses in scientific publications since the eighteenth century [1,31]. This study is the first to our knowledge to present final confirmation of the biological origin of jetsam ambergris samples as sperm whales, through DNA analysis. Beyond this, however, the present study lays out the potential of ambergris as a new source of genetic data related to sperm whales with a considerable longevity across time. Greater elucidation remains to be achieved through the study of the preservational conditions of DNA in ambrein and of the differential effects from multiple factors. However, the potential implications for aiding our understanding of past population dynamics in whales and their ecologically associated taxa may be profound. The oldest-known ambergris found within Pleistocene deposits feature permineralized squid beaks containing amino acids endogenous to squid, which the authors attribute to the preservational capacity of the local sediment [32,33]. Although it is unlikely DNA will be preserved for such an age (1.75 Ma), this finding might also be attributable to the effectiveness of ambergris and ambrein as preservational substrates. A great deal is still unknown about the ecology and adaptation of the marine giants formerly characterized as semi-mythical beasts, and ambergris may now prove a small but significant key to understanding some further aspects of them.
Ethics
All samples were obtained and analysed ethically in this study, part of the DNRF-128 PROTEIOS project. Of the four ambergris samples analysed in this study, three (S.01–S.03) were subsampled from material studied by S.J.R. in prior publications (see above), obtained as jetsam ambergris found on beaches (sample nos ACL102-S.01, ACL103-S.02 and ACL237-S.03). The fourth sample (TEXEL151212) was kindly provided as a subsample of the ambergris specimen TEXEL151212 by Dr A. Oosterbaan on behalf of the EcoMare Museum, Texel, The Netherlands, where the rest of this sample is still held. We are grateful for the museum's permission to use this sample.
Data accessibility
In support of data open-access, all data generated in this study has been made available on the Zenodo data repository [34]: https://doi.org/10.5281/zenodo.3528073.
Authors' contributions
R.M., M.-H.S.S., M.T.O., M.J.C. and S.J.R. participated in the design of the study and contributed to drafting the manuscript. M.-H.S.S. undertook laboratory work, R.M. undertook sequence alignment and computational analyses, and coordinated the preparation of the manuscript. All authors gave final approval for publication and agree to be held accountable for the work described therein.
Competing interests
We declare we have no competing interests.
Funding
Funding was provided by the Danish National Research Foundation grant PROTEIOS (DNRF128).
Acknowledgements
R.M. is grateful to Dr J. A. Samaniego for advice on genomic analyses, and also to M. McCarthy for advice and discussions. We are grateful to Dr A. Oosterbaan (Ecomare Museum, Texel) for sample S.04 and to anonymous collectors for samples S.01–S.03. The authors would like to acknowledge the assistance of the Danish National High-Throughput Sequencing Centre in Illumina data generation.
References
- 1.
Schwediawer D, Banks J . 1783 An account of ambergrise, by Dr. Schwediawer; presented by Sir Joseph Banks, PRS. Phil. Trans. R. Soc. Lond. 73, 226-241. (doi:10.1098/rstl.1783.0015) Link, Google Scholar - 2.
Lederer E . 1949 Chemistry and biochemistry of some mammalian secretions and excretions. J. Chem. Soc. 1949, 2115-2125. (doi:10.1039/jr9490002115) Crossref, ISI, Google Scholar - 3.
Dannenfeldt KH . 1982 Ambergris: the search for its origin. Isis 73, 382-397. (doi:10.1086/353040) Crossref, PubMed, ISI, Google Scholar - 4.
Clarke R . 2006 The origin of ambergris. Lat. Am. J. Aquat. Mamm. 5, 7-21. (doi:10.5597/lajam00087) Crossref, Google Scholar - 5.
Rowland SJ, Sutton PA, Belt ST, Fitzsimmons-Thoss V, Scarlett AG . 2018 Further spectral and chromatographic studies of ambergris. Nat. Prod. Res. 32, 2603-2609. (doi:10.1080/14786419.2018.1428599) Crossref, PubMed, ISI, Google Scholar - 6.
Brito C, Jordão VL, Pierce GJ . 2016 Ambergris as an overlooked historical marine resource: its biology and role as a global economic commodity. J. Mar. Biol. Assoc. UK 96, 585-596. (doi:10.1017/s0025315415000910) ISI, Google Scholar - 7.
Rowland SJ, Sutton PA . 2017 Chromatographic and spectral studies of jetsam and archived ambergris. Nat. Prod. Res. 31, 1752-1757. (doi:10.1080/14786419.2017.1290618) Crossref, PubMed, ISI, Google Scholar - 8.
Lambertsen RH, Kohn BA . 1987 Unusual multisystemic pathology in a sperm whale bull. J. Wildl. Dis. 23, 510-514. (doi:10.7589/0090-3558-23.3.510) Crossref, PubMed, ISI, Google Scholar - 9.
MacLeod CD, Santos MB, Pierce GJ . 2003 Review of data on diets of beaked whales: evidence of niche separation and geographic segregation. J. Mar. Biol. Assoc. UK 83, 651-665. (doi:10.1017/S0025315403007616h) Crossref, ISI, Google Scholar - 10.
Monteiro S, Ferreira M, Vingada JV, López A, Brownlow A, Méndez-Fernandez P . 2015 Application of stable isotopes to assess the feeding ecology of long-finned pilot whale (Globicephala melas) in the northeast Atlantic Ocean. J. Exp. Mar. Bio. Ecol. 465, 56-63. (doi:10.1016/j.jembe.2015.01.007) Crossref, ISI, Google Scholar - 11.
Abend AG, Smith TD . 1997 Differences in stable isotope ratios of carbon and nitrogen between long-finned pilot whales (Globicephala melas) and their primary prey in the western north Atlantic. ICES J. Mar. Sci. 54, 500-503. (doi:10.1006/jmsc.1996.0192) Crossref, ISI, Google Scholar - 12.
Ueda D, Hoshino T, Sato T . 2013 Cyclization of squalene from both termini: identification of an onoceroid synthase and enzymatic synthesis of ambrein. J. Am. Chem. Soc. 135, 18 335-18 338. (doi:10.1021/ja4107226) Crossref, ISI, Google Scholar - 13.
Rowland SJ, Sutton PA, Knowles TDJ . 2019 The age of ambergris. Nat. Prod. Res. 33, 3134-3142. (doi:10.1080/14786419.2018.1523163) Crossref, PubMed, ISI, Google Scholar - 14.
Gilbert MTP 2007 Whole-genome shotgun sequencing of mitochondria from ancient hair shafts. Science 317, 1927-1930. (doi:10.1126/science.1146971) Crossref, PubMed, ISI, Google Scholar - 15.
Carøe C, Gopalakrishnan S, Vinner L, Mak SST, Sinding MHS, Samaniego JA, Wales N, Sicheritz-Pontén T, Gilbert MTP . 2018 Single-tube library preparation for degraded DNA. Methods Ecol. Evol. 9, 410-419. (doi:10.1111/2041-210X.12871) Crossref, ISI, Google Scholar - 16.
Schubert M . 2014 Characterization of ancient and modern genomes by SNP detection and phylogenomic and metagenomic analysis using PALEOMIX. Nat. Protoc. 9, 1056-1082. (doi:10.1038/nprot.2014.063) Crossref, PubMed, ISI, Google Scholar - 17.
Brabraham Informatics. In press. FastQC: A quality control tool for high throughput sequence data. See https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed 3 November 2019). Google Scholar
- 18.
Schubert M, Lindgreen S, Orlando L . 2016 AdapterRemoval v2: rapid adapter trimming, identification, and read merging. BMC Res. Notes 9, 88. (doi:10.1186/s13104-016-1900-2) Crossref, PubMed, Google Scholar - 19.
Li H, Durbin R . 2009 Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754-1760. (doi:10.1093/bioinformatics/btp324) Crossref, PubMed, ISI, Google Scholar - 20.
Jónsson H, Ginolhac A, Schubert M, Johnson PLF, Orlando L . 2013 mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics 29, 1682-1684. (doi:10.1093/bioinformatics/btt193) Crossref, PubMed, ISI, Google Scholar - 21.
Korneliussen TS, Albrechtsen A, Nielsen R . 2014 ANGSD: analysis of next generation sequencing data. BMC Bioinform. 15, 356. (doi:10.1186/s12859-014-0356-4) Crossref, PubMed, ISI, Google Scholar - 22.
Katoh K, Standley DM . 2013 MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772-780. (doi:10.1093/molbev/mst010) Crossref, PubMed, ISI, Google Scholar - 23.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K . 2018 MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547-1549. (doi:10.1093/molbev/msy096) Crossref, PubMed, ISI, Google Scholar - 24.
Hasegawa M, Kishino H, Yano T . 1985 Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 22, 160-174. (doi:10.1007/BF02101694) Crossref, PubMed, ISI, Google Scholar - 25.
Arnason U, Gullberg A, Widegren B . 1993 Cetacean mitochondrial DNA control region: sequences of all extant baleen whales and two sperm whale species. Mol. Biol. Evol. 10, 960-970. (doi:10.1093/oxfordjournals.molbev.a040061) PubMed, ISI, Google Scholar - 26.
Nemoto T, Okiyama M, Iwasaki N, Kikuchi T . 1988 Squid as predators on krill (Euphausia superba) and prey for sperm whales in the Southern Ocean. In Antarctic Ocean and resources variability (ed.Sahrhage D ), pp. 292-296. Berlin, Germany: Springer. Crossref, Google Scholar - 27.
Graven HD . 2015 Impact of fossil fuel emissions on atmospheric radiocarbon and various applications of radiocarbon over this century. Proc. Natl Acad. Sci. USA 112, 9542-9545. (doi:10.1073/pnas.1504467112) Crossref, PubMed, ISI, Google Scholar - 28.
Wilson J, van Doorn NL, Collins MJ . 2012 Assessing the extent of bone degradation using glutamine deamidation in collagen. Anal. Chem. 84, 9041-9048. (doi:10.1021/ac301333t) Crossref, PubMed, ISI, Google Scholar - 29.
Collins MJ, Waite ER, Van Duin ACT . 1999 Predicting protein decomposition: the case of aspartic–acid racemization kinetics. Phil. Trans. R. Soc. Lond. B 354, 51-64. (doi:10.1098/rstb.1999.0359) Link, ISI, Google Scholar - 30.
Smith SC, Whitehead H . 2000 The diet of Galapagos sperm whales Physeter macrocephalus as indicated by fecal sample analysis. Mar. Mamm. Sci. 16, 315-325. (doi:10.1111/j.1748-7692.2000.tb00927.x) Crossref, ISI, Google Scholar - 31.
Dudley P . 1724 II. An essay upon the natural history of Whales, with a particular account of the ambergris found in the sperma ceti Whale. In a letter to the publisher, from the Honourable Paul Dudley, Esq; FRS. Phil. Trans. R. Soc. Lond. 33, 256-269. (doi:10.1098/rstl.1724.0053) Link, Google Scholar - 32.
Baldanza A, Bizzarri R, Famiani F, Monaco P, Pellegrino R, Sassi P . 2013 Enigmatic, biogenically induced structures in Pleistocene marine deposits: a first record of fossil ambergris. Geology 41, 1075-1078. (doi:10.1130/G34731.1) Crossref, ISI, Google Scholar - 33.
Monaco P, Baldanza A, Bizzarri R, Famiani F, Lezzerini M, Sciuto F . 2014 Ambergris cololites of Pleistocene sperm whales from central Italy and description of the new ichnogenus and ichnospecies Ambergrisichnus alleronae. Palaeontol. Electron. 17, 129A. ISI, Google Scholar - 34.
Macleod R, Sinding M-HS, Olsen MT, Collins MJ, Rowland SJ . 2020Data from: DNA preserved in jetsam whale ambergris . Zenodo Open Data Repository. (doi:10.5281/zenodo.3528073) Google Scholar