Large-scale Phylogenetic Analysis Research Articles

Uncovering the boundaries of Campylobacter species through large-scale phylogenetic and nucleotide identity analyses.

Campylobacter species are typically helical shaped, Gram-negative, and non-spore-forming bacteria. Species in this genus include established foodborne and animal pathogens as well as emerging pathogens. The accumulation of genomic data from the Campylobacter genus has increased exponentially in recent years, accompanied by the discovery of putative new species. At present, the lack of a standardized species boundary complicates distinguishing established and novel species. We defined the Campylobacter genus core genome (500 loci) using publicly available Campylobacter complete genomes (n = 498) and constructed a core genome phylogeny using 2,193 publicly available Campylobacter genomes to examine inter-species diversity and species boundaries. Utilizing 8,440 Campylobacter genomes representing 33 species and 8 subspecies, we found species delineation based on an average nucleotide identity (ANI) cutoff of 94.2% is consistent with the core genome phylogeny. We identified 60 ANI genomic species that delineated Campylobacter species in concordance with previous comparative genetic studies. All pairwise ANI genomic species pairs had in silico DNA-DNA hybridization scores of less than 70%, supporting their delineation as separate species. We provide the tool Campylobacter Genomic Species typer (CampyGStyper) that assigns ANI genomic species to query genomes based on ANI similarities to medoid genomes from each ANI genomic species with an accuracy of 99.96%. The ANI genomic species definitions proposed here allow consistent species definition in the Campylobacter genus and will facilitate the detection of novel species in the future.IMPORTANCEIn recent years, Campylobacter has gained recognition as the leading cause of bacterial gastroenteritis worldwide, leading to a substantial rise in the collection of genomic data of the Campylobacter genus in public databases. Currently, a standardized Campylobacter species boundary at the genomic level is absent, leading to challenges in detecting emerging pathogens and defining putative novel species within this genus. We used a comprehensive representation of genomes of the Campylobacter genus to construct a core genome phylogenetic tree. Furthermore, we found an average nucleotide identity (ANI) of 94.2% as the optimal cutoff to define the Campylobacter species. Using this cutoff, we identified 60 ANI genomic species which provided a standardized species definition and nomenclature. Importantly, we have developed Campylobacter Genomic Species typer (CampyGStyper), which can robustly and accurately assign these ANI genomic species to Campylobacter genomes, thereby aiding pathogen surveillance and facilitating evolutionary and epidemiological studies of existing and emerging pathogens in the genus Campylobacter.

Patchwork: Alignment-Based Retrieval and Concatenation of Phylogenetic Markers from Genomic Data.

Low-coverage whole-genome sequencing (also known as "genome skimming") is becoming an increasingly affordable approach to large-scale phylogenetic analyses. While already routinely used to recover organellar genomes, genome skimming is rather rarely utilized for recovering single-copy nuclear markers. One reason might be that only few tools exist to work with this data type within a phylogenomic context, especially to deal with fragmented genome assemblies. We here present a new software tool called Patchwork for mining phylogenetic markers from highly fragmented short-read assemblies as well as directly from sequence reads. Patchwork is an alignment-based tool that utilizes the sequence aligner DIAMOND and is written in the programming language Julia. Homologous regions are obtained via a sequence similarity search, followed by a "hit stitching" phase, in which adjacent or overlapping regions are merged into a single unit. The novel sliding window algorithm trims away any noncoding regions from the resulting sequence. We demonstrate the utility of Patchwork by recovering near-universal single-copy orthologs within a benchmarking study, and we additionally assess the performance of Patchwork in comparison with other programs. We find that Patchwork allows for accurate retrieval of (putatively) single-copy genes from genome skimming data sets at different sequencing depths with high computational speed, outperforming existing software targeting similar tasks. Patchwork is released under the GNU General Public License version 3. Installation instructions, additional documentation, and the source code itself are all available via GitHub at https://github.com/fethalen/Patchwork.

CHESS 3: an improved, comprehensive catalog of human genes and transcripts based on large-scale expression data, phylogenetic analysis, and protein structure

CHESS 3 represents an improved human gene catalog based on nearly 10,000 RNA-seq experiments across 54 body sites. It significantly improves current genome annotation by integrating the latest reference data and algorithms, machine learning techniques for noise filtering, and new protein structure prediction methods. CHESS 3 contains 41,356 genes, including 19,839 protein-coding genes and 158,377 transcripts, with 14,863 protein-coding transcripts not in other catalogs. It includes all MANE transcripts and at least one transcript for most RefSeq and GENCODE genes. On the CHM13 human genome, the CHESS 3 catalog contains an additional 129 protein-coding genes. CHESS 3 is available at http://ccb.jhu.edu/chess.

Large-scale phylogenetic analysis reveals genetic diversity and geographic distribution of rabies virus in South-East and South Asia.

South-East Asia (SEA) and South Asia (SA) are two important geographic regions with the most severe enzootic rabies in the world. In these regions, phylogenetic analysis of rabies virus (RABV) has been conducted only at a country level; the results obtained from different countries are scattered and unequal, with a non-uniform system to name RABV genotypes. Therefore, it is difficult to undertake origin-tracking and compare inter-country RABV evolution and transmission. To avoid the confusion in understanding and to generate a panoramic picture of RABV genetic diversity, distribution, and transmission in SEA and SA, the present study conducted a systematic phylogenetic analysis by combining all sequences representing 2368 RABV strains submitted to GenBank by 14 rabies endemic SEA and SA countries. The results showed that RABVs circulating in two regions were classified into four major clades and many subclades: the Asia clade is circulating only in SEA, the Indian subcontinent, and Arctic-like clades only in SA, while the Cosmopolitan clade has been detected in both regions. The results also showed a wide range of hosts were infected by divergent RABV subclades, with dogs being the major transmission source. However, wildlife rabies was also found to be an important issue with 6 wild carnivore species identified as potential sources of spillover risk for sylvatic rabies to humans, domestic animals, and other wild animals. Current findings indicate that the two regions have separate virus clades circulating thus indicating the absence of cross-transmission between the regions. The study emphasizes the importance of phylogenetic analysis in the regions using uniform genotyping and naming systems for rabies surveillance, to coordinate actions of member countries to eliminate dog-mediated human rabies by 2030.

Human-to-swine introductions and onward transmission of 2009 H1N1 pandemic influenza viruses in Brazil.

Once established in the human population, the 2009 H1N1 pandemic virus (H1N1pdm09) was repeatedly introduced into swine populations globally with subsequent onward transmission among pigs. To identify and characterize human-to-swine H1N1pdm09 introductions in Brazil, we conducted a large-scale phylogenetic analysis of 4,141 H1pdm09 hemagglutinin (HA) and 3,227 N1pdm09 neuraminidase (NA) gene sequences isolated globally from humans and swine between 2009 and 2022. Phylodynamic analysis revealed that during the period between 2009 and 2011, there was a rapid transmission of the H1N1pdm09 virus from humans to swine in Brazil. Multiple introductions of the virus were observed, but most of them resulted in self-limited infections in swine, with limited onward transmission. Only a few sustained transmission clusters were identified during this period. After 2012, there was a reduction in the number of human-to-swine H1N1pdm09 transmissions in Brazil. The virus underwent continuous antigenic drift, and a balance was established between swine-to-swine transmission and extinction, with minimal sustained onward transmission from humans to swine. These results emphasize the dynamic interplay between human-to-swine transmission, antigenic drift, and the establishment of swine-to-swine transmission in shaping the evolution and persistence of H1N1pdm09 in swine populations.

Functional divergence of CYP76AKs shapes the chemodiversity of abietane-type diterpenoids in genus Salvia

The genus Salvia L. (Lamiaceae) comprises myriad distinct medicinal herbs, with terpenoids as one of their major active chemical groups. Abietane-type diterpenoids (ATDs), such as tanshinones and carnosic acids, are specific to Salvia and exhibit taxonomic chemical diversity among lineages. To elucidate how ATD chemical diversity evolved, we carried out large-scale metabolic and phylogenetic analyses of 71 Salvia species, combined with enzyme function, ancestral sequence and chemical trait reconstruction, and comparative genomics experiments. This integrated approach showed that the lineage-wide ATD diversities in Salvia were induced by differences in the oxidation of the terpenoid skeleton at C-20, which was caused by the functional divergence of the cytochrome P450 subfamily CYP76AK. These findings present a unique pattern of chemical diversity in plants that was shaped by the loss of enzyme activity and associated catalytic pathways.

Identification of a novel lipoic acid biosynthesis pathway reveals the complex evolution of lipoate assembly in prokaryotes.

Lipoic acid is an essential biomolecule found in all domains of life and is involved in central carbon metabolism and dissimilatory sulfur oxidation. The machineries for lipoate assembly in mitochondria and chloroplasts of higher eukaryotes, as well as in the apicoplasts of some protozoa, are all of prokaryotic origin. Here, we provide experimental evidence for a novel lipoate assembly pathway in bacteria based on a sLpl(AB) lipoate:protein ligase, which attaches octanoate or lipoate to apo-proteins, and 2 radical SAM proteins, LipS1 and LipS2, which work together as lipoyl synthase and insert 2 sulfur atoms. Extensive homology searches combined with genomic context analyses allowed us to precisely distinguish between the new and established pathways and map them on the tree of life. This not only revealed a much wider distribution of lipoate biogenesis systems than expected, in particular, the novel sLpl(AB)-LipS1/S2 pathway, and indicated a highly modular nature of the enzymes involved, with unforeseen combinations, but also provided a new framework for the evolution of lipoate assembly. Our results show that dedicated machineries for both de novo lipoate biogenesis and scavenging from the environment were implemented early in evolution and that their distribution in the 2 prokaryotic domains was shaped by a complex network of horizontal gene transfers, acquisition of additional genes, fusions, and losses. Our large-scale phylogenetic analyses identify the bipartite archaeal LplAB ligase as the ancestor of the bacterial sLpl(AB) proteins, which were obtained by horizontal gene transfer. LipS1/S2 have a more complex evolutionary history with multiple of such events but probably also originated in the domain archaea.

Developmental strategies underlying gigantism and miniaturization in non-avialan theropod dinosaurs.

In amniotes, the predominant developmental strategy underlying body size evolution is thought to be adjustments to the rate of growth rather than its duration. However, most theoretical and experimental studies supporting this axiom focus on pairwise comparisons and/or lack an explicit phylogenetic framework. We present the first large-scale phylogenetic comparative analysis examining developmental strategies underlying the evolution of body size, focusing on non-avialan theropod dinosaurs. We reconstruct ancestral states of growth rate and body mass in a taxonomically rich dataset, finding that contrary to expectations, changes in the rate and duration of growth played nearly equal roles in the evolution of the vast body size disparity present in non-avialan theropods-and perhaps that of amniotes in general.

Global population structure and genomic surveillance framework of carbapenem-resistant Salmonella enterica

Due to the frequent international and intercontinental transmission of multidrug-resistant bacteria, it is imperative to understand the epidemiology, phylogeography, and population structure of carbapenem-resistant Salmonella enterica (CRSE) across the globe. During the period of 2015–2022, two blaNDM-carrying S. enterica strains were recovered from 3695 Salmonella strains in four hospitals in China. The global phylogenetic framework and geographical distribution of CRSE were defined by our recently updated bacterial whole genome sequence typing and source tracking database BacWGSTdb 2.0 to measure the diversity and evolutionary relatedness in context with epidemiological metadata. Phylogeny for all carbapenemase gene-harboring plasmids in S. enterica based on the pairwise Mash differences was also constructed to evaluate the potential transmission of these plasmids in a global context. A large-scale phylogenetic analysis grouped global CRSE into nine distinct clades. The small genetic distance (< 20 SNPs) between 198 pairs of CRSE suggested the presence of clonal transmission. Global CRSE have significant geographical variations, which was associated with the clonal lineages and carbapenemase genes. Carbapenemase gene-carrying plasmids with a high degree of similarity have surfaced in various hosts and countries. The widespread of multiple-replicon plasmids that offer a great capacity to accommodate multiple antimicrobial resistance genes is continuously enhancing the potential risk of CRSE isolates to propagate globally. Both clonal spread of strains and horizontal transfer of carbapenemase gene-harboring plasmids contribute to the global dissemination of CRSE. Our findings on the worldwide spread and transmission dynamics of this emerging bacterium has increased the knowledge of its global epidemics. Continued epidemiological surveillance is necessary to prevent global outbreak of multidrug-resistant Salmonella infections.

Variation in plant Toll/Interleukin-1 receptor domain protein dependence on ENHANCED DISEASE SUSCEPTIBILITY 1.

Toll/Interleukin-1 receptor (TIR) domains are integral to immune systems across all kingdoms. In plants, TIRs are present in nucleotide-binding leucine-rich repeat (NLR) immune receptors, NLR-like, and TIR-only proteins. Although TIR-NLR and TIR signaling in plants require the ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) protein family, TIRs persist in species that have no EDS1 members. To assess whether particular TIR groups evolved with EDS1, we searched for TIR-EDS1 co-occurrence patterns. Using a large-scale phylogenetic analysis of TIR domains from 39 algal and land plant species, we identified 4 TIR families that are shared by several plant orders. One group occurred in TIR-NLRs of eudicots and another in TIR-NLRs across eudicots and magnoliids. Two further groups were more widespread. A conserved TIR-only group co-occurred with EDS1 and members of this group elicit EDS1-dependent cell death. In contrast, a maize (Zea mays) representative of TIR proteins with tetratricopeptide repeats was also present in species without EDS1 and induced EDS1-independent cell death. Our data provide a phylogeny-based plant TIR classification and identify TIRs that appear to have evolved with and are dependent on EDS1, while others have EDS1-independent activity.

Towards the phylogenetic placement of the enigmatic African genus Prolabeops Schultz, 1941.

The small cyprinid genus Prolabeops Schultz, 1941 is restricted to the Nyong and Sanaga River systems in Cameroon. In the past, the genus had been suggested to be either a member of the Labeoninae, Torinae or the Smiliogastrinae mainly on the basis of morphological similarities, and it is nowadays considered as incertae sedis within the Cypriniformes. This study provides the first attempt to reveal the phylogenetic position of Prolabeops using molecular data. For this purpose, the authors sequenced a large fraction of the mitochondrial genome (c. 13,600 bp), including all mitochondrial protein coding genes, of two Prolabeops melanhypopterus specimens and an additional four Enteromius specimens. The large-scale phylogenetic analysis was based on an alignment including all mitochondrial protein coding genes of 902 specimens representing c. 899 cypriniform species. Prolabeops was clearly recovered within the African Smiliogastrinae, forming a weakly supported clade together with Enteromius jae, Enteromius hulstaerti and Barboides gracilis. The study data underline the urgent need of a thorough taxonomic revision of the small African barbs collectively placed in the genus Enteromius.

Taxonomic, palaeobiological and evolutionary implications of a phylogenetic hypothesis for Ornithischia (Archosauria: Dinosauria)

Abstract The origin and evolutionary relationships of ornithischian dinosaurs are topics that have undergone a series of substantial revisions. At present there are several competing hypotheses concerning the relationship between Ornithischia and the other principal clades of Dinosauria. Some hypotheses have posited a tree topology within Dinosauria that imply a ‘ghost-lineage’ for Ornithischia (whose representatives make their first unambiguous appearance in the Hettangian) that extends through a substantial portion of Triassic time. In contrast, other hypotheses have placed conventionally Triassic dinosauromorph (stem-lineage Dinosauria) taxa within the clade Ornithischia. Recently, a large-scale phylogenetic analysis recovered an array of taxa, known as ‘silesaurids’, as a paraphyletic assemblage of taxa (referred to in this article using the informal terms silesaurs or silesaurians) on the branch leading to the clade Ornithischia. This latter hypothesis of relationships would account for the apparent absence of Triassic ornithischians, because stem-lineage ornithischians (silesaurs in this article) are exclusively Triassic. However, the analysis that produced this novel topology used a dataset that, in its original form, did not include all early representatives of Ornithischia (sensu lato), and did not incorporate all the anatomical characters that have been suggested to unite Ornithischia with other dinosaurian clades (Theropoda and Sauropodomorpha). Nor did the initial study go on to expand upon some important taxonomic, palaeobiological and evolutionary implications of a topology that links a paraphyletic array of silesaurs to the clade Ornithischia. The present article addresses these latter issues by expansion and re-analysis of the original dataset. The results find further support for the hypothesis that silesaurs comprise a paraphyletic grouping of taxa on the stem of Ornithischia and that successive silesaur taxa acquire anatomical characters anagenetically in a process that culminates in the assembly of what may be described as a ‘traditional’ ornithischian. The overall topology of the consensus tree remains but little changed from the original analysis, despite the addition of new taxa and characters. To provide stability to this area of the tree and to preserve the most important of the relevant taxonomic names, we suggest a revised taxonomic framework for ornithischians that is consistent with this new topology. We retain the name Ornithischia for the total-group (traditional Ornithischia and its stem-lineage), while we resuscitate a name originally proposed by Richard Owen, Prionodontia (= ‘coarse edged teeth’) for the clade containing only the so-called traditional ornithischian (= ‘bird-hipped’) dinosaurs. We also erect Parapredentata as a more exclusive subclade in Ornithischia. This novel taxonomic framework is intended to provide phylogenetic clarity and a degree of stability in Ornithischia and Dinosauria as further analyses and new data continue to refine and re-shape the tree. The data presented in this study represent a stage in our attempt to establish an early dinosaur dataset in which character definitions and character scores are agreed upon and used consistently.

Large-Scale Phylogenetic Analysis Reveals a New Genetic Clade among Escherichia coli O26 Strains.

ABSTRACTShiga toxin-producing Escherichia coli (STEC) O26 is the predominant non-O157 serogroup causing hemolytic uremic syndrome worldwide. Moreover, the serogroup is highly dynamic and harbors several pathogenic clones. Here, we investigated the phylogenetic relationship of STEC O26 at a global level based on 1,367 strains from 20 countries deposited in NCBI and Enterobase databases. The whole-genome-based analysis identified a new genetic clade, called ST29C4. The new clade was unique in terms of multilocus sequence type (ST29), CRISPR (group Ia), and dominant plasmid gene profile (ehxA+/katP-/espP-/etpD-). Moreover, the combination of multiple typing methods (core genome single nucleotide polymorphism [SNP] typing, CRISPR typing, and virulence genes analysis) demonstrated that this new lineage ST29C4 was in the intermediate phylogenetic position between ST29C3 and other non-ST29C3 strains. Besides, we observed that ST29C4 harbored extraintestinal pathogenic E. coli (ExPEC)-related virulence gene (VG), tsh, and STEC-associated VG, stx2a, suggesting the emergence of a hybrid pathogen. The ST29C4 strains also exhibited high similarity in stx2a-prophage and integrase with the O104:H4 strain, further demonstrating its potential risk to human health. Collectively, the large-scale phylogenetic analysis extends the understanding of the clonal structure of O26 strains and provides new insights for O26 strain microevolution.IMPORTANCE Shiga toxin-producing Escherichia coli (STEC) O26 is the second prevalent STEC serogroup only to O157, which can cause a series of diseases ranging from mild diarrhea to life-threatening hemolytic uremic syndrome (HUS). The serogroup is highly diverse and multiple clones are characterized, including ST29C1-C3 and ST21C1-C2. However, the phylogenetic relationship of these clones remains fully unclear. In this study, we revealed a new genetic clade among O26 strains, ST29C4, which was unique in terms of CRISPR, multilocus sequence type (MLST), and plasmid gene profile (PGP). Moreover, the combination of multiple typing methods demonstrated that this new clone was located in the intermediate phylogenetic position between ST29C3 and other non-ST29C3 strains (i.e., ST29C1-C2 and ST21C1-C2). Overall, the large-scale phylogenetic analysis extends our current understanding of O26 microevolution.

Origin, evolution and diversification of plant ARGONAUTE proteins.

Argonaute (AGO) proteins are central players in RNA interference in eukaryotes. They associate with small RNAs (sRNA) and lead to transcriptional or posttranscriptional silencing of targets, thereby regulating diverse biological processes. The molecular and biological functions of AGO proteins have been extensively characterized, particularly in a few angiosperm species, leading to the recognition that the AGO family has expanded to accommodate diverse sRNAs thereby performing diverse biological functions. However, understanding of the expansion of AGO proteins in plants is still limited, due to a dearth of knowledge of AGO proteins in green algal groups. Here, we identified more than 2900 AGO proteins from 244 plant species, including green algae, and performed a large-scale phylogenetic analysis. The phylogeny shows that the plant AGO family gave rise to four clades after the emergence of hydrobiontic algae and prior to the emergence of land plants. Subsequent parallel expansion in ferns and angiosperms resulted in eight main clades in angiosperms: AGO2, AGO7, AGO6, AGO4, AGO1, AGO10a, AGO10b and AGO5. On the basis of this phylogeny, we identified two novel AGO4 orthologs that Arabidopsis does not have, and redefined AGO10, which is composed of AGO10a and AGO10b. Finally, we propose a hypothetical evolutionary model of AGO proteins in plants. Our studies provide a deeper understanding of the phylogenetic relationships of AGO family members in the green lineage, which would help to further reveal their roles as RNAi effectors.

Nucleocapsid mutations R203K/G204R increase the infectivity, fitness, and virulence of SARS-CoV-2

Previous work found that the co-occurring mutations R203K/G204R on the SARS-CoV-2 nucleocapsid (N) protein are increasing in frequency among emerging variants of concern or interest. Through a combination of in silico analyses, this study demonstrates that R203K/G204R are adaptive, while large-scale phylogenetic analyses indicate that R203K/G204R associate with the emergence of the high-transmissibility SARS-CoV-2 lineage B.1.1.7. Competition experiments suggest that the 203K/204R variants possess a replication advantage over the preceding R203/G204 variants, possibly related to ribonucleocapsid (RNP) assembly. Moreover, the 203K/204R virus shows increased infectivity in human lung cells and hamsters. Accordingly, we observe a positive association between increased COVID-19 severity and sample frequency of 203K/204R. Our work suggests that the 203K/204R mutations contribute to the increased transmission and virulence of select SARS-CoV-2 variants. In addition to mutations in the spike protein, mutations in the nucleocapsid protein are important for viral spreading during the pandemic.

New spinosaurids from the Wessex Formation (Early Cretaceous, UK) and the European origins of Spinosauridae

Spinosaurids are among the most distinctive and yet poorly-known of large-bodied theropod dinosaurs, a situation exacerbated by their mostly fragmentary fossil record and competing views regarding their palaeobiology. Here, we report two new Early Cretaceous spinosaurid specimens from the Wessex Formation (Barremian) of the Isle of Wight. Large-scale phylogenetic analyses using parsimony and Bayesian techniques recover the pair in a new clade within Baryonychinae that also includes the hypodigm of the African spinosaurid Suchomimus. Both specimens represent distinct and novel taxa, herein named Ceratosuchops inferodios gen. et sp. nov. and Riparovenator milnerae gen. et sp. nov. A palaeogeographic reconstruction suggests a European origin for Spinosauridae, with at least two dispersal events into Africa. These new finds provide welcome information on poorly sampled areas of spinosaurid anatomy, suggest that sympatry was present and potentially common in baryonychines and spinosaurids as a whole, and contribute to updated palaeobiogeographic reconstructions for the clade.

Cophylogenetic analyses of Trachymyrmex ant-fungal specificity: "One to one with some exceptions".

Over the past few decades, large-scale phylogenetic analyses of fungus-gardening ants and their symbiotic fungi have depicted strong concordance among major clades of ants and their symbiotic fungi, yet within clades, fungus sharing is widespread among unrelated ant lineages. Sharing has been explained using a diffuse coevolution model within major clades. Understanding horizontal exchange within clades has been limited by conventional genetic markers that lack both interspecific and geographic variation. To examine whether reports of horizontal exchange were indeed due to symbiont sharing or the result of employing relatively uninformative molecular markers, samples of Trachymyrmex arizonensis and Trachymyrmex pomonae and their fungi were collected from native populations in Arizona and genotyped using conventional marker genes and genome-wide single nucleotide polymorphisms (SNPs). Conventional markers of the fungal symbionts generally exhibited cophylogenetic patterns that were consistent with some symbiont sharing, but most fungal clades had low support. SNP analysis, in contrast, indicated that each ant species exhibited fidelity to its own fungal subclade with only one instance of a colony growing a fungus that was otherwise associated with a different ant species. This evidence supports a pattern of codivergence between Trachymyrmex species and their fungi, and thus a diffuse coevolutionary model may not accurately predict symbiont exchange. These results suggest that fungal sharing across host species in these symbioses may be less extensive than previously thought.

Origin and evolution of nonulosonic acid synthases and their relationship with bacterial pathogenicity revealed by a large-scale phylogenetic analysis.

Nonulosonic acids (NulOs) are a group of nine-carbon monosaccharides with different functions in nature. N-acetylneuraminic acid (Neu5Ac) is the most common NulO. It covers the membrane surface of all human cells and is a central molecule in the process of self-recognition via SIGLECS receptors. Some pathogenic bacteria escape the immune system by copying the sialylation of the host cell membrane. Neu5Ac production in these bacteria is catalysed by the enzyme NeuB. Some bacteria can also produce other NulOs named pseudaminic and legionaminic acids, through the NeuB homologues PseI and LegI, respectively. In Opisthokonta eukaryotes, the biosynthesis of Neu5Ac is catalysed by the enzyme NanS. In this study, we used publicly available data of sequences of NulOs synthases to investigate its distribution within the three domains of life and its relationship with pathogenic bacteria. We mined the KEGG database and found 425 NeuB sequences. Most NeuB sequences (58.74 %) from the KEGG orthology database were classified as from environmental bacteria; however, sequences from pathogenic bacteria showed higher conservation and prevalence of a specific domain named SAF. Using the HMM profile we identified 13 941 NulO synthase sequences in UniProt. Phylogenetic analysis of these sequences showed that the synthases were divided into three main groups that can be related to the lifestyle of these bacteria: (I) predominantly environmental, (II) intermediate and (III) predominantly pathogenic. NeuB was widely distributed in the groups. However, LegI and PseI were more concentrated in groups II and III, respectively. We also found that PseI appeared later in the evolutionary process, derived from NeuB. We use this same methodology to retrieve sialic acid synthase sequences from Archaea and Eukarya. A large-scale phylogenetic analysis showed that while the Archaea sequences are spread across the tree, the eukaryotic NanS sequences were grouped in a specific branch in group II. None of the bacterial NanS sequences grouped with the eukaryotic branch. The analysis of conserved residues showed that the synthases of Archaea and Eukarya present a mutation in one of the three catalytic residues, an E134D change, related to a Neisseria meningitidis reference sequence. We also found that the conservation profile is higher between NeuB of pathogenic bacteria and NanS of eukaryotes than between NeuB of environmental bacteria and NanS of eukaryotes. Our large-scale analysis brings new perspectives on the evolution of NulOs synthases, suggesting their presence in the last common universal ancestor.

Genetic evidence for sex-biased dispersal and cryptic diversity in the greater horseshoe bat, Rhinolophus ferrumequinum

Dispersal plays an important role in the ecological and evolutionary processes of natural populations. Mating behavior (or mating system) is a critical factor shaping dispersal patterns and extents in social mammals, sometimes driving the evolution of sex-biased dispersal. Using molecular markers with contrasting modes of inheritance (mitochondrial DNA and nuclear microsatellites), we determined the population genetic structure and evolutionary history of the great horseshoe bat, Rhinolophus ferrumequinum inhabiting eleven national parks of South Korea, being known as a biodiversity hotspot. Despite apparent matrilineal structure observed over space, there was weak nuclear geographic structure, suggesting female philopatry with male-biased dispersal. The analyses indicated the signal of nonrandom mating (i.e. inbreeding), which is at least partly due to female’s sedentary lifestyle. The large-scale phylogenetic analysis revealed unexpected deep divergence among three distinct clades (Southwest China, East China, and Northeast Asia including South Korea), suggesting these may possibly represent cryptic species complex in R. ferrumequinum. Our findings of sex differences in dispersal in R. ferrumequinum inform conservation managements to enhance the population connectivity among the national parks through promoting female dispersal. Our study also highlights cryptic species diversity in a temperate bat that should have a priority for conservation.

Establishment and lineage dynamics of the SARS-CoV-2 epidemic in the UK.

The UK’s COVID-19 epidemic during early 2020 was one of world’s largest and unusually well represented by virus genomic sampling. Here we reveal the fine-scale genetic lineage structure of this epidemic through analysis of 50,887 SARS-CoV-2 genomes, including 26,181 from the UK sampled throughout the country’s first wave of infection. Using large-scale phylogenetic analyses, combined with epidemiological and travel data, we quantify the size, spatio-temporal origins and persistence of genetically-distinct UK transmission lineages. Rapid fluctuations in virus importation rates resulted in >1000 lineages; those introduced prior to national lockdown tended to be larger and more dispersed. Lineage importation and regional lineage diversity declined after lockdown, while lineage elimination was size-dependent. We discuss the implications of our genetic perspective on transmission dynamics for COVID-19 epidemiology and control.