Search

Selecting the best model of sequence evolution for a multiple-sequence-alignment (MSA) constitutes the first step of phylogenetic tree reconstruction. Common approaches for inferring nucleotide models typically apply maximum likelihood (ML) methods, with discrimination between models determined by one of several information criteria. This requires tree reconstruction and optimisation which can be computationally expensive. We demonstrate that neural networks can be used to perform model selection, without the need to reconstruct trees, optimise parameters, or calculate likelihoods. We introduce ModelRevelator, a model selection tool underpinned by two deep neural networks. The first neural network, NNmodelfind, recommends one of six commonly used models of sequence evolution, ranging in complexity from Jukes and Cantor to General Time Reversible. The second, NNalphafind, recommends whether or not a Γ-distributed rate heterogeneous model should be incorporated, and if so, provides an estimate of the shape parameter, ɑ. Users can simply input an MSA into ModelRevelator, and swiftly receive output recommending the evolutionary model, inclusive of the presence or absence of rate heterogeneity, and an estimate of ɑ. We show that ModelRevelator performs comparably with likelihood-based methods and the recently published machine learning method ModelTeller over a wide range of parameter settings, with significant potential savings in computational effort. Further, we show that this performance is not restricted to the alignments on which the networks were trained, but is maintained even on unseen empirical data. We expect that ModelRevelator will provide a valuable alternative for phylogeneticists, especially where traditional methods of model selection are computationally prohibitive., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Efficient repair of DNA lesions is essential for the faithful transmission of genetic information between somatic cells and for genome integrity across generations. Plants have multiple, partially redundant, and overlapping DNA repair pathways, probably due to the less constricted germline and the inevitable exposure to light including higher energy wavelengths. Many proteins involved in DNA repair and their mode of actions are well described. In contrast, a role for DNA damage-associated RNA components, evident from many other organisms, is less well understood. Here, we have challenged young Arabidopsis thaliana plants with two different types of genotoxic stress and performed de novo assembly and transcriptome analysis. We identified three long noncoding RNAs (lncRNAs) that are lowly or not expressed under regular conditions but up-regulated or induced by DNA damage. We generated CRISPR/Cas deletion mutants and found that the absence of the lncRNAs impairs the recovery capacity of the plants from genotoxic stress. The genetic loci are highly conserved among world-wide distributed Arabidopsis accessions and within related species in the Brassicaceae group. Together, these results suggest that the lncRNAs have a conserved function in connection with DNA damage and provide a basis for mechanistic analysis of their role., Competing Interests: Conflicts of interest: The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

Transcription factors are crucial drivers of cellular differentiation during animal development and often share ancient evolutionary origins. The T-box transcription factor Brachyury plays a pivotal role as an early mesoderm determinant and neural repressor in vertebrates; yet, the ancestral function and key evolutionary transitions of the role of this transcription factor remain obscure. Here, we present a genome-wide target-gene screen using chromatin immunoprecipitation sequencing in the sea anemone Nematostella vectensis, an early branching non-bilaterian, and the sea urchin Strongylocentrotus purpuratus, a representative of the sister lineage of chordates. Our analysis reveals an ancestral gene regulatory feedback loop connecting Brachyury, FoxA and canonical Wnt signalling involved in axial patterning that predates the cnidarian-bilaterian split about 700 million years ago. Surprisingly, we also found that part of the gene regulatory network controlling the fate of neuromesodermal progenitors in vertebrates was already present in the common ancestor of cnidarians and bilaterians. However, while several endodermal and neuronal Brachyury target genes are ancestrally shared, hardly any of the key mesodermal downstream targets in vertebrates are found in the sea anemone or the sea urchin. Our study suggests that a limited number of target genes involved in mesoderm formation were newly acquired in the vertebrate lineage, leading to a dramatic shift in the function of this ancestral developmental regulator., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Maximum likelihood and maximum parsimony are two key methods for phylogenetic tree reconstruction. Under certain conditions, each of these two methods can perform more or less efficiently, resulting in unresolved or disputed phylogenies. We show that a neural network can distinguish between four-taxon alignments that were evolved under conditions susceptible to either long-branch attraction or long-branch repulsion. When likelihood and parsimony methods are discordant, the neural network can provide insight as to which tree reconstruction method is best suited to the alignment. When applied to the contentious case of Strepsiptera evolution, our method shows robust support for the current scientific view, that is, it places Strepsiptera with beetles, distant from flies., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

IQ-TREE (http://www.iqtree.org, last accessed February 6, 2020) is a user-friendly and widely used software package for phylogenetic inference using maximum likelihood. Since the release of version 1 in 2014, we have continuously expanded IQ-TREE to integrate a plethora of new models of sequence evolution and efficient computational approaches of phylogenetic inference to deal with genomic data. Here, we describe notable features of IQ-TREE version 2 and highlight the key advantages over other software., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The recent discovery and structure determination of a novel ubiquitin-like dimerization domain in protein kinase D (PKD) has significant implications for its activation. PKD is a serine/threonine kinase activated by the lipid second messenger diacylglycerol (DAG). It is an essential and highly conserved protein that is implicated in plasma membrane directed trafficking processes from the trans-Golgi network. However, many open questions surround its mechanism of activation, its localization, and its role in the biogenesis of cargo transport carriers. In reviewing this field, the focus is primarily on the mechanisms that control the activation of PKD at precise locations in the cell. In light of the new structural findings, the understanding of the mechanisms underlying PKD activation is critically evaluated, with particular emphasis on the role of dimerization in PKD autophosphorylation, and the provenance and recognition of the DAG that activates PKD., (© 2020 The Authors. BioEssays published by WILEY Periodicals, Inc.)

Climate change is causing arctic regions to warm disproportionally faster than those at lower latitudes, leading to alterations in carbon and nitrogen cycling, and potentially higher greenhouse gas emissions. It is thus increasingly important to better characterize the microorganisms driving arctic biogeochemical processes and their potential responses to changing conditions. Here, we describe a novel thaumarchaeon enriched from an arctic soil, Candidatus Nitrosocosmicus arcticus strain Kfb, which has been maintained for seven years in stable laboratory enrichment cultures as an aerobic ammonia oxidizer, with ammonium or urea as substrates. Genomic analyses show that this organism harbors all genes involved in ammonia oxidation and in carbon fixation via the 3-hydroxypropionate/4-hydroxybutyrate cycle, characteristic of all AOA, as well as the capability for urea utilization and potentially also for heterotrophic metabolism, similar to other AOA. Ca . N. arcticus oxidizes ammonia optimally between 20 and 28°C, well above average temperatures in its native high arctic environment (-13-4°C). Ammonia oxidation rates were nevertheless much lower than those of most cultivated mesophilic AOA (20-45°C). Intriguingly, we repeatedly observed apparent faster growth rates (based on marker gene counts) at lower temperatures (4-8°C) but without detectable nitrite production. Together with potential metabolisms predicted from its genome content, these observations indicate that Ca . N. arcticus is not a strict chemolithotrophic ammonia oxidizer and add to cumulating evidence for a greater metabolic and physiological versatility of AOA. The physiology of Ca . N. arcticus suggests that increasing temperatures might drastically affect nitrification in arctic soils by stimulating archaeal ammonia oxidation.

The conserved INO80 chromatin remodeling complex is involved in regulation of DNA damage repair, replication and transcription. It is commonly recruited to the transcription start region and contributes to the establishment of promoter-proximal nucleosomes. We find a substantial influence of INO80 on nucleosome dynamics and gene expression during stress induced transcription. Transcription induced by osmotic stress leads to genome-wide remodeling of promoter proximal nucleosomes. INO80 function is required for timely return of evicted nucleosomes to the 5΄ end of induced genes. Reduced INO80 function in Arp8-deficient cells leads to correlated prolonged transcription and nucleosome eviction. INO80 and the related complex SWR1 regulate incorporation of the H2A.Z isoform at promoter proximal nucleosomes. However, H2A.Z seems not to influence osmotic stress induced gene regulation. Furthermore, we show that high rates of transcription promote INO80 recruitment to promoter regions, suggesting a connection between active transcription and promoter proximal nucleosome remodeling. In addition, we find that absence of INO80 enhances bidirectional promoter activity at highly induced genes and expression of a number of stress induced transcripts. We suggest that INO80 has a direct repressive role via promoter proximal nucleosome remodeling to limit high levels of transcription in yeast., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The evolution of reptiles, birds, and mammals was associated with the origin of unique integumentary structures. Studies on lizards, chicken, and humans have suggested that the evolution of major structural proteins of the outermost, cornified layers of the epidermis was driven by the diversification of a gene cluster called Epidermal Differentiation Complex (EDC). Turtles have evolved unique defense mechanisms that depend on mechanically resilient modifications of the epidermis. To investigate whether the evolution of the integument in these reptiles was associated with specific adaptations of the sequences and expression patterns of EDC-related genes, we utilized newly available genome sequences to determine the epidermal differentiation gene complement of turtles. The EDC of the western painted turtle (Chrysemys picta bellii) comprises more than 100 genes, including at least 48 genes that encode proteins referred to as beta-keratins or corneous beta-proteins. Several EDC proteins have evolved cysteine/proline contents beyond 50% of total amino acid residues. Comparative genomics suggests that distinct subfamilies of EDC genes have been expanded and partly translocated to loci outside of the EDC in turtles. Gene expression analysis in the European pond turtle (Emys orbicularis) showed that EDC genes are differentially expressed in the skin of the various body sites and that a subset of beta-keratin genes within the EDC as well as those located outside of the EDC are expressed predominantly in the shell. Our findings give strong support to the hypothesis that the evolutionary innovation of the turtle shell involved specific molecular adaptations of epidermal differentiation., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Large phylogenomics data sets require fast tree inference methods, especially for maximum-likelihood (ML) phylogenies. Fast programs exist, but due to inherent heuristics to find optimal trees, it is not clear whether the best tree is found. Thus, there is need for additional approaches that employ different search strategies to find ML trees and that are at the same time as fast as currently available ML programs. We show that a combination of hill-climbing approaches and a stochastic perturbation method can be time-efficiently implemented. If we allow the same CPU time as RAxML and PhyML, then our software IQ-TREE found higher likelihoods between 62.2% and 87.1% of the studied alignments, thus efficiently exploring the tree-space. If we use the IQ-TREE stopping rule, RAxML and PhyML are faster in 75.7% and 47.1% of the DNA alignments and 42.2% and 100% of the protein alignments, respectively. However, the range of obtaining higher likelihoods with IQ-TREE improves to 73.3-97.1%. IQ-TREE is freely available at http://www.cibiv.at/software/iqtree., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

In bilaterians, three orthogonal body axes define the animal form, with distinct anterior-posterior, dorsal-ventral and left-right asymmetries. The key signalling factors are Wnt family proteins for the anterior-posterior axis, Bmp family proteins for the dorsal-ventral axis and Nodal for the left-right axis. Cnidarians, the sister group to bilaterians, are characterized by one oral-aboral body axis, which exhibits a distinct biradiality of unknown molecular nature. Here we analysed the biradial growth pattern in the radially symmetrical cnidarian polyp Hydra, and we report evidence of Nodal in a pre-bilaterian clade. We identified a Nodal-related gene (Ndr) in Hydra magnipapillata, and this gene is essential for setting up an axial asymmetry along the main body axis. This asymmetry defines a lateral signalling centre, inducing a new body axis of a budding polyp orthogonal to the mother polyp's axis. Ndr is expressed exclusively in the lateral bud anlage and induces Pitx, which encodes an evolutionarily conserved transcription factor that functions downstream of Nodal. Reminiscent of its function in vertebrates, Nodal acts downstream of β-Catenin signalling. Our data support an evolutionary scenario in which a 'core-signalling cassette' consisting of β-Catenin, Nodal and Pitx pre-dated the cnidarian-bilaterian split. We presume that this cassette was co-opted for various modes of axial patterning: for example, for lateral branching in cnidarians and left-right patterning in bilaterians.

Background: The availability of many gene alignments with overlapping taxon sets raises the question of which strategy is the best to infer species phylogenies from multiple gene information. Methods and programs abound that use the gene alignment in different ways to reconstruct the species tree. In particular, different methods combine the original data at different points along the way from the underlying sequences to the final tree. Accordingly, they are classified into superalignment, supertree and medium-level approaches. Here, we present a simulation study to compare different methods from each of these three approaches., Results: We observe that superalignment methods usually outperform the other approaches over a wide range of parameters including sparse data and gene-specific evolutionary parameters. In the presence of high incongruency among gene trees, however, other combination methods show better performance than the superalignment approach. Surprisingly, some supertree and medium-level methods exhibit, on average, worse results than a single gene phylogeny with complete taxon information., Conclusions: For some methods, using the reconstructed gene tree as an estimation of the species tree is superior to the combination of incomplete information. Superalignment usually performs best since it is less susceptible to stochastic error. Supertree methods can outperform superalignment in the presence of gene-tree conflict.

Regeneration is widespread throughout the animal kingdom, but our molecular understanding of this process in adult animals remains poorly understood. Wnt/β-catenin signaling plays crucial roles throughout animal life from early development to adulthood. In intact and regenerating planarians, the regulation of Wnt/β-catenin signaling functions to maintain and specify anterior/posterior (A/P) identity. Here, we explore the expression kinetics and RNAi phenotypes for secreted members of the Wnt signaling pathway in the planarian Schmidtea mediterranea. Smed-wnt and sFRP expression during regeneration is surprisingly dynamic and reveals fundamental aspects of planarian biology that have been previously unappreciated. We show that after amputation, a wounding response precedes rapid re-organization of the A/P axis. Furthermore, cells throughout the body plan can mount this response and reassess their new A/P location in the complete absence of stem cells. While initial stages of the amputation response are stem cell independent, tissue remodeling and the integration of a new A/P address with anatomy are stem cell dependent. We also show that WNT5 functions in a reciprocal manner with SLIT to pattern the planarian mediolateral axis, while WNT11-2 patterns the posterior midline. Moreover, we perform an extensive phylogenetic analysis on the Smed-wnt genes using a method that combines and integrates both sequence and structural alignments, enabling us to place all nine genes into Wnt subfamilies for the first time., (Copyright © 2010 Elsevier Inc. All rights reserved.)

Wnt genes and beta-catenin signaling are involved in axial patterning processes in vertebrate embryogenesis in setting up the Spemann-Mangold organizer in amphibian embryos. An organizer with a similar function is present in the hypostome of an adult Hydra polyp. Previously, a Hydra ortholog of Wnt3 (HyWnt3), which is expressed in the hypostome, has been described. Here, ten additional Hydra Wnt genes have been identified. Of these, six (HyWnt1, -7, -9/10a, -9/10c, -11, and -16) are expressed in the adult hypostome. And, as is HyWnt3, these six Wnt genes are also expressed when a new head organizer is formed during head regeneration and bud formation. The kinetics of Wnt gene expressions during head regeneration suggests that a cascade of consecutive Wnt activation accompanies regeneration, and HyWnt3 begins this cascade. Recombinant HyWnt3 protein induced body column tissue to undergo head formation. It also increased the head formation capacity in the head regeneration-deficient mutant strain reg-16 to that of wild-type strains. In addition our data reveal striking similarities in the molecular basis of the organizer in Hydra and axis polarization in chordates (e.g. Spemann's organizer) as well as it's role in regeneration suggesting a conserved function of Wnt signaling in setting up this ancient metazoan signaling center.

Until recently, cave bears were believed to have only inhabited Europe. However, recent morphological evidence suggests that cave bears' geographic range extended as far east as Transbaikalia, Eastern Siberia. These Asian cave bears were morphologically distinct from European cave bears. However, how they related to European lineages remains unclear, stressing the need to assess the phylogenetic and phylogeographic relationship between Asian cave bears and their European relatives. In this work, we address this issue using a 227 base-pair fragment of the mitochondrial control region obtained from nine fossil bone samples from eight sites from the Urals, Caucasus, Altai Mountains, Ukraine and Yana River region in Eastern Siberia. Results of the phylogenetic analyses indicate that (i) the cave bear from the Yana River is most closely related to cave bears from the Caucasus region; (ii) the Caucasus/Yana group of bears is genetically very distinct from both European cave bears and brown bears, suggesting that these bears could represent an independent species; and (iii) the Western European cave bear lineage reached at least temporarily to the Altai Mountains, 7000 km east of their known centre of distribution. These results suggest that the diversity of cave bears was greater than previously believed, and that they could survive in a much wider range of ecological conditions than previously assumed. They also agree with recent studies on other extinct and extant species, such as wolves, hyenas and steppe bison, which have also revealed higher genetic and ecological diversity in Pleistocene populations than previously known.

Background: Anomalous gene trees (AGTs) are gene trees with a topology different from a species tree that are more probable to observe than congruent gene trees. In this paper we propose a rooted triple approach to finding the correct species tree in the presence of AGTs., Results: Based on simulated data we show that our method outperforms the extended majority rule consensus strategy, while still resolving the species tree. Applying both methods to a metazoan data set of 216 genes, we tested whether AGTs substantially interfere with the reconstruction of the metazoan phylogeny., Conclusion: Evidence of AGTs was not found in this data set, suggesting that erroneously reconstructed gene trees are the most significant challenge in the reconstruction of phylogenetic relationships among species with current data. The new method does however rule out the erroneous reconstruction of deep or poorly resolved splits in the presence of lineage sorting.

pIPHULA is the parallel program to estimate the parameters of a realistic model of population growth.

TREE-PUZZLE provides a means to analyze and reconstruct evolutionary relationships and trees based on quartets, i.e., groups of four sequences. Basic Protocol 1 explains how to reconstruct trees based on the maximum-likelihood principle and quartet puzzling. Basic Protocol 2 discusses likelihood mapping, a method to visualize phylogenetic content in a multiple sequence alignment. Basic Protocol 3 explains how to compare tree topologies using different tests.

The dickkopf (dkk) gene family encodes secreted antagonists of Wnt signalling proteins, which have important functions in the control of cell fate, proliferation, and cell polarity during development. Here, we report the isolation, from a regeneration-specific signal peptide screen, of a novel dickkopf gene from the fresh water cnidarian Hydra. Comparative sequence analysis demonstrates that the Wnt antagonistic subfamily Dkk1/Dkk2/Dkk4 and the non-modulating subfamily Dkk3 separated prior to the divergence of cnidarians and bilaterians. In steady-state Hydra, hydkk1/2/4-expression is inversely related to that of hywnt3a. hydkk1/2/4 is an early injury and regeneration responsive gene, and hydkk1/2/4-expressing gland cells are essential for head regeneration in Hydra, although once the head has regenerated they are excluded from it. Activation of Wnt/beta-Catenin signalling leads to the complete downregulation of hydkk1/2/4 transcripts. When overexpressed in Xenopus, HyDkk1/2/4 has similar Wnt-antagonizing activity to the Xenopus gene Dkk1. Based on the corresponding expression patterns of hydkk1/2/4 and neuronal genes, we suggest that the body column of Hydra is a neurogenic environment suppressing Wnt signalling and facilitating neurogenesis.

Summary: IQPNNI is a program to infer maximum-likelihood phylogenetic trees from DNA or protein data with a large number of sequences. We present an improved and MPI-parallel implementation showing very good scaling and speed-up behavior.

The Wnt gene family encodes secreted signalling molecules that control cell fate in animal development and human diseases. Despite its significance, the evolution of this metazoan-specific protein family is unclear. In vertebrates, twelve Wnt subfamilies were defined, of which only six have counterparts in Ecdysozoa (for example, Drosophila and Caenorhabditis). Here, we report the isolation of twelve Wnt genes from the sea anemone Nematostella vectensis, a species representing the basal group within cnidarians. Cnidarians are diploblastic animals and the sister-group to bilaterian metazoans. Phylogenetic analyses of N. vectensis Wnt genes reveal a thus far unpredicted ancestral diversity within the Wnt family. Cnidarians and bilaterians have at least eleven of the twelve known Wnt gene subfamilies in common; five subfamilies appear to be lost in the protostome lineage. Expression patterns of Wnt genes during N. vectensis embryogenesis indicate distinct roles of Wnts in gastrulation, resulting in serial overlapping expression domains along the primary axis of the planula larva. This unexpectedly complex inventory of Wnt family signalling factors evolved in early multi-cellular animals about 650 million years (Myr) ago, predating the Cambrian explosion by at least 100 Myr (refs 5, 8). It emphasizes the crucial function of Wnt genes in the diversification of eumetazoan body plans.

Dobrava virus (DOBV) occurs in two different rodent species, Apodemus flavicollis (DOBV-Af) and A. agrarius (DOBV-Aa). We sequenced the S and M genomic segments from sympatric DOBV-Af and DOBV-Aa strains which fell into two distinct genetic lineages. Molecular phylogenetic analyses gave evidence for genetic reassortment between S and M segments of DOBV-Af and DOBV-Aa and indicated homologous recombination events in DOBV evolution. DOBV-Af and DOBV-Aa are distinct but also subject to genetic exchanges that affect their evolutionary trajectories.

Summary: TREE-PUZZLE is a program package for quartet-based maximum-likelihood phylogenetic analysis (formerly PUZZLE, Strimmer and von Haeseler, Mol. Biol. Evol., 13, 964-969, 1996) that provides methods for reconstruction, comparison, and testing of trees and models on DNA as well as protein sequences. To reduce waiting time for larger datasets the tree reconstruction part of the software has been parallelized using message passing that runs on clusters of workstations as well as parallel computers., Availability: http://www.tree-puzzle.de. The program is written in ANSI C. TREE-PUZZLE can be run on UNIX, Windows and Mac systems, including Mac OS X. To run the parallel version of PUZZLE, a Message Passing Interface (MPI) library has to be installed on the system. Free MPI implementations are available on the Web (cf. http://www.lam-mpi.org/mpi/implementations/).