The enigmatic mitochondrial genome of Rhabdopleura compacta (Pterobranchia) reveals insights into selection of an efficient tRNA system and supports monophyly of Ambulacraria

Zooid morphology and molecular phylogeny of the graptolite Rhabdopleura annulata (Hemichordata, Pterobranchia) from Heron Island, Australia

Rhabdopleura Allman, 1869 is one of the longest surviving animal genera. The five-known species are the only living Graptolithina, a group well known from their diverse Paleozoic fossil record. Here we add information on the soft-bodied zooids and molecular phylogenetics of Rhabdopleura annulata Norman, 1921, which was previously only known from its tubes. Tubes and zooids were collected from Heron Island, Queensland, Australia. Zooids have a single pair of tentaculated arms. Dark pigment granules are found throughout the body, and particularly dense in the pair of arms and the anterior lip of the cephalic shield. Colonies grow encrusted in and on coral debris. The tubes are either creeping or erect, but no stolon has been found. Inside of the coral matrix lacunae, the tube cortex formed a parchment-like wallpaper. Phylogenetic analysis based on combined 18S+16S rRNA sequences placed R. annulata as sister to the remaining rhabdopleurids, albeit with weak support. The biogeographic range of R. annulata extends from Indonesia to Tasmania, and New Zealand. Its occurrence on Heron Island does not extend this range, but highlights that rhabdopleurids may be more common, and in shallower waters, than previously appreciated, permitting further studies that may shed light on graptolite paleobiology.

New Insight Into Avian Papillomavirus Ecology and Evolution From Characterization of Novel Wild Bird Papillomaviruses

Viruses in the family Papillomaviridae have circular dsDNA genomes of approximately 5.7–8.6 kb that are packaged within non-enveloped, icosahedral capsids. The known papillomavirus (PV) representatives infect vertebrates, and there are currently more than 130 recognized PV species in more than 50 genera. We identified 12 novel avian papillomavirus (APV) types in wild birds that could represent five distinct species and two genera. Viruses were detected in paired oropharyngeal/cloacal swabs collected from six bird species, increasing the number of avian species known to harbor PVs by 40%. A new duck PV (DuPV-3) was found in mallard and American black duck (27.6% estimated prevalence) that was monophyletic with other known DuPVs. A single viral type was identified in Atlantic puffin (PuPV-1, 9.8% estimated prevalence), while a higher genetic diversity was found in other Charadriiformes. Specifically, three types [gull PV-1 (GuPV-1), -2, and -3] were identified in two gull species (estimated prevalence of 17% and 2.6% in American herring and great black-backed gull, respectively), and seven types [kittiwake PV-1 (KiPV-1) through -7] were found in black-legged kittiwake (81.3% estimated prevalence). Significantly higher DuPV-3 circulation was observed in spring compared to fall and in adults compared to juveniles. The studied host species’ tendencies to be in crowded environments likely affect infection rates and their migratory behaviors could explain the high viral diversity, illustrating how host behavior can influence viral ecology and distribution. For DuPV-3, GuPV-1, PuPV-1, and KiPV-2, we obtained the complete genomic sequences, which showed the same organization as other known APVs. Phylogenetic analyses showed evidence for virus–host co-divergence at the host taxonomic levels of family, order, and inter-order, but we also observed that host-specificity constraints are relaxed among highly related hosts as we found cross-species transmission within ducks and within gulls. Furthermore, the phylogeny of viruses infecting the Charadriiformes did not match the host phylogeny and gull viruses formed distinct monophyletic clades with kittiwake viruses, possibly reflecting past host-switching events. Considering the vast PV genotype diversity in other hosts and the large number of bird species, many more APVs likely remain to be discovered.

Mitogenomics Reveals a Novel Genetic Code in Hemichordata

The diverse array of codon reassignments demonstrate that the genetic code is not universal in nature. Exploring mechanisms underlying codon reassignment is critical for understanding the evolution of the genetic code during translation. Hemichordata, comprising worm-like Enteropneusta and colonial filter-feeding Pterobranchia, is the sister taxon of echinoderms and is more distantly related to chordates. However, only a few hemichordate mitochondrial genomes have been sequenced, hindering our understanding of mitochondrial genome evolution within Deuterostomia. In this study, we sequenced four mitochondrial genomes and two transcriptomes, including representatives of both major hemichordate lineages and analyzed together with public available data. Contrary to the current understanding of the mitochondrial genetic code in hemichordates, our comparative analyses suggest that UAA encodes Tyr instead of a "Stop" codon in the pterobranch lineage Cephalodiscidae. We also predict that AAA encodes Lys in pterobranch and enteropneust mitochondrial genomes, contradicting the previous assumption that hemichordates share the same genetic code with echinoderms for which AAA encodes Asn. Thus, we propose a new mitochondrial genetic code for Cephalodiscus and a revised code for enteropneusts. Moreover, our phylogenetic analyses are largely consistent with previous phylogenomic studies. The only exception is the phylogenetic position of the enteropneust Stereobalanus, whose placement as sister to all other described enteropneusts. With broader taxonomic sampling, we provide evidence that evolution of mitochondrial gene order and genetic codes in Hemichordata are more dynamic than previously thought and these findings provide insights into mitochondrial genome evolution within this clade.

Three new Diplozoidae mitogenomes expose unusual compositional biases within the Monogenea class: implications for phylogenetic studies

BACKGROUND

As the topologies produced by previous molecular and morphological studies were contradictory and unstable (polytomy), evolutionary relationships within the Diplozoidae family and the Monogenea class (controversial relationships among the Discocotylinea, Microcotylinea and Gastrocotylinea suborders) remain unresolved. Complete mitogenomes carry a relatively large amount of information, sufficient to provide a much higher phylogenetic resolution than traditionally used morphological traits and/or single molecular markers. However, their implementation is hampered by the scarcity of available monogenean mitogenomes. Therefore, we sequenced and characterized mitogenomes belonging to three Diplozoidae family species, and conducted comparative genomic and phylogenomic analyses for the entire Monogenea class.

RESULTS

Taxonomic identification was inconclusive, so two of the species were identified merely to the genus level. The complete mitogenomes of Sindiplozoon sp. and Eudiplozoon sp. are 14,334 bp and 15,239 bp in size, respectively. Paradiplozoon opsariichthydis (15,385 bp) is incomplete: an approximately 2000 bp-long gap within a non-coding region could not be sequenced. Each genome contains the standard 36 genes (atp8 is missing). G + T content and the degree of GC- and AT-skews of these three mitogenome (and their individual elements) were higher than in other monogeneans. nad2, atp6 and nad6 were the most variable PCGs, whereas cox1, nad1 and cytb were the most conserved. Mitochondrial phylogenomics analysis, conducted using concatenated amino acid sequences of all PCGs, indicates that evolutionary relationships of the three genera are: (Eudiplozoon, (Paradiplozoon, Sindiplozoon)); and of the three suborders: (Discocotylinea, (Microcotylinea, Gastrocotylinea)). These intergeneric relationships were also supported by the skewness and principal component analyses.

CONCLUSIONS

Our results show that nad2, atp6 and nad6 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in Diplozoidae. Nucleotide bias and codon and amino acid usage patterns of the three diplozoid mitogenomes are more similar to cestodes and trematodes than to other monogenean flatworms. This unusual mutational bias was reflected in disproportionately long branches in the phylogram. Our study offsets the scarcity of molecular data for the subclass Polyopisthocotylea to some extent, and might provide important new insights into the evolutionary history of the three genera and three suborders.

Hemichordate molecular phylogeny reveals a novel cold-water clade of harrimaniid acorn worms

Hemichordates are instrumental to understanding early deuterostome and chordate evolution, yet diversity and relationships within the group have been understudied. Recently, there has been renewed interest in hemichordate diversity and taxonomy, although current findings suggest that much hemichordate diversity remains to be discovered. Herein, we present a molecular phylogenetic study based on nuclear 18S rDNA sequence data, which includes 35 previously unsampled taxa and represents all recognized hemichordate families. We include mitochondrial 16S rDNA data from 66 enteropneust taxa and three pterobranch Rhabdopleura species, and recover colonial pterobranchs and solitary enteropneusts as reciprocally monophyletic taxa. Our phylogenetic results also reveal a previously unknown clade of at least four species of harrimaniid enteropneusts from cold waters, including Antarctica, the North Atlantic around Iceland and Norway, and the deep sea off Oregon. These small worms (1-5 mm in length), occur from 130 to 2950 m and are not closely related to other deep-sea harrimaniids, indicating that diversity of enteropneusts within the deep sea is broader than previously described in the literature. Discovery of this clade, as well as larger torquaratorids from Antarctica, strengthens hypotheses of close associations between Antarctic and deep-sea fauna.

The genome sequence of the colonial chordate, Botryllus schlosseri

Botryllus schlosseri is a colonial urochordate that follows the chordate plan of development following sexual reproduction, but invokes a stem cell-mediated budding program during subsequent rounds of asexual reproduction. As urochordates are considered to be the closest living invertebrate relatives of vertebrates, they are ideal subjects for whole genome sequence analyses. Using a novel method for high-throughput sequencing of eukaryotic genomes, we sequenced and assembled 580 Mbp of the B. schlosseri genome. The genome assembly is comprised of nearly 14,000 intron-containing predicted genes, and 13,500 intron-less predicted genes, 40% of which could be confidently parceled into 13 (of 16 haploid) chromosomes. A comparison of homologous genes between B. schlosseri and other diverse taxonomic groups revealed genomic events underlying the evolution of vertebrates and lymphoid-mediated immunity. The B. schlosseri genome is a community resource for studying alternative modes of reproduction, natural transplantation reactions, and stem cell-mediated regeneration. DOI:http://dx.doi.org/10.7554/eLife.00569.001.

Base-pairing versatility determines wobble sites in tRNA anticodons of vertebrate mitogenomes

BACKGROUND

Vertebrate mitochondrial genomes typically have one transfer RNA (tRNA) for each synonymous codon family. This limited anticodon repertoire implies that each tRNA anticodon needs to wobble (establish a non-Watson-Crick base pairing between two nucleotides in RNA molecules) to recognize one or more synonymous codons. Different hypotheses have been proposed to explain the factors that determine the nucleotide composition of wobble sites in vertebrate mitochondrial tRNA anticodons. Until now, the two major postulates--the "codon-anticodon adaptation hypothesis" and the "wobble versatility hypothesis"--have not been formally tested in vertebrate mitochondria because both make the same predictions regarding the composition of anticodon wobble sites. The same is true for the more recent "wobble cost hypothesis".

PRINCIPAL FINDINGS

In this study we have analyzed the occurrence of synonymous codons and tRNA anticodon wobble sites in 1553 complete vertebrate mitochondrial genomes, focusing on three fish species with mtDNA codon usage bias reversal (L-strand is GT-rich). These mitogenomes constitute an excellent opportunity to study the evolution of the wobble nucleotide composition of tRNA anticodons because due to the reversal the predictions for the anticodon wobble sites differ between the existing hypotheses. We observed that none of the wobble sites of tRNA anticodons in these unusual mitochondrial genomes coevolved to match the new overall codon usage bias, suggesting that nucleotides at the wobble sites of tRNA anticodons in vertebrate mitochondrial genomes are determined by wobble versatility.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that, at wobble sites of tRNA anticodons in vertebrate mitogenomes, selection favors the most versatile nucleotide in terms of wobble base-pairing stability and that wobble site composition is not influenced by codon usage. These results are in agreement with the "wobble versatility hypothesis".

The evolution of the mitochondrial genetic code in arthropods revisited

A variant of the invertebrate mitochondrial genetic code was previously identified in arthropods (Abascal et al. 2006a, PLoS Biol 4:e127) in which, instead of translating the AGG codon as serine, as in other invertebrates, some arthropods translate AGG as lysine. Here, we revisit the evolution of the genetic code in arthropods taking into account that (1) the number of arthropod mitochondrial genomes sequenced has triplicated since the original findings were published; (2) the phylogeny of arthropods has been recently resolved with confidence for many groups; and (3) sophisticated probabilistic methods can be applied to analyze the evolution of the genetic code in arthropod mitochondria. According to our analyses, evolutionary shifts in the genetic code have been more common than previously inferred, with many taxonomic groups displaying two alternative codes. Ancestral character-state reconstruction using probabilistic methods confirmed that the arthropod ancestor most likely translated AGG as lysine. Point mutations at tRNA-Lys and tRNA-Ser correlated with the meaning of the AGG codon. In addition, we identified three variables (GC content, number of AGG codons, and taxonomic information) that best explain the use of each of the two alternative genetic codes.