Rare genomic changes as a tool for phylogenetics

Plastome phylogenomics unravels the evolutionary relationships and biogeographic history of Chloranthaceae

BACKGROUND

Chloranthaceae is a pantropical family of flowering plants distributed mainly across Central and South America, East Asia, and the Pacific, encompassing approximately 73 species belonging to four extant genera-Ascarina, Chloranthus, Hedyosmum, and Sarcandra. As one of the most ancient lineages of extant angiosperms, Chloranthaceae holds substantial value in traditional medicine globally and offers major insights into the evolutionary history of flowering plants. However, phylogenetic relationships within this family remain partially resolved, and its origin continues to be debated. We here sequenced, assembled, annotated, and compared the chloroplast genomes (plastomes) of 22 Chloranthaceae species to investigate their plastome evolution and reconstructed the family's phylogeny by using both plastome and nuclear ribosomal DNA (nrDNA). Additionally, we employed the divergence-extinction-cladogenesis model to infer the ancestral area and historical migration patterns of Chloranthaceae.

RESULTS

Results The plastomes of Chloranthaceae exhibit a typical quadripartite structure, with genome sizes ranging from 157,449 bp to 159,218 bp and encoding 113 unique genes, including 17-18 genes duplicated in inverted repeat (IR) regions. Interestingly, the genus Hedyosmum demonstrated notable IR expansions, resulting in trnHGUG duplication in five species. A comparative analysis unveiled substantial variations in simple sequence repeats (SSRs) and tandem repeats across different species, whereas interspersed repeats remained relatively stable. Correlation analyses revealed positive correlations between plastome size and both IR length and SSR count. Furthermore, our phylogenetic reconstruction based on plastome and nrDNA (18 S-ITS1-5.8 S-ITS2-26 S) data provided robust support for the intergeneric relationships, confirming the monophyly of each genus and resolving long-standing uncertainties in the intrageneric relationships of Chloranthus and Hedyosmum. Notably, H. orientale, the sole Hedyosmum species found in Asia, was phylogenetically embedded within the genus, rather than occupying a basal position contrary to previous results. Our biogeographical reconstruction, incorporating both extant and fossil evidence, supports a Laurasian origin of Chloranthaceae and its subsequent dispersal to the tropics.

CONCLUSION

Conclusions This study provides the most comprehensive plastome-based phylogenomic analysis of Chloranthaceae to date, offering novel insights into its plastome evolution, resolving long-standing ambiguities in interspecific and phylogenetic relationships and reconstructing the probable biogeographic origin and historical migrations of this family. The proposed resources establish a robust framework for future evolutionary investigations within the realm of Chloranthaceae.

BAD2matrix: Phylogenomic matrix concatenation, indel coding, and more

Premise

Common steps in phylogenomic matrix production include biological sequence concatenation, morphological data concatenation, insertion/deletion (indel) coding, gene content (presence/absence) coding, removing uninformative characters for parsimony analysis, recording with reduced amino acid alphabets, and occupancy filtering. Existing software does not accomplish these tasks on a phylogenomic scale using a single program.

Methods and Results

BAD2matrix is a Python script that performs the above-mentioned steps in phylogenomic matrix construction for DNA or amino acid sequences as well as morphological data. The script works in UNIX-like environments (e.g., LINUX, MacOS, Windows Subsystem for LINUX).

Conclusions

BAD2matrix helps simplify phylogenomic pipelines and can be downloaded from https://github.com/dpl10/BAD2matrix/tree/master under a GNU General Public License v2.

Annelid Comparative Genomics and the Evolution of Massive Lineage-Specific Genome Rearrangement in Bilaterians

The organization of genomes into chromosomes is critical for processes such as genetic recombination, environmental adaptation, and speciation. All animals with bilateral symmetry inherited a genome structure from their last common ancestor that has been highly conserved in some taxa but seemingly unconstrained in others. However, the evolutionary forces driving these differences and the processes by which they emerge have remained largely uncharacterized. Here, we analyze genome organization across the phylum Annelida using 23 chromosome-level annelid genomes. We find that while many annelid lineages have maintained the conserved bilaterian genome structure, the Clitellata, a group containing leeches and earthworms, possesses completely scrambled genomes. We develop a rearrangement index to quantify the extent of genome structure evolution and show that, compared to the last common ancestor of bilaterians, leeches and earthworms have among the most highly rearranged genomes of any currently sampled species. We further show that bilaterian genomes can be classified into two distinct categories-high and low rearrangement-largely influenced by the presence or absence, respectively, of chromosome fission events. Our findings demonstrate that animal genome structure can be highly variable within a phylum and reveal that genome rearrangement can occur both in a gradual, stepwise fashion, or rapid, all-encompassing changes over short evolutionary timescales.

Insertions and Deletions: Computational Methods, Evolutionary Dynamics, and Biological Applications

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.

Acoelomorph flatworm monophyly is a long-branch attraction artefact obscuring a clade of Acoela and Xenoturbellida

Acoelomorpha is a broadly accepted clade of bilaterian animals made up of the fast-evolving, morphologically simple, mainly marine flatworm lineages Acoela and Nemertodermatida. Phylogenomic studies support Acoelomorpha's close relationship with the slowly evolving and similarly simplistic Xenoturbella, together forming the phylum Xenacoelomorpha. The phylogenetic placement of Xenacoelomorpha amongst bilaterians is controversial, with some studies supporting Xenacoelomorpha as the sister group to all other bilaterians, implying that their simplicity may be representative of early bilaterians. Others propose that this placement is an error resulting from the fast-evolving Acoelomorpha, and instead suggest that they are the degenerate sister group to Ambulacraria. Perhaps as a result of this debate, internal xenacoelomorph relationships have been somewhat overlooked at a phylogenomic scale. Here, I employ a highly targeted approach to detect and overcome possible phylogenomic error in the relationship between Xenoturbella and the fast-evolving acoelomorph flatworms. The results indicate that the subphylum Acoelomorpha is a long-branch attraction artefact obscuring a previously undiscovered clade comprising Xenoturbella and Acoela, which I name Xenacoela. The findings also suggest that Xenacoelomorpha is not the sister group to all other bilaterians. This study provides a template for future efforts aimed at discovering and correcting unrecognized long-branch attraction artefacts throughout the tree of life.

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.

A plurality of morphological characters need not equate with phylogenetic accuracy: A rare genomic change refutes the placement of Solifugae and Pseudoscorpiones in Haplocnemata

Recent advances in higher-level invertebrate phylogeny have leveraged shared features of genomic architecture to resolve contentious nodes across the tree of life. Yet, the interordinal relationships within Chelicerata have remained recalcitrant given competing topologies in recent molecular analyses. As such, relationships between topologically unstable orders remain supported primarily by morphological cladistic analyses. Solifugae, one such unstable chelicerate order, has long been thought to be the sister group of Pseudoscorpiones, forming the clade Haplocnemata, on the basis of eight putative morphological synapomorphies. The discovery, however, of a shared whole genome duplication placing Pseudoscorpiones in Arachnopulmonata provides the opportunity for a simple litmus test evaluating the validity of Haplocnemata. Here, we present the first developmental transcriptome of a solifuge (Titanopuga salinarum) and survey copy numbers of the homeobox genes for evidence of systemic duplication. We find that over 70% of the identified homeobox genes in T. salinarum are retained in a single copy, while representatives of the arachnopulmonates retain orthologs of those genes as two or more copies. Our results refute the placement of Solifugae in Haplocnemata. Subsequent reevaluation of putative interordinal morphological synapomorphies among chelicerates reveals a high incidence of homoplasy, reversals, and inaccurate coding within Haplocnemata and other small clades, as well as Arachnida more broadly, suggesting existing morphological character matrices are insufficient to resolve chelicerate phylogeny.

Whole genome resequencing unveils low-temperature stress tolerance specific genomic variations in jute (Corchorus sp.)

Jute (Corchorus sp.), a commercially important and eco-friendly crop, is widely cultivated in Bangladesh, India, and China. Some varieties of this tropical plant such as the Corchorus olitorius. Variety accession no. 2015 (acc. 2015) has been found to be low-temperature tolerant. The current study was designed to explore the genome-wide variations present in the tolerant plant acc. 2015 in comparison to the sensitive farmer popular variety Corchorus olitorius var. O9897 using the whole genome resequencing technique. Among different variations, intergenic Single Nucleotide Polymorphism (SNPs) and Insertion-Deletion (InDels) were found in the highest percentage whereas approximately 3% SNPs and 2% InDels were found in exonic regions in both plants. Gene enrichment analysis indicated the presence of acc. 2015 specific SNPs in the genes encoding peroxidase, ER lumen protein retaining receptor, and hexosyltransferase involved in stress response (GO:0006950) which were not present in sensitive variety O9897. Besides, distinctive copy number variation regions (CNVRs) comprising 120 gene loci were found in acc. 2015 with a gain of function from multiple copy numbers but absent in O9897. Gene ontology analysis revealed these gene loci to possess different receptors like kinases, helicases, phosphatases, transcription factors especially Myb transcription factors, regulatory proteins containing different binding domains, annexin, laccase, acyl carrier protein, potassium transporter, and vesicular transporter proteins that are responsible for low temperature induced adaptation pathways in plants. This work of identifying genomic variations linked to cold stress tolerance traits will help to develop successful markers that will pave the way to develop genetically modified cold-resistant jute lines for year-round cultivation to meet the demand for a sustainable fiber crop economy.

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.

Phylomitogenomics of two Neotropical species of long-legged crickets Endecous Saussure, 1878 (Orthoptera: Phalangopsidae)

Mitochondrial genomes have provided significant insights into the evolution of several insects. A typical mitogenome contains 37 genes, and variations in gene order can indicate evolutionary relationships between species. In this study, we have assembled the first complete mitogenomes of Endecous chape and E. onthophagus and analyzed the phylogenetic implications for the Gryllidea infraorder. We performed DNA extractions and genome sequencing for both Endecous species. Subsequently, we searched for raw data in the Sequence Read Archive (SRA) in NCBI. Using the SRA data, we assembled the partial mitogenome of Dianemobius nigrofasciatus and annotated the protein-coding genes (PCGs) for nine species. Phylogenomic relationships were reconstructed using Maximum Likelihood (ML) and Bayesian Inference (BI), utilizing the PCGs from 49 Gryllidea species. The mitogenome lengths of E. chape and E. onthophagus are 16,266 bp and 16,023 bp, respectively, while D. nigrofasciatus has a length of 15,359 bp. Our results indicate that species within the infraorder exhibit four types of gene order arrangements that align with their phylogenetic relationships. Both phylogenomic trees displayed strong support, and the ML corroborated with the literature. Gryllidea species have significantly contributed to various fields, and studying their mitogenomes can provide valuable insights into this infraorder evolution.

Phylogenomics resolves the higher-level phylogeny of herbivorous eriophyoid mites (Acariformes: Eriophyoidea)

BACKGROUND

Eriophyoid mites (Eriophyoidea) are among the largest groups in the Acariformes; they are strictly phytophagous. The higher-level phylogeny of eriophyoid mites, however, remains unresolved due to the limited number of available morphological characters-some of them are homoplastic. Nevertheless, the eriophyoid mites sequenced to date showed highly variable mitochondrial (mt) gene orders, which could potentially be useful for resolving the higher-level phylogenetic relationships.

RESULTS

Here, we sequenced and compared the complete mt genomes of 153 eriophyoid mite species, which showed 54 patterns of rearranged mt gene orders relative to that of the hypothetical ancestor of arthropods. The shared derived mt gene clusters support the monophyly of eriophyoid mites (Eriophyoidea) as a whole and the monophylies of six clades within Eriophyoidea. These monophyletic groups and their relationships were largely supported in the phylogenetic trees inferred from mt genome sequences as well. Our molecular dating results showed that Eriophyoidea originated in the Triassic and diversified in the Cretaceous, coinciding with the diversification of angiosperms.

CONCLUSIONS

This study reveals multiple molecular synapomorphies (i.e. shared derived mt gene clusters) at different levels (i.e. family, subfamily or tribe level) from the complete mt genomes of 153 eriophyoid mite species. We demonstrated the use of derived mt gene clusters in unveiling the higher-level phylogeny of eriophyoid mites, and underlines the origin of these mites and their co-diversification with angiosperms.

Phylogenomics studies and molecular markers reliably demarcate genus Pseudomonas sensu stricto and twelve other Pseudomonadaceae species clades representing novel and emended genera

Genus Pseudomonas is a large assemblage of diverse microorganisms, not sharing a common evolutionary history. To clarify their evolutionary relationships and classification, we have conducted comprehensive phylogenomic and comparative analyses on 388 Pseudomonadaceae genomes. In phylogenomic trees, Pseudomonas species formed 12 main clusters, apart from the "Aeruginosa clade" containing its type species, P. aeruginosa. In parallel, our detailed analyses on protein sequences from Pseudomonadaceae genomes have identified 98 novel conserved signature indels (CSIs), which are uniquely shared by the species from different observed clades/groups. Six CSIs, which are exclusively shared by species from the "Aeruginosa clade," provide reliable demarcation of this clade corresponding to the genus Pseudomonas sensu stricto in molecular terms. The remaining 92 identified CSIs are specific for nine other Pseudomonas species clades and the genera Azomonas and Azotobacter which branch in between them. The identified CSIs provide strong independent evidence of the genetic cohesiveness of these species clades and offer reliable means for their demarcation/circumscription. Based on the robust phylogenetic and molecular evidence presented here supporting the distinctness of the observed Pseudomonas species clades, we are proposing the transfer of species from the following clades into the indicated novel genera: Alcaligenes clade - Aquipseudomonas gen. nov.; Fluvialis clade - Caenipseudomonas gen. nov.; Linyingensis clade - Geopseudomonas gen. nov.; Oleovorans clade - Ectopseudomonas gen. nov.; Resinovorans clade - Metapseudomonas gen. nov.; Straminea clade - Phytopseudomonas gen. nov.; and Thermotolerans clade - Zestomonas gen. nov. In addition, descriptions of the genera Azomonas, Azotobacter, Chryseomonas, Serpens, and Stutzerimonas are emended to include information for the CSIs specific for them. The results presented here should aid in the development of a more reliable classification scheme for Pseudomonas species.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Characterization of the Complete Mitochondrial Genome of Agelas nakamurai from the South China Sea

The Agelas genus sponges are widely distributed and provide shelter for organisms that inhabit reefs. However, there is a lack of research on the genetic diversity of the Agelas sponges. Additionally, only one Agelas mitochondrial genome has been documented, leaving the characteristics of the Agelas genus's mitogenome in need of further clarification. To address this research gap, we utilized Illumina HiSeq4000 sequencing and de novo assembly to ascertain the complete mitochondrial genome of Agelas sp. specimens, sourced from the South China Sea. Our analysis of the cox1 barcoding similarity and phylogenetic relationship reveals that taxonomically, the Agelas sp. corresponds to Agelas nakamurai. The mitogenome of Agelas nakamurai is 20,885 bp in length, encoding 14 protein-coding genes, 24 transfer RNA genes, and 2 ribosomal RNA genes. Through a comparison of the mitochondrial genes, we discovered that both Agelas nakamurai and Agelas schmidti have an identical gene arrangement. Furthermore, we observed a deletion in the trnD gene and duplication and remodeling of the trnL gene in the Agelas nakamurai's mitogenome. Our evolutionary analysis also identified lineage-specific positive selection sites in the nad3 and nad5 genes of the Agelas sponges' mitogenome. These findings shed light on the gene rearrangement events and positive selection sites in the mitogenome of Agelas nakamurai, providing valuable molecular insights into the evolutionary processes of this genus.

Generation of somatic de novo structural variation as a hallmark of cellular senescence in human lung fibroblasts

Cellular senescence is characterized by replication arrest in response to stress stimuli. Senescent cells accumulate in aging tissues and can trigger organ-specific and possibly systemic dysfunction. Although senescent cell populations are heterogeneous, a key feature is that they exhibit epigenetic changes. Epigenetic changes such as loss of repressive constitutive heterochromatin could lead to subsequent LINE-1 derepression, a phenomenon often described in the context of senescence or somatic evolution. LINE-1 elements decode the retroposition machinery and reverse transcription generates cDNA from autonomous and non-autonomous TEs that can potentially reintegrate into genomes and cause structural variants. Another feature of cellular senescence is mitochondrial dysfunction caused by mitochondrial damage. In combination with impaired mitophagy, which is characteristic of senescent cells, this could lead to cytosolic mtDNA accumulation and, as a genomic consequence, integrations of mtDNA into nuclear DNA (nDNA), resulting in mitochondrial pseudogenes called numts. Thus, both phenomena could cause structural variants in aging genomes that go beyond epigenetic changes. We therefore compared proliferating and senescent IMR-90 cells in terms of somatic de novo numts and integrations of a non-autonomous composite retrotransposons - the so-called SVA elements-that hijack the retropositional machinery of LINE-1. We applied a subtractive and kinetic enrichment technique using proliferating cell DNA as a driver and senescent genomes as a tester for the detection of nuclear flanks of de novo SVA integrations. Coupled with deep sequencing we obtained a genomic readout for SVA retrotransposition possibly linked to cellular senescence in the IMR-90 model. Furthermore, we compared the genomes of proliferative and senescent IMR-90 cells by deep sequencing or after enrichment of nuclear DNA using AluScan technology. A total of 1,695 de novo SVA integrations were detected in senescent IMR-90 cells, of which 333 were unique. Moreover, we identified a total of 81 de novo numts with perfect identity to both mtDNA and nuclear hg38 flanks. In summary, we present evidence for possible age-dependent structural genomic changes by paralogization that go beyond epigenetic modifications. We hypothesize, that the structural variants we observe potentially impact processes associated with replicative aging of IMR-90 cells.

Phylogenomic and molecular markers based studies on Staphylococcaceae and Gemella species. Proposals for an emended family Staphylococcaceae and three new families (Abyssicoccaceae fam. nov., Salinicoccaceae fam. nov. and Gemellaceae fam. nov.) harboring four new genera, Lacicoccus gen. nov., Macrococcoides gen. nov., Gemelliphila gen. nov., and Phocicoccus gen. nov

The family Staphylococcacae and genus Gemella contain several organisms of clinical or biotechnological importance. We report here comprehensive phylogenomic and comparative analyses on 112 available genomes from species in these taxa to clarify their evolutionary relationships and classification. In a phylogenomic tree based on 678 core proteins, Gemella species were separated from Staphylococcacae by a long branch indicating that they constitute a distinct family (Gemellaceae fam. nov.). In this tree, Staphylococcacae species formed two main clades, one encompassing the genera Aliicoccus, Jeotgalicoccus, Nosocomiicoccus and Salinicoccus (Family "Salinicoccaceae"), while the other clade consisted of the genera Macrococcus, Mammaliicoccus and Staphylococcus (Family Staphylococcaceae emend.). In this tree, species from the genera Gemella, Jeotgalicoccus, Macrococcus and Salinicoccus each formed two distinct clades. Two species clades for these genera are also observed in 16S rRNA gene trees and supported by average amino acid identity analysis. We also report here detailed analyses on protein sequences from Staphylococcaceae and Gemella genomes to identify conserved signature indels (CSIs) which are specific for different genus and family-level clades. These analyses have identified 120 novel CSIs robustly demarcating different proposed families and genera. The identified CSIs provide independent evidence that the genera Gemella, Jeotgalicoccus, Macrococcus and Salinicoccus consist of two distinct clades, which can be reliably distinguished based on multiple exclusively shared CSIs. We are proposing transfers of the species from the novel clades of the above four genera into the genera Gemelliphila gen. nov., Phocicoccus gen. nov., Macrococcoides gen. nov. and Lacicoccus gen. nov., respectively. The identified CSIs also provide strong evidence for division of Staphylococcaceae into an emended family Staphylococcaceae and two new families, Abyssicoccaceae fam. nov. and Salinicoccaceae fam. nov. All of these families can be reliably demarcated based on several exclusively shared CSIs.

Complete mitochondrial genomes of the slugs Deroceraslaeve (Agriolimacidae) and Ambigolimaxvalentianus (Limacidae) provide insights into the phylogeny of Stylommatophora (Mollusca, Gastropoda)

In this study, we sequenced two complete mitogenomes from Deroceraslaeve and Ambigolimaxvalentianus. The mitogenome of Ambigolimaxvalentianus represented the first such data from the family Limacidae. The lengths of the mitogenomes of Deroceraslaeve and Ambigolimaxvalentianus were 14,773 bp and 15,195 bp, respectively. The entire set of 37 mitochondrial genes were identified for both mitogenomes. Compared with the mitogenome of Achatinafulica, the trnP_trnA tRNA cluster was rearranged in both Deroceraslaeve and Ambigolimaxvalentianus. The secondary structures of tRNA and rRNA genes for the two species were predicted. Phylogenetic analyses based on amino acid sequences supported (1) monophyly of Stylommatophora, (2) division of Stylommatophora into the 'achatinoid' clade (i.e., the suborder Achatinina) and the 'non-achatinoid' clade (i.e., the suborder Helicina), (3) placement of the Orthurethra in the 'non-achatinoid' clade, and (4) monophyly of each of the superfamilies Helicoidea, Urocoptoidea, Succineoidea, Arionoidea, Pupilloidea and Limacoidea. The exemplars of Helicidae, Philomycidae and Achatinellidae displayed many more mitochondrial gene rearrangements than other species of Stylommatophora.

Ancient gene linkages support ctenophores as sister to other animals

A central question in evolutionary biology is whether sponges or ctenophores (comb jellies) are the sister group to all other animals. These alternative phylogenetic hypotheses imply different scenarios for the evolution of complex neural systems and other animal-specific traits1-6. Conventional phylogenetic approaches based on morphological characters and increasingly extensive gene sequence collections have not been able to definitively answer this question7-11. Here we develop chromosome-scale gene linkage, also known as synteny, as a phylogenetic character for resolving this question12. We report new chromosome-scale genomes for a ctenophore and two marine sponges, and for three unicellular relatives of animals (a choanoflagellate, a filasterean amoeba and an ichthyosporean) that serve as outgroups for phylogenetic analysis. We find ancient syntenies that are conserved between animals and their close unicellular relatives. Ctenophores and unicellular eukaryotes share ancestral metazoan patterns, whereas sponges, bilaterians, and cnidarians share derived chromosomal rearrangements. Conserved syntenic characters unite sponges with bilaterians, cnidarians, and placozoans in a monophyletic clade to the exclusion of ctenophores, placing ctenophores as the sister group to all other animals. The patterns of synteny shared by sponges, bilaterians, and cnidarians are the result of rare and irreversible chromosome fusion-and-mixing events that provide robust and unambiguous phylogenetic support for the ctenophore-sister hypothesis. These findings provide a new framework for resolving deep, recalcitrant phylogenetic problems and have implications for our understanding of animal evolution.

Nine Mitochondrial Genomes of Phasmatodea with Two Novel Mitochondrial Gene Rearrangements and Phylogeny

The classification of stick and leaf insects (Order Phasmatodea) is flawed at various taxonomic ranks due to a lack of robust phylogenetic relationships and convergent morphological characteristics. In this study, we sequenced nine new mitogenomes that ranged from 15,011 bp to 17,761 bp in length. In the mitogenome of Carausis sp., we found a translocation of trnR and trnA, which can be explained by the tandem duplication/random loss (TDRL) model. In the Stheneboea repudiosa Brunner von Wattenwyl, 1907, a novel mitochondrial structure of 12S rRNA-CR1-trnI-CR2-trnQ-trnM was found for the first time in Phasmatodea. Due to the low homology of CR1 and CR2, we hypothesized that trnI was inverted through recombination and then translocated into the middle of the control region. Control region repeats were frequently detected in the newly sequenced mitogenomes. To explore phylogenetic relationships in Phasmatodea, mtPCGs from 56 Phasmatodean species (composed of 9 stick insects from this study, 31 GenBank data, and 16 data derived from transcriptome splicing) were used for Bayesian inference (BI), and maximum likelihood (ML) analyses. Both analyses supported the monophyly of Lonchodinae and Necrosciinae, but Lonchodidae was polyphyletic. Phasmatidae was monophyletic, and Clitumninae was paraphyletic. Phyllidae was located at the base of Neophasmatodea and formed a sister group with the remaining Neophasmatodea. Bacillidae and Pseudophasmatidae were recovered as a sister group. Heteroptergidae was monophyletic, and the Heteropteryginae sister to the clade (Obriminae + Dataminae) was supported by BI analysis and ML analysis.

Phylogeny and expression of tetraspanin CD9 paralogues in rainbow trout (Oncorhynchus mykiss)

CD9 is a member of the tetraspanin family, which is characterised by a unique domain structure and conserved motifs. In mammals, CD9 is found in tetraspanin-enriched microdomains (TEMs) on the surface of virtually every cell type. CD9 has a wide variety of roles, including functions within the immune system. Here we show the first in-depth analysis of the cd9 gene family in salmonids, showing that this gene has expanded to six paralogues in three groups (cd9a, cd9b, cd9c) through whole genome duplication events. We suggest that through genome duplications, cd9 has undergone subfunctionalisation in the paralogues and that cd9c1 and cd9c2 in particular are involved in antiviral responses in salmonid fish. We show that these paralogues are significantly upregulated in parallel to classic interferon-stimulated genes (ISGs) active in the antiviral response. Expression analysis of cd9 may therefore become an interesting target to assess teleost responses to viruses.

A genomic timescale for placental mammal evolution

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

Evolution of mitogenomic gene order in Orthoptera

Mitochondrial gene order has contributed to the elucidation of evolutionary relationships in several animal groups. It generally has found its application as a phylogenetic marker for deep nodes. Yet, in Orthoptera limited research has been performed on the gene order, although the group represents one of the oldest insect orders. We performed a comprehensive study on mitochondrial genome rearrangements (MTRs) within Orthoptera in the context of mitogenomic sequence-based phylogeny. We used 280 published mitogenome sequences from 256 species, including three outgroup species, to reconstruct a molecular phylogeny. Using a heuristic approach, we assigned MTR scenarios to the edges of the phylogenetic tree and reconstructed ancestral gene orders to identify possible synapomorphies in Orthoptera. We found all types of MTRs in our dataset: inversions, transpositions, inverse transpositions, and tandem-duplication/random loss events (TDRL). Most of the suggested MTRs were in single and unrelated species. Out of five MTRs which were unique in subgroups of Orthoptera, we suggest four of them to be synapomorphies; those were in the infraorder Acrididea, in the tribe Holochlorini, in the subfamily Pseudophyllinae, and in the two families Phalangopsidae and Gryllidae or their common ancestor (leading to the relationship ((Phalangopsidae + Gryllidae) + Trigonidiidae)). However, similar MTRs have been found in distant insect lineages. Our findings suggest convergent evolution of specific mitochondrial gene orders in several species, deviant from the evolution of the mitogenome DNA sequence. As most MTRs were detected at terminal nodes, a phylogenetic inference of deeper nodes based on MTRs is not supported. Hence, the marker does not seem to aid resolving the phylogeny of Orthoptera, but adds further evidence for the complex evolution of the whole group, especially at the genetic and genomic levels. The results indicate a high demand for more research on patterns and underlying mechanisms of MTR events in Orthoptera.

Reconstruction of hundreds of reference ancestral genomes across the eukaryotic kingdom

Ancestral sequence reconstruction is a fundamental aspect of molecular evolution studies and can trace small-scale sequence modifications through the evolution of genomes and species. In contrast, fine-grained reconstructions of ancestral genome organizations are still in their infancy, limiting our ability to draw comprehensive views of genome and karyotype evolution. Here we reconstruct the detailed gene contents and organizations of 624 ancestral vertebrate, plant, fungi, metazoan and protist genomes, 183 of which are near-complete chromosomal gene order reconstructions. Reconstructed ancestral genomes are similar to their descendants in terms of gene content as expected and agree precisely with reference cytogenetic and in silico reconstructions when available. By comparing successive ancestral genomes along the phylogenetic tree, we estimate the intra- and interchromosomal rearrangement history of all major vertebrate clades at high resolution. This freely available resource introduces the possibility to follow evolutionary processes at genomic scales in chronological order, across multiple clades and without relying on a single extant species as reference.

Chloroplast genomes of four Carex species: Long repetitive sequences trigger dramatic changes in chloroplast genome structure

The chloroplast genomes of angiosperms usually have a stable circular quadripartite structure that exhibits high consistency in genome size and gene order. As one of the most diverse genera of angiosperms, Carex is of great value for the study of evolutionary relationships and speciation within its genus, but the study of the structure of its chloroplast genome is limited due to its highly expanded and restructured genome with a large number of repeats. In this study, we provided a more detailed account of the chloroplast genomes of Carex using a hybrid assembly of second- and third-generation sequencing and examined structural variation within this genus. The study revealed that chloroplast genomes of four Carex species are significantly longer than that of most angiosperms and are characterized by high sequence rearrangement rates, low GC content and gene density, and increased repetitive sequences. The location of chloroplast genome structural variation in the species of Carex studied is closely related to the positions of long repeat sequences; this genus provides a typical example of chloroplast structural variation and expansion caused by long repeats. Phylogenetic relationships constructed based on the chloroplast protein-coding genes support the latest taxonomic system of Carex, while revealing that structural variation in the chloroplast genome of Carex may have some phylogenetic significance. Moreover, this study demonstrated a hybrid assembly approach based on long and short reads to analyze complex chloroplast genome assembly and also provided an important reference for the analysis of structural rearrangements of chloroplast genomes in other taxa.

Filtering artifactual signal increases support for Xenacoelomorpha and Ambulacraria sister relationship in the animal tree of life

Conflicting studies place a group of bilaterian invertebrates containing xenoturbellids and acoelomorphs, the Xenacoelomorpha, as either the primary emerging bilaterian phylum^1,2,3,4,5,6 or within Deuterostomia, sister to Ambulacraria.^7,8,9,10,11 Although their placement as sister to the rest of Bilateria supports relatively simple morphology in the ancestral bilaterian, their alternative placement within Deuterostomia suggests a morphologically complex ancestral bilaterian along with extensive loss of major phenotypic traits in the Xenacoelomorpha. Recent studies have questioned whether Deuterostomia should be considered monophyletic at all.^10,12,13 Hidden paralogy and poor phylogenetic signal present a major challenge for reconstructing species phylogenies.^{14,15,16,17,18} Here, we assess whether these issues have contributed to the conflict over the placement of Xenacoelomorpha. We reanalyzed published datasets, enriching for orthogroups whose gene trees support well-resolved clans elsewhere in the animal tree.¹⁶ We find that most genes in previously published datasets violate incontestable clans, suggesting that hidden paralogy and low phylogenetic signal affect the ability to reconstruct branching patterns at deep nodes in the animal tree. We demonstrate that removing orthogroups that cannot recapitulate incontestable relationships alters the final topology that is inferred, while simultaneously improving the fit of the model to the data. We discover increased, but ultimately not conclusive, support for the existence of Xenambulacraria in our set of filtered orthogroups. At a time when we are progressing toward sequencing all life on the planet, we argue that long-standing contentious issues in the tree of life will be resolved using smaller amounts of better quality data that can be modeled adequately.¹⁹.

Complete chloroplast genomes of Sorbus sensu stricto (Rosaceae): comparative analyses and phylogenetic relationships

BACKGROUND

Sorbus sensu stricto (Sorbus s.s.) is a genus with important economical values because of its beautiful leaves, and flowers and especially the colorful fruits. It belongs to the tribe Maleae of the family Rosaceae, and comprises about 90 species mainly distributed in China. There is on-going dispute about its infrageneric classification and species delimitation as the species are morphologically similar. With the aim of shedding light on the circumscription of taxa within the genus, phylogenetic analyses were performed using 29 Sorbus s.s. chloroplast (cp) genomes (16 newly sequenced) representing two subgenera and eight sections.

RESULTS

The 16 cp genomes newly sequenced range between 159,646 bp and 160,178 bp in length. All the samples examined and 22 taxa re-annotated in Sorbus sensu lato (Sorbus s.l.) contain 113 unique genes with 19 of these duplicated in the inverted repeat (IR). Six hypervariable regions including trnR-atpA, petN-psbM, rpl32-trnL, trnH-psbA, trnT-trnL and ndhC-trnV were screened and 44-53 SSRs and 14-31 dispersed repeats were identified as potential molecular markers. Phylogenetic analyses under ML/BI indicated that Sorbus s.l. is polyphyletic, but Sorbus s.s. and the other five segregate genera, Aria, Chamaemespilus, Cormus, Micromeles and Torminalis are monophyletic. Two major clades and four sub-clades resolved with full-support within Sorbus s.s. are not consistent with the existing infrageneric classification. Two subgenera, subg. Sorbus and subg. Albocarmesinae are supported as monophyletic when S. tianschanica is transferred to subg. Albocarmesinae from subg. Sorbus and S. hupehensis var. paucijuga transferred to subg. Sorbus from subg. Albocarmesinae, respectively. The current classification at sectional level is not supported by analysis of cp genome phylogeny.

CONCLUSION

Phylogenomic analyses of the cp genomes are useful for inferring phylogenetic relationships in Sorbus s.s. Though genome structure is highly conserved in the genus, hypervariable regions and repeat sequences used are the most promising molecule makers for population genetics, species delimitation and phylogenetic studies.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

BACKGROUND

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

RESULTS

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

CONCLUSIONS

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.

Phylogenomic Analyses of the Tenthredinoidea Support the Familial Rank of Athaliidae (Insecta, Tenthredinoidea)

The systematic status of the genus Athalia and related genera is a perennial controversy in sawfly taxonomy. Several authors have hypothesized that the placement of Athalia within the Tenthredinidae is artificial, but no studies have focused on this topic. If the hypothesis that Athalia does not belong to Tenthredinidae can be supported, the taxonomic framework of Tenthredinoidea needs revision. We present a comprehensive phylogenomic study of Tenthredinoidae, focusing on the positions of Athalia and related genera by sampling 80 representatives mainly of the Tenthredinoidea, including Heptamelinae and Blasticotomidae. Our phylogenetic reconstructions based on nuclear genes and mitochondrial (mt) sequences support Athalia and related genera as a distinct clade sister to Tenthredinidae + (Cimbicidae + Diprionidae). A comparison of symphytan mitochondrial genomes reveals an innovative gene rearrangement pattern in Athaliidae, in which Dentathalia demonstrates a more ancestral pattern than Athalia and Hypsathalia. The lineage specificity of mt rRNA secondary structures also provides sufficient support to consider Athaliidae as a separate family. In summary, the phylogeny and genomic structural changes unanimously support the taxonomic treatment of Athaliidae as a family and the re-establishment of Dentathalia as a valid genus.

Revealing the Complete Chloroplast Genome of an Andean Horticultural Crop, Sweet Cucumber (Solanum muricatum), and Its Comparison with Other Solanaceae Species

Sweet cucumber (Solanum muricatum) sect. Basarthrum is a neglected horticultural crop native to the Andean region. It is naturally distributed very close to other two Solanum crops of high importance, potatoes, and tomatoes. To date, molecular tools for this crop remain undetermined. In this study, the complete sweet cucumber chloroplast (cp) genome was obtained and compared with seven Solanaceae species. The cp genome of S. muricatum was 155,681 bp in length and included a large single copy (LSC) region of 86,182 bp and a small single-copy (SSC) region of 18,360 bp, separated by a pair of inverted repeats (IR) regions of 25,568 bp. The cp genome possessed 87 protein-coding genes (CDS), 37 transfer RNA (tRNA) genes, eight ribosomal RNA (rRNA) genes, and one pseudogene. Furthermore, 48 perfect microsatellites were identified. These repeats were mainly located in the noncoding regions. Whole cp genome comparative analysis revealed that the SSC and LSC regions showed more divergence than IR regions. Similar to previous studies, our phylogenetic analysis showed that S. muricatum is a sister species to members of sections Petota + Lycopersicum + Etuberosum. We expect that this first sweet cucumber chloroplast genome will provide potential molecular markers and genomic resources to shed light on the genetic diversity and population studies of S. muricatum, which will allow us to identify varieties and ecotypes. Finally, the features and the structural differentiation will provide us with information about the genes of interest, generating tools for the most precise selection of the best individuals of sweet cucumber, in less time and with fewer resources.

Comparative Mitogenomic Analyses of Hydropsychidae Revealing the Novel Rearrangement of Protein-Coding Gene and tRNA (Trichoptera: Annulipalpia)

Gene rearrangement of the mitochondrial genome of insects, especially the rearrangement of protein-coding genes, has long been a hot topic for entomologists. Although mitochondrial gene rearrangement is common within Annulipalpia, protein-coding gene rearrangement is relatively rare. As the largest family in Annulipalpia, the available mitogenomes from Hydropsychidae Curtis, 1835 are scarce, and thus restrict our interpretation of the mitogenome characteristic. In this study, we obtained 19 novel mitogenomes of Hydropsychidae, of which the mitogenomes of the genus Arctopsyche are published for the first time. Coupled with published hydropsychid mitogenome, we analyzed the nucleotide composition evolutionary rates and gene rearrangements of the mitogenomes among subfamilies. As a result, we found two novel gene rearrangement patterns within Hydropsychidae, including rearrangement of protein-coding genes. Meanwhile, our results consider that the protein-coding gene arrangement of Potamyia can be interpreted by the tandem duplication/random loss (TDRL) model. In addition, the phylogenetic relationships within Hydropsychidae constructed by two strategies (Bayesian inference and maximum likelihood) strongly support the monophyly of Arctopscychinae, Diplectroninae, Hydropsychinae, and Macronematinae. Our study provides new insights into the mechanisms and patterns of mitogenome rearrangements in Hydropsychidae.

Caryophylliids (Anthozoa, Scleractinia) and mitochondrial gene order: insights from mitochondrial and nuclear phylogenomics

Molecularly, the family Caryophylliidae is polyphyletic and different sets of genetic data converge towards a consensus that a taxonomic review of this family is necessary. Overall, the order of genes in the mitochondrial genome (mitogenome) together with DNA sequences have been used to successfully untangle evolutionary relationships in several groups of organisms. Published mitogenomes of two caryophylliid genera (Desmophyllum and Solenosmilia) present a transposition of the gene block containing cob, nad2, and nad6, which is located between nad5 5' exon and trnW, while that of Polycyathus chaishanensis presents the same gene order as the majority of scleractinian corals. In molecular-based evolutionary reconstructions, caryophylliids that have the mitochondrial gene rearrangement were recovered as a monophyletic lineage ("true" caryophylliids), while members of the genus Polycyathus were placed in a different position. In this study, additional mitogenomes of this family were assembled and included in evolutionary reconstructions of Scleractinia in order to improve our understanding on whether the mitogenome gene rearrangement is limited to and, therefore, could be a synapomorphy of the actual members of Caryophylliidae. Specimens of Caryophyllia scobinosa, Premocyathus sp., Heterocyathus sulcatus, and Trochocyathus caryophylloides, as well as Desmophyllum pertusum and Solenosmilia variabilis from the Southwest Atlantic were sequenced using Illumina platforms. Then, mitochondrial genomes were assembled and annotated, and nuclear datasets were recovered in-silico from assembled contigs using a previously published set of baits. Evolutionary reconstructions were performed using mitochondrial and nuclear datasets and based on Maximum Likelihood and Bayesian Inference. Obtained mitogenomes are circular and range between 15,816 and 18,225 bp in size and from 30.76% to 36.63% in GC content. The gene rearrangement is only seen in C. scobinosa, D. pertusum, Premocyathus sp., and S. variabilis, which were recovered as a monophyletic clade in both mitochondrial and nuclear phylogenies. On the other hand, the "caryophylliids" with the canonical mitogenome gene order were not recovered within this clade. Differences in features of the skeleton of "true" caryophylliids in comparison to traditional members of the family were observed and offer further support that the gene rearrangement might be seen as a synapomorphy of family Caryophylliidae.

Phylogenomic Analyses and Molecular Signatures Elucidating the Evolutionary Relationships amongst the Chlorobia and Ignavibacteria Species: Robust Demarcation of Two Family-Level Clades within the Order Chlorobiales and Proposal for the Family Chloroherpetonaceae fam. nov

Evolutionary relationships amongst Chlorobia and Ignavibacteria species/strains were examined using phylogenomic and comparative analyses of genome sequences. In a phylogenomic tree based on 282 conserved proteins, the named Chlorobia species formed a monophyletic clade containing two distinct subclades. One clade, encompassing the genera Chlorobaculum, Chlorobium, Pelodictyon, and Prosthecochloris, corresponds to the family Chlorobiaceae, whereas another clade, harboring Chloroherpeton thalassium, Candidatus Thermochlorobacter aerophilum, Candidatus Thermochlorobacteriaceae bacterium GBChlB, and Chlorobium sp. 445, is now proposed as a new family (Chloroherpetonaceae fam. nov). In parallel, our comparative genomic analyses have identified 47 conserved signature indels (CSIs) in diverse proteins that are exclusively present in members of the class Chlorobia or its two families, providing reliable means for identification. Two known Ignavibacteria species in our phylogenomic tree are found to group within a larger clade containing several Candidatus species and uncultured Chlorobi strains. A CSI in the SecY protein is uniquely shared by the species/strains from this “larger Ignavibacteria clade”. Two additional CSIs, which are commonly shared by Chlorobia species and the “larger Ignavibacteria clade”, support a specific relationship between these two groups. The newly identified molecular markers provide novel tools for genetic and biochemical studies and identification of these organisms.

Phylogenomic Coalescent Analyses of Avian Retroelements Infer Zero-Length Branches at the Base of Neoaves, Emergent Support for Controversial Clades, and Ancient Introgressive Hybridization in Afroaves

Retroelement insertions (RIs) are low-homoplasy characters that are ideal data for addressing deep evolutionary radiations, where gene tree reconstruction errors can severely hinder phylogenetic inference with DNA and protein sequence data. Phylogenomic studies of Neoaves, a large clade of birds (>9000 species) that first diversified near the Cretaceous−Paleogene boundary, have yielded an array of robustly supported, contradictory relationships among deep lineages. Here, we reanalyzed a large RI matrix for birds using recently proposed quartet-based coalescent methods that enable inference of large species trees including branch lengths in coalescent units, clade-support, statistical tests for gene flow, and combined analysis with DNA-sequence-based gene trees. Genome-scale coalescent analyses revealed extremely short branches at the base of Neoaves, meager branch support, and limited congruence with previous work at the most challenging nodes. Despite widespread topological conflicts with DNA-sequence-based trees, combined analyses of RIs with thousands of gene trees show emergent support for multiple higher-level clades (Columbea, Passerea, Columbimorphae, Otidimorphae, Phaethoquornithes). RIs express asymmetrical support for deep relationships within the subclade Afroaves that hints at ancient gene flow involving the owl lineage (Strigiformes). Because DNA-sequence data are challenged by gene tree-reconstruction error, analysis of RIs represents one approach for improving gene tree-based methods when divergences are deep, internodes are short, terminal branches are long, and introgressive hybridization further confounds species−tree inference.

Categorical edge-based analyses of phylogenomic data reveal conflicting signals for difficult relationships in the avian tree

Phylogenetic analyses fail to yield a satisfactory resolution of some relationships in the tree of life even with genome-scale datasets, so the failure is unlikely to reflect limitations in the amount of data. Gene tree conflicts are particularly notable in studies focused on these contentious nodes, and taxon sampling, different analytical methods, and/or data type effects can further confound analyses. Although many efforts have been made to incorporate biological conflicts, few studies have curated individual genes for their efficiency in phylogenomic studies. Here, we conduct an edge-based analysis of Neoavian evolution, examining the phylogenetic efficacy of two recent phylogenomic bird datasets and three datatypes (ultraconserved elements [UCEs], introns, and coding regions). We assess the potential causes for biases in signal-resolution for three difficult nodes: the earliest divergence of Neoaves, the position of the enigmatic Hoatzin (Opisthocomus hoazin), and the position of owls (Strigiformes). We observed extensive conflict among genes for all data types and datasets even after meticulous curation. Edge-based analyses (EBA) increased congruence and provided information about the impact of data type, GC content variation (GC_CV), and outlier genes on each of nodes we examined. First, outlier gene signals appeared to drive different patterns of support for the relationships among the earliest diverging Neoaves. Second, the placement of Hoatzin was highly variable, although our EBA did reveal a previously unappreciated data type effect with an impact on its position. It also revealed that the resolution with the most support here was Hoatzin + shorebirds. Finally, GC_CV, rather than data type (i.e., coding vs non-coding) per se, was correlated with a signal that supports monophyly of owls + Accipitriformes (hawks, eagles, and vultures). Eliminating high GC_CV loci increased the signal for owls + mousebirds. Categorical EBA was able to reveal the nature of each edge and provide a way to highlight especially problematic branches that warrant a further examination. The current study increases our understanding about the contentious parts of the avian tree, which show even greater conflicts than appreciated previously.

Genomic evidence for homoploid hybrid speciation between ancestors of two different genera

Homoploid hybrid speciation (HHS) has been increasingly recognized as occurring widely during species diversification of both plants and animals. However, previous studies on HHS have mostly focused on closely-related species while it has been rarely reported or tested between ancestors of different genera. Here, we explore the likely HHS origin of Carpinus sect. Distegocarpus between sect. Carpinus and Ostrya in the family Betulaceae. We generate a chromosome-level reference genome for C. viminea of sect. Carpinus and re-sequence genomes of 44 individuals from the genera Carpinus and Ostrya. Our integrated analyses of all genomic data suggest that sect. Distegocarpus, which has three species, likely originates through HHS during the early divergence between Carpinus and Ostrya. Our study highlights the likelihood of an HHS event between ancestors of the extant genera during their initial divergences, which may have led to reticulate phylogenies at higher taxonomic levels.

Champagne: Automated Whole-Genome Phylogenomic Character Matrix Method Using Large Genomic Indels for Homoplasy-Free Inference

We present Champagne, a whole-genome method for generating character matrices for phylogenomic analysis using large genomic indel events. By rigorously picking orthologous genes and locating large insertion and deletion events, Champagne delivers a character matrix that considerably reduces homoplasy compared with morphological and nucleotide-based matrices, on both established phylogenies and difficult-to-resolve nodes in the mammalian tree. Champagne provides ample evidence in the form of genomic structural variation to support incomplete lineage sorting and possible introgression in Paenungulata and human-chimp-gorilla which were previously inferred primarily through matrices composed of aligned single-nucleotide characters. Champagne also offers further evidence for Myomorpha as sister to Sciuridae and Hystricomorpha in the rodent tree. Champagne harbors distinct theoretical advantages as an automated method that produces nearly homoplasy-free character matrices on the whole-genome scale.

Phylogenomic and comparative genomic analyses of Leuconostocaceae species: identification of molecular signatures specific for the genera Leuconostoc, Fructobacillus and Oenococcus and proposal for a novel genus Periweissella gen. nov

The genera Convivina, Fructobacillus, Leuconostoc, Oenococcus and Weissella, which formed the family Leuconostocaceae, have recently been merged within the family Lactobacillaceae. Using genome sequences for 47 of the 52 named species from these genera, we report here comprehensive phylogenomic and comparative analyses on protein sequences from these species using multiple approaches. In a phylogenomic tree based on concatenated sequences of 498 core proteins from these five genera, and in a 16S rRNA gene tree, members of the genera Fructobacillus, Leuconostoc and Oenococcus formed distinct strongly supported clades. In contrast, Weissella species grouped into two distinct unrelated clades designated as the ‘Weissella main clade’ and ‘Weissella clade 2’. The presence of these clades is also seen in a matrix of pairwise average amino acid identity based on core protein sequences. In parallel, comparative genomic studies on protein sequences from Leuconostocaceae genomes have identified 46 conserved signature indels (CSIs) in diverse proteins that are unique characteristics of the different observed species clades. Of these identified CSIs, five, five and 13 CSIs are uniquely present in members of the genera Fructobacillus, Leuconostoc and Oenococcus, respectively. We also report here six and five CSIs that are exclusively present in the species from the Weissella main clade and Weissella clade 2, respectively, providing independent evidence supporting their distinctness from each other. The remaining 12 identified CSIs are commonly shared by some or all of the species from the genera Convivina, Fructobacillus and Leuconostoc, clarifying their interrelationships. The identified CSIs provide novel and reliable means for the identification/circumscription of members of the genera Fructobacillus, Leuconostoc and Oenococcus as well as the two Weissella species clades in molecular terms. Based on the strong phylogenetic and molecular evidence presented here, we propose that the genus Weissella be limited to only the species from the Weissella main clade, whereas the species forming Weissella clade 2 should be transferred to a new genus Periweissella gen. nov.

Conserved Signatures in Protein Sequences Reliably Demarcate Different Clades of Rodents/Glires Species and Consolidate Their Evolutionary Relationships

The grandorder Glires, consisting of the orders Rodentia and Lagomorpha, encompasses a significant portion of the extant mammalian species including Rat, Mouse, Squirrel, Guinea pig and Beaver. Glires species play an important role in the ecosystem and provide valuable animal models for genetic studies and animal testing. Thus, it is important to reliably determine their evolutionary relationships and identify molecular characteristics that are specific for different species groups within the Glires. In this work, we have constructed a phylogenetic tree for >30 genome sequenced Glires species based on concatenated sequences of 25 conserved proteins. In this tree, members of different orders, suborders, and families within Glires formed strongly supported clades, and their interrelationships were also generally reliably resolved. In parallel, we conducted comparative analyses on more than 1500 protein sequences from Glires species to identify highly conserved molecular markers. These markers were comprised of conserved signature indels (CSIs) in proteins, which are specific for different Rodentia/Glires clades. Of the 41 novel CSIs identified in this work, some are specific for the entire Glires, Rodentia, or Lagomorpha clades, whereas many others reliably demarcate different family/suborder level clades of Rodentia (viz. Myomorpha, Castorimorpha, Sciuromorpha, Hystricomorpha, and Muroidea). Additionally, some of the CSIs also provide information regarding the interrelationships among Rodentia subgroups. Our analysis has also identified one CSI that is commonly shared by the Glires and Scandentia species (tree shrew), however, its evolutionary significance is unclear. Several of the identifed rodents-specific CSIs are present in conserved disease-related proteins. Thus, they provide novel molecular markers for genetic and biochemical studies on the functions of these proteins.

A chromosome-level genome assembly of an alpine plant Crucihimalaya lasiocarpa provides insights into high-altitude adaptation

It remains largely unknown how plants adapt to high-altitude habitats. Crucihimalaya (Brassicaceae) is an alpine genus occurring in the Qinghai-Tibet Plateau characterized by cold temperatures and strong ultraviolet radiation. Here, we generated a chromosome-level genome for C. lasiocarpa with a total size of 255.8 Mb and a scaffold N50 size of 31.9 Mb. We first examined the karyotype origin of this species and found that the karyotype of five chromosomes resembled the ancestral karyotype of the Brassicaceae family, while the other three showed strong chromosomal structural variations. In combination with the rough genome sequence of another congener (C. himalaica), we found that the significantly expanded gene families and positively selected genes involved in alpine adaptation have occurred since the origin of this genus. Our new findings provide valuable information for the chromosomal karyotype evolution of Brassicaceae and investigations of high-altitude environment adaptation of the genus.

Chromosome Painting in Lonchorhina aurita Sheds Light onto the Controversial Phylogenetic Position of Sword-Nosed Bats (Chiroptera, Phyllostomidae)

The subfamily Lonchorhininae encompasses 6 species of sword-nosed bats (Lonchorhina) and is one of the most problematic lineages in the Neotropical leaf-nosed bats (Phyllostomidae) phylogeny. There are at least 5 different hypotheses to explain when the subfamily diverged from the remaining phyllostomids, but none with robust statistical support. Here, we generated a chromosome painting homology map of Lonchorhina aurita karyotype (2n = 32 and FN = 60) using whole-chromosome probes of Macrotus californicus (MCA; 2n = 40 and FN = 60). We placed the karyotype changes of L. aurita in a phylogenetic context to discuss the most likely branching position of Lonchorhininae based on karyotypic evolution. We show that L. aurita has a derived karyotype with 24 segments homologous to the 20 MCA chromosomes used as probes. Comparative analyses between 7 published painted bats species across 4 phyllostomid subfamilies (Macrotinae, Phyllostominae, Glossophaginae, and Lonchophyllinae) revealed that one inversion (MCA 4inv) and one fusion (MCA 17 + 18) are shared derived features between the karyotypes of L. aurita and species of Phyllostominae not yet observed in other bats. Our data show that chromosomal homology maps may contribute with new insights into a long-standing phylogenetic debate that has endured for decades.

First Report on Mitochondrial Gene Rearrangement in Non-Biting Midges, Revealing a Synapomorphy in Stenochironomus Kieffer (Diptera: Chironomidae)

Simple Summary

Gene rearrangement is an additional type of data to support relationships of taxa, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. The concept to use mitochondrial gene rearrangements as phylogenetic markers has been proposed since the mid-1980s, the synapomorphic gene rearrangements have been identified from many lineages. However, mitochondrial gene rearrangement has never been observed in the non-biting midges (Diptera: Chironomidae). Here, seven new mitogenomes of the genus Stenochironomus were sequenced and analyzed. Coupled with published data, phylogenetic analyses were performed within Chironominae. The present study showed that mitogenomes of Stenochironomus are showing a higher A+T bias than other chironomid species. A synapomorphic gene rearrangement that the gene order rearranges from trnI-trnQ-trnM to trnI-trnM-trnQ was identified within Stenochironomus, which is the first instance of mitochondrial gene rearrangement discovered in the Chironomidae. The monophyly of the genus Stenochironomus was strongly supported by mitogenomes. Our study provides new insights into the mitochondrial gene order of Chironomidae, and provides a valuable resource for understanding synapomorphic gene rearrangements.

Abstract

(1) Background: Gene rearrangement of mitochondrial genome, especially those with phylogenetic signals, has long fascinated evolutionary biologists. The synapomorphic gene rearrangements have been identified across multiple orders and at many different taxonomic levels, supporting the monophyletic or systematic relationships of related lineages. However, mitochondrial gene rearrangement has never been observed in the non-biting midges (Diptera: Chironomidae); (2) methods: in this study, the complete mitogenomes of seven Stenochironomus species were sequenced and analyzed for the first time; (3) results: each mitogenome of Stenochironomus contains 37 typical genes and a control region. The whole mitogenomes of Stenochironomus species exhibit a higher A+T bias than other published chironomid species. The gene order rearranges from trnI-trnQ-trnM to trnI-trnM-trnQ in all the seven mitogenomes of Stenochironomus, which might be act as a synapomorphy of the genus, supporting the monophyletic of Stenochironomus species. In addition, another derived gene cluster: trnA-trnG-ND3-trnR exists in Stenochironomus tobaduodecimus. The derived gene orders described above are the first case of mitochondrial gene rearrangement in Chironomidae. Coupled with published data, phylogenetic relationships were reconstructed within Chironominae, and strongly supported the monophyly of Stenochironomus; (4) conclusions: our study provides new insights into the mitochondrial gene order of Chironomidae, and provides a valuable resource for understanding the synapomorphic gene rearrangements.

OUP accepted manuscript

Homology of highly divergent genes often cannot be determined from sequence similarity alone. For example, we recently identified in the aphid Hormaphis cornu a family of rapidly evolving bicycle genes, which encode novel proteins implicated as plant gall effectors, and sequence similarity search methods yielded few putative bicycle homologs in other species. Coding sequence-independent features of genes, such as intron-exon boundaries, often evolve more slowly than coding sequences, however, and can provide complementary evidence for homology. We found that a linear logistic regression classifier using only structural features of bicycle genes identified many putative bicycle homologs in other species. Independent evidence from sequence features and intron locations supported homology assignments. To test the potential roles of bicycle genes in other aphids, we sequenced the genome of a second gall-forming aphid, Tetraneura nigriabdominalis and found that many bicycle genes are strongly expressed in the salivary glands of the gall forming foundress. In addition, bicycle genes are strongly overexpressed in the salivary glands of a non-gall forming aphid, Acyrthosiphon pisum, and in the non-gall forming generations of H. cornu. These observations suggest that Bicycle proteins may be used by multiple aphid species to manipulate plants in diverse ways. Incorporation of gene structural features into sequence search algorithms may aid identification of deeply divergent homologs, especially of rapidly evolving genes involved in host-parasite interactions.

Conserved molecular signatures in the spike protein provide evidence indicating the origin of SARS-CoV-2 and a Pangolin-CoV (MP789) by recombination(s) between specific lineages of Sarbecoviruses

Both SARS-CoV-2 and SARS coronaviruses (CoVs) are members of the subgenus Sarbecovirus. To understand the origin of SARS-CoV-2, sequences for the spike and nucleocapsid proteins from sarbecoviruses were analyzed to identify molecular markers consisting of conserved inserts or deletions (termed CSIs) that are specific for either a particular clade of Sarbecovirus or are commonly shared by two or more clades of these viruses. Three novel CSIs in the N-terminal domain (NTD) of the spike protein S1-subunit (S1-NTD) are uniquely shared by SARS-CoV-2, Bat-CoV-RaTG13 and most pangolin CoVs (SARS-CoV-2r clade). Three other sarbecoviruses viz. bat-CoVZXC21, -CoVZC45 and -PrC31 (forming CoVZC/PrC31 clade), and a pangolin-CoV_MP789 also contain related CSIs in the same positions. In contrast to the S1-NTD, both SARS and SARS-CoV-2r viruses contain two large CSIs in the S1-C-terminal domain (S1-CTD) that are absent in the CoVZC/PrC31 clade. One of these CSIs, consisting of a 12 aa insert, is also present in the RShSTT clade (Cambodia-CoV strains). Sequence similarity studies show that the S1-NTD of SARS-CoV-2r viruses is most similar to the CoVZC/PrC31 clade, whereas their S1-CTD exhibits highest similarity to the RShSTT- (and the SARS-related) CoVs. Results from the shared presence of CSIs and sequence similarity studies on different CoV lineages support the inference that the SARS-CoV-2r cluster of viruses has originated by a genetic recombination between the S1-NTD of the CoVZC/PrC31 clade of CoVs and the S1-CTD of RShSTT/SARS viruses, respectively. We also present compelling evidence, based on the shared presence of CSIs and sequence similarity studies, that the pangolin-CoV_MP789, whose receptor-binding domain is most similar to the SARS-CoV-2 virus, has resulted from another independent recombination event involving the S1-NTD of the CoVZC/PrC31 CoVs and the S1-CTD of an unidentified SARS-CoV-2r related virus. The SARS-CoV-2 virus involved in this latter recombination event is postulated to be most similar to the SARS-CoV-2. Several other CSIs reported here are specific for other clusters of sarbecoviruses including a clade consisting of bat-SARS-CoVs (BM48-31/BGR/2008 and SARS_BtKY72). Structural mapping studies show that the identified CSIs form distinct loops/patches on the surface of the spike protein. It is hypothesized that these novel loops/patches on the spike protein, through their interactions with other host components, should play important roles in the biology/pathology of SARS-CoV-2 virus. Lastly, the CSIs specific for different clades of sarbecoviruses including SARS-CoV-2r clade provide novel means for the identification of these viruses and other potential applications.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.

Phylogenomics and molecular signatures support division of the order Neisseriales into emended families Neisseriaceae and Chromobacteriaceae and three new families Aquaspirillaceae fam. nov., Chitinibacteraceae fam. nov., and Leeiaceae fam. nov

The order Neisseriales contains 37 genera harboring 122 species with validly published names, which are placed into two families, Neisseriaceae and Chromobacteriaceae. Genome sequences are now available for 35 of the 37 Neisseriales genera for reliably determining their evolutionary relationships and taxonomy. We report here comprehensive phylogenomic and comparative analyses on protein sequences from 110 Neisseriales genomes plus 3 Chitinimonas genomes using multiple approaches. In a phylogenomic tree based on 596 core proteins, Neisseriales species formed 5 strongly supported clades. In addition to the clades for Neisseriaceae and Chromobacteriaceae families, three novel species clades designated as the "Chitinibacteraceae", "Aquaspirillaceae", and "Leeiaceae" were observed. The genus Chitinimonas grouped reliably with members of the "Chitinibacteraceae" clade. The major clades within the order Neisseriales can also be distinguished based on average amino acid identity analysis. In parallel, our comparative genomic studies have identified 30 conserved signature indels (CSIs) that are specific for members of the order Neisseriales or its five main clades. One of these CSIs is uniquely shared by all Neisseriales, whereas 8, 4, 9, 3 and 5 CSIs are distinctive characteristics of the Neisseriaceae, Chromobacteriaceae, "Chitinibacteraceae", "Aquaspirillaceae" and "Leeiaceae" clades, respectively. Based on the strong phylogenetic and molecular evidence presented here, we are proposing that the three newly identified clades should be recognized as novel families (Chitinibacteraceae fam. nov., Aquaspirillaceae fam. nov. and Leeiaceae fam. nov.) within the order Neisseriales. In addition, we are also emending descriptions of the families Neisseriaceae and Chromobacteriaceae regarding their constituent genera and other distinguishing characteristics.

Comparative Analysis of Mitochondrial Genomes in Gryllidea (Insecta: Orthoptera): Implications for Adaptive Evolution in Ant-Loving Crickets

Species of infraorder Gryllidea, or crickets, are useful invertebrate models for studying developmental biology and neuroscience. They have also attracted attention as alternative protein sources for human food and animal feed. Mitochondrial genomic information on related invertebrates, such as katydids, and locusts, has recently become available in attempt to clarify the controversial classification schemes, although robust phylogenetic relationships with emphasis on crickets remain elusive. Here, we report newly sequenced complete mitochondrial genomes of crickets to study their phylogeny, genomic rearrangements, and adaptive evolution. First, we conducted de novo assembly of mitochondrial genomes from eight cricket species and annotated protein-coding genes and transfer and ribosomal RNAs using automatic annotations and manual curation. Next, by combining newly described protein-coding genes with public data of the complete Gryllidea genomes and gene annotations, we performed phylogenetic analysis and found gene order rearrangements in several branches. We further analyzed genetic signatures of selection in ant-loving crickets (Myrmecophilidae), which are small wingless crickets that inhabit ant nests. Three distinct approaches revealed two positively selected sites in the cox1 gene in these crickets. Protein 3D structural analyses suggested that these selected sites could influence the interaction of respiratory complex proteins, conferring benefits to ant-loving crickets with a unique ecological niche and morphology. These findings enhance our understanding of the genetic basis of cricket evolution without relying on estimates based on a limited number of molecular markers.

Complete Mitochondrial Genomes of Metcalfa pruinosa and Salurnis marginella (Hemiptera: Flatidae): Genomic Comparison and Phylogenetic Inference in Fulgoroidea

The complete mitochondrial genomes (mitogenomes) of two DNA barcode-defined haplotypes of Metcalfa pruinosa and one of Salurnis marginella (Hemiptera: Flatidae) were sequenced and compared to those of other Fulgoroidea species. Furthermore, the mitogenome sequences were used to reconstruct phylogenetic relationships among fulgoroid families. The three mitogenomes, including that of the available species of Flatidae, commonly possessed distinctive structures in the 1702-1836 bp A+T-rich region, such as two repeat regions at each end and a large centered nonrepeat region. All members of the superfamily Fulgoroidea, including the Flatidae, consistently possessed a motiflike sequence (TAGTA) at the ND1 and trnS₂ junction. The phylogenetic analyses consistently recovered the familial relationships of (((((Ricaniidae + Issidae) + Flatidae) + Fulgoridae) + Achilidae) + Derbidae) in the amino acid-based analysis, with the placement of Cixiidae and Delphacidae as the earliest-derived lineages of fulgoroid families, whereas the monophyly of Delphacidae was not congruent between tree-constructing algorithms.

Tunicates Illuminate the Enigmatic Evolution of Chordate Metallothioneins by Gene Gains and Losses, Independent Modular Expansions, and Functional Convergences

To investigate novel patterns and processes of protein evolution, we have focused in the metallothioneins (MTs), a singular group of metal-binding, cysteine-rich proteins that, due to their high degree of sequence diversity, still represents a "black hole" in Evolutionary Biology. We have identified and analyzed more than 160 new MTs in nonvertebrate chordates (especially in 37 species of ascidians, 4 thaliaceans, and 3 appendicularians) showing that prototypic tunicate MTs are mono-modular proteins with a pervasive preference for cadmium ions, whereas vertebrate and cephalochordate MTs are bimodular proteins with diverse metal preferences. These structural and functional differences imply a complex evolutionary history of chordate MTs-including de novo emergence of genes and domains, processes of convergent evolution, events of gene gains and losses, and recurrent amplifications of functional domains-that would stand for an unprecedented case in the field of protein evolution.

Phylogenomic and comparative genomic analyses of species of the family Pseudomonadaceae: Proposals for the genera Halopseudomonas gen. nov. and Atopomonas gen. nov., merger of the genus Oblitimonas with the genus Thiopseudomonas, and transfer of some misclassified species of the genus Pseudomonas into other genera

The evolutionary relationships among species of the family Pseudomonadaceae were examined based on 255 available genomes representing >85 % of the species from this family. In a phylogenetic tree based on concatenated sequences of 118 core proteins, most species of the genus Pseudomonas grouped within one large cluster which also included members of the genera Azotobacter and Azomonas. Within this large cluster 18-30 clades/subclades of species of the genus Pseudomonas consisting of between 1 and 36 species, were observed. However, a number of species of the genus Pseudomonas branched outside of this main cluster and were interspersed among other genera of the family Pseudomonadaceae. This included a strongly supported clade (Pertucinogena clade) consisting of 19 mainly halotolerant species. The distinctness of this clade from all other members of the family Pseudomonadaceae is strongly supported by 24 conserved signature indels (CSIs) in diverse proteins that are exclusively found in all members of this clade. Nine uncharacterized members of the genus Pseudomonas also shared these CSIs and they branched within the Pertucinogena clade, indicating their affiliation to this clade. On the basis of the strong evidence supporting the distinctness of the Pertucinogena clade, we are proposing transfer of species from this clade into a novel genus Halopseudomonas gen. nov. Pseudomonas caeni also branches outside of the main cluster and groups reliably with Oblitimonas alkaliphila and Thiopseudomonas denitrificans. Six identified CSIs are uniquely shared by these three species and we are proposing their integration into the emended genus Thiopseudomonas, which has priority over the name Oblitimonas. We are also proposing transfer of the deep-branching Pseudomonas hussainii, for which 22 exclusive CSIs have been identified, into the genus Atopomonas gen. nov. Lastly, we present strong evidence that the species Pseudomonas cissicola and Pseudomonas geniculata are misclassified into the genus Pseudomonas and that they are specifically related to the genera Xanthomonas and Stenotrophomonas, respectively. In addition, we are also reclassifying 'Pseudomonas acidophila' as Paraburkholderia acidicola sp. nov. (Type strain: G-6302=ATCC 31363=BCRC 13035).