When did the ancestor of true bugs become stinky? Disentangling the phylogenomics of Hemiptera–Heteroptera

The phylogeny of true bugs (Hemiptera: Heteroptera), one of the most diverse insect groups in terms of morphology and ecology, has been the focus of attention for decades with respect to several deep nodes between the suborders of Hemiptera and the infraorders of Heteroptera. Here, we assembled a phylogenomic data set of 53 taxa and 3102 orthologous genes to investigate the phylogeny of Hemiptera–Heteroptera, and both concatenation and coalescent methods were used. A binode-control approach for data filtering was introduced to reduce the incongruence between different genes, which can improve the performance of phylogenetic reconstruction. Both hypotheses (Coleorrhyncha + Heteroptera) and (Coleorrhyncha + Auchenorrhyncha) received support from various analyses, in which the former is more consistent with the morphological evidence. Based on a divergence time estimation performed on genes with a strong phylogenetic signal, the origin of true bugs was dated to 290–268 Ma in the Permian, the time in Earth's history with the highest concentration of atmospheric oxygen. During this time interval, at least 1007 apomorphic amino acids were retained in the common ancestor of the extant true bugs. These molecular apomorphies are located in 553 orthologous genes, which suggests the common ancestor of the extant true bugs may have experienced large-scale evolution at the genome level.     

Exploring tree-building methods and distinct molecular data to recover a known asymmetric phage phylogeny

An experimental phylogeny was constructed using bacteriophage T7 and a propagation protocol, in the presence of the mutagen N-methyl-N′-nitro-N′-nitrosoguanidine, based on Hillis et al. [Hillis, D.M., Bull, J.J., White, M.E., Badgett, M.R., Molineux, I.J., 1992. Experimental phylogenetics, generation of a known phylogeny. Science 255, 589–592]. The topology presented in this study has a considerable variation in branch lengths and is less symmetric than the one presented by Hillis et al. [Hillis, D.M., Bull, J.J., White, M.E., Badgett, M.R., Molineux, I.J., 1992. Experimental phylogenetics, generation of a known phylogeny. Science 255, 589–592]. These features are known to present additional difficulties to phylogenetic inference methods. The performance of several phylogenetic methods (conventional and less conventional) was tested using restriction site and nucleotide data. Only methods that encompassed a molecular clock or those based on sequence signatures recovered the true phylogeny. Nevertheless a likelihood ratio test rejected the hypothesis of the existence of a molecular clock when the whole sequence data set was considered. This fact or the particular substitution pattern (mainly G → A and C → T) may be related to the unexpected performance of distance methods based on sequence signatures. To test if the results could have been predicted by simulation studies we estimated the evolution parameters from the real phylogeny and used them to simulate evolution along the same tree (parametric bootstrap). We found that simulation could predict most but not all of the problems encountered by phylogenetic inference methods in the real phylogeny. Short interior branches may be more prone to error than predicted by theoretical studies.     

Brassicaceae phylogeny and trichome evolution

To estimate the evolutionary history of the mustard family (Brassicaceae or Cruciferae), we sampled 113 species, representing 101 of the roughly 350 genera and 17 of the 19 tribes of the family, for the chloroplast gene ndhF. The included accessions increase the number of genera sampled over previous phylogenetic studies by four-fold. Using parsimony, likelihood, and Bayesian methods, we reconstructed the phylogeny of the gene and used the Shimodaira-Hasegawa test (S-H test) to compare the phylogenetic results with the most recent tribal classification for the family. The resultant phylogeny allowed a critical assessment of variations in fruit morphology and seed anatomy, upon which the current classification is based. We also used the S-H test to examine the utility of trichome branching patterns for describing monophyletic groups in the ndhF phylogeny. Our phylogenetic results indicate that 97 of 114 ingroup accessions fall into one of 21 strongly supported clades. Some of these clades can themselves be grouped into strongly to moderately supported monophyletic groups. One of these lineages is a novel grouping overlooked in previous phylogenetic studies. Results comparing 30 different scenarios of evolution by the S-H test indicate that five of 12 tribes represented by two or more genera in the study are clearly polyphyletic, although a few tribes are not sampled well enough to establish para- or polyphyly. In addition, branched trichomes likely evolved independently several times in the Brassicaceae, although malpighiaceous and stellate trichomes may each have a single origin.     

FACTORS DETERMINING THE ACCURACY OF CLADOGRAM ESTIMATION: EVALUATION USING COMPUTER SIMULATION

We developed a simulation model of phylogenesis with which we generated a large number of phylogenies and associated data matrices. We examined the characteristics of these and evaluated the success of three taxonomic methods (Wagner parsimony, character compatibility, and UPGMA clustering) as estimators of phylogeny, paying particular attention to the consequences of changes in certain evolutionary assumptions: relative rate of evolution in three different evolutionary contexts (phyletic, parent lineage, and daughter lineage); relative rate of evolution in different directions (novel forward, convergent forward, or reverse); variation of evolutionary rates; and topology of the phylogenetic tree. Except for variation of evolutionary rates, all the evolutionary parameters that we controlled had significant effects on accuracy of phylogenetic reconstructions. Unexpectedly, the topology of the phylogeny was the most important single factor affecting accuracy; some phylogenies are more readily estimated than others for simply historical reasons. We conclude that none of the three estimation methods is very accurate, that the differences in accuracy among them are rather small, and that historical effects (the branching pattern of a phylogeny) may outweigh biological effects in determining the accuracy with which a phylogeny can be reconstructed.     

Identifying heterogeneity in rates of morphological evolution: discrete character change in the evolution of lungfish (Sarcopterygii; Dipnoi)

Quantifying rates of morphological evolution is important in many macroevolutionary studies, and critical when assessing possible adaptive radiations and episodes of punctuated equilibrium in the fossil record. However, studies of morphological rates of change have lagged behind those on taxonomic diversification, and most authors have focused on continuous characters and quantifying patterns of morphological rates over time. Here, we provide a phylogenetic approach, using discrete characters and three statistical tests to determine points on a cladogram (branches or entire clades) that are characterized by significantly high or low rates of change. These methods include a randomization approach that identifies branches with significantly high rates and likelihood ratio tests that pinpoint either branches or clades that have significantly higher or lower rates than the pooled rate of the remainder of the tree. As a test case for these methods, we analyze a discrete character dataset of lungfish, which have long been regarded as "living fossils" due to an apparent slowdown in rates since the Devonian. We find that morphological rates are highly heterogeneous across the phylogeny and recover a general pattern of decreasing rates along the phylogenetic backbone toward living taxa, from the Devonian until the present. Compared with previous work, we are able to report a more nuanced picture of lungfish evolution using these new methods.     

Measures of clade confidence do not correlate with accuracy of phylogenetic trees

Metrics of phylogenetic tree reliability, such as parametric bootstrap percentages or Bayesian posterior probabilities, represent internal measures of the topological reproducibility of a phylogenetic tree, while the recently introduced aLRT (approximate likelihood ratio test) assesses the likelihood that a branch exists on a maximum-likelihood tree. Although those values are often equated with phylogenetic tree accuracy, they do not necessarily estimate how well a reconstructed phylogeny represents cladistic relationships that actually exist in nature. The authors have therefore attempted to quantify how well bootstrap percentages, posterior probabilities, and aLRT measures reflect the probability that a deduced phylogenetic clade is present in a known phylogeny. The authors simulated the evolution of bacterial genes of varying lengths under biologically realistic conditions, and reconstructed those known phylogenies using both maximum likelihood and Bayesian methods. Then, they measured how frequently clades in the reconstructed trees exhibiting particular bootstrap percentages, aLRT values, or posterior probabilities were found in the true trees. The authors have observed that none of these values correlate with the probability that a given clade is present in the known phylogeny. The major conclusion is that none of the measures provide any information about the likelihood that an individual clade actually exists. It is also found that the mean of all clade support values on a tree closely reflects the average proportion of all clades that have been assigned correctly, and is thus a good representation of the overall accuracy of a phylogenetic tree.     

Misinformative characters and phylogeny shape

The discrepancy between theoretical and observed distributions of tree shapes in recent surveys of phylogeny estimates has lead to investigations of possible biological and methodological causes. I investigated whether the phylogenetic quality of characters is related to the tree shape on which they evolve. Simulated evolution revealed shape-related tendencies for characters to indicate correct cladistic relationships; these differences were measured by examining the characters directly, without deriving any phylogeny estimates. Tree stemminess indices correlated strongly with character quality when characters evolved either speciationally or phyletically. Tree balance was a significant correlate of character quality under speciational evolution but not under phyletic evolution. These results help explain the findings of other simulation studies. With additional study of the behavior of evolving characters and their interaction with phylogenetic methods, we might be able to increase the accuracy of tree estimation and compensate for potential biases related to shape. These results give further reason for caution in trusting phylogeny estimates.     

Molecular phylogeny and evolution of phenotype in silica‐scaled chrysophyte genus Mallomonas

The evolution of phenotypes is highly understudied in protists, due to the dearth of morphological characters, missing fossil record, and/or unresolved phylogeny in the majority of taxa. The chrysophyte genus Mallomonas (Stramenopiles) forms species‐specific silica scales with extraordinary diversity of their ornamentation. In this paper, we molecularly characterized three additional species to provide an updated phylogeny of 43 species, and combined this with evaluations of 24 morphological traits. Using phylogenetic comparative methods, we evaluated phylogenetic signal in traits, reconstructed the trait evolution, and compared the overall phylogenetic and morphological diversity. The majority of traits showed strong phylogenetic signal and mostly dynamic evolution. Phylogenetic relatedness was often reflected by the phenotypic similarity. Both V‐rib and dome are very conservative structures that are presumably involved in precise scale overlap and bristle attachment, respectively. Based on modern species, it seems the dome firstly appeared on apical and/or caudal scales, and only later emerged on body scales. Bristle was presumably present in the common ancestor and gradually elongated ever since. However, most other morphological traits readily changed during the evolution of Mallomonas.     

Phylogenetic Analysis of SARS-CoV-2 Data Is Difficult

Numerous studies covering some aspects of SARS-CoV-2 data analyses are being published on a daily basis, including a regularly updated phylogeny on nextstrain.org. Here, we review the difficulties of inferring reliable phylogenies by example of a data snapshot comprising a quality-filtered subset of 8,736 out of all 16,453 virus sequences available on May 5, 2020 from gisaid.org. We find that it is difficult to infer a reliable phylogeny on these data due to the large number of sequences in conjunction with the low number of mutations. We further find that rooting the inferred phylogeny with some degree of confidence either via the bat and pangolin outgroups or by applying novel computational methods on the ingroup phylogeny does not appear to be credible. Finally, an automatic classification of the current sequences into subclasses using the mPTP tool for molecular species delimitation is also, as might be expected, not possible, as the sequences are too closely related. We conclude that, although the application of phylogenetic methods to disentangle the evolution and spread of COVID-19 provides some insight, results of phylogenetic analyses, in particular those conducted under the default settings of current phylogenetic inference tools, as well as downstream analyses on the inferred phylogenies, should be considered and interpreted with extreme caution.     

Maximum Parsimony and the Skewness Test: A Simulation Study of the Limits of Applicability

The maximum parsimony (MP) method for inferring phylogenies is widely used, but little is known about its limitations in non-asymptotic situations. This study employs large-scale computations with simulated phylogenetic data to estimate the probability that MP succeeds in finding the true phylogeny for up to twelve taxa and 256 characters. The set of candidate phylogenies are taken to be unrooted binary trees; for each simulated data set, the tree lengths of all (2n − 5)!! candidates are computed to evaluate quantities related to the performance of MP, such as the probability of finding the true phylogeny, the probability that the tree with the shortest length is unique, the probability that the true phylogeny has the shortest tree length, and the expected inverse of the number of trees sharing the shortest length. The tree length distributions are also used to evaluate and extend the skewness test of Hillis for distinguishing between random and phylogenetic data. The results indicate, for example, that the critical point after which MP achieves a success probability of at least 0.9 is roughly around 128 characters. The skewness test is found to perform well on simulated data and the study extends its scope to up to twelve taxa.     

Character evolution in Hydrozoa (phylum Cnidaria)

The diversity of hydrozoan life cycles, as manifested in the wide range of polyp, colony, and medusa morphologies, has been appreciated for centuries. Unraveling the complex history of characters involved in this diversity is critical for understanding the processes driving hydrozoan evolution. In this study, we use a phylogenetic approach to investigate the evolution of morphological characters in Hydrozoa. A molecular phylogeny is reconstructed using ribosomal DNA sequence data. Several characters involving polyp, colony, and medusa morphology are coded in the terminal taxa. These characters are mapped onto the phylogeny and then the ancestral character states are reconstructed. This study confirms the complex evolutionary history of hydrozoan morphological characters. Many of the characters involving polyp, colony, and medusa morphology appear as synapomorphies for major hydrozoan clades, yet homoplasy is commonplace.     

The Effect of Recombination on the Accuracy of Phylogeny Estimation

Phylogenetic studies based on DNA sequences typically ignore the potential occurrence of recombination, which may produce different alignment regions with different evolutionary histories. Traditional phylogenetic methods assume that a single history underlies the data. If recombination is present, can we expect the inferred phylogeny to represent any of the underlying evolutionary histories? We examined this question by applying traditional phylogenetic reconstruction methods to simulated recombinant sequence alignments. The effect of recombination on phylogeny estimation depended on the relatedness of the sequences involved in the recombinational event and on the extent of the different regions with different phylogenetic histories. Given the topologies examined here, when the recombinational event was ancient, or when recombination occurred between closely related taxa, one of the two phylogenies underlying the data was generally inferred. In this scenario, the evolutionary history corresponding to the majority of the positions in the alignment was generally recovered. Very different results were obtained when recombination occurred recently among divergent taxa. In this case, when the recombinational breakpoint divided the alignment in two regions of similar length, a phylogeny that was different from any of the true phylogenies underlying the data was inferred.     

Phylogenetically structured variance in felid bite force: the role of phylogeny in the evolution of biting performance

A key question in evolution is the degree to which morphofunctional complexes are constrained by phylogeny. We investigated the role of phylogeny in the evolution of biting performance, quantified as bite forces, using phylogenetic eigenvector regression. Results indicate that there are strong phylogenetic signals in both absolute and size‐adjusted bite forces, although it is weaker in the latter. This indicates that elimination of size influences reduces the level of phylogenetic inertia and that the majority of the phylogenetic constraint is a result of size. Tracing the evolution of bite force through phylogeny by character optimization also supports this notion, in that relative bite force is randomly distributed across phylogeny whereas absolute bite force diverges according to clade. The nonphylogenetically structured variance in bite force could not be sufficiently explained by species‐unique morphology or by ecology. This study demonstrates the difficulties in identifying causes of nonphylogenetically structured variance in morphofunctional character complexes.     

How our view of animal phylogeny was reshaped by molecular approaches: lessons learned

In the late 1980s, researchers began applying molecular sequencing tools to questions of deep animal phylogeny. These advances in sequencing were accompanied with improvements in computation and phylogenetic methods, and served to significantly reshape our understanding of metazoan evolution. Prior to this time, researchers asserted phylogenetic hypotheses based on their experience with taxa and to some degree, their authority. Molecular phylogenetic tools provided discrete methods and objective characters for reconstructing phylogeny. Nonetheless, major changes to widely accepted views, such as animal phylogeny, take time to be accepted. Development and acceptance of our current understanding of animal evolution occurred in three main phases: initial hypotheses based on 18S data, confirmation with additional molecular markers, and continued refinement with phylogenomics. With the advent of ideas such as Lophotrochozoa and Ecdysozoa, flaws in the traditional view became apparent. We now understand that complex morphological and embryological features (e.g., segmentation, coelom formation, development of body cavities) are much more evolutionarily plastic than previously recognized. Here, I explore how the transition from the traditional to the modern phylogenetic understanding of animal phylogeny occurred and examine some implications of this change in understanding. As the field moves forward, the utility of morphological and embryological characters for reconstruction of deep animal phylogeny should be discouraged. Instead, these characters should be interpreted in the light of independent phylogeny.     

Does a tree-like phylogeny only exist at the tips in the prokaryotes?

The extent to which prokaryotic evolution has been influenced by horizontal gene transfer (HGT) and therefore might be more of a network than a tree is unclear. Here we use supertree methods to ask whether a definitive prokaryotic phylogenetic tree exists and whether it can be confidently inferred using orthologous genes. We analysed an 11-taxon dataset spanning the deepest divisions of prokaryotic relationships, a 10-taxon dataset spanning the relatively recent gamma-proteobacteria and a 61-taxon dataset spanning both, using species for which complete genomes are available. Congruence among gene trees spanning deep relationships is not better than random. By contrast, a strong, almost perfect phylogenetic signal exists in gamma-proteobacterial genes. Deep-level prokaryotic relationships are difficult to infer because of signal erosion, systematic bias, hidden paralogy and/or HGT. Our results do not preclude levels of HGT that would be inconsistent with the notion of a prokaryotic phylogeny. This approach will help decide the extent to which we can say that there is a prokaryotic phylogeny and where in the phylogeny a cohesive genomic signal exists.     

CONDUCTING PHYLOGENETIC COMPARATIVE STUDIES WHEN THE PHYLOGENY IS NOT KNOWN

A method is proposed to conduct phylogenetic analyses of comparative or interspecific data when the true phylogeny is not known. Standard models of speciation and/or extinction or other methods are used to generate a sample from the set of all possible phylogenies for the measured species. The comparative data are then analyzed on each of the possible trees to obtain a distribution of possible evolutionary statistics for these data. The mean of this distribution is proposed as a reasonable estimate of the true evolutionary statistic of interest. Ways of obtaining confidence intervals and of developing hypothesis tests for this mean statistic are also proposed. The method can be used with any comparative method or phylogenetic analysis technique when phylogenetic relationships among species are not known or when branch lengths for a phylogeny in units of expected character change (as required by most methods) are not available. Computer programs to conduct the analyses are available on request.     

Phylogenomics of Fargesia and Yushania reveals a history of reticulate evolution

Reticulate evolution is a common and important driving force in angiosperm evolution. In this study, we analyzed the phylogenetic signals of genomic regions with different inheritance patterns to understand the evolutionary process of organisms using species-rich Himalaya–Hengduan taxa of bamboos (Fargesia Franchet and Yushania Keng). We constructed phylogenetic trees using different sampling strategies and reconstruction methods based on genome skimming and double digest restriction-site-associated DNA sequencing data. We assessed the congruence of topologies generated from different datasets and employed several approaches to reveal the causes of phylogenetic incongruence, including the detection of hybridization and introgression using PhyloNetworks and the D-statistic test (ABBA-BABA test). We found that, in the plastome-based phylogeny, Fargesia bamboos can be clustered into three groups and Yushania was nested within one of them, which contradicts the nuclear–double digest restriction-site-associated DNA sequencing-based phylogeny. Moreover, the genetic variation of chloroplast DNA is significantly correlated with geographical distribution. The strong signal of incomplete lineage sorting, hybridization, introgression, and cytoplasmic gene flow found among genera and species suggests that reticulate evolution is the main cause for the phylogenetic incongruence between nuclear and chloroplast datasets. Our results add evidence that genomes with different inheritance patterns can reveal distinct evolutionary histories of species and suggest that reticulate evolution is prevalent in rapidly diversifying groups.     

The quest for orthologs: finding the corresponding gene across genomes

Orthology is a key evolutionary concept in many areas of genomic research. It provides a framework for subjects as diverse as the evolution of genomes, gene functions, cellular networks and functional genome annotation. Although orthologous proteins usually perform equivalent functions in different species, establishing true orthologous relationships requires a phylogenetic approach, which combines both trees and graphs (networks) using reliable species phylogeny and available genomic data from more than two species, and an insight into the processes of molecular evolution. Here, we evaluate the available bioinformatics tools and provide a set of guidelines to aid researchers in choosing the most appropriate tool for any situation.