Contentious relationships in phylogenomic studies can be driven by a handful of genes

Genome Sequences Reveal Cryptic Speciation in the Human Pathogen Histoplasma capsulatum

Histoplasma capsulatum is a pathogenic fungus that causes life-threatening lung infections. About 500,000 people are exposed to H. capsulatum each year in the United States, and over 60% of the U.S. population has been exposed to the fungus at some point in their life. We performed genome-wide population genetics and phylogenetic analyses with 30 Histoplasma isolates representing four recognized areas where histoplasmosis is endemic and show that the Histoplasma genus is composed of at least four species that are genetically isolated and rarely interbreed. Therefore, we propose a taxonomic rearrangement of the genus.IMPORTANCE The evolutionary processes that give rise to new pathogen lineages are critical to our understanding of how they adapt to new environments and how frequently they exchange genes with each other. The fungal pathogen Histoplasma capsulatum provides opportunities to precisely test hypotheses about the origin of new genetic variation. We find that H. capsulatum is composed of at least four different cryptic species that differ genetically and also in virulence. These results have implications for the epidemiology of histoplasmosis because not all Histoplasma species are equivalent in their geographic range and ability to cause disease.

Expanded Taxonomic Sampling Coupled with Gene Genealogy Interrogation Provides Unambiguous Resolution for the Evolutionary Root of Angiosperms

The branching order of major angiosperm lineages is a challenging phylogenetic question that has received substantial attention in recent years. Two main competing hypotheses place the New Caledonian Amborella as either sister to all other extant angiosperms (Amborella-sister) or to the water lilies (Amborella + Nymphaeales). Here, we revisit this question by expanding a transcriptomic data set of 310 genes previously assembled to include data from seven species comprising two major lineages of flowering plants that were poorly represented or missing from the original study. We also applied gene genealogy interrogation, a recent approach based on constrained tree searches in combination with topology tests, to account for gene tree estimation error and its downstream effects in coalescent analyses. In addition to gene genealogy interrogation, we conducted a large number of multilocus analyses, including concatenation and coalescent approaches (using both unconstrained and constrained gene trees), and based on different data sets (original and expanded) and data types (nucleotide and amino acid sequences). We show that the majority of gene trees favor Amborella-sister topology, and all multilocus analyses conducted (concatenation and coalescent) provide overwhelming support for this hypothesis regardless of data type. Beyond resolving the evolutionary root of angiosperms with confidence, our results highlight the importance of both broadening taxonomic sampling in phylogenomics and addressing the effects of gene tree error in summary coalescent inferences.

The maternal-zygotic transition and zygotic activation of the Mnemiopsis leidyi genome occurs within the first three cleavage cycles

The maternal-zygotic transition (MZT) describes the developmental reprogramming of gene expression marked by the degradation of maternally supplied gene products and activation of the zygotic genome. While the timing and duration of the MZT vary among taxa, little is known about early-stage transcriptional dynamics in the non-bilaterian phylum Ctenophora. We sought to better understand the extent of maternal mRNA loading and subsequent differential transcript abundance during the earliest stages of development by performing comprehensive RNA-sequencing-based analyses of mRNA abundance in single- and eight-cell stage embryos in the lobate ctenophore Mnemiopsis leidyi. We found 1,908 contigs with significant differential abundance between single- and eight-cell stages, of which 1,208 contigs were more abundant at the single-cell stage and 700 contigs were more abundant at the eight-cell stage. Of the differentially abundant contigs, 267 were exclusively present in the eight-cell samples, providing strong evidence that both the MZT and zygotic genome activation (ZGA) have commenced by the eight-cell stage. Many highly abundant transcripts encode genes involved in molecular mechanisms critical to the MZT, such as maternal transcript degradation, serine/threonine kinase activity, and chromatin remodeling. Our results suggest that chromosomal restructuring, which is critical to ZGA and the initiation of transcriptional regulation necessary for normal development, begins by the third cleavage within 1.5 hr post-fertilization in M. leidyi.

Maximizing Power in Phylogenetics and Phylogenomics: A Perspective Illuminated by Fungal Big Data

Since its original inception over 150 years ago by Darwin, we have made tremendous progress toward the reconstruction of the Tree of Life. In particular, the transition from analyzing datasets comprised of small numbers of loci to those comprised of hundreds of loci, if not entire genomes, has aided in resolving some of the most vexing of evolutionary problems while giving us a new perspective on biodiversity. Correspondingly, phylogenetic trees have taken a central role in fields that span ecology, conservation, and medicine. However, the rise of big data has also presented phylogenomicists with a new set of challenges to experimental design, quantitative analyses, and computation. The sequencing of a number of very first genomes presented significant challenges to phylogenetic inference, leading fungal phylogenomicists to begin addressing pitfalls and postulating solutions to the issues that arise from genome-scale analyses relevant to any lineage across the Tree of Life. Here we highlight insights from fungal phylogenomics for topics including systematics and species delimitation, ecological and phenotypic diversification, and biogeography while providing an overview of progress made on the reconstruction of the fungal Tree of Life. Finally, we provide a review of considerations to phylogenomic experimental design for robust tree inference. We hope that this special issue of Advances in Genetics not only excites the continued progress of fungal evolutionary biology but also motivates the interdisciplinary development of new theory and methods designed to maximize the power of genomic scale data in phylogenetic analyses.

Ctenophore relationships and their placement as the sister group to all other animals

Ctenophora, compromising approximately 200 described species, is an important lineage for understanding metazoan evolution and is of great ecological and economic importance. Ctenophore diversity includes species with unique colloblasts used for prey capture, smooth and striated muscles, benthic and pelagic lifestyles, and locomotion with ciliated paddles or muscular propulsion. However, ancestral states of traits are debated and relationships among many lineages are unresolved. Here, using 27 newly sequenced ctenophore transcriptomes, publicly available data, and methods to control systematic error we establish the placement of Ctenophora as the sister group to all other animals and refine phylogenetic relationships within ctenophores. Molecular clock analyses suggest modern ctenophore diversity originated approximately 350MYA ± 88 MY, conflicting with previous hypotheses of approximately 65 MYA. We recover Euplokamis dunlapae, a species with striated muscles, as the sister lineage to other sampled ctenophores. Ancestral state reconstruction shows the most recent common ancestor of extant ctenophores was pelagic, possessed tentacles, was bioluminescent, and did not have separate sexes. Our results imply at least two transitions from a pelagic to a benthic lifestyle within Ctenophora, suggesting such transitions were more common in animal diversification than appreciated.

Early-branching euteleost relationships: areas of congruence between concatenation and coalescent model inferences

Phylogenetic inference based on evidence from DNA sequences has led to significant strides in the development of a stable and robustly supported framework for the vertebrate tree of life. To date, the bulk of those advances have relied on sequence data from a small number of genome regions that have proven unable to produce satisfactory answers to consistently recalcitrant phylogenetic questions. Here, we re-examine phylogenetic relationships among early-branching euteleostean fish lineages classically grouped in the Protacanthopterygii using DNA sequence data surrounding ultraconserved elements. We report and examine a dataset of thirty-four OTUs with 17,957 aligned characters from fifty-three nuclear loci. Phylogenetic analysis is conducted in concatenated, joint gene trees and species tree estimation and summary coalescent frameworks. All analytical frameworks yield supporting evidence for existing hypotheses of relationship for the placement of Lepidogalaxias salamandroides, monophyly of the Stomiatii and the presence of an esociform + salmonid clade. Lepidogalaxias salamandroides and the Esociformes + Salmoniformes are successive sister lineages to all other euteleosts in the majority of analyses. The concatenated and joint gene trees and species tree analysis types produce high support values for this arrangement. However, inter-relationships of Argentiniformes, Stomiatii and Neoteleostei remain uncertain as they varied by analysis type while receiving strong and contradictory indices of support. Topological differences between analysis types are also apparent within the otomorph and the percomorph taxa in the data set. Our results identify concordant areas with strong support for relationships within and between early-branching euteleost lineages but they also reveal limitations in the ability of larger datasets to conclusively resolve other aspects of that phylogeny.

Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements

The phylogeny of eutherian mammals contains some of the most recalcitrant nodes in the tetrapod tree of life. We combined comprehensive taxon and character sampling to explore three of the most debated interordinal relationships among placental mammals. We performed in silico extraction of ultraconserved element loci from 72 published genomes and invitro enrichment and sequencing of ultraconserved elements from 28 additional mammals, resulting in alignments of 3,787 loci. We analyzed these data using concatenated and multispecies coalescent phylogenetic approaches, topological tests, and exploration of support among individual loci to identify the root of Eutheria and the sister groups of tree shrews (Scandentia) and horses (Perissodactyla). Individual loci provided weak, but often consistent support for topological hypotheses. Although many gene trees lacked accepted species-tree relationships, summary coalescent topologies were largely consistent with inferences from concatenation. At the root of Eutheria, we identified consistent support for a sister relationship between Xenarthra and Afrotheria (i.e., Atlantogenata). At the other nodes of interest, support was less consistent. We suggest Scandentia is the sister of Primatomorpha (Euarchonta), but we failed to reject a sister relationship between Scandentia and Glires. Similarly, we suggest Perissodactyla is sister to Cetartiodactyla (Euungulata), but a sister relationship between Perissodactyla and Chiroptera remains plausible.

Evolutionary origin of synapses and neurons - Bridging the gap

The evolutionary origin of synapses and neurons is an enigmatic subject that inspires much debate. Non-bilaterian metazoans, both with and without neurons and their closest relatives already contain many components of the molecular toolkits for synapse functions. The origin of these components and their assembly into ancient synaptic signaling machineries are particularly important in light of recent findings on the phylogeny of non-bilaterian metazoans. The evolution of synapses and neurons are often discussed only from a metazoan perspective leaving a considerable gap in our understanding. By taking an integrative approach we highlight the need to consider different, but extremely relevant phyla and to include the closest unicellular relatives of metazoans, the ichthyosporeans, filastereans and choanoflagellates, to fully understand the evolutionary origin of synapses and neurons. This approach allows for a detailed understanding of when and how the first pre- and postsynaptic signaling machineries evolved.

Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements

The phylogeny of eutherian mammals contains some of the most recalcitrant nodes in the tetrapod tree of life. We combined comprehensive taxon and character sampling to explore three of the most debated interordinal relationships among placental mammals. We performed in silico extraction of ultraconserved element loci from 72 published genomes and invitro enrichment and sequencing of ultraconserved elements from 28 additional mammals, resulting in alignments of 3,787 loci. We analyzed these data using concatenated and multispecies coalescent phylogenetic approaches, topological tests, and exploration of support among individual loci to identify the root of Eutheria and the sister groups of tree shrews (Scandentia) and horses (Perissodactyla). Individual loci provided weak, but often consistent support for topological hypotheses. Although many gene trees lacked accepted species-tree relationships, summary coalescent topologies were largely consistent with inferences from concatenation. At the root of Eutheria, we identified consistent support for a sister relationship between Xenarthra and Afrotheria (i.e., Atlantogenata). At the other nodes of interest, support was less consistent. We suggest Scandentia is the sister of Primatomorpha (Euarchonta), but we failed to reject a sister relationship between Scandentia and Glires. Similarly, we suggest Perissodactyla is sister to Cetartiodactyla (Euungulata), but a sister relationship between Perissodactyla and Chiroptera remains plausible.

Targeted NGS for species level phylogenomics: "made to measure" or "one size fits all"?

Targeted high-throughput sequencing using hybrid-enrichment offers a promising source of data for inferring multiple, meaningfully resolved, independent gene trees suitable to address challenging phylogenetic problems in species complexes and rapid radiations. The targets in question can either be adopted directly from more or less universal tools, or custom made for particular clades at considerably greater effort. We applied custom made scripts to select sets of homologous sequence markers from transcriptome and WGS data for use in the flowering plant genus Erica (Ericaceae). We compared the resulting targets to those that would be selected both using different available tools (Hyb-Seq; MarkerMiner), and when optimising for broader clades of more distantly related taxa (Ericales; eudicots). Approaches comparing more divergent genomes (including MarkerMiner, irrespective of input data) delivered fewer and shorter potential markers than those targeted for Erica. The latter may nevertheless be effective for sequence capture across the wider family Ericaceae. We tested the targets delivered by our scripts by obtaining an empirical dataset. The resulting sequence variation was lower than that of standard nuclear ribosomal markers (that in Erica fail to deliver a well resolved gene tree), confirming the importance of maximising the lengths of individual markers. We conclude that rather than searching for "one size fits all" universal markers, we should improve and make more accessible the tools necessary for developing "made to measure" ones.

Dating early animal evolution using phylogenomic data

Information about the geological timeframe during which animals radiated into their major subclades is crucial to understanding early animal ecology and evolution. Unfortunately, the pre-Cambrian fossil record is sparse and its interpretation controversial. Relaxed molecular-clock methods provide an alternative means of estimating the timing of cladogenesis deep in the metazoan tree of life. So far, thorough molecular clock studies focusing specifically on Metazoa as a whole have been based on relatively small datasets or incomplete representation of the main non-bilaterian lineages (such as sponges and ctenophores), which are fundamental for understanding early metazoan evolution. Here, we use a previously published phylogenomic dataset that includes a fair sampling of all relevant groups to estimate the timing of early animal evolution with Bayesian relaxed-clock methods. According to our results, all non-bilaterian phyla, as well as total-group Bilateria, evolved in an ancient radiation during a geologically relatively short time span, before the onset of long-term global glaciations (“Snowball Earth”; ~720–635 Ma). Importantly, this result appears robust to alterations of a number of important analytical variables, such as models of among-lineage rate variation and sets of fossil calibrations used.

The Understanding of the Metazoan Skeletal System, Based on the Initial Discoveries with Siliceous and Calcareous Sponges

Initiated by studies on the mechanism of formation of the skeletons of the evolutionary oldest still extant multicellular animals, the sponges (phylum Porifera) have provided new insights into the mechanism of formation of the Ca-phosphate/hydroxyapatite skeleton of vertebrate bone. Studies on the formation of the biomineral skeleton of sponges revealed that both the formation of the inorganic siliceous skeletons (sponges of the class of Hexactinellida and Demospongiae) and of the calcareous skeletons (class of Calcarea) is mediated by enzymes (silicatein: polymerization of biosilica; and carbonic anhydrase: deposition of Ca-carbonate). Detailed studies of the initial mineralization steps in human bone-forming cells showed that this process is also controlled by enzymes, starting with the deposition of Ca-carbonate bio-seeds, mediated by carbonic anhydrases-II and -IX, followed by non-enzymatic transformation of the formed amorphous Ca-carbonate deposits into amorphous Ca-phosphate and finally hydroxyapatite crystals. The required phosphate is provided by enzymatic (alkaline phosphatase-mediated) degradation of an inorganic polymer, polyphosphate (polyP), which also acts as a donor for chemically useful energy in this process. These new discoveries allow the development of novel biomimetic strategies for treatment of bone diseases and defects.

Lokiarchaea are close relatives of Euryarchaeota, not bridging the gap between prokaryotes and eukaryotes

The eocyte hypothesis, in which Eukarya emerged from within Archaea, has been boosted by the description of a new candidate archaeal phylum, "Lokiarchaeota", from metagenomic data. Eukarya branch within Lokiarchaeota in a tree reconstructed from the concatenation of 36 universal proteins. However, individual phylogenies revealed that lokiarchaeal proteins sequences have different evolutionary histories. The individual markers phylogenies revealed at least two subsets of proteins, either supporting the Woese or the Eocyte tree of life. Strikingly, removal of a single protein, the elongation factor EF2, is sufficient to break the Eukaryotes-Lokiarchaea affiliation. Our analysis suggests that the three lokiarchaeal EF2 proteins have a chimeric organization that could be due to contamination and/or homologous recombination with patches of eukaryotic sequences. A robust phylogenetic analysis of RNA polymerases with a new dataset indicates that Lokiarchaeota and related phyla of the Asgard superphylum are sister group to Euryarchaeota, not to Eukarya, and supports the monophyly of Archaea with their rooting in the branch leading to Thaumarchaeota.

Targeted NGS for species level phylogenomics: "made to measure" or "one size fits all"?

Targeted high-throughput sequencing using hybrid-enrichment offers a promising source of data for inferring multiple, meaningfully resolved, independent gene trees suitable to address challenging phylogenetic problems in species complexes and rapid radiations. The targets in question can either be adopted directly from more or less universal tools, or custom made for particular clades at considerably greater effort. We applied custom made scripts to select sets of homologous sequence markers from transcriptome and WGS data for use in the flowering plant genus Erica (Ericaceae). We compared the resulting targets to those that would be selected both using different available tools (Hyb-Seq; MarkerMiner), and when optimising for broader clades of more distantly related taxa (Ericales; eudicots). Approaches comparing more divergent genomes (including MarkerMiner, irrespective of input data) delivered fewer and shorter potential markers than those targeted for Erica. The latter may nevertheless be effective for sequence capture across the wider family Ericaceae. We tested the targets delivered by our scripts by obtaining an empirical dataset. The resulting sequence variation was lower than that of standard nuclear ribosomal markers (that in Erica fail to deliver a well resolved gene tree), confirming the importance of maximising the lengths of individual markers. We conclude that rather than searching for "one size fits all" universal markers, we should improve and make more accessible the tools necessary for developing "made to measure" ones.