Maximum Likelihood Estimation for Unrooted 3-Leaf Trees: An Analytic Solution for the CFN Model.

Maximum likelihood estimation is among the most widely-used methods for inferring phylogenetic trees from sequence data. This paper solves the problem of computing solutions to the maximum likelihood problem for 3-leaf trees under the 2-state symmetric mutation model (CFN model). Our main result is a closed-form solution to the maximum likelihood problem for unrooted 3-leaf trees, given generic data; this result characterizes all of the ways that a maximum likelihood estimate can fail to exist for generic data and provides theoretical validation for predictions made in Parks and Goldman (Syst Biol 63(5):798-811, 2014). Our proof makes use of both classical tools for studying group-based phylogenetic models such as Hadamard conjugation and reparameterization in terms of Fourier coordinates, as well as more recent results concerning the semi-algebraic constraints of the CFN model. To be able to put these into practice, we also give a complete characterization to test genericity.

An entropy-based study on the mutational landscape of SARS-CoV-2 in USA: Comparing different variants and revealing co-mutational behavior of proteins.

COVID-19 emergency has pushed the international scientific community to use every resource to combat the spread of the virus, to understand its biology and predict its possible evolution in terms of new variants. Since the first SARS-CoV-2 virus nucleotide and amino acid sequences were made available, information theory was used to study how viral information content was changing over time and then trace the evolution of its mutational landscape. In this work we analyzed SARS-CoV-2 sequences collected mainly in the USA in a period from March 2020 until December 2022 and computed mutation profiles of viral proteins over time through an entropy-based approach using Shannon Entropy and Hellinger distance. This representation allows an at-a-glance view of the mutational landscape of viral proteins over time and can provide new insights on the evolution of the virus from different points of view. Non-structural proteins typically showed flat mutation profiles, characterized by a very low Average mutation Entropy, while accessory and structural proteins showed mostly non uniform and high mutation profiles, often coupled with the predominance of variants. Interestingly NSP2 protein, whose function is currently still debated, falls in the same branch of NSP14 and NSP10 in the phylogenetic tree of mutations constructed through correlations of mutation profiles, suggesting a co-evolution of those proteins and a possible functional link with each other. To the best of our knowledge this is the first study based on a massive amount of data (n = 107,939,973) that analyzes from an entropy point of view the mutational landscape of SARS-CoV-2 over time and depicts a mutational temporal profile of each protein of the virus.

Ocean to Tree: Leveraging Single-Molecule RNA-Seq to Repair Genome Gene Models and Improve Phylogenomic Analysis of Gene and Species Evolution.

Understanding gene evolution across genomes and organisms, including ctenophores, can provide unexpected biological insights. It enables powerful integrative approaches that leverage sequence diversity to advance biomedicine. Sequencing and bioinformatic tools can be inexpensive and user-friendly, but numerous options and coding can intimidate new users. Distinct challenges exist in working with data from diverse species but may go unrecognized by researchers accustomed to gold-standard genomes. Here, we provide a high-level workflow and detailed pipeline to enable animal collection, single-molecule sequencing, and phylogenomic analysis of gene and species evolution. As a demonstration, we focus on (1) PacBio RNA-seq of the genome-sequenced ctenophore Mnemiopsis leidyi, (2) diversity and evolution of the mechanosensitive ion channel Piezo in genetic models and basal-branching animals, and (3) associated challenges and solutions to working with diverse species and genomes, including gene model updating and repair using single-molecule RNA-seq. We provide a Python Jupyter Notebook version of our pipeline (GitHub Repository: Ctenophore-Ocean-To-Tree-2023 https://github.com/000generic/Ctenophore-Ocean-To-Tree-2023 ) that can be run for free in the Google Colab cloud to replicate our findings or modified for specific or greater use. Our protocol enables users to design new sequencing projects in ctenophores, marine invertebrates, or other novel organisms. It provides a simple, comprehensive platform that can ease new user entry into running their evolutionary sequence analyses.

Detection and Phylogenetic Analysis of Caprine Arthritis Encephalitis Virus Using TaqMan-based qPCR in Eastern China.

Caprine arthritis encephalitis is an infectious disease caused by the caprine arthritis encephalitis virus that infects goats, sheep, and other small ruminants. An outbreak of CAEV could be extremely harmful to the goat farming industry and could cause severe economic losses. We designed specific primers and probes for the gag gene and established a TaqMan real-time quantitative polymerase chain reaction assay. This method's correlation coefficient (R2) was >0.999, and the sensitivity of the assay to the plasmid-carried partial gag gene was approximately 10 copies/µL, 1000 times higher than that of conventional PCR. No specific fluorescence was detected for other sheep viruses. Using this method, we tested 776 asymptomatic sheep blood samples and 4 neurodegenerative sheep brain samples from six farms in eastern China, and the positivity rate was 0.77% (6/780). The gag gene was partially sequenced in the three positive samples and compared with the sequences from other representative strains in GenBank. The results revealed that all three strains belonged to the B1 subtype and were most closely related to the strains from Shanxi and Gansu, previously isolated in China, with their homology ranging from 97.7% to 98.9%. These results suggest that the designed RT-qPCR assay can be used to detect subclinical CAEV in sheep and that the virus is still present in eastern China.

Choice of Metric Divergence in Genome Sequence Comparison.

The paper introduces a novel probability descriptor for genome sequence comparison, employing a generalized form of Jensen-Shannon divergence. This divergence metric stems from a one-parameter family, comprising fractions up to a maximum value of half. Utilizing this metric as a distance measure, a distance matrix is computed for the new probability descriptor, shaping Phylogenetic trees via the neighbor-joining method. Initial exploration involves setting the parameter at half for various species. Assessing the impact of parameter variation, trees drawn at different parameter values (half, one-fourth, one-eighth). However, measurement scales decrease with parameter value increments, with higher similarity accuracy corresponding to lower scale values. Ultimately, the highest accuracy aligns with the maximum parameter value of half. Comparative analyses against previous methods, evaluating via Symmetric Distance (SD) values and rationalized perception, consistently favor the present approach's results. Notably, outcomes at the maximum parameter value exhibit the most accuracy, validating the method's efficacy against earlier approaches.

Genomic and phylogenomic analysis of Fusobacteriaceae family and proposal to reclassify Fusobacterium naviforme Jungano 1909 into a novel genus as Zandiella naviformis gen. nov., comb. nov. and reclassification of Fusobacterium necrophorum subsp. funduliforme as later heterotypic synonym of Fusobacterium necrophorum subsp. necrophorum and Fusobacterium equinum as later heterotypic synonym of Fusobacterium gonidiaformans.

The family Fusobacteriaceae is a large family within the phylum Fusobacteriota. The reclassification of F. naviforme as Zandiella naviformis gen. nov., comb. nov. is proposed because of the separate and distinct phylogenetic situation on the basis of the results of 16S rRNA gene sequence analysis, the genetic and genomic differences from all other species and subspecies in the Fusobacteriaceae family. The type strain is ATCC 25832; CCUG 50052; NCTC 13121. In phylogenetic trees drawn using complete genome sequences and 16S rRNA gene sequences, F. necrophorum subsp. funduliforme and F. equinum were clades together with F. necrophorum subsp. necrophorum and F. gonidiaformans, respectively. The average nucleotide identity, average amino acid identity, and digital DNA-DNA hybridization values between themes exceeded the cut-off values for species delineation. Based on these results, F. necrophorum subsp. funduliforme and F. equinum should be reclassified as later heterotypic synonyms of F. necrophorum subsp. necrophorum and F. gonidiaformans, respectively.

NewickTreeModifier: A simple web tool to prune and modify Newick trees.

Large-scale selection analyses of protein-coding sequences and phylogenetic tree reconstructions require suitable trees in Newick format. We developed the NewickTreeModifier (NTM), a simple web-based tool to trim and modify Newick trees for such analyses. The users can choose provided master trees or upload a tree to prune it to selected species available in FASTA, NEXUS, or PHYLIP sequence format with an internal converter, a simple species list, or directly determined from a checklist interface of the master trees. Plant, insect, and vertebrate master trees comprise the maximum number of species in an up-to-date phylogenetic order directly transferable to the pruned Newick outfile. NTM is available at https://retrogenomics.uni-muenster.de/tools/ntm.

Molecular Characterization and Pathogenicity of an Infectious cDNA Clone of Youcai Mosaic Virus on Solanum nigrum.

Virus infections cause devastative economic losses for various plant species, and early diagnosis and prevention are the most effective strategies to avoid the losses. Exploring virus genomic evolution and constructing virus infectious cDNA clones is essential to achieve a deeper understanding of the interaction between host plant and virus. Therefore, this work aims to guide people to better prevent, control, and utilize the youcai mosaic virus (YoMV). Here, the YoMV was found to infect the Solanum nigrum under natural conditions. Then, an infectious cDNA clone of YoMV was successfully constructed using triple-shuttling vector-based yeast recombination. Furthermore, we established phylogenetic trees based on the complete genomic sequences, the replicase gene, movement protein gene, and coat protein gene using the corresponding deposited sequences in NCBI. Simultaneously, the evolutionary relationship of the YoMV discovered on S. nigrum to others was determined and analyzed. Moreover, the constructed cDNA infectious clone of YoMV from S. nigrum could systematically infect the Nicotiana benthamiana and S. nigrum by agrobacterium-mediated infiltration. Our investigation supplied a reverse genetic tool for YoMV study, which will also contribute to in-depth study and profound understanding of the interaction between YoMV and host plant.

TreeViewer: Flexible, modular software to visualise and manipulate phylogenetic trees.

Phylogenetic trees illustrate evolutionary relationships between taxa or genes. Tree figures are crucial when presenting results and data, and by creating clear and effective plots, researchers can describe many kinds of evolutionary patterns. However, producing tree plots can be a time-consuming task, especially as multiple different programs are often needed to adjust and illustrate all data associated with a tree. We present TreeViewer, a new software to draw phylogenetic trees. TreeViewer is flexible, modular, and user-friendly. Plots are produced as the result of a user-defined pipeline, which can be finely customised and easily applied to different trees. Every feature of the program is documented and easily accessible, either in the online manual or within the program's interface. We show how TreeViewer can be used to produce publication-ready figures, saving time by not requiring additional graphical post-processing tools. TreeViewer is freely available for Windows, macOS, and Linux operating systems and distributed under an AGPLv3 licence from https://treeviewer.org. It has a graphical user interface (GUI), as well as a command-line interface, which is useful to work with very large trees and for automated pipelines. A detailed user manual with examples and tutorials is also available. TreeViewer is mainly aimed at users wishing to produce highly customised, publication-quality tree figures using a single GUI software tool. Compared to other GUI tools, TreeViewer offers a richer feature set and a finer degree of customisation. Compared to command-line-based tools and software libraries, TreeViewer's graphical interface is more accessible. The flexibility of TreeViewer's approach to phylogenetic tree plotting enables the program to produce a wide variety of publication-ready figures. Users are encouraged to create their own custom modules to expand the functionalities of the program. This sets the scene for an ever-expanding and ever-adapting software framework that can easily adjust to respond to new challenges.

The k-Robinson-Foulds Dissimilarity Measures for Comparison of Labeled Trees.

Understanding the mutational history of tumor cells is a critical endeavor in unraveling the mechanisms that drive the onset and progression of cancer. Modeling tumor cell evolution with labeled trees motivates researchers to develop different measures to compare labeled trees. Although the Robinson-Foulds (RF) distance is widely used for comparing species trees, its applicability to labeled trees reveals certain limitations. This study introduces the k-RF dissimilarity measures, tailored to address the challenges of labeled tree comparison. The RF distance is succinctly expressed as n-RF in the space of labeled trees with n nodes. Like the RF distance, the k-RF is a pseudometric for multiset-labeled trees and becomes a metric in the space of 1-labeled trees. By setting k to a small value, the k-RF dissimilarity can capture analogous local regions in two labeled trees with different size or different labels.

Maximum likelihood estimation and natural pairwise estimating equations are identical for three sequences and a symmetric 2-state substitution model.

Consider the problem of estimating the branch lengths in a symmetric 2-state substitution model with a known topology and a general, clock-like or star-shaped tree with three leaves. We show that the maximum likelihood estimates are analytically tractable and can be obtained from pairwise sequence comparisons. Furthermore, we demonstrate that this property does not generalize to larger state spaces, more complex models or larger trees. Our arguments are based on an enumeration of the free parameters of the model and the dimension of the minimal sufficient data vector. Our interest in this problem arose from discussions with our former colleague Freddy Bugge Christiansen.

Predicting Impacts of Contact Tracing on Epidemiological Inference from Phylogenetic Data.

Robust sampling methods are foundational to many inference problems in the phylodynamic field, yet the impact of using contact tracing, a type of non-uniform sampling used in public health applications, is not well understood. To investigate and quantify how this non-uniform sampling method influences recovered phylogenetic tree structure, we developed a new simulation tool called SEEPS (Sequence Evolution and Epidemiological Process Simulator) that allows for the simulation of contact tracing and the resulting transmission tree, pathogen phylogeny, and corresponding virus genetic sequences. Importantly, SEEPS takes within-host evolution into account when generating pathogen phylogenies and sequences from transmission histories. Using SEEPS, we demonstrate that contact tracing can significantly impact the structure of the resulting tree as described by popular tree statistics. Contact tracing generates phylogenies that are less balanced than the underlying transmission process, less representative of the larger epidemiological process, and affects the internal/external branch length ratios that characterize specific epidemiological scenarios. We also examine a 2007-2008 Swedish HIV-1 outbreak and the broader 1998-2010 European HIV-1 epidemic to highlight the differences in contact tracing and expected phylogenies. Aided by SEEPS, we show that the Swedish outbreak was strongly influenced by contact tracing even after downsampling, while the broader European Union epidemic showed little evidence of universal contact tracing, agreeing with the known epidemiological information about sampling and spread. SEEPS is available at github.com/MolEvolEpid/SEEPS.

First Report of the Occurrence of Diaporthe nobilis Causing Postharvest Fruit Rot on 'Shine Muscat' Grape in Yuncheng, China.

In recent years, 'Shine Muscat' (Vitis labruscana × Vitis vinifera) has been the most popular table grape among consumers in Asian countries (Choi et al., 2021). In November 2020, fruit rot was observed in packing house 2 months after cold-stored in Yuncheng (35°16'N, 111°08'E), Shanxi Province, China. Of 38 boxes (5 kg per box) grape were selected randomly from the packing house with a stock of 7,000 kg. The incidence of grape rot was 71.5%, observed in 143 of the 200 clusters sampled. The diseased fruit showed symptoms including pedicel browning, and extending lesions and softening of berries. Samples (1 mm2) cut from the border between diseased and healthy tissue were surface sterilized with 70% ethanol for 1 min, then rinsed with sterile distilled water thrice, and placed on potato dextrose agar (PDA) plate 25 ° C in the dark for 5 days. After subculturing for three times, 26 isolates with a similar type of fungal colony were obtained. Colonies on PDA was initially white fluffy aerial, but gradually turned grey. Pycnidia were black and globes, producing α and ß conidia. α conidia abundant, aseptate, hyaline, ellipsoid, both ends obtuse, biguttulate, 6.52~9.05 µm×2.19~3.25 µm (n=50). ß conidia abundant, hyaline, aseptate, filiform, flexuous to J-shaped, 22.79~36.75 µm×1.24~1.92 µm (n=50). For further identification, phylogenetic analysis of combined internal transcribed spacer (ITS), ß-tubulin (TUB), translation elongation factor (EF1-α) and calmodulin (CAL) datasets were carried out. The representative isolate SS1 and SS2 were selected for identification. Primers ITS1 and ITS4 (White et al. 1990) were used to amplify the ITS region. The gene fragments of TUB, CAL and EF1-α were amplified by the primers Bt2a/Bt2b (Glass and Donaldson 1995), CAL228F/CAL737R and EF1-728F/EF1-986R, respectively (Carbone and Kohn 1999). The sequences were deposited in GenBank (SS1: MW644526 for ITS, OQ718912 - 4 for CAL, TEF-1α, and TUB, respectively; SS2: MW644527 for ITS, OQ718915 -7 for CAL, TEF-1α, and TUB, respectively). BLASTn analysis of the NCBI database indicated that the ITS, TUB, EF1-α, and CAL sequences had very high nucleotide homology (98 to 100%) to ex-type sequences for D. nobilis (CBS 116953, Gomes et al., 2013) in GenBank (ITS: KC343147, TUB: KC344115, EF1-α: KC343873, and CAL: KC343389). Phylogenetic tree based on the combined sequences revealed that the two isolates clustered well with D. nobilis. To confirm its pathogenicity, 30 fresh-harvested 'Shine Muscat' grape fruits were used to mock inoculate. Wounds (5 mm wide by 2 mm deep) were created uniformly at the equator of the grape fruit with a sterile puncher and inoculated with mycelial plugs (5 mm in diameter) prepared from 7-day-old fungal cultures grown on PDA (Wild 1994). The control fruits were inoculated with agar disks without mycelium. Inoculated fruits were placed in a chamber with the 80% relative humidity at 25°C for 8 days. After inoculation, all treated fruits turned brown, with significant necrotic lesions; fruits in control showed no symptoms. Fungal colonies reisolated from inoculated fruits had similar morphological characteristics as D. nobilis. Diaporthe spp. are responsible for diseases on a wide range of plants hosts, causing root and fruit rots, dieback, cankers, leaf spots, blights, decay and wilt (Gomes et al., 2013). D. nobilis causing postharvest grape rot can seriously restrict the development of local grape industry.

An Unsupervised Classifier for Whole-Genome Phylogenies, the Maxwell© Tool.

The development of phylogenetic trees based on RNA or DNA sequences generally requires a precise and limited choice of important RNAs, e.g., messenger RNAs of essential proteins or ribosomal RNAs (like 16S), but rarely complete genomes, making it possible to explain evolution and speciation. In this article, we propose revisiting a classic phylogeny of archaea from only the information on the succession of nucleotides of their entire genome. For this purpose, we use a new tool, the unsupervised classifier Maxwell, whose principle lies in the Burrows-Wheeler compression transform, and we show its efficiency in clustering whole archaeal genomes.

Contribution to the Knowledge of Gastrointestinal Nematodes in Roe Deer (Capreolus capreolus) from the Province of León, Spain: An Epidemiological and Molecular Study.

A study of gastrointestinal nematodes in roe deer was carried out in the regional hunting reserves of Riaño and Mampodre, Province of León, Spain, to provide information on their prevalence and intensity of infection in relation to the sampling areas, age of the animals, and body weight. Through a regulated necropsy of the animals, all of them harbored gastrointestinal nematodes in their digestive tract, with a mean intensity of parasitism of 638 ± 646.1 nematodes/infected animal. Eleven genera were found and 18 species of gastrointestinal nematodes were identified, three of them polymorphic: Trichostrongylus axei, Trichostrongylus vitrinus, Trichostrongylus capricola, Trichostrongylus colubriformis, Haemonchus contortus, Spiculopteragia spiculoptera/Spiculopteragia mathevossiani, Ostertagia leptospicularis/Ostertagia kolchida, Ostertagia (Grosspiculopteragia) occidentalis, Teladorsagia circumcincta/Teladorsagia trifurcate, Marshallagia marshalli, Nematodirus europaeus, Cooperia oncophora, Capillaria bovis, Oesophagostomum venulosum, and Trichuris ovis. All of them have already been cited in roe deer in Europe, but Marshallagia marshalli, Capillaria bovis, and Ostertagia (Grosspiculopteragia) occidentalis are reported for the first time in Spain in this host. The abomasum was the intestinal section, where the prevalence (98.9%) and mean intensity (x¯ = 370.7 ± 374.4 worms/roe deer; range 3-1762) were significantly higher, but no statistically significant differences were found when comparing the sampling areas and age of animals. The animals with lower body weight had a higher parasite load than those in better physical condition, finding, in this case, statistically significant differences (p = 0.0020). Seven genera and 14 species were identified. In the small intestine, 88% of the animals examined presented gastrointestinal nematodes, with an average intensity of x¯ = 131.7 ± 225.6 parasites/infected animal, ranging between 4-1254 worms. No statistically significant differences were found when the three parameters studied were compared. Four genera and seven species were identified. In the large intestine/cecum, 78.3% of the examined roe deer presented adult worms, with an average intensity of 6.3 ± 5.5 worms/infected animal; range 1-26 worms. Only statistically significant differences were observed when considering the mean intensity of parasitism and the sampling area (p = 0.0093). Two genera and two species were identified. Several of the species found in the study were studied molecularly, and with the sequences obtained compared with those deposited in GenBank, phylogenetic trees were prepared to determine their taxonomic status. Using coprological techniques, the existing correlation in the shedding of gastrointestinal nematode eggs in roe deer was investigated with that of semi-extensive sheep farms in the same study area to verify the existence of cross-transmission of these parasites between wild and domestic animals. The high values found in the studied parameters show that northern Spain is an area of high-intensity infection for roe deer.

Robustness of Felsenstein's Versus Transfer Bootstrap Supports With Respect to Taxon Sampling.

The bootstrap method is based on resampling sequence alignments and re-estimating trees. Felsenstein's bootstrap proportions (FBP) are the most common approach to assess the reliability and robustness of sequence-based phylogenies. However, when increasing taxon sampling (i.e., the number of sequences) to hundreds or thousands of taxa, FBP tend to return low support for deep branches. The transfer bootstrap expectation (TBE) has been recently suggested as an alternative to FBP. TBE is measured using a continuous transfer index in [0,1] for each bootstrap tree, instead of the binary {0,1} index used in FBP to measure the presence/absence of the branch of interest. TBE has been shown to yield higher and more informative supports while inducing a very low number of falsely supported branches. Nonetheless, it has been argued that TBE must be used with care due to sampling issues, especially in datasets with a high number of closely related taxa. In this study, we conduct multiple experiments by varying taxon sampling and comparing FBP and TBE support values on different phylogenetic depths, using empirical datasets. Our results show that the main critique of TBE stands in extreme cases with shallow branches and highly unbalanced sampling among clades, but that TBE is still robust in most cases, while FBP is inescapably negatively impacted by high taxon sampling. We suggest guidelines and good practices in TBE (and FBP) computing and interpretation.

The complete chloroplast genome of Aegle marmelos and its phylogenetic analysis.

Aegle marmelos (L.) Correa 1800, a plant belonging to the Rutaceae family, is extensively used in Tibetan medicine. We employed Illumina HiSeq reads to assemble the complete chloroplast (cp) genome of A. marmelos, which spans 144,538 bp. The genome comprises 114 genes, including 75 protein-coding genes, 31 tRNA genes, and 8 rRNA genes. It is characterized by four regions: The large single-copy (LSC) region (74,253 bp), the inverted repeat A (IRa) region (26,015 bp), the small single-copy (SSC) region (18,255 bp), and the inverted repeat B (IRb) region (26,015 bp). Phylogenomic analysis demonstrated a close relationship between A. marmelos and Citrus. The assembly of The cp genome in this study serves as a foundation for conservation efforts and phylogenetic investigations of A. marmelos, paving the way for future experimentation.

Sequencing and analysis of the complete mitochondrial genome of Eothenomys eleusis Thomas 1911 from China and its phylogenetic analysis.

The complete mitogenome sequence of Eothenomys eleusis Thomas 1911 was determined using PCR. A circular double-stranded structure makes up the mitochondrial genome of E. eleusis. The complete length of the mitochondrial genome is 16,419 bp. The mitochondrial genome of E. eleusis included 13 protein-coding genes, 1 control region, 22 tRNA genes, 2 rRNA genes and 1 origin of L strand replication. The total base composition of E. eleusis mitochondrial genome was A (32.6%), T (26.3%), G (13.6%) and C (27.5%). We found significant A-T skew in base composition, especially in control regions and protein-coding genes. E. eleusis was supported by bootstrap values of 100%. This study verifies the evolutionary status of E. eleusis in Myodini tribe of Cricetidae at the molecular level. The mitochondrial genome would be a significant supplement for the E. eleusis genetic background.

The bacterial origin of mitochondria: Incorrect phylogenies and the importance of metabolic traits.

This article provides an updated review on the evolution of mitochondria from bacteria, which were likely related to extant alphaproteobacteria. Particular attention is given to the timeline of oxygen history on Earth and the entwined phases of eukaryotic evolution that produced the animals that still populate our planet. Mitochondria of early-branching unicellular eukaryotes and plants appear to retain partial or vestigial traits that were directly inherited from the alphaproteobacterial ancestors of the organelles. Most of such traits define the current aerobic physiology of mitochondria. Conversely, the anaerobic traits that would be essential in the syntrophic associations postulated for the evolution of eukaryotic cells are scantly present in extant alphaproteobacteria, and therefore cannot help defining from which bacterial lineage the ancestors of mitochondria originated. This question has recently been addressed quantitatively, reaching the novel conclusion that marine bacteria related to Iodidimonas may be the living relatives of protomitochondria. Additional evidence is presented that either support or does not contrast this novel view of the bacterial origin of mitochondria.

Topological analysis of vertebrate phylogenetic trees utilizing Horton's first law.

Phylogenetic trees visually represent evolution and diversification. While many studies have focused on the number and length of edges (branches), topological properties, such as edge connection patterns, are also important. In this study, the topological properties of phylogenetic trees were quantified, focusing on edge connection patterns. Horton's first law was applied to quantify the overall, rather than local, topological properties of phylogenetic trees. The topological properties of vertebrate phylogenetic trees for spiny-rayed fishes, Amphibians, turtles, Squamata, Aves, and placental mammals were analyzed. The topological features discussed herein include the number of first-order edges, maximum order, and bifurcation ratio. The average bifurcation ratio of all trees was approximately 3, suggesting that phylogenetic trees for different taxa have a common mechanism of evolution. Vertebrate phylogenetic trees were compared with artificial branching objects created from neutral stochastic branching model simulations. The topological properties of the actual vertebrate phylogenetic trees agreed with those of the artificial branching objects. Our study suggests that evolutionary events do not change the overall topological properties of actual phylogenetic trees, even if the number and length of the edges change. Specifically, non-neutral events (e.g., environmental changes and mass extinction) are not main factors associated with topological properties. The results instead demonstrate a relationship between the bifurcation ratio and symmetricity in the context of temporal changes of topological properties. When the number of first-order edges increased and the maximum order remained constant, the bifurcation ratio increased and symmetricity decreased. When the number of first-order edges increased and the maximum order increased by one, the bifurcation ratio decreased and symmetricity increased.