Complete reference genome assemblies and annotations of three Escherichia coli MRE162 clones

Escherichia coli MRE162 was originally isolated from a toilet pan in 1949 and since been utilized in numerous studies. Here, we sequence, assemble, and annotate clones held at three laboratories providing reference-level assemblies. We show the uniqueness of MRE162 to strains in open databases and make the UK clone publically available.

Reference genome assembly and annotation of two Bacillus cereus sensu lato strains from Etosha National Park, Namibia

Bacillus cereus sensu lato (s.l.) poses health and security issues. Here, we report the reference genome assembly of two Bacillus cereus s.l. strains, isolated from Etosha National Park, Namibia (FFI_BCgr36 and FFI_BCgr46). These unique genomes open for better understanding of environmental diversity and improvements in detection of threatening species.

Paleozoic origins of cheilostome bryozoans and their parental care inferred by a new genome-skimmed phylogeny

Phylogenetic relationships and the timing of evolutionary events are essential for understanding evolution on longer time scales. Cheilostome bryozoans are a group of ubiquitous, species-rich, marine colonial organisms with an excellent fossil record but lack phylogenetic relationships inferred from molecular data. We present genome-skimmed data for 395 cheilostomes and combine these with 315 published sequences to infer relationships and the timing of key events among c. 500 cheilostome species. We find that named cheilostome genera and species are phylogenetically coherent, rendering fossil or contemporary specimens readily delimited using only skeletal morphology. Our phylogeny shows that parental care in the form of brooding evolved several times independently but was never lost in cheilostomes. Our fossil calibration, robust to varied assumptions, indicates that the cheilostome lineage and parental care therein could have Paleozoic origins, much older than the first known fossil record of cheilostomes in the Late Jurassic.

A broadly resolved molecular phylogeny of New Zealand cheilostome bryozoans as a framework for hypotheses of morphological evolution

Larger molecular phylogenies based on ever more genes are becoming commonplace with the advent of cheaper and more streamlined sequencing and bioinformatics pipelines. However, many groups of inconspicuous but no less evolutionarily or ecologically important marine invertebrates are still neglected in the quest for understanding species- and higher-level phylogenetic relationships. Here, we alleviate this issue by presenting the molecular sequences of 165 cheilostome bryozoan species from New Zealand waters. New Zealand is our geographic region of choice as its cheilostome fauna is taxonomically, functionally and ecologically diverse, and better characterized than many other such faunas in the world. Using this most taxonomically broadly-sampled and statistically-supported cheilostome phylogeny comprising 214 species, when including previously published sequences, and 17 genes (2 nuclear and 15 mitochondrial) we tested several existing systematic hypotheses based solely on morphological observations. We find that lower taxonomic level hypotheses (species and genera) are robust while our inferred trees did not reflect current higher-level systematics (family and above), illustrating a general need for the rethinking of current hypotheses. To illustrate the utility of our new phylogeny, we reconstruct the evolutionary history of frontal shields (i.e., a calcified body-wall layer in ascus-bearing cheilostomes) and ask if its presence has any bearing on the diversification rates of cheilostomes.

A molecular phylogeny of historical and contemporary specimens of an under-studied micro-invertebrate group

Resolution of relationships at lower taxonomic levels is crucial for answering many evolutionary questions, and as such, sufficiently varied species representation is vital. This latter goal is not always achievable with relatively fresh samples. To alleviate the difficulties in procuring rarer taxa, we have seen increasing utilization of historical specimens in building molecular phylogenies using high throughput sequencing. This effort, however, has mainly focused on large-bodied or well-studied groups, with small-bodied and under-studied taxa under-prioritized. Here, we utilize both historical and contemporary specimens, to increase the resolution of phylogenetic relationships among a group of under-studied and small-bodied metazoans, namely, cheilostome bryozoans. In this study, we pioneer the sequencing of air-dried cheilostomes, utilizing a recently developed library preparation method for low DNA input. We evaluate a de novo mitogenome assembly and two iterative methods, using the sequenced target specimen as a reference for mapping, for our sequences. In doing so, we present mitochondrial and ribosomal RNA sequences of 43 cheilostomes representing 37 species, including 14 from historical samples ranging from 50 to 149 years old. The inferred phylogenetic relationships of these samples, analyzed together with publicly available sequence data, are shown in a statistically well-supported 65 taxa and 17 genes cheilostome tree, which is also the most broadly sampled and largest to date. The robust phylogenetic placement of historical samples whose contemporary conspecifics and/or congenerics have been sequenced verifies the appropriateness of our workflow and gives confidence in the phylogenetic placement of those historical samples for which there are no close relatives sequenced. The success of our workflow is highlighted by the circularization of a total of 27 mitogenomes, seven from historical cheilostome samples. Our study highlights the potential of utilizing DNA from micro-invertebrate specimens stored in natural history collections for resolving phylogenetic relationships among species.

The Genome of the Great Gerbil Reveals Species-Specific Duplication of an MHCII Gene

The great gerbil (Rhombomys opimus) is a social rodent living in permanent, complex burrow systems distributed throughout Central Asia, where it serves as the main host of several important vector-borne infectious pathogens including the well-known plague bacterium (Yersinia pestis). Here, we present a continuous annotated genome assembly of the great gerbil, covering over 96% of the estimated 2.47-Gb genome. Taking advantage of the recent genome assemblies of the sand rat (Psammomys obesus) and the Mongolian gerbil (Meriones unguiculatus), comparative immunogenomic analyses reveal shared gene losses within TLR gene families (i.e., TLR8, TLR10, and the entire TLR11-subfamily) for Gerbillinae, accompanied with signs of diversifying selection of TLR7 and TLR9. Most notably, we find a great gerbil-specific duplication of the MHCII DRB locus. In silico analyses suggest that the duplicated gene provides high peptide binding affinity for Yersiniae epitopes as well as Leishmania and Leptospira epitopes, putatively leading to increased capability to withstand infections by these pathogens. Our study demonstrates the power of whole-genome sequencing combined with comparative genomic analyses to gain deeper insight into the immunogenomic landscape of the great gerbil and its close relatives.

A genome-skimmed phylogeny of a widespread bryozoan family, Adeonidae

BACKGROUND

Understanding the phylogenetic relationships among species is one of the main goals of systematic biology. Simultaneously, credible phylogenetic hypotheses are often the first requirement for unveiling the evolutionary history of traits and for modelling macroevolutionary processes. However, many non-model taxa have not yet been sequenced to an extent such that statistically well-supported molecular phylogenies can be constructed for these purposes. Here, we use a genome-skimming approach to extract sequence information for 15 mitochondrial and 2 ribosomal operon genes from the cheilostome bryozoan family, Adeonidae, Busk, 1884, whose current systematics is based purely on morphological traits. The members of the Adeonidae are, like all cheilostome bryozoans, benthic, colonial, marine organisms. Adeonids are also geographically widely-distributed, often locally common, and are sometimes important habitat-builders.

RESULTS

We successfully genome-skimmed 35 adeonid colonies representing 6 genera (Adeona, Adeonellopsis, Bracebridgia, Adeonella, Laminopora and Cucullipora). We also contributed 16 new, circularised mitochondrial genomes to the eight previously published for cheilostome bryozoans. Using the aforementioned mitochondrial and ribosomal genes, we inferred the relationships among these 35 samples. Contrary to some previous suggestions, the Adeonidae is a robustly supported monophyletic clade. However, the genera Adeonella and Laminopora are in need of revision: Adeonella is polyphyletic and Laminopora paraphyletically forms a clade with some Adeonella species. Additionally, we assign a sequence clustering identity using cox1 barcoding region of 99% at the species and 83% at the genus level.

CONCLUSIONS

We provide sequence data, obtained via genome-skimming, that greatly increases the resolution of the phylogenetic relationships within the adeonids. We present a highly-supported topology based on 17 genes and substantially increase availability of circularised cheilostome mitochondrial genomes, and highlight how we can extend our pipeline to other bryozoans.

Enigmatic Diphyllatea eukaryotes: culturing and targeted PacBio RS amplicon sequencing reveals a higher order taxonomic diversity and global distribution

BACKGROUND

The class Diphyllatea belongs to a group of enigmatic unicellular eukaryotes that play a key role in reconstructing the morphological innovation and diversification of early eukaryotic evolution. Despite its evolutionary significance, very little is known about the phylogeny and species diversity of Diphyllatea. Only three species have described morphology, being taxonomically divided by flagella number, two or four, and cell size. Currently, one 18S rRNA Diphyllatea sequence is available, with environmental sequencing surveys reporting only a single partial sequence from a Diphyllatea-like organism. Accordingly, geographical distribution of Diphyllatea based on molecular data is limited, despite morphological data suggesting the class has a global distribution. We here present a first attempt to understand species distribution, diversity and higher order structure of Diphyllatea.

RESULTS

We cultured 11 new strains, characterised these morphologically and amplified their rRNA for a combined 18S-28S rRNA phylogeny. We sampled environmental DNA from multiple sites and designed new Diphyllatea-specific PCR primers for long-read PacBio RSII technology. Near full-length 18S rRNA sequences from environmental DNA, in addition to supplementary Diphyllatea sequence data mined from public databases, resolved the phylogeny into three deeply branching and distinct clades (Diphy I - III). Of these, the Diphy III clade is entirely novel, and in congruence with Diphy II, composed of species morphologically consistent with the earlier described Collodictyon triciliatum. The phylogenetic split between the Diphy I and Diphy II + III clades corresponds with a morphological division of Diphyllatea into bi- and quadriflagellate cell forms.

CONCLUSIONS

This altered flagella composition must have occurred early in the diversification of Diphyllatea and may represent one of the earliest known morphological transitions among eukaryotes. Further, the substantial increase in molecular data presented here confirms Diphyllatea has a global distribution, seemingly restricted to freshwater habitats. Altogether, the results reveal the advantage of combining a group-specific PCR approach and long-read high-throughput amplicon sequencing in surveying enigmatic eukaryote lineages. Lastly, our study shows the capacity of PacBio RS when targeting a protist class for increasing phylogenetic resolution.

Single Cell Transcriptomics, Mega-Phylogeny, and the Genetic Basis of Morphological Innovations in Rhizaria

The innovation of the eukaryote cytoskeleton enabled phagocytosis, intracellular transport, and cytokinesis, and is largely responsible for the diversity of morphologies among eukaryotes. Still, the relationship between phenotypic innovations in the cytoskeleton and their underlying genotype is poorly understood. To explore the genetic mechanism of morphological evolution of the eukaryotic cytoskeleton, we provide the first single cell transcriptomes from uncultured, free-living unicellular eukaryotes: the polycystine radiolarian Lithomelissa setosa (Nassellaria) and Sticholonche zanclea (Taxopodida). A phylogenomic approach using 255 genes finds Radiolaria and Foraminifera as separate monophyletic groups (together as Retaria), while Cercozoa is shown to be paraphyletic where Endomyxa is sister to Retaria. Analysis of the genetic components of the cytoskeleton and mapping of the evolution of these on the revised phylogeny of Rhizaria reveal lineage-specific gene duplications and neofunctionalization of α and β tubulin in Retaria, actin in Retaria and Endomyxa, and Arp2/3 complex genes in Chlorarachniophyta. We show how genetic innovations have shaped cytoskeletal structures in Rhizaria, and how single cell transcriptomics can be applied for resolving deep phylogenies and studying gene evolution in uncultured protist species.

Conservation of the abscission signaling peptide IDA during Angiosperm evolution: withstanding genome duplications and gain and loss of the receptors HAE/HSL2

The peptide INFLORESCENCE DEFICIENT IN ABSCISSION (IDA), which signals through the leucine-rich repeat receptor-like kinases HAESA (HAE) and HAESA-LIKE2 (HSL2), controls different cell separation events in Arabidopsis thaliana. We hypothesize the involvement of this signaling module in abscission processes in other plant species even though they may shed other organs than A. thaliana. As the first step toward testing this hypothesis from an evolutionarily perspective we have identified genes encoding putative orthologs of IDA and its receptors by BLAST searches of publically available protein, nucleotide and genome databases for angiosperms. Genes encoding IDA or IDA-LIKE (IDL) peptides and HSL proteins were found in all investigated species, which were selected as to represent each angiosperm order with available genomic sequences. The 12 amino acids representing the bioactive peptide in A. thaliana have virtually been unchanged throughout the evolution of the angiosperms; however, the number of IDL and HSL genes varies between different orders and species. The phylogenetic analyses suggest that IDA, HSL2, and the related HSL1 gene, were present in the species that gave rise to the angiosperms. HAE has arisen from HSL1 after a genome duplication that took place after the monocot-eudicots split. HSL1 has also independently been duplicated in the monocots, while HSL2 has been lost in gingers (Zingiberales) and grasses (Poales). IDA has been duplicated in eudicots to give rise to functionally divergent IDL peptides. We postulate that the high number of IDL homologs present in the core eudicots is a result of multiple whole genome duplications (WGD). We substantiate the involvement of IDA and HAE/HSL2 homologs in abscission by providing gene expression data of different organ separation events from various species.

Gene duplication, loss and selection in the evolution of saxitoxin biosynthesis in alveolates

A group of marine dinoflagellates (Alveolata, Eukaryota), consisting of ∼10 species of the genus Alexandrium, Gymnodinium catenatum and Pyrodinium bahamense, produce the toxin saxitoxin and its analogues (STX), which can accumulate in shellfish, leading to ecosystem and human health impacts. The genes, sxt, putatively involved in STX biosynthesis, have recently been identified, however, the evolution of these genes within dinoflagellates is not clear. There are two reasons for this: uncertainty over the phylogeny of dinoflagellates; and that the sxt genes of many species of Alexandrium and other dinoflagellate genera are not known. Here, we determined the phylogeny of STX-producing and other dinoflagellates based on a concatenated eight-gene alignment. We determined the presence, diversity and phylogeny of sxtA, domains A1 and A4 and sxtG in 52 strains of Alexandrium, and a further 43 species of dinoflagellates and thirteen other alveolates. We confirmed the presence and high sequence conservation of sxtA, domain A4, in 40 strains (35 Alexandrium, 1 Pyrodinium, 4 Gymnodinium) of 8 species of STX-producing dinoflagellates, and absence from non-producing species. We found three paralogs of sxtA, domain A1, and a widespread distribution of sxtA1 in non-STX producing dinoflagellates, indicating duplication events in the evolution of this gene. One paralog, clade 2, of sxtA1 may be particularly related to STX biosynthesis. Similarly, sxtG appears to be generally restricted to STX-producing species, while three amidinotransferase gene paralogs were found in dinoflagellates. We investigated the role of positive (diversifying) selection following duplication in sxtA1 and sxtG, and found negative selection in clades of sxtG and sxtA1, clade 2, suggesting they were functionally constrained. Significant episodic diversifying selection was found in some strains in clade 3 of sxtA1, a clade that may not be involved in STX biosynthesis, indicating pressure for diversification of function.

Alexandrium diversaporum sp. nov., a new non-saxitoxin producing species: Phylogeny, morphology and sxtA genes

Species of the PST producing planktonic marine dinoflagellate genus Alexandrium have been intensively scrutinised, and it is therefore surprising that new taxa can still be found. Here we report a new species, Alexandrium diversaporum nov. sp., isolated from spherical cysts found at two sites in Tasmania, Australia. This species differs in its morphology from all previously reported Alexandrium species, possessing a unique combination of morphological features: the presence of 2 size classes of thecal pores on the cell surface, a medium cell size, the size and shape of the 6″, 1', 2⁗ and Sp plates, the lack of a ventral pore, a lack of anterior and posterior connecting pores, and a lack of chain formation. We determined the relationship of the two strains to other species of Alexandrium based on an alignment of concatenated SSU-ITS1, 5.8S, ITS2 and partial LSU ribosomal RNA sequences, and found A. diversaporum to be a sister group to Alexandrium leei with high support. A. leei shares several morphological features, including the relative size and shapes of the 6″, 1', 2⁗ and Sp plates and the fact that some strains of A. leei have two size classes of thecal pores. We examined A. diversaporum strains for saxitoxin production and found them to be non-toxic. The species lacked sequences for the domain A4 of sxtA, as has been previously found for non-saxitoxin producing species of Alexandrium.

Evolution and distribution of saxitoxin biosynthesis in dinoflagellates

Numerous species of marine dinoflagellates synthesize the potent environmental neurotoxic alkaloid, saxitoxin, the agent of the human illness, paralytic shellfish poisoning. In addition, certain freshwater species of cyanobacteria also synthesize the same toxic compound, with the biosynthetic pathway and genes responsible being recently reported. Three theories have been postulated to explain the origin of saxitoxin in dinoflagellates: The production of saxitoxin by co-cultured bacteria rather than the dinoflagellates themselves, convergent evolution within both dinoflagellates and bacteria and horizontal gene transfer between dinoflagellates and bacteria. The discovery of cyanobacterial saxitoxin homologs in dinoflagellates has enabled us for the first time to evaluate these theories. Here, we review the distribution of saxitoxin within the dinoflagellates and our knowledge of its genetic basis to determine the likely evolutionary origins of this potent neurotoxin.

When naked became armored: an eight-gene phylogeny reveals monophyletic origin of theca in dinoflagellates

The dinoflagellates are a diverse lineage of microbial eukaryotes. Dinoflagellate monophyly and their position within the group Alveolata are well established. However, phylogenetic relationships between dinoflagellate orders remain unresolved. To date, only a limited number of dinoflagellate studies have used a broad taxon sample with more than two concatenated markers. This lack of resolution makes it difficult to determine the evolution of major phenotypic characters such as morphological features or toxin production e.g. saxitoxin. Here we present an improved dinoflagellate phylogeny, based on eight genes, with the broadest taxon sampling to date. Fifty-five sequences for eight phylogenetic markers from nuclear and mitochondrial regions were amplified from 13 species, four orders, and concatenated phylogenetic inferences were conducted with orthologous sequences. Phylogenetic resolution is increased with addition of support for the deepest branches, though can be improved yet further. We show for the first time that the characteristic dinoflagellate thecal plates, cellulosic material that is present within the sub-cuticular alveoli, appears to have had a single origin. In addition, the monophyly of most dinoflagellate orders is confirmed: the Dinophysiales, the Gonyaulacales, the Prorocentrales, the Suessiales, and the Syndiniales. Our improved phylogeny, along with results of PCR to detect the sxtA gene in various lineages, allows us to suggest that this gene was probably acquired separately in Gymnodinium and the common ancestor of Alexandrium and Pyrodinium and subsequently lost in some descendent species of Alexandrium.

Discovery of nuclear-encoded genes for the neurotoxin saxitoxin in dinoflagellates

Saxitoxin is a potent neurotoxin that occurs in aquatic environments worldwide. Ingestion of vector species can lead to paralytic shellfish poisoning, a severe human illness that may lead to paralysis and death. In freshwaters, the toxin is produced by prokaryotic cyanobacteria; in marine waters, it is associated with eukaryotic dinoflagellates. However, several studies suggest that saxitoxin is not produced by dinoflagellates themselves, but by co-cultured bacteria. Here, we show that genes required for saxitoxin synthesis are encoded in the nuclear genomes of dinoflagellates. We sequenced >1.2×10⁶ mRNA transcripts from the two saxitoxin-producing dinoflagellate strains Alexandrium fundyense CCMP1719 and A. minutum CCMP113 using high-throughput sequencing technology. In addition, we used in silico transcriptome analyses, RACE, qPCR and conventional PCR coupled with Sanger sequencing. These approaches successfully identified genes required for saxitoxin-synthesis in the two transcriptomes. We focused on sxtA, the unique starting gene of saxitoxin synthesis, and show that the dinoflagellate transcripts of sxtA have the same domain structure as the cyanobacterial sxtA genes. But, in contrast to the bacterial homologs, the dinoflagellate transcripts are monocistronic, have a higher GC content, occur in multiple copies, contain typical dinoflagellate spliced-leader sequences and eukaryotic polyA-tails. Further, we investigated 28 saxitoxin-producing and non-producing dinoflagellate strains from six different genera for the presence of genomic sxtA homologs. Our results show very good agreement between the presence of sxtA and saxitoxin-synthesis, except in three strains of A. tamarense, for which we amplified sxtA, but did not detect the toxin. Our work opens for possibilities to develop molecular tools to detect saxitoxin-producing dinoflagellates in the environment.

Biosynthesis and molecular genetics of polyketides in marine dinoflagellates

Marine dinoflagellates are the single most important group of algae that produce toxins, which have a global impact on human activities. The toxins are chemically diverse, and include macrolides, cyclic polyethers, spirolides and purine alkaloids. Whereas there is a multitude of studies describing the pharmacology of these toxins, there is limited or no knowledge regarding the biochemistry and molecular genetics involved in their biosynthesis. Recently, however, exciting advances have been made. Expressed sequence tag sequencing studies have revealed important insights into the transcriptomes of dinoflagellates, whereas other studies have implicated polyketide synthase genes in the biosynthesis of cyclic polyether toxins, and the molecular genetic basis for the biosynthesis of paralytic shellfish toxins has been elucidated in cyanobacteria. This review summarises the recent progress that has been made regarding the unusual genomes of dinoflagellates, the biosynthesis and molecular genetics of dinoflagellate toxins. In addition, the evolution of these metabolic pathways will be discussed, and an outlook for future research and possible applications is provided.

Evolutionary position of breviate amoebae and the primary eukaryote divergence

Integration of ultrastructural and molecular sequence data has revealed six supergroups of eukaryote organisms (excavates, Rhizaria, chromalveolates, Plantae, Amoebozoa and opisthokonts), and the root of the eukaryote evolutionary tree is suggested to lie between unikonts (Amoebozoa, opisthokonts) and bikonts (the other supergroups). However, some smaller lineages remain of uncertain affinity. One of these unassigned taxa is the anaerobic, free-living, amoeboid flagellate Breviata anathema, which is of key significance as it is unclear whether it is a unikont (i.e. possibly the deepest branching amoebozoan) or a bikont. To establish its evolutionary position, we sequenced thousands of Breviata genes and calculated trees using 78 protein sequences. Our trees and specific substitutions in the 18S RNA sequence indicate that Breviata is related to other Amoebozoa, thereby significantly increasing the cellular diversity of this phylum and establishing Breviata as a deep-branching unikont. We discuss the implications of these results for the ancestral state of Amoebozoa and eukaryotes generally, demonstrating that phylogenomics of phylogenetically 'nomadic' species can elucidate key questions in eukaryote evolution. Furthermore, mitochondrial genes among the Breviata ESTs demonstrate that Breviata probably contains a modified anaerobic mitochondrion. With these findings, remnants of mitochondria have been detected in all putatively deep-branching amitochondriate organisms.

Seven gene phylogeny of heterokonts

Nucleotide ssu and lsu rDNA sequences of all major lineages of autotrophic (Ochrophyta) and heterotrophic (Bigyra and Pseudofungi) heterokonts were combined with amino acid sequences from four protein-coding genes (actin, beta-tubulin, cox1 and hsp90) in a multigene approach for resolving the relationship between heterokont lineages. Applying these multigene data in Bayesian and maximum likelihood analyses improved the heterokont tree compared to previous rDNA analyses by placing all plastid-lacking heterotrophic heterokonts sister to Ochrophyta with robust support, and divided the heterotrophic heterokonts into the previously recognized phyla, Bigyra and Pseudofungi. Our trees identified the heterotrophic heterokonts Bicosoecida, Blastocystis and Labyrinthulida (Bigyra) as the earliest diverging lineages. A separate analysis of the phototrophic lineages, by adding the rbcL gene, further resolved the Ochrophyta lineages by increased support for several important nodes. Except for the positioning of Chrysophyceae, Eustigmatophyceae, Raphidophyceae and Pinguiophyceae, all main branches of Ochrophyta were resolved. Our results support the transfer of classes Dictyochophyceae and Pelagophyceae from subphylum Phaeista to Khakista. Based on all our trees, in combination with current knowledge about ultrastructure of heterokonts we suggest that a more advanced flagellar apparatus originated at one occasion in the ancestor of Phaeista whereas, Khakista independently reduced their flagellar apparatus and gained chlorophyll c(3).