The revised classification of eukaryotes

Dynamic species classification of microorganisms across time, abiotic and biotic environments-A sliding window approach

The development of video-based monitoring methods allows for rapid, dynamic and accurate monitoring of individuals or communities, compared to slower traditional methods, with far reaching ecological and evolutionary applications. Large amounts of data are generated using video-based methods, which can be effectively processed using machine learning (ML) algorithms into meaningful ecological information. ML uses user defined classes (e.g. species), derived from a subset (i.e. training data) of video-observed quantitative features (e.g. phenotypic variation), to infer classes in subsequent observations. However, phenotypic variation often changes due to environmental conditions, which may lead to poor classification, if environmentally induced variation in phenotypes is not accounted for. Here we describe a framework for classifying species under changing environmental conditions based on the random forest classification. A sliding window approach was developed that restricts temporal and environmentally conditions to improve the classification. We tested our approach by applying the classification framework to experimental data. The experiment used a set of six ciliate species to monitor changes in community structure and behavior over hundreds of generations, in dozens of species combinations and across a temperature gradient. Differences in biotic and abiotic conditions caused simplistic classification approaches to be unsuccessful. In contrast, the sliding window approach allowed classification to be highly successful, as phenotypic differences driven by environmental change, could be captured by the classifier. Importantly, classification using the random forest algorithm showed comparable success when validated against traditional, slower, manual identification. Our framework allows for reliable classification in dynamic environments, and may help to improve strategies for long-term monitoring of species in changing environments. Our classification pipeline can be applied in fields assessing species community dynamics, such as eco-toxicology, ecology and evolutionary ecology.

Comparative interactomics provides evidence for functional specialization of the nuclear pore complex

The core architecture of the eukaryotic cell was established well over one billion years ago, and is largely retained in all extant lineages. However, eukaryotic cells also possess lineage-specific features, frequently keyed to specific functional requirements. One quintessential core eukaryotic structure is the nuclear pore complex (NPC), responsible for regulating exchange of macromolecules between the nucleus and cytoplasm as well as acting as a nuclear organizational hub. NPC architecture has been best documented in one eukaryotic supergroup, the Opisthokonts (e.g. Saccharomyces cerevisiae and Homo sapiens), which although compositionally similar, have significant variations in certain NPC subcomplex structures. The variation of NPC structure across other taxa in the eukaryotic kingdom however, remains poorly understood. We explored trypanosomes, highly divergent organisms, and mapped and assigned their NPC proteins to specific substructures to reveal their NPC architecture. We showed that the NPC central structural scaffold is conserved, likely across all eukaryotes, but more peripheral elements can exhibit very significant lineage-specific losses, duplications or other alterations in their components. Amazingly, trypanosomes lack the major components of the mRNA export platform that are asymmetrically localized within yeast and vertebrate NPCs. Concomitant with this, the trypanosome NPC is ALMOST completely symmetric with the nuclear basket being the only major source of asymmetry. We suggest these features point toward a stepwise evolution of the NPC in which a coating scaffold first stabilized the pore after which selective gating emerged and expanded, leading to the addition of peripheral remodeling machineries on the nucleoplasmic and cytoplasmic sides of the pore.

How Embryophytic is the Biosynthesis of Phenylpropanoids and their Derivatives in Streptophyte Algae?

The origin of land plants from algae is a long-standing question in evolutionary biology. It is becoming increasingly clear that many characters that were once assumed to be 'embryophyte specific' can in fact be found in their closest algal relatives, the streptophyte algae. One such case is the phenylpropanoid pathway. While biochemical data indicate that streptophyte algae harbor lignin-like components, the phenylpropanoid core pathway, which serves as the backbone of lignin biosynthesis, has been proposed to have arisen at the base of the land plants. Here we revisit this hypothesis using a wealth of new sequence data from streptophyte algae. Tracing the biochemical pathway towards lignin biogenesis, we show that most of the genes required for phenylpropanoid synthesis and the precursors for lignin production were already present in streptophyte algae. Nevertheless, phylogenetic analyses and protein structure predictions of one of the key enzyme classes in lignin production, cinnamyl alcohol dehydrogenase (CAD), suggest that CADs of streptophyte algae are more similar to sinapyl alcohol dehydrogenases (SADs). This suggests that the end-products of the pathway leading to lignin biosynthesis in streptophyte algae may facilitate the production of lignin-like compounds and defense molecules. We hypothesize that streptophyte algae already possessed the genetic toolkit from which the capacity to produce lignin later evolved in vascular plants.

UniEuk: Time to Speak a Common Language in Protistology!

Universal taxonomic frameworks have been critical tools to structure the fields of botany, zoology, mycology, and bacteriology as well as their large research communities. Animals, plants, and fungi have relatively solid, stable morpho‐taxonomies built over the last three centuries, while bacteria have been classified for the last three decades under a coherent molecular taxonomic framework. By contrast, no such common language exists for microbial eukaryotes, even though environmental ‘‐omics’ surveys suggest that protists make up most of the organismal and genetic complexity of our planet's ecosystems! With the current deluge of eukaryotic meta‐omics data, we urgently need to build up a universal eukaryotic taxonomy bridging the protist ‐omics age to the fragile, centuries‐old body of classical knowledge that has effectively linked protist taxa to morphological, physiological, and ecological information. UniEuk is an open, inclusive, community‐based and expert‐driven international initiative to build a flexible, adaptive universal taxonomic framework for eukaryotes. It unites three complementary modules, EukRef, EukBank, and EukMap, which use phylogenetic markers, environmental metabarcoding surveys, and expert knowledge to inform the taxonomic framework. The UniEuk taxonomy is directly implemented in the European Nucleotide Archive at EMBL‐EBI, ensuring its broad use and long‐term preservation as a reference taxonomy for eukaryotes.

The Protistan Microbiome of Grassland Soil: Diversity in the Mesoscale

Genomic data for less than one quarter of ∼1.8 million named species on earth exist in public databases like GenBank. Little information exists on the estimated one million small sized (1-100μm) heterotrophic nanoflagellates and ciliates and their taxa-area relationship. We analyzed environmental DNA from 150 geo-referenced grassland plots representing topographical and land-use ranges typical for Central Europe. High through-put barcoding allowed the identification of operational taxonomic units (OTUs) at species level, with high pairwise identity to reference sequences (≥99.7%), but also the identification of sequences at the genus (≥97%) and class (≥80%) taxonomic level. Species richness analyses revealed, on average, 100 genus level OTUs (332 unique individual read (UIR) and 56 class level OTUs per gram of soil sample in the mesoscale (1-1000km). Database shortfalls were highlighted by increased uncertain taxonomic lineages at lower resolution (≥80% sequence identity). No single barcode occurred ubiquitously across all sites. Taxa-area relationships indicated that OTUs spread over the entire mesoscale were more similar than in the local scale and increased land-use (fertilization, mowing and grazing) promoted taxa-area separation. Only a small fraction of sequences strictly matched reference library sequences, suggesting a large protistan "dark matter" in soil which warrants further research.

Origin and diversity of testate amoebae shell composition: Example of Bullinularia indica living in Sphagnum capillifolium

Testate amoebae are free-living shelled protists that build a wide range of shells with various sizes, shapes, and compositions. Recent studies showed that xenosomic testate amoebae shells could be indicators of atmospheric particulate matter (PM) deposition. However, no study has yet been conducted to assess the intra-specific mineral, organic, and biologic grain diversity of a single xenosomic species in a natural undisturbed environment. This study aims at providing new information about grain selection to develop the potential use of xenosomic testate amoebae shells as bioindicators of the multiple-origin mineral/organic diversity of their proximal environment. To fulfil these objectives, we analysed the shell content of 38 Bullinularia indica individuals, a single xenosomic testate amoeba species living in Sphagnum capillifolium, by scanning electron microscope (SEM) coupled with X-ray spectroscopy. The shells exhibited high diversities of mineral, organic, and biomineral grains, which confirms their capability to recycle xenosomes. Mineral grain diversity and size of B. indica matched those of the atmospheric natural mineral PM deposited in the peatbog. Calculation of grain size sorting revealed a discrete selection of grains agglutinated by B. indica. These results are a first step towards understanding the mechanisms of particle selection by xenosomic testate amoebae in natural conditions.

Well-positioned nucleosomes punctuate polycistronic pol II transcription units and flank silent VSG gene arrays in Trypanosoma brucei

Background

The compaction of DNA in chromatin in eukaryotes allowed the expansion of genome size and coincided with significant evolutionary diversification. However, chromatin generally represses DNA function, and mechanisms coevolved to regulate chromatin structure and its impact on DNA. This included the selection of specific nucleosome positions to modulate accessibility to the DNA molecule. Trypanosoma brucei, a member of the Excavates supergroup, falls in an ancient evolutionary branch of eukaryotes and provides valuable insight into the organization of chromatin in early genomes.

Results

We have mapped nucleosome positions in T. brucei and identified important differences compared to other eukaryotes: The RNA polymerase II initiation regions in T. brucei do not exhibit pronounced nucleosome depletion, and show little evidence for defined −1 and +1 nucleosomes. In contrast, a well-positioned nucleosome is present directly on the splice acceptor sites within the polycistronic transcription units. The RNA polyadenylation sites were depleted of nucleosomes, with a single well-positioned nucleosome present immediately downstream of the predicted sites. The regions flanking the silent variant surface glycoprotein (VSG) gene cassettes showed extensive arrays of well-positioned nucleosomes, which may repress cryptic transcription initiation. The silent VSG genes themselves exhibited a less regular nucleosomal pattern in both bloodstream and procyclic form trypanosomes. The DNA replication origins, when present within silent VSG gene cassettes, displayed a defined nucleosomal organization compared with replication origins in other chromosomal core regions.

Conclusions

Our results indicate that some organizational features of chromatin are evolutionarily ancient, and may already have been present in the last eukaryotic common ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0121-9) contains supplementary material, which is available to authorized users.

Nuclear codon reassignments in the genomics era and mechanisms behind their evolution

The canonical genetic code ubiquitously translates nucleotide into peptide sequence with several alterations known in viruses, bacteria, mitochondria, plastids, and single-celled eukaryotes. A new hypothesis to explain genetic code changes, termed tRNA loss driven codon reassignment, has been proposed recently when the polyphyly of the yeast codon reassignment events has been uncovered. According to this hypothesis, the driving force for genetic code changes are tRNA or translation termination factor loss-of-function mutations or loss-of-gene events. The free codon can subsequently be captured by all tRNAs that have an appropriately mutated anticodon and are efficiently charged. Thus, codon capture most likely happens by near-cognate tRNAs and tRNAs whose anticodons are not part of the recognition sites of the respective aminoacyl-tRNA-synthetases. This hypothesis comprehensively explains the CTG codon translation as alanine in Pachysolen yeast together with the long known translation of the same codon as serine in Candida albicans and related species, and can also be applied to most other known reassignments.

SILVA, RDP, Greengenes, NCBI and OTT - how do these taxonomies compare?

Background

A key step in microbiome sequencing analysis is read assignment to taxonomic units. This is often performed using one of four taxonomic classifications, namely SILVA, RDP, Greengenes or NCBI. It is unclear how similar these are and how to compare analysis results that are based on different taxonomies.

Results

We provide a method and software for mapping taxonomic entities from one taxonomy onto another. We use it to compare the four taxonomies and the Open Tree of life Taxonomy (OTT).

Conclusions

While we find that SILVA, RDP and Greengenes map well into NCBI, and all four map well into the OTT, mapping the two larger taxonomies on to the smaller ones is problematic.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3501-4) contains supplementary material, which is available to authorized users.

Characterization of Cytosine Methylation and the DNA Methyltransferases of Toxoplasma gondii

DNA methylation is a key epigenetic modification which confers phenotypic plasticity and adaptation. Cyst-forming strains of Toxoplasma gondii undergo tachyzoite to bradyzoite conversion after initial acute infection of a host, and the reverse conversion may occur in immune-suppressed hosts. The formation of m⁵C is catalyzed by DNA methyltransferase (DNMT). We identified two functional DNA methyltransferases, TgDNMTa and TgDNMTb, in T. gondii that may mediate DNA methylation. The recombinant proteins showed intrinsic methyltransferase activity; both have higher transcription levels in bradyzoites than that in tachyzoites. We performed genome-wide analysis of DNA methylation in tachyzoites and bradyzoites. The results showed more methylation sites in bradyzoites than that in tachyzoites. The most significantly enriched GO-terms of genes with DNA methylation were associated with basal cellular processes such as energy metabolism and parasite resistance to host immunity. Tachyzoite proliferation in parasitophorous vacuoles (PV) can be inhibited by the DNA methyltransferase inhibitor 5-azacytidine, a chemical analogue of the nucleotide cytosine that can inactivate DNA methyltransferases. These findings provide the first confirmation of DNA methylation in T. gondii.

Analysis of triglyceride synthesis unveils a green algal soluble diacylglycerol acyltransferase and provides clues to potential enzymatic components of the chloroplast pathway

Background

Microalgal triglyceride (TAG) synthesis has attracted considerable attention. Particular emphasis has been put towards characterizing the algal homologs of the canonical rate-limiting enzymes, diacylglycerol acyltransferase (DGAT) and phospholipid:diacylglycerol acyltransferase (PDAT). Less work has been done to analyze homologs from a phylogenetic perspective. In this work, we used HMMER iterative profiling and phylogenetic and functional analyses to determine the number and sequence characteristics of algal DGAT and PDAT, as well as related sequences that constitute their corresponding superfamilies. We included most algae with available genomes, as well as representative eukaryotic and prokaryotic species.

Results

Amongst our main findings, we identified a novel clade of DGAT1-like proteins exclusive to red algae and glaucophyta and a previously uncharacterized subclade of DGAT2 proteins with an unusual number of transmembrane segments. Our analysis also revealed the existence of a novel DGAT exclusive to green algae with moderate similarity to plant soluble DGAT3. The DGAT3 clade shares a most recent ancestor with a group of uncharacterized proteins from cyanobacteria. Subcellular targeting prediction suggests that most green algal DGAT3 proteins are imported to the chloroplast, evidencing that the green algal chloroplast might have a soluble pathway for the de novo synthesis of TAGs. Heterologous expression of C. reinhardtii DGAT3 produces an increase in the accumulation of TAG, as evidenced by thin layer chromatography.

Conclusions

Our analysis contributes to advance in the knowledge of complex superfamilies involved in lipid metabolism and provides clues to possible enzymatic players of chloroplast TAG synthesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3602-0) contains supplementary material, which is available to authorized users.

Discovering Protein-Coding Genes from the Environment: Time for the Eukaryotes?

Eukaryotic microorganisms from diverse environments encompass a large number of taxa, many of them still unknown to science. One strategy to mine these organisms for genes of biotechnological relevance is to use a pool of eukaryotic mRNA directly extracted from environmental samples. Recent reports demonstrate that the resulting metatranscriptomic cDNA libraries can be screened by expression in yeast for a wide range of genes and functions from many of the different eukaryotic taxa. In combination with novel emerging high-throughput technologies, we anticipate that this approach should contribute to exploring the functional diversity of the eukaryotic microbiota.

Structures related to attachment and motility in the marine eugregarine Cephaloidophora cf. communis (Apicomplexa)

Gregarines represent a highly diversified group of ancestral apicomplexans, with various modes of locomotion and host-parasite interactions. The eugregarine parasite of the barnacle Balanus balanus, Cephaloidophora cf. communis, exhibits interesting organisation of its attachment apparatus along with unique motility modes. The pellicle covered gregarine is arranged into longitudinal epicytic folds. The epimerite is separated from the protomerite by a septum consisting of tubulin-rich filamentous structures and both are packed with microneme-like structures suggestive of their function in the production of adhesives important for attachment and secreted through the abundant epimerite pores. Detached trophozoites and gamonts are capable of gliding motility, enriched by jumping and rotational movements with rapid changes in gliding direction and cell flexions. Actin in its polymerised form (F-actin) is distributed throughout the entire gregarine, while myosin, detected in the cortical region of the cell, follows the pattern of the epicytic folds. Various motility modes exhibited by individuals of C. cf. communis, together with significant changes in their cell shape during locomotion, are not concordant with the gliding mechanisms generally described in apicomplexan zoites and indicate that additional structures must be involved (e.g. two 12-nm filaments; the specific dentate appearance of internal lamina inside the epicytic folds).

A trophic framework for animal origins

Metazoans emerged in a microbial world and play a unique role in the biosphere as the only complex multicellular eukaryotes capable of phagocytosis. While the bodyplan and feeding mode of the last common metazoan ancestor remain unresolved, the earliest multicellular stem-metazoans likely subsisted on picoplankton (planktonic microbes 0.2-2 μm in diameter) and dissolved organic matter (DOM), similarly to modern sponges. Once multicellular stem-metazoans emerged, they conceivably modulated both the local availability of picoplankton, which they preferentially removed from the water column for feeding, and detrital particles 2-100 μm in diameter, which they expelled and deposited into the benthos as waste products. By influencing the availability of these heterotrophic food sources, the earliest multicellular stem-metazoans would have acted as ecosystem engineers, helping create the ecological conditions under which other metazoans, namely detritivores and non-sponge suspension feeders incapable of subsisting on picoplankton and DOM, could emerge and diversify. This early style of metazoan feeding, specifically the phagocytosis of small eukaryotic prey, could have also encouraged the evolution of larger, even multicellular, eukaryotic forms less prone to metazoan consumption. Therefore, the first multicellular stem-metazoans, through their feeding, arguably helped bridge the strictly microbial food webs of the Proterozoic Eon (2.5-0.541 billion years ago) to the more macroscopic, metazoan-sustaining food webs of the Phanerozoic Eon (0.541-0 billion years ago).

Key Ecological Roles for Zoosporic True Fungi in Aquatic Habitats

The diversity and abundance of zoosporic true fungi have been analyzed recently using fungal sequence libraries and advances in molecular methods, such as high-throughput sequencing. This review focuses on four evolutionary primitive true fungal phyla: the Aphelidea, Chytridiomycota, Neocallimastigomycota, and Rosellida (Cryptomycota), most species of which are not polycentric or mycelial (filamentous), rather they tend to be primarily monocentric (unicellular). Zoosporic fungi appear to be both abundant and diverse in many aquatic habitats around the world, with abundance often exceeding other fungal phyla in these habitats, and numerous novel genetic sequences identified. Zoosporic fungi are able to survive extreme conditions, such as high and extremely low pH; however, more work remains to be done. They appear to have important ecological roles as saprobes in decomposition of particulate organic substrates, pollen, plant litter, and dead animals; as parasites of zooplankton and algae; as parasites of vertebrate animals (such as frogs); and as symbionts in the digestive tracts of mammals. Some chytrids cause economically important diseases of plants and animals. They regulate sizes of phytoplankton populations. Further metagenomics surveys of aquatic ecosystems are expected to enlarge our knowledge of the diversity of true zoosporic fungi. Coupled with studies on their functional ecology, we are moving closer to unraveling the role of zoosporic fungi in carbon cycling and the impact of climate change on zoosporic fungal populations.

Parasitism as the main factor shaping peptide vocabularies in current organisms

Self/non-self-discrimination by vertebrate immune systems is based on the recognition of the presence of peptides in proteins of a parasite that are not contained in the proteins of a host. Therefore, a reduction of the number of 'words' in its own peptide vocabulary could be an efficient evolutionary strategy of parasites for escaping recognition. Here, we compared peptide vocabularies of 30 endoparasitic and 17 free-living unicellular organisms and also eight multicellular parasitic and 16 multicellular free-living organisms. We found that both unicellular and multicellular parasites used a significantly lower number of different pentapeptides than free-living controls. Impoverished pentapeptide vocabularies in parasites were observed across all five clades that contain both the parasitic and free-living species. The effect of parasitism on a number of peptides used in an organism's proteins is larger than effects of all other studied factors, including the size of a proteome, the number of encoded proteins, etc. This decrease of pentapeptide diversity was partly compensated for by an increased number of hexapeptides. Our results support the hypothesis of parasitism-associated reduction of peptide vocabulary and suggest that T-cell receptors mostly recognize the five amino acids-long part of peptides that are presented in the groove of major histocompatibility complex molecules.

Loss of GET pathway orthologs in Arabidopsis thaliana causes root hair growth defects and affects SNARE abundance

Soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins are key players in cellular trafficking and coordinate vital cellular processes, such as cytokinesis, pathogen defense, and ion transport regulation. With few exceptions, SNAREs are tail-anchored (TA) proteins, bearing a C-terminal hydrophobic domain that is essential for their membrane integration. Recently, the Guided Entry of Tail-anchored proteins (GET) pathway was described in mammalian and yeast cells that serve as a blueprint of TA protein insertion [Schuldiner M, et al. (2008) Cell 134(4):634-645; Stefanovic S, Hegde RS (2007) Cell 128(6):1147-1159]. This pathway consists of six proteins, with the cytosolic ATPase GET3 chaperoning the newly synthesized TA protein posttranslationally from the ribosome to the endoplasmic reticulum (ER) membrane. Structural and biochemical insights confirmed the potential of pathway components to facilitate membrane insertion, but the physiological significance in multicellular organisms remains to be resolved. Our phylogenetic analysis of 37 GET3 orthologs from 18 different species revealed the presence of two different GET3 clades. We identified and analyzed GET pathway components in Arabidopsis thaliana and found reduced root hair elongation in Atget lines, possibly as a result of reduced SNARE biogenesis. Overexpression of AtGET3a in a receptor knockout (KO) results in severe growth defects, suggesting presence of alternative insertion pathways while highlighting an intricate involvement for the GET pathway in cellular homeostasis of plants.

Deep Transcriptome Sequencing of Two Green Algae, Chara vulgaris and Chlamydomonas reinhardtii, Provides No Evidence of Organellar RNA Editing

Nearly all land plants post-transcriptionally modify specific nucleotides within RNAs, a process known as RNA editing. This adaptation allows the correction of deleterious mutations within the asexually reproducing and presumably non-recombinant chloroplast and mitochondrial genomes. There are no reports of RNA editing in any of the green algae so this phenomenon is presumed to have originated in embryophytes either after the invasion of land or in the now extinct algal ancestor of all land plants. This was challenged when a recent in silico screen for RNA edit sites based on genomic sequence homology predicted edit sites in the green alga Chara vulgaris, a multicellular alga found within the Streptophyta clade and one of the closest extant algal relatives of land plants. In this study, the organelle transcriptomes of C. vulgaris and Chlamydomonas reinhardtii were deep sequenced for a comprehensive assessment of RNA editing. Initial analyses based solely on sequence comparisons suggested potential edit sites in both species, but subsequent high-resolution melt analysis, RNase H-dependent PCR (rhPCR), and Sanger sequencing of DNA and complementary DNAs (cDNAs) from each of the putative edit sites revealed them to be either single-nucleotide polymorphisms (SNPs) or spurious deep sequencing results. The lack of RNA editing in these two lineages is consistent with the current hypothesis that RNA editing evolved after embryophytes split from its ancestral algal lineage.

Diversity of the Photosynthetic Paulinella Species, with the Description of Paulinella micropora sp. nov. and the Chromatophore Genome Sequence for strain KR01

The thecate filose amoeba Paulinella chromatophora is a good model organism for understanding plastid organellogenesis because its chromatophore was newly derived from an alpha-cyanobacterium. Paulinella chromatophora was the only known photosynthetic Paulinella species until recent studies that suggested a species level of diversity. Here, we described a new photosynthetic species P. micropora sp. nov. based on morphological and molecular evidence from a newly established strain KR01. The chromatophore genome of P. micropora KR01 was fully determined; the genome was 976,991bp in length, the GC content was 39.9%, and 908 genes were annotated. A pairwise comparison of chromatophore genome sequences between strains KR01 and FK01, representing two different natural populations of P. micropora, showed a 99.85% similarity. Differences between the two strains included single nucleotide polymorphisms (SNPs) in CDSs, which resulted in 357 synonymous and 280 nonsynonymous changes, along with 245 SNPs in non-coding regions. Indels (37) and microinversions (14) were also detected. Species diversity for photosynthetic Paulinella was surveyed using samples collected from around the world. We compared our new species to two photosynthetic species, P. chromatophora and P. longichromatophora. Phylogenetic analyses using four gene markers revealed three distinct lineages of photosynthetic Paulinella species including P. micropora sp. nov.

Lineage-Specific and Highly Derived Gene Sequences Among Amoebozoa, Revealed by the Comparative Analysis of Transcriptomes from Twelve Amoebozoan Species

Amoebozoa represent a difficult group for traditional morphology-based taxonomy. Molecular approaches, such as gene sequencing and DNA barcoding have greatly enhanced our knowledge of the diversity of these organisms. However, metagenomic studies of Amoebozoa still did not provide as impressive results as they did among some other groups of protists. In environmental DNA surveys done on fragments of SSU rDNA gene and other traditional DNA barcodes, Amoebozoa genes normally constitute a minor part of the total gene diversity and represent only the most abundant lineages. A potential way to resolve this problem is the usage of DNA barcodes based on genes, which are unique or highly derived in this group of organisms. In the present study, we attempted to find such genes and gene families with a low level of paralogy, potentially appropriate as Amoebozoa-specific DNA barcodes. For this we re-assembled transcriptomes of 12 amoebozoan species available from the public databases and performed gene annotation and identification of orthologous genes. In our analysis Amoebozoa-specific and highly derived sequences formed 53,182 clusters of orthologs, containing from 2 to 299 proteins each. Some of these genes may be a potential target for DNA barcoding of Amoebozoa.

DNA metabarcoding reveals that 200-μm-size-fractionated filtering is unable to discriminate between planktonic microbial and large eukaryotes

Microeukaryotic plankton (0.2-200 μm) are critical components of aquatic ecosystems and key players in global ecological processes. High-throughput sequencing is currently revolutionizing their study on an unprecedented scale. However, it is currently unclear whether we can accurately, effectively and quantitatively depict the microeukaryotic plankton communities using traditional size-fractionated filtering combined with molecular methods. To address this, we analysed the eukaryotic plankton communities both with, and without, prefiltering with a 200 μm pore-size sieve -by using SSU rDNA-based high-throughput sequencing on 16 samples with three replicates in each sample from two subtropical reservoirs sampled from January to October in 2013. We found that ~25% reads were classified as metazoan in both size groups. The species richness, alpha and beta diversity of plankton community and relative abundance of reads in 99.2% eukaryotic OTUs showed no significant changes after prefiltering with a 200 μm pore-size sieve. We further found that both >0.2 μm and 0.2-200 μm eukaryotic plankton communities, especially the abundant plankton subcommunities, exhibited very similar, and synchronous, spatiotemporal patterns and processes associated with almost identical environmental drivers. The lack of an effect on community structure from prefiltering suggests that environmental DNA from larger metazoa is introduced into the smaller size class. Therefore, size-fractionated filtering with 200 μm is insufficient to discriminate between the eukaryotic plankton size groups in metabarcoding approaches. Our results also highlight the importance of sequencing depth, and strict quality filtering of reads, when designing studies to characterize microeukaryotic plankton communities.

SDH6 and SDH7 Contribute to Anchoring Succinate Dehydrogenase to the Inner Mitochondrial Membrane in Arabidopsis thaliana

The succinate dehydrogenase complex (complex II) is a highly conserved protein complex composed of the SDH1 to SDH4 subunits in bacteria and in the mitochondria of animals and fungi. The reason for the occurrence of up to four additional subunits in complex II of plants, termed SDH5 to SDH8, so far is a mystery. Here, we present a biochemical approach to investigate the internal subunit arrangement of Arabidopsis (Arabidopsis thaliana) complex II. Using low-concentration detergent treatments, the holo complex is dissected into subcomplexes that are analyzed by a three-dimensional gel electrophoresis system. Protein identifications by mass spectrometry revealed that the largest subcomplex (IIa) represents the succinate dehydrogenase domain composed of SDH1 and SDH2. Another subcomplex (IIb) is composed of the SDH3, SDH4, SDH6, and SDH7 subunits. All four proteins include transmembrane helices and together form the membrane anchor of complex II. Sequence analysis revealed that SDH3 and SDH4 lack helices conserved in other organisms. Using homology modeling and phylogenetic analyses, we present evidence that SDH6 and SDH7 substitute missing sequence stretches of SDH3 and SDH4 in plants. Together with SDH5, which is liberated upon dissection of complex II into subcomplexes, SDH6 and SDH7 also add some hydrophilic mass to plant complex II, which possibly inserts further functions into this smallest protein complex of the oxidative phosphorylation system (which is not so small in plants).

Comprehensive transcriptome analysis provides new insights into nutritional strategies and phylogenetic relationships of chrysophytes

Background

Chrysophytes are protist model species in ecology and ecophysiology and important grazers of bacteria-sized microorganisms and primary producers. However, they have not yet been investigated in detail at the molecular level, and no genomic and only little transcriptomic information is available. Chrysophytes exhibit different trophic modes: while phototrophic chrysophytes perform only photosynthesis, mixotrophs can gain carbon from bacterial food as well as from photosynthesis, and heterotrophs solely feed on bacteria-sized microorganisms. Recent phylogenies and megasystematics demonstrate an immense complexity of eukaryotic diversity with numerous transitions between phototrophic and heterotrophic organisms. The question we aim to answer is how the diverse nutritional strategies, accompanied or brought about by a reduction of the plasmid and size reduction in heterotrophic strains, affect physiology and molecular processes.

Results

We sequenced the mRNA of 18 chrysophyte strains on the Illumina HiSeq platform and analysed the transcriptomes to determine relations between the trophic mode (mixotrophic vs. heterotrophic) and gene expression. We observed an enrichment of genes for photosynthesis, porphyrin and chlorophyll metabolism for phototrophic and mixotrophic strains that can perform photosynthesis. Genes involved in nutrient absorption, environmental information processing and various transporters (e.g., monosaccharide, peptide, lipid transporters) were present or highly expressed only in heterotrophic strains that have to sense, digest and absorb bacterial food. We furthermore present a transcriptome-based alignment-free phylogeny construction approach using transcripts assembled from short reads to determine the evolutionary relationships between the strains and the possible influence of nutritional strategies on the reconstructed phylogeny. We discuss the resulting phylogenies in comparison to those from established approaches based on ribosomal RNA and orthologous genes. Finally, we make functionally annotated reference transcriptomes of each strain available to the community, significantly enhancing publicly available data on Chrysophyceae.

Conclusions

Our study is the first comprehensive transcriptomic characterisation of a diverse set of Chrysophyceaen strains. In addition, we showcase the possibility of inferring phylogenies from assembled transcriptomes using an alignment-free approach. The raw and functionally annotated data we provide will prove beneficial for further examination of the diversity within this taxon. Our molecular characterisation of different trophic modes presents a first such example.

Exploring the evolution of the proteins of the plant nuclear envelope

In this study, we explore the plasticity during evolution of proteins of the higher plant nuclear envelope (NE) from the most ancestral plant species to advanced angiosperms. The higher plant NE contains a functional Linker of Nucleoskeleton and Cytoskeleton (LINC) complex based on conserved Sad1-Unc84 (SUN) domain proteins and plant specific Klarsicht/Anc1/Syne homology (KASH) domain proteins. Recent evidence suggests the presence of a plant lamina underneath the inner membrane and various coiled-coil proteins have been hypothesized to be associated with it including Crowded Nuclei (CRWN; also termed LINC and NMCP), Nuclear Envelope Associated Protein (NEAP) protein families as well as the CRWN binding protein KAKU4. SUN domain proteins appear throughout with a key role for mid-SUN proteins suggested. Evolution of KASH domain proteins has resulted in increasing complexity, with some appearing in all species considered, while other KASH proteins are progressively gained during evolution. Failure to identify CRWN homologs in unicellular organisms included in the study and their presence in plants leads us to speculate that convergent evolution may have occurred in the formation of the lamina with each kingdom having new proteins such as the Lamin B receptor (LBR) and Lamin-Emerin-Man1 (LEM) domain proteins (animals) or NEAPs and KAKU4 (plants). Our data support a model in which increasing complexity at the nuclear envelope occurred through the plant lineage and suggest a key role for mid-SUN proteins as an early and essential component of the nuclear envelope.

Evolutionary Origin of Euglena

Euglenids (Excavata, Discoba, Euglenozoa, Euglenida) is a group of free-living, single-celled flagellates living in the aquatic environments. The uniting and unique morphological feature of euglenids is the presence of a cell covering called the pellicle. The morphology and organization of the pellicle correlate well with the mode of nutrition and cell movement. Euglenids exhibit diverse modes of nutrition, including phagotrophy and photosynthesis. Photosynthetic species (Euglenophyceae) constitute a single subclade within euglenids. Their plastids embedded by three membranes arose as the result of a secondary endosymbiosis between phagotrophic eukaryovorous euglenid and the Pyramimonas-related green alga. Within photosynthetic euglenids three evolutionary lineages can be distinguished. The most basal lineage is formed by one mixotrophic species, Rapaza viridis. Other photosynthetic euglenids are split into two groups: predominantly marine Eutreptiales and freshwater Euglenales. Euglenales are divided into two families: Phacaceae, comprising three monophyletic genera (Discoplastis, Lepocinclis, Phacus) and Euglenaceae with seven monophyletic genera (Euglenaformis, Euglenaria, Colacium, Cryptoglena, Strombomonas, Trachelomonas, Monomorphina) and polyphyletic genus Euglena. For 150 years researchers have been studying Euglena based solely on morphological features what resulted in hundreds of descriptions of new taxa and many artificial intra-generic classification systems. In spite of the progress towards defining Euglena, it still remains polyphyletic and morphologically almost undistinguishable from members of the recently described genus Euglenaria; members of both genera have cells undergoing metaboly (dynamic changes in cell shape), large chloroplasts with pyrenoids and monomorphic paramylon grains. Model organisms Euglena gracilis Klebs, the species of choice for addressing fundamental questions in eukaryotic biochemistry, cell and molecular biology, is a representative of the genus Euglena.

Euglena gracilis Genome and Transcriptome: Organelles, Nuclear Genome Assembly Strategies and Initial Features

Euglena gracilis is a major component of the aquatic ecosystem and together with closely related species, is ubiquitous worldwide. Euglenoids are an important group of protists, possessing a secondarily acquired plastid and are relatives to the Kinetoplastidae, which themselves have global impact as disease agents. To understand the biology of E. gracilis, as well as to provide further insight into the evolution and origins of the Kinetoplastidae, we embarked on sequencing the nuclear genome; the plastid and mitochondrial genomes are already in the public domain. Earlier studies suggested an extensive nuclear DNA content, with likely a high degree of repetitive sequence, together with significant extrachromosomal elements. To produce a list of coding sequences we have combined transcriptome data from both published and new sources, as well as embarked on de novo sequencing using a combination of 454, Illumina paired end libraries and long PacBio reads. Preliminary analysis suggests a surprisingly large genome approaching 2 Gbp, with a highly fragmented architecture and extensive repeat composition. Over 80% of the RNAseq reads from E. gracilis maps to the assembled genome sequence, which is comparable with the well assembled genomes of T. brucei and T. cruzi. In order to achieve this level of assembly we employed multiple informatics pipelines, which are discussed here. Finally, as a preliminary view of the genome architecture, we discuss the tubulin and calmodulin genes, which highlight potential novel splicing mechanisms.

Medical Parasitology Taxonomy Update: January 2012 to December 2015

Parasites of medical importance have long been classified taxonomically by morphological characteristics. However, molecular-based techniques have been increasingly used and relied on to determine evolutionary distances for the basis of rational hierarchal classifications. This has resulted in several different classification schemes for parasites and changes in parasite taxonomy. The purpose of this Minireview is to provide a single reference for diagnostic laboratories that summarizes new and revised clinically relevant parasite taxonomy from January 2012 through December 2015.

Expansion of the molecular and morphological diversity of Acanthamoebidae (Centramoebida, Amoebozoa) and identification of a novel life cycle type within the group

BACKGROUND

Acanthamoebidae is a "family" level amoebozoan group composed of the genera Acanthamoeba, Protacanthamoeba, and very recently Luapeleamoeba. This clade of amoebozoans has received considerable attention from the broader scientific community as Acanthamoeba spp. represent both model organisms and human pathogens. While the classical composition of the group (Acanthamoeba + Protacanthamoeba) has been well accepted due to the morphological and ultrastructural similarities of its members, the Acanthamoebidae has never been highly statistically supported in single gene phylogenetic reconstructions of Amoebozoa either by maximum likelihood (ML) or Bayesian analyses.

RESULTS

Here we show using a phylogenomic approach that the Acanthamoebidae is a fully supported monophyletic group within Amoebozoa with both ML and Bayesian analyses. We also expand the known range of morphological and life cycle diversity found in the Acanthamoebidae by demonstrating that the amoebozoans "Protostelium" arachisporum, Dracoamoeba jormungandri n. g. n. sp., and Vacuolamoeba acanthoformis n.g. n.sp., belong within the group. We also found that "Protostelium" pyriformis is clearly a species of Acanthamoeba making it the first reported sporocarpic member of the genus, that is, an amoeba that individually forms a walled, dormant propagule elevated by a non-cellular stalk. Our phylogenetic analyses recover a fully supported Acanthamoebidae composed of five genera. Two of these genera (Acanthamoeba and Luapeleameoba) have members that are sporocarpic.

CONCLUSIONS

Our results provide high statistical support for an Acanthamoebidae that is composed of five distinct genera. This study increases the known morphological diversity of this group and shows that species of Acanthamoeba can include spore-bearing stages. This further illustrates the widespread nature of spore-bearing stages across the tree of Amoebozoa.

REVIEWERS

This article was reviewed by Drs. Eugene Koonin, Purificacion Lopez-Garcia and Sandra Baldauf. Sandra Baldauf was nominated by Purificacion Lopez-Garcia, an Editorial Board member.

A small mitochondrial protein present in myzozoans is essential for malaria transmission

Myzozoans (which include dinoflagellates, chromerids and apicomplexans) display notable divergence from their ciliate sister group, including a reduced mitochondrial genome and divergent metabolic processes. The factors contributing to these divergent processes are still poorly understood and could serve as potential drug targets in disease-causing protists. Here, we report the identification and characterization of a small mitochondrial protein from the rodent-infecting apicomplexan parasite Plasmodium berghei that is essential for development in its mosquito host. Parasites lacking the gene mitochondrial protein ookinete developmental defect (mpodd) showed malformed parasites that were unable to transmit to mosquitoes. Knockout parasites displayed reduced mitochondrial mass without affecting organelle integrity, indicating no role of the protein in mitochondrial biogenesis or morphology maintenance but a likely role in mitochondrial import or metabolism. Using genetic complementation experiments, we identified a previously unrecognized Plasmodium falciparum homologue that can rescue the mpodd(−) phenotype, thereby showing that the gene is functionally conserved. As far as can be detected, mpodd is found in myzozoans, has homologues in the phylum Apicomplexa and appears to have arisen in free-living dinoflagellates. This suggests that the MPODD protein has a conserved mitochondrial role that is important for myzozoans. While previous studies identified a number of essential proteins which are generally highly conserved evolutionarily, our study identifies, for the first time, a non-canonical protein fulfilling a crucial function in the mitochondrion during parasite transmission.

Concerted Up-regulation of Aldehyde/Alcohol Dehydrogenase (ADHE) and Starch in Chlamydomonas reinhardtii Increases Survival under Dark Anoxia

Aldehyde/alcohol dehydrogenases (ADHEs) are bifunctional enzymes that commonly produce ethanol from acetyl-CoA with acetaldehyde as intermediate and play a key role in anaerobic redox balance in many fermenting bacteria. ADHEs are also present in photosynthetic unicellular eukaryotes, where their physiological role and regulation are, however, largely unknown. Herein we provide the first molecular and enzymatic characterization of the ADHE from the photosynthetic microalga Chlamydomonas reinhardtii Purified recombinant ADHE catalyzed the reversible NADH-mediated interconversions of acetyl-CoA, acetaldehyde, and ethanol but seemed to be poised toward the production of ethanol from acetaldehyde. Phylogenetic analysis of the algal fermentative enzyme supports a vertical inheritance from a cyanobacterial-related ancestor. ADHE was located in the chloroplast, where it associated in dimers and higher order oligomers. Electron microscopy analysis of ADHE-enriched stromal fractions revealed fine spiral structures, similar to bacterial ADHE spirosomes. Protein blots showed that ADHE is regulated under oxic conditions. Up-regulation is observed in cells exposed to diverse physiological stresses, including zinc deficiency, nitrogen starvation, and inhibition of carbon concentration/fixation capacity. Analyses of the overall proteome and fermentation profiles revealed that cells with increased ADHE abundance exhibit better survival under dark anoxia. This likely relates to the fact that greater ADHE abundance appeared to coincide with enhanced starch accumulation, which might reflect ADHE-mediated anticipation of anaerobic survival.

Widespread Inter- and Intra-Domain Horizontal Gene Transfer of d-Amino Acid Metabolism Enzymes in Eukaryotes

Analysis of the growing number of available fully-sequenced genomes has shown that Horizontal Gene Transfer (HGT) in eukaryotes is more common than previously thought. It has been proposed that genes with certain functions may be more prone to HGT than others, but we still have a very poor understanding of the selective forces driving eukaryotic HGT. Recent work uncovered that d-amino acid racemases have been commonly transferred from bacteria to fungi, but their role in the receiving organisms is currently unknown. Here, we set out to assess whether d-amino acid racemases are commonly transferred to and between eukaryotic groups. For this we performed a global survey that used a novel automated phylogeny-based HGT-detection algorithm (Abaccus). Our results revealed that at least 7.0% of the total eukaryotic racemase repertoire is the result of inter- or intra-domain HGT. These transfers are significantly enriched in plant-associated fungi. For these, we hypothesize a possible role for the acquired racemases allowing to exploit minoritary nitrogen sources in plant biomass, a nitrogen-poor environment. Finally, we performed experiments on a transferred aspartate-glutamate racemase in the fungal human pathogen Candida glabrata, which however revealed no obvious biological role.

Blastocystis Mitochondrial Genomes Appear to Show Multiple Independent Gains and Losses of Start and Stop Codons

Complete mitochondrion-related organelle (MRO) genomes of several subtypes (STs) of the unicellular stramenopile Blastocystis are presented. Complete conservation of gene content and synteny in gene order is observed across all MRO genomes, comprising 27 protein coding genes, 2 ribosomal RNA genes, and 16 transfer RNA (tRNA) genes. Despite the synteny, differences in the degree of overlap between genes were observed between subtypes and also between isolates within the same subtype. Other notable features include unusual base-pairing mismatches in the predicted secondary structures of some tRNAs. Intriguingly, the rps4 gene in some MRO genomes is missing a start codon and, based on phylogenetic relationships among STs, this loss has happened twice independently. One unidentified open reading frame (orf160) is present in all MRO genomes. However, with the exception of ST4 where the feature has been lost secondarily, orf160 contains variously one or two in-frame stop codons. The overall evidence suggests that both the orf160 and rps4 genes are functional in all STs, but how they are expressed remains unclear.

Bothrosome Formation in Schizochytrium aggregatum (Labyrinthulomycetes, Stramenopiles) during Zoospore Settlement

Labyrinthulomycetes are characterized by the presence of ectoplasmic nets originating from an organelle known as the bothrosome, whose evolutionary origin is unclear. To address this issue, we investigated the developmental process from a zoospore to a vegetative cell in Schizochytrium aggregatum. After disappearance of the flagellum during zoospore settlement, the bothrosome emerged at the anterior-ventral pole of the cells. A new Golgi body also appeared at this stage, and the bothrosome was positioned close to both the new and the old Golgi bodies. This observation suggested that the Golgi body is related to the formation of the bothrosome. Actin appeared as a spot in the same location as the newly appeared bothrosome, as determined by immunofluorescence labeling. An immunoelectron microscopic analysis revealed that actin was present in the ectoplasmic nets and in the cytoplasm around the bothrosome, indicating that the electron-dense materials of the bothrosome are not the polar center of F-actin. This suggests that actin filaments pull the endoplasmic reticulum to the bothrosome and induce the membrane to become evaginated within ectoplasmic nets.

The interaction of a Trypanosoma brucei KH-domain protein with a ribonuclease is implicated in ribosome processing

Ribosomal RNA maturation is best understood in yeast. While substantial efforts have been made to explore parts of these essential pathways in animals, the similarities and uniquenesses of rRNA maturation factors in non-Opisthokonts remain largely unexplored. Eukaryotic ribosome synthesis requires the coordinated activities of hundreds of Assembly Factors (AFs) that transiently associate with pre-ribosomes, many of which are essential. Pno1 and Nob1 are two of six AFs that are required for the cytoplasmic maturation of the 20S pre-rRNA to 18S rRNA in yeast where it has been almost exclusively analyzed. Specifically, Nob1 ribonucleolytic activity generates the mature 3'-end of 18S rRNA. We identified putative Pno1 and Nob1 homologues in the protist Trypanosoma brucei, named TbPNO1 and TbNOB1, and set out to explore their rRNA maturation role further as they are both essential for normal growth. TbPNO1 is a nuclear protein with limited cytosolic localization relative to its yeast homologue. Like in yeast, it interacts directly with TbNOB1, with indications of associations with a larger AF-containing complex. Interestingly, in the absence of TbPNO1, TbNOB1 exhibits non-specific degradation activity on RNA substrates, and its cleavage activity becomes specific only in the presence of TbPNO1, suggesting that TbPNO1-TbNOB1 interaction is essential for regulation and site-specificity of TbNOB1 activity. These results highlight a conserved role of the TbPNO1-TbNOB1 complex in 18S rRNA maturation across eukaryotes; yet reveal a novel role of their interaction in regulation of TbNOB1 enzymatic activity.

An Interactome-Centered Protein Discovery Approach Reveals Novel Components Involved in Mitosome Function and Homeostasis in Giardia lamblia

Protozoan parasites of the genus Giardia are highly prevalent globally, and infect a wide range of vertebrate hosts including humans, with proliferation and pathology restricted to the small intestine. This narrow ecological specialization entailed extensive structural and functional adaptations during host-parasite co-evolution. An example is the streamlined mitosomal proteome with iron-sulphur protein maturation as the only biochemical pathway clearly associated with this organelle. Here, we applied techniques in microscopy and protein biochemistry to investigate the mitosomal membrane proteome in association to mitosome homeostasis. Live cell imaging revealed a highly immobilized array of 30–40 physically distinct mitosome organelles in trophozoites. We provide direct evidence for the single giardial dynamin-related protein as a contributor to mitosomal morphogenesis and homeostasis. To overcome inherent limitations that have hitherto severely hampered the characterization of these unique organelles we applied a novel interaction-based proteome discovery strategy using forward and reverse protein co-immunoprecipitation. This allowed generation of organelle proteome data strictly in a protein-protein interaction context. We built an initial Tom40-centered outer membrane interactome by co-immunoprecipitation experiments, identifying small GTPases, factors with dual mitosome and endoplasmic reticulum (ER) distribution, as well as novel matrix proteins. Through iterative expansion of this protein-protein interaction network, we were able to i) significantly extend this interaction-based mitosomal proteome to include other membrane-associated proteins with possible roles in mitosome morphogenesis and connection to other subcellular compartments, and ii) identify novel matrix proteins which may shed light on mitosome-associated metabolic functions other than Fe-S cluster biogenesis. Functional analysis also revealed conceptual conservation of protein translocation despite the massive divergence and reduction of protein import machinery in Giardia mitosomes.

Organelles with endosymbiotic origin are present in virtually all extant eukaryotes and have undergone considerable remodeling during > 1 billion years of evolution. Highly diverged organelles such as mitosomes or plastids in some parasitic protozoa are the product of extensive secondary reduction. They are sufficiently unique to generate interest as targets for pharmacological intervention, in addition to providing a rich ground for evolutionary cell biologists. The so-called mitochondria-related organelles (MROs) comprise mitosomes and hydrogenosomes, with the former having lost any role in energy metabolism along with the organelle genome. The mitosomes of the intestinal pathogen Giardia lamblia are the most highly reduced MROs known and have proven difficult to investigate because of their extreme divergence and their unique biophysical properties. Here, we implemented a novel strategy aimed at systematic analysis of the organelle proteome by iterative expansion of a protein-protein interaction network. We combined serial forward and reverse co-immunoprecipitations with mass spectrometry analysis, data mining, and validation by subcellular localization and/or functional analysis to generate an interactome network centered on a giardial Tom40 homolog. This iterative ab initio proteome reconstruction provided protein-protein interaction data in addition to identifying novel organelle proteins and functions. Building on this data we generated information on organelle replication, mitosome morphogenesis and organelle dynamics in living cells.

Multiciliated Cells in Animals

Many animal cells assemble single cilia involved in motile and/or sensory functions. In contrast, multiciliated cells (MCCs) assemble up to 300 motile cilia that beat in a coordinate fashion to generate a directional fluid flow. In the human airways, the brain, and the oviduct, MCCs allow mucus clearance, cerebrospinal fluid circulation, and egg transportation, respectively. Impairment of MCC function leads to chronic respiratory infections and increased risks of hydrocephalus and female infertility. MCC differentiation during development or repair involves the activation of a regulatory cascade triggered by the inhibition of Notch activity in MCC progenitors. The downstream events include the simultaneous assembly of a large number of basal bodies (BBs)-from which cilia are nucleated-in the cytoplasm of the differentiating MCCs, their migration and docking at the plasma membrane associated to an important remodeling of the actin cytoskeleton, and the assembly and polarization of motile cilia. The direction of ciliary beating is coordinated both within cells and at the tissue level by a combination of planar polarity cues affecting BB position and hydrodynamic forces that are both generated and sensed by the cilia. Herein, we review the mechanisms controlling the specification and differentiation of MCCs and BB assembly and organization at the apical surface, as well as ciliary assembly and coordination in MCCs.

Plant Proteins Are Smaller Because They Are Encoded by Fewer Exons than Animal Proteins

Protein size is an important biochemical feature since longer proteins can harbor more domains and therefore can display more biological functionalities than shorter proteins. We found remarkable differences in protein length, exon structure, and domain count among different phylogenetic lineages. While eukaryotic proteins have an average size of 472 amino acid residues (aa), average protein sizes in plant genomes are smaller than those of animals and fungi. Proteins unique to plants are ∼81aa shorter than plant proteins conserved among other eukaryotic lineages. The smaller average size of plant proteins could neither be explained by endosymbiosis nor subcellular compartmentation nor exon size, but rather due to exon number. Metazoan proteins are encoded on average by ∼10 exons of small size [∼176 nucleotides (nt)]. Streptophyta have on average only ∼5.7 exons of medium size (∼230nt). Multicellular species code for large proteins by increasing the exon number, while most unicellular organisms employ rather larger exons (>400nt). Among subcellular compartments, membrane proteins are the largest (∼520aa), whereas the smallest proteins correspond to the gene ontology group of ribosome (∼240aa). Plant genes are encoded by half the number of exons and also contain fewer domains than animal proteins on average. Interestingly, endosymbiotic proteins that migrated to the plant nucleus became larger than their cyanobacterial orthologs. We thus conclude that plants have proteins larger than bacteria but smaller than animals or fungi. Compared to the average of eukaryotic species, plants have ∼34% more but ∼20% smaller proteins. This suggests that photosynthetic organisms are unique and deserve therefore special attention with regard to the evolutionary forces acting on their genomes and proteomes.

Recapitulating phylogenies using k-mers: from trees to networks

Ernst Haeckel based his landmark Tree of Life on the supposed ontogenic recapitulation of phylogeny, i.e. that successive embryonic stages during the development of an organism re-trace the morphological forms of its ancestors over the course of evolution. Much of this idea has since been discredited. Today, phylogenies are often based on families of molecular sequences. The standard approach starts with a multiple sequence alignment, in which the sequences are arranged relative to each other in a way that maximises a measure of similarity position-by-position along their entire length. A tree (or sometimes a network) is then inferred. Rigorous multiple sequence alignment is computationally demanding, and evolutionary processes that shape the genomes of many microbes (bacteria, archaea and some morphologically simple eukaryotes) can add further complications. In particular, recombination, genome rearrangement and lateral genetic transfer undermine the assumptions that underlie multiple sequence alignment, and imply that a tree-like structure may be too simplistic. Here, using genome sequences of 143 bacterial and archaeal genomes, we construct a network of phylogenetic relatedness based on the number of shared k-mers (subsequences at fixed length k). Our findings suggest that the network captures not only key aspects of microbial genome evolution as inferred from a tree, but also features that are not treelike. The method is highly scalable, allowing for investigation of genome evolution across a large number of genomes. Instead of using specific regions or sequences from genome sequences, or indeed Haeckel’s idea of ontogeny, we argue that genome phylogenies can be inferred using k-mers from whole-genome sequences. Representing these networks dynamically allows biological questions of interest to be formulated and addressed quickly and in a visually intuitive manner.

MyTH4-FERM myosins have an ancient and conserved role in filopod formation

The formation of filopodia in Metazoa and Amoebozoa requires the activity of myosin 10 (Myo10) in mammalian cells and of Dictyostelium unconventional myosin 7 (DdMyo7) in the social amoeba Dictyostelium However, the exact roles of these MyTH4-FERM myosins (myosin tail homology 4-band 4.1, ezrin, radixin, moesin; MF) in the initiation and elongation of filopodia are not well defined and may reflect conserved functions among phylogenetically diverse MF myosins. Phylogenetic analysis of MF myosin domains suggests that a single ancestral MF myosin existed with a structure similar to DdMyo7, which has two MF domains, and that subsequent duplications in the metazoan lineage produced its functional homolog Myo10. The essential functional features of the DdMyo7 myosin were identified using quantitative live-cell imaging to characterize the ability of various mutants to rescue filopod formation in myo7-null cells. The two MF domains were found to function redundantly in filopod formation with the C-terminal FERM domain regulating both the number of filopodia and their elongation velocity. DdMyo7 mutants consisting solely of the motor plus a single MyTH4 domain were found to be capable of rescuing the formation of filopodia, establishing the minimal elements necessary for the function of this myosin. Interestingly, a chimeric myosin with the Myo10 MF domain fused to the DdMyo7 motor also was capable of rescuing filopod formation in the myo7-null mutant, supporting fundamental functional conservation between these two distant myosins. Together, these findings reveal that MF myosins have an ancient and conserved role in filopod formation.

Evolution in the Cycles of Life

The life cycles of eukaryotes alternate between haploid and diploid phases, which are initiated by meiosis and gamete fusion, respectively. In both ascomycete and basidiomycete fungi and chlorophyte algae, the haploid-to-diploid transition is regulated by a pair of paralogous homeodomain protein encoding genes. That a common genetic program controls the haploid-to-diploid transition in phylogenetically disparate eukaryotic lineages suggests this may be the ancestral function for homeodomain proteins. Multicellularity has evolved independently in many eukaryotic lineages in either one or both phases of the life cycle. Organisms, such as land plants, exhibiting a life cycle whereby multicellular bodies develop in both the haploid and diploid phases are often referred to as possessing an alternation of generations. We review recent progress on understanding the genetic basis for the land plant alternation of generations and highlight the roles that homeodomain-encoding genes may have played in the evolution of complex multicellularity in this lineage.

Integrated omics study of lipid droplets from Plasmodiophora brassicae

Plasmodiophora brassicae causes clubroot disease in cruciferous. In this report, lipid droplets were observed in the resting spores of P. brassicae. 295 lipid droplet-associated proteins were identified and categorized into nine groups. Transcriptome analysis of these proteins during three different zoosporic stages revealed differences in gene expression pattern. GO enrichment analysis revealed that these proteins associated with lipid droplets were mainly linked to biosynthesis and metabolism. GC-MS analysis revealed that lipid droplets contain seven types of free fatty acids: saturated fatty acids C16:0 and C18:0, and unsaturated fatty acids C18:1Δ9, C18:1Δ11, C18:2, C20:4 and C20:5. P. brassicae accumulated a large amount of triacylglycerols (TAGs). We systematically analyzed the putative proteins involved in TAG biosynthesis and its metabolic pathway. KEGG pathway analysis defined 3390 genes, including 167 genes involved in lipid metabolism. Transcriptome analysis revealed that 162 candidate enzymes involved in lipid metabolism were differential expressed. Our omics studies are the first to investigate the lipid droplet organelles in P. brassicae, providing a reference resource to study protist lipid droplets.

Polymetallic nodules, sediments, and deep waters in the equatorial North Pacific exhibit highly diverse and distinct bacterial, archaeal, and microeukaryotic communities

Concentrated seabed deposits of polymetallic nodules, which are rich in economically valuable metals (e.g., copper, nickel, cobalt, manganese), occur over vast areas of the abyssal Pacific Ocean floor. Little is currently known about the diversity of microorganisms inhabiting abyssal habitats. In this study, sediment, nodule, and water column samples were collected from the Clarion-Clipperton Zone of the Eastern North Pacific. The diversities of prokaryote and microeukaryote communities associated with these habitats were examined. Microbial community composition and diversity varied with habitat type, water column depth, and sediment horizon. Thaumarchaeota were relatively enriched in the sediments and nodules compared to the water column, whereas Gammaproteobacteria were the most abundant sequences associated with nodules. Among the Eukaryota, rRNA genes belonging to the Cryptomonadales were relatively most abundant among organisms associated with nodules, whereas rRNA gene sequences deriving from members of the Alveolata were relatively enriched in sediments and the water column. Nine operational taxonomic unit (OTU)s were identified that occur in all nodules in this dataset, as well as all nodules found in a study 3000-9000 km from our site. Microbial communities in the sediments had the highest diversity, followed by nodules, and then by the water column with <1/3 the number of OTUs as in the sediments.