Tree Reconstruction Research Articles

Reproducible Model Selection Using Bagged Posteriors.

Bayesian model selection is premised on the assumption that the data are generated from one of the postulated models. However, in many applications, all of these models are incorrect (that is, there is misspecification). When the models are misspecified, two or more models can provide a nearly equally good fit to the data, in which case Bayesian model selection can be highly unstable, potentially leading to self-contradictory findings. To remedy this instability, we propose to use bagging on the posterior distribution ("BayesBag") - that is, to average the posterior model probabilities over many bootstrapped datasets. We provide theoretical results characterizing the asymptotic behavior of the posterior and the bagged posterior in the (misspecified) model selection setting. We empirically assess the BayesBag approach on synthetic and real-world data in (i) feature selection for linear regression and (ii) phylogenetic tree reconstruction. Our theory and experiments show that, when all models are misspecified, BayesBag (a) provides greater reproducibility and (b) places posterior mass on optimal models more reliably, compared to the usual Bayesian posterior; on the other hand, under correct specification, BayesBag is slightly more conservative than the usual posterior, in the sense that BayesBag posterior probabilities tend to be slightly farther from the extremes of zero and one. Overall, our results demonstrate that BayesBag provides an easy-to-use and widely applicable approach that improves upon Bayesian model selection by making it more stable and reproducible.

Serratia silvae sp. nov., Isolated from Forest Soil.

The Gram-negative, oxidase-negative, catalase-positive, rod-shaped strain Arafor3T was isolated from forest soil (France). Comparative 16S rRNA gene analysis and phylogenetic analysis based on (1) multilocus sequence analysis (MLSA) with four housekeeping genes (atpD, gyrB, infB and rpoB) and (2) genomes indicated that strain Arafor3T shared 98.83% 16S rRNA gene sequence similarity with the type strain of Serratia fonticola DSM 4576T and was closely related to this same strain in the MLSA and in the phylogenomic tree reconstruction. Average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) comparisons of strain Arafor3T with its nearest neighbor S. fonticola DSM 4576T showed 93.5% identity and 55.7% sequence similarity, respectively, and were lower than the 96% and 70% species-level cut-off values relating to these analyses (Logan et al. in Int J Syst Evol Microbiol 59:2114-21, 2009, https://doi.org/10.1099/ijs.0.013649-0 ). The strain differed from S. fonticola in that it was urease and arginine dihydrolase negative. The major fatty acids of strain Arafor3T are C16:0, C16:1 ω7c/C16:1 ω6c, C14:0, C14:0 3-OH/16:1 isoI, and C18:1 ω7c. The major respiratory quinone is Q8. The polar lipid profile consisted of diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol, and 6 unknown lipids. The mol G + C% content of the genomic DNA of strain Arafor3T was 53.49%. Hence, Arafor3T represents a novel species within the genus Serratia, for which the name Serratia silvae sp. nov. is proposed. The type strain is Arafor3T (=LMG 32338T = CIP 111939T).

Reconstructing clonal tree for phylo-phenotypic characterization of cancer using single-cell transcriptomics

Functional characterization of the cancer clones can shed light on the evolutionary mechanisms driving cancer’s proliferation and relapse mechanisms. Single-cell RNA sequencing data provide grounds for understanding the functional state of cancer as a whole; however, much research remains to identify and reconstruct clonal relationships toward characterizing the changes in functions of individual clones. We present PhylEx that integrates bulk genomics data with co-occurrences of mutations from single-cell RNA sequencing data to reconstruct high-fidelity clonal trees. We evaluate PhylEx on synthetic and well-characterized high-grade serous ovarian cancer cell line datasets. PhylEx outperforms the state-of-the-art methods both when comparing capacity for clonal tree reconstruction and for identifying clones. We analyze high-grade serous ovarian cancer and breast cancer data to show that PhylEx exploits clonal expression profiles beyond what is possible with expression-based clustering methods and clear the way for accurate inference of clonal trees and robust phylo-phenotypic analysis of cancer.

Genome-scale angiosperm phylogenies based on nuclear, plastome, and mitochondrial datasets.

Angiosperms dominate the Earth's ecosystems and provide most of the basic necessities for human life. The major angiosperm clades comprise 64 orders, as recognized by the APG IV classification. However, the phylogenetic relationships of angiosperms remain unclear, as phylogenetic trees with different topologies have been reconstructed depending on the sequence datasets utilized, from targeted genes to transcriptomes. Here, we used currently available de novo genome data to reconstruct the phylogenies of 366 angiosperm species from 241 genera belonging to 97 families across 43 of the 64 orders based on orthologous genes from the nuclear, plastid, and mitochondrial genomes of the same species with compatible datasets. The phylogenetic relationships were largely consistent with previously constructed phylogenies based on sequence variations in each genome type. However, there were major inconsistencies in the phylogenetic relationships of the five Mesangiospermae lineages when different genomes were examined. We discuss ways to address these inconsistencies, which could ultimately lead to the reconstruction of a comprehensive angiosperm tree of life. The angiosperm phylogenies presented here provide a basic framework for further updates and comparisons. These phylogenies can also be used as guides to examine the evolutionary trajectories among the three genome types during lineage radiation.

A phylogenetic study of dengue virus in urban Vietnam shows long-term persistence of endemic strains.

Dengue virus (DENV) causes repeated outbreaks of disease in endemic areas, with patterns of local transmission strongly influenced by seasonality, importation via human movement, immunity, and vector control efforts. An understanding of how each of these interacts to enable endemic transmission (continual circulation of local virus strains) is largely unknown. There are times of the year when no cases are reported, often for extended periods of time, perhaps wrongly implying the successful eradication of a local strain from that area. Individuals who presented at a clinic or hospital in four communes in Nha Trang, Vietnam, were initially tested for DENV antigen presence. Enrolled positive individuals then had their corresponding household members invited to participate, and those who enrolled were tested for DENV. The presence of viral nucleic acid in all samples was confirmed using quantitative polymerase chain reaction, and positive samples were then whole-genome sequenced using an amplicon and target enrichment library preparation techniques and Illumina MiSeq sequencing technology. Generated consensus genome sequences were then analysed using phylogenetic tree reconstruction to categorise sequences into clades with a common ancestor, enabling investigations of both viral clade persistence and introductions. Hypothetical introduction dates were additionally assessed using a molecular clock model that calculated the time to the most recent common ancestor (TMRCA). We obtained 511 DENV whole-genome sequences covering four serotypes and more than ten distinct viral clades. For five of these clades, we had sufficient data to show that the same viral lineage persisted for at least several months. We noted that some clades persisted longer than others during the sampling time, and by comparison with other published sequences from elsewhere in Vietnam and around the world, we saw that at least two different viral lineages were introduced into the population during the study period (April 2017-2019). Next, by inferring the TMRCA from the construction of molecular clock phylogenies, we predicted that two of the viral lineages had been present in the study population for over a decade. We observed five viral lineages co-circulating in Nha Trang from three DENV serotypes, with two likely to have remained as uninterrupted transmission chains for a decade. This suggests clade cryptic persistence in the area, even during periods of low reported incidence.

Phocaeicola oris sp. nov., an anaerobic bacterium isolated from the saliva of a patient with oral squamous cell carcinoma.

A Gram-stain-negative or -positive, strictly anaerobic, non-spore-forming and pleomorphic bacterium (designated 14-104T) was isolated from the saliva sample of a patient with oral squamous cell carcinoma. It was an acid-tolerant neutralophilic mesophile, growing at between 20 and 40 °C (with optimum growth at 30 °C) and pH between pH 3.0 and 7.0 (with optimum growth at pH 6.0-7.0). It contained anteiso-C15 : 0 and C15 : 0 as the major fatty acids. The genome size of strain 14-104T was 2.98 Mbp, and the G+C content was 39.6 mol%. It shared <87 % 16S rRNA sequence similarity, <71 % orthologous average nucleotide identity, <76 % average amino acid identity and <68 %% of conserved proteins with its closest relative, Phocaeicola abscessus CCUG 55929T. Reconstruction of phylogenetic and phylogenomic trees revealed that strain 14-104T and P. abscessus CCUG 55929T were clustered as a distinct clade without any other terminal node. The phylogenetic and phylogenomic analyses along with physiological and chemotaxonomic data indicated that strain 14-104T represents a novel species in the genus Phocaeicola, for which the name Phocaeicola oris sp. nov. is proposed. The type strain is 14-104T (=BCRC 81305T= NBRC 115041T).

Aristophania vespae gen. nov., sp. nov., isolated from wasps, is related to Bombella and Oecophyllibacter, isolated from bees and ants.

Acetic acid bacteria (family Acetobacteraceae) are found in the gut of most insects. Two clades are currently recognized: Commensalibacter-Entomobacter and Bombella-Oecophyllibacter. The latter group is only found in hymenopteran insects and the described species have been isolated from bees and ants. In this study, two new strains DDB2-T1T (=KACC 21507T=LMG 31759T) and DM15PD (=CCM 9165=DSM 112731=KACC 22353=LMG 32454) were isolated from wasps collected in the Republic of Korea and Germany, respectively. Molecular and phenotypic analysis revealed that the strains are closely related, with 16S rRNA gene sequences showing 100 % identity and genomic average nucleotide identity (ANI) values ≥99 %. The closest related species based on type strain 16S rRNA gene sequences are Swingsia samuiensis, Acetobacter peroxydans, Bombella favorum and Bombella intestini (94.8-94.7% identity), whereas the closest related species based on type strain genome analysis are Saccharibacter floricola and Bombella intestini (ANI values of 68.8 and 68.2 %, respectively). The reconstruction of a phylogenomic tree based on 107 core proteins revealed that the branch leading to DDB2-T1T and DM15PD is localized between Oecophyllibacter and Saccharibacter-Bombella. Further genomic distance metrics such as ANI, percentage of conserved proteins and alignment fraction values were consistent with these strains belonging to a new genus. The key phenotypic characteristics were one MALDI-TOF-MS peak (m/z=4601.9±2.0) and the ability to produce acid from d-arabinose. Based on this polyphasic approach, including phylogenetics, phylogenomics, genome distance calculations, ecology and phenotypic characteristics, we propose to name the novel strains Aristophania vespae gen. nov., sp. nov., with the type strain DDB2-T1T (=KACC 21507T=LMG 31759T).

DNA barcode of seven species coral from Sepulu, Madura Island, Indonesia

Abstract. Insafitri, Nursalim N, Kholilah N, Kurniasih EM, Cahyani NKD, Nugraha WA, Ambariyanto A. 2023. DNA barcode of seven species coral from Sepulu, Madura Island, Indonesia. Biodiversitas 24: 317-323. Diversity of coral species is needed because it affects the disruption of coral reef ecosystems, nutrient recycling, coral reef growth and the habitat of biota that live in these ecosystems. The main constituents of coral reef ecosystems are coral animals that belong to the Anthozoa class. Sepulu waters (Madura Island), Bangkalan District, East Java Province, Indonesia, has a percentage of coral cover that is in the bad category, management of coral reefs in this area has not been the main focus and it is feared that the situation will get worse so efforts are needed to maintain coral reefs to be sustainable and better. This management effort requires an analysis of species identification. DNA barcoding of corals has never been done on corals of Sepulu waters. The purpose of this study is to find out the DNA of barcoding coral species on Sepulu waters. Genetic marker Cytochrome Oxidase I of the mitochondrial genome DNA (mtDNA) was used to analyze genetic diversity. Reconstruction of phylogenetic tree and genetic diversity were done by using software MEGA X. Research results showed that sample DBP011101 closely related to Dipsatraea maxima 99.55%, DBP011102 closely related to Porites rus 99.4%, DBP011103 closely related to Acropora hyacinthus 99.11%, DBP011104 closely related to Porites cf lichen 98.66%, DBP011106 closely related to Dipsastraea rotumana 95.37%, DBP011107 closely related to Porites horrisoni 98.81%, and DBP011108 closely related to Acropora valida 98.38%. The results of the study are the first coral barcode DNA on Sepulu waters that can help in determining coral species that are useful in managing the sustainability of coral resources, especially on Sepulu waters.

Individual Tree Reconstruction Based on Circular Truncated Cones From Portable LiDAR Scanner Data

The LiDAR scanning approach captures a large amount of tree point cloud information, including the accurate coordinate of points, local topology, and overall geometry. However, the complex tree structures, e.g., curvature and occlusion of branches, bring challenges to the 3-D tree reconstruction. In this letter, we propose an innovative solution for obtaining the complete skeletons of individual trees and 3-D structures for modeling using point clouds. First, we obtain individual trees from input street scene and segment tree into small successive pieces, with the centers of each piece serving as skeleton candidate points. Second, we conduct the interpolation based on the Euclidean distance and orientation yields the entire skeleton completely. Finally, a high-precision 3-D model of trees is constructed by cylindrically fitting the skeleton relying on the optimized circular truncated cones depending on the branch orientation and cylindrical curvature. Experiments on various trees demonstrate the high efficiency and effectiveness of our method. Our method achieves 98% accuracy and takes less than 1 min in the reconstruction. Compared with other methods, our method reduces the time by more than 95%.

Spatial evolution of human cultures inferred through Bayesian phylogenetic analysis.

Spatial distribution of human culture reflects both descent from the common ancestor and horizontal transmission among neighbouring populations. To analyse empirically documented geographical variations in cultural repertoire, we will describe a framework for Bayesian statistics in a spatially explicit model. To consider both horizontal transmission and mutation of the cultural trait in question, our method employs a network model in which populations are represented by nodes. Using algorithms borrowed from Bayesian phylogenetic analysis, we will perform a Markov chain Monte Carlo (MCMC) method to compute the posterior distributions of parameters, such as the rate of horizontal transmission and the mutation rates among trait variants, as well as the identity of trait variants in unobserved populations. Besides the inference of model parameters, our method enables the reconstruction of the genealogical tree of the focal trait, provided that the mutation rate is sufficiently small. We will also describe a heuristic algorithm to reduce the dimension of the parameter space explored in the MCMC method, where we simulate the coalescent process in the network of populations. Numerical examples show that our algorithms compute the posterior distribution of model parameters within a practical computation time, although the posterior distribution tends to be broad if we use uninformative priors.

Unraveling the secrets of Eclipta alba (L.) Hassk.: a comprehensive study of morpho-anatomy and DNA barcoding.

Safety regarding herbal products is very necessary; therefore, routine identification of raw materials should be performed to ensure that the raw materials used in pharmaceutical products are suitable for their intended use. In order for the identification-related data obtained to be accurate, the identification of various kinds of markers is also very necessary. The purpose of this study was to describe the characteristics of Eclipta alba (L.) Hassk. based on qualitative morpho-anatomical markers and quantitative DNA coding. The morphology of this plant has herbaceous habit with a taproot and a stem with branches that appear from the middle. Leaves are single type imperfectly arranged oppositely, lanceolatus, finely serrated on the edges, tapered at the base, pointed at the end, and have a pinnate and hairy leaf surface. The flowers consist of ray flowers and tube flowers with a cup shape. Meanwhile, in terms of anatomy, E. alba has aerenchyma, which are scattered in the cortex of the root and stem. In addition, there are anisocytic stomata, glandular trichomes, and non-glandural trichomes with an elongated shape accompanied by ornamentation found on the leaf epidermis. The results of sequence alignment and phylogenetic tree reconstruction show that the sample plants are closely related to species in the genus Eclipta.

Detection of Vibrio sp. on vaname shrimp (Litopanaeus vannamei) cultivation pond in Gorontalo, Indonesia

Abstract. Ulfiani P, Hasan AM, Retnowati Y. 2022. Detection of Vibrio sp. on vaname shrimp (Litopanaeus vannamei) cultivation pond in Gorontalo, Indonesia. Nusantara Bioscience 14: 340-346. Vibrio is a pathogenic bacteria that can attack cultured shrimp. Therefore, it is necessary to increase the intensity of water changes and efforts to control water quality to maintain dissolved oxygen levels and manage organic matter content in ponds to prevent the increase of pathogenic bacteria, especially Vibrio. In addition, water physicochemical factors can influence the presence and population of Vibrio bacteria in ponds. This study aims to determine the presence of Vibrio sp. in vaname shrimp culture ponds in Gorontalo, Indonesia. Descriptive quantitative analysis was used to analyze the data collected from this study. Shrimp pond water samples were obtained from shrimp ponds in Molosipat Village, West Popayato Sub-district, Pohuwato District, Gorontalo, Indonesia, using a purposive sampling technique with an interval of three days from the first water change to the last before the water change in one period (approximately 10 days). The water samples were subjected to physicochemical analysis, including temperature, pH, salinity, ammonia, nitrate, and dissolved oxygen. The results showed that there were Vibrio bacteria in the Thiosulfate Citrate BileSalt Sucrose Agar (TCBSA) medium with a population that tends to increase in the range of 1.0 × 103 - 2.7 × 104 CFU/mL. The isolation results showed that two isolates were named ISL1 and ISL2, with the colony morphology characters of ISL1 isolates being green and ISL2 isolates being yellow. The results of microscopic observations of the two isolates showed comma-shaped and gram-negative cells. Molecular identification using the 16S rRNA gene and reconstruction of the phylogenetic tree using the Neighbor-Joining algorithm showed that ISL1 isolate is very closely related to Vibrio parahaemolyticus strain 11-1 (relationship level 99%) and ISL2 is closely related to V. parahaemolyticus strain B2-1 (relationship level 86%).

The GEN-ERA toolbox: unified and reproducible workflows for research in microbial genomics.

Microbial culture collections play a key role in taxonomy by studying the diversity of their strains and providing well-characterized biological material to the scientific community for fundamental and applied research. These microbial resource centers thus need to implement new standards in species delineation, including whole-genome sequencing and phylogenomics. In this context, the genomic needs of the Belgian Coordinated Collections of Microorganisms were studied, resulting in the GEN-ERA toolbox. The latter is a unified cluster of bioinformatic workflows dedicated to both bacteria and small eukaryotes (e.g., yeasts). This public toolbox allows researchers without a specific training in bioinformatics to perform robust phylogenomic analyses. Hence, it facilitates all steps from genome downloading and quality assessment, including genomic contamination estimation, to tree reconstruction. It also offers workflows for average nucleotide identity comparisons and metabolic modeling. Nextflow workflows are launched by a single command and are available on the GEN-ERA GitHub repository (https://github.com/Lcornet/GENERA). All the workflows are based on Singularity containers to increase reproducibility. The toolbox was developed for a diversity of microorganisms, including bacteria and fungi. It was further tested on an empirical dataset of 18 (meta)genomes of early branching Cyanobacteria, providing the most up-to-date phylogenomic analysis of the Gloeobacterales order, the first group to diverge in the evolutionary tree of Cyanobacteria. The GEN-ERA toolbox can be used to infer completely reproducible comparative genomic and metabolic analyses on prokaryotes and small eukaryotes. Although designed for routine bioinformatics of culture collections, it can also be used by all researchers interested in microbial taxonomy, as exemplified by our case study on Gloeobacterales.

Tree Reconstruction Using Topology Optimisation

Generating accurate digital tree models from scanned environments is invaluable for forestry, agriculture, and other outdoor industries in tasks such as identifying fall hazards, estimating trees’ biomass and calculating traversability. Existing methods for tree reconstruction rely on sparse feature identification to segment a forest into individual trees and generate a branch structure graph, limiting their application to easily separable trees and uniform forests. However, the natural world is a messy place in which trees present with significant heterogeneity and are frequently encroached upon by the surrounding environment. We present a general method for extracting the branch structure of trees from point cloud data, which estimates the structure of trees by adapting the methods of structural topology optimisation to find the optimal material distribution to interpolate the input data. We present the results of this optimisation over a wide variety of scans, and discuss the benefits and drawbacks of this novel approach to tree structure reconstruction. Our method generates detailed and accurate tree structures, with a mean Surface Error (SE) of 15 cm over 13 diverse tree datasets.

A Method Based on Improved iForest for Trunk Extraction and Denoising of Individual Street Trees

Currently, the street tree resource survey using Mobile laser scanning (MLS) represents a hot spot around the world. Refined trunk extraction is an essential step for 3D reconstruction of street trees. However, due to scanning errors and the effects of occlusion by various types of features in the urban environment, street tree point cloud data processing has the problem of excessive noise. For the noise points that are difficult to remove using statistical methods in close proximity to the tree trunk, we propose an adaptive trunk extraction and denoising method for street trees based on an improved iForest (Isolation Forest) algorithm. Firstly, to extract the individual tree trunk points, the trunk and the crown are distinguished from the individual tree point cloud through point cloud slicing. Next, the iForest algorithm is improved by conducting automatic calculation of the contamination and further used to denoise the tree trunk point cloud. Finally, the method is validated with five datasets of different scenes. The results indicate that our method is robust and effective in extracting and denoising tree trunks. Compared with the traditional Statistical Outlier Removal (SOR) filter and Radius filter denoising methods, the denoising accuracy of the proposed method can be improved by approximately 30% for noise points close to tree trunks. Compared to iForest, the proposed method automatically calculates the contamination, improving the automation of the algorithm. Our method can provide more precise trunk point clouds for 3D reconstruction of street trees.

Description and Genome Characterization of Three Novel Fungal Strains Isolated from Mars 2020 Mission-Associated Spacecraft Assembly Facility Surfaces-Recommendations for Two New Genera and One Species.

National Aeronautics and Space Administration's (NASA) spacecraft assembly facilities are monitored for the presence of any bacteria or fungi that might conceivably survive a transfer to an extraterrestrial environment. Fungi present a broad and diverse range of phenotypic and functional traits to adapt to extreme conditions, hence the detection of fungi and subsequent eradication of them are needed to prevent forward contamination for future NASA missions. During the construction and assembly for the Mars 2020 mission, three fungal strains with unique morphological and phylogenetic properties were isolated from spacecraft assembly facilities. The reconstruction of phylogenetic trees based on several gene loci (ITS, LSU, SSU, RPB,&amp;nbsp;TUB, TEF1) using multi-locus sequence typing (MLST) and whole genome sequencing (WGS) analyses supported the hypothesis that these were novel species. Here we report the genus or species-level classification of these three novel strains via a polyphasic approach using phylogenetic analysis, colony and cell morphology, and comparative analysis of WGS. The strain FJI-L9-BK-P1 isolated from the Jet Propulsion Laboratory Spacecraft Assembly Facility (JPL-SAF) exhibited a putative phylogenetic relationship with the strain Aaosphaeria arxii CBS175.79 but showed distinct morphology and microscopic features. Another JPL-SAF strain, FJII-L3-CM-DR1, was phylogenetically distinct from members of the family Trichomeriaceae and exhibited morphologically different features from the genera Lithohypha and Strelitziana. The strain FKI-L1-BK-DR1 isolated from the Kennedy Space Center facility was identified as a member of Dothideomycetes&amp;nbsp;incertae sedis and is closely related to the family Kirschsteiniotheliaceae according to a phylogenetic analysis. The polyphasic taxonomic approach supported the recommendation for establishing two novel genera and one novel species. The names Aaosphaeria pasadenensis (FJI-L9-BK-P1 = NRRL 64424 = DSM 114621), Pasadenomyces melaninifex (FJII-L3-CM-DR1 = NRRL 64433 = DSM 114623), and Floridaphiala radiotolerans (FKI-L1-BK-DR1 = NRRL 64434 = DSM 114624) are proposed as type species. Furthermore, resistance to ultraviolet-C and presence of specific biosynthetic gene cluster(s) coding for metabolically active compounds are unique to these strains.

Stochastic Safety Radius on UPGMA

Unweighted Pair Group Method with Arithmetic Mean (UPGMA) is one of the most popular distance-based methods to reconstruct an equidistant phylogenetic tree from a distance matrix computed from an alignment of sequences. Since we use equidistant trees as gene trees for phylogenomic analyses under the multi-species coalescent model and since an input distance matrix computed from an alignment of each gene in a genome is estimated via the maximum likelihood estimators, it is important to conduct a robust analysis on UPGMA. Stochastic safety radius, introduced by Steel and Gascuel, provides a lower bound for the probability that a phylogenetic tree reconstruction method returns the true tree topology from a given distance matrix. In this article, we compute the stochastic safety radius of UPGMA for a phylogenetic tree with n leaves. Computational experiments show an improved gap between empirical probabilities estimated from random samples and the true tree topology from UPGMA, increasing confidence in phylogenic results.

Defining Binary Phylogenetic Trees Using Parsimony

Phylogenetic (i.e., leaf-labeled) trees play a fundamental role in evolutionary research. A typical problem is to reconstruct such trees from data like DNA alignments (whose columns are often referred to as characters), and a simple optimization criterion for such reconstructions is maximum parsimony. It is generally assumed that this criterion works well for data in which state changes are rare. In the present manuscript, we prove that each binary phylogenetic tree T with nge 20 k leaves is uniquely defined by the set A_k(T), which consists of all characters with parsimony score k on T. This can be considered as a promising first step toward showing that maximum parsimony as a tree reconstruction criterion is justified when the number of changes in the data is relatively small.

SPLACE: A tool to automatically SPLit, Align, and ConcatenatE genes for phylogenomic inference of several organisms.

The reconstruction of phylogenomic trees containing multiple genes is best achieved by using a supermatrix. The advent of NGS technology made it easier and cheaper to obtain multiple gene data in one sequencing run. When numerous genes and organisms are used in the phylogenomic analysis, it is difficult to organize all information and manually align the gene sequences to further concatenate them. This study describes SPLACE, a tool to automatically SPLit, Align, and ConcatenatE the genes of all species of interest to generate a supermatrix file, and consequently, a phylogenetic tree, while handling possible missing data. In our findings, SPLACE was the only tool that could automatically align gene sequences and also handle missing data; and, it required only a few minutes to produce a supermatrix FASTA file containing 83 aligned and concatenated genes from the chloroplast genomes of 270 plant species. It is an open-source tool and is publicly available at https://github.com/reinator/splace.

A scalable model for simulating multi-round antibody evolution and benchmarking of clonal tree reconstruction methods.

Affinity maturation (AM) of B cells through somatic hypermutations (SHMs) enables the immune system to evolve to recognize diverse pathogens. The accumulation of SHMs leads to the formation of clonal lineages of antibody-secreting b cells that have evolved from a common naïve B cell. Advances in high-throughput sequencing have enabled deep scans of B cell receptor repertoires, paving the way for reconstructing clonal trees. However, it is not clear if clonal trees, which capture microevolutionary time scales, can be reconstructed using traditional phylogenetic reconstruction methods with adequate accuracy. In fact, several clonal tree reconstruction methods have been developed to fix supposed shortcomings of phylogenetic methods. Nevertheless, no consensus has been reached regarding the relative accuracy of these methods, partially because evaluation is challenging. Benchmarking the performance of existing methods and developing better methods would both benefit from realistic models of clonal lineage evolution specifically designed for emulating B cell evolution. In this paper, we propose a model for modeling B cell clonal lineage evolution and use this model to benchmark several existing clonal tree reconstruction methods. Our model, designed to be extensible, has several features: by evolving the clonal tree and sequences simultaneously, it allows modeling selective pressure due to changes in affinity binding; it enables scalable simulations of large numbers of cells; it enables several rounds of infection by an evolving pathogen; and, it models building of memory. In addition, we also suggest a set of metrics for comparing clonal trees and measuring their properties. Our results show that while maximum likelihood phylogenetic reconstruction methods can fail to capture key features of clonal tree expansion if applied naively, a simple post-processing of their results, where short branches are contracted, leads to inferences that are better than alternative methods.