Excluding loci with substitution saturation improves inferences from phylogenomic data

PhyIN: trimming alignments by phylogenetic incompatibilities among neighbouring sites

In phylogenomics, regions of low alignment reliability and high noise are typically trimmed from multiple sequence alignments before they are used in phylogenetic inference. I introduce a new trimming tool, PhyIN, which deletes regions in which a large proportion of sites (characters) have conflicting phylogenetic signal. It does not require inference of a phylogenetic tree, as it finds neighbouring characters that cannot agree on any possible tree. In phylogenomic data of ultraconserved elements (UCE), PhyIN effectively finds the boundaries between chaotic (conflicted) and orderly regions of alignments with data for only a single locus. Its ability to work on individual loci allows it to preserve discord between gene trees and species trees.

Phylogenomic analyses of Blattodea combining traditional methods, incremental tree-building, and quality-aware support

Despite the many advances of the genomic era, there is a persistent problem in assessing the uncertainty of phylogenomic hypotheses. We see this in the recent history of phylogenetics for cockroaches and termites (Blattodea), where huge advances have been made, but there are still major inconsistencies between studies. To address this, we present a phylogenetic analysis of Blattodea that emphasizes identification and quantification of uncertainty. We analyze 1183 gene domains using three methods (multi-species coalescent inference, concatenation, and a supermatrix-supertree hybrid approach) and assess support for controversial relationships while considering data quality. The hybrid approach-here dubbed "tiered phylogenetic inference"-incorporates information about data quality into an incremental tree building framework. Leveraging this method, we are able to identify cases of low or misleading support that would not be possible otherwise, and explore them more thoroughly with follow-up tests. In particular, quality annotations pointed towards nodes with high bootstrap support that later turned out to have large ambiguities, sometimes resulting from low-quality data. We also clarify issues related to some recalcitrant nodes: Anaplectidae's placement lacks unbiased signal, Ectobiidae s.s. and Anaplectoideini need greater taxon sampling, the deepest relationships among most Blaberidae lack signal. As a result, several previous phylogenetic uncertainties are now closer to being resolved (e.g., African and Malagasy "Rhabdoblatta" spp. are the sister to all other Blaberidae, and Oxyhaloinae is sister to the remaining Blaberidae). Overall, we argue for more approaches to quantifying support that take data quality into account to uncover the nature of recalcitrant nodes.

Dating in the Dark: Elevated Substitution Rates in Cave Cockroaches (Blattodea: Nocticolidae) Have Negative Impacts on Molecular Date Estimates

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analyzed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all 13 major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales.

Reliable estimation of tree branch lengths using deep neural networks

A phylogenetic tree represents hypothesized evolutionary history for a set of taxa. Besides the branching patterns (i.e., tree topology), phylogenies contain information about the evolutionary distances (i.e. branch lengths) between all taxa in the tree, which include extant taxa (external nodes) and their last common ancestors (internal nodes). During phylogenetic tree inference, the branch lengths are typically co-estimated along with other phylogenetic parameters during tree topology space exploration. There are well-known regions of the branch length parameter space where accurate estimation of phylogenetic trees is especially difficult. Several novel studies have recently demonstrated that machine learning approaches have the potential to help solve phylogenetic problems with greater accuracy and computational efficiency. In this study, as a proof of concept, we sought to explore the possibility of machine learning models to predict branch lengths. To that end, we designed several deep learning frameworks to estimate branch lengths on fixed tree topologies from multiple sequence alignments or its representations. Our results show that deep learning methods can exhibit superior performance in some difficult regions of branch length parameter space. For example, in contrast to maximum likelihood inference, which is typically used for estimating branch lengths, deep learning methods are more efficient and accurate. In general, we find that our neural networks achieve similar accuracy to a Bayesian approach and are the best-performing methods when inferring long branches that are associated with distantly related taxa. Together, our findings represent a next step toward accurate, fast, and reliable phylogenetic inference with machine learning approaches.

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

Accurate Recapitulation of Chikungunya Virus Complete Coding Sequence Phylogeny Using Variable Genome Regions for Genomic Surveillance

Chikungunya virus (CHIKV) is transmitted by mosquito bites and causes chikungunya fever (CHIKF). CHIKV has a single-stranded RNA genome and belongs to a single serotype with three genotypes. The Asian lineage has recently emerged in the Western Hemisphere, likely due to travel-associated introduction. Genetic variation accumulates in the CHIKV genome as the virus replicates, creating new lineages. Whole genome sequencing is ideal for studying virus evolution and spread but is expensive and complex. This study investigated whether specific, highly variable regions of the CHIKV genome could recapitulate the phylogeny obtained with a complete coding sequence (CDS). Our results revealed that concatenated highly variable regions accurately reconstructed CHIKV phylogeny, exhibiting statistically indistinguishable branch lengths and tree confidence compared to CDS. In addition, these regions adequately inferred the evolutionary relationships among CHIKV isolates from the American outbreak with similar results to the CDS. This finding suggests that highly variable regions can effectively capture the evolutionary relationships among CHIKV isolates, offering a simpler approach for future studies. This approach could be particularly valuable for large-scale surveillance efforts.

Complexity of avian evolution revealed by family-level genomes

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Evolution of bioluminescence in Anthozoa with emphasis on Octocorallia

Bioluminescence is a widespread phenomenon that has evolved multiple times across the tree of life, converging among diverse fauna and habitat types. The ubiquity of bioluminescence, particularly in marine environments where it is commonly used for communication and defense, highlights the adaptive value of this trait, though the evolutionary origins and timing of emergence remain elusive for a majority of luminous organisms. Anthozoan cnidarians are a diverse group of animals with numerous bioluminescent species found throughout the world's oceans, from shallow waters to the light-limited deep sea where bioluminescence is particularly prominent. This study documents the presence of bioluminescent Anthozoa across depth and explores the diversity and evolutionary origins of bioluminescence among Octocorallia-a major anthozoan group of marine luminous organisms. Using a phylogenomic approach and ancestral state reconstruction, we provide evidence for a single origin of bioluminescence in Octocorallia and infer the age of occurrence to around the Cambrian era, approximately 540 Ma-setting a new record for the earliest timing of emergence of bioluminescence in the marine environment. Our results further suggest this trait was largely maintained in descendants of a deep-water ancestor and bioluminescent capabilities may have facilitated anthozoan diversification in the deep sea.

ClockstaRX: Testing Molecular Clock Hypotheses With Genomic Data

Phylogenomic data provide valuable opportunities for studying evolutionary rates and timescales. These analyses require theoretical and statistical tools based on molecular clocks. We present ClockstaRX, a flexible platform for exploring and testing evolutionary rate signals in phylogenomic data. Here, information about evolutionary rates in branches across gene trees is placed in Euclidean space, allowing data transformation, visualization, and hypothesis testing. ClockstaRX implements formal tests for identifying groups of loci and branches that make a large contribution to patterns of rate variation. This information can then be used to test for drivers of genomic evolutionary rates or to inform models for molecular dating. Drawing on the results of a simulation study, we recommend forms of data exploration and filtering that might be useful prior to molecular-clock analyses.

300 million years apart: the extreme case of macromorphological skeletal convergence between deltocyathids and a turbinoliid coral (Anthozoa, Scleractinia)

The integration of morphological and molecular lines of evidence has enabled the family Deltocyathidae to be erected to accommodate Deltocyathus species that were previously ascribed to the family Caryophylliidae. However, although displaying the same morphological characteristics as other species of Deltocyathus , molecular data suggested that D. magnificus was phylogenetically distant from Deltocyathidae, falling within the family Turbinoliidae instead. To elucidate the enigmatic evolutionary history of this species and skeletal microstructural features, the phylogenetic relationships of Deltocyathidae and Turbinoliidae were investigated using nuclear ultraconserved and exon loci and complete mitochondrial genomes. Both nuclear and mitochondrial phylogenomic reconstructions confirmed the position of D. magnificus within turbinolids. Furthermore, a novel mitochondrial gene order was uncovered for Deltocyathidae species. This gene order was not present in Turbinoliidae or in D. magnificus that both have the scleractinian canonical gene order, further indicating the taxonomic utility of mitochondrial gene order. D. magnificus is therefore formally moved to the family Turbinoliidae and accommodated in a new genus (Dennantotrochus Kitahara, Vaga & Stolarski, gen. nov.). Surprisingly, turbinolids and deltocyathids do not differ in microstructural organisation of the skeleton that consists of densely packed, individualised rapid accretion deposits and thickening deposits composed of fibres perpendicular to the skeleton surface. Therefore, although both families are clearly evolutionarily divergent, macromorphological features indicate a case of skeletal convergence while these may still share conservative biomineralisation mechanisms. ZooBank: urn:lsid:zoobank.org:pub:5F1C0E25-3CC6-4D1F-B1F0-CD9D0014678E.

The Neotropical endemic liverwort subfamily Micropterygioideae had circum-Antarctic links to the rest of the Lepidoziaceae during the early Cretaceous

Lepidoziaceae are the third-largest family of liverworts, with about 860 species distributed on all continents. The evolutionary history of this family has not been satisfactorily resolved, with taxa such as Micropterygioideae yet to be included in phylogenetic analyses. We inferred a dated phylogeny of Lepidoziaceae using a data set consisting of 13 genetic markers, sampled from 147 species. Based on our phylogenetic estimate, we used statistical dispersal-vicariance analysis to reconstruct the biogeographic history of the family. We inferred a crown age of 197 Ma (95% credible interval 157-240 Ma) for the family in the Australian region, with most major lineages also originating in the same region. Micropterygioideae are placed as the sister group to Lembidioideae, with these two lineages diverging from each other about 132 Ma in the South American-Australian region. With South America and Australia being connected through Antarctica at the time, our results suggest a circum-Antarctic link between Micropterygioideae and the rest of the family. Crown Micropterygioideae were inferred to have arisen 45 Ma in South America before the continent separated from Antarctica. Extinction from southern temperate regions might explain the present-day restriction of Micropterygioideae to the Neotropics. Our study reveals the influence of past geological events, such as continental drift, on the evolution and distribution of a widespread and diverse family of liverworts.

The Chalcidoidea bush of life: evolutionary history of a massive radiation of minute wasps

Chalcidoidea are mostly parasitoid wasps that include as many as 500 000 estimated species. Capturing phylogenetic signal from such a massive radiation can be daunting. Chalcidoidea is an excellent example of a hyperdiverse group that has remained recalcitrant to phylogenetic resolution. We combined 1007 exons obtained with Anchored Hybrid Enrichment with 1048 ultra-conserved elements (UCEs) for 433 taxa including all extant families, >95% of all subfamilies, and 356 genera chosen to represent the vast diversity of the superfamily. Going back and forth between the molecular results and our collective knowledge of morphology and biology, we detected bias in the analyses that was driven by the saturation of nucleotide data. Our final results are based on a concatenated analysis of the least saturated exons and UCE datasets (2054 loci, 284 106 sites). Our analyses support an expected sister relationship with Mymarommatoidea. Seven previously recognized families were not monophyletic, so support for a new classification is discussed. Natural history in some cases would appear to be more informative than morphology, as illustrated by the elucidation of a clade of plant gall associates and a clade of taxa with planidial first-instar larvae. The phylogeny suggests a transition from smaller soft-bodied wasps to larger and more heavily sclerotized wasps, with egg parasitism as potentially ancestral for the entire superfamily. Deep divergences in Chalcidoidea coincide with an increase in insect families in the fossil record, and an early shift to phytophagy corresponds with the beginning of the "Angiosperm Terrestrial Revolution". Our dating analyses suggest a middle Jurassic origin of 174 Ma (167.3-180.5 Ma) and a crown age of 162.2 Ma (153.9-169.8 Ma) for Chalcidoidea. During the Cretaceous, Chalcidoidea may have undergone a rapid radiation in southern Gondwana with subsequent dispersals to the Northern Hemisphere. This scenario is discussed with regard to knowledge about the host taxa of chalcid wasps, their fossil record and Earth's palaeogeographic history.

A three-genome ultraconserved element phylogeny of cryptophytes

Cryptophytes are single celled protists found in all aquatic environments. They are composed of a heterotrophic genus, Goniomonas, and a largely autotrophic group comprising many genera. Cryptophytes evolved through secondary endosymbiosis between a host eukaryotic heterotroph and a symbiont red alga. This merger resulted in a four-genome system that includes the nuclear and mitochondrial genomes from the host and a second nuclear genome (nucleomorph) and plastid genome inherited from the symbiont. Here, we make use of different genomes (with potentially distinct evolutionary histories) to perform a phylogenomic study of the early history of cryptophytes. Using ultraconserved elements from the host nuclear genome and symbiont nucleomorph and plastid genomes, we produce a three-genome phylogeny of 91 strains of cryptophytes. Our phylogenetic analyses find that that there are three major cryptophyte clades: Clade 1 comprises Chroomonas and Hemiselmis species, Clade 2, a taxonomically rich clade, comprises at least twelve genera, and Clade 3, comprises the heterotrophic Goniomonas species. Each of these major clades include both freshwater and marine species, but subclades within these clades differ in degrees of niche conservatism. Finally, we discuss priorities for taxonomic revision to Cryptophyceae based on previous studies and in light of these phylogenomic analyses.

Mito-nuclear discordance within Anthozoa, with notes on unique properties of their mitochondrial genomes

Whole mitochondrial genomes are often used in phylogenetic reconstruction. However, discordant patterns in species relationships between mitochondrial and nuclear phylogenies are commonly observed. Within Anthozoa (Phylum Cnidaria), mitochondrial (mt)-nuclear discordance has not yet been examined using a large and comparable dataset. Here, we used data obtained from target-capture enrichment sequencing to assemble and annotate mt genomes and reconstruct phylogenies for comparisons to phylogenies inferred from hundreds of nuclear loci obtained from the same samples. The datasets comprised 108 hexacorals and 94 octocorals representing all orders and > 50% of extant families. Results indicated rampant discordance between datasets at every taxonomic level. This discordance is not attributable to substitution saturation, but rather likely caused by introgressive hybridization and unique properties of mt genomes, including slow rates of evolution driven by strong purifying selection and substitution rate variation. Strong purifying selection across the mt genomes caution their use in analyses that rely on assumptions of neutrality. Furthermore, unique properties of the mt genomes were noted, including genome rearrangements and the presence of nad5 introns. Specifically, we note the presence of the homing endonuclease in ceriantharians. This large dataset of mitochondrial genomes further demonstrates the utility of off-target reads generated from target-capture data for mt genome assembly and adds to the growing knowledge of anthozoan evolution.

Using PhyloSuite for molecular phylogeny and tree-based analyses

Phylogenetic analysis has entered the genomics (multilocus) era. For less experienced researchers, conquering the large number of software programs required for a multilocus-based phylogenetic reconstruction can be somewhat daunting and time-consuming. PhyloSuite, a software with a user-friendly GUI, was designed to make this process more accessible by integrating multiple software programs needed for multilocus and single-gene phylogenies and further streamlining the whole process. In this protocol, we aim to explain how to conduct each step of the phylogenetic pipeline and tree-based analyses in PhyloSuite. We also present a new version of PhyloSuite (v1.2.3), wherein we fixed some bugs, made some optimizations, and introduced some new functions, including a number of tree-based analyses, such as signal-to-noise calculation, saturation analysis, spurious species identification, and etc. The step-by-step protocol includes background information (i.e., what the step does), reasons (i.e., why do the step), and operations (i.e., how to do it). This protocol will help researchers quick-start their way through the multilocus phylogenetic analysis, especially those interested in conducting organelle-based analyses.

Placing human gene families into their evolutionary context

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Comparative genomics of the Western Hemisphere soft tick-borne relapsing fever borreliae highlights extensive plasmid diversity

BACKGROUND

Tick-borne relapsing fever (TBRF) is a globally prevalent, yet under-studied vector-borne disease transmitted by soft and hard bodied ticks. While soft TBRF (sTBRF) spirochetes have been described for over a century, our understanding of the molecular mechanisms facilitating vector and host adaptation is poorly understood. This is due to the complexity of their small (~ 1.5 Mb) but fragmented genomes that typically consist of a linear chromosome and both linear and circular plasmids. A majority of sTBRF spirochete genomes' plasmid sequences are either missing or are deposited as unassembled sequences. Consequently, our goal was to generate complete, plasmid-resolved genomes for a comparative analysis of sTBRF species of the Western Hemisphere.

RESULTS

Utilizing a Borrelia specific pipeline, genomes of sTBRF spirochetes from the Western Hemisphere were sequenced and assembled using a combination of short- and long-read sequencing technologies. Included in the analysis were the two recently isolated species from Central and South America, Borrelia puertoricensis n. sp. and Borrelia venezuelensis, respectively. Plasmid analyses identified diverse sequences that clustered plasmids into 30 families; however, only three families were conserved and syntenic across all species. We also compared two species, B. venezuelensis and Borrelia turicatae, which were isolated ~ 6,800 km apart and from different tick vector species but were previously reported to be genetically similar.

CONCLUSIONS

To truly understand the biological differences observed between species of TBRF spirochetes, complete chromosome and plasmid sequences are needed. This comparative genomic analysis highlights high chromosomal synteny across the species yet diverse plasmid composition. This was particularly true for B. turicatae and B. venezuelensis, which had high average nucleotide identity yet extensive plasmid diversity. These findings are foundational for future endeavors to evaluate the role of plasmids in vector and host adaptation.

Systematic position of the Clicking Frog (Kassinula Laurent, 1940), the problem of chimeric sequences and the revised classification of the family Hyperoliidae

The systematics of the African frog family Hyperoliidae has undergone turbulent changes in last decades. Representatives of several genera have not been genetically investigated or with only limited data, and their phylogenetic positions are thus still not reliably known. This is the case of the De Witte's Clicking Frog (Kassinula wittei) which belongs to a monotypic genus. This miniature frog occurs in a poorly studied region, southeastern Democratic Republic of the Congo, northern Zambia, Angola. So far it is not settled whether this genus belongs to the subfamily Kassininae as a relative of the genus Kassina, or to the subfamily Hyperoliinae as a relative of the genus Afrixalus. Here we present for the first time a multilocus phylogenetic reconstruction (using five nuclear and one mitochondrial marker) of the family Hyperoliidae, including Kassinula. We demonstrate with high confidence that Kassinula is a member of Hyperoliinae belonging to a clade also containing Afrixalus (sub-Saharan Africa), Heterixalus (Madagascar) and Tachycnemis (Seychelles). We find that Kassinula represents a divergent lineage (17-25 Mya), which supports its separate genus-level status, but its exact systematic position remains uncertain. We propose to name the clade to which the above four genera belong as the tribe Tachycnemini Channing, 1989. A new taxonomy of the family Hyperoliidae was recently proposed by Dubois et al. (2021: Megataxa 5, 1-738). We demonstrate here that the new taxonomy was based on a partially erroneous phylogenetic reconstruction resulting from a supermatrix analysis of chimeric DNA sequences combining data from two families, Hyperoliidae and Arthroleptidae (the case of Cryptothylax). We therefore correct the erroneous part and propose a new, revised suprageneric taxonomy of the family Hyperoliidae. We also emphasize the importance of inspecting individual genetic markers before their concatenation or coalescent-based tree reconstructions to avoid analyses of chimeric DNA sequences producing incorrect phylogenetic reconstructions. Especially when phylogenetic reconstructions are used to propose taxonomies and systematic classifications.