Genomic insights into the secondary aquatic transition of penguins.

Penguins lost the ability to fly more than 60 million years ago, subsequently evolving a hyper-specialized marine body plan. Within the framework of a genome-scale, fossil-inclusive phylogeny, we identify key geological events that shaped penguin diversification and genomic signatures consistent with widespread refugia/recolonization during major climate oscillations. We further identify a suite of genes potentially underpinning adaptations related to thermoregulation, oxygenation, diving, vision, diet, immunity and body size, which might have facilitated their remarkable secondary transition to an aquatic ecology. Our analyses indicate that penguins and their sister group (Procellariiformes) have the lowest evolutionary rates yet detected in birds. Together, these findings help improve our understanding of how penguins have transitioned to the marine environment, successfully colonizing some of the most extreme environments on Earth.

Genomic signatures of inbreeding in a critically endangered parrot, the kakapo.

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kakapo (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kakapo, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Early sexual dimorphism in the developing gut microbiome of northern elephant seals.

The gut microbiome is an integral part of a species' ecology, but we know little about how host characteristics impact its development in wild populations. Here, we explored the role of such intrinsic factors in shaping the gut microbiome of northern elephant seals (Mirounga angustirostris) during a critical developmental window of 6 weeks after weaning, when the pups stay ashore without feeding. We found substantial sex differences in the early-life gut microbiome, even though males and females could not yet be distinguished morphologically. Sex and age both explained around 15% of the variation in gut microbial beta diversity, while microbial communities sampled from the same individual showed high levels of similarity across time, explaining another 40% of the variation. Only a small proportion of the variation in beta diversity was explained by health status, assessed by full blood counts, but clinically healthy individuals had a greater microbial alpha diversity than their clinically abnormal peers. Across the post-weaning period, the northern elephant seal gut microbiome was highly dynamic. We found evidence for several colonization and extinction events as well as a decline in Bacteroides and an increase in Prevotella, a pattern that has previously been associated with the transition from nursing to solid food. Lastly, we show that genetic relatedness was correlated with gut microbiome similarity in males but not females, again reflecting early sex differences. Our study represents a naturally diet-controlled and longitudinal investigation of how intrinsic factors shape the early gut microbiome in a species with extreme sex differences in morphology and life history.

High-coverage genomes to elucidate the evolution of penguins.

BACKGROUND: Penguins (Sphenisciformes) are a remarkable order of flightless wing-propelled diving seabirds distributed widely across the southern hemisphere. They share a volant common ancestor with Procellariiformes close to the Cretaceous-Paleogene boundary (66 million years ago) and subsequently lost the ability to fly but enhanced their diving capabilities. With ∼20 species among 6 genera, penguins range from the tropical Galápagos Islands to the oceanic temperate forests of New Zealand, the rocky coastlines of the sub-Antarctic islands, and the sea ice around Antarctica. To inhabit such diverse and extreme environments, penguins evolved many physiological and morphological adaptations. However, they are also highly sensitive to climate change. Therefore, penguins provide an exciting target system for understanding the evolutionary processes of speciation, adaptation, and demography. Genomic data are an emerging resource for addressing questions about such processes. RESULTS: Here we present a novel dataset of 19 high-coverage genomes that, together with 2 previously published genomes, encompass all extant penguin species. We also present a well-supported phylogeny to clarify the relationships among penguins. In contrast to recent studies, our results demonstrate that the genus Aptenodytes is basal and sister to all other extant penguin genera, providing intriguing new insights into the adaptation of penguins to Antarctica. As such, our dataset provides a novel resource for understanding the evolutionary history of penguins as a clade, as well as the fine-scale relationships of individual penguin lineages. Against this background, we introduce a major consortium of international scientists dedicated to studying these genomes. Moreover, we highlight emerging issues regarding ensuring legal and respectful indigenous consultation, particularly for genomic data originating from New Zealand Taonga species. CONCLUSIONS: We believe that our dataset and project will be important for understanding evolution, increasing cultural heritage and guiding the conservation of this iconic southern hemisphere species assemblage.

Fur seal microbiota are shaped by the social and physical environment, show mother-offspring similarities and are associated with host genetic quality.

Despite an increasing appreciation of the importance of host-microbe interactions in ecological and evolutionary processes, the factors shaping microbial communities in wild populations remain poorly understood. We therefore exploited a natural experiment provided by two adjacent Antarctic fur seal (Arctocephalus gazella) colonies of high and low social density and combined 16S rRNA metabarcoding with microsatellite profiling of mother-offspring pairs to investigate environmental and genetic influences on skin microbial communities. Seal-associated bacterial communities differed profoundly between the two colonies, despite the host populations themselves being genetically undifferentiated. Consistent with the hypothesis that social stress depresses bacterial diversity, we found that microbial alpha diversity was significantly lower in the high-density colony. Seals from one of the colonies that contained a stream also carried a subset of freshwater-associated bacteria, indicative of an influence of the physical environment. Furthermore, mothers and their offspring shared similar microbial communities, in support of the notion that microbes may facilitate mother-offspring recognition. Finally, a significant negative association was found between bacterial diversity and heterozygosity, a measure of host genetic quality. Our study thus reveals a complex interplay between environmental and host genetic effects, while also providing empirical support for the leash model of host control, which posits that bacterial communities are driven not only by bottom-up species interactions, but also by top-down host regulation. Taken together, our findings have broad implications for understanding host-microbe interactions as well as prokaryotic diversity in general.

Estimating the biodiversity of terrestrial invertebrates on a forested island using DNA barcodes and metabarcoding data.

Invertebrates are a major component of terrestrial ecosystems, however, estimating their biodiversity is challenging. We compiled an inventory of invertebrate biodiversity along an elevation gradient on the temperate forested island of Hauturu, New Zealand, by DNA barcoding of specimens obtained from leaf litter samples and pitfall traps. We compared the barcodes and biodiversity estimates from this data set with those from a parallel DNA metabarcoding analysis of soil from the same locations, and with pre-existing sequences in reference databases, before exploring the use of combined data sets as a basis for estimating total invertebrate biodiversity. We obtained 1,282 28S and 1,610 COI barcodes from a total of 1,947 invertebrate specimens, which were clustered into 247 (28S) and 366 (COI) OTUs, of which ≤ 10% were represented in GenBank. Coleoptera were most abundant (730 sequenced specimens), followed by Hymenoptera, Diptera, Lepidoptera, and Amphipoda. The most abundant OTU from both the 28S (153 sequences) and COI (140 sequences) data sets was an undescribed beetle from the family Salpingidae. Based on the occurrences of COI OTUs along the elevation gradient, we estimated there are ~1,000 arthropod species (excluding mites) on Hauturu, including 770 insects, of which 344 are beetles. A DNA metabarcoding analysis of soil DNA from the same sites resulted in the identification of similar numbers of OTUs in most invertebrate groups compared with the DNA barcoding, but less than 10% of the DNA barcoding COI OTUs were also detected by the metabarcoding analysis of soil DNA. A mark-recapture analysis based on the overlap between these data sets estimated the presence of approximately 6,800 arthropod species (excluding mites) on the island, including ~3,900 insects. Estimates of New Zealand-wide biodiversity for selected arthropod groups based on matching of the COI DNA barcodes with pre-existing reference sequences suggested over 13,200 insect species are present, including 4,000 Coleoptera, 2,200 Diptera, and 2,700 Hymenoptera species, and 1,000 arachnid species (excluding mites). These results confirm that metabarcoding analyses of soil DNA tends to recover different components of terrestrial invertebrate biodiversity compared to traditional invertebrate sampling, but the combined methods provide a novel basis for estimating invertebrate biodiversity.

Genetic sex assignment in wild populations using genotyping-by-sequencing data: A statistical threshold approach.

Establishing the sex of individuals in wild systems can be challenging and often requires genetic testing. Genotyping-by-sequencing (GBS) and other reduced-representation DNA sequencing (RRS) protocols (e.g., RADseq, ddRAD) have enabled the analysis of genetic data on an unprecedented scale. Here, we present a novel approach for the discovery and statistical validation of sex-specific loci in GBS data sets. We used GBS to genotype 166 New Zealand fur seals (NZFS, Arctocephalus forsteri) of known sex. We retained monomorphic loci as potential sex-specific markers in the locus discovery phase. We then used (i) a sex-specific locus threshold (SSLT) to identify significantly male-specific loci within our data set; and (ii) a significant sex-assignment threshold (SSAT) to confidently assign sex in silico the presence or absence of significantly male-specific loci to individuals in our data set treated as unknowns (98.9% accuracy for females; 95.8% for males, estimated via cross-validation). Furthermore, we assigned sex to 86 individuals of true unknown sex using our SSAT and assessed the effect of SSLT adjustments on these assignments. From 90 verified sex-specific loci, we developed a panel of three sex-specific PCR primers that we used to ascertain sex independently of our GBS data, which we show amplify reliably in at least two other pinniped species. Using monomorphic loci normally discarded from large SNP data sets is an effective way to identify robust sex-linked markers for nonmodel species. Our novel pipeline can be used to identify and statistically validate monomorphic and polymorphic sex-specific markers across a range of species and RRS data sets.

Managing shifting species: Ancient DNA reveals conservation conundrums in a dynamic world.

The spread of exotic species represents a major driver of biological change across the planet. While dispersal and colonization are natural biological processes, we suggest that the failure to recognize increasing rates of human-facilitated self-introductions may represent a threat to native lineages. Notably, recent biogeographic analyses have revealed numerous cases of biological range shifts in response to anthropogenic impacts and climate change. In particular, ancient DNA analyses have revealed several cases in which lineages traditionally thought to be long-established "natives" are in fact recent colonizers. Such range expansion events have apparently occurred in response to human-mediated native biodiversity declines and ecosystem change, particularly in recently colonized, isolated ecosystems such as New Zealand. While such events can potentially boost local biodiversity, the spread of exotic lineages may also hasten the decline of indigenous species, so it is essential that conservation managers recognize these rapid biotic shifts.​.

Invader or resident? Ancient-DNA reveals rapid species turnover in New Zealand little penguins.

The expansion of humans into previously unoccupied parts of the globe is thought to have driven the decline and extinction of numerous vertebrate species. In New Zealand, human settlement in the late thirteenth century AD led to the rapid demise of a distinctive vertebrate fauna, and also a number of 'turnover' events where extinct lineages were subsequently replaced by closely related taxa. The recent genetic detection of an Australian little penguin (Eudyptula novaehollandiae) in southeastern New Zealand may potentially represent an additional 'cryptic' invasion. Here we use ancient-DNA (aDNA) analysis and radiocarbon dating of pre-human, archaeological and historical Eudyptula remains to reveal that the arrival of E. novaehollandiae in New Zealand probably occurred between AD 1500 and 1900, following the anthropogenic decline of its sister taxon, the endemic Eudyptula minor. This rapid turnover event, revealed by aDNA, suggests that native species decline can be masked by invasive taxa, and highlights the potential for human-mediated biodiversity shifts.

Coalescent Modelling Suggests Recent Secondary-Contact of Cryptic Penguin Species.

Molecular genetic analyses present powerful tools for elucidating demographic and biogeographic histories of taxa. Here we present genetic evidence showing a dynamic history for two cryptic lineages within Eudyptula, the world's smallest penguin. Specifically, we use a suite of genetic markers to reveal that two congeneric taxa ('Australia' and 'New Zealand') co-occur in southern New Zealand, with only low levels of hybridization. Coalescent modelling suggests that the Australian little penguin only recently expanded into southern New Zealand. Analyses conducted under time-dependent molecular evolutionary rates lend support to the hypothesis of recent anthropogenic turnover, consistent with shifts detected in several other New Zealand coastal vertebrate taxa. This apparent turnover event highlights the dynamic nature of the region's coastal ecosystem.

Evaluating a multigene environmental DNA approach for biodiversity assessment.

BACKGROUND: There is an increasing demand for rapid biodiversity assessment tools that have a broad taxonomic coverage. Here we evaluate a suite of environmental DNA (eDNA) markers coupled with next generation sequencing (NGS) that span the tree of life, comparing them with traditional biodiversity monitoring tools within ten 20×20 meter plots along a 700 meter elevational gradient. RESULTS: From six eDNA datasets (one from each of 16S, 18S, ITS, trnL and two from COI) we identified sequences from 109 NCBI taxonomy-defined phyla or equivalent, ranging from 31 to 60 for a given eDNA marker. Estimates of alpha and gamma diversity were sensitive to the number of sequence reads, whereas beta diversity estimates were less sensitive. The average within-plot beta diversity was lower than between plots for all markers. The soil beta diversity of COI and 18S markers showed the strongest response to the elevational variation of the eDNA markers (COI: r=0.49, p<0.001; 18S: r=0.48, p<0.001). Furthermore pairwise beta diversities for these two markers were strongly correlated with those calculated from traditional vegetation and invertebrate biodiversity measures. CONCLUSIONS: Using a soil-based eDNA approach, we demonstrate that standard phylogenetic markers are capable of recovering sequences from a broad diversity of eukaryotes, in addition to prokaryotes by 16S. The COI and 18S eDNA markers are the best proxies for aboveground biodiversity based on the high correlation between the pairwise beta diversities of these markers and those obtained using traditional methods.

Analysis of a deep transcriptome from the mantle tissue of Patella vulgata Linnaeus (Mollusca: Gastropoda: Patellidae) reveals candidate biomineralising genes.

The gastropod Patella vulgata is abundant on rocky shores in Northern Europe and a significant grazer of intertidal algae. Here we report the application of Illumina sequencing to develop a transcriptome from the adult mantle tissue of P. vulgata. We obtained 47,237,104 paired-end reads of 51 bp, trialled de novo assembly methods and settled on the additive multiple K method followed by redundancy removal as resulting in the most comprehensive assembly. This yielded 29,489 contigs of at least 500 bp in length. We then used three methods to search for candidate genes relevant to biomineralisation: searches via BLAST and Hidden Markov Models for homologues of biomineralising genes from other molluscs, searches for predicted proteins containing tandem repeats and searches for secreted proteins that lacked a transmembrane domain. From the results of these searches we selected 15 contigs for verification by RT-PCR, of which 14 were successfully amplified and cloned. These included homologues of Pif-177/BSMP, Perlustrin, SPARC, AP24, Follistatin-like and Carbonic anhydrase, as well as three containing extensive G-X-Y repeats as found in nacrein. We selected two for further verification by in situ hybridisation, demonstrating expression in the larval shell field. We conclude that de novo assembly of Illumina data offers a cheap and rapid route to a predicted transcriptome that can be used as a resource for further biological study.