Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia

A taxon-rich and genome-scale phylogeny of Opisthokonta

Ancient divergences within Opisthokonta-a major lineage that includes organisms in the kingdoms Animalia, Fungi, and their unicellular relatives-remain contentious. To assess progress toward a genome-scale Opisthokonta phylogeny, we conducted the most taxon rich phylogenomic analysis using sets of genes inferred with different orthology inference methods and established the geological timeline of Opisthokonta diversification. We also conducted sensitivity analysis by subsampling genes or taxa from the full data matrix based on filtering criteria previously shown to improve phylogenomic inference. We found that approximately 85% of internal branches were congruent across data matrices and the approaches used. Notably, the use of different orthology inference methods was a substantial contributor to the observed incongruence: analyses using the same set of orthologs showed high congruence of 97% to 98%, whereas different sets of orthologs resulted in somewhat lower congruence (87% to 91%). Examination of unicellular Holozoa relationships suggests that the instability observed across varying gene sets may stem from weak phylogenetic signals. Our results provide a comprehensive Opisthokonta phylogenomic framework that will be useful for illuminating ancient evolutionary episodes concerning the origin and diversification of the 2 major eukaryotic kingdoms and emphasize the importance of investigating effects of orthology inference on phylogenetic analyses to resolve ancient divergences.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Inference of phylogenetic trees directly from raw sequencing reads using Read2Tree

Current methods for inference of phylogenetic trees require running complex pipelines at substantial computational and labor costs, with additional constraints in sequencing coverage, assembly and annotation quality, especially for large datasets. To overcome these challenges, we present Read2Tree, which directly processes raw sequencing reads into groups of corresponding genes and bypasses traditional steps in phylogeny inference, such as genome assembly, annotation and all-versus-all sequence comparisons, while retaining accuracy. In a benchmark encompassing a broad variety of datasets, Read2Tree is 10-100 times faster than assembly-based approaches and in most cases more accurate-the exception being when sequencing coverage is high and reference species very distant. Here, to illustrate the broad applicability of the tool, we reconstruct a yeast tree of life of 435 species spanning 590 million years of evolution. We also apply Read2Tree to >10,000 Coronaviridae samples, accurately classifying highly diverse animal samples and near-identical severe acute respiratory syndrome coronavirus 2 sequences on a single tree. The speed, accuracy and versatility of Read2Tree enable comparative genomics at scale.

Interkingdom network analyses reveal microalgae and protostomes as keystone taxa involved in nutrient cycling in large freshwater lake sediment

Few studies have explored the role of interkingdom interactions between bacteria and microeukaryotes in nutrient cycling in lake ecosystems. We conducted sediment sampling from 40 locations covering Hongze Lake and analyzed their chemical properties. Intra- and interkingdom networks were constructed using 16S and 18S rRNA gene amplicon sequencing. Microeukaryotic intranetworks were more complex in spring than in autumn, while no clear variation in the complexity of bacterial intranetworks was found between autumn and spring. Larger and more complex bacterial-microeukaryotic bipartite networks emerged in spring than in autumn, correlated with lower carbon, nitrogen, and phosphorus levels in spring, likely resulting in intense microbial competition. Bacteria and microeukaryotes played different topological roles in interkingdom networks, with microeukaryotes contributing to the networks' greater complexity. Seven keystone modules were identified in spring and autumn nutrient cycling. Importantly, keystone taxa in these modules belonged to photoautotrophic microalgae or predatory protostomes, indicating that these organisms are key drivers in lake sediment nutrient cycling. Our results suggested that nutrient content variation in autumn and spring changes interkingdom networks' topological structure between bacteria and microeukaryotes. Microalgae and protostomes are essential in freshwater lake nutrient cycling and may be targeted to modulate nutrient cycling in large freshwater ecosystems.

Compositionally Constrained Sites Drive Long-Branch Attraction

Accurate phylogenies are fundamental to our understanding of the pattern and process of evolution. Yet, phylogenies at deep evolutionary timescales, with correspondingly long branches, have been fraught with controversy resulting from conflicting estimates from models with varying complexity and goodness of fit. Analyses of historical as well as current empirical datasets, such as alignments including Microsporidia, Nematoda, or Platyhelminthes, have demonstrated that inadequate modeling of across-site compositional heterogeneity, which is the result of biochemical constraints that lead to varying patterns of accepted amino acids along sequences, can lead to erroneous topologies that are strongly supported. Unfortunately, models that adequately account for across-site compositional heterogeneity remain computationally challenging or intractable for an increasing fraction of contemporary datasets. Here, we introduce "compositional constraint analysis," a method to investigate the effect of site-specific constraints on amino acid composition on phylogenetic inference. We show that more constrained sites with lower diversity and less constrained sites with higher diversity exhibit ostensibly conflicting signals under models ignoring across-site compositional heterogeneity that lead to long-branch attraction artifacts and demonstrate that more complex models accounting for across-site compositional heterogeneity can ameliorate this bias. We present CAT-posterior mean site frequencies (PMSF), a pipeline for diagnosing and resolving phylogenetic bias resulting from inadequate modeling of across-site compositional heterogeneity based on the CAT model. CAT-PMSF is robust against long-branch attraction in all alignments we have examined. We suggest using CAT-PMSF when convergence of the CAT model cannot be assured. We find evidence that compositionally constrained sites are driving long-branch attraction in two metazoan datasets and recover evidence for Porifera as the sister group to all other animals. [Animal phylogeny; cross-site heterogeneity; long-branch attraction; phylogenomics.].

Identifying the Best Approximating Model in Bayesian Phylogenetics: Bayes Factors, Cross-Validation or wAIC?

There is still no consensus as to how to select models in Bayesian phylogenetics, and more generally in applied Bayesian statistics. Bayes factors are often presented as the method of choice, yet other approaches have been proposed, such as cross-validation or information criteria. Each of these paradigms raises specific computational challenges, but they also differ in their statistical meaning, being motivated by different objectives: either testing hypotheses or finding the best-approximating model. These alternative goals entail different compromises, and as a result, Bayes factors, cross-validation, and information criteria may be valid for addressing different questions. Here, the question of Bayesian model selection is revisited, with a focus on the problem of finding the best-approximating model. Several model selection approaches were re-implemented, numerically assessed and compared: Bayes factors, cross-validation (CV), in its different forms (k-fold or leave-one-out), and the widely applicable information criterion (wAIC), which is asymptotically equivalent to leave-one-out cross-validation (LOO-CV). Using a combination of analytical results and empirical and simulation analyses, it is shown that Bayes factors are unduly conservative. In contrast, CV represents a more adequate formalism for selecting the model returning the best approximation of the data-generating process and the most accurate estimates of the parameters of interest. Among alternative CV schemes, LOO-CV and its asymptotic equivalent represented by the wAIC, stand out as the best choices, conceptually and computationally, given that both can be simultaneously computed based on standard Markov chain Monte Carlo runs under the posterior distribution. [Bayes factor; cross-validation; marginal likelihood; model comparison; wAIC.].

Read2Tree: scalable and accurate phylogenetic trees from raw reads

The inference of phylogenetic trees is foundational to biology. However, state-of-the-art phylogenomics requires running complex pipelines, at significant computational and labour costs, with additional constraints in sequencing coverage, assembly and annotation quality. To overcome these challenges, we present Read2Tree, which directly processes raw sequencing reads into groups of corresponding genes. In a benchmark encompassing a broad variety of datasets, our assembly-free approach was 10-100x faster than conventional approaches, and in most cases more accurate-the exception being when sequencing coverage was high and reference species very distant. To illustrate the broad applicability of the tool, we reconstructed a yeast tree of life of 435 species spanning 590 million years of evolution. Applied to Coronaviridae samples, Read2Tree accurately classified highly diverse animal samples and near-identical SARS-CoV-2 sequences on a single tree-thereby exhibiting remarkable breadth and depth. The speed, accuracy, and versatility of Read2Tree enables comparative genomics at scale.

Host identity is the dominant factor in the assembly of nematode and tardigrade gut microbiomes in Antarctic Dry Valley streams

Recent work examining nematode and tardigrade gut microbiomes has identified species-specific relationships between host and gut community composition. However, only a handful of species from either phylum have been examined. How microbiomes differ among species and what factors contribute to their assembly remains unexplored. Cyanobacterial mats within Antarctic Dry Valley streams host a simple and tractable natural ecosystem of identifiable microinvertebrates to address these questions. We sampled 2 types of coexisting mats (i.e., black and orange) across four spatially isolated streams, hand-picked single individuals of two nematode species (i.e., Eudorylaimus antarcticus and Plectus murrayi) and tardigrades, to examine their gut microbiomes using 16S and 18S rRNA metabarcoding. All gut microbiomes (bacterial and eukaryotic) were significantly less diverse than the mats they were isolated from. In contrast to mats, microinvertebrates' guts were depleted of Cyanobacteria and differentially enriched in taxa of Bacteroidetes, Proteobacteria, and Fungi. Among factors investigated, gut microbiome composition was most influenced by host identity while environmental factors (e.g., mats and streams) were less important. The importance of host identity in predicting gut microbiome composition suggests functional value to the host, similar to other organisms with strong host selected microbiomes.

Out of chaos: Phylogenomics of Asian Sonerileae

Sonerileae is a diverse Melastomataceae lineage comprising ca. 1000 species in 44 genera, with >70% of genera and species distributed in Asia. Asian Sonerileae are taxonomically intractable with obscure generic circumscriptions. The backbone phylogeny of this group remains poorly resolved, possibly due to complexity caused by rapid species radiation in early and middle Miocene, which hampers further systematic study. Here, we used genome resequencing data to reconstruct the phylogeny of Asian Sonerileae. Three parallel datasets, viz. single-copy ortholog (SCO), genomic SNPs, and whole plastome, were assembled from genome resequencing data of 205 species for this purpose. Based on these genome-scale data, we provided the first well resolved phylogeny of Asian Sonerileae, with 34 major clades identified and 74% of the interclade relationships consistently resolved by both SCO and genomic data. Meanwhile, widespread phylogenetic discordance was detected among SCO gene trees as well as species trees reconstructed using different tree estimation methods (concatenation/site-based coalescent method/summary method) or different datasets (SCO/genomic/plastome). We explored sources of discordance using multiple approaches and found that the observed discordance in Asian Sonerileae was mainly caused by a combination of biased distribution of missing data, random noise from uninformative genes, incomplete lineage sorting, and hybridization/introgression. Exploration of these sources can enable us to generate hypotheses for future testing, which is the first step towards understanding the evolution of Asian Sonerileae. We also detected high levels of homoplasy for some characters traditionally used in taxonomy, which explains current chaotic generic delimitations. The backbone phylogeny of Asian Sonerileae revealed in this study offers a solid basis for future taxonomic revision at the generic level.

Evolutionary origin and sequence signatures of the heterodimeric ABCG5/ABCG8 transporter

ATP-binding cassette (ABC) systems, characterized by ABC-type nucleotide-binding domains (NBDs), play crucial roles in various aspects of human physiology. Human ABCG5 and ABCG8 form a heterodimeric transporter that functions in the efflux of sterols. We used sequence similarity search, multiple sequence alignment, phylogenetic analysis, and structure comparison to study the evolutionary origin and sequence signatures of ABCG5 and ABCG8. Orthologs of ABCG5 and ABCG8, supported by phylogenetic analysis and signature residues, were identified in bilaterian animals, Filasterea, Fungi, and Amoebozoa. Such a phylogenetic distribution suggests that ABCG5 and ABCG8 could have originated in the last common ancestor of Amorphea (the unikonts), the eukaryotic group including Amoebozoa and Opisthokonta. ABCG5 and ABCG8 were missing in genomes of various lineages such as snakes, jawless vertebrates, non-vertebrate chordates, echinoderms, and basal metazoan groups. Amino-acid changes in key positions in ABCG8 Walker A motif and/or ABCG5 C-loop were observed in most tetrapod organisms, likely resulted in the loss of ATPase activity at one nucleotide-binding site. ABCG5 and ABCG8 in Ecdysozoa (such as insects) exhibit elevated evolutionary rates and accumulate various changes in their NBD functional motifs. Alignment inspection revealed several residue positions that show different amino-acid usages in ABCG5/ABCG8 compared to other ABCG subfamily proteins. These residues were mapped to the structural cores of transmembrane domains (TMDs), the NBD-TMD interface, and the interface between TMDs. They serve as sequence signatures to differentiate ABCG5/ABCG8 from other ABCG subfamily proteins, and some of them may contribute to substrate specificity of the ABCG5/ABCG8 transporter.

Novel Dicistroviruses in an Unexpected Wide Range of Invertebrates

Dicistroviruses are members of a rapidly growing family of small RNA viruses. Related sequences have been discovered in many environmental samples, indicating that our knowledge about dicistrovirus diversity and host range is still limited. In this study, we performed a systematic search against the publicly available transcriptome database, and identified large numbers of dicistrovirus-like sequences in a wide variety of eukaryotic species. The origins of these sequences were 108 invertebrates (including 77 insect species belonging to 18 orders) and 11 plants, revealing new associations between dicistroviruses and hosts. Finally, 83 transcripts corresponding to nearly-complete viral genomes were retrieved from the RNA-seq data, of which most sequences showed limited similarity to known dicistroviruses and might present previously unreported virus species. Phylogenetic analysis suggested that horizontal virus transfer has occurred between diverse hosts and has important implications for dicistrovirus evolution. The results will provide new insight into the hidden diversity of the Dicistroviridae, and help us to better understand the viral evolution, host range and the possible way of transmission.

Concordance-based approaches for the inference of relationships and molecular rates with phylogenomic datasets

Gene tree conflict is common and finding methods to analyze and alleviate the negative effects that conflict has on species tree analysis is a crucial part of phylogenomics. This study aims to expand the discussion of inferring species trees and molecular branch lengths when conflict is present. Conflict is typically examined in two ways: inferring its prevalence, and inferring the influence of the individual genes (how strongly one gene supports any given topology compared to an alternative topology). Here, we examine a procedure for incorporating both conflict and the influence of genes in order to infer evolutionary relationships. All supported relationships in the gene trees are analyzed and the likelihood of the genes constrained to these relationships is summed to provide a likelihood for the relationship. Consensus tree assembly is conducted based on the sum of likelihoods for a given relationship and choosing relationships based on the most likely relationship assuming it does not conflict with a relationship that has a higher likelihood score. If it is not possible for all most likely relationships to be combined into a single bifurcating tree then multiple trees are produced and a consensus tree with a polytomy is created. This procedure allows for more influential genes to have greater influence on an inferred relationship, does not assume conflict has arisen from any one source, and does not force the dataset to produce a single bifurcating tree. Using this approach on three empirical datasets, we examine and discuss the relationship between influence and prevalence of gene tree conflict. We find that in one of the datasets, assembling a bifurcating consensus tree solely composed of the most likely relationships is impossible. To account for conflict in molecular rate analysis we also introduce a concordance-based approach to the summary and estimation of branch lengths suitable for downstream comparative analyses. We demonstrate through simulation that even under high levels of stochastic conflict, the mean and median of the concordant rates recapitulate the true molecular rate better than using a supermatrix approach. Using a large phylogenomic dataset, we examine rate heterogeneity across concordant genes with a focus on the branch subtending crown angiosperms. Notably, we find highly variable rates of evolution along the branch subtending crown angiosperms. The approaches outlined here have several limitations, but they also represent some alternative methods for harnessing the complexity of phylogenomic datasets and enrich our inferences of both species' relationships and evolutionary processes.

Human SIRT1 Multispecificity Is Modulated by Active-Site Vicinity Substitutions during Natural Evolution

Many enzymes that catalyze protein post-translational modifications can specifically modify multiple target proteins. However, little is known regarding the molecular basis and evolution of multispecificity in these enzymes. Here, we used a combined bioinformatics and experimental approaches to investigate the evolution of multispecificity in the sirtuin-1 (SIRT1) deacetylase. Guided by bioinformatics analysis of SIRT1 orthologs and substrates, we identified and examined important amino acid substitutions that have occurred during the evolution of sirtuins in Metazoa and Fungi. We found that mutation of human SIRT1 at these positions, based on sirtuin orthologs from Fungi, could alter its substrate specificity. These substitutions lead to reduced activity toward K382 acetylated p53 protein, which is only present in Metazoa, without affecting the high activity toward the conserved histone substrates. Results from ancestral sequence reconstruction are consistent with a model in which ancestral sirtuin proteins exhibited multispecificity, suggesting that the multispecificity of some metazoan sirtuins, such as hSIRT1, could be a relatively ancient trait.

Scalable Empirical Mixture Models That Account for Across-Site Compositional Heterogeneity

Biochemical demands constrain the range of amino acids acceptable at specific sites resulting in across-site compositional heterogeneity of the amino acid replacement process. Phylogenetic models that disregard this heterogeneity are prone to systematic errors, which can lead to severe long-branch attraction artifacts. State-of-the-art models accounting for across-site compositional heterogeneity include the CAT model, which is computationally expensive, and empirical distribution mixture models estimated via maximum likelihood (C10–C60 models). Here, we present a new, scalable method EDCluster for finding empirical distribution mixture models involving a simple cluster analysis. The cluster analysis utilizes specific coordinate transformations which allow the detection of specialized amino acid distributions either from curated databases or from the alignment at hand. We apply EDCluster to the HOGENOM and HSSP databases in order to provide universal distribution mixture (UDM) models comprising up to 4,096 components. Detailed analyses of the UDM models demonstrate the removal of various long-branch attraction artifacts and improved performance compared with the C10–C60 models. Ready-to-use implementations of the UDM models are provided for three established software packages (IQ-TREE, Phylobayes, and RevBayes).

Evidence for sponges as sister to all other animals from partitioned phylogenomics with mixture models and recoding

Resolving the relationships between the major lineages in the animal tree of life is necessary to understand the origin and evolution of key animal traits. Sponges, characterized by their simple body plan, were traditionally considered the sister group of all other animal lineages, implying a gradual increase in animal complexity from unicellularity to complex multicellularity. However, the availability of genomic data has sparked tremendous controversy as some phylogenomic studies support comb jellies taking this position, requiring secondary loss or independent origins of complex traits. Here we show that incorporating site-heterogeneous mixture models and recoding into partitioned phylogenomics alleviates systematic errors that hamper commonly-applied phylogenetic models. Testing on real datasets, we show a great improvement in model-fit that attenuates branching artefacts induced by systematic error. We reanalyse key datasets and show that partitioned phylogenomics does not support comb jellies as sister to other animals at either the supermatrix or partition-specific level.

Resolving the Early Divergence Pattern of Teleost Fish Using Genome-Scale Data

Regarding the phylogenetic relationship of the three primary groups of teleost fishes, Osteoglossomorpha (bonytongues and others), Elopomorpha (eels and relatives), Clupeocephala (the remaining teleost fish), early morphological studies hypothesized the first divergence of Osteoglossomorpha, whereas the recent prevailing view is the first divergence of Elopomorpha. Molecular studies supported all the possible relationships of the three primary groups. This study analyzed genome-scale data from four previous studies: 1) 412 genes from 12 species, 2) 772 genes from 15 species, 3) 1,062 genes from 30 species, and 4) 491 UCE loci from 27 species. The effects of the species, loci, and models used on the constructed tree topologies were investigated. In the analyses of the data sets (1)–(3), although the first divergence of Clupeocephala that left the other two groups in a sister relationship was supported by concatenated sequences and gene trees of all the species and genes, the first divergence of Elopomorpha among the three groups was supported using species and/or genes with low divergence of sequence and amino-acid frequencies. This result corresponded to that of the UCE data set (4), whose sequence divergence was low, which supported the first divergence of Elopomorpha with high statistical significance. The increase in accuracy of the phylogenetic construction by using species and genes with low sequence divergence was predicted by a phylogenetic informativeness approach and confirmed by computer simulation. These results supported that Elopomorpha was the first basal group of teleost fish to have diverged, consistent with the prevailing view of recent morphological studies.

Mitochondrial gene order of the freshwater bryozoan Cristatella mucedo retains ancestral lophotrochozoan features

Bryozoans are aquatic colonial suspension-feeders abundant in many marine and freshwater benthic communities. At the same time, the phylum is under studied on both morphological and molecular levels, and its position on the metazoan tree of life is still disputed. Bryozoa include the exclusively marine Stenolaemata, predominantly marine Gymnolaemata and exclusively freshwater Phylactolaemata. Here we report the mitochondrial genome of the phylactolaemate bryozoan Cristatella mucedo. This species has the largest (21,008 bp) of all currently known bryozoan mitogenomes, containing a typical metazoan gene compendium as well as a number of non-coding regions, three of which are longer than 1500 bp. The trnS1/trnG/nad3 region is presumably duplicated in this species. Comparative analysis of the gene order in C. mucedo and another phylactolaemate bryozoan, Pectinatella magnifica, confirmed their close relationships, and revealed a stronger similarity to mitogenomes of phoronids and other lophotrochozoan species than to marine bryozoans, indicating the ancestral nature of their gene arrangement. We suggest that the ancestral gene order underwent substantial changes in different bryozoan cladesshowing mosaic distribution of conservative gene blocks regardless of their phylogenetic position. Altogether, our results support the early divergence of Phylactolaemata from the rest of Bryozoa.

cognac: rapid generation of concatenated gene alignments for phylogenetic inference from large, bacterial whole genome sequencing datasets

Background

The quantity of genomic data is expanding at an increasing rate. Tools for phylogenetic analysis which scale to the quantity of available data are required. To address this need, we present cognac, a user-friendly software package to rapidly generate concatenated gene alignments for phylogenetic analysis.

Results

We illustrate that cognac is able to rapidly identify phylogenetic marker genes using a data driven approach and efficiently generate concatenated gene alignments for very large genomic datasets. To benchmark our tool, we generated core gene alignments for eight unique genera of bacteria, including a dataset of over 11,000 genomes from the genus Escherichia producing an alignment with 1353 genes, which was constructed in less than 17 h.

Conclusions

We demonstrate that cognac presents an efficient method for generating concatenated gene alignments for phylogenetic analysis. We have released cognac as an R package (https://github.com/rdcrawford/cognac) with customizable parameters for adaptation to diverse applications.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

The phylogeny of the African wood mice (Muridae, Hylomyscus) based on complete mitochondrial genomes and five nuclear genes reveals their evolutionary history and undescribed diversity

Wood mice of the genus Hylomyscus, are small-sized rodents widely distributed in lowland and montane rainforests in tropical Africa, where they can be locally abundant. Recent morphological and molecular studies have increased the number of recognized species from 8 to 18 during the last 15 years. We used complete mitochondrial genomes and five nuclear genes to infer the number of candidate species within this genus and depict its evolutionary history. In terms of gene sampling and geographical and taxonomic coverage, this is the most comprehensive review of the genus Hylomyscus to date. The six species groups (aeta, alleni, anselli, baeri, denniae and parvus) defined on morphological grounds are monophyletic. Species delimitation analyses highlight undescribed diversity within this genus: perhaps up to 10 taxa need description or elevation from synonymy, pending review of type specimens. Our divergence dating and biogeographical analyses show that diversification of the genus occurred after the end of the Miocene and is closely linked to the history of the African forest. The formation of the Rift Valley combined with the declining global temperatures during the Late Miocene caused the fragmentation of the forests and explains the first split between the denniae group and remaining lineages. Subsequently, periods of increased climatic instability during Plio-Pleistocene probably resulted in elevated diversification in both lowland and montane forest taxa.

Origin and Evolution of the Neuroendocrine Control of Reproduction in Vertebrates, With Special Focus on Genome and Gene Duplications

In humans, as in the other mammals, the neuroendocrine control of reproduction is ensured by the brain-pituitary gonadotropic axis. Multiple internal and environmental cues are integrated via brain neuronal networks, ultimately leading to the modulation of the activity of gonadotropin-releasing hormone (GnRH) neurons. The decapeptide GnRH is released into the hypothalamic-hypophysial portal blood system and stimulates the production of pituitary glycoprotein hormones, the two gonadotropins luteinizing hormone and follicle-stimulating hormone. A novel actor, the neuropeptide kisspeptin, acting upstream of GnRH, has attracted increasing attention in recent years. Other neuropeptides, such as gonadotropin-inhibiting hormone/RF-amide related peptide, and other members of the RF-amide peptide superfamily, as well as various nonpeptidic neuromediators such as dopamine and serotonin also provide a large panel of stimulatory or inhibitory regulators. This paper addresses the origin and evolution of the vertebrate gonadotropic axis. Brain-pituitary neuroendocrine axes are typical of vertebrates, the pituitary gland, mediator and amplifier of brain control on peripheral organs, being a vertebrate innovation. The paper reviews, from molecular and functional perspectives, the evolution across vertebrate radiation of some key actors of the vertebrate neuroendocrine control of reproduction and traces back their origin along the vertebrate lineage and in other metazoa before the emergence of vertebrates. A focus is given on how gene duplications, resulting from either local events or from whole genome duplication events, and followed by paralogous gene loss or conservation, might have shaped the evolutionary scenarios of current families of key actors of the gonadotropic axis.

Homology modeling identified for purported drug targets to the neuroprotective effects of levodopa and asiaticoside-D in degenerated cerebral ganglions of Lumbricus terrestris

CONTEXT:

Homology modeling plays role in determining the therapeutic targets dreadful for condition such as neurodegenerative diseases (NDD), which pose challenge in achieving the effective managements. The structures of the serotonin transporter (SERT), aquaporin (AQP), and tropomyosin receptor kinase (TrkA) which are implicated in NDD pathology are still unknown for Lumbricus terrestris, but the three-dimensional (3D) structure of the human counterpart for modeling.

AIM:

This study aims to generate and evaluate the 3D structure of TrkA, SERT, and AQP proteins and their interaction with the ligands, namely Asiaticoside-D (AD) and levodopa (L-DOPA) the anti-NDD agents.

SUBJECTS AND METHODS:

Homology modeling for SERT, AQP, and TrkA proteins of Lumbricus terrestris using SWISS-MODEL Server and the modeled structure was validated using Rampage Server. Wet-lab analysis of their correspondent m-RNA levels was also done to validate the in silico data.

RESULTS:

It was found that TrkA had moderately high homology (67%) to human while SERT and AQP could exhibit 58% and 42%, respectively. The reliability of the model was assessed by Ramachandran plot analysis. Interactions of AD with the SERT, AQP-4, and TrkA showed the binding energies as −9.93, 8.88, and −7.58 of Kcal/mol, respectively, while for L-DOPA did show −3.93, −5.13, and −6.0 Kcal/mol, respectively. The levels of SERT, TrkA, and AQP-4 were significantly reduced (P < 0.001) on ROT induced when compared to those of control worms. On ROT + AD supplementation group (III), m-RNA levels were significantly increased (P < 0.05) when compared to those of ROT induced worms (group II).

CONCLUSION:

Our pioneering docking data propose the possible of target which is proved useful for therapeutic investigations against the unconquered better of NDD.

Unexpected rDNA divergence between two morphologically minimalistic nematodes with description of a new species (Tylenchomorpha: Tylenchidae)

Species of the family Tylenchidae are encountered in large numbers in soils. The genus Labrys was recently described characterised by a remarkable lip pattern that differs from all other known Tylenchidae genera. Here we describe a curious new species, Labrys fujianensis sp. n., that morphologically fits the genus Labrys but which is genetically divergent. The phylogeny was inferred based on 18S and 28S rDNA and light and scanning electron microscopy were used to extract detailed morphologies. The phylogenetic position of this species and its phenotypic convergence are discussed. The possibility of a long-branch attraction artefact was inspected both by removal of variable nucleotide sites and monophyletic testing of topologies. The results confirmed the divergent positioning of the presented species and it is demonstrated that the genetic diversity in Tylenchidae may be much higher than expected due to morphological homoplasy.

Noise and biases in genomic data may underlie radically different hypotheses for the position of Iguania within Squamata

Squamate reptiles are a major component of vertebrate biodiversity whose crown-clade traces its origin to a narrow window of time in the Mesozoic during which the main subclades diverged in rapid succession. Deciphering phylogenetic relationships among these lineages has proven challenging given the conflicting signals provided by genomic and phenomic data. Most notably, the placement of Iguania has routinely differed between data sources, with morphological evidence supporting a sister relationship to the remaining squamates (Scleroglossa hypothesis) and molecular data favoring a highly nested position alongside snakes and anguimorphs (Toxicofera hypothesis). We provide novel insights by generating an expanded morphological dataset and exploring the presence of phylogenetic signal, noise, and biases in molecular data. Our analyses confirm the presence of strong conflicting signals for the position of Iguania between morphological and molecular datasets. However, we also find that molecular data behave highly erratically when inferring the deepest branches of the squamate tree, a consequence of limited phylogenetic signal to resolve this ancient radiation with confidence. This, in turn, seems to result from a rate of evolution that is too high for historical signals to survive to the present. Finally, we detect significant systematic biases, with iguanians and snakes sharing faster rates of molecular evolution and a similarly biased nucleotide composition. A combination of scant phylogenetic signal, high levels of noise, and the presence of systematic biases could result in the misplacement of Iguania. We regard this explanation to be at least as plausible as the complex scenario of convergence and reversals required for morphological data to be misleading. We further evaluate and discuss the utility of morphological data to resolve ancient radiations, as well as its impact in combined-evidence phylogenomic analyses, with results relevant for the assessment of evidence and conflict across the Tree of Life.

Genetic and molecular basis of the immune system in the brachiopod Lingula anatina

The extension of comparative immunology to non-model systems, such as mollusks and annelids, has revealed an unexpected diversity in the complement of immune receptors and effectors among evolutionary lineages. However, several lophotrochozoan phyla remain unexplored mainly due to the lack of genomic resources. The increasing accessibility of high-throughput sequencing technologies offers unique opportunities for extending genome-wide studies to non-model systems. As a result, the genome-based study of the immune system in brachiopods allows a better understanding of the alternative survival strategies developed by these immunologically neglected phyla. Here we present a detailed overview of the molecular components of the immune system identified in the genome of the brachiopod Lingula anatina. Our findings reveal conserved intracellular signaling pathways as well as unique strategies for pathogen detection and killing in brachiopods.

Phylogenomics of Lophotrochozoa with Consideration of Systematic Error

Phylogenomic studies have improved understanding of deep metazoan phylogeny and show promise for resolving incongruences among analyses based on limited numbers of loci. One region of the animal tree that has been especially difficult to resolve, even with phylogenomic approaches, is relationships within Lophotrochozoa (the animal clade that includes molluscs, annelids, and flatworms among others). Lack of resolution in phylogenomic analyses could be due to insufficient phylogenetic signal, limitations in taxon and/or gene sampling, or systematic error. Here, we investigated why lophotrochozoan phylogeny has been such a difficult question to answer by identifying and reducing sources of systematic error. We supplemented existing data with 32 new transcriptomes spanning the diversity of Lophotrochozoa and constructed a new set of Lophotrochozoa-specific core orthologs. Of these, 638 orthologous groups (OGs) passed strict screening for paralogy using a tree-based approach. In order to reduce possible sources of systematic error, we calculated branch-length heterogeneity, evolutionary rate, percent missing data, compositional bias, and saturation for each OG and analyzed increasingly stricter subsets of only the most stringent (best) OGs for these five variables. Principal component analysis of the values for each factor examined for each OG revealed that compositional heterogeneity and average patristic distance contributed most to the variance observed along the first principal component while branch-length heterogeneity and, to a lesser extent, saturation contributed most to the variance observed along the second. Missing data did not strongly contribute to either. Additional sensitivity analyses examined effects of removing taxa with heterogeneous branch lengths, large amounts of missing data, and compositional heterogeneity. Although our analyses do not unambiguously resolve lophotrochozoan phylogeny, we advance the field by reducing the list of viable hypotheses. Moreover, our systematic approach for dissection of phylogenomic data can be applied to explore sources of incongruence and poor support in any phylogenomic data set. [Annelida; Brachiopoda; Bryozoa; Entoprocta; Mollusca; Nemertea; Phoronida; Platyzoa; Polyzoa; Spiralia; Trochozoa.].

The male sexual apparatus in the order Scorpiones (Arachnida): a comparative study of functional morphology as a tool to define hypotheses of homology

BACKGROUND

Insemination in scorpions is carried out by means of a partly sclerotized structure, the spermatophore, which is composed of two separate halves, the hemispermatophores. In most genera these reproductive structures can be used to differentiate species. However, many taxa such as the genus Euscorpius and the family Diplocentridae lack the morphological diversity observed in the copulatory organs of many other arthropods, rendering them useless for species level taxonomy. Such structural stasis, however, suggests that hemispermatophores have evolved relatively slowly and may thus provide a stronger phylogenetic signal for recognizing supra-generic ranks than previously thought. Based on the postulate that the phenotypic stability observed in some groups is the consequence of functional constraint, the most comprehensive comparative study of the male sexual apparatus to date was conducted for a complete reassessment of the morphology, phylogenetic value and hypotheses of homology of these structures.

RESULTS

Hemispermatophores, pre- and post-insemination spermatophores, as well as the inherent mechanisms of insemination, were studied across the whole order, allowing the recognition and description of a series of five basic bauplans for the capsular region. For the most part, these patterns appear to be consistent within each major taxonomic group, but several cases of incongruence between spermatophore morphology and taxonomy raises questions about the monophyly of some clades. The Bothriuridae are traditionally regarded as a basal scorpionoid family. However, except for the genus Lisposoma, bothriurid hemispermatophores and spermatophores are morphologically more similar to those of the Chactoidea than to those of scorpionoids. On the other hand, the male copulatory structures of the hormurid clade (Hormiops (Hormurus + Liocheles)) are more akin to those of Diplocentridae and Heteroscorpionidae than to those of other hormurids.

CONCLUSIONS

Spermatophore capsular patterns appears to be congruent with a recent phylogeny of the order Scorpiones based on phylogenomic data that placed Bothriuridae outside of Scorpionoidea and Liocheles outside of Hormuridae, in contradicton with earlier phylogenetic reconstructions based on morphology. This raises questions about the potential use of functionally constrained traits to assess the reliability of contradicting phylogenetic hypotheses and emphasizes the need for a thorough reassessment of the scorpion phylogenetic relationships.

Metagenomics reshapes the concepts of RNA virus evolution by revealing extensive horizontal virus transfer

Virus metagenomics is a young research filed but it has already transformed our understanding of virus diversity and evolution, and illuminated at a new level the connections between virus evolution and the evolution and ecology of the hosts. In this review article, we examine the new picture of the evolution of RNA viruses, the dominant component of the eukaryotic virome, that is emerging from metagenomic data analysis. The major expansion of many groups of RNA viruses through metagenomics allowed the construction of substantially improved phylogenetic trees for the conserved virus genes, primarily, the RNA-dependent RNA polymerases (RdRp). In particular, a new superfamily of widespread, small positive-strand RNA viruses was delineated that unites tombus-like and noda-like viruses. Comparison of the genome architectures of RNA viruses discovered by metagenomics and by traditional methods reveals an extent of gene module shuffling among diverse virus genomes that far exceeds the previous appreciation of this evolutionary phenomenon. Most dramatically, inclusion of the metagenomic data in phylogenetic analyses of the RdRp resulted in the identification of numerous, strongly supported groups that encompass RNA viruses from diverse hosts including different groups of protists, animals and plants. Notwithstanding potential caveats, in particular, incomplete and uneven sampling of eukaryotic taxa, these highly unexpected findings reveal horizontal virus transfer (HVT) between diverse hosts as the central aspect of RNA virus evolution. The vast and diverse virome of invertebrates, particularly nematodes and arthropods, appears to be the reservoir, from which the viromes of plants and vertebrates evolved via multiple HVT events.

Reconstructing the muscular ground pattern of phylactolaemate bryozoans: first data from gelatinous representatives

BACKGROUND

Phylactolaemata is commonly regarded the earliest branch within Bryozoa and thus the sister group to the other bryozoan taxa, Cyclostomata and Gymnolaemata. Therefore, the taxon is important for the reconstruction of the bryozoan morphological ground pattern. In this study the myoanatomy of Pectinatella magnifica, Cristatella mucedo and Hyalinella punctata was analysed by means of histology, f-actin staining and confocal laser-scanning microscopy in order to fill gaps in knowledge concerning the myoanatomy of Phylactolaemata.

RESULTS

The retractor muscles and muscles of the aperture, gut, body wall, tentacle sheath, lophophore constitute the most prominent muscular subsets in these species. The lophophore shows longitudinal muscle bands in the tentacles, lophophoral arm muscles, epistome musculature and hitherto undescribed muscles of the ring canal. In general the muscular system of the three species is very similar with differences mainly in the body wall, tentacle sheath and epistome. The body wall contains an orthogonal grid of musculature. The epistome exhibits either a muscular meshwork in the epistomal wall or muscle fibers traversing the epistomal cavity. The whole tentacle sheath possesses a regular mesh of muscles in Pectinatella and Cristatella, whereas circular muscles are limited to the tentacle sheath base in Hyalinella.

CONCLUSION

This study is the first to describe muscles of the ring canal and contributes to reconstructing muscular features for the last common ancestor of all bryozoans. The data available suggest that two longitudinal muscle bands in the tentacles, as well as retractor muscles and longitudinal and circular muscles in the tentacle sheath, were present in the last common bryozoan ancestor. Comparisons among bryozoans shows that several apomorphies are present in the myoanatomy of each class- level taxon such as the epistomal musculature and musculature of the lophophoral arms in phylactolaemates, annular muscles in cyclostomes and parietal muscles in gymnolaemates.

The Composite 259-kb Plasmid of Martelella mediterranea DSM 17316T-A Natural Replicon with Functional RepABC Modules from Rhodobacteraceae and Rhizobiaceae

A multipartite genome organization with a chromosome and many extrachromosomal replicons (ECRs) is characteristic for Alphaproteobacteria. The best investigated ECRs of terrestrial rhizobia are the symbiotic plasmids for legume root nodulation and the tumor-inducing (Ti) plasmid of Agrobacterium tumefaciens. RepABC plasmids represent the most abundant alphaproteobacterial replicon type. The currently known homologous replication modules of rhizobia and Rhodobacteraceae are phylogenetically distinct. In this study, we surveyed type-strain genomes from the One Thousand Microbial Genomes (KMG-I) project and identified a roseobacter-specific RepABC-type operon in the draft genome of the marine rhizobium Martelella mediterranea DSM 17316^T. PacBio genome sequencing demonstrated the presence of three circular ECRs with sizes of 593, 259, and 170-kb. The rhodobacteral RepABC module is located together with a rhizobial equivalent on the intermediate sized plasmid pMM259, which likely originated in the fusion of a pre-existing rhizobial ECR with a conjugated roseobacter plasmid. Further evidence for horizontal gene transfer (HGT) is given by the presence of a roseobacter-specific type IV secretion system on the 259-kb plasmid and the rhodobacteracean origin of 62% of the genes on this plasmid. Functionality tests documented that the genuine rhizobial RepABC module from the Martelella 259-kb plasmid is only maintained in A. tumefaciens C58 (Rhizobiaceae) but not in Phaeobacter inhibens DSM 17395 (Rhodobacteraceae). Unexpectedly, the roseobacter-like replication system is functional and stably maintained in both host strains, thus providing evidence for a broader host range than previously proposed. In conclusion, pMM259 is the first example of a natural plasmid that likely mediates genetic exchange between roseobacters and rhizobia.

First fossil of an oestroid fly (Diptera: Calyptratae: Oestroidea) and the dating of oestroid divergences

Calyptrate flies include about 22,000 extant species currently classified into Hippoboscoidea (tsetse, louse, and bat flies), the muscoid grade (house flies and relatives) and the Oestroidea (blow flies, bot flies, flesh flies, and relatives). Calyptrates are abundant in nearly all terrestrial ecosystems, often playing key roles as decomposers, parasites, parasitoids, vectors of pathogens, and pollinators. For oestroids, the most diverse group within calyptrates, definitive fossils have been lacking. The first unambiguous fossil of Oestroidea is described based on a specimen discovered in amber from the Dominican Republic. The specimen was identified through digital dissection by CT scans, which provided morphological data for a cladistic analysis of its phylogenetic position among extant oestroids. The few known calyptrate fossils were used as calibration points for a molecular phylogeny (16S, 28S, CAD) to estimate the timing of major diversification events among the Oestroidea. Results indicate that: (a) the fossil belongs to the family Mesembrinellidae, and it is identified and described as Mesembrinella caenozoica sp. nov.; (b) the mesembrinellids form a sister clade to the Australian endemic Ulurumyia macalpinei (Ulurumyiidae) (McAlpine's fly), which in turn is sister to all remaining oestroids; (c) the most recent common ancestor of extant Calyptratae lived just before the K-Pg boundary (ca. 70 mya); and (d) the radiation of oestroids began in the Eocene (ca. 50 mya), with the origin of the family Mesembrinellidae dated at ca. 40 mya. These results provide new insight into the timing and rate of oestroid diversification and highlight the rapid radiation of some of the most diverse and ecologically important families of flies. ZooBank accession number-urn:lsid:zoobank.org:pub:0DC5170B-1D16-407A-889E-56EED3FE3627.

Identification of introns harboring functional sequence elements through positional conservation

Many human introns carry out a function, in the sense that they are critical to maintain normal cellular activity. Their identification is fundamental to understanding cellular processes and disease. However, being noncoding elements, such functional introns are poorly predicted based on traditional approaches of sequence and structure conservation. Here, we generated a dataset of human functional introns that carry out different types of functions. We showed that functional introns share common characteristics, such as higher positional conservation along the coding sequence and reduced loss rates, regardless of their specific function. A unique property of the data is that if an intron is unknown to be functional, it still does not mean that it is indeed non-functional. We developed a probabilistic framework that explicitly accounts for this unique property, and predicts which specific human introns are functional. We show that we successfully predict function even when the algorithm is trained on introns with a different type of function. This ability has many implications in studying regulatory networks, gene regulation, the effect of mutations outside exons on human disease, and on our general understanding of intron evolution and their functional exaptation in mammals.

Segmentation in Tardigrada and diversification of segmental patterns in Panarthropoda

The origin and diversification of segmented metazoan body plans has fascinated biologists for over a century. The superphylum Panarthropoda includes three phyla of segmented animals-Euarthropoda, Onychophora, and Tardigrada. This superphylum includes representatives with relatively simple and representatives with relatively complex segmented body plans. At one extreme of this continuum, euarthropods exhibit an incredible diversity of serially homologous segments. Furthermore, distinct tagmosis patterns are exhibited by different classes of euarthropods. At the other extreme, all tardigrades share a simple segmented body plan that consists of a head and four leg-bearing segments. The modular body plans of panarthropods make them a tractable model for understanding diversification of animal body plans more generally. Here we review results of recent morphological and developmental studies of tardigrade segmentation. These results complement investigations of segmentation processes in other panarthropods and paleontological studies to illuminate the earliest steps in the evolution of panarthropod body plans.

The impact of GC bias on phylogenetic accuracy using targeted enrichment phylogenomic data

The field of sequence based phylogenetic analyses is currently being transformed by novel hybrid-based targeted enrichment methods, such as the use of ultraconserved elements (UCEs). Rather than analyzing relationships among organisms using a small number of genes, these methods now allow us to evaluate relationships with many hundreds to thousands of individual gene loci. However, the inclusion of thousands of loci does not necessarily overcome the long-standing challenge of incongruence among phylogenetic trees derived from different genes or gene regions. One factor that impacts the level of incongruence in phylogenomic data sets is the level of GC bias. GC rich gene regions are prone to higher recombination rates than AT rich regions, driven by a process referred to as "GC biased gene conversion". As a result, high GC content can be negatively associated with phylogenetic accuracy, but the extent to which this impacts incongruence among UCEs is currently unstudied. We investigated the impact of GC content on phylogeny reconstruction using in silico captured UCE data for the corbiculate bees (Hymenoptera: Apidae). The phylogeny of this group has been the subject of extensive study, and incongruence among gene trees is thought to be a source of phylogenetic error. We conducted coalescent- and concatenation-based analyses of 810 individual gene loci from all 13 currently available bee genomes, including 8 corbiculate taxa. Both coalescent- and concatenation-based methods converged on a single topology for the corbiculate tribes. In contrast to concatenation, the coalescent-based methods revealed significant topological conflict at nodes involving the orchid bees (Euglossini) and honeybees (Apini). Partitioning the loci by GC content reveals decreasing support for the inferred topology with increasing GC bias. Based on the results of this study, we report the first evidence that GC biased gene conversion may contribute to topological incongruence in studies based on ultraconserved elements.

A Large and Consistent Phylogenomic Dataset Supports Sponges as the Sister Group to All Other Animals

Resolving the early diversification of animal lineages has proven difficult, even using genome-scale datasets. Several phylogenomic studies have supported the classical scenario in which sponges (Porifera) are the sister group to all other animals ("Porifera-sister" hypothesis), consistent with a single origin of the gut, nerve cells, and muscle cells in the stem lineage of eumetazoans (bilaterians + ctenophores + cnidarians). In contrast, several other studies have recovered an alternative topology in which ctenophores are the sister group to all other animals (including sponges). The "Ctenophora-sister" hypothesis implies that eumetazoan-specific traits, such as neurons and muscle cells, either evolved once along the metazoan stem lineage and were then lost in sponges and placozoans or evolved at least twice independently in Ctenophora and in Cnidaria + Bilateria. Here, we report on our reconstruction of deep metazoan relationships using a 1,719-gene dataset with dense taxonomic sampling of non-bilaterian animals that was assembled using a semi-automated procedure, designed to reduce known error sources. Our dataset outperforms previous metazoan gene superalignments in terms of data quality and quantity. Analyses with a best-fitting site-heterogeneous evolutionary model provide strong statistical support for placing sponges as the sister-group to all other metazoans, with ctenophores emerging as the second-earliest branching animal lineage. Only those methodological settings that exacerbated long-branch attraction artifacts yielded Ctenophora-sister. These results show that methodological issues must be carefully addressed to tackle difficult phylogenetic questions and pave the road to a better understanding of how fundamental features of animal body plans have emerged.

Cross-talk among immune and neuroendocrine systems in molluscs and other invertebrate models

The comparison between immune and neuroendocrine systems in vertebrates and invertebrates suggest an ancient origin and a high degree of conservation for the mechanisms underlying the integration between immune and stress responses. This suggests that in both vertebrates and invertebrates the stress response involves the integrated network of soluble mediators (e.g., neurotransmitters, hormones and cytokines) and cell functions (e.g., chemotaxis and phagocytosis), that interact with a common objective, i.e., the maintenance of body homeostasis. During evolution, several changes observed in the stress response of more complex taxa could be the result of new roles of ancestral molecules, such as ancient immune mediators may have been recruited as neurotransmitters and hormones, or vice versa. We review older and recent evidence suggesting that immune and neuro-endocrine functions during the stress response were deeply intertwined already at the dawn of multicellular organisms. These observations found relevant reflections in the demonstration that immune cells can transdifferentiate in olfactory neurons in crayfish and the recently re-proposed neural transdifferentiation in humans.

Multiple measures could alleviate long-branch attraction in phylogenomic reconstruction of Cupressoideae (Cupressaceae)

Long-branch attraction (LBA) is a major obstacle in phylogenetic reconstruction. The phylogenetic relationships among Juniperus (J), Cupressus (C) and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) subclades of Cupressoideae are controversial. Our initial analyses of plastid protein-coding gene matrix revealed both J and C with much longer stem branches than those of HCX, so their sister relationships may be attributed to LBA. We used multiple measures including data filtering and modifying, evolutionary model selection and coalescent phylogenetic reconstruction to alleviate the LBA artifact. Data filtering by strictly removing unreliable aligned regions and removing substitution saturation genes and rapidly evolving sites could significantly reduce branch lengths of subclades J and C and recovered a relationship of J (C, HCX). In addition, using coalescent phylogenetic reconstruction could elucidate the LBA artifact and recovered J (C, HCX). However, some valid methods for other taxa were inefficient in alleviating the LBA artifact in J-C-HCX. Different strategies should be carefully considered and justified to reduce LBA in phylogenetic reconstruction of different groups. Three subclades of J-C-HCX were estimated to have experienced ancient rapid divergence within a short period, which could be another major obstacle in resolving relationships. Furthermore, our plastid phylogenomic analyses fully resolved the intergeneric relationships of Cupressoideae.

Phylogenetic analysis of ionotropic L-glutamate receptor genes in the Bilateria, with special notes on Aplysia californica

Background

The neurotransmitter L-Glutamate (L-Glu) acting at ionotropic L-Glu receptors (iGluR) conveys fast excitatory signal transmission in the nervous systems of all animals. iGluR-dependent neurotransmission is a key component of the synaptic plasticity that underlies learning and memory. During learning, two subtypes of iGluR, α-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPAR) and N-methyl-D-aspartate receptors (NMDAR), are dynamically regulated postsynaptically in vertebrates. Invertebrate organisms such as Aplysia californica (Aplysia) are well-studied models for iGluR-mediated function, yet no studies to date have analyzed the evolutionary relationships between iGluR genes in these species and those in vertebrates, to identify genes that may mediate plasticity. We conducted a thorough phylogenetic analysis spanning Bilateria to elucidate these relationships. The expression status of iGluR genes in the Aplysia nervous system was also examined.

Results

Our analysis shows that ancestral genes for both NMDAR and AMPAR subtypes were present in the common bilaterian ancestor. NMDAR genes show very high conservation in motifs responsible for forming the conductance pore of the ion channel. The number of NMDAR subunits is greater in vertebrates due to an increased number of splice variants and an increased number of genes, likely due to gene duplication events. AMPAR subunits form an orthologous group, and there is high variability in the number of AMPAR genes in each species due to extensive taxon specific gene gain and loss. qPCR results show that all 12 Aplysia iGluR subunits are expressed in all nervous system ganglia.

Conclusions

Orthologous NMDAR subunits in all species studied suggests conserved function across Bilateria, and potentially a conserved mechanism of neuroplasticity and learning. Vertebrates display an increased number of NMDAR genes and splice variants, which may play a role in their greater diversity of physiological responses. Extensive gene gain and loss of AMPAR genes may result in different physiological properties that are taxon specific. Our results suggest a significant role for L-Glu mediated responses throughout the Aplysia nervous system, consistent with L-Glu’s role as the primary excitatory neurotransmitter.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0871-1) contains supplementary material, which is available to authorized users.

Outgroup, alignment and modelling improvements indicate that two TNFSF13-like genes existed in the vertebrate ancestor

The molecular machinery required for lymphocyte development and differentiation appears to have emerged concomitantly with distinct B- and T-like lymphocyte subsets in the ancestor of all vertebrates. The TNFSF superfamily (TNFSF) members BAFF (TNFSF13/Blys) and APRIL (TNFSF13) are key regulators of B cell development survival, and activation in mammals, but the temporal emergence of these molecules, and their precise relationship to the newly identified TNFSF gene BALM (BAFF and APRIL-like molecule), have not yet been elucidated. Here, to resolve the early evolutionary history of this family, we improved outgroup sampling and alignment quality, and applied better fitting substitution models compared to past studies. Our analyses reveal that BALM is a definitive TNFSF13 family member, which split from BAFF in the gnathostome (jawed vertebrate) ancestor. Most importantly, however, we show that both the APRIL and BAFF lineages existed in the ancestors of all extant vertebrates. This implies that APRIL has been lost, or is yet to be found, in cyclostomes (jawless vertebrates). Our results suggest that lineage-specific gene duplication and loss events have caused lymphocyte regulation, despite shared origins, to become secondarily distinct between gnathostomes and cyclostomes. Finally, the structure of lamprey BAFF-like, and its phylogenetic placement as sister to BAFF and BALM, but not the more slowly evolving APRIL, indicates that the primordial lymphocyte regulator was more APRIL-like than BAFF-like.

Redefining the invertebrate RNA virosphere

Current knowledge of RNA virus biodiversity is both biased and fragmentary, reflecting a focus on culturable or disease-causing agents. Here we profile the transcriptomes of over 220 invertebrate species sampled across nine animal phyla and report the discovery of 1,445 RNA viruses, including some that are sufficiently divergent to comprise new families. The identified viruses fill major gaps in the RNA virus phylogeny and reveal an evolutionary history that is characterized by both host switching and co-divergence. The invertebrate virome also reveals remarkable genomic flexibility that includes frequent recombination, lateral gene transfer among viruses and hosts, gene gain and loss, and complex genomic rearrangements. Together, these data present a view of the RNA virosphere that is more phylogenetically and genomically diverse than that depicted in current classification schemes and provide a more solid foundation for studies in virus ecology and evolution.

Evolutionary switches between two serine codon sets are driven by selection

Serine is the only amino acid that is encoded by two disjoint codon sets so that a tandem substitution of two nucleotides is required to switch between the two sets. Previously published evidence suggests that, for the most evolutionarily conserved serines, the codon set switch occurs by simultaneous substitution of two nucleotides. Here we report a genome-wide reconstruction of the evolution of serine codons in triplets of closely related species from diverse prokaryotes and eukaryotes. The results indicate that the great majority of codon set switches proceed by two consecutive nucleotide substitutions, via a threonine or cysteine intermediate, and are driven by selection. These findings imply a strong pressure of purifying selection in protein evolution, which in the case of serine codon set switches occurs via an initial deleterious substitution quickly followed by a second, compensatory substitution. The result is frequent reversal of amino acid replacements and, at short evolutionary distances, pervasive homoplasy.

Phylogenetic Resolution of Deep Eukaryotic and Fungal Relationships Using Highly Conserved Low-Copy Nuclear Genes

A comprehensive and reliable eukaryotic tree of life is important for many aspects of biological studies from comparative developmental and physiological analyses to translational medicine and agriculture. Both gene-rich and taxon-rich approaches are effective strategies to improve phylogenetic accuracy and are greatly facilitated by marker genes that are universally distributed, well conserved, and orthologous among divergent eukaryotes. In this article, we report the identification of 943 low-copy eukaryotic genes and we show that many of these genes are promising tools in resolving eukaryotic phylogenies, despite the challenges of determining deep eukaryotic relationships. As a case study, we demonstrate that smaller subsets of ∼20 and 52 genes could resolve controversial relationships among widely divergent taxa and provide strong support for deep relationships such as the monophyly and branching order of several eukaryotic supergroups. In addition, the use of these genes resulted in fungal phylogenies that are congruent with previous phylogenomic studies that used much larger datasets, and successfully resolved several difficult relationships (e.g., forming a highly supported clade with Microsporidia, Mitosporidium and Rozella sister to other fungi). We propose that these genes are excellent for both gene-rich and taxon-rich analyses and can be applied at multiple taxonomic levels and facilitate a more complete understanding of the eukaryotic tree of life.