The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights.

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Marine environmental DNA biomonitoring reveals seasonal patterns in biodiversity and identifies ecosystem responses to anomalous climatic events.

Marine ecosystems are changing rapidly as the oceans warm and become more acidic. The physical factors and the changes to ocean chemistry that they drive can all be measured with great precision. Changes in the biological composition of communities in different ocean regions are far more challenging to measure because most biological monitoring methods focus on a limited taxonomic or size range. Environmental DNA (eDNA) analysis has the potential to solve this problem in biological oceanography, as it is capable of identifying a huge phylogenetic range of organisms to species level. Here we develop and apply a novel multi-gene molecular toolkit to eDNA isolated from bulk plankton samples collected over a five-year period from a single site. This temporal scale and level of detail is unprecedented in eDNA studies. We identified consistent seasonal assemblages of zooplankton species, which demonstrates the ability of our toolkit to audit community composition. We were also able to detect clear departures from the regular seasonal patterns that occurred during an extreme marine heatwave. The integration of eDNA analyses with existing biotic and abiotic surveys delivers a powerful new long-term approach to monitoring the health of our world's oceans in the context of a rapidly changing climate.

Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture.

The evolution of cetaceans, from their early transition to an aquatic lifestyle to their subsequent diversification, has been the subject of numerous studies. However, although the higher-level relationships among cetacean families have been largely settled, several aspects of the systematics within these groups remain unresolved. Problematic clades include the oceanic dolphins (37 spp.), which have experienced a recent rapid radiation, and the beaked whales (22 spp.), which have not been investigated in detail using nuclear loci. The combined application of high-throughput sequencing with techniques that target specific genomic sequences provide a powerful means of rapidly generating large volumes of orthologous sequence data for use in phylogenomic studies. To elucidate the phylogenetic relationships within the Cetacea, we combined sequence capture with Illumina sequencing to generate data for $\sim $3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from $>$38,000 exons with existing sequences from 11 cetaceans and seven outgroup taxa, we produced the first comprehensive comparative genomic data set for cetaceans, spanning 6,527,596 aligned base pairs (bp) and 89 taxa. Phylogenetic trees reconstructed with maximum likelihood and Bayesian inference of concatenated loci, as well as with coalescence analyses of individual gene trees, produced mostly concordant and well-supported trees. Our results completely resolve the relationships among beaked whales as well as the contentious relationships among oceanic dolphins, especially the problematic subfamily Delphinidae. We carried out Bayesian estimation of species divergence times using MCMCTree and compared our complete data set to a subset of clocklike genes. Analyses using the complete data set consistently showed less variance in divergence times than the reduced data set. In addition, integration of new fossils (e.g., Mystacodon selenensis) indicates that the diversification of Crown Cetacea began before the Late Eocene and the divergence of Crown Delphinidae as early as the Middle Miocene. [Cetaceans; phylogenomics; Delphinidae; Ziphiidae; dolphins; whales.].

How much is enough? Effects of technical and biological replication on metabarcoding dietary analysis.

DNA metabarcoding is increasingly used in dietary studies to estimate diversity, composition and frequency of occurrence of prey items. However, few studies have assessed how technical and biological replication affect the accuracy of diet estimates. This study addresses these issues using the European free-tailed bat Tadarida teniotis, involving high-throughput sequencing of a small fragment of the COI gene in 15 separate faecal pellets and a 15-pellet pool per each of 20 bats. We investigated how diet descriptors were affected by variability among (a) individuals, (b) pellets of each individual and (c) PCRs of each pellet. In addition, we investigated the impact of (d) analysing separate pellets vs. pellet pools. We found that diet diversity estimates increased steadily with the number of pellets analysed per individual, with seven pellets required to detect ~80% of prey species. Most variation in diet composition was associated with differences among individual bats, followed by pellets per individual and PCRs per pellet. The accuracy of frequency of occurrence estimates increased with the number of pellets analysed per bat, with the highest error rates recorded for prey consumed infrequently by many individuals. Pools provided poor estimates of diet diversity and frequency of occurrence, which were comparable to analysing a single pellet per individual, and consistently missed the less common prey items. Overall, our results stress that maximizing biological replication is critical in dietary metabarcoding studies and emphasize that analysing several samples per individual rather than pooled samples produce more accurate results.

Challenges for assessing vertebrate diversity in turbid Saharan water-bodies using environmental DNA.

The Sahara desert is the largest warm desert in the world and a poorly explored area. Small water-bodies occur across the desert and are crucial habitats for vertebrate biodiversity. Environmental DNA (eDNA) is a powerful tool for species detection and is being increasingly used to conduct biodiversity assessments. However, there are a number of difficulties with sampling eDNA from such turbid water-bodies and it is often not feasible to rely on electrical tools in remote desert environments. We trialled a manually powered filtering method in Mauritania, using pre-filtration to circumvent problems posed by turbid water in remote arid areas. From nine vertebrate species expected in the water-bodies, four were detected visually, two via metabarcoding, and one via both methods. Difficulties filtering turbid water led to severe constraints, limiting the sampling protocol to only one sampling point per study site, which alone may largely explain why many of the expected vertebrate species were not detected. The amplification of human DNA using general vertebrate primers is also likely to have contributed to the low number of taxa identified. Here we highlight a number of challenges that need to be overcome to successfully conduct metabarcoding eDNA studies for vertebrates in desert environments in Africa.

High occurrence of jellyfish predation by black-browed and Campbell albatross identified by DNA metabarcoding.

Gelatinous zooplankton are a large component of the animal biomass in all marine environments, but are considered to be uncommon in the diet of most marine top predators. However, the diets of key predator groups like seabirds have conventionally been assessed from stomach content analyses, which cannot detect most gelatinous prey. As marine top predators are used to identify changes in the overall species composition of marine ecosystems, such biases in dietary assessment may impact our detection of important ecosystem regime shifts. We investigated albatross diet using DNA metabarcoding of scats to assess the prevalence of gelatinous zooplankton consumption by two albatross species, one of which is used as an indicator species for ecosystem monitoring. Black-browed and Campbell albatross scats were collected from eight breeding colonies covering the circumpolar range of these birds over two consecutive breeding seasons. Fish was the main dietary item at most sites; however, cnidarian DNA, primarily from scyphozoan jellyfish, was present in 42% of samples overall and up to 80% of samples at some sites. Jellyfish was detected during all breeding stages and consumed by adults and chicks. Trawl fishery catches of jellyfish near the Falkland Islands indicate a similar frequency of jellyfish occurrence in albatross diets in years of high and low jellyfish availability, suggesting jellyfish consumption may be selective rather than opportunistic. Warmer oceans and overfishing of finfish are predicted to favour jellyfish population increases, and we demonstrate here that dietary DNA metabarcoding enables measurements of the contribution of gelatinous zooplankton to the diet of marine predators.

Antarctic krill population genomics: apparent panmixia, but genome complexity and large population size muddy the water.

Antarctic krill (Euphausia superba; hereafter krill) are an incredibly abundant pelagic crustacean which has a wide, but patchy, distribution in the Southern Ocean. Several studies have examined the potential for population genetic structuring in krill, but DNA-based analyses have focused on a limited number of markers and have covered only part of their circum-Antarctic range. We used mitochondrial DNA and restriction site-associated DNA sequencing (RAD-seq) to investigate genetic differences between krill from five sites, including two from East Antarctica. Our mtDNA results show no discernible genetic structuring between sites separated by thousands of kilometres, which is consistent with previous studies. Using standard RAD-seq methodology, we obtained over a billion sequences from >140 krill, and thousands of variable nucleotides were identified at hundreds of loci. However, downstream analysis found that markers with sufficient coverage were primarily from multicopy genomic regions. Careful examination of these data highlights the complexity of the RAD-seq approach in organisms with very large genomes. To characterize the multicopy markers, we recorded sequence counts from variable nucleotide sites rather than the derived genotypes; we also examined a small number of manually curated genotypes. Although these analyses effectively fingerprinted individuals, and uncovered a minor laboratory batch effect, no population structuring was observed. Overall, our results are consistent with panmixia of krill throughout their distribution. This result may indicate ongoing gene flow. However, krill's enormous population size creates substantial panmictic inertia, so genetic differentiation may not occur on an ecologically relevant timescale even if demographically separate populations exist.

Molecular biomarkers for chronological age in animal ecology.

The chronological age of an individual animal predicts many of its biological characteristics, and these in turn influence population-level ecological processes. Animal age information can therefore be valuable in ecological research, but many species have no external features that allow age to be reliably determined. Molecular age biomarkers provide a potential solution to this problem. Research in this area of molecular ecology has so far focused on a limited range of age biomarkers. The most commonly tested molecular age biomarker is change in average telomere length, which predicts age well in a small number of species and tissues, but performs poorly in many other situations. Epigenetic regulation of gene expression has recently been shown to cause age-related modifications to DNA and to cause changes in abundance of several RNA types throughout animal lifespans. Age biomarkers based on these epigenetic changes, and other new DNA-based assays, have already been applied to model organisms, humans and a limited number of wild animals. There is clear potential to apply these marker types more widely in ecological studies. For many species, these new approaches will produce age estimates where this was previously impractical. They will also enable age information to be gathered in cross-sectional studies and expand the range of demographic characteristics that can be quantified with molecular methods. We describe the range of molecular age biomarkers that have been investigated to date and suggest approaches for developing the newer marker types as age assays in nonmodel animal species.

Improving accuracy of DNA diet estimates using food tissue control materials and an evaluation of proxies for digestion bias.

Ecologists are increasingly interested in quantifying consumer diets based on food DNA in dietary samples and high-throughput sequencing of marker genes. It is tempting to assume that food DNA sequence proportions recovered from diet samples are representative of consumer's diet proportions, despite the fact that captive feeding studies do not support that assumption. Here, we examine the idea of sequencing control materials of known composition along with dietary samples in order to correct for technical biases introduced during amplicon sequencing and biological biases such as variable gene copy number. Using the Ion Torrent PGM(©) , we sequenced prey DNA amplified from scats of captive harbour seals (Phoca vitulina) fed a constant diet including three fish species in known proportions. Alongside, we sequenced a prey tissue mix matching the seals' diet to generate tissue correction factors (TCFs). TCFs improved the diet estimates (based on sequence proportions) for all species and reduced the average estimate error from 28 ± 15% (uncorrected) to 14 ± 9% (TCF-corrected). The experimental design also allowed us to infer the magnitude of prey-specific digestion biases and calculate digestion correction factors (DCFs). The DCFs were compared with possible proxies for differential digestion (e.g. fish protein%, fish lipid%) revealing a strong relationship between the DCFs and percent lipid of the fish prey, suggesting prey-specific corrections based on lipid content would produce accurate diet estimates in this study system. These findings demonstrate the value of parallel sequencing of food tissue mixtures in diet studies and offer new directions for future research in quantitative DNA diet analysis.

DNA metabarcoding and the cytochrome c oxidase subunit I marker: not a perfect match.

DNA metabarcoding enables efficient characterization of species composition in environmental DNA or bulk biodiversity samples, and this approach is making significant and unique contributions in the field of ecology. In metabarcoding of animals, the cytochrome c oxidase subunit I (COI) gene is frequently used as the marker of choice because no other genetic region can be found in taxonomically verified databases with sequences covering so many taxa. However, the accuracy of metabarcoding datasets is dependent on recovery of the targeted taxa using conserved amplification primers. We argue that COI does not contain suitably conserved regions for most amplicon-based metabarcoding applications. Marker selection deserves increased scrutiny and available marker choices should be broadened in order to maximize potential in this exciting field of research.

A genomic predictor for age at sexual maturity for mammalian species.

Age at sexual maturity is a key life history trait that can be used to predict population growth rates and develop life history models. In many wild animal species, the age at sexual maturity is not accurately quantified. This results in a reduced ability to accurately model demography of wild populations. Recent studies have indicated the potential for CpG density within gene promoters to be predictive of other life history traits, specifically maximum lifespan. Here, we have developed a machine learning model using gene promoter CpG density to predict the mean age at sexual maturity in mammalian species. In total, 91 genomes were used to identify 101 unique gene promoters predictive of age at sexual maturity across males and females. We found these gene promoters to be most predictive of age at sexual maturity in females (R 2 = 0.881) compared to males (R 2 = 0.758). The median absolute error rate was also found to be lower in females (0.427 years) compared to males (0.785 years). This model provides a novel method for species-level age at sexual maturity prediction without the need for long-term monitoring. This study also highlights a potential epigenetic mechanism for the onset of sexual maturity, indicating the possibility of using epigenetic biomarkers for this important life history trait.

TICI: a taxon-independent community index for eDNA-based ecological health assessment.

Global biodiversity is declining at an ever-increasing rate. Yet effective policies to mitigate or reverse these declines require ecosystem condition data that are rarely available. Morphology-based bioassessment methods are difficult to scale, limited in scope, suffer prohibitive costs, require skilled taxonomists, and can be applied inconsistently between practitioners. Environmental DNA (eDNA) metabarcoding offers a powerful, reproducible and scalable solution that can survey across the tree-of-life with relatively low cost and minimal expertise for sample collection. However, there remains a need to condense the complex, multidimensional community information into simple, interpretable metrics of ecological health for environmental management purposes. We developed a riverine taxon-independent community index (TICI) that objectively assigns indicator values to amplicon sequence variants (ASVs), and significantly improves the statistical power and utility of eDNA-based bioassessments. The TICI model training step uses the Chessman iterative learning algorithm to assign health indicator scores to a large number of ASVs that are commonly encountered across a wide geographic range. New sites can then be evaluated for ecological health by averaging the indicator value of the ASVs present at the site. We trained a TICI model on an eDNA dataset from 53 well-studied riverine monitoring sites across New Zealand, each sampled with a high level of biological replication (n = 16). Eight short-amplicon metabarcoding assays were used to generate data from a broad taxonomic range, including bacteria, microeukaryotes, fungi, plants, and animals. Site-specific TICI scores were strongly correlated with historical stream condition scores from macroinvertebrate assessments (macroinvertebrate community index or MCI; R2 = 0.82), and TICI variation between sample replicates was minimal (CV = 0.013). Taken together, this demonstrates the potential for taxon-independent eDNA analysis to provide a reliable, robust and low-cost assessment of ecological health that is accessible to environmental managers, decision makers, and the wider community.

Marine eDNA sampling from submerged surfaces with paint rollers.

Environmental DNA (eDNA) analyses of species present in marine environments is the most effective biological diversity measurement tool currently available. eDNA sampling methods are an intrinsically important part of the eDNA biodiversity analysis process. Identification and development of eDNA sampling methods that are as rapid, affordable, versatile and practical as possible will improve rates of detection of marine species. Optimal outcomes of eDNA biodiversity surveys come from studies employing high levels of sampling replication, so any methods that make sampling faster and cheaper will improve scientific outcomes. eDNA sampling methods that can be applied more widely will also enable sampling from a greater range of marine surface micro-habitats, resulting in detection of a wider range of organisms. In this study, we compared diversity detection by several methods for sampling eDNA from submerged marine surfaces: polyurethane foam, nylon swabs, microfibre paint rollers, and sediment scoops. All of the methods produced a diverse range of species identifications, with >250 multicellular species represented by eDNA at the study site. We found that widely-available small paint rollers were an effective, readily available and affordable method for sampling eDNA from underwater marine surfaces. This approach enables the sampling of marine eDNA using extended poles, or potentially by remotely operated vehicles, where surface sampling by hand is impractical.

Stranded dolphin stomach contents represent the free-ranging population's diet.

Diet is a fundamental aspect of animal ecology. Cetacean prey species are generally identified by examining stomach contents of stranded individuals. Critical uncertainty in these studies is whether samples from stranded animals are representative of the diet of free-ranging animals. Over two summers, we collected faecal and gastric samples from healthy free-ranging individuals of an extensively studied bottlenose dolphin population. These samples were analysed by molecular prey detection and these data compared with stomach contents data derived from stranded dolphins from the same population collected over 22 years. There was a remarkable consistency in the prey species composition and relative amounts between the two datasets. The conclusions of past stomach contents studies regarding dolphin habitat associations, prey selection and proposed foraging mechanisms are supported by molecular data from live animals and the combined dataset. This is the first explicit test of the validity of stomach contents analysis for accurate population-scale diet determination of an inshore cetacean.

Calibrating epigenetic clocks with training data error.

Animal age data are valuable for management of wildlife populations. Yet, for most species, there is no practical method for determining the age of unknown individuals. However, epigenetic clocks, a molecular-based method, are capable of age prediction by sampling specific tissue types and measuring DNA methylation levels at specific loci. Developing an epigenetic clock requires a large number of samples from animals of known ages. For most species, there are no individuals whose exact ages are known, making epigenetic clock calibration inaccurate or impossible. For many epigenetic clocks, calibration samples with inaccurate age estimates introduce a degree of error to epigenetic clock calibration. In this study, we investigated how much error in the training data set of an epigenetic clock can be tolerated before it resulted in an unacceptable increase in error for age prediction. Using four publicly available data sets, we artificially increased the training data age error by iterations of 1% and then tested the model against an independent set of known ages. A small effect size increase (Cohen's d >0.2) was detected when the error in age was higher than 22%. The effect size increased linearly with age error. This threshold was independent of sample size. Downstream applications for age data may have a more important role in deciding how much error can be tolerated for age prediction. If highly precise age estimates are required, then it may be futile to embark on the development of an epigenetic clock when there is no accurately aged calibration population to work with. However, for other problems, such as determining the relative age order of pairs of individuals, a lower-quality calibration data set may be adequate.

Fish species lifespan prediction from promoter cytosine-phosphate-guanine density.

Lifespan is a key attribute of a species' life cycle and varies extensively among major lineages of animals. In fish, lifespan varies by several orders of magnitude, with reported values ranging from less than 1 year to approximately 400 years. Lifespan information is particularly useful for species management, as it can be used to estimate invasion potential, extinction risk and sustainable harvest rates. Despite its utility, lifespan is unknown for most fish species. This is due to the difficulties associated with accurately identifying the oldest individual(s) of a given species, and/or deriving lifespan estimates that are representative for an entire species. Recently it has been shown that CpG density in gene promoter regions can be used to predict lifespan in mammals and other vertebrates, with variable accuracy across taxa. To improve accuracy of lifespan prediction in a non-mammalian vertebrate group, here we develop a fish-specific genomic lifespan predictor. Our new model includes more than eight times the number of fish species included in the previous vertebrate model (n = 442) and uses fish-specific gene promoters as reference sequences. The model predicts fish lifespan from genomic CpG density alone (measured as CpG observed/expected ratio), explaining 64% of the variance between known and predicted lifespans. The predictions are highly robust to variation in genome quality and are applicable to all classes of fish; a taxonomically diverse and speciose group. The results demonstrate the value of promoter CpG density as a universal predictor of fish lifespan that can applied where empirical data are unavailable, or impracticable to obtain.

Who is eating what: diet assessment using next generation sequencing.

The analysis of food webs and their dynamics facilitates understanding of the mechanistic processes behind community ecology and ecosystem functions. Having accurate techniques for determining dietary ranges and components is critical for this endeavour. While visual analyses and early molecular approaches are highly labour intensive and often lack resolution, recent DNA-based approaches potentially provide more accurate methods for dietary studies. A suite of approaches have been used based on the identification of consumed species by characterization of DNA present in gut or faecal samples. In one approach, a standardized DNA region (DNA barcode) is PCR amplified, amplicons are sequenced and then compared to a reference database for identification. Initially, this involved sequencing clones from PCR products, and studies were limited in scale because of the costs and effort required. The recent development of next generation sequencing (NGS) has made this approach much more powerful, by allowing the direct characterization of dozens of samples with several thousand sequences per PCR product, and has the potential to reveal many consumed species simultaneously (DNA metabarcoding). Continual improvement of NGS technologies, on-going decreases in costs and current massive expansion of reference databases make this approach promising. Here we review the power and pitfalls of NGS diet methods. We present the critical factors to take into account when choosing or designing a suitable barcode. Then, we consider both technical and analytical aspects of NGS diet studies. Finally, we discuss the validation of data accuracy including the viability of producing quantitative data.

Variation in the tyrosinase gene associated with a white humpback whale (Megaptera novaeangliae).

Tyrosinase-negative oculocutaneous albinism (OCA1A) is characterized by lifelong white hair and skin, a phenotype that has been described in most mammalian species worldwide. Tyrosinase is the key enzyme in melanin biosynthesis, and mutations in the tyrosinase gene result in OCA1A. We examined sequence variation at exon 1 of the tyrosinase gene in 66 humpback whale samples collected from the east coast of Australia, including an anomalously white humpback whale known as "Migaloo." We identified 3 novel variants, including a cytosine deletion that results in a premature stop codon in exon 1. The deletion truncates the tyrosinase protein including the putative catalytic domains that are essential for tyrosinase enzymatic activity. Migaloo was homozygous for this deletion, suggesting that the albino phenotype is a consequence of inactive tyrosinase caused by the frameshift in the tyrosinase gene.

A review of applications of environmental DNA for reptile conservation and management.

Reptile populations are in decline globally, with total reptile abundance halving in the past half century, and approximately a fifth of species currently threatened with extinction. Research on reptile distributions, population trends, and trophic interactions can greatly improve the accuracy of conservation listings and planning for species recovery, but data deficiency is an impediment for many species. Environmental DNA (eDNA) can detect species and measure community diversity at diverse spatio-temporal scales, and is especially useful for detection of elusive, cryptic, or rare species, making it potentially very valuable in herpetology. We aim to summarize the utility of eDNA as a tool for informing reptile conservation and management and discuss the benefits and limitations of this approach. A literature review was conducted to collect all studies that used eDNA and focus on reptile ecology, conservation, or management. Results of the literature search are summarized into key discussion points, and the review also draws on eDNA studies from other taxa to highlight methodological challenges and to identify future research directions. eDNA has had limited application to reptiles, relative to other vertebrate groups, and little use in regions with high species richness. eDNA techniques have been more successfully applied to aquatic reptiles than to terrestrial reptiles, and most (64%) of studies focused on aquatic habitats. Two of the four reptilian orders dominate the existing eDNA studies (56% Testudines, 49% Squamata, 5% Crocodilia, 0% Rhynchocephalia). Our review provides direction for the application of eDNA as an emerging tool in reptile ecology and conservation, especially when it can be paired with traditional monitoring approaches. Technologies associated with eDNA are rapidly advancing, and as techniques become more sensitive and accessible, we expect eDNA will be increasingly valuable for addressing key knowledge gaps for reptiles.