Search

Environmental DNA (eDNA) metabarcoding is emerging as a novel, objective tool for monitoring marine metazoan biodiversity. Zooplankton biodiversity in the vast open ocean is currently monitored through continuous plankton recorder (CPR) surveys, using ship‐based bulk plankton sampling and morphological identification. We assessed whether eDNA metabarcoding (2 L filtered seawater) could capture similar Southern Ocean zooplankton biodiversity as conventional CPR bulk sampling (~1,500 L filtered seawater per CPR sample). We directly compared eDNA metabarcoding with (a) conventional morphological CPR sampling and (b) bulk DNA metabarcoding of CPR collected plankton (two transects for each comparison, 40 and 44 paired samples, respectively). A metazoan‐targeted cytochrome c oxidase I (COI) marker was used to characterize species‐level diversity. In the 2 L seawater eDNA samples, this marker amplified large amounts of non‐metazoan picoplanktonic algae, but eDNA metabarcoding still detected up to 1.6 times more zooplankton species than morphologically analysed bulk CPR samples. COI metabarcoding of bulk DNA samples mostly avoided nonmetazoan amplifications and recovered more zooplankton species than eDNA metabarcoding. However, eDNA metabarcoding detected roughly two thirds of metazoan species and identified similar taxa contributing to community differentiation across the subtropical front separating transects. We observed a diurnal pattern in eDNA data for copepods which perform diel vertical migrations, indicating a surprisingly short temporal eDNA signal. Compared to COI, a eukaryote‐targeted 18S ribosomal RNA marker detected a higher proportion, but lower diversity, of metazoans in eDNA. With refinement and standardization of methodology, eDNA metabarcoding could become an efficient tool for monitoring open ocean biodiversity. [ABSTRACT FROM AUTHOR]

Antarctic krill (Euphausia superba) is a highly abundant keystone species of the Southern Ocean ecosystem, directly connecting primary producers to high-trophic level predators. Sex ratios of krill vary remarkably between swarms and this phenomenon is poorly understood, as identification of krill sex relies on external morphological differences that appear late during development. Sex determination mechanisms in krill are unknown, but could include genetic, environmental or parasitic mechanisms. Similarly, virtually nothing is known about molecular sex differentiation. The krill genome has to date not been sequenced, and due to its enormous size and large amount of repetitive elements, it is currently not feasible to develop sex-specific DNA markers. To produce a reliable molecular marker for sex in krill and to investigate molecular sex differentiation we therefore focused on identifying sex-specific transcriptomic differences. Through transcriptomic analysis, we found large gene expression differences between testes and ovaries and identified three genes exclusively expressed in female whole krill from early juvenile stages onwards. The sex-specific expression of these three genes persisted through sexual regression, although our regressed samples originated from a krill aquarium and may differ from wild-regressed krill. Two slightly male-biased genes did not display sufficient expression differences to clearly differentiate sexes. Based on the expression of the three female-specific genes we developed a molecular test that for the first time allows the unambiguous sex determination of krill samples lacking external sex-specific features from juvenile stages onwards, including the sexually regressed krill we examined. [ABSTRACT FROM AUTHOR]

Most seabirds do not have any outward identifiers of their chronological age, so estimation of seabird population age structure generally requires expensive, long-term banding studies. We investigated the potential to use a molecular age biomarker to estimate age in short-tailed shearwaters (Ardenna tenuirostris). We quantified DNA methylation in several A. tenuirostris genes that have shown age-related methylation changes in mammals. In birds ranging from chicks to 21 years of age, bisulphite treated blood and feather DNA was sequenced and methylation levels analysed in 67 CpG sites in 13 target gene regions. From blood samples, five of the top relationships with age were identified in KCNC3 loci (CpG66: R2 = 0.325, p = 0.019). In feather samples ELOVL2 (CpG42: R2 = 0.285, p = 0.00048) and EDARADD (CpG46: R2 = 0.168, p = 0.0067) were also weakly correlated with age. However, the majority of markers had no clear association with age (of 131 comparisons only 12 had a p-value < 0.05) and statistical analysis using a penalised lasso approach did not produce an accurate ageing model. Our data indicate that some age-related signatures identified in orthologous mammalian genes are not conserved in the long-lived short tailed shearwater. Alternative molecular approaches will be required to identify a reliable biomarker of chronological age in these seabirds. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Age is a fundamental aspect of animal ecology, but is difficult to determine in many species. Humpback whales exemplify this as they have a lifespan comparable to humans, mature sexually as early as 4 years and have no reliable visual age indicators after their first year. Current methods for estimating humpback age cannot be applied to all individuals and populations. Assays for human age have recently been developed based on age-induced changes in DNA methylation of specific genes. We used information on age-associated DNA methylation in human and mouse genes to identify homologous gene regions in humpbacks. Humpback skin samples were obtained from individuals with a known year of birth and employed to calibrate relationships between cytosine methylation and age. Seven of 37 cytosines assayed for methylation level in humpback skin had significant age-related profiles. The three most age-informative cytosine markers were selected for a humpback epigenetic age assay. The assay has an R2 of 0.787 ( P = 3.04e−16) and predicts age from skin samples with a standard deviation of 2.991 years. The epigenetic method correctly determined which of parent-offspring pairs is the parent in more than 93% of cases. To demonstrate the potential of this technique, we constructed the first modern age profile of humpback whales off eastern Australia and compared the results to population structure 5 decades earlier. This is the first epigenetic age estimation method for a wild animal species and the approach we took for developing it can be applied to many other nonmodel organisms. [ABSTRACT FROM AUTHOR]

Recent advances in nanofluidic technologies have enabled the use of Integrated Fluidic Circuits (IFCs) for high-throughput Single Nucleotide Polymorphism (SNP) genotyping (GT). In this study, we implemented and validated a relatively low cost nanofluidic system for SNP-GT with and without Specific Target Amplification (STA). As proof of principle, we first validated the effect of input DNA copy number on genotype call rate using well characterised, digital PCR (dPCR) quantified human genomic DNA samples and then implemented the validated method to genotype 45 SNPs in the humpback whale, Megaptera novaeangliae, nuclear genome. When STA was not incorporated, for a homozygous human DNA sample, reaction chambers containing, on average 9 to 97 copies, showed 100% call rate and accuracy. Below 9 copies, the call rate decreased, and at one copy it was 40%. For a heterozygous human DNA sample, the call rate decreased from 100% to 21% when predicted copies per reaction chamber decreased from 38 copies to one copy. The tightness of genotype clusters on a scatter plot also decreased. In contrast, when the same samples were subjected to STA prior to genotyping a call rate and a call accuracy of 100% were achieved. Our results demonstrate that low input DNA copy number affects the quality of data generated, in particular for a heterozygous sample. Similar to human genomic DNA, a call rate and a call accuracy of 100% was achieved with whale genomic DNA samples following multiplex STA using either 15 or 45 SNP-GT assays. These calls were 100% concordant with their true genotypes determined by an independent method, suggesting that the nanofluidic system is a reliable platform for executing call rates with high accuracy and concordance in genomic sequences derived from biological tissue. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM PUBLISHER]

Aim: Delimiting recently diverged species is challenging. During speciation, genetic differentiation may be distributed unevenly across the genome, as different genomic regions can be subject to different selective pressures and evolutionary histories. Reliance on limited numbers of genetic markers that may be underpowered can make species delimitation even more challenging, potentially resulting in taxonomic inconsistencies. Rockhopper penguins of the genus Eudyptes comprise three broadly recognized taxa: northern (E. moseleyi), southern (E. chrysocome) and eastern rockhopper (E. filholi). Their taxonomic status has been controversial for decades, with researchers disagreeing about whether E. chrysocome and E. filholi are distinct species or conspecific. Our goal is to evaluate genome‐wide patterns of divergence to evaluate genetic differentiation and species delimitation in rockhopper penguins, and to assess which mechanisms may underlie previous discordance among nuclear versus mitochondrial analyses. Location: Sub‐Antarctic and temperate coastal regions of the Southern Hemisphere. Methods: We generated reduced‐representation genomic libraries using double digest restriction‐site associated DNA (ddRAD) sequencing to evaluate genetic differentiation, contemporary migration rates and admixture among colonies of rockhopper penguins. Results: The extent of genetic differentiation among the three taxa was consistently higher than population‐level genetic differentiation found within these and other penguin species. There was no evidence of admixture among the three taxa, suggesting the absence of ongoing gene flow among them. Species delimitation analyses based on molecular data, along with other lines of evidence, provide strong support for the taxonomic distinction of three species of rockhopper penguins. Main conclusions: Our results provide strong support for the existence of three distinct species of rockhopper penguins. The recognition of this taxonomic diversity is crucial for the management and conservation of this widely distributed species group. This study illustrates that widespread dispersive seabird lineages lacking obvious morphological differences may nevertheless have complex evolutionary histories and comprise cryptic species diversity. [ABSTRACT FROM AUTHOR]

Species inventories and biodiversity assessments are critical to conservation. Yet cryptic species or recolonizing species can be challenging to detect. DNA metabarcoding provides an alternative tool to identify species that can be difficult to observe during field surveys. We test the efficacy of DNA analysis to identify burrowing petrel species in a rapidly changing landscape, on a remote sub‐Antarctic island following pest eradication. Discarded feathers and scats provided high quality DNA for species identification, assisting in detection of new species arrivals and new breeding sites across Macquarie Island. We highlight how DNA metabarcoding informs species inventories and is a valuable tool to complement seabird field surveys. [ABSTRACT FROM AUTHOR]

AIM: Management of competition with predators is an important consideration for fisheries, particularly within marine protected areas (MPAs) where conservation is a primary objective. We aimed to test whether static no‐take zones within a large, sustainable‐use MPA prevented overlap between gentoo penguins and a krill fishery during two winters with contrasting prey fields. LOCATION: South Georgia, Southwest Atlantic Ocean. METHOD: We used satellite tracking (N = 16, June–September 2018) to describe gentoo penguin movements and distribution and quantified their overlap with the MPA's no‐take zone (NTZ) and the krill fishing grounds. DNA metabarcoding of scats (N = 220, April–September 2018) was used to quantify diet. RESULTS: When krill were at moderate densities and evenly distributed in 2001, gentoo penguins would have spent all of their time within the 12 NM NTZ, but when availability was low in 2018, they spent 46.3% of their time outside the NTZ and 9.6% within the krill fishing grounds. The extension of the NTZ to 30 km in response to this finding would have produced a 14.9% increase in protection for penguins and displaced 4% of fishery hauls. Gentoo penguin diet comprised 25.8% krill, which is lower than in the late 1980s but more than in 2009. MAIN CONCLUSIONS: Gentoo penguins extend their foraging range when krill is scarce, which increases the potential for spatial overlap with the krill fishery during periods of nutritional stress. Current regulations allow for expansion of both extent and catches by the krill fishery and, should this occur, gentoo penguins may face heightened risks from competition. A dynamic ocean management framework, that extends closed areas in response to near real‐time data on penguin movements and krill density estimates, may reduce the potential for competition in this sustainable‐use MPA while allowing a profitable krill fishery. [ABSTRACT FROM AUTHOR]

Aim: The conservation of biodiversity is hampered by data deficiencies, with many new species and subspecies awaiting description or reclassification. Population genomics and ecological niche modelling offer complementary new tools for uncovering functional units of phylogenetic diversity. We hypothesize that phylogenetically delineated lineages of gentoo penguins (Pygoscelis papua) distributed across Antarctica and sub‐Antarctic Islands are subject to spatially explicit ecological conditions that have limited gene flow, facilitating genetic differentiation, and thereby speciation processes. Location: Antarctica and sub‐Antarctic area. Methods: We identify divergent lineages for gentoo penguins using ddRAD‐seq and mtDNA, and generated species distribution models (SDMs) based on terrestrial and marine parameters. Results: Analyses of our genomic data supports the existence of four major lineages of gentoo penguin: (i) spanning the sub‐Antarctic archipelagos north of the Antarctic Polar Front (APF); (ii) Kerguelen Island; (iii) South America; and (iv) across maritime Antarctic and the Scotia Arc archipelagos. The APF, a major current system around Antarctica, acts as the most important barrier separating regional sister lineages. Our ecological analyses spanning both the terrestrial (breeding sites) and marine (feeding sites) realms recover limited niche overlap among the major lineages of gentoo penguin. We observe this pattern to correspond more closely with regional differentiation of marine conditions than to terrestrial macroenvironmental features. Main conclusions: Recognition of regional genetic lineages as discrete evolutionary entities that occupy distinct ecological niches and also differ morphologically should be considered a priority for conservation. Gentoo penguins provide a good example of how conservation policy can be directly impacted by new insights obtained through the integration of larger genomic datasets with novel approaches to ecological modelling. This is particularly pertinent to polar environments that are among the most rapidly changing environments on earth. [ABSTRACT FROM AUTHOR]

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.] [ABSTRACT FROM AUTHOR]

Age structure is a fundamental aspect of animal population biology. Age is strongly related to individual physiological condition, reproductive potential and mortality rate. Currently, there are no robust molecular methods for age estimation in birds. Instead, individuals must be ringed as chicks to establish known‐age populations, which is a labour‐intensive and expensive process. The estimation of chronological age using DNA methylation (DNAm) is emerging as a robust approach in mammals including humans, mice and some non‐model species. Here, we quantified DNAm in whole blood samples from a total of 71 known‐age Short‐tailed shearwaters (Ardenna tenuirostris) using digital restriction enzyme analysis of methylation (DREAM). The DREAM method measures DNAm levels at thousands of CpG dinucleotides throughout the genome. We identified seven CpG sites with DNAm levels that correlated with age. A model based on these relationships estimated age with a mean difference of 2.8 years to known age, based on validation estimates from models created by repeated sampling of training and validation data subsets. Longitudinal observation of individuals re‐sampled over 1 or 2 years generally showed an increase in estimated age (6/7 cases). For the first time, we have shown that epigenetic changes with age can be detected in a wild bird. This approach should be of broad interest to researchers studying age biomarkers in non‐model species and will allow identification of markers that can be assessed using targeted techniques for accurate age estimation in large population studies. [ABSTRACT FROM AUTHOR]

Abstract: DNA metabarcoding is an efficient method for measuring biodiversity, but the process of initiating long‐term DNA‐based monitoring programmes, or integrating with conventional programs, is only starting. In marine ecosystems, plankton surveys using the continuous plankton recorder (CPR) have characterized biodiversity along transects covering millions of kilometres with time‐series spanning decades. We investigated the potential for use of metabarcoding in CPR surveys. Samples (n = 53) were collected in two Southern Ocean transects and metazoans identified using standard microscopic methods and by high‐throughput sequencing of a cytochrome c oxidase subunit I marker. DNA increased the number of metazoan species identified and provided high‐resolution taxonomy of groups problematic in conventional surveys (e.g., larval echinoderms and hydrozoans). Metabarcoding also generally produced more detections than microscopy, but this sensitivity may make cross‐contamination during sampling a problem. In some samples, the prevalence of DNA from large plankton such as krill masked the presence of smaller species. We investigated adding a fixed amount of exogenous DNA to samples as an internal control to allow determination of relative plankton biomass. Overall, the metabarcoding data represent a substantial shift in perspective, making direct integration into current long‐term time‐series challenging. We discuss a number of hurdles that exist for progressing DNA metabarcoding from the current snapshot studies to the requirements of a long‐term monitoring programme. Given the power and continually increasing efficiency of metabarcoding, it is almost certain this approach will play an important role in future plankton monitoring. [ABSTRACT FROM AUTHOR]

In understanding the impact of commercial whaling, it is important to estimate the mixing of low latitude breeding populations on Antarctic feeding grounds, particularly the endangered humpback whale populations of Oceania. This paper estimates the degree of genetic differentiation among the putative populations of Oceania (New Caledonia, Tonga, the Cook Islands and French Polynesia) and Australia (western Australia and eastern Australia) using ten microsatellite loci and mtDNA, assesses the power of the data for a mixed-stock analysis, determines ways to improve statistical power for future studies and estimates the population composition of Antarctic samples collected in 2010 south of New Zealand and eastern Australia. A large proportion of individuals could not be assigned to a population of origin (> 52%) using a posterior probability threshold of > 0.90. The mixed-stock analysis simulations however, produced accurate results with humpback whales reapportioned to their population of origin above the 90% threshold for western Australia, New Caledonia and Oceania grouped using a combined mtDNA and microsatellite dataset. Removing the Cook Islands, considered a transient region for humpback whales, from the simulation analysis increased the ability to reapportion Tonga from 86% to 89% and French Polynesia from 89% to 92%. Breeding ground sample size was found to be a factor influencing the accuracy of population reapportionment whereas increasing the mixture or feeding ground sample size improved the precision of results. The mixed-stock analysis of our Antarctic samples revealed substantial contributions from both eastern Australia (53.2%, 6.8% SE) and New Caledonia (43.7%, 5.5% SE) [with Oceania contributing 46.8% (5.9% SE)] but not western Australia. Despite the need for more samples to improve estimates of population allocation, our study strengthens the emerging genetic and non-genetic evidence that Antarctic waters south of New Zealand and eastern Australia are used by humpback whales from both eastern Australia and the more vulnerable breeding population of New Caledonia, representing Oceania. [ABSTRACT FROM AUTHOR]

The Adélie penguin is the most important animal currently used for ecosystem monitoring in the Southern Ocean. The diet of this species is generally studied by visual analysis of stomach contents; or ratios of isotopes of carbon and nitrogen incorporated into the penguin from its food. There are significant limitations to the information that can be gained from these methods. We evaluated population diet assessment by analysis of food DNA in scats as an alternative method for ecosystem monitoring with Adélie penguins as an indicator species. Scats were collected at four locations, three phases of the breeding cycle, and in four different years. A novel molecular diet assay and bioinformatics pipeline based on nuclear small subunit ribosomal RNA gene (SSU rDNA) sequencing was used to identify prey DNA in 389 scats. Analysis of the twelve population sample sets identified spatial and temporal dietary change in Adélie penguin population diet. Prey diversity was found to be greater than previously thought. Krill, fish, copepods and amphipods were the most important food groups, in general agreement with other Adélie penguin dietary studies based on hard part or stable isotope analysis. However, our DNA analysis estimated that a substantial portion of the diet was gelatinous groups such as jellyfish and comb jellies. A range of other prey not previously identified in the diet of this species were also discovered. The diverse prey identified by this DNA-based scat analysis confirms that the generalist feeding of Adélie penguins makes them a useful indicator species for prey community composition in the coastal zone of the Southern Ocean. Scat collection is a simple and non-invasive field sampling method that allows DNA-based estimation of prey community differences at many temporal and spatial scales and provides significant advantages over alternative diet analysis approaches. [ABSTRACT FROM AUTHOR]

Eastern Rockhopper penguins from Marion Island in the sub-Antarctic region. The cover image relates to the Research Article https://doi.org/10.1111/ddi.13399 "Taxonomy based on limited genomic markers may underestimate species diversity of rockhopper penguins and threaten their conservation" by Frugone et al. [Extracted from the article]

Genetic heterogeneity is a difficulty frequently encountered in the search for genes conferring susceptibility to prostate cancer. To circumvent this issue, we selected a large prostate cancer pedigree for genome-wide linkage analysis from a population that is genetically homogeneous. Selected cases and first-degree relatives were genotyped with Affymetrix 10K SNP arrays, identifying a 14 Mb haplotype on chromosome 5 (5p13–q12) inherited identical-by-descent (IBD) by multiple cases. Microsatellite genotyping of additional deceased case samples confirmed that a total of eight cases inherited the common haplotype (P=0.0017). Re-sequencing of eight prioritised candidate genes in the region in six selected individuals identified 15 SNPs segregating with the IBD haplotype, located within the ITGA2 gene. Three of these polymorphisms were selected for genotyping in an independent Tasmanian data set comprising 127 cases with familial prostate cancer, 412 sporadic cases and 319 unaffected controls. Two were associated with prostate cancer risk: rs3212649 (OR=1.67 (1.07–2.6), P=0.0009) and rs1126643 (OR=1.52 (1.01–2.28), P=0.0088). Significant association was observed in both familial and sporadic prostate cancer. Although the functional SNP remains to be identified, considerable circumstantial evidence, provided by in vivo and in vitro studies, supports a role for ITGA2 in tumour development.European Journal of Human Genetics (2009) 17, 368–377; doi:10.1038/ejhg.2008.171; published online 1 October 2008 [ABSTRACT FROM AUTHOR]

Mesopelagic fish form an important link between zooplankton and higher trophic levels in Southern Ocean food webs, however their diets are poorly known. Most of the dietary information available comes from morphological analysis of stomach contents and to a lesser extent fatty acid and stable isotopes. DNA sequencing could substantially improve our knowledge of mesopelagic fish diets, but has not previously been applied. We used high-throughput DNA sequencing (HTS) of the 18S ribosomal DNA and mitochondrial cytochrome oxidase I (COI) to characterise stomach contents of four myctophid and one bathylagid species collected at the southern extension of the Kerguelen Plateau (southern Kerguelen Axis), one of the most productive regions in the Indian sector of the Southern Ocean. Diets of the four myctophid species were dominated by amphipods, euphausiids and copepods, whereas radiolarians and siphonophores contributed a much greater proportion of HTS reads for Bathylagus sp. Analysis of mitochondrial COI showed that all species preyed on Thysanoessa macrura , but Euphausia superba was only detected in the stomach contents of myctophids. Size-based shifts in diet were apparent, with larger individuals of both bathylagid and myctophid species more likely to consume euphausiids, but we found little evidence for regional differences in diet composition for each species over the survey area. The presence of DNA from coelenterates and other gelatinous prey in the stomach contents of all five species suggests the importance of these taxa in the diet of Southern Ocean mesopelagics has been underestimated to date. Our study demonstrates the use of DNA-based diet assessment to determine the role of mesopelagic fish and their trophic position in the Southern Ocean and inform the development of ecosystem models. [ABSTRACT FROM AUTHOR]

• Scarce genetic differentiation between Royal and macaroni penguins. • Population differentiation associated to distance between colonies and temperatures. • Macaroni penguins from Antarctic and sub-Antarctic exhibited genetic differentiation. The study of systematics in wide-ranging seabirds can be challenging due to the vast geographic scales involved, as well as the possible discordance between molecular, morphological and behavioral data. In the Southern Ocean, macaroni penguins (Eudyptes chrysolophus) are distributed over a circumpolar range including populations in Antarctic and sub-Antarctic areas. Macquarie Island, in its relative isolation, is home to a closely related endemic taxon — the royal penguin (Eudyptes schlegeli), which is distinguishable from E. chrysolophus mainly by facial coloration. Although these sister taxa are widely accepted as representing distinct species based on morphological grounds, the extent of their genome-wide differentiation remains uncertain. In this study, we use genome-wide Single Nucleotide Polymorphisms to test genetic differentiation between these geographically isolated taxa and evaluate the main drivers of population structure among breeding colonies of macaroni/royal penguins. Genetic similarity observed between macaroni and royal penguins suggests they constitute a single evolutionary unit. Nevertheless, royal penguins exhibited a tendency to cluster only with macaroni individuals from Kerguelen Island, suggesting that dispersal occurs mainly between these neighboring colonies. A stepping stone model of differentiation of macaroni/royal populations was further supported by a strong pattern of isolation by distance detected across its whole distribution range, possibly driven by large geographic distances between colonies as well as natal philopatry. However, we also detected intraspecific genomic differentiation between Antarctic and sub-Antarctic populations of macaroni penguins, highlighting the role of environmental factors together with geographic distance in the processes of genetic differentiation between Antarctic and sub-Antarctic waters. [ABSTRACT FROM AUTHOR]