Ron P. Dirks

The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system

Significance Snake venoms are toxic protein cocktails used for prey capture. To investigate the evolution of these complex biological weapon systems, we sequenced the genome of a venomous snake, the king cobra, and assessed the composition of venom gland expressed genes, small RNAs, and secreted venom proteins. We show that regulatory components of the venom secretory system may have evolved from a pancreatic origin and that venom toxin genes were co-opted by distinct genomic mechanisms. After co-option, toxin genes important for prey capture have massively expanded by gene duplication and evolved under positive selection, resulting in protein neofunctionalization. This diverse and dramatic venom-related genomic response seemingly occurs in response to a coevolutionary arms race between venomous snakes and their prey. Snakes are limbless predators, and many species use venom to help overpower relatively large, agile prey. Snake venoms are complex protein mixtures encoded by several multilocus gene families that function synergistically to cause incapacitation. To examine venom evolution, we sequenced and interrogated the genome of a venomous snake, the king cobra (Ophiophagus hannah), and compared it, together with our unique transcriptome, microRNA, and proteome datasets from this species, with data from other vertebrates. In contrast to the platypus, the only other venomous vertebrate with a sequenced genome, we find that snake toxin genes evolve through several distinct co-option mechanisms and exhibit surprisingly variable levels of gene duplication and directional selection that correlate with their functional importance in prey capture. The enigmatic accessory venom gland shows a very different pattern of toxin gene expression from the main venom gland and seems to have recruited toxin-like lectin genes repeatedly for new nontoxic functions. In addition, tissue-specific microRNA analyses suggested the co-option of core genetic regulatory components of the venom secretory system from a pancreatic origin. Although the king cobra is limbless, we recovered coding sequences for all Hox genes involved in amniote limb development, with the exception of Hoxd12. Our results provide a unique view of the origin and evolution of snake venom and reveal multiple genome-level adaptive responses to natural selection in this complex biological weapon system. More generally, they provide insight into mechanisms of protein evolution under strong selection.

Primitive duplicate Hox clusters in the European eel's genome.

The enigmatic life cycle and elongated body of the European eel (Anguilla anguilla L., 1758) have long motivated scientific enquiry. Recently, eel research has gained in urgency, as the population has dwindled to the point of critical endangerment. We have assembled a draft genome in order to facilitate advances in all provinces of eel biology. Here, we use the genome to investigate the eels complement of the Hox developmental transcription factors. We show that unlike any other teleost fish, the eel retains fully populated, duplicate Hox clusters, which originated at the teleost-specific genome duplication. Using mRNA-sequencing and in situ hybridizations, we demonstrate that all copies are expressed in early embryos. Theories of vertebrate evolution predict that the retention of functional, duplicate Hox genes can give rise to additional developmental complexity, which is not immediately apparent in the adult. However, the key morphological innovation elsewhere in the eels life history coincides with the evolutionary origin of its Hox repertoire.

A High-Throughput Screen for Tuberculosis Progression

One-third of the world population is infected with Mycobacterium tuberculosis and multi-drug resistant strains are rapidly evolving. The noticeable absence of a whole organism high-throughput screening system for studying the progression of tuberculosis is fast becoming the bottleneck in tuberculosis research. We successfully developed such a system using the zebrafish Mycobacterium marinum infection model, which is a well-characterized model for tuberculosis progression with biomedical significance, mimicking hallmarks of human tuberculosis pathology. Importantly, we demonstrate the suitability of our system to directly study M. tuberculosis, showing for the first time that the human pathogen can propagate in this vertebrate model, resulting in similar early disease symptoms to those observed upon M. marinum infection. Our system is capable of screening for disease progression via robotic yolk injection of early embryos and visual flow screening of late-stage larvae. We also show that this system can reliably recapitulate the standard caudal vein injection method with a throughput level of 2,000 embryos per hour. We additionally demonstrate the possibility of studying signal transduction leading to disease progression using reverse genetics at high-throughput levels. Importantly, we use reference compounds to validate our system in the testing of molecules that prevent tuberculosis progression, making it highly suited for investigating novel anti-tuberculosis compounds in vivo.

First draft genome sequence of the Japanese eel, Anguilla japonica

The Japanese eel is a much appreciated research object and very important for Asian aquaculture; however, its genomic resources are still limited. We have used a streamlined bioinformatics pipeline for the de novo assembly of the genome sequence of the Japanese eel from raw Illumina sequence reads. The total assembled genome has a size of 1.15 Gbp, which is divided over 323,776 scaffolds with an N50 of 52,849 bp, a minimum scaffold size of 200 bp and a maximum scaffold size of 1.14 Mbp. Direct comparison of a representative set of scaffolds revealed that all the Hox genes and their intergenic distances are almost perfectly conserved between the European and the Japanese eel. The first draft genome sequence of an organism strongly catalyzes research progress in multiple fields. Therefore, the Japanese eel genome sequence will provide a rich resource of data for all scientists working on this important fish species.

Comparison of the Exomes of Common Carp (Cyprinus carpio) and Zebrafish (Danio rerio)

Research on common carp, Cyprinus carpio, is beneficial for zebrafish research because of resources available owing to its large body size, such as the availability of sufficient organ material for transcriptomics, proteomics, and metabolomics. Here we describe the shot gun sequencing of a clonal double-haploid common carp line. The assembly consists of 511891 scaffolds with an N50 of 17 kb, predicting a total genome size of 1.4-1.5 Gb. A detailed analysis of the ten largest scaffolds indicates that the carp genome has a considerably lower repeat coverage than zebrafish, whilst the average intron size is significantly smaller, making it comparable to the fugu genome. The quality of the scaffolding was confirmed by comparisons with RNA deep sequencing data sets and a manual analysis for synteny with the zebrafish, especially the Hox gene clusters. In the ten largest scaffolds analyzed, the synteny of genes is almost complete. Comparisons of predicted exons of common carp with those of the zebrafish revealed only few genes specific for either zebrafish or carp, most of these being of unknown function. This supports the hypothesis of an additional genome duplication event in the carp evolutionary history, which--due to a higher degree of compactness--did not result in a genome larger than that of zebrafish.

Robotic injection of zebrafish embryos for high-throughput screening in disease models

The increasing use of zebrafish larvae for biomedical research applications is resulting in versatile models for a variety of human diseases. These models exploit the optical transparency of zebrafish larvae and the availability of a large genetic tool box. Here we present detailed protocols for the robotic injection of zebrafish embryos at very high accuracy with a speed of up to 2000 embryos per hour. These protocols are benchmarked for several applications: (1) the injection of DNA for obtaining transgenic animals, (2) the injection of antisense morpholinos that can be used for gene knock-down, (3) the injection of microbes for studying infectious disease, and (4) the injection of human cancer cells as a model for tumor progression. We show examples of how the injected embryos can be screened at high-throughput level using fluorescence analysis. Our methods open up new avenues for the use of zebrafish larvae for large compound screens in the search for new medicines.

Multiple kisspeptin receptors in early osteichthyans provide new insights into the evolution of this receptor family.

Deorphanization of GPR54 receptor a decade ago led to the characterization of the kisspeptin receptor (Kissr) in mammals and the discovery of its major role in the brain control of reproduction. While a single gene encodes for Kissr in eutherian mammals including human, other vertebrates present a variable number of Kissr genes, from none in birds, one or two in teleosts, to three in an amphibian, xenopus. In order to get more insight into the evolution of Kissr gene family, we investigated the presence of Kissr in osteichthyans of key-phylogenetical positions: the coelacanth, a representative of early sarcopterygians, the spotted gar, a non-teleost actinopterygian, and the European eel, a member of an early group of teleosts (elopomorphs). We report the occurrence of three Kissr for the first time in a teleost, the eel. As measured by quantitative RT-PCR, the three eel Kissr were differentially expressed in the brain-pituitary-gonadal axis, and differentially regulated in experimentally matured eels, as compared to prepubertal controls. Subfunctionalisation, as shown by these differences in tissue distribution and regulation, may have represented significant evolutionary constraints for the conservation of multiple Kissr paralogs in this species. Furthermore, we identified four Kissr in both coelacanth and spotted gar genomes, providing the first evidence for the presence of four Kissr in vertebrates. Phylogenetic and syntenic analyses supported the existence of four Kissr paralogs in osteichthyans and allowed to propose a clarified nomenclature of Kissr (Kissr-1 to -4) based on these paralogs. Syntenic analysis suggested that the four Kissr paralogs arose through the two rounds of whole genome duplication (1R and 2R) in early vertebrates, followed by multiple gene loss events in the actinopterygian and sarcopterygian lineages. Due to gene loss there was no impact of the teleost-specific whole genome duplication (3R) on the number of Kissr paralogs in current teleosts.

Rapid de novo assembly of the European eel genome from nanopore sequencing reads

We have sequenced the genome of the endangered European eel using the MinION by Oxford Nanopore, and assembled these data using a novel algorithm specifically designed for large eukaryotic genomes. For this 860 Mbp genome, the entire computational process takes two days on a single CPU. The resulting genome assembly significantly improves on a previous draft based on short reads only, both in terms of contiguity (N50 1.2 Mbp) and structural quality. This combination of affordable nanopore sequencing and light weight assembly promises to make high-quality genomic resources accessible for many non-model plants and animals.

First artificial hybrid of the eel species Anguilla australis and Anguilla anguilla

BackgroundStudies on artificial hybridization of different Anguilla species were conducted recently, i.e. female A. australis with male A. dieffenbachii, and female A. japonica with male A. anguilla. The existence of these artificial hybrids was however not demonstrated by independent genetic methods. Two species - A. anguilla and A. australis - that are phylogenetically close but have different sexual maturation times (12-25 weeks and 6-8 weeks, respectively), were expected to produce favourable hybrids for reproduction studies.ResultsA modification of the protocol for the reproduction of Anguilla japonica was used to produce eight-day Anguilla australis larvae, with a success rate of 71.4%. Thus ten out of 14 females produced eggs that could be fertilized, and three batches resulted in mass hatching. Hybrid larvae from female A. australis x male A. Anguilla survived for up to seven days post fertilization (dpf). The early development of the hybrid showed typical characteristics of A. anguilla tail pigmentation at 50 hours post fertilization (hpf), indicating expression of genes derived from the father.ConclusionsIn this paper we describe the first production of hybrid larvae from male A. anguilla and female A. australis and their survival for up to 7 dpf. A species-specific nucleotide difference in the 18 S rDNA gene confirmed that genes from both A. australis and A. anguilla were present in the hybrids. The developmental stages of the hybrid eel embryos and larvae are described using high resolution images. Video footage also indicated a heart beat in 5-dpf larva.

Testing Tuberculosis Drug Efficacy in a Zebrafish High-Throughput Translational Medicine Screen

ABSTRACT The translational value of zebrafish high-throughput screens can be improved when more knowledge is available on uptake characteristics of potential drugs. We investigated reference antibiotics and 15 preclinical compounds in a translational zebrafish-rodent screening system for tuberculosis. As a major advance, we have developed a new tool for testing drug uptake in the zebrafish model. This is important, because despite the many applications of assessing drug efficacy in zebrafish research, the current methods for measuring uptake using mass spectrometry do not take into account the possible adherence of drugs to the larval surface. Our approach combines nanoliter sampling from the yolk using a microneedle, followed by mass spectrometric analysis. To date, no single physicochemical property has been identified to accurately predict compound uptake; our method offers a great possibility to monitor how any novel compound behaves within the system. We have correlated the uptake data with high-throughput drug-screening data from Mycobacterium marinum-infected zebrafish larvae. As a result, we present an improved zebrafish larva drug-screening platform which offers new insights into drug efficacy and identifies potential false negatives and drugs that are effective in zebrafish and rodents. We demonstrate that this improved zebrafish drug-screening platform can complement conventional models of in vivo Mycobacterium tuberculosis-infected rodent assays. The detailed comparison of two vertebrate systems, fish and rodent, may give more predictive value for efficacy of drugs in humans.