Molecular Phylogenetics and Evolution最新文献

Poales, as one of the largest orders of angiosperm, holds crucial economic and ecological importance. Nevertheless, achieving a consensus topology has been challenging in previous studies due to limited molecular data and sparse taxon sampling. The uneven distribution of species diversity among families and the factors leading to elevated species richness in certain lineages have also been subjects of ongoing discussion and investigation. In this study, we conducted a comprehensive sampling, including representatives from all 14 families and 85 taxa of Poales, along with five additional outgroups. To reconstruct the phylogeny of Poales, we employed a combination of coalescent and concatenation methods on three nuclear gene sets (1093, 491, 143) and one plastid gene set (53), which were inferenced from genomic data. We also conducted phylogenetic hypothesis analyses to evaluate two major conflicting nodes detected in phylogenetic analyses. As a result, we successfully resolved the backbone of Poales and provided a timeline for its evolutionary history. We recovered the sister relationship between Typhaceae and Bromeliaceae as the earliest diverging families within Poales. The clade consisting of Ecdeiocoleaceae and Joinvilleaceae was recovered as the sister group of Poaceae. Within the xyrid clade, Mayacaceae and Erioaculaceae + Xyridaceae successively diverged along the backbone of Poales. The topology of [Aristidoideae, ((Micrairoideae, Panicoideae), (Arundinoideae, (Chloridoideae, Danthonioideae)))] within the PACMAD clade has received strong support from multiple findings. We also delved into the underlying biological factors that contributed to the conflicting nodes observed in the phylogenetic analysis. Apart from the uncertainty regarding the sister group of Poaceae caused by cytonuclear discordance, frequent hybridization and polyploidy may have contributed to other conflicting nodes. We identified 26 putative whole-genome duplication (WGD) events within Poales. However, apart from the σ-WGD and the ρ-WGD, we did not observe any potential polyploid events that could be directly linked to the species diversification in specific lineages. Furthermore, there was a significant increase in the net diversification rate of Poales following the K-Pg boundary.

Phylogenetic tree reconstruction with molecular data is important in many fields of life science research. The gold standard in this discipline is the phylogenetic tree reconstruction based on the Maximum Likelihood method. In this study, we present neural networks to predict the best model of sequence evolution and the correct topology for four sequence alignments of nucleotide or amino acid sequence data. We trained neural networks with different architectures using simulated alignments for a wide range of evolutionary models, model parameters and branch lengths. By comparing the accuracy of model and topology prediction of the trained neural networks with Maximum Likelihood and Neighbour Joining methods, we show that for quartet trees, the neural network classifier outperforms the Neighbour Joining method and is in most cases as good as the Maximum Likelihood method to infer the best model of sequence evolution and the best tree topology. These results are consistent for nucleotide and amino acid sequence data. We also show that our method is superior for model selection than previously published methods based on convolutionary networks. Furthermore, we found that neural network classifiers are much faster than the IQ-TREE implementation of the Maximum Likelihood method. Our results show that neural networks could become a true competitor for the Maximum Likelihood method in phylogenetic reconstructions.

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

Silkmoths (Bombycidae) have a disjunct distribution predominantly in the Southern Hemisphere and Asia. Here we reconstruct the phylogenetic history of the family to test competing hypotheses on their origin and assess how vicariance and long-distance dispersal shaped their current distribution. We sequenced up to 5,074 base pairs from six loci (COI, EF1-α, wgl, CAD, GAPDH, and RpS5) to infer the historical biogeography of Bombycidae. The multilocus dataset covering 20 genera (80 %) of the family, including 17 genera (94 %) of Bombycinae and 3 genera (43 %) of Epiinae, was used to estimate phylogenetic patterns, divergence times and biogeographic reconstruction. Dating estimates extrapolated from secondary calibration sources indicate the Bombycidae stem-group originated approximately 64 Mya. The subfamilies Epiinae (South America) and Bombycinae (Australia, Asia, East Palaearctic, and Africa) were reciprocally monophyletic, diverging at c. 56 Mya (95 % credibility interval: 66–46 Mya). The ‘basal’ lineage of Bombycinae — Gastridiota + Elachyophtalma — split from the rest of Bombycinae c. 53 Mya (95 % credibility interval: 63–43 Mya). Gastridiota is a monobasic genus with a relictual distribution in subtropical forests of eastern Australia. The Oriental and African genera comprised a monophyletic group: the Oriental region was inferred to have been colonized from a long-distance dispersal event from Australia to South-East Asia c. 53 Mya or possibly later (c. 36–26 Mya); Africa was subsequently colonized by dispersal from Asia c. 16 Mya (95 % credibility interval: 21–12 Mya). Based on the strongly supported phylogenetic relationships and estimates of divergence times, we conclude that Bombycidae had its origin in the fragment of Southern Gondwana consisting of Australia, Antarctica and South America during the Paleocene. The disjunction between South America (Epiinae) and Australia (Bombycinae) is best explained by vicariance in the Eocene, whereas the disjunct distribution in Asia and Africa is best explained by more recent dispersal events.

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

Springtails (Collembola) stand as one of the most abundant, widespread, and ancient terrestrial arthropods on earth. However, their evolutionary history and deep phylogenetic relationships remain elusive. In this study, we employed phylogenomic approaches to elucidate the basal relationships among Collembola. We sampled whole-genome data representing all major collembolan lineages in proportion to their known diversity. To account for potential phylogenomic biases, we implemented various data extraction, locus sampling, and signal filtering strategies to generate matrices. Subsequently, we applied a diverse array of tree-searching and rate-modelling methods to reconstruct the phylogeny. Our analyses, utilizing different matrices and methods, converged on the same unrooted relationships among collembolan ingroups, supporting the current ordinal classification and challenging the monophyly of Arthropleona and Symphypleona s.l. However, discrepancies across analyses existed in the root of Collembola. Among various root positions, those based on more informative matrices and biologically realistic models, favoring a basal topology of Entomobryomorpha + (Symphypleona s.s. + (Neelipleona + Poduromorpha)), were supported by subsequent methodological assessment, topology tests, and rooting analyses. This optimal topology suggests multiple independent reduction of the pronotum in non-poduromorph orders and aligns with the plesiomorphic status of neuroendocrine organs and epicuticular structure of Entomobryomorpha. Fossil-calibrated dating analyses based on the optimal topology indicated late-Paleozoic to mid-Mesozoic origins of the crown Collembola and four orders. In addition, our results questioned the monophyly of Isotomidae and Neanuridae, underscoring the need for further attention to the systematics of these families. Overall, this study provides novel insights into the phylogenetic backbone of Collembola, which will inform future studies on the systematics, ecology, and evolution of this significant arthropod lineage.

Phylogenetic trees are believed to contain a wealth of information on diversification processes. However, comparing phylogenetic trees is not straightforward due to their high dimensionality. Researchers have therefore defined a wide range of low-dimensional summary statistics. Currently, it remains unexplored to what extent these summary statistics cover the same underlying information and what summary statistics best explain observed variation across phylogenies. Furthermore, a large subset of available summary statistics focusses on measuring the topological features of a phylogenetic tree, but are often only explored at the extreme edge cases of the fully balanced or imbalanced tree and not for trees of intermediate balance.

Here, we introduce a new R package called ‘treestats’, that provides speed optimized code to compute 70 summary statistics. We study correlations between summary statistics on empirical trees and on trees simulated using several diversification models. Furthermore, we introduce an algorithm to create intermediately balanced trees in a well-defined manner, in order to explore variation in summary statistics across a balance gradient.

We find that almost all summary statistics are correlated with tree size, and find that it is difficult, if not impossible, to correct for tree size, unless the tree generating model is known. Furthermore, we find that across empirical and simulated trees, at least three large clusters of correlated summary statistics can be found, where statistics group together based on information used (topology or branching times). However, the finer grained correlation structure appears to depend strongly on either the taxonomic group studied (in empirical studies) or the tree generating model (in simulation studies).

Amongst statistics describing the (im)balance of a tree, we find that almost all statistics vary non-linearly, and sometimes even non-monotonically, with our generated balance gradient. This indicates that balance is perhaps a more complex property of a tree than previously thought. Furthermore, using our new imbalancing algorithm, we devise a numerical test to identify balance statistics, and identify several statistics as balance statistics that were not previously considered as such. Lastly, our results lead to several recommendations on which statistics to select when analyzing and comparing phylogenetic trees.

Green algae usually assigned to the genus Oophila are known to colonize egg capsules of amphibian egg masses across the Nearctic and Palearctic regions. We study the phylogenetic relationships of these algae using a phylotranscriptomic data set of 76 protein-coding single-copy nuclear genes. Our data set includes novel RNAseq data for six amphibian-associated and five free-living green algae, and draft genomes of two of the latter. Within the Oophila clade (nested within Moewusinia), we find samples from two European frogs (Rana dalmatina and R. temporaria) closely related to those of the North American frog R. aurora (Oophila subclade III). An isolate from the North American R. sylvatica (subclade IV) appears to be sister to the Japanese isolate from the salamander Hynobius nigrescens (subclade J1), and subclade I algae from Ambystoma maculatum are sister to all other lineages in the Oophila clade. Two free-living algae (Chlamydomonas nasuta and Cd. pseudogloeogama) are nested within the Oophila clade, and a strain of the type species of Chlorococcum (Cc. infusionum) is related to this assemblage. Our phylotranscriptomic tree suggests that recognition of different species within the Oophila clade (“clade B” of earlier studies) is warranted, and calls for a comprehensive taxonomic revision of Moewusinia.

Disentangling the genomic intricacies underlying speciation and the causes of discordance between sources of evidence can offer remarkable insights into evolutionary dynamics. The ant-eating spider Zodarion nitidum, found across the Middle East and Egypt, displays yellowish and blackish morphs that co-occur sympatrically. These morphs additionally differ in behavioral and physiological features and show complete pre-mating reproductive isolation. In contrast, they possess similar sexual features and lack distinct differences in their mitochondrial DNA. We analyzed both Z. nitidum morphs and outgroups using genome-wide and additional mitochondrial DNA data. The genomic evidence indicated that Yellow and Black are reciprocally independent lineages without signs of recent admixture. Interestingly, the sister group of Yellow is not Black but Z. luctuosum, a morphologically distinct species. Genomic gene flow analyses pinpointed an asymmetric nuclear introgression event, with Yellow contributing nearly 5 % of its genome to Black roughly 320,000 years ago, intriguingly aligning with the independently estimated origin of the mitochondrial DNA of Black. We conclude that the blackish and yellowish morphs of Z. nitidum are long-diverged distinct species, and that the ancient and modest genomic introgression event registered resulted in a complete mitochondrial takeover of Black by Yellow. This investigation underscores the profound long-term effects that even modest hybridization events can have on the genome of organisms. It also exemplifies the utility of phylogenetic networks for estimating historical events and how integrating independent lines of evidence can increase the reliability of such estimations.