New chromosome-scale genomes provide insights into marine adaptations of sea snakes (Hydrophis: Elapidae).

BACKGROUND: Sea snakes underwent a complete transition from land to sea within the last ~ 15 million years, yet they remain a conspicuous gap in molecular studies of marine adaptation in vertebrates. RESULTS: Here, we generate four new annotated sea snake genomes, three of these at chromosome-scale (Hydrophis major, H. ornatus and H. curtus), and perform detailed comparative genomic analyses of sea snakes and their closest terrestrial relatives. Phylogenomic analyses highlight the possibility of near-simultaneous speciation at the root of Hydrophis, and synteny maps show intra-chromosomal variations that will be important targets for future adaptation and speciation genomic studies of this system. We then used a strict screen for positive selection in sea snakes (against a background of seven terrestrial snake genomes) to identify genes over-represented in hypoxia adaptation, sensory perception, immune response and morphological development. CONCLUSIONS: We provide the best reference genomes currently available for the prolific and medically important elapid snake radiation. Our analyses highlight the phylogenetic complexity and conserved genome structure within Hydrophis. Positively selected marine-associated genes provide promising candidates for future, functional studies linking genetic signatures to the marine phenotypes of sea snakes and other vertebrates.

Functional Duplication of the Short-Wavelength-Sensitive Opsin in Sea Snakes: Evidence for Reexpanded Color Sensitivity Following Ancestral Regression.

Color vision is mediated by ancient and spectrally distinct cone opsins. Yet, while there have been multiple losses of opsin genes during the evolution of tetrapods, evidence for opsin gains via functional duplication is extremely scarce. Previous studies have shown that some secondarily marine elapid snakes have acquired expanded "UV-blue" sensitivity via changes at key spectral tuning amino acid sites of the Short-Wavelength Opsin 1 (SWS1) gene. Here, we use elapid reference genomes to show that the molecular origin of this adaptation involved repeated, proximal duplications of the SWS1 gene in the fully marine Hydrophis cyanocinctus. This species possesses four intact SWS1 genes; two of these genes have the ancestral UV sensitivity, and two have a derived sensitivity to the longer wavelengths that dominate marine habitats. We suggest that this remarkable expansion of the opsin repertoire of sea snakes functionally compensates for the ancestral losses of two middle-wavelength opsins in the earliest (dim-light adapted) snakes. This provides a striking contrast to the evolution of opsins during ecological transitions in mammals. Like snakes, early mammals lost two cone photopigments; however, lineages such as bats and cetaceans underwent further opsin losses during their adaptation to dim-light environments.

Horizontal Transposon Transfer and Its Implications for the Ancestral Ecology of Hydrophiine Snakes.

Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids ~30 Mya and have since rapidly diversified into six amphibious, ~60 marine and ~100 terrestrial species that fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route.

Horizontal transfer and subsequent explosive expansion of a DNA transposon in sea kraits (Laticauda).

Transposable elements (TEs) are self-replicating genetic sequences and are often described as important 'drivers of evolution'. This driving force is because TEs promote genomic novelty by enabling rearrangement, and through exaptation as coding and regulatory elements. However, most TE insertions potentially lead to neutral or harmful outcomes, therefore host genomes have evolved machinery to suppress TE expansion. Through horizontal transposon transfer (HTT) TEs can colonize new genomes, and since new hosts may not be able to regulate subsequent replication, these TEs may proliferate rapidly. Here, we describe HTT of the Harbinger-Snek DNA transposon into sea kraits (Laticauda), and its subsequent explosive expansion within Laticauda genomes. This HTT occurred following the divergence of Laticauda from terrestrial Australian elapids approximately 15-25 Mya. This has resulted in numerous insertions into introns and regulatory regions, with some insertions into exons which appear to have altered UTRs or added sequence to coding exons. Harbinger-Snek has rapidly expanded to make up 8-12% of Laticauda spp. genomes; this is the fastest known expansion of TEs in amniotes following HTT. Genomic changes caused by this rapid expansion may have contributed to adaptation to the amphibious-marine habitat.

Demographic analyses of marine and terrestrial snakes (Elapidae) using whole genome sequences.

The question of whether spatial aspects of evolution differ in marine versus terrestrial realms has endured since Ernst Mayr's 1954 essay on marine speciation. Marine systems are often suggested to support larger and more highly connected populations, but quantitative comparisons with terrestrial systems have been lacking. Here, we compared the population histories of marine and terrestrial elapid snakes using the pairwise sequentially Markovian coalescent (PSMC) model to track historical fluctuations in species' effective population sizes (Ne ) from individual whole-genome sequences. To do this we generated a draft genome for the olive sea snake (Aiysurus laevis) and analysed this alongside six published elapid genomes and their sequence reads (marine species Hydrophis curtus, H. melanocephalus and Laticauda laticaudata; terrestrial species Pseudonaja textilis, Naja Naja and Notechis scutatus). Counter to the expectation that marine species should show higher overall Ne and less pronounced fluctuations in Ne , our analyses reveal demographic patterns that are highly variable among species and do not clearly correspond to major ecological divisions. At deeper time intervals, the four marine elapids appear to have experienced relatively stable Ne , while each terrestrial species shows a prominent upturn in Ne starting at ~4 million years ago (Ma) followed by an equally strong decline. However, over the last million years, all seven species show strong and divergent fluctuations. Estimates of Ne in the most recent intervals (~10 kya) are lowest in two of four marine species (H. melanocephalus and Laticauda), and do not correspond to contemporary range sizes in marine or terrestrial taxa.

New Environment, New Invaders-Repeated Horizontal Transfer of LINEs to Sea Snakes.

Although numerous studies have found horizontal transposon transfer (HTT) to be widespread across metazoans, few have focused on HTT in marine ecosystems. To investigate potential recent HTTs into marine species, we searched for novel repetitive elements in sea snakes, a group of elapids which transitioned to a marine habitat at most 18 Ma. Our analysis uncovered repeated HTTs into sea snakes following their marine transition. The seven subfamilies of horizontally transferred LINE retrotransposons we identified in the olive sea snake (Aipysurus laevis) are transcribed, and hence are likely still active and expanding across the genome. A search of 600 metazoan genomes found all seven were absent from other amniotes, including terrestrial elapids, with the most similar LINEs present in fish and marine invertebrates. The one exception was a similar LINE found in sea kraits, a lineage of amphibious elapids which independently transitioned to a marine environment 25 Ma. Our finding of repeated horizontal transfer events into marine snakes greatly expands past findings that the marine environment promotes the transfer of transposons. Transposons are drivers of evolution as sources of genomic sequence and hence genomic novelty. We identified 13 candidate genes for HTT-induced adaptive change based on internal or neighboring HTT LINE insertions. One of these, ADCY4, is of particular interest as a part of the KEGG adaptation pathway "Circadian Entrainment." This provides evidence of the ecological interactions between species influencing evolution of metazoans not only through specific selection pressures, but also by contributing novel genomic material.

Metal resistant bacteria on gold particles: Implications of how anthropogenic contaminants could affect natural gold biogeochemical cycling.

In Earth's near-surface environments, gold biogeochemical cycling involves gold dissolution and precipitation processes, which are partly attributed to bacteria. These biogeochemical processes as well as abrasion (via physical transport) are known to act upon gold particles, thereby resulting in particle transformation including the development of pure secondary gold and altered morphology, respectively. While previous studies have inferred gold biogeochemical cycling from gold particles obtained from natural environments, little is known about how metal contamination in an environment could impact this cycle. Therefore, this study aims to infer how potentially toxic metal contaminants could affect the structure and chemistry of gold particles and therefore the biogeochemical cycling of gold. In doing so, river sediments and gold particles from the De Kaap Valley, South Africa, were analysed using both microanalytical and molecular techniques. Of the metal contaminants detected in the sediment, mercury can chemically interact with gold particles thereby directly altering particle morphology and "erasing" textural evidence indicative of particle transformation. Other metal contaminants (including mercury) indirectly affect gold cycling by exerting a selective pressure on bacteria living on the surface of gold particles. Particles harbouring gold-tolerant bacteria with diverse metal resistant genes, such as Arthrobacter sp. and Pseudomonas sp., contained nearly two times more secondary gold relative to particles harbouring bacteria with less gold-tolerance. In conclusion, metal contaminants can have a direct or indirect effect on gold biogeochemical cycling in natural environments impacted by anthropogenic activity.

Accelerated brain aging towards transcriptional inversion in a zebrafish model of the K115fs mutation of human PSEN2.

BACKGROUND: The molecular changes involved in Alzheimer's disease (AD) progression remain unclear since we cannot easily access antemortem human brains. Some non-mammalian vertebrates such as the zebrafish preserve AD-relevant transcript isoforms of the PRESENILIN genes lost from mice and rats. One example is PS2V, the alternative transcript isoform of the PSEN2 gene. PS2V is induced by hypoxia/oxidative stress and shows increased expression in late onset, sporadic AD brains. A unique, early onset familial AD mutation of PSEN2, K115fs, mimics the PS2V coding sequence suggesting that forced, early expression of PS2V-like isoforms may contribute to AD pathogenesis. Here we use zebrafish to model the K115fs mutation to investigate the effects of forced PS2V-like expression on the transcriptomes of young adult and aged adult brains. METHODS: We edited the zebrafish genome to model the K115fs mutation. To explore its effects at the molecular level, we analysed the brain transcriptome and proteome of young (6-month-old) and aged (24-month-old) wild type and heterozygous mutant female sibling zebrafish. Finally, we used gene co-expression network analysis (WGCNA) to compare molecular changes in the brains of these fish to human AD. RESULTS: Young heterozygous mutant fish show transcriptional changes suggesting accelerated brain aging and increased glucocorticoid signalling. These early changes precede a transcriptional 'inversion' that leads to glucocorticoid resistance and other likely pathological changes in aged heterozygous mutant fish. Notably, microglia-associated immune responses regulated by the ETS transcription factor family are altered in both our zebrafish mutant model and in human AD. The molecular changes we observe in aged heterozygous mutant fish occur without obvious histopathology and possibly in the absence of Aß. CONCLUSIONS: Our results suggest that forced expression of a PS2V-like isoform contributes to immune and stress responses favouring AD pathogenesis. This highlights the value of our zebrafish genetic model for exploring molecular mechanisms involved in AD pathogenesis.

Phototactic tails: Evolution and molecular basis of a novel sensory trait in sea snakes.

Dermal phototaxis has been reported in a few aquatic vertebrate lineages spanning fish, amphibians and reptiles. These taxa respond to light on the skin of their elongate hind-bodies and tails by withdrawing under cover to avoid detection by predators. Here, we investigated tail phototaxis in sea snakes (Hydrophiinae), the only reptiles reported to exhibit this sensory behaviour. We conducted behavioural tests in 17 wild-caught sea snakes of eight species by illuminating the dorsal surface of the tail and midbody skin using cold white, violet, blue, green and red light. Our results confirmed phototactic tail withdrawal in the previously studied Aipysurus laevis, revealed this trait for the first time in A. duboisii and A. tenuis, and suggested that tail photoreceptors have peak spectral sensitivities between blue and green light (457-514 nm). Based on these results, and an absence of photoresponses in five Aipysurus and Hydrophis species, we tentatively infer that tail phototaxis evolved in the ancestor of a clade of six Aipysurus species (comprising 10% of all sea snakes). Quantifying tail damage, we found that the probability of sustaining tail injuries was not influenced by tail phototactic ability in snakes. Gene profiling showed that transcriptomes of both tail skin and body skin lacked visual opsins but contained melanopsin (opn4x) in addition to key genes of the retinal regeneration and phototransduction cascades. This work suggests that a nonvisual photoreceptor (e.g., Gq rhabdomeric) signalling pathway underlies tail phototaxis, and provides candidate gene targets for future studies of this unusual sensory innovation in reptiles.