Outside the pattern: Evolution of the genital asymmetry in Saicinae (Hemiptera: Heteroptera: Reduviidae)

Despite genital structures in insects being consistently important as systematic and taxonomy evidence, within assassin bugs (Heteroptera, Reduviidae) at least, the male and female genitalic structures of several subfamilies are poorly or totally unknown. The genital structure is mostly symmetric within Saicinae genera, but male genital asymmetry has been recorded in Gallobelgicus, Polytoxus, and recently Pseudosaica. The subfamily has been considered as closely related to Emesinae and Visayanocorinae but this hypothesis has never been tested using a comprehensive taxa sampling, being a constraint to test morphological traits changes or exploring characters' evolution hypotheses. Here, we compiled a morphological data set of 170 characters that includes external morphological characters and genitalia of both sexes of Saicinae which was analyzed cladistically including 55 terminals, comprising 16 genera (64% of the generic diversity), 43 species of Saicinae and 12 outgroups. Saicinae was recovered as polyphyletic, Saicireta correntina is recovered as sister-species of Empicoris armatus + Collartida (Emesinae), Oncerotrachelus, Carayonia (Visayanocorinae), and the Clade Saicinae sensu stricto. Carayonia orientalis is recovered as sister-species of Saicinae sensu stricto. Ancestral state reconstruction of symmetry of the male genitalia shows an ancestor with symmetric male genitalia, two independent emergences of asymmetrical male genitalia within Saicinae sensu stricto, and the asymmetrical endosomal sclerites appearing before the other asymmetric traits.

Modeling Drug Resistance Emergence and Transmission in HIV-1 in the UK

A deeper understanding of HIV-1 transmission and drug resistance mechanisms can lead to improvements in current treatment policies. However, the rates at which HIV-1 drug resistance mutations (DRMs) are acquired and which transmitted DRMs persist are multi-factorial and vary considerably between different mutations. We develop a method for the estimation of drug resistance acquisition and transmission patterns. The method uses maximum likelihood ancestral character reconstruction informed by treatment roll-out dates and allows for the analysis of very large datasets. We apply our method to transmission trees reconstructed on the data obtained from the UK HIV Drug Resistance Database to make predictions for known DRMs. Our results show important differences between DRMs, in particular between polymorphic and non-polymorphic DRMs and between the B and C subtypes. Our estimates of reversion times, based on a very large number of sequences, are compatible but more accurate than those already available in the literature, with narrower confidence intervals. We consistently find that large resistance clusters are associated with polymorphic DRMs and DRMs with long loss times, which require special surveillance. As in other high-income countries (e.g., Switzerland), the prevalence of sequences with DRMs is decreasing, but among these, the fraction of transmitted resistance is clearly increasing compared to the fraction of acquired resistance mutations. All this indicates that efforts to monitor these mutations and the emergence of resistance clusters in the population must be maintained in the long term.

Characterization of the Complete Mitochondrial Genome of Eight Diurnal Hawkmoths (Lepidoptera: Sphingidae): New Insights into the Origin and Evolution of Diurnalism in Sphingids

In this study, the mitochondrial genomes of 22 species from three subfamilies in the Sphingidae were sequenced, assembled, and annotated. Eight diurnal hawkmoths were included, of which six were newly sequenced (Hemaris radians, Macroglossum bombylans, M. fritzei, M. pyrrhosticta, Neogurelca himachala, and Sataspes xylocoparis) and two were previously published (Cephonodes hylas and Macroglossum stellatarum). The mitochondrial genomes of these eight diurnal hawkmoths were comparatively analyzed in terms of sequence length, nucleotide composition, relative synonymous codon usage, non-synonymous/synonymous substitution ratio, gene spacing, and repeat sequences. The mitogenomes of the eight species, ranging in length from 15,201 to 15,461 bp, encode the complete set of 37 genes usually found in animal mitogenomes. The base composition of the mitochondrial genomes showed A+T bias. The most commonly used codons were UUA (Leu), AUU (Ile), UUU (Phe), AUA (Met), and AAU (Asn), whereas GCG (Ala) and CCG (Pro) were rarely used. A phylogenetic tree of Sphingidae was constructed based on both maximum likelihood and Bayesian methods. We verified the monophyly of the four current subfamilies of Sphingidae, all of which had high support. In addition, we performed divergence time estimation and ancestral character reconstruction analyses. Diurnal behavior in hawkmoths originated 29.19 million years ago (Mya). It may have been influenced by the combination of herbaceous flourishing, which occurred 26-28 Mya, the uplift of the Tibetan Plateau, and the large-scale evolution of bats in the Oligocene to Pre-Miocene. Moreover, diurnalism in hawkmoths had multiple independent origins in Sphingidae.

Untargeted Metabolomics Sheds Light on the Diversity of Major Classes of Secondary Metabolites in the Malpighiaceae Botanical Family

Natural products produced by plants are one of the most investigated natural sources, which substantially contributed to the development of the natural products field. Even though these compounds are widely explored, the literature still lacks comprehensive investigations aiming to explore the evolution of secondary metabolites produced by plants, especially if classical methodologies are employed. The development of sensitive hyphenated techniques and computational tools for data processing has enabled the study of large datasets, being valuable assets for chemosystematic studies. Here, we describe a strategy for chemotaxonomic investigations using the Malpighiaceae botanical family as a model. Our workflow was based on MS/MS untargeted metabolomics, spectral searches, and recently described in silico classification tools, which were mapped into the latest molecular phylogeny accepted for this family. The metabolomic analysis revealed that different ionization modes and extraction protocols significantly impacted the chemical profiles, influencing the chemotaxonomic results. Spectral searches within public databases revealed several clades or genera-specific molecular families, being potential chemical markers for these taxa, while the in silico classification tools were able to expand the Malpighiaceae chemical space. The classes putatively annotated were used for ancestral character reconstructions, which recovered several classes of metabolites as homoplasies (i.e., non-exclusive) or synapomorphies (i.e., exclusive) for all sampled clades and genera. Our workflow combines several approaches to perform a comprehensive evolutionary chemical study. We expect it to be used on further chemotaxonomic investigations to expand chemical knowledge and reveal biological insights for compounds classes in different biological groups.

Is evolution faster at ecotones? A test using rates and tempo of diet transitions in Neotropical Sigmodontinae (Rodentia, Cricetidae)

We evaluated whether evolution is faster at ecotones as niche shifts may be needed to persist under unstable environment. We mapped diet evolution along the evolutionary history of 350 sigmodontine species. Mapping was used in three new tip-based metrics of trait evolution - Transition Rates, Stasis Time, and Last Transition Time - which were spatialized at the assemblage level (aTR, aST, aTL). Assemblages were obtained by superimposing range maps on points located at core and ecotone of the 93 South American ecoregions. Using Linear Mixed Models, we tested whether ecotones have species with more changes from the ancestral diet (higher aTR), have maintained the current diet for a shorter time (lower aST), and have more recent transitions to the current diet (lower aLT) than cores. We found lower aTR, and higher aST and aLT at ecotones than at cores. Although ecotones are more heterogeneous, both environmentally and in relation to selection pressures they exert on organisms, ecotone species change little from the ancestral diet as generalist habits are necessary toward feeding in ephemeral environments. The need to incorporate phylogenetic uncertainty in tip-based metrics was evident from large uncertainty detected. Our study integrates ecology and evolution by analyzing how fast trait evolution is across space.

African Swine Fever Virus Circulation between Tanzania and Neighboring Countries: A Systematic Review and Meta-Analysis

For over 100 years after the description of the first case of African swine fever (ASF) in Kenya, ASF virus (ASFV) cross-border spread in eastern and southern Africa has not been fully investigated. In this manuscript, we reviewed systematically the available literature on molecular epidemiology of ASF in Tanzania and its eight neighboring countries in order to establish the transmission dynamics of ASFV between these countries. Data were retrieved from World Animal Health Information System (WAHIS), Google Scholar, PubMed, Scopus, and CrossRef databases, using the recommendations of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines and reviewed to document ASF outbreaks and ASFV genotypes distribution. Using phylogeographic approach applied to ASFV p72 sequence dataset, the evolutionary history and the dispersal pattern of the ASFV strains were assessed. From 2005 to 2019, a total of 1588 ASF outbreaks affecting 341,742 cases that led to 302,739 domestic pig deaths were reported. The case fatality rates (CFR) varied from 15.41% to 98.95% with an overall CFR of 88.58%. Fifteen different p72 ASFV genotypes were reported and the time to the most recent common ancestor (TMRCA) for ASFV strains dated back to 1652.233 (1626.473, 1667.735) with an evolutionary rate of 4.805 × 10⁻⁵ (2.5857 × 10⁻⁵, 9.7789 × 10⁻⁵). Phylogeographic dispersal analysis revealed several transboundary spread events of ASFV strains between these countries. These results suggest persistent circulation of ASFV in these countries and advocate for more research to improve our understanding of the transmission dynamics of the virus and for a regional approach to mitigate the spread of ASFV.

Niche partitioning among three snail-eating snakes revealed by dentition asymmetry and prey specialisation

The level of dentition asymmetry in snail-eating snakes may reflect their prey choice and feeding efficiency on asymmetric land snails. The three species of Pareas snakes (Squamata: Pareidae) in Taiwan, which form partially sympatric distribution on the island, provide a potential case to test the hypothesis of niche partitioning and character displacement with regard to dentition asymmetry and specialisation in feeding behaviour. In this study, behavioural experiments confirmed that P. formosensis feeds exclusively on slugs, whereas P. atayal and P. komaii consumed both. However, P. atayal more efficiently preys on land snails than P. komaii, exhibiting a shorter handling time and fewer mandibular retractions. Micro-CT and ancestral character reconstruction demonstrated the lowest asymmetry in P. formosensis (the slug specialist), the highest dentition asymmetry in P. atayal (the land snail specialist) and flexibility in P. komaii (the niche switcher): increased dentition asymmetry when sympatrically distributed with the slug eater (character displacement), and decreased asymmetry when living alone (ecological niche release). Ecological niche modelling showed that the distribution of P. formosensis is associated with the presence of slugs, while that of P. atayal could be explained by the land snails. Combining the results from morphology, phylogeny, behavioural experiments and ecological niche modelling, we showed that competition in the sympatric region might have facilitated character displacement among congeners, while the absence of competition in allopatric region has led to ecological niche release.

Phylogeography of the marine pathogen, Vibrio vulnificus, revealed the ancestral scenarios of its evolution

Vibrio vulnificus is the leading cause of seafood‐associated deaths worldwide. Despite the growing knowledge about the population structure of V. vulnificus, the evolutionary history and the ancestral relationships of strains isolated from various regions around the world have not been determined. Using the largest collection of sequence and isolate data of V. vulnificus to date, we applied ancestral character reconstruction to study the phylogeography of V. vulnificus. Multilocus sequence typing data from 10 housekeeping genes were used for the inference of ancestral states and reconstruction of the evolutionary history. The findings showed that the common ancestor of all V. vulnificus populations originated from East Asia, and later evolved into two main clusters that spread with time and eventually evolved into distinct populations in different parts of the world. While we found no meaningful insights concerning the evolution of V. vulnificus populations in the Middle East; however, we were able to reconstruct the ancestral scenarios of its evolution in East Asia, North America, and Western Europe.

A Fast Likelihood Method to Reconstruct and Visualize Ancestral Scenarios

The reconstruction of ancestral scenarios is widely used to study the evolution of characters along phylogenetic trees. One commonly uses the marginal posterior probabilities of the character states, or the joint reconstruction of the most likely scenario. However, marginal reconstructions provide users with state probabilities, which are difficult to interpret and visualize, whereas joint reconstructions select a unique state for every tree node and thus do not reflect the uncertainty of inferences. We propose a simple and fast approach, which is in between these two extremes. We use decision-theory concepts (namely, the Brier score) to associate each node in the tree to a set of likely states. A unique state is predicted in tree regions with low uncertainty, whereas several states are predicted in uncertain regions, typically around the tree root. To visualize the results, we cluster the neighboring nodes associated with the same states and use graph visualization tools. The method is implemented in the PastML program and web server. The results on simulated data demonstrate the accuracy and robustness of the approach. PastML was applied to the phylogeography of Dengue serotype 2 (DENV2), and the evolution of drug resistances in a large HIV data set. These analyses took a few minutes and provided convincing results. PastML retrieved the main transmission routes of human DENV2 and showed the uncertainty of the human-sylvatic DENV2 geographic origin. With HIV, the results show that resistance mutations mostly emerge independently under treatment pressure, but resistance clusters are found, corresponding to transmissions among untreated patients.

Divergence times and the evolution of morphological complexity in an early land plant lineage (Marchantiopsida) with a slow molecular rate

We present a complete generic-level phylogeny of the complex thalloid liverworts, a lineage that includes the model system Marchantia polymorpha. The complex thalloids are remarkable for their slow rate of molecular evolution and for being the only extant plant lineage to differentiate gas exchange tissues in the gametophyte generation. We estimated the divergence times and analyzed the evolutionary trends of morphological traits, including air chambers, rhizoids and specialized reproductive structures. A multilocus dataset was analyzed using maximum likelihood and Bayesian approaches. Relative rates were estimated using local clocks. Our phylogeny cements the early branching in complex thalloids. Marchantia is supported in one of the earliest divergent lineages. The rate of evolution in organellar loci is slower than for other liverwort lineages, except for two annual lineages. Most genera diverged in the Cretaceous. Marchantia polymorpha diversified in the Late Miocene, giving a minimum age estimate for the evolution of its sex chromosomes. The complex thalloid ancestor, excluding Blasiales, is reconstructed as a plant with a carpocephalum, with filament-less air chambers opening via compound pores, and without pegged rhizoids. Our comprehensive study of the group provides a temporal framework for the analysis of the evolution of critical traits essential for plants during land colonization.

Frequent but asymmetric niche shifts in Bulbophyllum orchids support environmental and climatic instability in Madagascar over Quaternary time scales

BACKGROUND

Species or clades may retain or shift their environmental niche space over evolutionary time. Understanding these processes offers insights into the environmental processes fuelling lineage diversification and might also provide information on past range dynamics of ecosystems. However, little is known about the relative contributions of niche conservatism versus niche divergence to species diversification in the tropics. Here, we examined broad-scale patterns of niche evolution within a Pliocene-Pleistocene clade of epiphytic Bulbophyllum orchids (30 spp.) whose collective distribution covers the northwest and eastern forest ecosystems of Madagascar.

RESULTS

Using species occurrence data, ecological niche models, and multivariate analyses of contributing variables, we identified a three-state niche distribution character for the entire clade, coinciding with three major forest biomes viz. phytogeographical provinces in Madagascar: A, Northwest 'Sambirano'; B, 'Eastern Lowlands'; and C, 'Central Highlands'. A time-calibrated phylogeny and Bayesian models of niche evolution were then used to detect general trends in the direction of niche change over the clade's history (≤5.3 Ma). We found highest transitions rates between lowlands (A and B) and (mostly from B) into the highland (C), with extremely low rates out of the latter. Lowland-to-highland transitions occurred frequently during the Quaternary, suggesting that climate-induced vegetational shifts promoted niche transitions and ecological speciation at this time.

CONCLUSIONS

Our results reveal that niche transitions occurred frequently and asymmetrically within this Madagascan orchid clade, and in particular over Quaternary time scales. Intrinsic features germane to Bulbophyllum (e.g., high dispersal ability, drought tolerance, multiple photosynthetic pathways) as well as extrinsic factors (ecological, historical) likely interacted to generate the niche transition patterns observed. In sum, our results support the emerging idea of dramatic environmental and climatic fluctuations in Madagascar during the recent geological past, which overturns the long-held paradigm of long-term stability in tropical forest settings. The generality of the patterns and timings reported here awaits the availability of additional comparative studies in other Madagascan endemics.

Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi

We studied the evolutionary history of the Parmeliaceae (Lecanoromycetes, Ascomycota), one of the largest families of lichen-forming fungi with complex and variable morphologies, also including several lichenicolous fungi. We assembled a six-locus data set including nuclear, mitochondrial and low-copy protein-coding genes from 293 operational taxonomic units (OTUs). The lichenicolous lifestyle originated independently three times in lichenized ancestors within Parmeliaceae, and a new generic name is introduced for one of these fungi. In all cases, the independent origins occurred c. 24 million yr ago. Further, we show that the Paleocene, Eocene and Oligocene were key periods when diversification of major lineages within Parmeliaceae occurred, with subsequent radiations occurring primarily during the Oligocene and Miocene. Our phylogenetic hypothesis supports the independent origin of lichenicolous fungi associated with climatic shifts at the Oligocene-Miocene boundary. Moreover, diversification bursts at different times may be crucial factors driving the diversification of Parmeliaceae. Additionally, our study provides novel insight into evolutionary relationships in this large and diverse family of lichen-forming ascomycetes.

Resolution of Brassicaceae Phylogeny Using Nuclear Genes Uncovers Nested Radiations and Supports Convergent Morphological Evolution

Brassicaceae is one of the most diverse and economically valuable angiosperm families with widely cultivated vegetable crops and scientifically important model plants, such as Arabidopsis thaliana. The evolutionary history, ecological, morphological, and genetic diversity, and abundant resources and knowledge of Brassicaceae make it an excellent model family for evolutionary studies. Recent phylogenetic analyses of the family revealed three major lineages (I, II, and III), but relationships among and within these lineages remain largely unclear. Here, we present a highly supported phylogeny with six major clades using nuclear markers from newly sequenced transcriptomes of 32 Brassicaceae species and large data sets from additional taxa for a total of 55 species spanning 29 out of 51 tribes. Clade A consisting of Lineage I and Macropodium nivale is sister to combined Clade B (with Lineage II and others) and a new Clade C. The ABC clade is sister to Clade D with species previously weakly associated with Lineage II and Clade E (Lineage III) is sister to the ABCD clade. Clade F (the tribe Aethionemeae) is sister to the remainder of the entire family. Molecular clock estimation reveals an early radiation of major clades near or shortly after the Eocene–Oligocene boundary and subsequent nested divergences of several tribes of the previously polytomous Expanded Lineage II. Reconstruction of ancestral morphological states during the Brassicaceae evolution indicates prevalent parallel (convergent) evolution of several traits over deep times across the entire family. These results form a foundation for future evolutionary analyses of structures and functions across Brassicaceae.

Settling down of seasonal migrants promotes bird diversification

How seasonal migration originated and impacted diversification in birds remains largely unknown. Although migratory behaviour is likely to affect bird diversification, previous studies have not detected any effect. Here, we infer ancestral migratory behaviour and the effect of seasonal migration on speciation and extinction dynamics using a complete bird tree of life. Our analyses infer that sedentary behaviour is ancestral, and that migratory behaviour evolved independently multiple times during the evolutionary history of birds. Speciation of a sedentary species into two sedentary daughter species is more frequent than speciation of a migratory species into two migratory daughter species. However, migratory species often diversify by generating a sedentary daughter species in addition to the ancestral migratory one. This leads to an overall higher migratory speciation rate. Migratory species also experience lower extinction rates. Hence, although migratory species represent a minority (18.5%) of all extant birds, they have a higher net diversification rate than sedentary species. These results suggest that the evolution of seasonal migration in birds has facilitated diversification through the divergence of migratory subpopulations that become sedentary, and illustrate asymmetrical diversification as a mechanism by which diversification rates are decoupled from species richness.

Reconstruction of ancestral 16S rRNA reveals mutation bias in the evolution of optimal growth temperature in the Thermotogae phylum

Optimal growth temperature is a complex trait involving many cellular components, and its physiology is not yet fully understood. Evolution of continuous characters, such as optimal growth temperature, is often modeled as a one-dimensional random walk, but such a model may be an oversimplification given the complex processes underlying the evolution of continuous characters. Recent articles have used ancestral sequence reconstruction to infer the optimal growth temperature of ancient organisms from the guanine and cytosine content of the stem regions of ribosomal RNA, allowing inferences about the evolution of optimal growth temperature. Here, we investigate the optimal growth temperature of the bacterial phylum Thermotogae. Ancestral sequence reconstruction using a nonhomogeneous model was used to reconstruct the stem guanine and cytosine content of 16S rRNA sequences. We compare this sequence reconstruction method with other ancestral character reconstruction methods, and show that sequence reconstruction generates smaller confidence intervals and different ancestral values than other reconstruction methods. Unbiased random walk simulation indicates that the lower temperature members of the Thermotogales have been under directional selection; however, when a simulation is performed that takes possible mutations into account, it is the high temperature lineages that are, in fact, under directional selection. We find that the evolution of Thermotogales optimal growth temperatures is best fit by a biased random walk model. These findings suggest that it may be easier to evolve from a high optimal growth temperature to a lower one than vice versa.