Taming the BEAST-A Community Teaching Material Resource for BEAST 2

A clear distinction and presence of Acropora aff. divaricata within Acropora cf. solitaryensis species complex along their biogeographic distribution in East Asia

In the Anthropocene, scleractinian corals face unprecedented threats from synergistic stressors, including rising seawater temperatures that surpass critical thresholds that lead to global coral reef degradation. With over 1,698 coral species in the order Scleractinia, their conservation is increasingly complex due to their morphological plasticity and the challenge of accurate species identification. The genus Acropora, with approximately 400 nominal species, exemplifies these challenges, as morphological traits often vary within and among species, complicating taxonomic efforts. Traditional methods based on skeletal characteristics are insufficient for delineating Acropora species, prompting the use of integrative approaches combining morphology, reproduction, and molecular data. In this study, we employ multi-locus phylogenetic analyses and morphometric measurements to distinguish different growth forms of Acropora cf. solitaryensis as distinct species and delineate the species range boundaries of A. aff. divaricata and A. cf. solitaryensis in East Asian coral ecosystems. We identify arborescent and intermediate morphotypes belonging to A. aff. divaricata, which is distributed from tropical reefs in southeastern Taiwan to temperate non-reefal communities in Shikoku, Japan. Conversely, the solid-plate morphotype aligns with A. cf. solitaryensis from the holotype locality at Solitary Island, Australia, found primarily in subtropical non-reefal regions in northern Taiwan and Japan. The distinct distribution patterns of A. aff. divaricata and A. cf. solitaryensis underscore the necessity for biogeographic sampling in Acropora taxonomy, considering the Kuroshio Current's impact on coral distributions, and a re-evaluation of poleward coral migration or expansion due to climate change. Our findings challenge the traditional taxonomy of A. aff. divaricata and A. cf. solitaryensis, revealing that they may instead encompass multiple species. This has significant implications for coral conservation strategies, as accurate species identification is crucial for understanding coral responses to environmental changes and informing conservation efforts in the face of climate change.

PhyloJunction: A Computational Framework for Simulating, Developing, and Teaching Evolutionary Models

We introduce PhyloJunction, a computational framework designed to facilitate the prototyping, testing, and characterization of evolutionary models. PhyloJunction is distributed as an open-source Python library that can be used to implement a variety of models, thanks to its flexible graphical modeling architecture and dedicated model specification language. Model design and use are exposed to users via command-line and graphical interfaces, which integrate the steps of simulating, summarizing, and visualizing data. This article describes the features of PhyloJunction-which include, but are not limited to, a general implementation of a popular family of phylogenetic diversification models-and, moving forward, how it may be expanded to not only include new models, but to also become a platform for conducting and teaching statistical learning.

Museum specimens shedding light on the evolutionary history and cryptic diversity of the hedgehog family Erinaceidae

The family Erinaceidae encompasses 27 extant species in two subfamilies: Erinaceinae, which includes spiny hedgehogs, and Galericinae, which comprises silky-furred gymnures and moonrats. Although they are commonly recognized by the general public, their phylogenetic history remains incompletely understood, and several species have never been included in any molecular analyses. Additionally, previous research suggested that the species diversity of Erinaceidae might be underestimated. In this study, we sequenced the mitochondrial genomes of 29 individuals representing 18 erinaceid species using 18 freshly collected tissue and 11 historical museum specimens. We also integrated previously published data for a concatenated analysis. We aimed to elucidate the evolutionary relationships within Erinaceidae, estimate divergence times, and uncover potential underestimated species diversity. Our data finely resolved intergeneric and interspecific relationships and presented the first molecular evidence for the phylogenetic position of Mesechinus wangi, Paraechinus micropus, and P. nudiventris. Our results revealed a sister relationship between Neotetracus and Neohylomys gymnures, as well as a sister relationship between Hemiechinus and Mesechinus, supporting previous hypotheses. Additionally, our findings provided a novel phylogenetic position for Paraechinus aethiopicus, placing it in a basal position within the genus. Furthermore, our study uncovered cryptic species diversity within Hylomys suillus as well as in Neotetracus sinensis, Atelerix albiventris, P. aethiopicus, and Hemiechinus auratus, most of which have been previously overlooked.

Mechanistic phylodynamic models do not provide conclusive evidence that non-avian dinosaurs were in decline before their final extinction

Phylodynamic models can be used to estimate diversification trajectories from time-calibrated phylogenies. Here we apply two such models to phylogenies of non-avian dinosaurs, a clade whose evolutionary history has been widely debated. Although some authors have suggested that the clade experienced a decline in diversity, potentially starting millions of years before the end-Cretaceous mass extinction, others have suggested that the group remained highly diverse right up until the Cretaceous-Paleogene (K-Pg) boundary. Our results show that model assumptions, likely with respect to incomplete sampling, have a large impact on whether dinosaurs appear to have experienced a long-term decline or not. The results are also highly sensitive to the topology and branch lengths of the phylogeny used. Developing comprehensive models of sampling bias, and building larger and more accurate phylogenies, are likely to be necessary steps for us to determine whether dinosaur diversity was or was not in decline before the end-Cretaceous mass extinction.

Young evolutionary origins of dioecy in the genus Asparagus

PREMISE

Dioecy (separate sexes) has independently evolved numerous times across the angiosperm phylogeny and is recently derived in many lineages. However, our understanding is limited regarding the evolutionary mechanisms that drive the origins of dioecy in plants. The recent and repeated evolution of dioecy across angiosperms offers an opportunity to make strong inferences about the ecological, developmental, and molecular factors influencing the evolution of dioecy, and thus sex chromosomes. The genus Asparagus (Asparagaceae) is an emerging model taxon for studying dioecy and sex chromosome evolution, yet estimates for the age and origin of dioecy in the genus are lacking.

METHODS

We use plastome sequences and fossil time calibrations in phylogenetic analyses to investigate the age and origin of dioecy in the genus Asparagus. We also review the diversity of sexual systems present across the genus to address contradicting reports in the literature.

RESULTS

We estimate that dioecy evolved once or twice approximately 2.78-3.78 million years ago in Asparagus, of which roughly 27% of the species are dioecious and the remaining are hermaphroditic with monoclinous flowers.

CONCLUSIONS

Our findings support previous work implicating a young age and the possibility of two origins of dioecy in Asparagus, which appear to be associated with rapid radiations and range expansion out of Africa. Lastly, we speculate that paleoclimatic oscillations throughout northern Africa may have helped set the stage for the origin(s) of dioecy in Asparagus approximately 2.78-3.78 million years ago.

Equine herpesvirus 4 infected domestic horses associated with Sintashta spoke-wheeled chariots around 4,000 years ago

Equine viral outbreaks have disrupted the socio-economic life of past human societies up until the late 19th century and continue to be of major concern to the horse industry today. With a seroprevalence of 60-80 per cent, equine herpesvirus 4 (EHV-4) is the most common horse pathogen on the planet. Yet, its evolutionary history remains understudied. Here, we screen the sequenced data of 264 archaeological horse remains to detect the presence of EHV-4. We recover the first ancient EHV-4 genome with 4.2× average depth-of-coverage from a specimen excavated in the Southeastern Urals and dated to the Early Bronze Age period, approximately 3,900 years ago. The recovery of an EHV-4 virus outside the upper respiratory tract not only points to an animal particularly infected but also highlights the importance of post-cranial bones in pathogen characterisation. Bayesian phylogenetic reconstruction provides a minimal time estimate for EHV-4 diversification to around 4,000 years ago, a time when modern domestic horses spread across the Central Asian steppes together with spoke-wheeled Sintashta chariots, or earlier. The analyses also considerably revise the diversification time of the two EHV-4 subclades from the 16th century based solely on modern data to nearly a thousand years ago. Our study paves the way for a robust reconstruction of the history of non-human pathogens and their impact on animal health.

PhyloJunction: a computational framework for simulating, developing, and teaching evolutionary models

We introduce PhyloJunction, a computational framework designed to facilitate the prototyping, testing, and characterization of evolutionary models. PhyloJunction is distributed as an open-source Python library that can be used to implement a variety of models, through its flexible graphical modeling architecture and dedicated model specification language. Model design and use are exposed to users via command-line and graphical interfaces, which integrate the steps of simulating, summarizing, and visualizing data. This paper describes the features of PhyloJunction - which include, but are not limited to, a general implementation of a popular family of phylogenetic diversification models - and, moving forward, how it may be expanded to not only include new models, but to also become a platform for conducting and teaching statistical learning.

Reconstructing the phylodynamic history and geographic spread of the CRF01_AE-predominant HIV-1 epidemic in the Philippines from PR/RT sequences sampled from 2008 to 2018

The Philippines has had a rapidly growing human immunodeficiency virus (HIV) epidemic with a shift in the prevalent subtype from B to CRF01_AE. However, the phylodynamic history of CRF01_AE in the Philippines has yet to be reconstructed. We conducted a descriptive retrospective study reconstructing the history of HIV-1 CRF01_AE transmissions in the Philippines through molecular epidemiology. Partial polymerase sequences (n = 1144) collected between 2008 and 2018 from three island groups were collated from the Research Institute for Tropical Medicine drug resistance genotyping database. Estimation of the time to the most recent common ancestor (tMRCA), effective reproductive number (Re), effective viral population size (Ne), relative migration rates, and geographic spread of CRF01_AE was performed with BEAST. Re and Ne were compared between CRF01_AE and B. Most CRF01_AE sequences formed a single clade with a tMRCA of June 1996 [95 per cent highest posterior density (HPD): December 1991, October 1999]. An increasing CRF01_AE Ne was observed from the tMRCA to 2013. The CRF01_AE Re reached peaks of 2.46 [95 per cent HPD: 1.76, 3.27] in 2007 and 2.52 [95 per cent HPD: 1.83, 3.34] in 2015. A decrease of CRF01_AE Re occurred in the intervening years of 2007 to 2011, reaching as low as 1.43 [95 per cent HPD: 1.06, 1.90] in 2011, followed by a rebound. The CRF01_AE epidemic most likely started in Luzon and then spread to the other island groups of the country. Both CRF01_AE and Subtype B exhibited similar patterns of Re fluctuation over time. These results characterize the subtype-specific phylodynamic history of the largest CRF01_AE cluster in the Philippines, which contextualizes and may inform past, present, and future public health measures toward controlling the HIV epidemic in the Philippines.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

Highly structured populations of deep-sea copepods associated with hydrothermal vents across the Southwest Pacific, despite contrasting life history traits

Hydrothermal vents are extreme environments, where abundant communities of copepods with contrasting life history traits co-exist along hydrothermal gradients. Here, we discuss how these traits may contribute to the observed differences in molecular diversity and population genetic structure. Samples were collected from vent locations across the globe including active ridges and back-arc basins and compared to existing deep-sea hydrothermal vent and shallow water data, covering a total of 22 vents and 3 non-vent sites. A total of 806 sequences of mtDNA from the Cox1 gene were used to reconstruct the phylogeny, haplotypic relationship and demography within vent endemic copepods (Dirivultidae, Stygiopontius spp.) and non-vent-endemic copepods (Ameiridae, Miraciidae and Laophontidae). A species complex within Stygiopontius lauensis was studied across five pacific back-arc basins at eight hydrothermal vent fields, with cryptic species being restricted to the basins they were sampled from. Copepod populations from the Lau, North Fiji and Woodlark basins are undergoing demographic expansion, possibly linked to an increase in hydrothermal activity in the last 10 kya. Highly structured populations of Amphiascus aff. varians 2 were also observed from the Lau to the Woodlark basins with populations also undergoing expansion. Less abundant harpacticoids exhibit little to no population structure and stable populations. This study suggests that similarities in genetic structure and demography may arise in vent-associated copepods despite having different life history traits. As structured meta-populations may be at risk of local extinction should major anthropogenic impacts, such as deep-sea mining, occur, we highlight the importance of incorporating a trait-based approach to investigate patterns of genetic connectivity and demography, particularly regarding area-based management tools and environmental management plans.

Phylogenetic identification of influenza virus candidates for seasonal vaccines

The seasonal influenza (flu) vaccine is designed to protect against those influenza viruses predicted to circulate during the upcoming flu season, but identifying which viruses are likely to circulate is challenging. We use features from phylogenetic trees reconstructed from hemagglutinin (HA) and neuraminidase (NA) sequences, together with a support vector machine, to predict future circulation. We obtain accuracies of 0.75 to 0.89 (AUC 0.83 to 0.91) over 2016-2020. We explore ways to select potential candidates for a seasonal vaccine and find that the machine learning model has a moderate ability to select strains that are close to future populations. However, consensus sequences among the most recent 3 years also do well at this task. We identify similar candidate strains to those proposed by the World Health Organization, suggesting that this approach can help inform vaccine strain selection.

Phylogenetic evidence reveals early Kra-Dai divergence and dispersal in the late Holocene

Studying language evolution brings a crucial perspective to bear on questions of human prehistory. As the most linguistically diverse region on earth, East and Southeast Asia have witnessed extensive sociocultural and ethnic contacts among different language communities. Especially, the Kra-Dai language family exhibits tremendous socio-cultural importance in these regions. Due to limited historical accounts, however, there are several controversies on their linguistic relatedness, ambiguities regarding the divergence time, and uncertainties on the dispersal patterns. To address these issues, here we apply Bayesian phylogenetic methods to analyze the largest lexical dataset containing 646 cognate sets compiled for 100 Kra-Dai languages. Our dated phylogenetic tree showed their initial divergence occurring approximately 4000 years BP. Phylogeographic results supported the early Kra-Dai language dispersal from the Guangxi-Guangdong area of South China towards Mainland Southeast Asia. Coupled with genetic, archaeological, paleoecologic, and paleoclimatic data, we demonstrated that the Kra-Dai language diversification could have coincided with their demic diffusion and agricultural spread shaped by the global climate change in the late Holocene. The interdisciplinary alignments shed light on reconstructing the prehistory of Kra-Dai languages and provide an indispensable piece of the puzzle for further studying prehistoric human activities in East and Southeast Asia.

Presumptive First Record of Myotis aurascens (Chiroptera, Vespertilionidae) from China with a Phylogenetic Analysis

Bat groups have a high degree of species diversity, and the taxonomic status and phylogenetic relationships among bat species have always been research hotspots. Due to the fact that morphological characteristics do not always reflect the evolutionary relationships among species, mitochondrial DNA has been widely used in the study of species relationships due to its maternal inheritance pattern. Myotis aurascens has been suggested as a possible synonym for M. davidii. However, the status of this classification has been controversial. In this study, the morphological and molecular characteristics of a M. aurascens captured from Inner Mongolia, China, were analyzed to determine its taxonomic status. In terms of morphological features, the body weight was 6.33 g, the head and body length were 45.10 mm, the forearm length was 35.87 mm, and the tragus length was 7.51 mm. These values all fell within the species signature data range. Nucleotide skew analysis of the protein-coding genes (PCGs) suggested that only five PCGs (ND1, ND2, COX2, ATP8, and ND4) showed AT-skew value within the mitogenome of M. aurascens. Except for ND6, the GC-skew values of the other PCGs were negative, reflecting the preference for C and T bases compared to G and A bases. Molecular phylogenetic analyses based on mitochondrial PCGs indicated that M. aurascens was a distinct species from M. davidii and phylogenetically closer to M. ikonnikovi, M. alcathoe, and M. mystacinus. Genetic distance analysis also showed that M. aurascens and M. davidii were distantly related. Therefore, the integrated analysis demonstrated that M. aurascens should be considered a distinct species rather than a synonym of M. davidii. Our study could provide a reference for enriching species diversity and research on conservation in China.

Role of paleoclimatic and paleohydrological processes in lineage divergence in freshwater organisms: A snippet from lentic genus Pila

The Indian subcontinent is extremely diverse in terms of its flora and fauna. However, only a handful of studies have aimed to understand the diversity of freshwater invertebrates using multiple lines of evidence in recent times. Here we aimed to estimate the cryptic diversity of two widespread freshwater snail species within the genus Pila (Röding, 1798) and uncover the processes behind lineage diversification in these species. We sequenced mitochondrial and nuclear markers from a comprehensive sampling of specimens from different river basins in India. We implemented an integrative taxonomy approach to delimit the lineages in these groups, employing phylogenetic, geometric morphometric and niche modelling-based methods. Then, we investigated the drivers of lineage divergence in these species using population genetic tools in conjunction with divergence time estimation. We found that both species consist of several genetically and ecologically distinct lineages. The genetic data showed that several of these lineages are restricted to a single or a few river basins. The divergence time estimation analyses indicated that the time frame of divergence within the species coincided with paleohydrological and paleoclimatic events in the Miocene. The diversification was primarily driven by allopatric isolation into different river basins. To conclude, the study sheds light on the complex interaction between the habitat preference of the species and the environment in shaping the diversification patterns in this group.

Propagating uncertainty about molecular evolution models and prior distributions to phylogenetic trees

Phylogenetic trees constructed from molecular sequence data rely on largely arbitrary assumptions about the substitution model, the distribution of substitution rates across sites, the version of the molecular clock, and, in the case of Bayesian inference, the prior distribution. Those assumptions affect results reported in the form of clade probabilities and error bars on divergence times and substitution rates. Overlooking the uncertainty in the assumptions leads to overly confident conclusions in the form of inflated clade probabilities and short confidence intervals or credible intervals. This paper demonstrates how to propagate that uncertainty by combining the models considered along with all of their assumptions, including their prior distributions. The combined models incorporate much more of the uncertainty than Bayesian model averages since the latter tend to settle on a single model due to the higher-level assumption that one of the models is true. Nucleotide sequence data illustrates the proposed model combination method.

Population structure and evolutionary history of the greater cane rat (Thryonomys swinderianus) from the Guinean Forests of West Africa

Grasscutter (Thryonomys swinderianus) is a large-body old world rodent found in sub-Saharan Africa. The body size and the unique taste of the meat of this major crop pest have made it a target of intense hunting and a potential consideration as a micro-livestock. However, there is insufficient knowledge on the genetic diversity of its populations across African Guinean forests. Herein, we investigated the genetic diversity, population structures and evolutionary history of seven Nigerian wild grasscutter populations together with individuals from Cameroon, Republic of Benin, and Ghana, using five mitochondrial fragments, including D-loop and cytochrome b (CYTB). D-loop haplotype diversity ranged from 0.571 (± 0.149) in Republic of Benin to 0.921 (± 0.013) in Ghana. Within Nigeria, the haplotype diversity ranged from 0.659 (± 0.059) in Cross River to 0.837 (± 0.075) in Ondo subpopulation. The fixation index (FST), haplotype frequency distribution and analysis of molecular variance revealed varying levels of population structures across populations. No significant signature of population contraction was detected in the grasscutter populations. Evolutionary analyses of CYTB suggests that South African population might have diverged from other populations about 6.1 (2.6-10.18, 95% CI) MYA. Taken together, this study reveals the population status and evolutionary history of grasscutter populations in the region.

Comparing the transmission potential from sequence and surveillance data of 2009 North American influenza pandemic waves

Technological advancements in phylodynamic modeling coupled with the accessibility of real-time pathogen genetic data are increasingly important for understanding the infectious disease transmission dynamics. In this study, we compare the transmission potentials of North American influenza A(H1N1)pdm09 derived from sequence data to that derived from surveillance data. The impact of the choice of tree-priors, informative epidemiological priors, and evolutionary parameters on the transmission potential estimation is evaluated. North American Influenza A(H1N1)pdm09 hemagglutinin (HA) gene sequences are analyzed using the coalescent and birth-death tree prior models to estimate the basic reproduction number (R ₀ ). Epidemiological priors gathered from published literature are used to simulate the birth-death skyline models. Path-sampling marginal likelihood estimation is conducted to assess model fit. A bibliographic search to gather surveillance-based R ₀ values were consistently lower (mean ≤ 1.2) when estimated by coalescent models than by the birth-death models with informative priors on the duration of infectiousness (mean ≥ 1.3 to ≤2.88 days). The user-defined informative priors for use in the birth-death model shift the directionality of epidemiological and evolutionary parameters compared to non-informative estimates. While there was no certain impact of clock rate and tree height on the R ₀ estimation, an opposite relationship was observed between coalescent and birth-death tree priors. There was no significant difference (p = 0.46) between the birth-death model and surveillance R ₀ estimates. This study concludes that tree-prior methodological differences may have a substantial impact on the transmission potential estimation as well as the evolutionary parameters. The study also reports a consensus between the sequence-based R ₀ estimation and surveillance-based R ₀ estimates. Altogether, these outcomes shed light on the potential role of phylodynamic modeling to augment existing surveillance and epidemiological activities to better assess and respond to emerging infectious diseases.

HIV-1 subtype B spread through cross-border clusters in the Balkans: a molecular analysis in view of incidence trends

OBJECTIVES

To analyze phylogenetic relations and assess the role of cross-border clusters in the spread of HIV-1 subtype B across the Balkans, given the general trends of new HIV diagnoses in seven Balkan countries.

DESIGN

Retrospective phylogenetic and trend analysis.

METHODS

In-depth phylogenetic, phylodynamic and phylogeographic analysis performed on 2415 HIV-1 subtype B sequences from 1999 to 2019 using maximal likelihood and Bayesian methods. The joinpoint regression analysis of new HIV diagnoses by country and modes of transmission using 2004-2019 ECDC data.

RESULTS

Ninety-three HIV-1 Subtype B transmission clusters (68% of studied sequences) were detected of which four cross-border clusters (11% of studied sequences). Phylodynamic analysis showed activity of cross-border clusters up until the mid-2000s, with a subsequent stationary growth phase. Phylogeography analyses revealed reciprocal spread patterns between Serbia, Slovenia and Montenegro and several introductions to Romania from these countries and Croatia. The joinpoint analysis revealed a reduction in new HIV diagnoses in Romania, Greece and Slovenia, whereas an increase in Serbia, Bulgaria, Croatia and Montenegro, predominantly among MSM.

CONCLUSION

Differing trends of new HIV diagnoses in the Balkans mirror differences in preventive policies implemented in participating countries. Regional spread of HIV within the countries of former Yugoslavia has continued to play an important role even after country break-up, whereas the spread of subtype B through multiple introductions to Romania suggested the changing pattern of travel and migration linked to European integration of Balkan countries in the early 2000s.

Genetic Diversity Evaluation and Conservation of Topmouth Culter (Culter alburnus) Germplasm in Five River Basins in China

To study the genetic diversity of Culter alburnus (C. alburnus) populations, we analyzed the genetic diversity of five C. alburnus populations from Songhua Lake (SH), Huaihe River (HH), Changjiang River (CJ), Taihu Lake (TH), and Gehu Lake (GH) based on mitochondrial COI gene sequences. The results showed that the average contents of bases T, C, A, and G in the 526 bp COI gene sequence were 25.3%, 18.1%, 28.1%, and 28.6%, respectively, which showed AT bias. A total of 115 polymorphic sites were detected in the five populations, and 11 haplotypes (Hap) were defined. The nucleotide diversity (P_i) of the five populations ranged from 0.00053 to 0.01834, and the haplotype diversity (H_d) ranged from 0.280 to 0.746, with the highest genetic diversity in the TH population, followed by the SH population, with lower genetic diversity in the HH, CJ and GH populations. The analysis of the fixation index (Fst) and the genetic distance between populations showed that there was significant genetic differentiation between the SH population and the other populations, and the genetic distances between all of them were far; the genetic diversity within populations was higher than that between populations. Neutral tests, mismatch distributions, and Bayesian skyline plot (BSP) analyses showed that the C. alburnus populations have not experienced population expansion and are relatively stable in historical dynamics.

Phylogeography of the mesic-adapted striped mouse, Rhabdomys dilectus chakae (Rodentia: Muridae) in forest margins of the Eastern Cape and southern KwaZulu-Natal provinces of South Africa

Rhabdomys is a genus that occupies a variety of habitats, including forest margins. Among the Rhabdomys taxa, Rhabdomys dilectus chakae has a distribution that covers the eastern seaboard of South Africa, with a poorly defined divergence date from its sister taxon Rhabdomys dilectus dilectus. Here, we study three mitochondrial markers (cytochrome b, cytochrome c oxidase I and partial control region) of R. d. chakae across the Eastern Cape and KwaZulu-Natal forests of South Africa, aiming to determine the cladogenesis (molecular dating) and effective population size of this subspecies through time, in addition to its cladogenesis in relationship to other species within the genus. A phylogenetic analysis revealed six clades within Rhabdomys, confirming that our study area is occupied solely by R. d. chakae, to the exclusion of other recognized sibling species. A fossil-calibrated Bayesian relaxed molecular clock estimated a recent split between R. d. chakae and R. d. dilectus ~1.4 ± 0.35 Mya and between two Rhabdomys pumilio groups, coastal A and B, at 1.16 ± 0.44 Mya. Coalescent Bayesian skyline plots revealed a stable population of R. d. chakae in the study area that was in slow decline until 2500 years ago, when there was an expansion in the late Holocene. Radiation within Rhabdomys dates as far back as 4.27 Mya, and subsequent demographic fluctuations primarily reflect palaeoclimatic changes.

Geography vs. past climate: the drivers of population genetic structure of the Himalayan langur

Background

Contemporary species distribution, genetic diversity and evolutionary history in many taxa are shaped by both historical and current climate as well as topography. The Himalayas show a huge variation in topography and climatic conditions across its entire range, and have experienced major climatic fluctuations in the past. However, very little is known regarding how this heterogenous landscape has moulded the distribution of Himalayan fauna. A recent study examined the effect of these historical events on the genetic diversity of the Himalayan langurs in Nepal Himalaya. However, this study did not include the samples from the Indian Himalayan region (IHR). Therefore, here we revisit the questions addressed in the previous study with a near complete sampling from the IHR, along with the samples from the Nepal Himalaya. We used the mitochondrial Cytochrome-b (Cyt-b, 746 bp) region combined with multiple phylogeographic analyses and palaeodistribution modelling.

Results

Our dataset contained 144 sequences from the IHR as well as the Nepal Himalaya. Phylogenetic analysis showed a low divergent western clade nested within high divergent group of eastern lineages and in the network analysis we identified 22 haplotypes over the entire distribution range of the Himalayan langurs. Samples from the Nepal Himalaya showed geographically structured haplotypes corresponding to different river barriers, whereas samples from IHR showed star-like topology with no structure. Our statistical phylogeography analysis using diyABC supported the model of east to west colonisation of these langurs with founder event during colonisation. Analysis of demographic history showed that the effective population size of the Himalayan langurs decreased at the onset of last glacial maximum (LGM) and started increasing post LGM. The palaeodistribution modelling showed that the extent of suitable habitat shifted from low elevation central Nepal, and adjoining parts of north India, during LGM to the western Himalaya at present.

Conclusion

The current genetic diversity and distribution of Himalayan langurs in the Nepal Himalaya has been shaped by river barriers, whereas the rivers in the IHR had relatively less time to act as a strong genetic barrier after the recent colonisation event. Further, the post LGM expansion could have had confounding effect on Himalayan langur population structure in both Nepal Himalaya and IHR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02054-1.

The Impact of Fast Radiation on the Phylogeny of Bactrocera Fruit Flies as Revealed by Multiple Evolutionary Models and Mutation Rate-Calibrated Clock

Several true fruit flies (Tephritidae) cause major damage to agriculture worldwide. Among them, species of the genus Bactrocera are extensively studied to understand the traits associated with their invasiveness and ecology. Comparative approaches based on a reliable phylogenetic framework are particularly effective, but several nodes of the Bactrocera phylogeny are still controversial, especially concerning the reciprocal affinities of the two major pests B. dorsalis and B. tryoni. Here, we analyzed a newly assembled genomic-scaled dataset using different models of evolution to infer a phylogenomic backbone of ten representative Bactrocera species and two outgroups. We further provide the first genome-scaled inference of their divergence by calibrating the clock using fossil records and the spontaneous mutation rate. The results reveal a closer relationship of B. dorsalis with B. latifrons than to B. tryoni, contrary to what was previously supported by mitochondrial-based phylogenies. By employing coalescent-aware and heterogeneous evolutionary models, we show that this incongruence likely derives from a hitherto undetected systematic error, exacerbated by incomplete lineage sorting and possibly hybridization. This agrees with our clock analysis, which supports a rapid and recent radiation of the clade to which B. dorsalis, B. latifrons and B. tryoni belong. These results provide a new picture of Bactrocera phylogeny that can serve as the basis for future comparative analyses.

Antibiotic-resistant organisms establish reservoirs in new hospital built environments and are related to patient blood infection isolates

Background

Healthcare-associated infections due to antibiotic-resistant organisms pose an acute and rising threat to critically ill and immunocompromised patients. To evaluate reservoirs of antibiotic-resistant organisms as a source of transmission to patients, we interrogated isolates from environmental surfaces, patient feces, and patient blood infections from an established and a newly built intensive care unit.

Methods

We used selective culture to recover 829 antibiotic-resistant organisms from 1594 environmental and 72 patient fecal samples, in addition to 81 isolates from blood cultures. We conducted antibiotic susceptibility testing and short- and long-read whole genome sequencing on recovered isolates.

Results

Antibiotic-resistant organism burden is highest in sink drains compared to other surfaces. Pseudomonas aeruginosa is the most frequently cultured organism from surfaces in both intensive care units. From whole genome sequencing, different lineages of P. aeruginosa dominate in each unit; one P. aeruginosa lineage of ST1894 is found in multiple sink drains in the new intensive care unit and 3.7% of blood isolates analyzed, suggesting movement of this clone between the environment and patients.

Conclusions

These results highlight antibiotic-resistant organism reservoirs in hospital built environments as an important target for infection prevention in hospitalized patients.

Mutation in Hemagglutinin Antigenic Sites in Influenza A pH1N1 Viruses from 2015-2019 in the United States Mountain West, Europe, and the Northern Hemisphere

H1N1 influenza A virus is a respiratory pathogen that undergoes antigenic shift and antigenic drift to improve viral fitness. Tracking the evolutionary trends of H1N1 aids with the current detection and the future response to new viral strains as they emerge. Here, we characterize antigenic drift events observed in the hemagglutinin (HA) sequence of the pandemic H1N1 lineage from 2015-2019. We observed the substitutions S200P, K147N, and P154S, together with other mutations in structural, functional, and/or epitope regions in 2015-2019 HA protein sequences from the Mountain West region of the United States, the larger United States, Europe, and other Northern Hemisphere countries. We reconstructed multiple phylogenetic trees to track the relationships and spread of these mutations and tested for evidence of selection pressure on HA. We found that the prevalence of amino acid substitutions at positions 147, 154, 159, 200, and 233 significantly changed throughout the studied geographical regions between 2015 and 2019. We also found evidence of coevolution among a subset of these amino acid substitutions. The results from this study could be relevant for future epidemiological tracking and vaccine prediction efforts. Similar analyses in the future could identify additional sequence changes that could affect the pathogenicity and/or infectivity of this virus in its human host.

Phylogeographical Analyses of a Relict Fern of Palaeotropical Flora (Vandenboschia speciosa): Distribution and Diversity Model in Relation to the Geological and Climate Events of the Late Miocene and Early Pliocene

Fern phylogeographic studies have mostly focused on the influence of the Pleistocene climate on fern distributions and the prevalence of long-distance dispersal. The effect of pre-Pleistocene events on the distributions of fern species is largely unexplored. Here, we elucidate a hypothetical scenario for the evolutionary history of Vandenboschia speciosa, hypothesised to be of Tertiary palaeotropical flora with a peculiar perennial gametophyte. We sequenced 40 populations across the species range in one plastid region and two variants of the nuclear gapCp gene and conducted time-calibrated phylogenetic, phylogeographical, and species distribution modelling analyses. Vandenboschia speciosa is an allopolyploid and had a Tertiary origin. Late Miocene aridification possibly caused the long persistence in independent refugia on the Eurosiberian Atlantic and Mediterranean coasts, with the independent evolution of gene pools resulting in two evolutionary units. The Cantabrian Cornice, a major refugium, could also be a secondary contact zone during Quaternary glacial cycles. Central European populations resulted from multiple post-glacial, long-distance dispersals. Vandenboschia speciosa reached Macaronesia during the Pliocene–Pleistocene, with a phylogeographical link between the Canary Islands, Madeira, and southern Iberia, and between the Azores and northwestern Europe. Our results support the idea that the geological and climate events of the Late Miocene/Early Pliocene shifted Tertiary fern distribution patterns in Europe.

Phylogeography and evolutionary dynamics analysis of porcine delta-coronavirus with host expansion to humans

From 2003 onwards, three pandemics have been caused by coronaviruses: severe acute respiratory syndrome coronavirus (SARS-CoV); middle east respiratory syndrome coronavirus (MERS-CoV); and, most recently, SARS-CoV-2. Notably, all three were transmitted from animals to humans. This would suggest that animals are potential sources of epidemics for humans. The emerging porcine delta-coronavirus was reported to infect children. This is a red flag that marks the ability of PDCoV to break barriers of cross-species transmission to humans. Therefore, we conducted molecular genetic analysis of global clade PDCoV to characterize spatio-temporal patterns of viral diffusion and genetic diversity. PDCoV was classified into three major lineages, according to distribution and phylogenetic analysis of PDCoV. It can be inferred based on the analysis results of the currently known PDCoV strains that PDCoV might originate in Asia. We also selected six special spike amino acid sequences to align and analyze to find seven significant mutation sites. The accumulation of these mutations may enhance dynamic movements, accelerating spike protein membrane fusion events and transmission. Altogether, our study offers a novel insight into the diversification, evolution, and interspecies transmission and origin of PDCoV and emphasizes the need to study the zoonotic potential of the PDCoV and comprehensive surveillance and enhanced biosecurity precautions for PDCoV. This article is protected by copyright. All rights reserved.

Multiple Recent Colonizations of the Australian Region by the Chydorus sphaericus Group (Crustacea: Cladocera)

Biotic introductions are an ongoing disruption for many ecosystems. For passively dispersed freshwater zooplankton, transcontinental introductions have been common but are poorly studied in the southern hemisphere. Here we assess the hypothesis of recent introduction for populations of the Chydorus sphaericus group (Crustacea: Cladocera) in Australia. We analyzed 254 sequences (63 original sequences) from the cytochrome oxidase I region of mitochondrial DNA of Chydorus sp., which included global representation. Three Australian populations were connected with separate clades in the northern hemisphere, suggesting multiple colonization events for Australia. The timescale of the divergences was consistent with recent (Quaternary) dispersal. As Australian populations are exposed to migrating birds from the northern hemisphere, both avian and anthropogenic sources are candidates for dispersal vectors. We concluded that recent cross-hemisphere dispersal in the Chydorus sphaericus group is more common than previously believed.

Variation in global distribution, population structures, and demographic history for four Trichiurus cutlassfishes

Background

Species-specific information on distribution and demographic patterns provides important implications for conservation and fisheries management. However, such information is often lacking for morphologically-similar species, which may lead to biases in the assessments of these species and even decrease effort towards sustainable management. Here, we aimed to uncover the distribution range, population structure and demographic history for four exploited Trichiurus cutlassfishes using genetics. These cutlassfishes contribute substantial global fisheries catch, with a high proportion of catch harvested from the NW Pacific.

Methods

We chose the widely available mitochondrial 16S ribosomal RNA (16S) as the genetic marker for cutlassfishes. We compiled the 16S sequence data from both the GenBank and a survey of trawler catch samples along the NW Pacific coasts 22-39°N. Genealogical relationships within each species was visualized with haplotype networks and potential population differentiations were further evaluated with AMOVA. Demographic histories were estimated using neutrality test, mismatch analysis, and the Bayesian skyline plot. The reconstructed phylogenetic trees were used to delimit and estimate the divergence time of species and included populations.

Results

In each of two cosmopolitan species, T. lepturus and T. nanhaiensis, we observed distinct populations along the coasts of warm oceans; such population differentiation might result from historical geographic barriers in the Pleistocene. In the NW Pacific, four Trichiurus species vary in their distribution habitats, which reflect differential ecological niches among these species. The small-sized T. brevis are primarily found in nearshore habitats; the warm-affiliated T. nanhaiensis are present along the path of the Kuroshio Current; the cold-affiliated T. japonicus spatially diverged from the widely-distributed T. lepturus, with the latter mainly occupy in warmer regions. Despite these differences, a single well-mixing fish stock, thus one management unit, was identified in each of the four species, presumably due to expansion of their population sizes predated the Last Glacial Maximum and a lack of distribution barrier. The most dominant T. japonicus, which have at least one magnitude higher effective population size than the others, show a unique abrupt size expansion event at 75 to 50-kilo years ago when the low sea level occurred during the ice age.

Main conclusions

The demographic history revealed by our genetic analyses advances understanding of the current distribution and population structure for these congeneric species. Moreover, the uncovered population structure provides insight into the assessment and management of these species. Such information complements contemporary knowledge about these species and enables us to forecast their ability to resist future environmental and anthropogenic disturbances.

Phylogenetic revision of the lichenized family Gomphillaceae (Ascomycota: Graphidales) suggests post-K-Pg boundary diversification and phylogenetic signal in asexual reproductive structures

We present the first broad molecular-phylogenetic revision of the lichenized family Gomphillaceae, based on 408 newly generated sequences of the mitochondrial SSU rDNA and nuclear LSU rDNA, representing 342 OTUs. The phylogenetic analysis of 20 out of the 28 currently accepted genera resulted in 48 clades. Twelve genera were resolved as monophyletic: Actinoplaca, Arthotheliopsis, Bullatina, Caleniopsis, Corticifraga, Gomphillus, Gyalectidium, Gyalidea, Jamesiella, Rolueckia, Rubrotricha, and Taitaia. Two genera resulted paraphyletic, namely Aulaxina (including Caleniopsis) and Asterothyrium (including Linhartia). Six genera were in part highly polyphyletic: Aderkomyces, Calenia, Echinoplaca, Gyalideopsis, Psorotheciopsis, and Tricharia. While ascoma morphology and anatomy has traditionally been considered as main character complex to distinguish genera, our study supported the notion that the characteristic asexual anamorph of Gomphillaceae, the so-called hyphophores, are diagnostic for most of the newly recognized clades. As a result, we recognize 26 new genus-level clades, three of which have names available (Microxyphiomyces, Psathyromyces, Spinomyces) and 23 that will require formal description as new genera. We also tested monophyly for 53 species-level names for which two or more specimens were sequenced: 27 were supported as monophyletic and representing a single species, 13 as monophyletic but with an internal topology suggesting cryptic speciation, four as paraphyletic, and nine as polyphyletic. These data suggest that species richness in the family is higher than indicated by the number of accepted names (currently 425); they also confirm that recently refined species concepts reflect species richness better than the broad concepts applied in Santesson's monograph. A divergence time analysis revealed that foliicolous Gomphillaceae diversified after the K-Pg-boundary and largely during the Miocene, a notion supported by limited data available for other common foliicolous lineages such as Chroodiscus (Graphidaceae), Pilocarpaceae, and Porinaceae. This contradicts recent studies suggesting that only macrofoliose Lecanoromycetes exhibit increased diversification rates in the Cenozoic.

Sympatric and independently evolving lineages in the Thoropa miliaris - T. taophora species complex (Anura: Cycloramphidae)

Species delimitation can be challenging and affected by subjectivity. Sibling lineages that occur in sympatry constitute good candidates for species delimitation regardless of the adopted species concept. The Thoropa miliaris + T. taophora species complex exhibits high genetic diversity distributed in several lineages that occur sympatrically in the southeastern Atlantic Forest of Brazil. We used 414 loci obtained by anchored hybrid enrichment to characterize genetic variation in the Thoropa miliaris species group (T. saxatilis, T megatympanum, T. miliaris, and T. taophora), combining assignment analyses with traditional and coalescent phylogeny reconstruction. We also investigated evolutionary independence in co-occurring lineages by estimating gene flow, and validated lineages under the multispecies coalescent. We recovered most previously described lineages as unique populations in assignment analyses; exceptions include two lineages within T. miliaris that are further substructured, and the merging of all T. taophora lineages. We found very low probabilities of gene flow between sympatric lineages, suggesting independent evolution. Species tree inferences and species delimitation yielded resolved relationships and indicate that all lineages constitute putative species that diverged during the Pliocene and Pleistocene, later than previously estimated.

Diversity and Paleodemography of the Addax (Addax nasomaculatus), a Saharan Antelope on the Verge of Extinction

Since the 19th century, the addax (Addax nasomaculatus) has lost approximately 99% of its former range. Along with its close relatives, the blue antelope (Hippotragus leucophaeus) and the scimitar-horned oryx (Oryx dammah), the addax may be the third large African mammal species to go extinct in the wild in recent times. Despite this, the evolutionary history of this critically endangered species remains virtually unknown. To gain insight into the population history of the addax, we used hybridization capture to generate ten complete mitochondrial genomes from historical samples and assembled a nuclear genome. We found that both mitochondrial and nuclear diversity are low compared to other African bovids. Analysis of mitochondrial genomes revealed a most recent common ancestor ~32 kya (95% CI 11–58 kya) and weak phylogeographic structure, indicating that the addax likely existed as a highly mobile, panmictic population across its Sahelo–Saharan range in the past. PSMC analysis revealed a continuous decline in effective population size since ~2 Ma, with short intermediate increases at ~500 and ~44 kya. Our results suggest that the addax went through a major bottleneck in the Late Pleistocene, remaining at low population size prior to the human disturbances of the last few centuries.

Nuclear phylogenies and genomics of a contact zone establish the species rank of Podarcis lusitanicus (Squamata, Lacertidae)

Unravelling when divergent lineages constitute distinct species can be challenging, particularly in complex scenarios combining cryptic diversity and phylogenetic discordances between different types of molecular markers. Combining a phylogenetic approach with the study of contact zones can help to overcome such difficulties. The Podarcis hispanicus species complex has proven to be prosperous in independent evolutionary units, sometimes associated with cryptic diversity. Previous studies have revealed that one of the species of this complex, P. guadarramae, comprises two deeply divergent yet morphologically indistinguishable evolutionary units, currently regarded as subspecies (P. g. guadarramae and P. g. lusitanicus). In this study we used molecular data to address the systematics of the two lineages of Podarcis guadarramae and the closely related P. bocagei. Firstly, we reconstructed the species tree of these three and two additional taxa based on 30 nuclear loci using the multispecies coalescent with and without gene flow. Secondly, we used SNPs obtained from RADseq data to analyze the population structure across the distribution limits P. g. lusitanicus and P. g. guadarramae, and for comparison, a contact zone between P. bocagei and P. g. lusitanicus. Nuclear phylogenetic relationships between these three taxa are clearly difficult to determine due to the influence of gene flow, but our results give little support to the monophyly of P. guadarramae, potentially due to a nearly simultaneous divergence between them. Genetic structure and geographic cline analysis revealed that the two lineages of P. guadarramae replace each other abruptly across the sampled region and that gene flow is geographically restricted, implying the existence of strong reproductive isolation. Podarcis bocagei and P. g. lusitanicus show a similar degree of genetic differentiation and reproductive isolation, with very low levels of admixture in syntopy. These results support that all three forms are equally differentiated and reproductively isolated. In consequence, we conclude that the two former subspecies of Podarcis guadarramae constitute valid, yet cryptic species, that should be referred to as P. lusitanicus and P. guadarramae.

Gene flow in phylogenomics: Sequence capture resolves species limits and biogeography of Afromontane forest endemic frogs from the Cameroon Highlands

Puddle frogs of the Phrynobatrachus steindachneri species complex are a useful group for investigating speciation and phylogeography in Afromontane forests of the Cameroon Volcanic Line, western Central Africa. The species complex is represented by six morphologically relatively cryptic mitochondrial DNA lineages, only two of which are distinguished at the species level - southern P. jimzimkusi and Lake Oku endemic P. njiomock, leaving the remaining four lineages identified as 'P. steindachneri'. In this study, the six mtDNA lineages are subjected to genomic sequence capture analyses and morphological examination to delimit species and to study biogeography. The nuclear DNA data (387 loci; 571,936 aligned base pairs) distinguished all six mtDNA lineages, but the topological pattern and divergence depths supported only four main clades: P. jimzimkusi, P. njiomock, and only two divergent evolutionary lineages within the four 'P. steindachneri' mtDNA lineages. One of the two lineages is herein described as a new species, P. amieti sp. nov. Reticulate evolution (hybridization) was detected within the species complex with morphologically intermediate hybrid individuals placed between the parental species in phylogenomic analyses, forming a ladder-like phylogenetic pattern. The presence of hybrids is undesirable in standard phylogenetic analyses but is essential and beneficial in the network multispecies coalescent. This latter approach provided insight into the reticulate evolutionary history of these endemic frogs. Introgressions likely occurred during the Middle and Late Pleistocene climatic oscillations, due to the cyclic connections (likely dominating during cold glacials) and separations (during warm interglacials) of montane forests. The genomic phylogeographic pattern supports the separation of the southern (Mt. Manengouba to Mt. Oku) and northern mountains at the onset of the Pleistocene. Further subdivisions occurred in the Early Pleistocene, separating populations from the northernmost (Tchabal Mbabo, Gotel Mts.) and middle mountains (Mt. Mbam, Mt. Oku, Mambilla Plateau), as well as the microendemic lineage restricted to Lake Oku (Mt. Oku). This unique model system is highly threatened as all the species within the complex have exhibited severe population declines in the past decade, placing them on the brink of extinction. In addition, Mount Oku is identified to be of particular conservation importance because it harbors three species of this complex. We, therefore, urge for conservation actions in the Cameroon Highlands to preserve their diversity before it is too late.

Total evidence phylogenetic analysis reveals polyphyly of Anostomoides and uncovers an unexpectedly ancient genus of Anostomidae fishes (Characiformes)

The nearly 150 species of Anostomidae comprise one of the most diverse and taxonomically dynamic families of Neotropical freshwater fishes. A recent revision of the enigmatic and poorly diagnosed genus Anostomoides demonstrated that it contains two valid species, each with complicated taxonomic histories; however, that study did not address their phylogenetic placement. Herein, we integrate molecular and morphological data to demonstrate their distant evolutionary relationship, and thus the polyphyly of Anostomoides. While we reconstruct one of the species in a previously hypothesized placement within a clade also containing Laemolyta, Rhytiodus and Schizodon, the other represents a morphologically and genetically distinctive lineage that diverged early in the history of the family. We describe and illustrate the osteology of this remarkable species, discuss the evolutionary implications of its unique suite of features, and use those characteristics to diagnose a new genus that evolved independently of all other known members of the family for approximately 37 Myr.

Phylogenetics and historical biogeography of Solanum section Brevantherum (Solanaceae)

Patterns of diversification in Neotropical plants have been studied intensively over the past decades. Most studies have focused on groups that migrated to and radiated into the Neotropics, however, with little focus on understanding diversification patterns in indigenous Neotropical groups. This study focuses on Solanum section Brevantherum Seithe (Solanaceae), a group of Neotropical nightshade shrubs or treelets defined mostly by terminal inflorescences with long peduncles, plurifoliate sympodial units and porrect-stellate, dendritic-echinoid, or lepidote trichomes. We generated sequences from two nuclear (ITS, waxy) and one plastid marker (trnT-F) to infer phylogenetic relationships under Bayesian and Maximum likelihood approaches. We reconstructed a time-calibrated tree to estimate both the ages of main splits and the ancestral ranges of the lineages. Finally, we carried out a biogeographic stochastic mapping (BSM) analysis to determine the main processes driving current distributions of the group. Results show the non-monophyly of the section as previously recognized and the homoplasy of morphological characters traditionally used to circumscribe it. Two main clades that encompass most species formerly recognized in section Brevantherum are recovered and named as the Erianthum and Abutiloides clades. Divergence time estimates suggest that the Erianthum and Abutiloides clades split around 5.7 Mya in the upper Miocene. Two main dispersal events from the Atlantic rainforest are supported in the Erianthum clade: one dispersal to Mesoamerica and a second dispersal to the Northern Andes. Within the Abutiloides clade, cladogenetic events were restricted to the Andean region. Our BSM analysis suggests within-area speciation and range expansion as the main processes shaping the extant distribution of species of both clades. As no putative morphological synapomorphies can yet be assigned to what could correspond to a new circumscription of Solanum section Brevantherum (with the exclusion of Solanum bullatum Vell. and inclusion of S. inelegans Rusby and four species described since the group last revision) we discourage the continued use of what would be an ambiguous sectional nomenclature.

Chloroplast genomes in Populus (Salicaceae): comparisons from an intensively sampled genus reveal dynamic patterns of evolution

The chloroplast is one of two organelles containing a separate genome that codes for essential and distinct cellular functions such as photosynthesis. Given the importance of chloroplasts in plant metabolism, the genomic architecture and gene content have been strongly conserved through long periods of time and as such are useful molecular tools for evolutionary inferences. At present, complete chloroplast genomes from over 4000 species have been deposited into publicly accessible databases. Despite the large number of complete chloroplast genomes, comprehensive analyses regarding genome architecture and gene content have not been conducted for many lineages with complete species sampling. In this study, we employed the genus Populus to assess how more comprehensively sampled chloroplast genome analyses can be used in understanding chloroplast evolution in a broadly studied lineage of angiosperms. We conducted comparative analyses across Populus in order to elucidate variation in key genome features such as genome size, gene number, gene content, repeat type and number, SSR (Simple Sequence Repeat) abundance, and boundary positioning between the four main units of the genome. We found that some genome annotations were variable across the genus owing in part from errors in assembly or data checking and from this provided corrected annotations. We also employed complete chloroplast genomes for phylogenetic analyses including the dating of divergence times throughout the genus. Lastly, we utilized re-sequencing data to describe the variations of pan-chloroplast genomes at the population level for P. euphratica. The analyses used in this paper provide a blueprint for the types of analyses that can be conducted with publicly available chloroplast genomes as well as methods for building upon existing datasets to improve evolutionary inference.

Myoglobin primary structure reveals multiple convergent transitions to semi-aquatic life in the world's smallest mammalian divers

The speciose mammalian order Eulipotyphla (moles, shrews, hedgehogs, solenodons) combines an unusual diversity of semi-aquatic, semi-fossorial, and fossorial forms that arose from terrestrial forbearers. However, our understanding of the ecomorphological pathways leading to these lifestyles has been confounded by a fragmentary fossil record, unresolved phylogenetic relationships, and potential morphological convergence, calling for novel approaches. The net surface charge of the oxygen-storing muscle protein myoglobin (ZMb), which can be readily determined from its primary structure, provides an objective target to address this question due to mechanistic linkages with myoglobin concentration. Here, we generate a comprehensive 71 species molecular phylogeny that resolves previously intractable intra-family relationships and then ancestrally reconstruct ZMb evolution to identify ancient lifestyle transitions based on protein sequence alone. Our phylogenetically informed analyses confidently resolve fossorial habits having evolved twice in talpid moles and reveal five independent secondary aquatic transitions in the order housing the world's smallest endothermic divers.

Bosminopsis deitersi (Crustacea: Cladocera) as an ancient species group: a revision

Water fleas (Crustacea: Cladocera) of the Family Bosminidae have been studied since the founding of paleolimnology and freshwater ecology. However, one species, Bosminopsis deitersi, stands out for its exceptional multicontinental range and broad ecological requirements. Here we use an integrated morphological and multilocus genetic approach to address the species problem in B. deitersi. We analyzed 32 populations of B. deitersi s. lat. Two nuclear and two mitochondrial loci were used to carry out the bGMYC, mPTP and STACEY algorithms for species delimitation. Detailed morphological study was also carried out across continents. The evidence indicated a widely distributed cryptic species in the Old World (Bosminopsis zernowi) that is genetically divergent from B. deitersi s.str. We revised the taxonomy and redescribed the species in this complex. Our sampling indicated that B. zernowi had weak genetic differentiation across its range. A molecular clock and biogeographic analysis with fossil calibrations suggested a Mesozoic origin for the Bosminopsis deitersi group. Our evidence rejects the single species hypothesis for B. deitersi and is consistent with an ancient species group (potentially Mesozoic) that shows marked morphological conservation. The family Bosminidae, then, has examples of both rapid morphological evolution (Holocene Bosmina), and morphological stasis (Bosminopsis).

Genetic evidence for the origin of Aedes aegypti, the yellow fever mosquito, in the southwestern Indian Ocean

Aedes aegypti is among the best-studied mosquitoes due to its critical role as a vector of human pathogens and ease of laboratory rearing. Until now, this species was thought to have originated in continental Africa, and subsequently colonized much of the world following the establishment of global trade routes. However, populations of this mosquito on the islands in the southwestern Indian Ocean (SWIO), where the species occurs with its nearest relatives referred to as the Aegypti Group, have received little study. We re-evaluated the evolutionary history of Ae. aegypti and these relatives, using three data sets: nucleotide sequence data, 18,489 SNPs and 12 microsatellites. We found that: (a) the Aegypti Group diverged 16 MYA (95% HPD: 7-28 MYA) from its nearest African/Asian ancestor; (b) SWIO populations of Ae. aegypti are basal to continental African populations; (c) after diverging 7 MYA (95% HPD: 4-15 MYA) from its nearest formally described relative (Ae. mascarensis), Ae. aegypti moved to continental Africa less than 85,000 years ago, where it recently (<1,000 years ago) split into two recognized subspecies Ae. aegypti formosus and a human commensal, Ae. aegypti aegypti; (d) the Madagascar samples form a clade more distant from all other Ae. aegypti than the named species Ae. mascarensis, implying that Madagascar may harbour a new cryptic species; and (e) there is evidence of introgression between Ae. mascarensis and Ae. aegypti on Réunion, and between the two subspecies elsewhere in the SWIO, a likely consequence of recent introductions of domestic Ae. aegypti aegypti from Asia.

Adaptive Metropolis-coupled MCMC for BEAST 2

With ever more complex models used to study evolutionary patterns, approaches that facilitate efficient inference under such models are needed. Metropolis-coupled Markov chain Monte Carlo (MCMC) has long been used to speed up phylogenetic analyses and to make use of multi-core CPUs. Metropolis-coupled MCMC essentially runs multiple MCMC chains in parallel. All chains are heated except for one cold chain that explores the posterior probability space like a regular MCMC chain. This heating allows chains to make bigger jumps in phylogenetic state space. The heated chains can then be used to propose new states for other chains, including the cold chain. One of the practical challenges using this approach, is to find optimal temperatures of the heated chains to efficiently explore state spaces. We here provide an adaptive Metropolis-coupled MCMC scheme to Bayesian phylogenetics, where the temperature difference between heated chains is automatically tuned to achieve a target acceptance probability of states being exchanged between individual chains. We first show the validity of this approach by comparing inferences of adaptive Metropolis-coupled MCMC to MCMC on several datasets. We then explore where Metropolis-coupled MCMC provides benefits over MCMC. We implemented this adaptive Metropolis-coupled MCMC approach as an open source package licenced under GPL 3.0 to the Bayesian phylogenetics software BEAST 2, available from https://github.com/nicfel/CoupledMCMC.

Cryptic species in White Cloud Mountain minnow, Tanichthys albonubes: Taxonomic and conservation implications

Cryptic species describe two or more species that had mistakenly been considered to be a single species, a phenomenon that has been found throughout the tree of life. Recognizing cryptic species is key to estimating the real biodiversity of the world and understanding evolutionary processes. Molecular methods present an unprecedented opportunity for biologists to question whether morphologically similar populations are actually cryptic species. The minnow Tanichthys albonubes is a critically endangered freshwater fish and was classified as a second-class state-protected animal in China. Previous studies have revealed highly divergent lineages with similar morphological characters in this species. Herein, we tested for cryptic species across the ranges of all known wild populations of this minnow. Using multilocus molecular (one mitochondrial gene, two nuclear genes and 13 microsatellite loci) and morphological data for 230 individuals from eight populations, we found deep genetic divergence among these populations with subtle morphological disparity. Morphological examination found variance among these populations in the number of branched anal-fin rays. Based on genetic data, we inferred eight monophyletic groups that were well supported by haplotype network and population clustering analyses. Species delimitation methods suggested eight putative species in the T. albonubes complex. Molecular dating suggested that these cryptic species diverged in the period from the Pliocene to the Pleistocene. Based on these findings, we propose the existence of seven cryptic species in the T. albonubes complex. Our results highlight the need for a taxonomic revision of Tanichthys. What is more, the conservation status of and conservation strategies for the T. albonubes complex should be reassessed as soon as possible.

Transmission route and introduction of pandemic SARS-CoV-2 between China, Italy, and Spain

We present a phylodynamic and phylogeographic analysis of this new severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) virus in this report. A tree of maximum credibility was constructed using the 72 entire genome sequences of this virus, from the three countries (China, Italy, and Spain) available as of 26 March 2020 on the GISAID reference frame. To schematize the current SARS‐CoV‐2 migration scenario between and within the three countries chosen, using the multitype bearth‐death model implemented in BEAST2. Bayesian phylogeographic reconstruction shows that SARS‐CoV‐2 has a rate of evolution of 2.11 × 10⁻³ per sites per year (95% highest posterior density: 1.56 × 10⁻³ to 3.89 × 10⁻³), and a geographic origin in Shanghai, where time until the most recent common ancestor (tMRCA) emerged, according to the analysis of the molecular clock, around 13 November 2019. While for Italy and Spain, there are two tMRCA for each country, which agree with the assumption of several introductions for these countries. That explains also this very short period of subepidermal circulation before the recent events. A total of 8 (median) migration events occurred during this short period, the largest proportion of which (6 events [75%]) occurred from Shanghai (China) to Spain and from Italy to Spain. Such events are marked by speeds of migration that are comparatively lower as compared with that from Shanghai to Italy. Shanghai's R₀ and Italy's are closer to each other, though Spain's is slightly higher. All these results allow us to conclude the need for an automatic system of mixed, molecular and classical epidemiological surveillance, which could play a role in this global surveillance of public health and decision‐making.

The multitype bearth‐death model implemented in BEAST2 is used.

The aim of this study was to estimate the migration rate of the SRAS‐CoV‐2 between and within the three countries chosen (CHINA, ITALY and SPAIN).

Identifying the R0 for each subpopulation.

Analyzing the number of transmission between these three geographic locations.

Bayesian inference of reassortment networks reveals fitness benefits of reassortment in human influenza viruses

Genetic recombination processes, such as reassortment, make it complex or impossible to use standard phylogenetic and phylodynamic methods. This is due to the fact that the shared evolutionary history of individuals has to be represented by a phylogenetic network instead of a tree. We therefore require novel approaches that allow us to coherently model these processes and that allow us to perform inference in the presence of such processes. Here, we introduce an approach to infer reassortment networks of segmented viruses using a Markov chain Monte Carlo approach. Our approach allows us to study different aspects of the reassortment process and allows us to show fitness benefits of reassortment events in seasonal human influenza viruses.

Reassortment is an important source of genetic diversity in segmented viruses and is the main source of novel pathogenic influenza viruses. Despite this, studying the reassortment process has been constrained by the lack of a coherent, model-based inference framework. Here, we introduce a coalescent-based model that allows us to explicitly model the joint coalescent and reassortment process. In order to perform inference under this model, we present an efficient Markov chain Monte Carlo algorithm to sample rooted networks and the embedding of phylogenetic trees within networks. This algorithm provides the means to jointly infer coalescent and reassortment rates with the reassortment network and the embedding of segments in that network from full-genome sequence data. Studying reassortment patterns of different human influenza datasets, we find large differences in reassortment rates across different human influenza viruses. Additionally, we find that reassortment events predominantly occur on selectively fitter parts of reassortment networks showing that on a population level, reassortment positively contributes to the fitness of human influenza viruses.

Predicting the short-term success of human influenza virus variants with machine learning

Seasonal influenza viruses are constantly changing and produce a different set of circulating strains each season. Small genetic changes can accumulate over time and result in antigenically different viruses; this may prevent the body's immune system from recognizing those viruses. Due to rapid mutations, in particular, in the haemagglutinin (HA) gene, seasonal influenza vaccines must be updated frequently. This requires choosing strains to include in the updates to maximize the vaccines' benefits, according to estimates of which strains will be circulating in upcoming seasons. This is a challenging prediction task. In this paper, we use longitudinally sampled phylogenetic trees based on HA sequences from human influenza viruses, together with counts of epitope site polymorphisms in HA, to predict which influenza virus strains are likely to be successful. We extract small groups of taxa (subtrees) and use a suite of features of these subtrees as key inputs to the machine learning tools. Using a range of training and testing strategies, including training on H3N2 and testing on H1N1, we find that successful prediction of future expansion of small subtrees is possible from these data, with accuracies of 0.71-0.85 and a classifier 'area under the curve' 0.75-0.9.

Disappearance of Icelandic Walruses Coincided with Norse Settlement

There is a growing body of evidence demonstrating the impacts of human arrival in new “pristine” environments, including terrestrial habitat alterations and species extinctions. However, the effects of marine resource utilization prior to industrialized whaling, sealing, and fishing have largely remained understudied. The expansion of the Norse across the North Atlantic offers a rare opportunity to study the effects of human arrival and early exploitation of marine resources. Today, there is no local population of walruses on Iceland, however, skeletal remains, place names, and written sources suggest that walruses existed, and were hunted by the Norse during the Settlement and Commonwealth periods (870–1262 AD). This study investigates the timing, geographic distribution, and genetic identity of walruses in Iceland by combining historical information, place names, radiocarbon dating, and genomic analyses. The results support a genetically distinct, local population of walruses that went extinct shortly after Norse settlement. The high value of walrus products such as ivory on international markets likely led to intense hunting pressure, which—potentially exacerbated by a warming climate and volcanism—resulted in the extinction of walrus on Iceland. We show that commercial hunting, economic incentives, and trade networks as early as the Viking Age were of sufficient scale and intensity to result in significant, irreversible ecological impacts on the marine environment. This is to one of the earliest examples of local extinction of a marine species following human arrival, during the very beginning of commercial marine exploitation.

Recent advances in Bayesian inference of isolation-with-migration models

Isolation-with-migration (IM) models have become popular for explaining population divergence in the presence of migrations. Bayesian methods are commonly used to estimate IM models, but they are limited to small data analysis or simple model inference. Recently three methods, IMa3, MIST, and AIM, resolved these limitations. Here, we describe the major problems addressed by these three software and compare differences among their inference methods, despite their use of the same standard likelihood function.