New approaches for unravelling reassortment pathways

Quantifying the influence of genetic context on duplicated mammalian genes

Gene duplication is a fundamental part of evolutionary innovation. While single-gene duplications frequently exhibit asymmetric evolutionary rates between paralogs, the extent to which this applies to multi-gene duplications remains unclear. In this study, we investigate the role of genetic context in shaping evolutionary divergence within multi-gene duplications, leveraging microsynteny to differentiate source and target copies. Using a dataset of 193 mammalian genome assemblies and a bird outgroup, we systematically analyze patterns of sequence divergence between duplicated genes and reference orthologs. We find that target copies, those relocated to new genomic environments, exhibit elevated evolutionary rates compared to source copies in the ancestral location. This asymmetry is influenced by the distance between copies and the size of the target copy. We also demonstrate that the polarization of rate asymmetry in paralogs, the "choice" of the slowly evolving copy, is biased towards collective, block-wise polarization in multi-gene duplications. Our findings highlight the importance of genetic context in modulating post-duplication divergence, where differences in cis-regulatory elements and co-expressed gene clusters between source and target copies may be responsible. This study presents a large-scale test of asymmetric evolution in multi-gene duplications, offering new insight into how genome architecture shapes functional diversification of paralogs.

Recombination and retroprocessing in broomrapes reveal a universal roadmap for mitochondrial evolution in heterotrophic plants

The altered life history strategies of heterotrophic organisms often leave a profound genetic footprint on energy metabolism related functions. In parasitic plants, the reliance on host-derived nutrients and loss of photosynthesis in holoparasites have led to highly degraded to absent plastid genomes, but its impact on mitochondrial genome (mitogenome) evolution has remained controversial. By examining mitogenomes from 45 Orobanchaceae species including three independent transitions to holoparasitism and key evolutionary intermediates, we identified measurable and predictable genetic alterations in genomic shuffling, RNA editing, and intracellular (IGT) and horizontal gene transfer (HGT) en route to a nonphotosynthetic lifestyle. In-depth comparative analyses revealed DNA recombination and repair processes, especially RNA-mediated retroprocessing, as significant drivers for genome structure evolution. In particular, we identified a novel RNA-mediated IGT and HGT mechanism, which has not been demonstrated in cross-species and inter-organelle transfers. Based on this, we propose a generalized dosage effect mechanism to explain the biased transferability of plastid DNA to mitochondria across green plants, especially in heterotrophic lineages like parasites and mycoheterotrophs. Evolutionary rates scaled with these genomic changes, but the direction and strength of selection varied substantially among genes and clades, resulting in high contingency in mitochondrial genome evolution. Finally, we describe a universal roadmap for mitochondrial evolution in heterotrophic plants where increased recombination and repair activities, rather than relaxed selection alone, lead to differentiated genome structure compared to free-living species.

Concepts and Methods for Predicting Viral Evolution

The seasonal human influenza virus undergoes rapid evolution, leading to significant changes in circulating viral strains from year to year. These changes are typically driven by adaptive mutations, particularly in the antigenic epitopes, the regions of the viral surface protein hemagglutinin targeted by human antibodies. Here, we describe a consistent set of methods for data-driven predictive analysis of viral evolution. Our pipeline integrates four types of data: (1) sequence data of viral isolates collected on a worldwide scale, (2) epidemiological data on incidences, (3) antigenic characterization of circulating viruses, and (4) intrinsic viral phenotypes. From the combined analysis of these data, we obtain estimates of relative fitness for circulating strains and predictions of clade frequencies for periods of up to 1 year. Furthermore, we obtain comparative estimates of protection against future viral populations for candidate vaccine strains, providing a basis for pre-emptive vaccine strain selection. Continuously updated predictions obtained from the prediction pipeline for influenza and SARS-CoV-2 are available at https://previr.app .

Human lncRNAs harbor conserved modules embedded in different sequence contexts

We analyzed the structure of human long non-coding RNA (lncRNAs) genes to investigate whether the non-coding transcriptome is organized in modular domains, as is the case for protein-coding genes. To this aim, we compared all known human lncRNA exons and identified 340 pairs of exons with high sequence and/or secondary structure similarity but embedded in a dissimilar sequence context. We grouped these pairs in 106 clusters based on their reciprocal similarities. These shared modules are highly conserved between humans and the four great ape species, display evidence of purifying selection and likely arose as a result of recent segmental duplications. Our analysis contributes to the understanding of the mechanisms driving the evolution of the non-coding genome and suggests additional strategies towards deciphering the functional complexity of this class of molecules.

Concepts and methods for predicting viral evolution

The seasonal human influenza virus undergoes rapid evolution, leading to significant changes in circulating viral strains from year to year. These changes are typically driven by adaptive mutations, particularly in the antigenic epitopes, the regions of the viral surface protein haemagglutinin targeted by human antibodies. Here we describe a consistent set of methods for data-driven predictive analysis of viral evolution. Our pipeline integrates four types of data: (1) sequence data of viral isolates collected on a worldwide scale, (2) epidemiological data on incidences, (3) antigenic characterization of circulating viruses, and (4) intrinsic viral phenotypes. From the combined analysis of these data, we obtain estimates of relative fitness for circulating strains and predictions of clade frequencies for periods of up to one year. Furthermore, we obtain comparative estimates of protection against future viral populations for candidate vaccine strains, providing a basis for pre-emptive vaccine strain selection. Continuously updated predictions obtained from the prediction pipeline for influenza and SARS-CoV-2 are available on the website previr.app.

Concepts and methods for predicting viral evolution

The seasonal human influenza virus undergoes rapid evolution, leading to significant changes in circulating viral strains from year to year. These changes are typically driven by adaptive mutations, particularly in the antigenic epitopes, the regions of the viral surface protein haemagglutinin targeted by human antibodies. Here we describe a consistent set of methods for data-driven predictive analysis of viral evolution. Our pipeline integrates four types of data: (1) sequence data of viral isolates collected on a worldwide scale, (2) epidemiological data on incidences, (3) antigenic characterization of circulating viruses, and (4) intrinsic viral phenotypes. From the combined analysis of these data, we obtain estimates of relative fitness for circulating strains and predictions of clade frequencies for periods of up to one year. Furthermore, we obtain comparative estimates of protection against future viral populations for candidate vaccine strains, providing a basis for pre-emptive vaccine strain selection. Continuously updated predictions obtained from the prediction pipeline for influenza and SARS-CoV-2 are available on the website previr.app.

Phylogenetic Insights into H7Nx Influenza Viruses: Uncovering Reassortment Patterns and Geographic Variability

Influenza A viruses (IAVs), which belong to the Orthomyxoviridae family, are RNA viruses characterized by a segmented genome that allows them to evolve and adapt rapidly. These viruses are mainly transmitted by wild waterfowl. In this study, we investigated the evolutionary processes of H7Nx (H7N1, H7N2, H7N3, H7N4, H7N5, H7N6, H7N7, H7N8, H7N9) viruses, which pose a significant pandemic risk due to the known cases of human infection and their potential for rapid genetic evolution and reassortment. The complete genome sequences of H7Nx influenza viruses (n = 3239) were compared between each other to investigate their phylogenetic relationships and reassortment patterns. For the selected viruses, phylogenetic trees were constructed for eight genome segments (PB2, PB1, PA, HA, NP, NA, M, NS) to assess the genetic diversity and geographic distribution of these viruses. Distinct phylogenetic clades with remarkable geographic patterns were found for the different segments. While the viruses were consistently grouped by subtype based on the NA segment sequences, the phylogeny of the other segment sequences, with the exception of the NS segment, showed distinct grouping patterns based on geographic origin rather than formal subtype assignment. Reassortment events leading to complex phylogenetic relationships were frequently observed. In addition, multiple cases of previously undescribed reassortments between subtypes were detected, emphasizing the fluidity of H7Nx virus populations. These results indicate a high degree of genetic diversity and reassortment within H7Nx influenza viruses. In other words, H7Nx viruses exist as constantly changing combinations of gene pools rather than stable genetic lineages.

Snake species assemblages across habitat types in four departments of the Republic of Congo, with emphasis on medically-relevant venomous species

Snakebite is a neglected public health crisis in sub-Saharan Africa. There is a particular lack of data (on snakes and snakebite) from the central African region. This study was conducted in the departments of Likouala, Sangha, Cuvette-Ouest and Kouilou, in the Republic of Congo. The objective was to inventory snakes in the four localities with particular emphasis on medically relevant venomous snakes in order to improve knowledge of snakes in order to minimize the risks of snakebite envenomation to humans. Two methods (active and passive) were used to collect specimens from different habitats. Fifty-one (51)snake specimens including 14 medically relevant snake specimens representing 3 families, 3 subfamilies, 5 genera, and 6 species, in addition to 22 harmless species, were collected. We found a high number of medically important venomous species in Okoyo, Mokéko and Mvouti districts with 3 species each. The highest number of medically important venomous species was recorded in natural forests and human habitations, five and four species respectively. The species obtained (Atheris squamigera, Bitis arietans, Bitis gabonica, Dendroaspis jamesoni, Naja melanoleuca and Toxicodryas blandingii) are medically relevant toxic species according to the WHO classification. Further studies would be necessary to assess the epidemiology of bite risks snakes and educate the public to minimize accidental human-snake contact.

The Number and Pattern of Viral Genomic Reassortments are not Necessarily Identifiable from Segment Trees

Reassortment is an evolutionary process common in viruses with segmented genomes. These viruses can swap whole genomic segments during cellular co-infection, giving rise to novel progeny formed from the mixture of parental segments. Since large-scale genome rearrangements have the potential to generate new phenotypes, reassortment is important to both evolutionary biology and public health research. However, statistical inference of the pattern of reassortment events from phylogenetic data is exceptionally difficult, potentially involving inference of general graphs in which individual segment trees are embedded. In this paper, we argue that, in general, the number and pattern of reassortment events are not identifiable from segment trees alone, even with theoretically ideal data. We call this fact the fundamental problem of reassortment, which we illustrate using the concept of the "first-infection tree," a potentially counterfactual genealogy that would have been observed in the segment trees had no reassortment occurred. Further, we illustrate four additional problems that can arise logically in the inference of reassortment events and show, using simulated data, that these problems are not rare and can potentially distort our observation of reassortment even in small data sets. Finally, we discuss how existing methods can be augmented or adapted to account for not only the fundamental problem of reassortment, but also the four additional situations that can complicate the inference of reassortment.

IAVCP (Influenza A Virus Consensus and Phylogeny): Automatic Identification of the Genomic Sequence of the Influenza A Virus from High-Throughput Sequencing Data

Recently, high-throughput sequencing of influenza A viruses has become a routine test. It should be noted that the extremely high diversity of the influenza A virus complicates the task of determining the sequences of all eight genome segments. For a fast and accurate analysis, it is necessary to select the most suitable reference for each segment. At the same time, there is no standardized method in the field of decoding sequencing results that allows the user to update the sequence databases to which the reads obtained by virus sequencing are compared. The IAVCP (influenza A virus consensus and phylogeny) was developed with the goal of automatically analyzing high-throughput sequencing data of influenza A viruses. Its goals include the extraction of a consensus genome directly from paired raw reads. In addition, the pipeline enables the identification of potential reassortment events in the evolutionary history of the virus of interest by analyzing the topological structure of phylogenetic trees that are automatically reconstructed.

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.

Modeling the evolution of segment trees reveals deficiencies in current inferential methods for genomic reassortment

Reassortment is an evolutionary process common in viruses with segmented genomes. These viruses can swap whole genomic segments during cellular co-infection, giving rise to new viral variants. Large-scale genome rearrangements, such as reassortment, have the potential to quickly generate new phenotypes, making the understanding of viral reassortment important to both evolutionary biology and public health research. In this paper, we argue that reassortment cannot be reliably inferred from incongruities between segment phylogenies using the established remove-and-rejoin or coalescent approaches. We instead show that reassortment must be considered in the context of a broader population process that includes the dynamics of the infected hosts. Using illustrative examples and simulation we identify four types of evolutionary events that are difficult or impossible to reconstruct with incongruence-based methods. Further, we show that these specific situations are very common and will likely occur even in small samples. Finally, we argue that existing methods can be augmented or modified to account for all the problematic situations that we identify in this paper. Robust assessment of the role of reassortment in viral evolution is difficult, and we hope to provide conceptual clarity on some important methodological issues that can arise in the development of the next generation of tools for studying reassortment.

Phoenix dactylifera (date palm; Arecaceae) putative lectin homologs: Genome-wide search, architecture analysis, and evolutionary relationship

The date palm, Phoenix dactylifera, is a vital crop in nations in the Middle East and North Africa. The date palm was thought to have outstanding traditional medicinal value because it was abundant in phytochemicals with diverse chemical structures. The date palm's ability to withstand harsh environments could be partly attributed to a class of proteins known as lectins, which are carbohydrate-binding proteins that can bind sugar moieties reversibly and without changing their chemical structures. After scanning the genome of P. dactylifera (GCF 009389715.1), this in silico study discovered 196 possible lectin homologs from 11 different families, some specific to plants. At the same time, others could also be found in other kingdoms of life. Their domain architectures and functional amino acid residues were investigated, and they yielded a 40% true-lectin with known conserved carbohydrate-binding residues. Further, their probable subcellular localization, physiochemical and phylogenetic analyses were also performed. Scanning all putative lectin homologs against the anticancer peptide (ACP) dataset found in the AntiCP2.0 webpage identified 26 genes with protein kinase receptors (Lec-KRs) belonging to 5 lectin families, which are reported to have at least one ACP motif. Our study offers the first account of Phoenix-lectins and their organization that can be used for further structural and functional analysis and investigating their potential as anticancer proteins.

Cross-Hemispheric Genetic Diversity and Spatial Genetic Structure of Callinectes sapidus Reovirus 1 (CsRV1)

The movement of viruses in aquatic systems is rarely studied over large geographic scales. Oceanic currents, host migration, latitude-based variation in climate, and resulting changes in host life history are all potential drivers of virus connectivity, adaptation, and genetic structure. To expand our understanding of the genetic diversity of Callinectes sapidus reovirus 1 (CsRV1) across a broad spatial and host life history range of its blue crab host (Callinectes sapidus), we obtained 22 complete and 96 partial genomic sequences for CsRV1 strains from the US Atlantic coast, Gulf of Mexico, Caribbean Sea, and the Atlantic coast of South America. Phylogenetic analyses of CsRV1 genomes revealed that virus genotypes were divided into four major genogroups consistent with their host geographic origins. However, some CsRV1 sequences from the US mid-Atlantic shared high genetic similarity with the Gulf of Mexico genotypes, suggesting potential human-mediated movement of CsRV1 between the US mid-Atlantic and Gulf coasts. This study advances our understanding of how climate, coastal geography, host life history, and human activity drive patterns of genetic structure and diversity of viruses in marine animals and contributes to the capacity to infer broadscale host population connectivity in marine ecosystems from virus population genetic data.

Characterization of the complete plastid genome of Korean endemic, Ajuga spectabilis Nakai (Lamiaceae)

Ajuga spectabilis Nakai is a Korean endemic species in Lamiaceae. In spite of its importance, genomic studies are not performed on this species. Here, we report the complete plastid genome sequences of A. spectabilis, which will provide valuable information for its natural conservation and future studies for the plastid genome evolution. The plastid genome is 150,417 bp in length, containing a large single-copy region (LSC) of 82,140 bp and a small single-copy (SSC) region of 17,165 bp which are separated by a pair of inverted repeats (IR) of 25,556 bp. It encodes 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The overall GC content is 38.3%, and those in the LSC, SSC, and IR regions are 36.4%, 32.2%, and 43.3%, respectively, which is consistent with other Ajuga species. Our phylogenetic analysis revealed that A. spectabilis formed a close relationship with A. ciliata and A. decumbens.

Plastome characteristics and species identification of Chinese medicinal wintergreens (Gaultheria, Ericaceae)

Wintergreen oil is a folk medicine widely used in foods, pesticides, cosmetics and drugs. In China, nine out of 47 species within Gaultheria (Ericaceae) are traditionally used as Chinese medicinal wintergreens; however, phylogenetic approaches currently used to discriminating these species remain unsatisfactory. In this study, we sequenced and characterized plastomes from nine Chinese wintergreen species and identified candidate DNA barcoding regions for Gaultheria. Each Gaultheria plastome contained 110 unique genes (76 protein-coding, 30 tRNA, and four rRNA genes). Duplication of trnfM, rps14, and rpl23 genes were detected, while all plastomes lacked ycf1 and ycf2 genes. Gaultheria plastomes shared substantially contracted SSC regions that contained only the ndhF gene. Moreover, plastomes of Gaultheria leucocarpa var. yunnanensis contained an inversion in the LSC region and an IR expansion to cover the ndhF gene. Multiple rearrangement events apparently occurred between the Gaultheria plastomes and those from several previously reported families in Ericales. Our phylogenetic reconstruction using 42 plastomes revealed well-supported relationships within all nine Gaultheria species. Additionally, seven mutational hotspot regions were identified as potential DNA barcodes for Chinese medicinal wintergreens. Our study is the first to generate complete plastomes and describe the structural variations of the complicated genus Gaultheria. In addition, our findings provide important resources for identification of Chinese medicinal wintergreens.

Cranial ecomorphology of turtles and neck retraction as a possible trigger of ecological diversification

Turtles have a highly modified body plan, including a rigid shell that constrains postcranial anatomy. Skull morphology and neck mobility may therefore be key to ecological specialization in turtles. However, the ecological signal of turtle skull morphologies has not been rigorously evaluated, leaving uncertainties about the roles of ecological adaptation and convergence. We evaluate turtle cranial ecomorphology using three-dimensional geometric morphometrics and phylogenetic comparative methods. Skull shape correlates with allometry, neck retraction capability, and different aquatic feeding ecologies. We find that ecological variables influence skull shape only, whereas a key functional variable (the capacity for neck retraction) influences both shape and size. Ecology and functional predictions from three-dimensional shape are validated by high success rates for extant species, outperforming previous two-dimensional approaches. We use this to infer ecological and functional traits of extinct species. Neck retraction evolved among crownward stem-turtles by the Late Jurassic, signaling functional decoupling of the skull and neck from the shell, possibly linked to a major episode of ecomorphological diversification. We also find strong evidence for convergent ecological adaptations among marine groups. This includes parallel loss of neck retraction, evidence for active hunting, possible grazing, and suction feeding in extinct marine groups. Our large-scale assessment of dietary and functional adaptation throughout turtle evolution reveals the timing and origin of their distinct ecomorphologies, and highlights the potential for ecology and function to have distinct effects on skull form.

TreeKnit: Inferring ancestral reassortment graphs of influenza viruses

When two influenza viruses co-infect the same cell, they can exchange genome segments in a process known as reassortment. Reassortment is an important source of genetic diversity and is known to have been involved in the emergence of most pandemic influenza strains. However, because of the difficulty in identifying reassortment events from viral sequence data, little is known about their role in the evolution of the seasonal influenza viruses. Here we introduce TreeKnit, a method that infers ancestral reassortment graphs (ARG) from two segment trees. It is based on topological differences between trees, and proceeds in a greedy fashion by finding regions that are compatible in the two trees. Using simulated genealogies with reassortments, we show that TreeKnit performs well in a wide range of settings and that it is as accurate as a more principled bayesian method, while being orders of magnitude faster. Finally, we show that it is possible to use the inferred ARG to better resolve segment trees and to construct more informative visualizations of reassortments.

Influenza viruses evolve quickly and escape immune defenses which requires frequent update of vaccines. Understanding this evolution is key to an effective public health response. The genome of influenza viruses is made up of 8 pieces called segments, each coding for different viral proteins. Within each segment, evolution is an asexual process in which genetic diversity is generated by mutations. But influenza also diversifies through reassortment which can occur when two different viruses infect the same cell: offsprings can then contain a combination of segments from both viruses. Reassortment is akin to sexual reproduction and can generate viruses that combine segments from diverged viral lineages. Reassortment is a crucial component of viral evolution, but it is challenging to reconstruct where reassortments happened and which segments share history. Here, we develop a method called TreeKnit to detect reassortment events. TreeKnit is based on genealogical trees of single segments that can be reconstructed using standard bioinformatics tools. Inconsistencies between these trees are then used as signs of reassortment. We show that TreeKnit is as accurate as other recent methods, but runs much faster. Our method will facilitate the study of reassortment and its consequences for influenza evolution.

A Bayesian approach to infer recombination patterns in coronaviruses

As shown during the SARS-CoV-2 pandemic, phylogenetic and phylodynamic methods are essential tools to study the spread and evolution of pathogens. One of the central assumptions of these methods is that the shared history of pathogens isolated from different hosts can be described by a branching phylogenetic tree. Recombination breaks this assumption. This makes it problematic to apply phylogenetic methods to study recombining pathogens, including, for example, coronaviruses. Here, we introduce a Markov chain Monte Carlo approach that allows inference of recombination networks from genetic sequence data under a template switching model of recombination. Using this method, we first show that recombination is extremely common in the evolutionary history of SARS-like coronaviruses. We then show how recombination rates across the genome of the human seasonal coronaviruses 229E, OC43 and NL63 vary with rates of adaptation. This suggests that recombination could be beneficial to fitness of human seasonal coronaviruses. Additionally, this work sets the stage for Bayesian phylogenetic tracking of the spread and evolution of SARS-CoV-2 in the future, even as recombinant viruses become prevalent.

The Jurassic rise of squamates as supported by lepidosaur disparity and evolutionary rates

The squamates (lizards, snakes, and relatives) today comprise more than 10,000 species, and yet their sister group, the Rhynchocephalia, is represented by a single species today, the tuatara. The explosion in squamate diversity has been tracked back to the Cretaceous Terrestrial Revolution, 100 million years ago (Ma), the time when flowering plants began their takeover of terrestrial ecosystems, associated with diversification of coevolving insects and insect-eating predators such as lizards, birds, and mammals. Squamates arose much earlier, but their long pre-Cretaceous history of some 150 million years (Myr) is documented by sparse fossils. Here, we provide evidence for an initial radiation of squamate morphology in the Middle and Late Jurassic (174-145 Ma), and show that they established their key ecological roles much earlier than had been assumed, and they have not changed them much since.

Plastid RNA editing reduction accompanied with genetic variations in Cymbidium, a genus with diverse lifestyle modes

Recent sequencing efforts have broadly uncovered the evolutionary trajectory of plastid genomes (plastomes) of flowering plants in diverse habitats, yet our knowledge of the evolution of plastid posttranscriptional modifications is limited. In this study, we generated 11 complete plastomes and performed ultra-deep transcriptome sequencing to investigate the co-evolution of plastid RNA editing and genetic variation in Cymbidium, a genus with diverse trophic lifestyles. Genome size and gene content is reduced in terrestrial and green mycoheterotrophic orchids relative to their epiphytic relatives. This could be partly due to extensive losses and pseudogenization of ndh genes for the plastid NADH dehydrogenase-like complex, but independent pseudogenization of ndh genes has also occurred in the epiphyte C. mannii, which was reported to use strong crassulacean acid metabolism photosynthesis. RNA editing sites are abundant but variable in number among Cymbidium plastomes. The nearly twofold variation in editing abundance is mainly due to extensive reduction of ancestral editing sites in ndh transcripts of terrestrial, mycoheterotrophic, and C. mannii plastomes. The co-occurrence of editing reduction and pseudogenization in ndh genes suggests functional constraints on editing machinery may be relaxed, leading to nonrandom loss of ancestral edited sites via reduced editing efficiency. This study represents the first systematic examination of RNA editing evolution linked to plastid genome variation in a single genus. We also propose an explanation for how genomic and posttranscriptional variations might be affected by lifestyle-associated ecological adaptation strategies in Cymbidium.

Phylogenomic and comparative analyses of Coffeeae alliance (Rubiaceae): deep insights into phylogenetic relationships and plastome evolution

BACKGROUND

The large and diverse Coffeeae alliance clade of subfamily Ixoroideae (Rubiaceae) consists of 10 tribes, > 90 genera, and > 2000 species. Previous molecular phylogenetics using limited numbers of markers were often unable to fully resolve the phylogenetic relationships at tribal and generic levels. Also, the structural variations of plastomes (PSVs) within the Coffeeae alliance tribes have been poorly investigated in previous studies. To fully understand the phylogenetic relationships and PSVs within the clade, highly reliable and sufficient sampling with superior next-generation analysis techniques is required. In this study, 71 plastomes (40 newly sequenced and assembled and the rest from the GenBank) were comparatively analyzed to decipher the PSVs and resolve the phylogenetic relationships of the Coffeeae alliance using four molecular data matrices.

RESULTS

All plastomes are typically quadripartite with the size ranging from 153,055 to 155,908 bp and contained 111 unique genes. The inverted repeat (IR) regions experienced multiple contraction and expansion; five repeat types were detected but the most abundant was SSR. The size of the Coffeeae alliance clade plastomes and its elements are affected by the IR boundary shifts and the repeat types. However, the emerging PSVs had no taxonomic and phylogenetic implications. Eight highly divergent regions were identified within the plastome regions ndhF, ccsA, ndhD, ndhA, ndhH, ycf1, rps16-trnQ-UUG, and psbM-trnD. These highly variable regions may be potential molecular markers for further species delimitation and population genetic analyses for the clade. Our plastome phylogenomic analyses yielded a well-resolved phylogeny tree with well-support at the tribal and generic levels within the Coffeeae alliance.

CONCLUSIONS

Plastome data could be indispensable in resolving the phylogenetic relationships of the Coffeeae alliance tribes. Therefore, this study provides deep insights into the PSVs and phylogenetic relationships of the Coffeeae alliance and the Rubiaceae family as a whole.

Intracellular translocation of HMGB1 is important for Zika virus replication in Huh7 cells

Neonatal microcephaly and adult Guillain-Barré syndrome are severe complications of Zika virus (ZIKV) infection. The robustly induced inflammatory cytokine expressions in ZIKV-infected patients may constitute a hallmark for severe disease. In the present study, the potential role of high mobility group box 1 protein (HMGB1) in ZIKV infection was investigated. HMGB1 protein expression was determined by the enzyme-linked immunosorbent assay (ELISA) and immunoblot assay. HMGB1's role in ZIKV infection was also explored using treatment with dexamethasone, an immunomodulatory drug, and HMGB1-knockdown (shHMGB1) Huh7 cells. Results showed that the Huh7 cells were highly susceptible to ZIKV infection. The infection was found to induce HMGB1 nuclear-to-cytoplasmic translocation, resulting in a > 99% increase in the cytosolic HMGB1 expression at 72-h post-infection (h.p.i). The extracellular HMGB1 level was elevated in a time- and multiplicity of infection (MOI)-dependent manner. Treatment of the ZIKV-infected cells with dexamethasone (150 µM) reduced HMGB1 extracellular release in a dose-dependent manner, with a maximum reduction of 71 ± 5.84% (P < 0.01). The treatment also reduced virus titers by over 83 ± 0.50% (P < 0.01). The antiviral effects, however, were not observed in the dexamethasone-treated shHMGB1 cells. These results suggest that translocation of HMGB1 occurred during ZIKV infection and inhibition of the translocation by dexamethasone coincided with a reduction in ZIKV replication. These findings highlight the potential of targeting the localization of HMGB1 in affecting ZIKV infection.

Phenotypic differentiation of the slow worm lizards (Squamata: Anguis) across their contact zone in Central Europe

Background

The application of molecular-phylogenetic approaches to taxonomy has had a dramatic effect on our understanding of the diversity of reptiles. These approaches have allowed researchers to reveal previously hidden lineages as well as taxonomic overestimation in morphologically plastic taxa. Slow worms, legless lizards of the genus Anguis (Squamata: Anguidae), were previously considered to comprise either one or two species, and morphology-based intraspecific taxonomy of Anguis fragilis remained controversial throughout the 20th century. After the discovery of deep genetic divergences within the genus, its taxonomy was reconsidered, and as a result, five extant species have been recognized. In order to better understand the patterns of their interspecific differentiation, here we studied phenotypic differences between the two most widespread of them—A. fragilis and A. colchica, and their putative hybrids across the contact zone of both species in Central Europe.

Methods

We used multivariate and univariate statistics and analyzed ten metric, eleven meristic, and six categorical phenotypic variables in material comprising a total of 326 individuals. We also genotyped individuals from the contact zone for one mitochondrial and two nuclear DNA fragments in order to delineate the distribution of individuals of hybrid and non-hybrid origin. The clines in morphological traits were studied using HZAR.

Results

We show that the two species are morphologically differentiated. Anguis fragilis has a less robust head, fewer scales covering the body, lower frequency of the external ear opening presence, lower frequency of separated prefrontal scales, higher frequency of prefrontal scales in contact with each other, and body coloration more similar to the juvenile coloration than A. colchica. Slow worms from the contact/hybrid zone are characterized by an intermediate morphology, with more similarities to A. fragilis than to A. colchica.

Discussion

None of the analyzed characters alone proved to be fully diagnostic, although more than 90% of all individuals could be successfully assigned to one or another species based on numbers of scales around the body. Our results indicate concordant, coincident, and steep clines in character states change. We present several hypotheses on the origin and evolutionary maintenance of the morphological divergence between both species and suggest that different evolutionary histories of the taxa rather than recently acting selection explain the observed morphological variation.

High diversity of small insectivorous mammals on Qinghai-Tibet Plateau and first description of karyotype for four endemics of China

Among seven species of the order Eulipotyphla (from southern Gansu and northern Sichuan Provinces, Central China) studied cytogenetically, karyotypes of one talpid species, Uropsilus aff. soricipes (2n = 36, NFa = 54), and three soricid species, Chodsigoa hypsibia (2n = 65, NFa = 66), Sorex cansulus (2n = 42, NFa = 64) and Sorex thibetanus (2n = 42, NFa = 60), are described cytogenetically for the first time. All four species are endemic to China with distribution ranges restricted to the Qinghai–Tibet Plateau and adjacent mountain ranges. The Ch. hypsibia karyotype consists of mostly acrocentric autosomes and one metacentric pair of autosomes; besides, a B chromosome was identified. No polymorphism was detected among karyotypes of other species, including shrews Sorex bedfordiae (2n = 26, NFa = 44), Anourosorex squamipes (2n = 48, NFa = 92) and Crocidura suaveolens (2n = 40, NFa = 44). The Chinese shrew mole U. aff. soricipes and three shrew species (S. bedfordiae, Ch. hypsibia and A. squamipes) represent autochthonous fauna of Central/Western China, whereas S. thibetanus, S. cansulus and C. suaveolens belong to phylogenetic groups occurring mostly to the north and west from China; therefore, they should be considered relatively recent colonisers. Thus, considering the relationships of the species within phylogenetic groups, our results on karyotypes are in good agreement with molecular genetic data.

New insights into the evolution of SPX gene family from algae to legumes; a focus on soybean

BACKGROUND

SPX-containing proteins have been known as key players in phosphate signaling and homeostasis. In Arabidopsis and rice, functions of some SPXs have been characterized, but little is known about their function in other plants, especially in the legumes.

RESULTS

We analyzed SPX gene family evolution in legumes and in a number of key species from algae to angiosperms. We found that SPX harboring proteins showed fluctuations in domain fusions from algae to the angiosperms with, finally, four classes appearing and being retained in the land plants. Despite these fluctuations, Lysine Surface Cluster (KSC), and the third residue of Phosphate Binding Sites (PBS) showed complete conservation in almost all of SPXs except few proteins in Selaginella moellendorffii and Papaver sumniferum, suggesting they might have different ligand preferences. In addition, we found that the WGD/segmentally or dispersed duplication types were the most frequent contributors to the SPX expansion, and that there is a positive correlation between the amount of WGD contribution to the SPX expansion in individual species and its number of EXS genes. We could also reveal that except SPX class genes, other classes lost the collinearity relationships among Arabidopsis and legume genomes. The sub- or neo-functionalization of the duplicated genes in the legumes makes it difficult to find the functional orthologous genes. Therefore, we used two different methods to identify functional orthologs in soybean and Medicago. High variance in the dynamic and spatial expression pattern of GmSPXs proved the new or sub-functionalization in the paralogs.

CONCLUSION

This comprehensive analysis revealed how SPX gene family evolved from algae to legumes and also discovered several new domains fused to SPX domain in algae. In addition, we hypothesized that there different phosphate sensing mechanisms might occur in S. moellendorffii and P. sumniferum. Finally, we predicted putative functional orthologs of AtSPXs in the legumes, especially, orthologs of AtPHO1, involved in long-distance Pi transportation. These findings help to understand evolution of phosphate signaling and might underpin development of new legume varieties with improved phosphate use efficiency.

Reassortment Network of Influenza A Virus

Influenza A virus (IAV) genomes are composed of eight single-stranded RNA segments. Genetic exchange through reassortment of the segmented genomes often endows IAVs with new genetic characteristics, which may affect transmissibility and pathogenicity of the viruses. However, a comprehensive understanding of the reassortment history of IAVs remains lacking. To this end, we assembled 40,296 whole-genome sequences of IAVs for analysis. Using a new clustering method based on Mean Pairwise Distances in the phylogenetic trees, we classified each segment of IAVs into clades. Correspondingly, reassortment events among IAVs were detected by checking the segment clade compositions of related genomes under specific environment factors and time period. We systematically identified 1,927 possible reassortment events of IAVs and constructed their reassortment network. Interestingly, minimum spanning tree of the reassortment network reproved that swine act as an intermediate host in the reassortment history of IAVs between avian species and humans. Moreover, reassortment patterns among related subtypes constructed in this study are consistent with previous studies. Taken together, our genome-wide reassortment analysis of all the IAVs offers an overview of the leaping evolution of the virus and a comprehensive network representing the relationships of IAVs.

SARS-CoV-2 transmission dynamics in Belarus in 2020 revealed by genomic and incidence data analysis

Background

Non-pharmaceutical interventions (NPIs) have been implemented worldwide to curb COVID-19 spread. Belarus is a rare case of a country with a relatively modern healthcare system, where highly limited NPIs have been enacted. Thus, investigation of Belarusian COVID-19 dynamics is essential for the local and global assessment of the impact of NPI strategies.

Methods

We integrate genomic epidemiology and surveillance methods to investigate the spread of SARS-CoV-2 in Belarus in 2020. We utilize phylodynamics, phylogeography, and probabilistic bias inference to study the virus import and export routes, the dynamics of the effective reproduction number, and the incidence of SARS-CoV-2 infection.

Results

Here we show that the estimated cumulative number of infections by June 2020 exceeds the confirmed case number by a factor of ~4 (95% confidence interval (2; 9)). Intra-country SARS-CoV-2 genomic diversity originates from at least 18 introductions from different regions, with a high proportion of regional transmissions. Phylodynamic analysis indicates a moderate reduction of the effective reproductive number after the introduction of limited NPIs, but its magnitude is lower than for developed countries with large-scale NPIs. On the other hand, the effective reproduction number estimate is comparable with that for the neighboring Ukraine, where NPIs were broader.

Conclusions

The example of Belarus demonstrates how countries with relatively low outward population mobility continue to be integral parts of the global epidemiological environment. Comparison of the effective reproduction number dynamics for Belarus and other countries reveals the effect of different NPI strategies but also emphasizes the role of regional Eastern European sociodemographic factors in the virus spread.

Belarus is one of few European countries that has enacted limited measures to contain SARS-CoV-2, the virus that causes COVID-19. We study the genetic sequences of the SARS-CoV-2 virus circulating in Belarus and other countries in 2020 to investigate how it might have been imported into the country and spread there. We show that the virus was repeatedly imported from and exported to different regions, including a large portion of regional transmissions that occurred despite stricter measures implemented by Belarus’ neighbors. There was a moderate reduction of the virus reproductive number—a measure of virus transmission speed—after April 2020, but its magnitude was lower than for developed countries with more stringent epidemiological interventions. These findings shed light on the COVID-19 spread in Eastern Europe and highlight the impact of public health policies and of regional factors on this spread.

Nemira et al. study the genomic epidemiology and phylodynamics of SARS-CoV-2 in Belarus. They identify potential introduction routes of the virus from other countries, determine that during the first wave of the pandemic the number of infections was likely several times higher than reported case numbers, and estimate the impact of early non-pharmaceutical interventions on SARS-CoV-2 transmission.

New interpretation of the cranial osteology of the Early Cretaceous turtle Arundelemys dardeni (Paracryptodira) based on a CT-based re-evaluation of the holotype

Arundelemys dardeni is an Early Cretaceous paracryptodire known from a single, incomplete, but generally well-preserved skull. Phylogenetic hypotheses of paracryptodires often find Arundelemys dardeni as an early branching baenid. As such, it has a central role in understanding the early evolution of the successful clade Baenidae, which survived the Cretaceous-Paleogene mass extinction, as well as the diversification of Paracryptodira into its subclades, which recent research suggests to perhaps include helochelydrids, compsemydids, pleurosternids, and baenids. Computer tomography scans of the holotype material that were produced for the initial description of Arundelemeys dardeni reveal several errors in the initial anatomical description of the species, which we correct based on element-by-element segmentation. In addition, we provide entirely novel anatomical information, including descriptions of several previously undescribed cranial bones, the endosseous labyrinth, and the cranial scutes, the latter of which are unknown for most paracryptodires. We provide an interpretation of cranial scutes which homologizes the scutes of Arundelemys dardeni with those of other stem turtles.

Progress and Challenge in Computational Identification of Influenza Virus Reassortment

Genomic reassortment is an important evolutionary mechanism for influenza viruses. In this process, the novel viruses acquire new characteristics by the exchange of the intact gene segments among multiple influenza virus genomes, which may cause flu endemics and epidemics within or even across hosts. Due to the safety and ethical limitations of the experimental studies on influenza virus reassortment, numerous computational researches on the influenza virus reassortment have been done with the explosion of the influenza virus genomic data. A great amount of computational methods and bioinformatics databases were developed to facilitate the identification of influenza virus reassortments. In this review, we summarized the progress and challenge of the bioinformatics research on influenza virus reassortment, which can guide the researchers to investigate the influenza virus reassortment events reasonably and provide valuable insight to develop the related computational identification tools.

Plastome phylogenomics of the East Asian endemic genus Dobinea

Dobinea is a dioecious genus endemic to East Asia that consists of two extant species: Dobinea delavayi and Dobinea vulgaris. Although the genus is morphologically distinct, its phylogenetic position remains controversial. In this study, we investigated the phylogenetic relationships between Dobinea and related taxa by sequencing the whole plastome DNA sequences for both extant species of Dobinea and comparing them to published plastomes within Sapindales. The complete plastomes of D. vulgaris and D. delavayi were 160,683 and 160, 154 base pairs (bp) in length, including a pair of inverted repeat regions (IRs, 26,889 and 26,759 bp) divided by the large single-copy region (LSC, 87,962 and 87,555 bp) and small single-copy region (SSC, 18,943 and 19,081 bp), and identically encoded 113 unique genes (79 protein-coding genes, 30 tRNAs, and 4 rRNA genes). Plastid phylogenomic analyses showed that Dobinea was a well-supported monophyletic unit and sister to the clade including tribes Anacardieae and Rhoideae, which suggests that Dobinea is a member of Anacardiaceae. In addition, molecular dating inferred D. delavayi and D. vulgaris diverged approximately 10.76 Ma, suggesting the divergence between these two species may have been driven by the intensification of the Asian summer monsoon and the establishment of distinct monsoon regimes in East Asia.

Genomes reveal genetic diversity of Piscine orthoreovirus in farmed and free-ranging salmonids from Canada and USA

Piscine orthoreovirus (PRV-1) is a segmented RNA virus, which is commonly found in salmonids in the Atlantic and Pacific Oceans. PRV-1 causes the heart and skeletal muscle inflammation disease in Atlantic salmon and is associated with several other disease conditions. Previous phylogenetic studies of genome segment 1 (S1) identified four main genogroups of PRV-1 (S1 genogroups I–IV). The goal of the present study was to use Bayesian phylogenetic inference to expand our understanding of the spatial, temporal, and host patterns of PRV-1 from the waters of the northeast Pacific. To that end, we determined the coding genome sequences of fourteen PRV-1 samples that were selected to improve our knowledge of genetic diversity across a broader temporal, geographic, and host range, including the first reported genome sequences from the northwest Atlantic (Eastern Canada). Nucleotide and amino acid sequences of the concatenated genomes and their individual segments revealed that established sequences from the northeast Pacific were monophyletic in all analyses. Bayesian inference phylogenetic trees of S1 sequences using BEAST and MrBayes also found that sequences from the northeast Pacific grouped separately from sequences from other areas. One PRV-1 sample (WCAN_BC17_AS_2017) from an escaped Atlantic salmon, collected in British Columbia but derived from Icelandic broodstock, grouped with other S1 sequences from Iceland. Our concatenated genome and S1 analysis demonstrated that PRV-1 from the northeast Pacific is genetically distinct but descended from PRV-1 from the North Atlantic. However, the analyses were inconclusive as to the timing and exact source of introduction into the northeast Pacific, either from eastern North America or from European waters of the North Atlantic. There was no evidence that PRV-1 was evolving differently between free-ranging Pacific Salmon and farmed Atlantic Salmon. The northeast Pacific PRV-1 sequences fall within genogroup II based on the classification of Garseth, Ekrem, and Biering (Garseth, A. H., Ekrem, T., and Biering, E. (2013) ‘Phylogenetic Evidence of Long Distance Dispersal and Transmission of Piscine Reovirus (PRV) between Farmed and Wild Atlantic Salmon’, PLoS One, 8: e82202.), which also includes North Atlantic sequences from Eastern Canada, Iceland, and Norway. The additional full-genome sequences herein strengthen our understanding of phylogeographical patterns related to the northeast Pacific, but a more balanced representation of full PRV-1 genomes from across its range, as well additional sequencing of archived samples, is still needed to better understand global relationships including potential transmission links among regions.

Cranial osteology of the Early Cretaceous turtle Pleurosternon bullockii (Paracryptodira: Pleurosternidae)

Pleurosternon bullockii is a turtle from the Early Cretaceous of Europe known from numerous postcranial remains. Only one skull has so far been referred to the species. Pleurosternon bullockii belongs to a group of turtles called pleurosternids, which is thought to include several poorly known taxa from the Late Jurassic and Early Cretaceous of Europe and North America. Pleurosternids and baenids, a group of North American turtles that lived from the Late Cretaceous to the Eocene, define a clade called Paracryptodira. Additionally, Paracryptodira likely includes compsemydids, and, potentially, helochelydrids. Character support for Paracryptodira is relatively weak, and many global phylogenetic studies fail to support paracryptodiran monophyly altogether. Proposed paracryptodiran synapomorphies are largely cranial, despite the poor characterization of pleurosternid cranial material. In addition to their questionable monophyly, the global position of paracryptodires is debated. Early studies suggest crown-turtle affinities, but most phylogenies find them as stem-turtles, irrespective of their monophyly. Here, we document the cranial osteology of Pleurosternon bullockii with the use of three-dimensional models derived from segmenting high-resolution X-ray micro-computed tomography (CT) scans. Pleurosternon bullockii has a primitive basipterygoid region of the skull, but a cryptodire-like acustico-jugular region. A surprising number of similarities with pleurodires exist, particularly in the laterally expanded external process of the pterygoid and in the posterior orbital wall. Our observations constitute an important step toward a phylogenetic re-evaluation of Paracryptodira.

Stable species boundaries despite ten million years of hybridization in tropical eels

Genomic evidence is increasingly underpinning that hybridization between taxa is commonplace, challenging our views on the mechanisms that maintain their boundaries. Here, we focus on seven catadromous eel species (genus Anguilla) and use genome-wide sequence data from more than 450 individuals sampled across the tropical Indo-Pacific, morphological information, and three newly assembled draft genomes to compare contemporary patterns of hybridization with signatures of past introgression across a time-calibrated phylogeny. We show that the seven species have remained distinct for up to 10 million years and find that the current frequencies of hybridization across species pairs contrast with genomic signatures of past introgression. Based on near-complete asymmetry in the directionality of hybridization and decreasing frequencies of later-generation hybrids, we suggest cytonuclear incompatibilities, hybrid breakdown, and purifying selection as mechanisms that can support species cohesion even when hybridization has been pervasive throughout the evolutionary history of clades.

HopPER: an adaptive model for probability estimation of influenza reassortment through host prediction

Background

Influenza reassortment, a mechanism where influenza viruses exchange their RNA segments by co-infecting a single cell, has been implicated in several major pandemics since 19th century. Owing to the significant impact on public health and social stability, great attention has been received on the identification of influenza reassortment.

Methods

We proposed a novel computational method named HopPER (Host-prediction-based Probability Estimation of Reassortment), that sturdily estimates reassortment probabilities through host tropism prediction using 147 new features generated from seven physicochemical properties of amino acids. We conducted the experiments on a range of real and synthetic datasets and compared HopPER with several state-of-the-art methods.

Results

It is shown that 280 out of 318 candidate reassortants have been successfully identified. Additionally, not only can HopPER be applied to complete genomes but its effectiveness on incomplete genomes is also demonstrated. The analysis of evolutionary success of avian, human and swine viruses generated through reassortment across different years using HopPER further revealed the reassortment history of the influenza viruses.

Conclusions

Our study presents a novel method for the prediction of influenza reassortment. We hope this method could facilitate rapid reassortment detection and provide novel insights into the evolutionary patterns of influenza viruses.

The complete chloroplast genome of Microcycas calocoma (Miq.) A. DC. (Zamiaceae, Cycadales) and evolution in Cycadales

Cycadales is an extant group of seed plants occurring in subtropical and tropical regions comprising putatively three families and 10 genera. At least one complete plastid genome sequence has been reported for all of the 10 genera except Microcycas, making it an ideal plant group to conduct comprehensive plastome comparisons at the genus level. This article reports for the first time the plastid genome of Microcycas calocoma. The plastid genome has a length of 165,688 bp with 134 annotated genes including 86 protein-coding genes, 47 non-coding RNA genes (39 tRNA and eight rRNA) and one pseudogene. Using global sequence variation analysis, the results showed that all cycad genomes share highly similar genomic profiles indicating significant slow evolution and little variation. However, identity matrices coinciding with the inverted repeat regions showed fewer similarities indicating that higher polymorphic events occur at those sites. Conserved non-coding regions also appear to be more divergent whereas variations in the exons were less discernible indicating that the latter comprises more conserved sequences. Phylogenetic analysis using 81 concatenated protein-coding genes of chloroplast (cp) genomes, obtained using maximum likelihood and Bayesian inference with high support values (>70% ML and = 1.0 BPP), confirms that Microcycas is closest to Zamia and forms a monophyletic clade with Ceratozamia and Stangeria. While Stangeria joined the Neotropical cycads Ceratozamia, Zamia and Microcyas, Bowenia grouped with the Southern Hemisphere cycads Encephalartos, Lepidozamia and Macrozamia. All Cycas species formed a distinct clade separated from the other genera. Dioon, on the other hand, was outlying from the rest of Zamiaceae encompassing two major clades—the Southern Hemisphere cycads and the Neotropical cycads. Analysis of the whole cp genomes in phylogeny also supports that the previously recognized family—Stangeriaceae—which contained Bowenia and Stangeria, is not monophyletic. Thus, the cp genome topology obtained in our study is congruent with other molecular phylogenies recognizing only a two-family classification (Cycadaceae and Zamiaceae) within extant Cycadales.

An update to the inventory of shore-fishes from the Parque Nacional Sistema Arrecifal Veracruzano, Veracruz, México

Data on marine and brackish-water fishes recorded in the area of the Parque Nacional Sistema Arrecifal Veracruzano in the southwest Gulf of Mexico were extracted from online aggregators of georeferenced location records, the recent ichthyological literature reviewed, and collections and observations made to provide a more complete faunal inventory for that park. Those actions added 95 species to a comprehensive inventory published in 2013, and brought the total to 472 species, an increase of 22%. Seventy-four percent of the additions came from online aggregators of georeferenced species records, which clearly demonstrates the value of reviewing and incorporating such data into species inventories. However, different aggregators recorded different sets of species, and some of their data were linked to outdated taxonomy or included identification errors. Hence individual records from multiple aggregators need to be obtained and reviewed for such issues when using such data to compile and revise faunal inventories. Existing lists also need to be carefully reviewed to ensure that errors are not perpetuated during updates.

The evolution of parental care diversity in amphibians

Parental care is extremely diverse across species, ranging from simple behaviours to complex adaptations, varying in duration and in which sex cares. Surprisingly, we know little about how such diversity has evolved. Here, using phylogenetic comparative methods and data for over 1300 amphibian species, we show that egg attendance, arguably one of the simplest care behaviours, is gained and lost faster than any other care form, while complex adaptations, like brooding and viviparity, are lost at very low rates, if at all. Prolonged care from the egg to later developmental stages evolves from temporally limited care, but it is as easily lost as it is gained. Finally, biparental care is evolutionarily unstable regardless of whether the parents perform complementary or similar care duties. By considering the full spectrum of parental care adaptations, our study reveals a more complex and nuanced picture of how care evolves, is maintained, or is lost.

Parental care can take many forms but how this diversity arises is not well understood. Here, the authors compile data for over 1300 amphibian species and show that different forms of care evolve at different rates, prolonged care can be easily reduced, and biparental care is evolutionarily unstable.

Evolution of the Piscine orthoreovirus Genome Linked to Emergence of Heart and Skeletal Muscle Inflammation in Farmed Atlantic Salmon (Salmo salar)

Heart and skeletal muscle inflammation (HSMI) in farmed Atlantic salmon (Salmo salar) was first diagnosed in Norway in 1999. The disease is caused by Piscine orthoreovirus-1 (PRV-1). The virus is prevalent in farmed Atlantic salmon, but not always associated with disease. Phylogeny and sequence analyses of 31 PRV-1 genomes collected over a 30-year period from fish with or without HSMI, grouped the viral sequences into two main monophylogenetic clusters, one associated with HSMI and the other with low virulent PRV-1 isolates. A PRV-1 strain from Norway sampled in 1988, a decade before the emergence of HSMI, grouped with the low virulent HSMI cluster. The two distinct monophylogenetic clusters were particularly evident for segments S1 and M2. Only a limited number of amino acids were unique to the association with HSMI, and they all located to S1 and M2 encoded proteins. The observed co-evolution of the S1-M2 pair coincided in time with the emergence of HSMI in Norway, and may have evolved through accumulation of mutations and/or segment reassortment. Sequences of S1-M2 suggest selection of the HSMI associated pair, and that this segment pair has remained almost unchanged in Norwegian salmon aquaculture since 1997. PRV-1 strains from the North American Pacific Coast and Faroe Islands have not undergone this evolution, and are more closely related to the PRV-1 precursor strains not associated with clinical HSMI.

The complete plastome of Panax stipuleanatus: Comparative and phylogenetic analyses of the genus Panax (Araliaceae)

Panax stipuleanatus (Araliaceae) is an endangered and medicinally important plant endemic to China. However, phylogenetic relationships within the genus Panax have remained unclear. In this study, we sequenced the complete plastome of P. stipuleanatus and included previously reported Panax plastomes to better understand the relationships between species and plastome evolution within the genus Panax. The plastome of P. stipuleanatus is 156,069 base pairs (bp) in length, consisting of a pair of inverted repeats (IRs, each 25,887 bp) that divide the plastome into a large single copy region (LSC, 86,126 bp) and a small single copy region (SSC, 8169 bp). The plastome contains 114 unigenes (80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes). Comparative analyses indicated that the plastome gene content and order, as well as the expansion/contraction of the IR regions, are all highly conserved within Panax. No significant positive selection in the plastid protein-coding genes was observed across the eight Panax species, suggesting the Panax plastomes may have undergone a strong purifying selection. Our phylogenomic analyses resulted in a phylogeny with high resolution and supports for Panax. Nine protein-coding genes and 10 non-coding regions presented high sequence divergence, which could be useful for identifying different Panax species.

Notes on recombination and reassortment in multipartite/segmented viruses

Besides evolving through nucleotide substitution, viruses frequently also evolve by genetic recombination which can occur when related viral variants co-infect the same cells. Viruses with segmented or multipartite genomes can additionally evolve via the reassortment of genomic components. Various computational techniques are now available for identifying and characterizing recombination and reassortment. While these techniques have revealed both that all well studied segmented and multipartite virus species show some capacity for reassortment, and that recombination is common in many multipartite species, they have indicated that recombination is either rare or does not occur in species with segmented genomes. Reassortment and recombination can make it very difficult to study segmented/multipartite viruses using metagenomics-based approaches. Notable challenges include, both the accurate identification and assignment of genomic components to individual genomes, and the differentiation between natural 'real' recombination events and artifactual 'fake' recombination events arising from the inaccurate de novo assembly of genome component sequences determined using short read sequencing.

Binary-state speciation and extinction method is conditionally robust to realistic violations of its assumptions

Background

Phylogenetic comparative methods allow us to test evolutionary hypotheses without the benefit of an extensive fossil record. These methods, however, make simplifying assumptions, among them that clades are always increasing or stable in diversity, an assumption we know to be false. This study simulates hypothetical clades to test whether the Binary State Speciation and Extinction (BiSSE) method can be used to correctly detect relative differences in diversification rate between ancestral and derived character states even as net diversification rates are declining overall. We simulate clades with declining but positive diversification rates, as well those in which speciation rates decline below extinction rates so that they are losing richness for part of their history. We run these analyses both with simulated symmetric and asymmetric speciation rates to test whether BiSSE can be used to detect them correctly.

Results

For simulations with a neutral character, the fit for a BiSSE model with a neutral character is better than alternative models so long as net diversification rates remain positive. Once net diversification rates become negative, the BiSSE model with the greatest likelihood often has a non-neutral character, even though there is no such character in the simulation. BiSSE’s usefulness in detecting real asymmetry in speciation rates improves with clade age, even well after net diversification rates have become negative.

Conclusions

BiSSE is most useful in analyzing clades of intermediate age, before they have reached peak diversity and gone into decline. After this point, users of BiSSE risk incorrectly inferring differential evolutionary rates when none exist. Fortunately, most studies using BiSSE and similar models focus on rapid, recent diversifications, and are less likely to encounter the biases BiSSE models are subject to for older clades. For extant groups that were once more diverse than now, however, caution should be taken in inferring past diversification patterns without fossil data.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1174-5) contains supplementary material, which is available to authorized users.

Fitness cost of reassortment in human influenza

Reassortment, which is the exchange of genome sequence between viruses co-infecting a host cell, plays an important role in the evolution of segmented viruses. In the human influenza virus, reassortment happens most frequently between co-existing variants within the same lineage. This process breaks genetic linkage and fitness correlations between viral genome segments, but the resulting net effect on viral fitness has remained unclear. In this paper, we determine rate and average selective effect of reassortment processes in the human influenza lineage A/H3N2. For the surface proteins hemagglutinin and neuraminidase, reassortant variants with a mean distance of at least 3 nucleotides to their parent strains get established at a rate of about 10⁻² in units of the neutral point mutation rate. Our inference is based on a new method to map reassortment events from joint genealogies of multiple genome segments, which is tested by extensive simulations. We show that intra-lineage reassortment processes are, on average, under substantial negative selection that increases in strength with increasing sequence distance between the parent strains. The deleterious effects of reassortment manifest themselves in two ways: there are fewer reassortment events than expected from a null model of neutral reassortment, and reassortant strains have fewer descendants than their non-reassortant counterparts. Our results suggest that influenza evolves under ubiquitous epistasis across proteins, which produces fitness barriers against reassortment even between co-circulating strains within one lineage.

The genome of the human influenza virus consists of 8 disjoint RNA polymer segments. These segments can undergo reassortment: when two viruses co-infect a host cell, they can produce viral offspring with a new combination of segments. In this paper, we show that reassortment within a given influenza lineage induces a fitness cost that increases in strength with increasing genetic distance of the parent viruses. Our finding suggests that evolution continuously produces viral proteins whose fitness depends on each other; reassortment reduces fitness by breaking up successful combinations of proteins. Thus, selection across proteins constrains viral evolution within a given lineage, and it may be an important factor in defining a viral species.

Frugivory-related traits promote speciation of tropical palms

Animal-mediated seed dispersal by frugivorous birds and mammals is central to the ecology and functioning of ecosystems, but whether and how frugivory-related traits have affected plant speciation remains little explored. Fruit size is directly linked to plant dispersal capacity and therefore influences gene flow and genetic divergence of plant populations. Using a global species-level phylogeny with comprehensive data on fruit sizes and plant species distributions, we test whether fruit size has affected speciation rates of palms (Arecaceae), a plant family characteristic of tropical rainforests. Globally, the results reveal that palms with small fruit sizes have increased speciation rates compared with those with large (megafaunal) fruits. Speciation of small-fruited palms is particularly high in the understory of tropical rainforests in the New World, and on islands in the Old World. This suggests that frugivory-related traits in combination with geography and the movement behaviour of frugivores can influence the speciation of fleshy-fruited plants.

Characterization of molecular diversity and genome-wide mapping of loci associated with resistance to stripe rust and stem rust in Ethiopian bread wheat accessions

BACKGROUND

The narrow genetic basis of resistance in modern wheat cultivars and the strong selection response of pathogen populations have been responsible for periodic and devastating epidemics of the wheat rust diseases. Characterizing new sources of resistance and incorporating multiple genes into elite cultivars is the most widely accepted current mechanism to achieve durable varietal performance against changes in pathogen virulence. Here, we report a high-density molecular characterization and genome-wide association study (GWAS) of stripe rust and stem rust resistance in 190 Ethiopian bread wheat lines based on phenotypic data from multi-environment field trials and seedling resistance screening experiments. A total of 24,281 single nucleotide polymorphism (SNP) markers filtered from the wheat 90 K iSelect genotyping assay was used to survey Ethiopian germplasm for population structure, genetic diversity and marker-trait associations.

RESULTS

Upon screening for field resistance to stripe rust in the Pacific Northwest of the United States and Ethiopia over multiple growing seasons, and against multiple races of stripe rust and stem rust at seedling stage, eight accessions displayed resistance to all tested races of stem rust and field resistance to stripe rust in all environments. Our GWAS results show 15 loci were significantly associated with seedling and adult plant resistance to stripe rust at false discovery rate (FDR)-adjusted probability (P) <0.10. GWAS also detected 9 additional genomic regions significantly associated (FDR-adjusted P < 0.10) with seedling resistance to stem rust in the Ethiopian wheat accessions. Many of the identified resistance loci were mapped close to previously identified rust resistance genes; however, three loci on the short arms of chromosomes 5A and 7B for stripe rust resistance and two on chromosomes 3B and 7B for stem rust resistance may be novel.

CONCLUSION

Our results demonstrate that considerable genetic variation resides within the landrace accessions that can be utilized to broaden the genetic base of rust resistance in wheat breeding germplasm. The molecular markers identified in this study should be useful in efficiently targeting the associated resistance loci in marker-assisted breeding for rust resistance in Ethiopia and other countries.

Hostplant change and paleoclimatic events explain diversification shifts in skipper butterflies (Family: Hesperiidae)

BACKGROUND

Skippers (Family: Hesperiidae) are a large group of butterflies with ca. 4000 species under 567 genera. The lack of a time-calibrated higher-level phylogeny of the group has precluded understanding of its evolutionary past. We here use a 10-gene dataset to reconstruct the most comprehensive time-calibrated phylogeny of the group, and explore factors that affected the diversification of these butterflies.

RESULTS

Ancestral state reconstructions show that the early hesperiid lineages utilized dicots as larval hostplants. The ability to feed on monocots evolved once at the K-Pg boundary (ca. 65 million years ago (Mya)), and allowed monocot-feeders to diversify much faster on average than dicot-feeders. The increased diversification rate of the monocot-feeding clade is specifically attributed to rate shifts in two of its descendant lineages. The first rate shift, a four-fold increase compared to background rates, happened ca. 50 Mya, soon after the Paleocene-Eocene thermal maximum, in a lineage of the subfamily Hesperiinae that mostly fed on forest monocots. The second rate shift happened ca. 40 Mya in a grass-feeding lineage of Hesperiinae when open-habitat grasslands appeared in the Neotropics owing to gradual cooling of the atmospheric temperature.

CONCLUSIONS

The evolution of monocot feeding strongly influenced diversification of skippers. We hypothesize that although monocot feeding was an intrinsic trait that allowed exploration of novel niches, the lack of extensive availability of monocots comprised an extrinsic limitation for niche exploration. The shifts in diversification rate coincided with paleoclimatic events during which grasses and forest monocots were diversified.