Towards complete and error-free genome assemblies of all vertebrate species.

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Embryonic development in the bonnethead (Sphyrna tiburo), a viviparous hammerhead shark.

BACKGROUND: The hammerhead sharks (family Sphyrnidae) are an immediately recognizable group of sharks due to their unique head shape. Though there has long been an interest in hammerhead development, there are currently no explicit staging tables published for any members of the group. The bonnethead Sphyrna tiburo is the smallest member of Sphyrnidae and is abundant in estuarine and nearshore waters in the Gulf of Mexico and Western North Atlantic Ocean. Due to their relative abundance, close proximity to shore, and brief gestation period, it has been possible to collect and document multiple embryonic specimens at progressive stages of development. RESULTS: We present the first comprehensive embryonic staging series for the Bonnethead, a viviparous hammerhead shark. Our stage series covers a period of development from stages that match the vertebrate phylotypic period, from Stage 23, through stages of morphological divergence to complete development at birth-Stage 35). Notably, we use a variety of techniques to document crucial stages that lead to their extreme craniofacial diversity, resulting in the formation of one of the most distinctive characters of any shark species, the cephalofoil or hammer-like head. CONCLUSION: Documenting the development of hard-to-access vertebrates, like this viviparous shark species, offers important information about how new and diverse morphologies arise that otherwise may remain poorly studied. This work will serve as a platform for future comparative developmental research both within sharks and across the phylogeny of vertebrates, underpinning the extreme potential of craniofacial development and morphological diversity in vertebrate animals.

White shark comparison reveals a slender body for the extinct megatooth shark, Otodus megalodon (Lamniformes: Otodontidae).

The megatooth shark, †Otodus megalodon, which likely reached at least 15 m in total length, is an iconic extinct shark represented primarily by its gigantic teeth in the Neogene fossil record. As one of the largest marine carnivores to ever exist, understanding the biology, evolution, and extinction of †O. megalodon is important because it had a significant impact on the ecology and evolution of marine ecosystems that shaped the present-day oceans. Some attempts inferring the body form of †O. megalodon have been carried out, but they are all speculative due to the lack of any complete skeleton. Here we highlight the fact that the previous total body length estimated from vertebral diameters of the extant white shark (Carcharodon carcharias) for an †O. megalodon individual represented by an incomplete vertebral column is much shorter than the sum of anteroposterior lengths of those fossil vertebrae. This factual evidence indicates that †O. megalodon had an elongated body relative to the body of the modern white shark. Although its exact body form remains unknown, this proposition represents the most parsimonious empirical evidence, which is a significant step towards deciphering the body form of †O. megalodon.

Reallocation of the magnificent catshark Proscyllium magnificum Last & Vongpanich, 2004 to the genus Ctenacis Compagno, 1973 (Carcharhiniformes: Proscylliidae).

The magnificent catshark Proscyllium magnificum was described in 2004 based off five specimens collected in the Andaman Sea off Myanmar. It was originally allocated to the genus Proscyllium, but recent molecular analyses suggested it was more closely related to the harlequin catshark Ctenacis fehlmanni from the western Indian Ocean. This study incorporated meristics and external and internal morphology, together with molecular data to reclassify the magnificent catshark as Ctenacis magnificum and provides revised diagnoses for the genera Ctenacis and Proscyllium. Ctenacis consists of two allopatric Indian Ocean species, while Proscyllium is monotypic genus confined to the northwest Pacific. The revised Ctenacis can be distinguished from Proscyllium in having a broader and longer head (head length 21%-23% vs. 16%-18% of total length), distance between pectoral and pelvic bases shorter than head length (vs. greater than head length), more teeth (upper jaw with 80-86 vs. 46-62 tooth files), and a complex colour pattern of dark reddish-brown blotches and saddles (vs. colour pattern of small black spots). A revised key to the genera of proscylliids and species of Ctenacis is provided.

Mitogenomic evidence of population differentiation of thorny skate, Amblyraja radiata, in the North Atlantic.

Management of thorny skate (Amblyraja radiata) in the Northwest Atlantic has posed a conservation dilemma for several decades due to the species' lack of response to strong conservation efforts in the US Gulf of Maine and the Canadian Scotian Shelf, confusion over the relationship between two reproductive size morphs of differing life histories that are sympatric in the Northwest Atlantic, and conflicting data on regional population connectivity throughout the species' broader range. To better assess potential A. radiata regional population differentiation and genetic links to life-history variation, we analysed complete mitochondrial genome sequences from 527 specimens collected across the species' North Atlantic geographic range, with particular emphasis on the Northwest Atlantic region. A high level of genetic diversity was evident across the North Atlantic, but significant genetic differentiation was identified between specimens inhabiting the Northwest (Gulf of Maine and Newfoundland) and Northeast (Greenland, Iceland, North Sea, and Arctic Circle) Atlantic. In the Northwest Atlantic, significant differentiation between the Gulf of Maine and Newfoundland regions was revealed; however, the overall level of differentiation was very low. No genetic difference was identified between the large and small reproductive morphs. The results of this study advance our understanding of A. radiata population structure in the North Atlantic but do not resolve all the questions confounding our understanding of the species' biology and evolutionary history.

Deciphering reticulate evolution of the largest group of polyploid vertebrates, the subfamily cyprininae (Teleostei: Cypriniformes).

Despite the rarity of polyploidy in animals, some groups with polyploid species exhibit complicated and interesting patterns of reticulate evolution. Here we focus on fishes in the subfamily Cyprininae, the largest polyploid group of vertebrates. The large number of polyploid taxa poses significant challenges for phylogenetic and evolutionary studies on this subfamily. In this study, we cloned and sequenced three single-copy nuclear loci to investigate the evolution of polyploidy in the Cyprininae. Topologies of nuclear gene trees were compared with a newly reconstructed mitochondrial tree. The data provided herein corroborate the hybrid origins of the tribes Torini, Cyprinini, Spinibarbini, Barbini, and also Probarbini. Based on results from this study and previous studies, we hypothesize that at least 13 independent polyploidization events have occurred during the evolution of the Cyprininae. We offer hypotheses on the origin of each polyploid group and show that a diploid group or the diploid ancestor of a polyploid group might have served as progenitor of one or two other polyploid groups. The evolutionary history of Cyprinine (since its first divergence) can be divided into three stages: the "Diploid stage" (69.2-43.4 Ma or million years ago), the "Tetraploidization stage" (43.4-18.9 Ma), and the "Hexaploidization stage" (18.9 Ma to present). The second stage is when all tetraploidization events happened, while the last stage is when all hexaploidization events and most genus- or species-specific polyploidization events occurred. The post-polyploidization dynamics of polyploid groups are complicated and warrant more genomic level studies. We showed that the subfamily Cyprininae may represent a more complicated polyploid system than most, if not all, other vertebrates and some plants, if one or more of the following factors are considered: number of polyploid species, number of different ploidy levels, and number and type of independent polyploidization events.

Integrative taxonomy identifies a new stingray species of the genus Hypanus Rafinesque, 1818 (Dasyatidae, Myliobatiformes), from the Tropical Southwestern Atlantic.

An integrative approach by the congruence of genetics, morphology and ecological niche modelling (ENM) was used to delimit a new species of Hypanus (Rafinesque, 1818), a recently resurrected genus of marine stingrays comprising eight species, five of which occur in the western Atlantic. The species with the widest distribution, Hypanus americanus (Hildebrand and Schroeder, 1928), from the northeastern coast of the United States to southeastern Brazil, was demonstrated to be paraphyletic based on protein-coding mitochondrial genome analyses. This data set also indicates that the genetic distance between the new species Hypanus berthalutzae sp. nov. and its three closely related species (H. americanus, H. longus and H. rudis) varies from 0.82% to 3.14%. In addition, Bayesian Analysis of Population Similarity using the mitochondrial gene mt-nd2 supports the separation of H. berthalutzae sp. nov. (southwestern Atlantic) from its sister species H. rudis (eastern Atlantic). Similarly, morphological and morphometric analyses corroborated four morphotypes within the H. americanus species group and indicated the ventral caudal fold height and length and interspiracular and interorbital lengths as useful measurements to distinguish among them. Claspers of adult males also exhibit morphological differences among species. The ENM agreed with molecular and morphological analyses and delimits the distribution of H. berthalutzae sp. nov. to shallow areas close to shore along the Brazilian coast, from the mouth of the Amazon River to São Paulo State, including the northeastern oceanic islands, suggesting that the great outflow of fresh water and sediments and the Mid-Atlantic Ridge might act as barriers. The integration of these data to describe a new species provides information relevant to their conservation status, because all species of the H. americanus species group are under the "data-deficient" status.

Demographic inferences after a range expansion can be biased: the test case of the blacktip reef shark (Carcharhinus melanopterus).

The evolutionary history of species is a dynamic process as they modify, expand, and contract their spatial distributions over time. Range expansions (REs) occur through a series of founder events that are followed by migration among neighboring demes. The process usually results in structured metapopulations and leaves a distinct signature in the genetic variability of species. Explicitly modeling the consequences of complex demographic events such as REs is computationally very intensive. Here we propose an an alternative approach that requires less computational effort than a comprehensive RE model, but that can recover the demography of species undergoing a RE, by combining spatially explicit modelling with simplified but realistic metapopulation models. We examine the demographic and colonization history of Carcharhinus melanopterus, an abundant reef-associated shark, as a test case. We first used a population genomics approach to statistically confirm the occurrence of a RE in C. melanopterus, and identify its origin in the Indo-Australian Archipelago. Spatial genetic modelling identified two waves of stepping-stone colonization: an eastward wave moving through the Pacific and a westward one moving through the Indian Ocean. We show that metapopulation models best describe the demographic history of this species and that not accounting for this may lead to incorrectly interpreting the observed genetic variation as signals of widespread population bottlenecks. Our study highlights insights that can be gained about demography by coupling metapopulation models with spatial modeling and underscores the need for cautious interpretation of population genetic data when advancing conservation priorities.

Phylogenomic analysis on the exceptionally diverse fish clade Gobioidei (Actinopterygii: Gobiiformes) and data-filtering based on molecular clocklikeness.

The use of genome-scale data to infer phylogenetic relationships has gained in popularity in recent years due to the progress made in target-gene capture and sequencing techniques. Data filtering, the approach of excluding data inconsistent with the model from analyses, presumably could alleviate problems caused by systematic errors in phylogenetic inference. Different data filtering criteria, such as those based on evolutionary rate and molecular clocklikeness as well as others have been proposed for selecting useful phylogenetic markers, yet few studies have tested these criteria using phylogenomic data. We developed a novel set of single-copy nuclear coding markers to capture thousands of target genes in gobioid fishes, a species-rich lineages of vertebrates, and tested the effects of data-filtering methods based on substitution rate and molecular clocklikeness while attempting to control for the compounding effects of missing data and variation in locus length. We found that molecular clocklikeness was a better predictor than overall substitution rate for phylogenetic usefulness of molecular markers in our study. In addition, when the 100 best ranked loci for our predictors were concatenated and analyzed using maximum likelihood, or combined in a coalescent-based species-tree analysis, the resulting trees showed a well-resolved topology of Gobioidei that mostly agrees with previous studies. However, trees generated from the 100 least clocklike frequently recovered conflicting, and in some cases clearly erroneous topologies with strong support, thus indicating strong systematic biases in those datasets. Collectively these results suggest that data filtering has the potential improve the performance of phylogenetic inference when using both a concatenation approach as well as methods that rely on input from individual gene trees (i.e. coalescent species-tree approaches), which may be preferred in scenarios where incomplete lineage sorting is likely to be an issue.

Cranial morphology in Mollisquama sp. (Squaliformes; Dalatiidae) and patterns of cranial evolution in dalatiid sharks.

Dalatiid sharks are members of a family of predominantly small, midwater meso- and bathypelagic chondrichthyans. The family is notable for both its number of monotypic genera and high morphological disparity. Three of the seven dalatiid genera are known only from holotype specimens (Mollisquama parini) or from only a handful of specimens (Euprotomicroides zantedeschia, Heteroscymnoides marleyi), with the only detailed anatomical work consistent across all taxa being studies of dentition. Here, we present detailed anatomical description of the second-ever specimen of Mollisquama (Mollisquama sp.) covering chondrocranial, jaw, dental, and muscular anatomy, derived from a phase-contrast synchrotron microtomographic scan. Mollisquama sp. is unique among dalatiids in possessing a deep carinal process, extending ventrally from the bar between the subethmoid region and basal angle in squaloid sharks, containing a large fenestra infiltrated by the suborbitalis muscle. Mollisquama sp. also exhibits additional possibly diagnostic features, including a planar configuration of the labial cartilages and the absence of labial folds; a pad-like orbital process on the palatoquadrate; and the origination of the suborbitalis muscle solely on the carina, rather than the intraorbital wall. Character optimization of anatomical data onto a phylogeny of dalatiid sharks suggests Mollisquama sp. to be among the most specialized in the family, expanding the existing dalatiid morphospace. However, the functional significance of such transformations remains unclear. Synchrotron-derived data, which do not require chemical pretreatment of specimens, may elucidate soft-tissue functional correlates in future studies of undersampled taxa, such as dalatiids.

DNA capture reveals transoceanic gene flow in endangered river sharks.

For over a hundred years, the "river sharks" of the genus Glyphis were only known from the type specimens of species that had been collected in the 19th century. They were widely considered extinct until populations of Glyphis-like sharks were rediscovered in remote regions of Borneo and Northern Australia at the end of the 20th century. However, the genetic affinities between the newly discovered Glyphis-like populations and the poorly preserved, original museum-type specimens have never been established. Here, we present the first (to our knowledge) fully resolved, complete phylogeny of Glyphis that includes both archival-type specimens and modern material. We used a sensitive DNA hybridization capture method to obtain complete mitochondrial genomes from all of our samples and show that three of the five described river shark species are probably conspecific and widely distributed in Southeast Asia. Furthermore we show that there has been recent gene flow between locations that are separated by large oceanic expanses. Our data strongly suggest marine dispersal in these species, overturning the widely held notion that river sharks are restricted to freshwater. It seems that species in the genus Glyphis are euryhaline with an ecology similar to the bull shark, in which adult individuals live in the ocean while the young grow up in river habitats with reduced predation pressure. Finally, we discovered a previously unidentified species within the genus Glyphis that is deeply divergent from all other lineages, underscoring the current lack of knowledge about the biodiversity and ecology of these mysterious sharks.

Historical introgression drives pervasive mitochondrial admixture between two species of pelagic sharks.

We use a genomic sampling of both nuclear and mitochondrial DNA markers to examine a pattern of genetic admixture between Carcharhinus galapagensis (Galapagos sharks) and Carcharhinus obscurus (dusky sharks), two well-known and closely related sharks that have been recognized as valid species for more than 100years. We describe widespread mitochondrial-nuclear discordance in which these species are readily distinguishable based on 2152 nuclear single nucleotide polymorphisms from 910 independent autosomal regions, but show pervasive mitochondrial admixture. The species are superficially morphologically cryptic as adults but show marked differences in internal anatomy, as well as niche separation. There was no indication of ongoing hybridization between the species. We conclude that the observed mitochondrial-nuclear discordance is likely due to historical mitochondrial introgression following a range expansion.

Conservation genomics of natural and managed populations: building a conceptual and practical framework.

The boom of massive parallel sequencing (MPS) technology and its applications in conservation of natural and managed populations brings new opportunities and challenges to meet the scientific questions that can be addressed. Genomic conservation offers a wide range of approaches and analytical techniques, with their respective strengths and weaknesses that rely on several implicit assumptions. However, finding the most suitable approaches and analysis regarding our scientific question are often difficult and time-consuming. To address this gap, a recent workshop entitled 'ConGen 2015' was held at Montana University in order to bring together the knowledge accumulated in this field and to provide training in conceptual and practical aspects of data analysis applied to the field of conservation and evolutionary genomics. Here, we summarize the expertise yield by each instructor that has led us to consider the importance of keeping in mind the scientific question from sampling to management practices along with the selection of appropriate genomics tools and bioinformatics challenges.

Can We Identify Genes with Increased Phylogenetic Reliability?

Topological heterogeneity among gene trees is widely observed in phylogenomic analyses and some of this variation is likely caused by systematic error in gene tree estimation. Systematic error can be mitigated by improving models of sequence evolution to account for all evolutionary processes relevant to each gene or identifying those genes whose evolution best conforms to existing models. However, the best method for identifying such genes is not well established. Here, we ask if filtering genes according to their clock-likeness or posterior predictive effect size (PPES, an inference-based measure of model violation) improves phylogenetic reliability and congruence. We compared these approaches to each other, and to the common practice of filtering based on rate of evolution, using two different metrics. First, we compared gene-tree topologies to accepted reference topologies. Second, we examined topological similarity among gene trees in filtered sets. Our results suggest that filtering genes based on clock-likeness and PPES can yield a collection of genes with more reliable phylogenetic signal. For the two exemplar data sets we explored, from yeast and amniotes, clock-likeness and PPES outperformed rate-based filtering in both congruence and reliability.

Molecular phylogeny of Squaliformes and first occurrence of bioluminescence in sharks.

BACKGROUND: Squaliform sharks represent approximately 27 % of extant shark diversity, comprising more than 130 species with a predominantly deep-dwelling lifestyle. Many Squaliform species are highly specialized, including some that are bioluminescent, a character that is reported exclusively from Squaliform sharks within Chondrichthyes. The interfamiliar relationships within the order are still not satisfactorily resolved. Herein we estimate the phylogenetic interrelationships of a generic level sampling of "squaloid" sharks and closely related taxa using aligned sequences derived from a targeted gene capture approach. The resulting phylogenetic estimate is further used to evaluate the age of first occurrence of bioluminescence in Squaliformes. RESULTS: Our dataset comprised 172 putative ortholog exon sequences. Phylogenetic estimates result in a fully resolved tree supporting a monophyletic lineage of Squaliformes excluding Echinorhinus. Non-luminous Squalidae are inferred to be the sister to a clade comprising all remaining Squaliform families. Our results suggest that the origin of photophores is coincident with an elevated diversification rate and the splitting of families Dalatiidae, Etmopteridae, Oxynotidae and Somniosidae at the transition of the Lower to the Upper Cretaceous. The presence of luminous organs was confirmed for the Sleeper shark genus Zameus. These results indicate that bioluminescence in sharks is not restricted solely to the families Etmopteridae and Dalatiidae as previously believed. CONCLUSIONS: The sister-clade to non-luminous Squalidae comprises five families. The presence of photophores is reported for extant members of three out of these five families based on results of this study, i.e. Lantern sharks (Etmopteridae), Kitefin sharks (Dalatiidae) and Sleeper sharks (Somniosidae). Our results suggest that the origin of luminous organs arose during the rapid diversification event that gave rise to the extant Squaliform families. These inferences are consistent with the idea of diversification of Squaliform sharks being associated with the emergence of new deep-sea habitats in the Lower Cretaceous, which may have been facilitated by the evolution of bioluminescence.

Mitogenomics of the Speartooth Shark challenges ten years of control region sequencing.

BACKGROUND: Mitochondrial DNA markers have long been used to identify population boundaries and are now a standard tool in conservation biology. In elasmobranchs, evolutionary rates of mitochondrial genes are low and variation between distinct populations can be hard to detect with commonly used control region sequencing or other single gene approaches. In this study we sequenced the whole mitogenome of 93 Critically Endangered Speartooth Shark Glyphis glyphis from the last three river drainages they inhabit in northern Australia. RESULTS: Genetic diversity was extremely low (π =0.00019) but sufficient to demonstrate the existence of barriers to gene flow among river drainages (AMOVA Φ ST =0.28283, P <0.00001). Surprisingly, the comparison with single gene sub-datasets revealed that ND5 and 12S were the only ones carrying enough information to detect similar levels of genetic structure. The control region exhibited only one mutation, which was not sufficient to detect any structure among river drainages. CONCLUSIONS: This study strongly supports the use of single river drainages as discrete management units for the conservation of G. glyphis. Furthermore when genetic diversity is low, as is often the case in elasmobranchs, our results demonstrate a clear advantage of using the whole mitogenome to inform population structure compared to single gene approaches. More specifically, this study questions the extensive use of the control region as the preferential marker for elasmobranch population genetic studies and whole mitogenome sequencing will probably uncover a large amount of cryptic population structure in future studies.

Evolutionary trends in the elasmobranch neurocranium.

The neurocranium (braincase) is one of the defining vertebrate characters. Housing the brain and other key sensory organs, articulating with the jaws and contributing to the shape of the anteriormost portion of the body, the braincase is undoubtedly of great functional importance. Through studying relationships between braincase shape and ecology we can gain an improved understanding of form-function relationships in extant and fossil taxa. Elasmobranchii (sharks and rays) represent an important case study of vertebrate braincase diversity as their neurocranium is simplified and somewhat decoupled from other components of the cranium relative to other vertebrates. Little is known about the associations between ecology and braincase shape in this clade. In this study we report patterns of mosaic cranial evolution in Elasmobranchii that differ significantly from those present in other clades. The degree of evolutionary modularity also differs between Selachii and Batoidea. In both cases innovation in the jaw suspension appears to have driven shifts in patterns of integration and modularity, subsequently facilitating ecological diversification. Our results confirm the importance of water depth and biogeography as drivers of elasmobranch cranial diversity and indicate that skeletal articulation between the neurocranium and jaws represents a major constraint upon the evolution of braincase shape in vertebrates.

The appendicular skeleton of the enigmatic shark Leptocharias smithii in comparison with other sharks of the order Carcharhiniformes (Elasmobranchii: Leptochariidae).

Leptocharias smithii has been poorly explored in anatomical terms. This species bears a mosaic of morphological characters and is considered to represent an intermediate condition between other carcharhiniform clades. In the present paper, the anatomy of the appendicular skeleton of the species is thoroughly investigated and compared with other representatives of the order Carcharhiniformes. Leptocharias bears exclusive characteristics, such as the visible separation of the pro- and mesopterygia but it also has an aplesodic pectoral fin, a condition shared with carcharhiniforms placed at the base of the phylogenetic tree and at the same time a chevron-shaped coracoid bar, a condition characteristic of charcharhiniforms placed at the apex of the phylogenetic tree. Additionally, in an attempt to understand the evolution of its appendicular skeleton and of other carcharhiniforms, 20 characters of the paired fins and girdles are explored and discussed in light of two recent phylogenetic hypotheses. Most of these characters were not previously explored and support not only the monophyly of Carcharhiniformes, such as the mesopterygium overlapping the metapterygium in ventral view, but also the monophyly of the less inclusive clade Hemigaleidae + (Galeocerdonidae + (Carcharhinidae+Sphyrnidae)), such as the morphology and arrangement of the distal radials, which are pointed and spaced.