The tuatara genome reveals ancient features of amniote evolution.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Bisulfite probing reveals DNA structural intricacies.

In recent decades, study of DNA structure has largely been focused on the interrelationships between nucleotides at the level of nearest neighbours. A little-utilized approach to probing structure on a larger scale is non-denaturing bisulfite modification of genomic DNA in conjunction with high-throughput sequencing. This technique revealed a marked gradient in reactivity increasing towards the 5' end of poly-dC:dG mononucleotide repeats as short as two base pairs, suggesting that access of the anion may be greater at these points due to positive-roll bending not predicted by existing models. Consistent with this, the 5' ends of these repeats are strikingly enriched at positions relative to the nucleosome dyad that bend towards the major groove, while their 3' ends tend to sit outside these areas. Mutation rates are also higher at the 5' ends of poly-dC:dG when CpG dinucleotides are excluded. These findings shed light on the mechanisms underlying bending/flexibility of the DNA double helix as well as the sequences that facilitate DNA packaging.

Host population crashes disrupt the diversity of associated marine microbiomes.

Host-associated microbial communities are shaped by myriad factors ranging from host conditions, environmental conditions and other microbes. Disentangling the ecological impact of each of these factors can be particularly difficult as many variables are correlated. Here, we leveraged earthquake-induced changes in host population structure to assess the influence of population crashes on marine microbial ecosystems. A large (7.8 magnitude) earthquake in New Zealand in 2016 led to widespread coastal uplift of up to ~6 m, sufficient to locally extirpate some intertidal southern bull kelp populations. These uplifted populations are slowly recovering, but remain at much lower densities than at nearby, less-uplifted sites. By comparing the microbial communities of the hosts from disturbed and relatively undisturbed populations using 16S rRNA gene amplicon sequencing, we observed that disturbed host populations supported higher functional, taxonomic and phylogenetic microbial beta diversity than non-disturbed host populations. Our findings shed light on microbiome ecological assembly processes, particularly highlighting that large-scale disturbances that affect host populations can dramatically influence microbiome structure. We suggest that disturbance-induced changes in host density limit the dispersal opportunities of microbes, with host community connectivity declining with the density of host populations.

Macroalgal microbiome biogeography is shaped by environmental drivers rather than geographical distance.

BACKGROUND AND AIMS: Contrasting patterns of host and microbiome biogeography can provide insight into the drivers of microbial community assembly. Distance-decay relationships are a classic biogeographical pattern shaped by interactions between selective and non-selective processes. Joint biogeography of microbiomes and their hosts is of increasing interest owing to the potential for microbiome-facilitated adaptation. METHODS: In this study, we examine the coupled biogeography of the model macroalga Durvillaea and its microbiome using a combination of genotyping by sequencing (host) and 16S rRNA amplicon sequencing (microbiome). Alongside these approaches, we use environmental data to characterize the relationship between the microbiome, the host, and the environment. KEY RESULTS: We show that although the host and microbiome exhibit shared biogeographical structure, these arise from different processes, with host biogeography showing classic signs of geographical distance decay, but with the microbiome showing environmental distance decay. Examination of microbial subcommunities, defined by abundance, revealed that the abundance of microbes is linked to environmental selection. As microbes become less common, the dominant ecological processes shift away from selective processes and towards neutral processes. Contrary to expectations, we found that ecological drift does not promote structuring of the microbiome. CONCLUSIONS: Our results suggest that although host macroalgae exhibit a relatively 'typical' biogeographical pattern of declining similarity with increasing geographical distance, the microbiome is more variable and is shaped primarily by environmental conditions. Our findings suggest that the Baas Becking hypothesis of 'everything is everywhere, the environment selects' might be a useful hypothesis to understand the biogeography of macroalgal microbiomes. As environmental conditions change in response to anthropogenic influences, the processes structuring the microbiome of macroalgae might shift, whereas those governing the host biogeography are less likely to change. As a result, increasingly decoupled host-microbe biogeography might be observed in response to such human influences.

Comparison of Reptilian Genomes Reveals Deletions Associated with the Natural Loss of γδ T Cells in Squamates.

T lymphocytes or T cells are key components of the vertebrate response to pathogens and cancer. There are two T cell classes based on their TCRs, αß T cells and γδ T cells, and each plays a critical role in immune responses. The squamate reptiles may be unique among the vertebrate lineages by lacking an entire class of T cells, the γδ T cells. In this study, we investigated the basis of the loss of the γδ T cells in squamates. The genome and transcriptome of a sleepy lizard, the skink Tiliqua rugosa, were compared with those of tuatara, Sphenodon punctatus, the last living member of the Rhynchocephalian reptiles. We demonstrate that the lack of TCRγ and TCRδ transcripts in the skink are due to large deletions in the T. rugosa genome. We also show that tuataras are on a growing list of species, including sharks, frogs, birds, alligators, and platypus, that can use an atypical TCRδ that appears to be a chimera of a TCR chain with an Ab-like Ag-binding domain. Tuatara represents the nearest living relative to squamates that retain γδ T cells. The loss of γδTCR in the skink is due to genomic deletions that appear to be conserved in other squamates. The genes encoding the αßTCR chains in the skink do not appear to have increased in complexity to compensate for the loss of γδ T cells.

Paternal hypoxia exposure primes offspring for increased hypoxia resistance.

BACKGROUND: In a time of rapid environmental change, understanding how the challenges experienced by one generation can influence the fitness of future generations is critically needed. Using tolerance assays and transcriptomic and methylome approaches, we use zebrafish as a model to investigate cross-generational acclimation to hypoxia. RESULTS: We show that short-term paternal exposure to hypoxia endows offspring with greater tolerance to acute hypoxia. We detected two hemoglobin genes that are significantly upregulated by more than 6-fold in the offspring of hypoxia exposed males. Moreover, the offspring which maintained equilibrium the longest showed greatest upregulation in hemoglobin expression. We did not detect differential methylation at any of the differentially expressed genes, suggesting that other epigenetic mechanisms are responsible for alterations in gene expression. CONCLUSIONS: Overall, our findings suggest that an epigenetic memory of past hypoxia exposure is maintained and that this environmentally induced information is transferred to subsequent generations, pre-acclimating progeny to cope with hypoxic conditions.

The association between personalities, alternative breeding strategies and reproductive success in dunnocks.

Although consistent between-individual differences in behaviour (i.e. animal personality) are ubiquitous in natural populations, relatively few studies have examined how personalities influence the formation of social relationships. Yet, behavioural characteristics of both sexes might be key when it comes to pair-bond formation, and cooperation with partners to successfully rear offspring. We here use a wild population of dunnocks (Prunella modularis) to first investigate whether individuals mate nonrandomly (i.e. assortative mating) with regard to four behavioural traits-flight-initiation distance (FID), provisioning, activity and vigilance-that differ in repeatability and have previously been associated with mating patterns and fitness in other species. Second, we test whether an individual's FID is associated with variability in the dunnocks' mating system (i.e. monogamous pairs vs. polygamous groups). Finally, we determine whether FID and provisioning of males and females associate with their reproductive success. We found no statistical support for assortative mating in FID between males and females. Interestingly, in polygamous groups, co-breeding males differed in their FIDs with dominant alpha males having significantly shorter FIDs compared with subordinate beta-males. Moreover, there was evidence for assortative mating in provisioning for alpha males and females in polygamous groups. We also found that male provisioning influenced reproductive success of both sexes, whereas female provisioning rates only positively correlated with her own but not their partner(s) reproductive output. Our results suggest that personality differences may have important implications for social relationships, the emergence of different mating patterns and ultimately reproductive success within populations.

Population Genomics of New Zealand Pouched Lamprey (kanakana; piharau; Geotria australis).

Pouched lamprey (Geotria australis) or kanakana/piharau is a culturally and ecologically significant jawless fish that is distributed throughout Aotearoa New Zealand. Despite its importance, much remains unknown about historical relationships and gene flow between populations of this enigmatic species within New Zealand. To help inform management, we assembled a draft G. australis genome and completed the first comprehensive population genomics analysis of pouched lamprey within New Zealand using targeted gene sequencing (Cyt-b and COI) and restriction site-associated DNA sequencing (RADSeq) methods. Employing 16 000 genome-wide single nucleotide polymorphisms (SNPs) derived from RADSeq (n = 186) and sequence data from Cyt-b (766 bp, n = 94) and COI (589 bp, n = 20), we reveal low levels of structure across 10 sampling locations spanning the species range within New Zealand. F-statistics, outlier analyses, and STRUCTURE suggest a single panmictic population, and Mantel and EEMS tests reveal no significant isolation by distance. This implies either ongoing gene flow among populations or recent shared ancestry among New Zealand pouched lamprey. We can now use the information gained from these genetic tools to assist managers with monitoring effective population size, managing potential diseases, and conservation measures such as artificial propagation programs. We further demonstrate the general utility of these genetic tools for acquiring information about elusive species.

Purifying Selection in Corvids Is Less Efficient on Islands.

Theory predicts that deleterious mutations accumulate more readily in small populations. As a consequence, mutation load is expected to be elevated in species where life-history strategies and geographic or historical contingencies reduce the number of reproducing individuals. Yet, few studies have empirically tested this prediction using genome-wide data in a comparative framework. We collected whole-genome sequencing data for 147 individuals across seven crow species (Corvus spp.). For each species, we estimated the distribution of fitness effects of deleterious mutations and compared it with proxies of the effective population size Ne. Island species with comparatively smaller geographic range sizes had a significantly increased mutation load. These results support the view that small populations have an elevated risk of mutational meltdown, which may contribute to the higher extinction rates observed in island species.

Evolution of the "world's only alpine parrot": Genomic adaptation or phenotypic plasticity, behaviour and ecology?

Climate warming, in particular in island environments, where opportunities for species to disperse are limited, may become a serious threat to cold adapted alpine species. In order to understand how alpine species may respond to a warming world, we need to understand the drivers that have shaped their habitat specialisation and the evolutionary adaptations that allow them to utilize alpine habitats. The endemic, endangered New Zealand kea (Nestor notabilis) is considered the only alpine parrot in the world. As a species commonly found in the alpine zone it may be highly susceptible to climate warming. But is it a true alpine specialist? Is its evolution driven by adaptation to the alpine zone, or is the kea an open habitat generalist that simply uses the alpine zone to, for example, avoid lower lying anthropogenic landscapes? We use whole genome data of the kea and its close, forest adapted sister species, the kaka (Nestor meridionalis) to reconstruct the evolutionary history of both species and identify the functional genomic differences that underlie their habitat specialisations. Our analyses do not identify major functional genomic differences between kea and kaka in pathways associated with high-altitude. Rather, we found evidence that selective pressures on adaptations commonly found in alpine species are present in both Nestor species, suggesting that selection for alpine adaptations has not driven their divergence. Strongly divergent demographic responses to past climate warming between the species nevertheless highlight potential future threats to kea survival in a warming world.

A genome-wide investigation of adaptive signatures in protein-coding genes related to tool behaviour in New Caledonian and Hawaiian crows.

Very few animals habitually manufacture and use tools. It has been suggested that advanced tool behaviour co-evolves with a suite of behavioural, morphological and life history traits. In fact, there are indications for such an adaptive complex in tool-using crows (genus Corvus species). Here, we sequenced the genomes of two habitually tool-using and ten non-tool-using crow species to search for genomic signatures associated with a tool-using lifestyle. Using comparative genomic and population genetic approaches, we screened for signals of selection in protein-coding genes in the tool-using New Caledonian and Hawaiian crows. While we detected signals of recent selection in New Caledonian crows near genes associated with bill morphology, our data indicate that genetic changes in these two lineages are surprisingly subtle, with little evidence at present for convergence. We explore the biological explanations for these findings, such as the relative roles of gene regulation and protein-coding changes, as well as the possibility that statistical power to detect selection in recently diverged lineages may have been insufficient. Our study contributes to a growing body of literature aiming to decipher the genetic basis of recently evolved complex behaviour.

Ovarian fluid proteome variation associates with sperm swimming speed in an externally fertilizing fish.

Sperm velocity is a key trait that predicts the outcome of sperm competition. By promoting or impeding sperm velocity, females can control fertilization via postcopulatory cryptic female choice. In Chinook salmon, ovarian fluid (OF), which surrounds the ova, mediates sperm velocity according to male and female identity, biasing the outcome of sperm competition towards males with faster sperm. Past investigations have revealed proteome variation in OF, but the specific components of OF that differentially mediate sperm velocity have yet to be characterized. Here we use quantitative proteomics to investigate whether OF protein composition explains variation in sperm velocity and fertilization success. We found that OF proteomes from six females robustly clustered into two groups and that these groups are distinguished by the abundance of a restricted set of proteins significantly associated with sperm velocity. Exposure of sperm to OF from females in group I had faster sperm compared to sperm exposed to the OF of group II females. Overall, OF proteins that distinguished between these groups were enriched for vitellogenin and calcium ion interactions. Our findings suggest that these proteins may form the functional basis for cryptic female choice via the biochemical and physiological mediation of sperm velocity.

Dunnock social status correlates with sperm speed, but fast sperm does not always equal high fitness.

Sperm competition theory predicts that males should modulate sperm investment according to their social status. Sperm speed, one proxy of sperm quality, also influences the outcome of sperm competition because fast sperm cells may fertilize eggs before slow sperm cells. We evaluated whether the social status of males predicted their sperm speed in a wild population of dunnocks (Prunella modularis). In addition to the traditional analysis of the average speed of sperm cells per sample, we also analysed subsamples of the fastest sperm cells per sample. In other words, we systematically evaluated the effects of including different numbers of the fastest sperm in our analyses, ranging from the 5-fastest sperm cells to the 100-fastest sperm cells in a sample. We further evaluated whether fitness, defined here as the number of chicks sired per male per breeding season, relates to the sperm speed in the same population. We found that males in monogamous pairings (i.e. low levels of sperm competition), produced the slowest sperm cells, whereas subordinate males in polyandrous male-male coalitions (i.e. high levels of sperm competition) produced the fastest sperm cells. This result was consistent regardless of the number of fastest sperm included in our analyses, but statistical support was conditional on the number of sperm cells included in the analysis. Interestingly, we found no significant relationship between fitness and sperm speed, which suggests that it is possible that the differential mating opportunities across social status levelled out any possible difference. Our study also suggests that it is important to identify biologically meaningful subsets of fastest sperm and cut-offs for inclusions for assessing sperm competition via sperm speed.

Resistance to natural and synthetic gene drive systems.

Scientists are rapidly developing synthetic gene drive elements intended for release into natural populations. These are intended to control or eradicate disease vectors and pests, or to spread useful traits through wild populations for disease control or conservation purposes. However, a crucial problem for gene drives is the evolution of resistance against them, preventing their spread. Understanding the mechanisms by which populations might evolve resistance is essential for engineering effective gene drive systems. This review summarizes our current knowledge of drive resistance in both natural and synthetic gene drives. We explore how insights from naturally occurring and synthetic drive systems can be integrated to improve the design of gene drives, better predict the outcome of releases and understand genomic conflict in general.

Conservation demands safe gene drive.

Interest in developing gene drive systems to control invasive species is growing, with New Zealand reportedly considering the nascent technology as a way to locally eliminate the mammalian pests that threaten its unique flora and fauna. If gene drives successfully eradicated these invasive populations, many would rejoice, but what are the possible consequences? Here, we explore the risk of accidental spread posed by self-propagating gene drive technologies, highlight new gene drive designs that might achieve better outcomes, and explain why we need open and international discussions concerning a technology that could have global ramifications.

The effects of transcription and recombination on mutational dynamics of short tandem repeats.

Short tandem repeats (STR) are ubiquitous components of the genomic architecture of most living organisms. Recent work has highlighted the widespread functional significance of such repeats, particularly around gene regulation, but the mutational processes underlying the evolution of these highly abundant and highly variable sequences are not fully understood. Traditional models assume that strand misalignment during replication is the predominant mechanism, but empirical data suggest the involvement of other processes including recombination and transcription. Despite this evidence, the relative influences of these processes have not previously been tested experimentally on a genome-wide scale. Using deep sequencing, we identify mutations at >200 microsatellites, across 700 generations in replicated populations of two otherwise identical sexual and asexual Saccharomyces cerevisiae strains. Using generalized linear models, we investigate correlates of STR mutability including the nature of the mutation, STR composition and contextual factors including recombination, transcription and replication origins. Sexual capability was not a significant predictor of microsatellite mutability, but, intriguingly, we identify transcription as a significant positive predictor. We also find that STR density is substantially increased in regions neighboring, but not within, recombination hotspots.

Molecular structure of sauropsid ß-keratins from tuatara (Sphenodon punctatus).

The birds and reptiles, collectively known as the sauropsids, can be subdivided phylogenetically into the archosaurs (birds, crocodiles), the testudines (turtles), the squamates (lizards, snakes) and the rhynchocephalia (tuatara). The structural framework of the epidermal appendages from the sauropsids, which include feathers, claws and scales, has previously been characterised by electron microscopy, infrared spectroscopy and X-ray diffraction analyses, as well as by studies of the amino acid sequences of the constituent ß-keratin proteins (also referred to as the corneous ß-proteins). An important omission in this work, however, was the lack of sequence and structural data relating to the epidermal appendages of the rhynchocephalia (tuatara), one of the two branches of the lepidosaurs. Considerable effort has gone into sequencing the tuatara genome and while this is not yet complete, there are now sufficient sequence data for conclusions to be drawn on the similarity of the ß-keratins from the tuatara to those of other members of the sauropsids. These results, together with a comparison of the X-ray diffraction pattern of tuatara claw with those from seagull feather and goanna claw, confirm that there is a common structural plan in the ß-keratins of all of the sauropsids, and not just those that comprise the archosaurs (birds and crocodiles), the testudines (turtles) and the squamates (lizards and snakes).