BRIEF RESEARCH REPORT article

Front. Environ. Sci., 21 June 2022

Sec. Freshwater Science

Volume 10 - 2022 | https://doi.org/10.3389/fenvs.2022.931467

Tell Us a Story Granddad: Age and Origin of an Iconic Australian Lungfish

  • 1. Environomics Future Science Platform, Commonwealth Scientific and Industrial Research Organisation (CSIRO), Indian Ocean Marine Research Centre, Crawley, WA, Australia

  • 2. Department of Regional Development, Manufacturing and Water, Bundaberg, QLD, Australia

  • 3. Seqwater, Ipswich, QLD, Australia

Abstract

The modern discovery of the Australian lungfish (Neoceratodus forsteri) by European settlers in 1870 was considered one of the most important events in natural history by leading international scientists and naturalists of that time. Its distinct evolutionary lineage and unique extant morphological characteristics fostered the romantic zoological moniker “living fossil”. Although N. forsteri were suspected of being long-lived, a reliable estimate of maximum age has remained elusive. Maximum age is critical to inform wildlife management and conservation efforts, including the use of population viability models. To estimate the maximum age for N. forsteri, we sourced DNA from “Granddad”, the presumed longest-living lungfish known in a zoological park and utilised an epigenetic ageing clock developed for N. forsteri. This lungfish specimen was gifted to the Chicago John G. Shedd Aquarium from Australia in 1933 and lived there for 84 years until death in 2017. We estimated the age of Granddad at death to be 109 years (±6 years), confirming N. forsteri as a true centenarian species. Genotyping also revealed the natal origin of Granddad to be the Burnett River, Queensland, Australia, the location of the species’ original discovery in 1870. We demonstrate the application of novel molecular techniques to a unique long-lived and captive-raised specimen, to improve estimates of maximum age for the species, and to identify natal origin. This information will support future conservation efforts for this iconic yet endangered species.

Introduction

The Australian lungfish (Neoceratodus forsteri), known to Australia’s First Nations people by various names including Dhal’la, Djelleh or Theebine (), was brought to the attention of the scientific community in 1870 as the “Burnett Salmon” (), after the river where it was first captured and its characteristic pink flesh. Leading international ichthyologists and naturalists venerated the discovery of this mysterious lungfish with its unique phylogeny that bridged the evolutionary gaps between fish, amphibians and tetrapods (; ). “Living fossil” was afforded to N. forsteri based on the discovery of an extant species that connected zoological orders with divergent evolutionary taxonomy (); for which fossil remains were discovered before living forms (); and for which a general morphological stasis was observed compared to the fossil record (). Classified in the genus Neoceratodus (meaning “new” and “horned tooth”) due to unique dentition that was well represented in fossil records dating back 140 million years (), N. forsteri sits within the air-breathing Dipnoi order, which flourished throughout the planet during the Devonian period (416–359 Mya) (). Today, only 6 lungfish species remain (4 African, 1 South American, and 1 Australian) (; ).

International interest in N. forsteri grew rapidly after their discovery, driving demand for living specimens to be shipped internationally (, , , , ; ). One specimen, which became known as “Granddad” (Figure 1), made the 20-day sea voyage to the United States of America in April 1933 for the opening of Chicago’s World Fair, and remained at the John G. Shedd Aquarium for 84 years until death in 2017. Neoceratodus forsteri have long been suspected of long life span due to their slow growth rates and large maximum adult sizes attained in the wild, with individuals of 100 cm in length shown to range in age from 10 to 75 years (; ). These highly variable juvenile and adult growth rates combined with a lack of accurate and cost-effective ageing methods, have hindered reliable maximum age estimates for the species. At time of shipment, Granddad was already of adult size, being approximately 109 cm and at death measured 120 cm (Mary-Anne Rodgers Pers Comm) thus providing a rare opportunity to investigate the potential maximum age of N. forsteri.

FIGURE 1

The maximum age, longevity, lifespan or life expectancy () of wildlife is often unknown or difficult to estimate (), however, it is critical to modern conservation planning and management approaches such as population viability modelling (). Estimates of maximum age can be affected by sampling error and can vary within populations of a species due to many genomic and environmental factors (; ). Novel molecular methods have the potential to complement quantitative population modelling tools for conservation planning, including epigenetic clocks that measure levels of DNA methylation to predict age (). Epigenetic clocks work by measuring DNA methylation at cytosine-phosphate-guanine (CpG) sites (). As a cell undergoes mitosis, DNA methylation at CpG sites is supposed to be maintained. However, cells can make mistakes and one of the hallmarks of ageing and cell senescence is a decrease in the maintenance of epigenetic modifications including DNA methylation. Several studies use this biological process to make epigenetic clocks to predict the age of a wide variety of vertebrate species (). This is particularly important for species such as N. forsteri that cannot be aged using traditional methods (). The most common form of DNA methylation in vertebrates is the addition of a methyl group to CpG sites (). DNA methylation has recently been used to develop a “Lungfish clock” to predict the age of N. forsteri ().

In Australia, N. forsteri are found in the Burnett, Mary and Brisbane River catchments of south-eastern Queensland (). Genomic studies indicate extant N. forsteri populations originated from the Burnett and Mary rivers with evidence of genetic inputs from the Burnett to the Mary River from the late Pleistocene to ∼12,000 years ago when sea levels rose (; ). In 1896, N. forsteri specimens were translocated by The Royal Society of Queensland into various catchments south of their presumed native range, including the Brisbane River, in response to an apparent rarity of juveniles and a perceived high extinction risk (; ). N. forsteri currently face multiple anthropogenic threats throughout its range (), which have elevated risk of population decline and possible extinction, resulting in recent relisting as “endangered” under the IUCN Red List of Threatened Species (). Concomitantly, numerous studies have focused on N. forsteri, aiming to address key knowledge gaps to inform conservation efforts (; ; ; ; ; ; ; ).

In this study, we aimed to build on conservation efforts by applying novel genomic methods to a unique long-lived specimen of N. forsteri to 1. Produce an updated maximum age estimate for N. forsteri, 2. Assess the maximum age of N. forsteri within the animal kingdom, and 3. Demonstrate methods that reveal the river of origin for this iconic flagship specimen. Results substantially revise the maximum age estimate for N. forsteri, contributing towards conservation and management outcomes for this iconic species in the wild.

Methods

Sample Collection and DNA Extraction

A sample of fin tissue from Granddad was preserved in 70% ethanol after euthanasia. DNA was extracted using the DNeasy Blood and Tissue Kit (QIAGEN) and RNAse A treated (QIAGEN) using the manufacturer’s protocol.

Age Prediction by DNA Methylation

To estimate the age of Granddad we used a previously developed Lungfish clock (). The Lungfish clock uses the DNA methylation levels of 31 cytosine-phosphate-guanine (CpG) sites to determine age. DNA (150 ng) is bisulphite treated as previously described (). Bisulphite treated DNA is then used for multiplex PCR targeting 31 CpG sites (Supplementary Tables S1, S2). The multiplex PCR product is barcoded with Fluidigm barcodes (Supplementary Tables S3, S4) and is sequenced on an Illumina MiSeq. The fastq file was trimmed 15 bp at the 5′ and 3′ ends to remove adaptor sequences (). Trimmed reads were aligned to the zebrafish reference genome (danRer10, Illumina iGenomes) with bismark version 0.20.0 and bowtie2 v2.3.4 with the following options: -bowtie2-N1-L15–bam-p2–score L, −0.6, −0.6 –non_directional (; ). DNA methylation levels were called using the bismark_methylation_extractor function with default parameters. The Lungfish clock was developed and validated using a data set of 141 known age Australian lungfish samples ranging from 1 to 77 years old. The method to determine the age of the Australian lungfish is provided in Supplementary Table S5. The revised maximum lifespan estimate for N. forsteri were compared to longevity records maintained in the Animal Ageing and Longevity Database ().

River Origin Genotyping

Restriction site-associated DNA sequencing (RAD-seq) data from 92 wild-caught Australian lungfish from a previous research project was used to identify the minimum number of single nucleotide polymorphisms (SNPs) for determining the origin of an individual within the species’ distribution (BioProject Accession# PRJNA477902) (). This data set included samples from three populations; the Brisbane River (n = 30); Burnett River (n = 32); and the Mary River (n = 30). Paired-end reads were quality trimmed using trimmomatic v 0.38 with the following options: - PE phred33 ILLUMINACLIP:TruSeq3-PE:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36 (). Trimmed reads were aligned to the Australian lungfish reference genome (neoFor_v3) with bowtie2 v2.3.4 default settings (; ). SNP calling was carried out using the Genome Analysis Toolkit (GATK) v4.1.9.0 best practices pipeline as it has been shown to perform better at differentiating populations compared to other software (; ; ).

SNP Markers for Population Differentiation

Although the three populations can be differentiated by 15,201 loci as previously described, it is not cost-effective to run RAD-seq on one individual (). Here, the challenge is to identify the minimum number of SNPs to identify the natal origin of Granddad. We used a combination of genome-wide association study (GWAS) and machine learning as previously described to identify the minimum SNPs for population differentiation ().

PLINK v1.9 was used to perform association tests between populations (). SNPs with Bonferroni correction p-value < 0.05 were considered significant. The significant SNPs were used as part of a Random Forest model to identify the minimum number of SNPs for population differentiation (). The caret R package was used to carry out the Random Forest modelling with a 10-fold cross validation (). 70% of the samples were randomly assigned to a training data set, with the remaining assigned to an independent testing data set. The performance of the model was assessed with an area under curve (AUC). A principal component analysis (PCA) was also used to assess the differentiation of the populations by the SNPs returned from the Random Forest modelling. All analyses were carried out using R version 4.0.4 ().

Targeted SNP Genotyping

The SNPs returned from the Random Forest modelling were collectively used to design an assay for the MassARRAY System (). The Agena Bioscience MassARRAY System enables genotyping of multiple targeted SNPs on small scale experiments. The MassARRAY genotyping was carried out at the Australian Genome Research Facility (AGRF). Multiplex PCR conditions were carried out as previously described ().

Results

Granddad Age at Death

In total 32,328 reads were aligned to a partitioned genome with an alignment success rate of 97.5%. The level of methylation at each CpG site is given in Supplementary Table S5. Using the model previously developed for Neoceratodus forsteri, the age at death of Granddad was found to be 109 years (103–115 years, within a 95% confidence interval).

The revised maximum age estimate for N. forsteri places the species as the 21st oldest living animal, out of 3,761 longevity entries in the AnAge database (). Of the 952 fish (aquatic vertebrates with gills) records in the database, N. forsteri ranks as the 12th oldest living fish species in the world (Figure 2). This “centenarian guild” includes well-known marine species such as the Greenland shark (Somniosus microcephalus) which is the current record holder at 392 years; several Rock fish species (Sebastes sp) ranging in maximum age between 116 and 205 years; and the Orange Roughy (Hoplostethus atlanticus) at 149 years. Out of the 11 centenarian fish older than N. forsteri, all except one, the Lake Sturgeon (Acipenser fulvescens), inhabit marine waters, or utilise marine habitats for part of their life cycle. At 152 years maximum age, A. fulvescens is the only other obligate freshwater fish within the centenarian guild, however, N. forsteri now holds the record for the longest living freshwater sub-tropical fish species in the world.

FIGURE 2

)] demonstrating the position of N. forsteri as a rare centenarian, and within the top 12 longest living fish, and the second oldest freshwater fish, in the world.

Granddad River of Origin

On average per sample, 1,994,292 reads aligned to the Australian lungfish reference genome with an alignment rate of 96.1%. The GATK returned a total of 48,386 SNPs which were quality checked and passed the default parameters in PLINK. In total, 2,610 SNPs were found to be associated with either one of the three extant river populations (Supplementary Figure S1). These associated SNPs were added to the random forest model to rank SNPs in terms of importance for river identification. The random forest model resulted in 35 SNPs that could be used to differentiate populations (AUC = 1, Supplementary Figure S2, and Supplementary Table S6). The performance in both the training and testing data set of the model was 100%. The random forest model predicted Granddad as originating from the Burnett River (Figure 3). This is also demonstrated with an unsupervised PCA with the 35 SNPs (Figure 3). This analysis highlights three distinctive extant river populations of the Australian lungfish with 35 SNPs.

FIGURE 3

Discussion

Australian lungfish (Neoceratodus forsteri) are considered the oldest surviving species from an ancient lineage of air-breathing, lobe-finned fish but face many modern anthropogenic threats which have led to contemporary population declines and increased extinction risk. Addressing key biological knowledge gaps to inform quantitative ecological risk assessment is critical to conservation efforts. Here we apply a novel epigenetic clock to determine the age of Granddad as 109 (±6) years and identify its natal origin as the Burnett River, Australia. At time of death, Granddad was likely the longest-lived N. forsteri known in a zoological park and confirms the species as a true centenarian. These findings provide a revised maximum age to aid in the conservation management of this endangered species. We also show that genotyping can effectively determine the natal origin of lungfish specimens, demonstrating a useful tool for forensic studies in the potential illegal trade of a CITES (Convention on International Trade in Endangered Species) species.

Recent population viability analyses (PVA) developed by the authors (unpublished data) to inform management of wild N. forsteri populations, used a previous maximum age estimate of 77 years (). Our findings extend the maximum age estimate for N. forsteri by 30% over previous estimates, to 109 (±6) years of age. This increase can affect many vital rates within quantitative population modelling tools, such as mortality probabilities, net reproduction, and fecundity potential, by extending the population age at first maturity through to reproductive senescence. N. forsteri can take at least 15 years to mature (), therefore extending the maximum age estimate from 77 to 109 years effectively adds two additional generations into a population model. Maximum age, however, can be artificially augmented in captive settings through improved husbandry, medical treatment and attenuation of key biological processes such as predation and competition (). In addition, species maximum age estimates derived from single records also come with caution (; ). Although captive animals such as Granddad can provide invaluable insight into species biology and behaviour, maximum age estimates for Australian lungfish, as with any long lived species, should utilise a broad sample of both captive and wild individuals to inform population modelling and other wildlife assessments (; ).

Epigenetic age clocks are also limited by the range of ages used in calibration, with decreasing confidence in estimates for animals outside the calibrated age range (). The Lungfish clock used here included “known-age” and “radiocarbon-aged” N. forsteri up to 77 years (). However, epigenetic clocks still offer a valuable tool for predicting ages beyond the calibrated range, albeit with lower certainty. Epigenetic clocks can also be improved as more samples within target age ranges become available and are included into the calibration data. Improving epigenetic clocks for older age classes is an important scientific activity that zoos and aquariums can support, by establishing an age for living specimens in exhibits that are within the calibrated age range and providing ongoing known-age samples beyond the calibrated range of existing clocks. This is particularly important for N. forsteri that exhibit highly variable growth rates in the wild and have populations heavily skewed towards adult fish (). In addition, the Lungfish clock provides a cost-effective ageing technique for effective management of wild populations through monitoring of sporadic recruitment events ().

Our findings have broader implications beyond new scientific knowledge, by enhancing social context for promoting conservation and awareness of N. forsteri. The benefits of human connection for compelling action to conserve a species cannot be underestimated. Studies show that zoo and aquarium exhibits can improve conservation ethics of visitors which can be further enhanced with specific information about animals that appeal to people such as animal age, longevity, place of origin or extinction status (; ). Non-human charisma is generally afforded to more conspicuous megafauna such as mammals (; ; ), however, certain flagship species, like N. forsteri (), are also powerful attractions for capturing the public interest. Granddad, was arguably the world’s first iconic freshwater flagship species connecting to over 100 million visitors at Chicago’s John G. Shedd Aquarium and rightfully held the mantle for the oldest living fish in a zoological park (). This study identifies the important contribution that public aquariums and zoos can make to both scientific research and community conservation outcomes, while providing an improved educational experience for visitors.

With freshwater biodiversity in decline and recognition that threats are inadequately addressed by policy and management actions (), global biodiversity advocates have called for specific actions to support freshwater biodiversity conservation, including the need to identify and promote freshwater flagship species, such as N. forsteri, in public education programs (; ; ). Granddad the lungfish was a pivotal ambassador for N. forsteri, a species under threat of extinction. Given the centenarian status of N. forsteri, living specimens on exhibit in zoological parks could continue roles as ambassadors for global freshwater biodiversity conservation. Other living specimens of N. forsteri, such as “Methuselah” at the California Academy of Sciences’ Steinhart Aquarium in San Francisco, could also exceed the current record set by Granddad and contribute to science and education.

In life, the N. forsteri known as Granddad was seen by over 100 million visitors at the John G. Shedd Aquarium in Chicago (). Even after death, Granddad’s DNA has contributed to conservation efforts and raising public awareness by revealing new biological insights into the species. When alive, Granddad secured his place as the oldest freshwater fish in a zoological park and possibly the world, and propelled N. forsteri into the rare class of true centenarians within the animal kingdom. In reference to a famous quote by the late biologist Dr Stephen Wainwright, “Nature reveals great secrets in her extreme forms”, Granddad represents one of nature’s most extreme living forms, and our application of new non-lethal molecular-based techniques has revealed some of N. forsteri’s greatest secrets including maximum age and geographic origin of Granddad. These findings will not only facilitate research, management actions and future policy to underpin conservation of wild and captive-raised lungfish, but also inform global conservation of freshwater biodiversity efforts.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/bioproject/, PRJNA477902.

Author contributions

BM performed the epigenetic age, river origin prediction, and reviewed the manuscript. TE and DR prepared the manuscript and sourced biological samples and funding for the genotyping of Granddad.

Funding

The authors would also like to thank the Burnett Mary Regional Group (BMRG) for funding the genotyping of Granddad. BM is supported by the CSIRO Environomics Future Science Platform and the North West Shelf Flatback Turtle Conservation Program.

Acknowledgments

The authors would like to thank all authors and contributors to the Australian lungfish RAD-seq data (BioProject Accession# PRJNA477902) which was used in this study. The authors acknowledge the contribution of staff at the John G. Shedd Aquarium for kindly donating the DNA material to enable this research.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenvs.2022.931467/full#supplementary-material

References

Summary

Keywords

maximum age, long-lived, epigenetics, natal origin, lifespan, longevity, centenarian

Citation

Mayne B, Espinoza T and Roberts D (2022) Tell Us a Story Granddad: Age and Origin of an Iconic Australian Lungfish. Front. Environ. Sci. 10:931467. doi: 10.3389/fenvs.2022.931467

Received

29 April 2022

Accepted

25 May 2022

Published

21 June 2022

Volume

10 - 2022

Edited by

Benjamin Cook, FRC Environmental, Australia

Reviewed by

Peter John Unmack, University of Canberra, Australia

Lucian Pârvulescu, West University of Timișoara, Romania

Updates

Copyright

*Correspondence: David Roberts,

This article was submitted to Freshwater Science, a section of the journal Frontiers in Environmental Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics