Comparative Time–Space Evolution of Bundibugyo Ebolavirus Glycoprotein Across the 2007, 2012, and 2026 Outbreaks

Technical Summary: Comparative Time–Space Evolution of Bundibugyo Ebolavirus Glycoprotein

Problem Statement
Bundibugyo ebolavirus (BDBV) has caused three documented outbreaks: Uganda (2007), the Democratic Republic of the Congo (DRC) (2012), and a cross-border DRC/Uganda event (2026). While complete glycoprotein (GP) sequences are available for these events, previous research has largely focused on individual outbreaks, phylogenetic reconstruction, or structural characterization in isolation. There is a lack of comprehensive comparative analyses that integrate temporal, geographical, and functional evolutionary characteristics across all three outbreaks within a unified positional framework. Consequently, the temporal continuity of amino acid substitutions, the persistence or reversion of specific mutations, and their distribution across functional domains have not been systematically characterized.

Methodology
The study analyzed 54 complete Bundibugyo ebolavirus GP amino acid sequences (676 residues each) representing the 2007, 2012, and 2026 outbreaks. The core methodological approach involved:

• Unified Coordinate Framework: All sequences were aligned without trimming or coordinate renumbering, establishing a single amino acid coordinate system (residues 1–676) to enable direct residue-by-residue comparison across all isolates.
• Substitution Identification: The aligned sequences were systematically examined to identify reproducible amino acid differences, distinguishing core outbreak-associated substitutions from sporadic variation.
• Evolutionary Classification: Identified substitutions were categorized into five distinct patterns based on their temporal distribution:
  1. Fixation: Substitutions established in later outbreaks.
  2. Reversion: Substitutions returning to amino acid states observed in earlier outbreaks.
  3. Emergence: Substitutions first detected in the 2026 outbreak.
  4. Within-outbreak polymorphism: Heterogeneous residues observed among isolates from the same outbreak period.
  5. Rare variant: Substitutions detected in only a single isolate.
• Functional Mapping: Substitution sites were mapped onto established functional domains of the GP, including the receptor-binding domain, mucin-like domain (MLD), furin cleavage site, internal fusion loop, immunosuppressive motif, and transmembrane domain.

Key Results

• Structural Conservation: The GP architecture remained highly conserved across all three outbreaks. All isolates encoded a protein of identical length (676 amino acids), with no outbreak-specific insertions or deletions affecting overall protein length. Evolutionary changes occurred exclusively through amino acid substitutions.
• Identification of Core Sites: The analysis identified 21 core amino acid substitution sites distributed throughout the GP sequence.
• Evolutionary Patterns: The 21 sites were classified as follows:
  • Fixation: 10 sites
  • Reversion: 4 sites
  • Emergence: 2 sites
  • Within-outbreak polymorphism: 4 sites
  • Rare variant: 1 site
• Domain-Specific Distribution: Substitutions were not uniformly distributed. They were concentrated primarily within the mucin-like domain (MLD) and the GP1 C-terminal region. Conversely, functionally critical regions—the furin cleavage site, internal fusion loop, immunosuppressive motif, and transmembrane domain—remained completely conserved across all three outbreak periods.

Key Contributions

• Standardized Comparative Framework: The study establishes a unified amino acid coordinate system that allows for direct comparison of homologous residues across successive outbreaks without the need for coordinate conversion. This facilitates the continuous tracking of individual substitution sites over time.
• Evolutionary Classification System: By categorizing substitutions into fixation, reversion, emergence, polymorphism, and rare variants, the study provides a temporal perspective on viral evolution that distinguishes established changes from transient or newly emerging variations.
• Integrated Time-Space Overview: The research offers the first integrated view of BDBV GP evolution across the three documented outbreaks, linking specific genetic changes to their geographical and temporal origins within a functional context.

Significance and Claims
The authors claim that this study provides a reproducible reference dataset and analytical framework for future genomic surveillance of Bundibugyo ebolavirus. By demonstrating that the virus evolves primarily through substitutions rather than insertions/deletions, and that functional constraints preserve essential domains while allowing variability in flexible regions (like the MLD), the study aids in interpreting antigenic variation and potential therapeutic targets. The unified framework is presented as a tool to facilitate the integration of temporal, geographical, and functional information, supporting comparative investigations of ebolavirus evolution and the monitoring of newly emerging lineages. The authors note that their findings are based on publicly available sequences and represent a structural and comparative analysis rather than experimental validation of biological consequences.