Gotree/Goalign: toolkit and Go API to facilitate the development of phylogenetic workflows

Repeated loss of function at HD mating-type genes and of recombination in anther-smut fungi

Basidiomycete fungi typically have two mating-type loci controlling mating compatibility, HD and PR, residing on different chromosomes. Loss-of-function in mating compatibility has been reported at the PR genes in a few heterothallic basidiomycetes, but not for the HD genes. In Microbotryum anther-smut fungi, there have been repeated linkage events between the HD and PR loci through chromosome fusions, leading to non-recombining regions. Here, we found that two sister Microbotryum species parasitizing Dianthus plants, M. superbum and M. shykoffianum, as well as the distantly related M. scorzonerae, have their HD and PR loci on different chromosomes, but with the PR chromosome fused with a part of the ancestral HD chromosome. In addition, recombination suppression has extended stepwise, generating evolutionary strata. In all three species, the HD genes lost their function in mating compatibility, natural diploid strains being often homozygous at the HD locus. Strains could be homozygous for a disrupted HD2 gene, that was hardly expressed during mating. Mating tests confirmed that a single genetic factor controlled mating compatibility and that haploid strains with identical HD alleles could mate and produce hyphae. This study shows that a unifactorial mating-type determinism can evolve, repeatedly, from a bifactorial system, by different mechanisms.

Unique trajectory of gene family evolution from genomic analysis of nearly all known species in an ancient yeast lineage

Gene gains and losses are a major driver of genome evolution; their precise characterization can provide insights into the origin and diversification of major lineages. Here, we examined gene family evolution of 1154 genomes from nearly all known species in the medically and technologically important yeast subphylum Saccharomycotina. We found that yeast gene family evolution differs from that of plants, animals, and filamentous ascomycetes, and is characterized by smaller overall gene numbers yet larger gene family sizes for a given gene number. Faster-evolving lineages (FELs) in yeasts experienced significantly higher rates of gene losses-commensurate with a narrowing of metabolic niche breadth-but higher speciation rates than their slower-evolving sister lineages (SELs). Gene families most often lost are those involved in mRNA splicing, carbohydrate metabolism, and cell division and are likely associated with intron loss, metabolic breadth, and non-canonical cell cycle processes. Our results highlight the significant role of gene family contractions in the evolution of yeast metabolism, genome function, and speciation, and suggest that gene family evolutionary trajectories have differed markedly across major eukaryotic lineages.

A multi-hospital, clinician-initiated bacterial genomics programme to investigate treatment failure in severe Staphylococcus aureus infections

Bacterial genomics is increasingly used for infectious diseases surveillance, outbreak detection and prediction of antibiotic resistance. With expanding availability of rapid whole-genome sequencing, bacterial genomics data could become a valuable tool for clinicians managing bacterial infections, driving precision medicine strategies. Here, we present a clinician-driven bacterial genomics framework that applies within-patient evolutionary analysis to identify in real-time microbial genetic changes that have an impact on treatment outcomes of severe Staphylococcus aureus infections, a strategy that is increasingly used in cancer genomics. Our approach uses a combination of bacterial genomics and antibiotic susceptibility testing to identify and track bacterial adaptive mutations that underlie microbiologically documented treatment failure (i.e. ongoing positive cultures [persistent infection] or new positive cultures after initial response [recurrent infection]). We show the potential added value of our approach to clinicians and propose a roadmap for the use of bacterial genomics to advance the management of severe bacterial infections.

Too Far From Relatives? Impact of the Genetic Distance on the Success of Exon Capture in Phylogenomics

The exon capture approach allows for sequencing a large number of loci to reconstruct phylogenetic relationships at varying taxonomic levels. In order to efficiently recover the targeted loci, the probes designed to capture the exons need to be genetically sufficiently similar to bind to the DNA, with a proposed limit of 10% of divergence. However, this threshold has never been tested with a specific protocol. We have designed a set of probes to capture 1125 exons in the Neogastropoda (Mollusca, Gastropoda), processed with the same protocol from the field to the DNA sequencing to control for potential bias in DNA quantity and quality. We sequenced 30 different species, including 14 species of Neogastropoda and 16 species of Caenogastropoda non-Neogastropoda. Each species includes five specimens, for a total of 150 specimens, and for four specimens among the 150, DNA extracts were aliquoted in four samples, sequenced separately, to estimate the intraspecific and intraspecimen variability in capture success. Our results confirm the impact of genetic distance on the success of exon capture with a negative linear correlation between the genetic distance and the number of exons captured for each sample. Consequently, designing new capture probes would allow for capturing exons in genetically more distant groups without the need to redesign a new set of exons.

The genomic consequences and persistence of sociality in spiders

In cooperatively breeding social animals, a few individuals account for all reproduction. In some taxa, sociality is accompanied by a transition from outcrossing to inbreeding. In concert, these traits reduce effective population size, potentially rendering transitions to sociality "evolutionarily dead-ends." We addressed this hypothesis in a comparative genomic study in spiders, in which sociality has evolved independently at least 23 times, but social branches are recent and short. We present genomic evidence for the evolutionary dead-end hypothesis in a spider genus with three independent transitions to sociality. We assembled and annotated high-quality, chromosome-level reference genomes from three pairs of closely related social and subsocial Stegodyphus species. We timed the divergence between the social and subsocial species pairs to be from 1.3 million to 1.8 million years. Social evolution in spiders involves a shift from outcrossing to inbreeding and from an equal to a female-biased sex ratio, causing severe reductions in effective population size and decreased efficacy of selection. We show that transitions to sociality only had full effect on purifying selection at 119, 260, and 279 kya, respectively, and follow similar convergent trajectories of progressive loss of diversity and shifts to an increasingly female-biased sex ratio. This almost deterministic genomic response to sociality may explain why social spider lineages do not persist. What causes species extinction is not clear, but either could be selfish meiotic drive eliminating the production of males or could be an inability to retain genome integrity in the face of extremely reduced efficacy of selection.

Molecular characterization and prevalence assessment of Durettenema sp. (Nematoda: Trichostrongyloidea) in the great leaf-nosed bats (Hipposideros armiger) in South China

The viruses associated with bats have generated significant concern; however, there is limited knowledge regarding the endoparasites that affect these mammals. This study involved the collection of seven nematode specimens (three males and four females) from the intestines of Hipposideros armiger in Shaoguan City, Guangdong, China. Next-generation sequencing was employed to obtain the mitochondrial DNA (mtDNA) genome, which was determined to be 14,130 base pairs in length. The mitochondrial genome comprised 12 protein-coding genes, 21 tRNA genes, 2 rRNA genes, and an AT-rich non-coding region. Phylogenetic analyses based on mtDNA sequences indicated that the nematode forms a sister clade to Nematodirus, exhibiting only 74% nucleotide identity. In contrast, the nuclear ITS1 gene demonstrated a high degree of nucleotide identity (98.6%-98.8%) with Durettenema guangdongense. Consequently, the parasitic nematode identified from H. armiger is likely to belong to the genus Durettenema and has been designated as Durettenema sp. 888. Furthermore, an epidemiological investigation revealed the presence of the parasitic nematode infections in H. armiger collected from Guangdong, Guangxi, and Guizhou Provinces. Given the widespread distribution of H. armiger and their tendency to inhabit areas in close proximity to human dwellings, the influence of parasite prevalence on bat population numbers and potential for human and domestic animal transmission of this pathogen warrants further investigation.

Whole-genome automated assembly pipeline for Chlamydia trachomatis strains from reference, in vitro and clinical samples using the integrated CtGAP pipeline

Whole genome sequencing (WGS) is pivotal for the molecular characterization of Chlamydia trachomatis (Ct)-the leading bacterial cause of sexually transmitted infections and infectious blindness worldwide. Ct WGS can inform epidemiologic, public health and outbreak investigations of these human-restricted pathogens. However, challenges persist in generating high-quality genomes for downstream analyses given its obligate intracellular nature and difficulty with in vitro propagation. No single tool exists for the entirety of Ct genome assembly, necessitating the adaptation of multiple programs with varying success. Compounding this issue is the absence of reliable Ct reference strain genomes. We, therefore, developed CtGAP-Chlamydia trachomatisGenome Assembly Pipeline-as an integrated 'one-stop-shop' pipeline for assembly and characterization of Ct genome sequencing data from various sources including isolates, in vitro samples, clinical swabs and urine. CtGAP, written in Snakemake, enables read quality statistics output, adapter and quality trimming, host read removal, de novo and reference-guided assembly, contig scaffolding, selective ompA, multi-locus-sequence and plasmid typing, phylogenetic tree construction, and recombinant genome identification. Twenty Ct reference genomes were also generated. Successfully validated on a diverse collection of 363 samples containing Ct, CtGAP represents a novel pipeline requiring minimal bioinformatics expertise with easy adaptation for use with other bacterial species.

Reduced evolutionary constraint accompanies ongoing radiation in deep-sea anglerfishes

Colonization of a novel habitat is often followed by phenotypic diversification in the wake of ecological opportunity. However, some habitats should be inherently more constraining than others if the challenges of that environment offer few evolutionary solutions. We examined this push-and-pull on macroevolutionary diversification following habitat transitions in the anglerfishes (Lophiiformes). We constructed a phylogeny with extensive sampling (1,092 loci and ~38% of species), combined with three-dimensional phenotypic data from museum specimens. We used these datasets to examine the tempo and mode of phenotypic diversification. The deep-sea pelagic anglerfishes originated from a benthic ancestor and shortly after experienced rapid lineage diversification rates. This transition incurred shifts towards larger jaws, smaller eyes and a more laterally compressed body plan. Despite these directional trends, this lineage still evolved high phenotypic disparity in body, skull and jaw shapes. In particular, bathypelagic anglerfishes show high variability in body elongation, while benthic anglerfishes are constrained around optimal shapes. Within this radiation, phenotypic evolution was concentrated among recently diverged lineages, notably those that deviated from the archetypical globose body plan. Taken together, these results demonstrate that spectacular evolutionary radiations can unfold even within environments with few ecological resources and demanding physiological challenges.

Convergent evolution of oxidized sugar metabolism in commensal and pathogenic microbes in the inflamed gut

Inflammation-associated perturbations of the gut microbiome are well characterized, but poorly understood. Here, we demonstrate that disparate taxa recapitulate the metabolism of the oxidized sugars glucarate and galactarate, utilizing enzymatically divergent, yet functionally equivalent, gud/gar pathways. The divergent pathway in commensals includes a putative 5-KDG aldolase (GudL) and an uncharacterized ABC transporter (GarABC) that recapitulate the function of their non-homologous counterparts in pathogens. A systematic bioinformatic search for the gud/gar pathway in gut microbes identified 887 species putatively capable of metabolizing oxidized sugars. Previous studies showed that inflammation-derived nitrate, formed by nitric oxide reacting with superoxide, promotes pathogen growth. Our findings reveal a parallel phenomenon: oxidized sugars, also produced from reactions with nitric oxide, serve as alternative carbon sources for commensal microbes. Previously considered a pathogen virulence factor, oxidized sugar metabolism is also present in specific commensals and may contribute to their increased relative abundance in gastrointestinal inflammation.

Genome sequences of four Ixodes species expands understanding of tick evolution

BACKGROUND

Ticks, hematophagous Acari, pose a significant threat by transmitting various pathogens to their vertebrate hosts during feeding. Despite advances in tick genomics, high-quality genomes were lacking until recently, particularly in the genus Ixodes, which includes the main vectors of Lyme disease.

RESULTS

Here, we present the genome sequences of four tick species, derived from a single female individual, with a particular focus on the European species Ixodes ricinus, achieving a chromosome-level assembly. Additionally, draft assemblies were generated for the three other Ixodes species, I. persulcatus, I. pacificus, and I. hexagonus. The quality of the four genomes and extensive annotation of several important gene families have allowed us to study the evolution of gene repertoires at the level of the genus Ixodes and of the tick group. We have determined gene families that have undergone major amplifications during the evolution of ticks, while an expression atlas obtained for I. ricinus reveals striking patterns of specialization both between and within gene families. Notably, several gene family amplifications are associated with a proliferation of single-exon genes-most strikingly for fatty acid elongases and sulfotransferases.

CONCLUSIONS

The integration of our data with existing genomes establishes a solid framework for the study of gene evolution, improving our understanding of tick biology. In addition, our work lays the foundations for applied research and innovative control targeting these organisms.

Vibrio cidicii genomes recovered from Baltic Sea samples in Denmark

We report the genomic characteristics of the human pathogen Vibrio cidicii isolated from seawater and green algae in the Baltic Sea. Initially misidentified as Vibrio vulnificus through culture and MALDI-TOF, whole-genome sequencing (WGS) confirmed them as V. cidicii, highlighting the importance of WGS analysis in accurate classification of emerging pathogens.

Exploring SNP filtering strategies: the influence of strict vs soft core

Phylogenetic analyses are crucial for understanding microbial evolution and infectious disease transmission. Bacterial phylogenies are often inferred from SNP alignments, with SNPs as the fundamental signal within these data. SNP alignments can be reduced to a 'strict core' by removing those sites that do not have data present in every sample. However, as sample size and genome diversity increase, a strict core can shrink markedly, discarding potentially informative data. Here, we propose and provide evidence to support the use of a 'soft core' that tolerates some missing data, preserving more information for phylogenetic analysis. Using large datasets of Neisseria gonorrhoeae and Salmonella enterica serovar Typhi, we assess different core thresholds. Our results show that strict cores can drastically reduce informative sites compared to soft cores. In a 10 000-genome alignment of Salmonella enterica serovar Typhi, a 95% soft core yielded ten times more informative sites than a 100% strict core. Similar patterns were observed in N. gonorrhoeae. We further evaluated the accuracy of phylogenies built from strict- and soft-core alignments using datasets with strong temporal signals. Soft-core alignments generally outperformed strict cores in producing trees displaying clock-like behaviour; for instance, the N. gonorrhoeae 95% soft-core phylogeny had a root-to-tip regression R 2 of 0.50 compared to 0.21 for the strict-core phylogeny. This study suggests that soft-core strategies are preferable for large, diverse microbial datasets. To facilitate this, we developed Core-SNP-filter (https://github.com/rrwick/Core-SNP-filter), an open-source software tool for generating soft-core alignments from whole-genome alignments based on user-defined thresholds.

Viral niche-partitioning: comparative genomics of giant viruses across environmental gradients in a high Arctic freshwater-saltwater lake

Giant viruses (GVs; Nucleocytoviricota) impact the biology and ecology of a wide range of eukaryotic hosts, with implications for global biogeochemical cycles. Here, we investigated GV niche separation in highly stratified Lake A at the northern coast of Ellesmere Island, Nunavut, Canada. This lake is composed of a layer of ice-covered freshwater that overlies saltwater derived from the ancient Arctic Ocean, and it therefore provides a broad gradient of environmental conditions and ecological habitats, each with a distinct protist community and rich assemblages of associated GVs. The upper layer (mixolimnion) had measurable light and oxygen, and contained diverse GVs linked to photosynthetic protists, indicating adaptation to surface biotic and abiotic conditions. In contrast, the saline lower layer (monimolimnion), lacking oxygen and light, hosted GVs associated with predicted heterotrophic protists, some of which are known for a predatory lifestyle, and with several viral genes suggesting adaptation to deep-water anaerobic conditions. Our observations underscore the coupling between physical and chemical gradients, microeukaryotes and their associated GVs in Lake A, and provide insight into the potential for GVs to directly and indirectly impact host metabolism. There were similarities between the genetic composition of GVs and the metabolic processes of their potential hosts, implying co-evolution and niche-adaptation within the lake habitats. Notably, we found a greater presence of viral rhodopsins in deeper water layers, suggesting an evolutionary relationship with potential hosts capable of supplementing their energetic needs to thrive in low energy, anoxic conditions.

Group A streptococcal infections in Alberta, Canada 2018-2023

Group A streptococcal or Streptococcus pyogenes infections have been increasing post-COVID-19 pandemic. We describe the epidemiology of S. pyogenes pharyngitis and invasive disease in Alberta, Canada 2018-2023. Positive pharyngitis specimens were identified from throat swabs collected from pharyngitis patients. Invasive S. pyogenes was defined as the isolation of S. pyogenes from a normally sterile site or severe skin infection. S. pyogenes isolates were emm typed. Pharyngitis and invasive disease displayed seasonal trends preceding the COVID-19 pandemic followed by a sharp decrease during COVID-19 intervention measures. After the lifting of interventions, rates of pharyngitis and invasive disease rose. There were 182 983 positive pharyngitis specimens between 2018 and 2023 for a positivity rate of 17.6%. The highest rates occurred in the 0-9 age group in 2023 (41.5%). Invasive disease increased in 2022-2023 driven by emm1 and 12 types. M1UK strain was the most frequent M1 type associated with invasive disease (59% of M1 isolates sequenced). Notably, out of 182 983 pharyngitis cases, there were 111 cases of invasive S. pyogenes detected for an invasive disease rate of 0.06%. This descriptive epidemiology of S. pyogenes pharyngitis and invasive S. pyogenes disease highlights the rapid increase in cases of S. pyogenes occurring in western Canada and illustrates the critical need for a vaccine.

SARS-CoV-2 introductions to the island of Ireland: a phylogenetic and geospatiotemporal study of infection dynamics

BACKGROUND

Ireland's COVID-19 response combined extensive SARS-CoV-2 testing to estimate incidence, with whole genome sequencing (WGS) for genome surveillance. As an island with two political jurisdictions-Northern Ireland (NI) and Republic of Ireland (RoI)-and access to detailed passenger travel data, Ireland provides a unique setting to study virus introductions and evaluate public health measures. Using a substantial Irish genomic dataset alongside global data from GISAID, this study aimed to trace the introduction and spread of SARS-CoV-2 across the island.

METHODS

We recursively searched for 29,518 SARS-CoV-2 genome sequences collected in Ireland from March 2020 to June 2022 within the global SARS-CoV-2 phylogenetic tree and identified clusters based on shared last common non-Irish ancestors. A maximum parsimony approach was used to assign a likely country of origin to each cluster. The geographic locations and collection dates of the samples in each introduction cluster were used to map the spread of the virus across Ireland. Downsampling was used to model the impact of varying levels of sequencing and normalisation for population permitted comparison between jurisdictions.

RESULTS

Six periods spanning the early introductions and the emergence of Alpha, Delta, and Omicron variants were studied in detail. Among 4439 SARS-CoV-2 introductions to Ireland, 2535 originated in England, with additional cases largely from the rest of Great Britain, United States of America, and Northwestern Europe. Introduction clusters ranged in size from a single to thousands of cases. Introductions were concentrated in the densely populated Dublin and Belfast areas, with many clusters spreading islandwide. Genetic phylogeny was able to effectively trace localised transmission patterns. Introduction rates were similar in NI and RoI for most variants, except for Delta, which was more frequently introduced to NI.

CONCLUSIONS

Tracking individual introduction events enables detailed modelling of virus spread patterns and clearer assessment of the effectiveness of control measures. Stricter travel restrictions in RoI likely reduced Delta introductions but not infection rates, which were similar across jurisdictions. Local and global sequencing levels influence the information available from phylogenomic analyses and we describe an approach to assess the ability of a chosen WGS level to detect virus introductions.

Genome-wide phylogenetic analysis and expansion of gene families involved in detoxification in Smittia aterrima (Meigen)and Smittia pratorum (Goetghebuer) (Diptera, Chironomidae)

Smittia aterrima (Meigen, 1818) and Smittia pratorum (Goetghebuer, 1927) are important indicator insects for aquatic environments, showing extensive tolerance to the environment. However, the genome-wide phylogenetic relationships and characteristics of the detoxification mechanisms in S. aterrima and S. pratorum remain unclear. Based on the genomes of the two species obtained in our preliminary studies and nine genomes from the NCBI database, we found that chironomids diverged from other mosquitoes approximately 200 million years ago (MYA), and S. aterrima and S. pratorum diverged about 30 MYA according to phylogenetic analysis. Gene family evolution analysis showed significant expansion of 43 and 15 gene families in S. aterrima and S. pratorum, respectively, particularly those related to detoxification pathways. Positive selection analysis reveals that genes under positive selection are crucial for promoting environmental adaptation. Additionally, the detoxification-associated gene families including Cytochrome P450 (CYP), Glutathione S-transferases (GST), ATP-binding cassette (ABC), carboxylesterase (CCE), and UDP-glucuronosyltransferase (UGT) were annotated. Our analysis results show that these five detoxification gene families have significantly expanded in the chironomid genomes. This study highlights the genome evolution of chironomids and their responses to mechanisms of tolerance to environmental challenges.

Evolution of sequence, structural and functional diversity of the ubiquitous DNA/RNA-binding Alba domain

The DNA/RNA-binding Alba domain is prevalent across all kingdoms of life. First discovered in archaea, this protein domain has evolved from RNA- to DNA-binding, with a concomitant expansion in the range of cellular processes that it regulates. Despite its widespread presence, the full extent of its sequence, structural, and functional diversity remains unexplored. In this study, we employed iterative searches in PSI-BLAST to identify 15,161 unique Alba domain-containing proteins from the NCBI non-redundant protein database. Sequence similarity network (SSN) analysis clustered them into 13 distinct subgroups, including the archaeal Alba and eukaryotic Rpp20/Pop7 and Rpp25/Pop6 groups, as well as novel fungal and Plasmodium-specific Albas. Sequence and structural conservation analysis of the subgroups indicated high preservation of the dimer interface, with Alba domains from unicellular eukaryotes notably exhibiting structural deviations towards their C-terminal end. Finally, phylogenetic analysis, while supporting SSN clustering, revealed the evolutionary branchpoint at which the eukaryotic Rpp20- and Rpp25-like clades emerged from archaeal Albas, and the subsequent taxonomic lineage-based divergence within each clade. Taken together, this comprehensive analysis enhances our understanding of the evolutionary history of Alba domain-containing proteins across diverse organisms.

Fpa (YlaN) is an iron(II) binding protein that functions to relieve Fur-mediated repression of gene expression in Staphylococcus aureus

Iron (Fe) is a trace nutrient required by nearly all organisms. As a result of the demand for Fe and the toxicity of non-chelated cytosolic ionic Fe, regulatory systems have evolved to tightly balance Fe acquisition and usage while limiting overload. In most bacteria, including the mammalian pathogen Staphylococcus aureus, the ferric uptake regulator (Fur) is the primary transcriptional regulator controlling the transcription of genes that code for Fe uptake and utilization proteins. Fpa (formerly YlaN) was demonstrated to be essential in Bacillus subtilis unless excess Fe is added to the growth medium, suggesting a role in Fe homeostasis. Here, we demonstrate that Fpa is essential in S. aureus upon Fe deprivation. Null fur alleles bypassed the essentiality of Fpa. The absence of Fpa abolished the derepression of Fur-regulated genes during Fe limitation. Bioinformatic analyses suggest that fpa was recruited to Gram-positive bacteria and, once acquired, was maintained in the genome as it co-evolved with Fur. Consistent with a role for Fpa in alleviating Fur-dependent repression, Fpa and Fur interacted in vivo, and Fpa decreased the DNA-binding ability of Fur in vitro. Fpa bound Fe(II) in vitro using oxygen or nitrogen ligands with an association constant that is consistent with a physiological role in Fe homeostasis. These findings have led to a model wherein Fpa is an Fe(II) binding protein that influences Fur-dependent regulation through direct interaction.IMPORTANCEIron (Fe) is an essential nutrient for nearly all organisms. If Fe homeostasis is not maintained, Fe may accumulate in the cytosol, which can be toxic. Questions remain about how cells efficiently balance Fe uptake and usage to prevent overload. Iron uptake and proper metalation of proteins are essential processes in the mammalian bacterial pathogen Staphylococcus aureus. Understanding the gene products involved in the genetic regulation of Fe uptake and usage and the physiological adaptations that S. aureus uses to survive in Fe-depleted conditions provides insight into pathogenesis. Herein, we demonstrate that the DNA-binding activity of the ferric uptake regulator transcriptional repressor is alleviated under Fe limitation, but uniquely, in S. aureus, alleviation requires the presence of Fpa.

Gut Microbial Utilization of the Alternative Sweetener, D-Allulose, via AlsE

D-allulose, a rare sugar with emerging potential as a low-calorie sweetener, has garnered attention as an alternative to other commercially available alternative sweeteners, such as sugar alcohols, which often cause severe gastrointestinal discomfort. D-allulose-6-phosphate 3-epimerase (AlsE) is a prokaryotic enzyme that converts D-allulose-6-phosphate into D-fructose-6-phopshate, enabling its use as a carbon source. However, the taxonomic breadth of AlsE across gut bacteria remains poorly understood, hindering insights into the utilization of D-allulose by microbial communities. In this study, we provide experimental evidence showing that Clostridium innocuum is capable of D-allulose metabolism via a homologous AlsE. A bioinformatics search of 85,202 bacterial genomes identified 116 bacterial species with AlsE homologs, suggesting a limited distribution of AlsE in bacteria. Additionally, Escherichia coli contains a copy of alsE , but it does not grow on D-allulose as a sole carbon source unless alsE is heterologously expressed. A metagenomic analysis revealed that 15.8% of 3,079 adult healthy human metagenomic samples that we analyzed contained alsE , suggesting a limited prevalence of the enzyme in the gut microbiome. These results suggest that the gut microbiome has limited capacity to metabolize D-allulose via alsE , supporting its use as an alternative sweetener with minimal impact on microbial composition and gastrointestinal symptoms. This finding also enables personalized nutrition, allowing diabetic individuals to assess their gut microbiota for alsE , and manage glycemic response while reducing gastrointestinal distress.

Eukaryotic viruses encode the ribosomal protein eL40

Viruses in the phylum Nucleocytoviricota are large, complex and have an exceptionally diverse metabolic repertoire. Some encode hundreds of products involved in the translation of mRNA into protein, but none was known to encode any of the proteins in ribosomes, the central engines of translation. With the discovery of the eL40 gene in FloV-SA2, we report the first example of a eukaryotic virus encoding a ribosomal protein and show that this gene is also present and expressed in other uncultivated marine giant viruses. FloV-SA2 also encodes a "group II" viral rhodopsin, a viral light-activated protein of unknown function previously only reported in metagenomes. FloV-SA2 is thus a valuable model system for investigating new mechanisms by which viruses manipulate eukaryotic cell metabolism.

Genomic and biogeographic characterisation of the novel prasinovirus Mantoniella tinhauana virus 1

Mamiellophyceae are a ubiquitous class of unicellular green algae in the global ocean. Their ecological importance is highlighted in studies focused on the prominent genera Micromonas, Ostreococcus, and Bathycoccus. Mamiellophyceae are susceptible to prasinoviruses, double-stranded DNA viruses belonging to the nucleocytoplasmic large DNA virus group. Our study represents the first isolation of a prasinovirus in the South China Sea and the only one to infect the globally distributed genus Mantoniella. We conducted a comparative analysis with previously identified viral relatives, encompassing morphological characteristics, host specificity, marker-based phylogenetic placement, and whole-genome sequence comparisons. Although it shares morphological and genetic similarities with established prasinoviruses, this novel virus showed distinct genetic traits, confining its infection to the species Mantoniella tinhauana. We also explored the global biogeography of this prasinovirus and its host by mapping metagenomic data and analysing their relationship with various environmental parameters. Our results demonstrate a pronounced link between the virus and its host, both found predominantly in higher latitudes in the surface ocean. By gaining an increased understanding of the relationships between viruses, hosts, and environments, researchers can better make predictions and potentially implement measures to mitigate the consequences of climate change on oceanic processes.

Next-generation bioinformatics: An ultrafast and user-friendly tool for phylogenomic data exploration

With increasingly large genomic datasets, even routine bioinformatic tasks can be arduous, computationally demanding, and time-consuming. Additionally, most bioinformatic programs do not have a graphical user interface (GUI) and thus, require users to be minimally competent in command-line. These impediments present significant economic and technological barriers for practitioners who do not have access to advanced computational resources and support. In this issue of Molecular Ecology Resources, Handika and Esselstyn (2024) developed an ultrafast and memory-efficient bioinformatic tool, SEGUL, that performs common manipulations and calculations of summary statistics on phylogenomic datasets. SEGUL has two main features that distinguish it from other bioinformatic programs: (1) it is based on the recently emergent, high-performance programming language Rust, and (2) it has a GUI written using Flutter, a cross-platform programming framework that also supports mobile operating systems (mobile iOS, iPadOS and Android). By leveraging and combining the power of Rust and Flutter, SEGUL achieves significantly faster computation times and lower memory usage across different platforms and CPU architectures compared to similar programs. The unique and innovative features of SEGUL pave the way for a new era of bioinformatics that can be more energy-efficient, inclusive, and available to a broader swathe of users.

SEGUL: Ultrafast, memory-efficient and mobile-friendly software for manipulating and summarizing phylogenomic datasets

Phylogenetic studies now routinely require manipulating and summarizing thousands of data files. For most of these tasks, currently available software requires considerable computing resources and substantial knowledge of command-line applications. We develop an ultrafast and memory-efficient software, SEGUL, that performs common phylogenomic dataset manipulations and calculates statistics summarizing essential data features. Our software is available as standalone command-line interface (CLI) and graphical user interface (GUI) applications, and as a library for Rust, R and Python, with possible support of other languages. The CLI and library versions run native on Windows, Linux and macOS, including Apple ARM Macs. The GUI version extends support to include mobile iOS, iPadOS and Android operating systems. SEGUL leverages the high performance of the Rust programming language to offer fast execution times and low memory footprints regardless of dataset size and platform choice. The inclusion of a GUI minimizes bioinformatics barriers to phylogenomics while SEGUL's efficiency reduces economic barriers by allowing analysis on inexpensive hardware. Our support for mobile operating systems further enables teaching phylogenomics where access to computing power is limited.

Ancestral Sequence Reconstruction as a Tool to Detect and Study De Novo Gene Emergence

New protein-coding genes can evolve from previously noncoding genomic regions through a process known as de novo gene emergence. Evidence suggests that this process has likely occurred throughout evolution and across the tree of life. Yet, confidently identifying de novo emerged genes remains challenging. Ancestral sequence reconstruction is a promising approach for inferring whether a gene has emerged de novo or not, as it allows us to inspect whether a given genomic locus ancestrally harbored protein-coding capacity. However, the use of ancestral sequence reconstruction in the context of de novo emergence is still in its infancy and its capabilities, limitations, and overall potential are largely unknown. Notably, it is difficult to formally evaluate the protein-coding capacity of ancestral sequences, particularly when new gene candidates are short. How well-suited is ancestral sequence reconstruction as a tool for the detection and study of de novo genes? Here, we address this question by designing an ancestral sequence reconstruction workflow incorporating different tools and sets of parameters and by introducing a formal criterion that allows to estimate, within a desired level of confidence, when protein-coding capacity originated at a particular locus. Applying this workflow on ∼2,600 short, annotated budding yeast genes (<1,000 nucleotides), we found that ancestral sequence reconstruction robustly predicts an ancient origin for the most widely conserved genes, which constitute "easy" cases. For less robust cases, we calculated a randomization-based empirical P-value estimating whether the observed conservation between the extant and ancestral reading frame could be attributed to chance. This formal criterion allowed us to pinpoint a branch of origin for most of the less robust cases, identifying 49 genes that can unequivocally be considered de novo originated since the split of the Saccharomyces genus, including 37 Saccharomyces cerevisiae-specific genes. We find that for the remaining equivocal cases we cannot rule out different evolutionary scenarios including rapid evolution, multiple gene losses, or a recent de novo origin. Overall, our findings suggest that ancestral sequence reconstruction is a valuable tool to study de novo gene emergence but should be applied with caution and awareness of its limitations.

Chromosomal inversions harbour excess mutational load in the coral, Acropora kenti, on the Great Barrier Reef

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future.

Large-scale investigation of species-specific orphan genes in the human gut microbiome elucidates their evolutionary origins

Species-specific genes, also known as orphans, are ubiquitous across life's domains. In prokaryotes, species-specific orphan genes (SSOGs) are mostly thought to originate in external elements such as viruses followed by horizontal gene transfer, whereas the scenario of native origination, through rapid divergence or de novo, is mostly dismissed. However, quantitative evidence supporting either scenario is lacking. Here, we systematically analyzed genomes from 4644 human gut microbiome species and identified more than 600,000 unique SSOGs, representing an average of 2.6% of a given species' pangenome. These sequences are mostly rare within each species yet show signs of purifying selection. Overall, SSOGs use optimal codons less frequently, and their proteins are more disordered than those of conserved genes (i.e., non-SSOGs). Importantly, across species, the GC content of SSOGs closely matches that of conserved ones. In contrast, the ∼5% of SSOGs that share similarity to known viral sequences have distinct characteristics, including lower GC content. Thus, SSOGs with similarity to viruses differ from the remaining SSOGs, contrasting an external origination scenario for most of them. By examining the orthologous genomic region in closely related species, we show that a small subset of SSOGs likely evolved natively de novo and find that these genes also differ in their properties from the remaining SSOGs. Our results challenge the notion that external elements are the dominant source of prokaryotic genetic novelty and will enable future studies into the biological role and relevance of species-specific genes in the human gut.

NeMu: a comprehensive pipeline for accurate reconstruction of neutral mutation spectra from evolutionary data

The recognized importance of mutational spectra in molecular evolution is yet to be fully exploited beyond human cancer studies and model organisms. The wealth of intraspecific polymorphism data in the GenBank repository, covering a broad spectrum of genes and species, presents an untapped opportunity for detailed mutational spectrum analysis. Existing methods fall short by ignoring intermediate substitutions on the inner branches of phylogenetic trees and lacking the capability for cross-species mutational comparisons. To address these challenges, we present the NeMu pipeline, available at https://nemu-pipeline.com, a tool grounded in phylogenetic principles designed to provide comprehensive and scalable analysis of mutational spectra. Utilizing extensive sequence data from numerous available genome projects, NeMu rapidly and accurately reconstructs the neutral mutational spectrum. This tool, facilitating the reconstruction of gene- and species-specific mutational spectra, contributes to a deeper understanding of evolutionary mechanisms across the broad spectrum of known species.

Unique trajectory of gene family evolution from genomic analysis of nearly all known species in an ancient yeast lineage

Gene gains and losses are a major driver of genome evolution; their precise characterization can provide insights into the origin and diversification of major lineages. Here, we examined gene family evolution of 1,154 genomes from nearly all known species in the medically and technologically important yeast subphylum Saccharomycotina. We found that yeast gene family and genome evolution are distinct from plants, animals, and filamentous ascomycetes and are characterized by small genome sizes and smaller gene numbers but larger gene family sizes. Faster-evolving lineages (FELs) in yeasts experienced significantly higher rates of gene losses-commensurate with a narrowing of metabolic niche breadth-but higher speciation rates than their slower-evolving sister lineages (SELs). Gene families most often lost are those involved in mRNA splicing, carbohydrate metabolism, and cell division and are likely associated with intron loss, metabolic breadth, and non-canonical cell cycle processes. Our results highlight the significant role of gene family contractions in the evolution of yeast metabolism, genome function, and speciation, and suggest that gene family evolutionary trajectories have differed markedly across major eukaryotic lineages.

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.

Deep Learning and Likelihood Approaches for Viral Phylogeography Converge on the Same Answers Whether the Inference Model Is Right or Wrong

Analysis of phylogenetic trees has become an essential tool in epidemiology. Likelihood-based methods fit models to phylogenies to draw inferences about the phylodynamics and history of viral transmission. However, these methods are often computationally expensive, which limits the complexity and realism of phylodynamic models and makes them ill-suited for informing policy decisions in real-time during rapidly developing outbreaks. Likelihood-free methods using deep learning are pushing the boundaries of inference beyond these constraints. In this paper, we extend, compare, and contrast a recently developed deep learning method for likelihood-free inference from trees. We trained multiple deep neural networks using phylogenies from simulated outbreaks that spread among 5 locations and found they achieve close to the same levels of accuracy as Bayesian inference under the true simulation model. We compared robustness to model misspecification of a trained neural network to that of a Bayesian method. We found that both models had comparable performance, converging on similar biases. We also implemented a method of uncertainty quantification called conformalized quantile regression that we demonstrate has similar patterns of sensitivity to model misspecification as Bayesian highest posterior density (HPD) and greatly overlap with HPDs, but have lower precision (more conservative). Finally, we trained and tested a neural network against phylogeographic data from a recent study of the SARS-Cov-2 pandemic in Europe and obtained similar estimates of region-specific epidemiological parameters and the location of the common ancestor in Europe. Along with being as accurate and robust as likelihood-based methods, our trained neural networks are on average over 3 orders of magnitude faster after training. Our results support the notion that neural networks can be trained with simulated data to accurately mimic the good and bad statistical properties of the likelihood functions of generative phylogenetic models.

Genomic characterisation and ecological distribution of Mantoniella tinhauana: a novel Mamiellophycean green alga from the Western Pacific

Mamiellophyceae are dominant marine algae in much of the ocean, the most prevalent genera belonging to the order Mamiellales: Micromonas, Ostreococcus and Bathycoccus, whose genetics and global distributions have been extensively studied. Conversely, the genus Mantoniella, despite its potential ecological importance, remains relatively under-characterised. In this study, we isolated and characterised a novel species of Mamiellophyceae, Mantoniella tinhauana, from subtropical coastal waters in the South China Sea. Morphologically, it resembles other Mantoniella species; however, a comparative analysis of the 18S and ITS2 marker genes revealed its genetic distinctiveness. Furthermore, we sequenced and assembled the first genome of Mantoniella tinhauana, uncovering significant differences from previously studied Mamiellophyceae species. Notably, the genome lacked any detectable outlier chromosomes and exhibited numerous unique orthogroups. We explored gene groups associated with meiosis, scale and flagella formation, shedding light on species divergence, yet further investigation is warranted. To elucidate the biogeography of Mantoniella tinhauana, we conducted a comprehensive analysis using global metagenomic read mapping to the newly sequenced genome. Our findings indicate this species exhibits a cosmopolitan distribution with a low-level prevalence worldwide. Understanding the intricate dynamics between Mamiellophyceae and the environment is crucial for comprehending their impact on the ocean ecosystem and accurately predicting their response to forthcoming environmental changes.

Frequent nonhomologous replacement of replicative helicase loaders by viruses in Vibrionaceae

Several microbial genomes lack textbook-defined essential genes. If an essential gene is absent from a genome, then an evolutionarily independent gene of unknown function complements its function. Here, we identified frequent nonhomologous replacement of an essential component of DNA replication initiation, a replicative helicase loader gene, in Vibrionaceae. Our analysis of Vibrionaceae genomes revealed two genes with unknown function, named vdhL1 and vdhL2, that were substantially enriched in genomes without the known helicase-loader genes. These genes showed no sequence similarities to genes with known function but encoded proteins structurally similar with a viral helicase loader. Analyses of genomic syntenies and coevolution with helicase genes suggested that vdhL1/2 encodes a helicase loader. The in vitro assay showed that Vibrio harveyi VdhL1 and Vibrio ezurae VdhL2 promote the helicase activity of DnaB. Furthermore, molecular phylogenetics suggested that vdhL1/2 were derived from phages and replaced an intrinsic helicase loader gene of Vibrionaceae over 20 times. This high replacement frequency implies the host's advantage in acquiring a viral helicase loader gene.

Adaptations to nitrogen availability drive ecological divergence of chemosynthetic symbionts

Bacterial symbionts, with their shorter generation times and capacity for horizontal gene transfer (HGT), play a critical role in allowing marine organisms to cope with environmental change. The closure of the Isthmus of Panama created distinct environmental conditions in the Tropical Eastern Pacific (TEP) and Caribbean, offering a "natural experiment" for studying how closely related animals evolve and adapt under environmental change. However, the role of bacterial symbionts in this process is often overlooked. We sequenced the genomes of endosymbiotic bacteria in two sets of sister species of chemosymbiotic bivalves from the genera Codakia and Ctena (family Lucinidae) collected on either side of the Isthmus, to investigate how differing environmental conditions have influenced the selection of symbionts and their metabolic capabilities. The lucinid sister species hosted different Candidatus Thiodiazotropha symbionts and only those from the Caribbean had the genetic potential for nitrogen fixation, while those from the TEP did not. Interestingly, this nitrogen-fixing ability did not correspond to symbiont phylogeny, suggesting convergent evolution of nitrogen fixation potential under nutrient-poor conditions. Reconstructing the evolutionary history of the nifHDKT operon by including other lucinid symbiont genomes from around the world further revealed that the last common ancestor (LCA) of Ca. Thiodiazotropha lacked nif genes, and populations in oligotrophic habitats later re-acquired the nif operon through HGT from the Sedimenticola symbiont lineage. Our study suggests that HGT of the nif operon has facilitated niche diversification of the globally distributed Ca. Thiodiazotropha endolucinida species clade. It highlights the importance of nitrogen availability in driving the ecological diversification of chemosynthetic symbiont species and the role that bacterial symbionts may play in the adaptation of marine organisms to changing environmental conditions.

Accurate Detection of Convergent Mutations in Large Protein Alignments With ConDor

Evolutionary convergences are observed at all levels, from phenotype to DNA and protein sequences, and changes at these different levels tend to be correlated. Notably, convergent mutations can lead to convergent changes in phenotype, such as changes in metabolism, drug resistance, and other adaptations to changing environments. We propose a two-component approach to detect mutations subject to convergent evolution in protein alignments. The "Emergence" component selects mutations that emerge more often than expected, while the "Correlation" component selects mutations that correlate with the convergent phenotype under study. With regard to Emergence, a phylogeny deduced from the alignment is provided by the user and is used to simulate the evolution of each alignment position. These simulations allow us to estimate the expected number of mutations in a neutral model, which is compared to the observed number of mutations in the data studied. In Correlation, a comparative phylogenetic approach, is used to measure whether the presence of each of the observed mutations is correlated with the convergent phenotype. Each component can be used on its own, for example Emergence when no phenotype is available. Our method is implemented in a standalone workflow and a webserver, called ConDor. We evaluate the properties of ConDor using simulated data, and we apply it to three real datasets: sedge PEPC proteins, HIV reverse transcriptase, and fish rhodopsin. The results show that the two components of ConDor complement each other, with an overall accuracy that compares favorably to other available tools, especially on large datasets.

Evolutionary trajectory of pattern recognition receptors in plants

Cell-surface receptors play pivotal roles in many biological processes, including immunity, development, and reproduction, across diverse organisms. How cell-surface receptors evolve to become specialised in different biological processes remains elusive. To shed light on the immune-specificity of cell-surface receptors, we analyzed more than 200,000 genes encoding cell-surface receptors from 350 genomes and traced the evolutionary origin of immune-specific leucine-rich repeat receptor-like proteins (LRR-RLPs) in plants. Surprisingly, we discovered that the motifs crucial for co-receptor interaction in LRR-RLPs are closely related to those of the LRR-receptor-like kinase (RLK) subgroup Xb, which perceives phytohormones and primarily governs growth and development. Functional characterisation further reveals that LRR-RLPs initiate immune responses through their juxtamembrane and transmembrane regions, while LRR-RLK-Xb members regulate development through their cytosolic kinase domains. Our data suggest that the cell-surface receptors involved in immunity and development share a common origin. After diversification, their ectodomains, juxtamembrane, transmembrane, and cytosolic regions have either diversified or stabilised to recognise diverse ligands and activate differential downstream responses. Our work reveals a mechanism by which plants evolve to perceive diverse signals to activate the appropriate responses in a rapidly changing environment.

FicD genes in invertebrates: A tale of transposons, pathogenic and integrated viruses

Many gene families are shared across the tree of life between distantly related species because of horizontal gene transfers (HGTs). However, the frequency of HGTs varies strongly between gene families and biotic realms suggesting differential selection pressures and functional bias. One gene family with a wide distribution are FIC-domain containing enzymes (FicDs). FicDs catalyze AMPylation, a post-translational protein modification consisting in the addition of adenosine monophosphate to accessible residues of target proteins. Beside the well-known conservation of FicDs in deuterostomes, we report the presence of a conserved FicD gene ortholog in a large number of protostomes and microbial eukaryotes. We also reported additional FicD gene copies in the genomes of some rotifers, parasitic worms and bivalves. A few dsDNA viruses of these invertebrates, including White spot syndrome virus, Cherax quadricarinatus iridovirus, Ostreid herpesvirus-1 and the beetle nudivirus, carry copies of FicDs, with phylogenetic analysis suggesting a common origin of these FicD copies and the duplicated FicDs of their invertebrate hosts. HGTs and gene duplications possibly mediated by endogenous viruses or genetic mobile elements seem to have contributed to the transfer of AMPylation ability from bacteria and eukaryotes to pathogenic viruses, where this pathway could have been hijacked to promote viral infection.

Unveiling Prasinovirus diversity and host specificity through targeted enrichment in the South China Sea

Unicellular green picophytoplankton from the Mamiellales order are pervasive in marine ecosystems and susceptible to infections by prasinoviruses, large double-stranded DNA viruses within the Nucleocytoviricota phylum. We developed a double-stranded DNA virus enrichment and shotgun sequencing method, and successfully assembled 80 prasinovirus genomes from 43 samples in the South China Sea. Our research delivered the first direct estimation of 94% accuracy in correlating genome similarity to host range. Stirkingly, our analyses uncovered unexpected host-switching across diverse algal lineages, challenging the existing paradigms of host-virus co-speciation and revealing the dynamic nature of viral evolution. We also detected six instances of horizontal gene transfer between prasinoviruses and their hosts, including a novel alternative oxidase. Additionally, diversifying selection on a major capsid protein suggests an ongoing co-evolutionary arms race. These insights not only expand our understanding of prasinovirus genomic diversity but also highlight the intricate evolutionary mechanisms driving their ecological success and shaping broader virus-host interactions in marine environments.

Robustness of Felsenstein's Versus Transfer Bootstrap Supports With Respect to Taxon Sampling

The bootstrap method is based on resampling sequence alignments and re-estimating trees. Felsenstein's bootstrap proportions (FBP) are the most common approach to assess the reliability and robustness of sequence-based phylogenies. However, when increasing taxon sampling (i.e., the number of sequences) to hundreds or thousands of taxa, FBP tend to return low support for deep branches. The transfer bootstrap expectation (TBE) has been recently suggested as an alternative to FBP. TBE is measured using a continuous transfer index in [0,1] for each bootstrap tree, instead of the binary {0,1} index used in FBP to measure the presence/absence of the branch of interest. TBE has been shown to yield higher and more informative supports while inducing a very low number of falsely supported branches. Nonetheless, it has been argued that TBE must be used with care due to sampling issues, especially in datasets with a high number of closely related taxa. In this study, we conduct multiple experiments by varying taxon sampling and comparing FBP and TBE support values on different phylogenetic depths, using empirical datasets. Our results show that the main critique of TBE stands in extreme cases with shallow branches and highly unbalanced sampling among clades, but that TBE is still robust in most cases, while FBP is inescapably negatively impacted by high taxon sampling. We suggest guidelines and good practices in TBE (and FBP) computing and interpretation.

Genomic surveillance reveals dynamic shifts in the connectivity of COVID-19 epidemics

The maturation of genomic surveillance in the past decade has enabled tracking of the emergence and spread of epidemics at an unprecedented level. During the COVID-19 pandemic, for example, genomic data revealed that local epidemics varied considerably in the frequency of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage importation and persistence, likely due to a combination of COVID-19 restrictions and changing connectivity. Here, we show that local COVID-19 epidemics are driven by regional transmission, including across international boundaries, but can become increasingly connected to distant locations following the relaxation of public health interventions. By integrating genomic, mobility, and epidemiological data, we find abundant transmission occurring between both adjacent and distant locations, supported by dynamic mobility patterns. We find that changing connectivity significantly influences local COVID-19 incidence. Our findings demonstrate a complex meaning of "local" when investigating connected epidemics and emphasize the importance of collaborative interventions for pandemic prevention and mitigation.

APOBEC3F Is a Mutational Driver of the Human Monkeypox Virus Identified in the 2022 Outbreak

BACKGROUND

On May 6, 2022, a powerful outbreak of monkeypox virus (MPXV) had been reported outside of Africa, with many continuing new cases being reported around the world. Analysis of mutations among the 2 different lineages present in the 2021 and 2022 outbreaks revealed the presence of G->A mutations occurring in the 5'GpA context, indicative of APOBEC3 cytidine deaminase activity.

METHODS

By using a sensitive polymerase chain reaction (differential DNA denaturation PCR) method allowing differential amplification of AT-rich DNA, we analyzed the level of APOBEC3-induced MPXV editing in infected cells and in patients.

RESULTS

We demonstrate that G->A hypermutated MPXV genomes can be recovered experimentally from APOBEC3 transfection followed by MPXV infection. Here, among the 7 human APOBEC3 cytidine deaminases (A3A-A3C, A3DE, A3F-A3H), only APOBEC3F was capable of extensively deaminating cytidine residues in MPXV genomes. Hyperedited genomes were also recovered in ∼42% of analyzed patients. Moreover, we demonstrate that substantial repair of these mutations occurs. Upon selection, corrected G->A mutations escaping drift loss contribute to the MPXV evolution observed in the current epidemic.

CONCLUSIONS

Stochastic or transient overexpression of the APOBEC3F gene exposes the MPXV genome to a broad spectrum of mutations that may be modeling the mutational landscape after multiple cycles of viral replication.

Bioinformatic investigation of discordant sequence data for SARS-CoV-2: insights for robust genomic analysis during pandemic surveillance

The COVID-19 pandemic has necessitated the rapid development and implementation of whole-genome sequencing (WGS) and bioinformatic methods for managing the pandemic. However, variability in methods and capabilities between laboratories has posed challenges in ensuring data accuracy. A national working group comprising 18 laboratory scientists and bioinformaticians from Australia and New Zealand was formed to improve data concordance across public health laboratories (PHLs). One effort, presented in this study, sought to understand the impact of the methodology on consensus genome concordance and interpretation. SARS-CoV-2 WGS proficiency testing programme (PTP) data were retrospectively obtained from the 2021 Royal College of Pathologists of Australasia Quality Assurance Programmes (RCPAQAP), which included 11 participating Australian laboratories. The submitted consensus genomes and reads from eight contrived specimens were investigated, focusing on discordant sequence data and findings were presented to the working group to inform best practices. Despite using a variety of laboratory and bioinformatic methods for SARS-CoV-2 WGS, participants largely produced concordant genomes. Two participants returned five discordant sites in a high-Cτ replicate, which could be resolved with reasonable bioinformatic quality thresholds. We noted ten discrepancies in genome assessment that arose from nucleotide heterogeneity at three different sites in three cell-culture-derived control specimens. While these sites were ultimately accurate after considering the participants' bioinformatic parameters, it presented an interesting challenge for developing standards to account for intrahost single nucleotide variation (iSNV). Observed differences had little to no impact on key surveillance metrics, lineage assignment and phylogenetic clustering, while genome coverage <90 % affected both. We recommend PHLs bioinformatically generate two consensus genomes with and without ambiguity thresholds for quality control and downstream analysis, respectively, and adhere to a minimum 90 % genome coverage threshold for inclusion in surveillance interpretations. We also suggest additional PTP assessment criteria, including primer efficiency, detection of iSNVs and minimum genome coverage of 90 %. This study underscores the importance of multidisciplinary national working groups in informing guidelines in real time for bioinformatic quality acceptance criteria. It demonstrates the potential for enhancing public health responses through improved data concordance and quality control in SARS-CoV-2 genomic analysis during pandemic surveillance.

BioConvert: a comprehensive format converter for life sciences

Bioinformatics is a field known for the numerous standards and formats that have been developed over the years. This plethora of formats, sometimes complementary, and often redundant, poses many challenges to bioinformatics data analysts. They constantly need to find the best tool to convert their data into the suitable format, which is often a complex, technical and time consuming task. Moreover, these small yet important tasks are often difficult to make reproducible. To overcome these difficulties, we initiated BioConvert, a collaborative project to facilitate the conversion of life science data from one format to another. BioConvert aggregates existing software within a single framework and complemented them with original code when needed. It provides a common interface to make the user experience more streamlined instead of having to learn tens of them. Currently, BioConvert supports about 50 formats and 100 direct conversions in areas such as alignment, sequencing, phylogeny, and variant calling. In addition to being useful for end-users, BioConvert can also be utilized by developers as a universal benchmarking framework for evaluating and comparing numerous conversion tools. Additionally, we provide a web server implementing an online user-friendly interface to BioConvert, hence allowing direct use for the community.

Genomic epidemiology offers high resolution estimates of serial intervals for COVID-19

Serial intervals - the time between symptom onset in infector and infectee - are a fundamental quantity in infectious disease control. However, their estimation requires knowledge of individuals' exposures, typically obtained through resource-intensive contact tracing efforts. We introduce an alternate framework using virus sequences to inform who infected whom and thereby estimate serial intervals. We apply our technique to SARS-CoV-2 sequences from case clusters in the first two COVID-19 waves in Victoria, Australia. We find that our approach offers high resolution, cluster-specific serial interval estimates that are comparable with those obtained from contact data, despite requiring no knowledge of who infected whom and relying on incompletely-sampled data. Compared to a published serial interval, cluster-specific serial intervals can vary estimates of the effective reproduction number by a factor of 2-3. We find that serial interval estimates in settings such as schools and meat processing/packing plants are shorter than those in healthcare facilities.

Weighted ASTRID: fast and accurate species trees from weighted internode distances

BACKGROUND

Species tree estimation is a basic step in many biological research projects, but is complicated by the fact that gene trees can differ from the species tree due to processes such as incomplete lineage sorting (ILS), gene duplication and loss (GDL), and horizontal gene transfer (HGT), which can cause different regions within the genome to have different evolutionary histories (i.e., "gene tree heterogeneity"). One approach to estimating species trees in the presence of gene tree heterogeneity resulting from ILS operates by computing trees on each genomic region (i.e., computing "gene trees") and then using these gene trees to define a matrix of average internode distances, where the internode distance in a tree T between two species x and y is the number of nodes in T between the leaves corresponding to x and y. Given such a matrix, a tree can then be computed using methods such as neighbor joining. Methods such as ASTRID and NJst (which use this basic approach) are provably statistically consistent, very fast (low degree polynomial time) and have had high accuracy under many conditions that makes them competitive with other popular species tree estimation methods. In this study, inspired by the very recent work of weighted ASTRAL, we present weighted ASTRID, a variant of ASTRID that takes the branch uncertainty on the gene trees into account in the internode distance.

RESULTS

Our experimental study evaluating weighted ASTRID typically shows improvements in accuracy compared to the original (unweighted) ASTRID, and shows competitive accuracy against weighted ASTRAL, the state of the art. Our re-implementation of ASTRID also improves the runtime, with marked improvements on large datasets.

CONCLUSIONS

Weighted ASTRID is a new and very fast method for species tree estimation that typically improves upon ASTRID and has comparable accuracy to weighted ASTRAL, while remaining much faster. Weighted ASTRID is available at https://github.com/RuneBlaze/internode .

Integrating full and partial genome sequences to decipher the global spread of canine rabies virus

Despite the rapid growth in viral genome sequencing, statistical methods face challenges in handling historical viral endemic diseases with large amounts of underutilized partial sequence data. We propose a phylogenetic pipeline that harnesses both full and partial viral genome sequences to investigate historical pathogen spread between countries. Its application to rabies virus (RABV) yields precise dating and confident estimates of its geographic dispersal. By using full genomes and partial sequences, we reduce both geographic and genetic biases that often hinder studies that focus on specific genes. Our pipeline reveals an emergence of the present canine-mediated RABV between years 1301 and 1403 and reveals regional introductions over a 700-year period. This geographic reconstruction enables us to locate episodes of human-mediated introductions of RABV and examine the role that European colonization played in its spread. Our approach enables phylogeographic analysis of large and genetically diverse data sets for many viral pathogens.

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

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

Enhanced neutralization escape to therapeutic monoclonal antibodies by SARS-CoV-2 omicron sub-lineages

The landscape of SARS-CoV-2 variants dramatically diversified with the simultaneous appearance of multiple subvariants originating from BA.2, BA.4, and BA.5 Omicron sub-lineages. They harbor a specific set of mutations in the spike that can make them more evasive to therapeutic monoclonal antibodies. In this study, we compared the neutralizing potential of monoclonal antibodies against the Omicron BA.2.75.2, BQ.1, BQ.1.1, and XBB variants, with a pre-Omicron Delta variant as a reference. Sotrovimab retains some activity against BA.2.75.2, BQ.1, and XBB as it did against BA.2/BA.5, but is less active against BQ.1.1. Within the Evusheld/AZD7442 cocktail, Cilgavimab lost all activity against all subvariants studied, resulting in loss of Evusheld activity. Finally, Bebtelovimab, while still active against BA.2.75, also lost all neutralizing activity against BQ.1, BQ.1.1, and XBB variants.

In-Depth Characterization of Full-Length Archived Viral Genomes after Nine Years of Posttreatment HIV Control

In the search for control of human immunodeficiency virus type 1 (HIV-1) infection without antiretroviral therapy, posttreatment controllers (PTCs) are models of HIV remission. To better understand their mechanisms of control, we characterized the HIV blood reservoirs of 8 PTCs (median of 9.4 years after treatment interruption) in comparison with those of 13 natural HIV infection controllers (HICs) (median of 18 years of infection) and with those of individuals receiving efficient antiretroviral therapy initiated during either primary HIV infection (PHIs; n = 8) or chronic HIV infection (CHIs; n = 6). This characterization was performed with single-genome amplification and deep sequencing. The proviral diversity, which reflects the history of past viral replication, was lower in the PTCs, PHIs, and aviremic HICs than in the blipper HICs and CHIs. The proportions of intact and defective proviruses among the proviral pool in PTCs were not significantly different from those of other groups. When looking at the quantities of proviruses per million peripheral blood mononuclear cells (PBMCs), they had similar amounts of intact proviruses as other groups but smaller amounts of defective proviruses than CHIs, suggesting a role of these forms in HIV pathogenesis. Two HICs but none of the PTCs harbored only proviruses with deletion in nef; these attenuated strains could contribute to viral control in these participants. We show, for the first time, the presence of intact proviruses and low viral diversity in PTCs long after treatment interruption, as well as the absence of evolution of the proviral quasispecies in subsequent samples. This reflects low residual replication over time. Further data are necessary to confirm these results. IMPORTANCE Most people living with HIV need antiretroviral therapy to control their infection and experience viral relapse in case of treatment interruption, because of viral reservoir (proviruses) persistence. Knowing that proviruses are very diverse and most of them are defective in treated individuals, we aimed to characterize the HIV blood reservoirs of posttreatment controllers (PTCs), rare models of drug-free remission, in comparison with spontaneous controllers and treated individuals. At a median time of 9 years after treatment interruption, which is unprecedented in the literature, we showed that the proportions and quantities of intact proviruses were similar between PTCs and other individuals. Unlike 2/7 spontaneous controllers who harbored only nef-deleted proviruses, which are attenuated strains, which could contribute to their control, no such case was observed in PTCs. Furthermore, PTCs displayed low viral genetic diversity and no evolution of their reservoirs, indicating very low residual replication, despite the presence of intact proviruses.

De novo birth of functional microproteins in the human lineage

Small open reading frames (sORFs) can encode functional "microproteins" that perform crucial biological tasks. However, their size makes them less amenable to genomic analysis, and their origins and conservation are poorly understood. Given their short length, it is plausible that some of these functional microproteins have recently originated entirely de novo from noncoding sequences. Here we sought to identify such cases in the human lineage by reconstructing the evolutionary origins of human microproteins previously found to have measurable, statistically significant fitness effects. By tracing the formation of each ORF and its transcriptional activation, we show that novel microproteins with significant phenotypic effects have emerged de novo throughout animal evolution, including two after the human-chimpanzee split. Notably, traditional methods for assessing coding potential would miss most of these cases. This evidence demonstrates that the functional potential intrinsic to sORFs can be relatively rapidly and frequently realized through de novo gene emergence.

Genome-level analyses resolve an ancient lineage of symbiotic ascomycetes

Ascomycota account for about two-thirds of named fungal species.1 Over 98% of known Ascomycota belong to the Pezizomycotina, including many economically important species as well as diverse pathogens, decomposers, and mutualistic symbionts.2 Our understanding of Pezizomycotina evolution has until now been based on sampling traditionally well-defined taxonomic classes.3,4,5 However, considerable diversity exists in undersampled and uncultured, putatively early-diverging lineages, and the effect of these on evolutionary models has seldom been tested. We obtained genomes from 30 putative early-diverging lineages not included in recent phylogenomic analyses and analyzed these together with 451 genomes covering all available ascomycete genera. We show that 22 of these lineages, collectively representing over 600 species, trace back to a single origin that diverged from the common ancestor of Eurotiomycetes and Lecanoromycetes over 300 million years BP. The new clade, which we recognize as a more broadly defined Lichinomycetes, includes lichen and insect symbionts, endophytes, and putative mycorrhizae and encompasses a range of morphologies so disparate that they have recently been placed in six different taxonomic classes. To test for shared hidden features within this group, we analyzed genome content and compared gene repertoires to related groups in Ascomycota. Regardless of their lifestyle, Lichinomycetes have smaller genomes than most filamentous Ascomycota, with reduced arsenals of carbohydrate-degrading enzymes and secondary metabolite gene clusters. Our expanded genome sample resolves the relationships of numerous "orphan" ascomycetes and establishes the independent evolutionary origins of multiple mutualistic lifestyles within a single, morphologically hyperdiverse clade of fungi.