Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “phylogenomics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Phylogenomic inference of protein molecular function: advances and challenges.

MOTIVATION: Protein families evolve a multiplicity of functions through gene duplication, speciation and other processes. As a number of studies have shown, standard methods of protein function prediction produce systematic errors on these data. Phylogenomic analysis--combining phylogenetic tree construction, integration of experimental data and differentiation of orthologs and paralogs--has been proposed to address these errors and improve the accuracy of functional classification. The explicit integration of structure prediction and analysis in this framework, which we call structural phylogenomics, provides additional insights into protein superfamily evolution. RESULTS: Results of protein functional classification using phylogenomic analysis show fewer expected false positives overall than when pairwise methods of functional classification are employed. We present an overview of the motivations and fundamental principles of phylogenomic analysis, new methods developed for the key tasks, benchmark datasets for these tasks (when available) and suggest procedures to increase accuracy. We also discuss some of the methods used in the Celera Genomics high-throughput phylogenomic classification of the human genome. AVAILABILITY: Software tools from the Berkeley Phylogenomics Group are available at http://phylogenomics.berkeley.edu

Algorithms↗

Phydbac (phylogenomic display of bacterial genes): An interactive resource for the annotation of bacterial genomes.

Phydbac is a web interactive resource based on phylogenomic profiling, designed to help microbiologists to annotate bacterial proteins. Phylogenomic annotation is based on the assumption that functionally linked protein-coding genes must evolve in a coordinated manner. The detection of subsets of co-evolving genes within a given genome involves the computation of protein sequence conservation profiles across a spectrum of microbial species, followed by the identification of significant pairwise correlations between them. Many ongoing studies are devoted to the problem of computing the most biologically significant phylogenomic profiles and how best identifying clusters of 'functionally interacting' genes. Here we introduce a web tool, Phydbac, allowing the dynamic construction of phylogenomic profiles of protein sequences of interest and their interactive display. In addition, Phydbac can identify Escherichia coli proteins exhibiting the evolution pattern most similar to arbitrary query protein sequences, hence providing functional hints for open reading frames (ORFs) of hypothetical or unknown function. The phylogenomic profiles of all E.coli K-12 protein-coding genes are pre-computed, allowing queries about E.coli genes to be answered instantaneously. The profiles and phylogenomic neighborhoods are computed using an original method shown to perform better than previous ones. An extension of Phydbac, including precomputed profiles for all available bacterial genomes (including major pathogens) will soon be available. Phydbac can be accessed at: http://igs-server.cnrs-mrs.fr/phydbac/.

Bacterial Proteins↗

Phylogenomic subsampling and upsampling for efficient evolutionary analyses of big data.

Long runtimes, high memory demands, and reliance on high-performance computing impede phylogenomic analyses. We review a scalable phylogenomic subsampling with upsampling (PSU) framework to address this challenge, which reduces runtime and memory requirements by orders of magnitude. In PSU, small subsamples of sites from a concatenated alignment are analyzed, which are expanded by upsampling before inference, and the resulting inferences are aggregated to obtain evolutionary estimates. PSU harnesses the fact that the computational cost of maximum likelihood analysis is strongly influenced by the number of distinct site patterns in the concatenated alignment, whereas statistical power depends primarily on the amount of evolutionary information represented by the total number of sites and substitutions. By reducing the former while restoring the latter through upsampling, PSU can approximate many full-alignment analyses at substantially lower computational cost. Analysis of simulated and empirical datasets shows that PSU can accurately estimate bootstrap support values, select the optimal substitution model, test evolutionary hypotheses, and infer branch lengths, divergence times, and associated uncertainty measures. PSU also provides distributions of inferred clade support across independent subsamples, enabling detection of conflicting phylogenetic signals that may remain hidden in conventional bootstrap analysis of concatenated alignments. Automated tuning of subsample size, the number of subsamples, and the number of upsampling replicates make PSU practical. We suggest that PSU is a general approach for scalable phylogenomic inference using a broad range of statistical methods. By enabling analyses of genome-scale alignments on commodity hardware, PSU broadens research access and reduces environmental and infrastructural costs of big-data phylogenomics.

Phylogeny↗

Phylogenomic subsampling and upsampling for efficient evolutionary analyses of big data.

Long runtimes, high memory demands, and reliance on high-performance computing impede phylogenomic analyses. We review a scalable phylogenomic subsampling with upsampling (PSU) framework, in which small subsamples of sites from a concatenated alignment are expanded by upsampling before inference, and the resulting analyses are then aggregated to obtain evolutionary estimates. PSU harnesses the fact that the computational cost of maximum likelihood analysis is strongly influenced by the number of distinct site patterns in the concatenated alignment, whereas statistical power depends primarily on the amount of evolutionary information represented by the total number of sites and substitutions. By reducing the former while restoring the latter through upsampling, PSU can approximate many full-data analyses at substantially lower computational cost. Analysis of simulated and empirical datasets shows that PSU can accurately estimate bootstrap support values, select the optimal substitution model, test evolutionary hypotheses, and infer branch lengths, divergence times, and associated uncertainty measures, while reducing runtime and memory requirements by orders of magnitude. PSU also provides distributions of inferred clade support across independent subsamples, enabling detection of conflicting phylogenetic signals that may remain hidden in conventional bootstrap analysis. Automated tuning of subsample size, the number of subsamples, and the number of upsampling replicates make PSU practical across diverse datasets. We suggest that PSU is a general strategy for scalable phylogenomic inference using a broad range of statistical methods. By enabling analyses of genome-scale alignments on commodity hardware, PSU broadens research access and reduces environmental and infrastructural costs of big-data phylogenomics.

confidence limits↗

Moraxella species isolated from blood cultures in Europe (MORAXEu): a multicentre study of epidemiology and antimicrobial susceptibility with complementary phylogenomic analysis of publicly available genomes.

INTRODUCTION: Moraxella species are fastidious Gram-negative bacteria capable of causing opportunistic infections, including bloodstream infections, especially in immunocompromised patients. Data on their epidemiology, antimicrobial susceptibility, and phylogenomics in Europe remains limited. METHODS: We conducted a multicentre, retrospective, observational study across 56 European hospital centres between January 1st 2020 and December 31st 2024. All Moraxella species isolated from blood cultures (BCs) were included. Species distribution and antimicrobial susceptibility profiles were analysed. We also performed a phylogenomic analysis of Moraxella genomes deposited in GenBank. RESULTS: A total of 709 Moraxella isolates were included. Moraxella osloensis (61.1%; n = 433/709) and Moraxella catarrhalis (20.4%; n = 145/709) were the most frequently identified species, followed by Moraxella nonliquefaciens (6.5%; n = 46/709) and Moraxella atlantae (4.8%; n = 34/709). Species distribution differed by age. M. catarrhalis was predominant in paediatric patients, whereas M. atlantae was more common in adults. Most isolates showed > 90% susceptibility to amoxicillin/clavulanate, cefotaxime, fluoroquinolones, and trimethoprim/sulfamethoxazole. Cefotaxime resistance in M. osloensis was more frequent in adults than in children (49% vs. 8%, p = 0.009). Phylogenomic analysis demonstrated the distinction of a core Moraxella group from the divergent Faucicola lineage, with M. osloensis and M. atlantae clustering within the latter. CONCLUSIONS: The epidemiology of Moraxella species from BCs in Europe showed age-group-specific differences in species distribution and generally favourable antimicrobial susceptibility patterns. Phylogenomic data corroborated recent taxonomic revisions, highlighting the need for improved diagnostics, harmonized nomenclature and sustained surveillance to inform management and stewardship of Moraxella bacteraemia.

Faucicola↗

Phydbac2: improved inference of gene function using interactive phylogenomic profiling and chromosomal location analysis.

Phydbac (phylogenomic display of bacterial genes) implemented a method of phylogenomic profiling using a distance measure based on normalized BLAST scores. This method was able to increase the predictive power of phylogenomic profiling by about 25% when compared to the classical approach based on Hamming distances. Here we present a major extension of Phydbac (named here Phydbac2), that extends both the concept and the functionality of the original web-service. While phylogenomic profiles remain the central focus of Phydbac2, it now integrates chromosomal proximity and gene fusion analyses as two additional non-similarity-based indicators for inferring pairwise gene functional relationships. Moreover, all presently available (January 2004) fully sequenced bacterial genomes and those of three lower eukaryotes are now included in the profiling process, thus increasing the initial number of reference genomes (71 in Phydbac) to 150 in Phydbac2. Using the KEGG metabolic pathway database as a benchmark, we show that the predictive power of Phydbac2 is improved by 27% over the previous version. This gain is accounted for on one hand, by the increased number of reference genomes (11%) and on the other hand, as a result of including chromosomal proximity into the distance measure (16%). The expanded functionality of Phydbac2 now allows the user to query more than 50 different genomes, including at least one member of each major bacterial group, most major pathogens and potential bio-terrorism agents. The search for co-evolving genes based on consensus profiles from multiple organisms, the display of Phydbac2 profiles side by side with COG information, the inclusion of KEGG metabolic pathway maps the production of chromosomal proximity maps, and the possibility of collecting and processing results from different Phydbac queries in a common shopping cart are the main new features of Phydbac2. The Phydbac2 web server is available at http://igs-server.cnrs-mrs.fr/phydbac/.

Artificial Gene Fusion↗

Novel relationships among ten fish model species revealed based on a phylogenomic analysis using ESTs.

The power of comparative phylogenomic analyses also depends on the amount of data that are included in such studies. We used expressed sequence tags (ESTs) from fish model species as a proof of principle approach in order to test the reliability of using ESTs for phylogenetic inference. As expected, the robustness increases with the amount of sequences. Although some progress has been made in the elucidation of the phylogeny of teleosts, relationships among the main lineages of the derived fish (Euteleostei) remain poorly defined and are still debated. We performed a phylogenomic analysis of a set of 42 of orthologous genes from 10 available fish model systems from seven different orders (Salmoniformes, Siluriformes, Cypriniformes, Tetraodontiformes, Cyprinodontiformes, Beloniformes, and Perciformes) of euteleostean fish to estimate divergence times and evolutionary relationships among those lineages. All 10 fish species serve as models for developmental, aquaculture, genomic, and comparative genetic studies. The phylogenetic signal and the strength of the contribution of each of the 42 orthologous genes were estimated with randomly chosen data subsets. Our study revealed a molecular phylogeny of higher-level relationships of derived teleosts, which indicates that the use of multiple genes produces robust phylogenies, a finding that is expected to apply to other phylogenetic issues among distantly related taxa. Our phylogenomic analyses confirm that the euteleostean superorders Ostariophysi and Acanthopterygii are monophyletic and the Protacanthopterygii and Ostariophysi are sister clades. In addition, and contrary to the traditional phylogenetic hypothesis, our analyses determine that killifish (Cyprinodontiformes), medaka (Beloniformes), and cichlids (Perciformes) appear to be more closely related to each other than either of them is to pufferfish (Tetraodontiformes). All 10 lineages split before or during the fragmentation of the supercontinent Pangea in the Jurassic.

Animals↗

Phylomat: an automated protein motif analysis tool for phylogenomics.

Recent progress in genomics, proteomics, and bioinformatics enables unprecedented opportunities to examine the evolutionary history of molecular, cellular, and developmental pathways through phylogenomics. Accordingly, we have developed a motif analysis tool for phylogenomics (Phylomat, http://alg.ncsa.uiuc.edu/pmat) that scans predicted proteome sets for proteins containing highly conserved amino acid motifs or domains for in silico analysis of the evolutionary history of these motifs/domains. Phylomat enables the user to download results as full protein or extracted motif/domain sequences from each protein. Tables containing the percent distribution of a motif/domain in organisms normalized to proteome size are displayed. Phylomat can also align the set of full protein or extracted motif/domain sequences and predict a neighbor-joining tree from relative sequence similarity. Together, Phylomat serves as a user-friendly data-mining tool for the phylogenomic analysis of conserved sequence motifs/domains in annotated proteomes from the three domains of life.

Algorithms↗

Functional genome annotation through phylogenomic mapping.

Accurate determination of functional interactions among proteins at the genome level remains a challenge for genomic research. Here we introduce a genome-scale approach to functional protein annotation--phylogenomic mapping--that requires only sequence data, can be applied equally well to both finished and unfinished genomes, and can be extended beyond single genomes to annotate multiple genomes simultaneously. We have developed and applied it to more than 200 sequenced bacterial genomes. Proteins with similar evolutionary histories were grouped together, placed on a three dimensional map and visualized as a topographical landscape. The resulting phylogenomic maps display thousands of proteins clustered in mountains on the basis of coinheritance, a strong indicator of shared function. In addition to systematic computational validation, we have experimentally confirmed the ability of phylogenomic maps to predict both mutant phenotype and gene function in the delta proteobacterium Myxococcus xanthus.

Bacterial Proteins↗

Annotation of bacterial genomes using improved phylogenomic profiles.

MOTIVATION: Phylogenomic profiling is a large-scale comparative genomic method used to infer protein function from evolutionary information first described in a binary form by Pellegrini et al. (1999). Here, we propose improvements of this approach including the use of normalized Blastp bit scores, a normalization of the matrix of profiles to take into account the evolutionary distances between bacteria, the definition of a phylogenomic neighborhood based on continuous pairwise distances between genes and an original annotation procedure including the computation of a p-value for each functional assignment. RESULTS: The method presented here increases the number of Ecocyc enzymes identified as being evolutionarily related by about 25% with respect to the original binary form (absent/present) method. The fraction of 'false' positives is shown to be smaller than 20%. Based on their phylogenomic relationships, genes of unknown function can then be automatically related to annotated genes. Each gene annotation predicted is associated with a p-value, i.e. its probability to be obtained by chance. The validity of this method was extensively tested on a large set of genes of known function using the MultiFun database. We find that 50% of 3122 function attributions that can be made at a p-value level of 10(-11) correspond to the actual gene annotation. The method can be readily applied to any newly sequenced microbial genome. In contrast to earlier work on the same topic, our approach avoids the use of arbitrary cut-off values, and provides a reliability estimate of the functional predictions in form of p-values.

Algorithms↗

Comparative phylogenomics and transcriptional regulatory networks of AQPs, HSPs, and LEA proteins in salt-stressed Portulaca oleracea.

Soil salinization severely threatens global food security, necessitating systematic investigations of halophytes like Portulaca oleracea to decode the molecular mechanisms of environmental resilience. Utilizing an integrated framework of deep learning-based genome annotation (58,817 predicted genes; 96.5% BUSCO completeness), multi-tissue RNA-Seq, phylogenomics, and gene regulatory network (GRN) inference, the synergistic orchestration of 78 aquaporins (AQPs), 525 heat shock proteins (HSPs), and 119 late embryogenesis abundant (LEA) proteins was elucidated. The active transcriptome, encompassing 39,065 expressed loci, revealed a systemic growth-defense trade-off. Tissues displayed distinct adaptive mechanisms: leaves modulated intracellular water balance via specialized AQPs, whereas adult roots maintained proteostasis through robust HSP20/HSP70 induction. Phylogenomic clustering across 154 species demonstrated that salinity tolerance constitutes an evolutionary mosaic, identifying 81 halophyte-exclusive orthogroups and 1129 species-specific clusters. Comparative topology across six independent GRNs (4.2M-5.3 M edges) unmasked a highly modular transcriptional reprogramming strategy governed by a core apparatus of 22 stress-exclusive regulators, with functional enrichment heavily prioritizing protein dimerization and chromatin remodeling. Theoretically, the distinct convergence of Trihelix transcription factors with guard cell differentiation pathways offers a candidate transcriptomic framework to explain the plant's characteristic C4-CAM photosynthetic plasticity under severe osmotic pressure. Practically, these evolutionary blueprints and specific master switches transcend single-gene transgenic limitations. Utilizing these root-sustained and stress-inducible targets under localized promoters provides a naturally optimized, network-level precision engineering roadmap to transfer robust, compartmentalized halotolerance to sensitive glycophytic crops.

Gene Regulatory Networks↗

Phylogenomics of Desulfuromonadia supports reclassification of Geobacter psychrophilus as Irobacter psychrophilus comb. nov. and proposal of Geosyntrophus gen. nov.

Genome-resolved phylogenomics reveals widespread misclassification of metal-reducing bacteria historically assigned to Geobacter based on 16S rRNA gene phylogeny, and highlights species that persist only as 16S rRNA entries without genomes for robust taxonomic resolution. Here, we resolve two such lineages by integrating whole-genome phylogeny with average amino acid identity (AAI) and percentage of conserved proteins (POCP) across 418 dereplicated genomes of Desulfuromonadia. We report a draft genome of the psychrophilic iron-reducing bacterium Geobacter psychrophilus (100% completeness). Phylogenomic analyses place both Geobacter psychrophilus and the GTDB placeholder genus g__JACRCG01 within the family 'Pseudopelobacteraceae', outside Geobacteraceae sensu stricto. Within this framework, G. psychrophilus forms a distinct, well-supported lineage separated from neighbouring genera by discontinuities in AAI and POCP, supporting its reclassification as Irobacter psychrophilus comb. nov. Additionally, we show that Geosyntrophus acetoxidans, a non-axenic syntrophic bacterium, forms a coherent genus with 51 other environmental genomes (placeholder genus g__JACRCG01), for which we propose the replacement name Geosyntrophus gen. nov. Comparative genome analysis revealed conserved family-level metabolic traits together with genus-specific differences in respiratory metabolism, while ANI-based clustering identified substantial species-level diversity within both proposed genera. Metagenome and 16S rRNA-gene survey data further show that Geosyntrophus and Irobacter occur in broadly similar aquatic and subsurface habitats spanning from the Arctic to the Antarctic. Together, these results resolve the taxonomy of two previously ambiguous Desulfuromonadales lineages and shed light on their environmental distribution.

AAI↗

Phylogenomics and female reproductive morphology reframe the classification of the Halymeniales (Rhodophyta).

The red algal order Halymeniales (Rhodophyta) exhibits remarkable morphological and taxonomic diversity but its higher-level relationships remain poorly resolved. Here, we present a comprehensive phylogenomic analysis based on newly generated plastid (170 protein-coding genes), mitochondrial (23 genes), and complete nuclear ribosomal cistron sequences from 56 taxa, complemented with an expanded rbcL dataset encompassing 334 sequences. Our results provide a robust phylogenomic framework for the Halymeniales, offering a taxonomic backbone for future systematic studies. The analyses consistently recover six early-diverging lineages (Acrodiscus, Isabbottia, Norrissia, Pachymenia, Zymurgia, and Tsengia) and two strongly supported larger clades (Halymenia s.l. and Grateloupia s.l.). While most small and recently described genera are monophyletic, several traditional genera (e.g., Halymenia, Cryptonemia, Grateloupia) are poly- or paraphyletic, requiring considerable taxonomic revision. At the family level, the data indicate that reinstatement of the Grateloupiaceae sensu Kim et al. (2021) would entail a revised circumscription of the Halymeniaceae and the recognition of at least five small families to accommodate the early-diverging lineages. Although such a revised classification would result in monophyletic families, it is not supported by morpho-anatomical characters. Instead, we propose a more stable two-family system, recognizing a broadly circumscribed Halymeniaceae that is sister to the Tsengiaceae. Female reproductive characters, particularly the structure of carpogonial and auxiliary cell ampullae, support this two-family system and further characterize many genus-level clades, although substantial convergence across lineages exists.

Phylogeny↗

Phylogenomics and evolution of the Lauraceae based on targeted capture data.

The family Lauraceae, a hyper-diverse magnoliid family comprising approximately 63 genera and over 3,000 species, plays a key ecological role in tropical and subtropical forests. Yet deep relationships among its nine tribes remain unresolved, likely due to limited sampling and complex evolutionary processes such as incomplete lineage sorting (ILS) and gene flow. To address these challenges, we generated datasets of 255 single-copy nuclear genes and chloroplast genomes using a newly designed Lauraceae-specific probe set, achieving the most comprehensive genus-level sampling (84%) to date. Phylogenomic analyses reconstructed a robust nuclear tree, which resolved the Neocinnamomeae as sister to the Caryodaphnopsideae and revealed pronounced gene tree conflict and pervasive cytonuclear discordance. To investigate the evolutionary processes underlying these patterns, comprehensive analyses were conducted. The results indicate that conflicting nuclear gene trees reflect the combined effects of ILS, gene tree estimation error, and gene flow, with ILS dominating across the core Lauraceae, whereas cytonuclear discordance is primarily driven by extensive gene flow. Diversification analyses further indicate that episodes of rapid lineage accumulation coincide with major gene flow events, suggesting a potential role of gene flow in the diversification of Lauraceae. Overall, this study provides a robust nuclear phylogenomic framework for Lauraceae and demonstrates that gene flow had profound effects on its evolutionary history, shedding light on the contribution of gene flow to the diversification of hyper-diverse tropical plant lineages.

Cytonuclear discordance↗

Phylogenomic analysis of the alpha proteasome gene family from early-diverging eukaryotes.

We employed a phylogenomic approach to study the evolution of alpha subunits of the proteasome gene family from early diverging eukaryotes. BLAST similarity searches of the Giardia lamblia genome identified all seven alpha proteasome genes characteristic of eukaryotes from the crown group. In addition, a PCR strategy for the amplification of multiple alpha subunit sequences generated single alpha proteasome products for representatives of the Kinetoplastida (Leishmania major), the Parabasalia (Trichomonas vaginalis), and the Microsporidia (Vairimorpha sp., Nosema sp., Endoreticulata sp., and Spraguea lophii). The kinetoplastid Trypanosoma cruzi and the eukaryote crown group Acanthamoeba castellanii yielded two distinct alpha proteasome genes each. The presence of seven distinct alpha proteasome genes in G. lamblia, one of the earliest-diverging eukaryotes, indicates that the alpha proteasome gene family evolved rapidly from a minimum of one gene in Archaea to seven or more in Eukarya. Results from the phylogenomic analysis are consistent with the idea that the Diplomonida (as represented by G. lamblia), the Kinetoplastida, the Parabasalia, and the Microsporidia diverged after the duplication events that originated the alpha proteasome gene family. A model for the early origin and evolution of the proteasome gene family is presented.

Base Sequence↗

Identification of the mammalian Not gene via a phylogenomic approach.

Despite the great morphological diversity of early embryos, the underlying mechanisms of gastrulation are known to be broadly conserved in vertebrates. However, a number of genes characterized as fulfilling an essential function in this process in several model organisms display no clear ortholog in mammalian genomes. We have devised an in silico phylogenomic approach, based on exhaustive similarity searches in vertebrate genomes and subsequent bayesian phylogenetic analyses, to identify such missing genes, presumed to be highly divergent. This approach has been used to identify mammalian orthologs of Not, an homeodomain containing gene previously characterized in Xenopus, chick and zebrafish as playing a critical role in the formation of the notochord. This attempt led to the identification of a highly divergent mammalian Not-related gene in the mouse, human and rat. The results from phylogenetic reconstructions, synteny analyses, expression pattern analyses in wild-type and mutant mouse embryos, and overexpression experiments in Xenopus embryos converge to confirm these genes as representatives of the Not family in mammals. The identification of the mammalian Not gene delivers an important component for the understanding of the genetics underlying notochord formation in mammals and its evolution among vertebrates. The phylogenomic method used to retrieve this gene thus provides a tool, which can complement or validate genome annotations in situations when they are weakly supported.

Amino Acid Sequence↗

Beyond species trees: pervasive gene flow limits phylogenomic resolution in the diversification of Juniperus from the Qinghai-Tibet Plateau.

Understanding how lineages diversify despite persistent ancestral polymorphism and recurrent gene flow remains a central challenge in evolutionary biology. Juniperus distributed across the Qinghai-Tibet Plateau provide an ideal system for addressing this question because repeated geological uplift and climatic oscillations have likely promoted cycles of lineage divergence, range shifts, and secondary contact. Here, we combined approximately 1.08 million genome-wide SNPs from 164 individuals representing thirteen Juniperus lineages with phylogenomic datasets comprising 3,381 nuclear single-copy genes and nearly complete plastomes. We detected extensive phylogenomic discordance and cytonuclear incongruence across genomic datasets. Topology weighting, coalescent simulations, quartet-based tests, and analyses of gene flow and reticulation collectively support the interpretation that these patterns were shaped by the combined effects of prolonged incomplete lineage sorting and gene flow during lineage diversification. Ecological niche analyses further provide a spatial and climatic context in which environmentally similar lineages may have had greater opportunities for secondary contact during historical range shifts. Collectively, our results reveal that the evolutionary history of Qinghai-Tibet Plateau Juniperus is characterized by reticulate diversification rather than strictly bifurcating evolution, and demonstrate how genome-wide discordance can provide biological insights into the evolutionary processes underlying lineage diversification.

Gene Flow↗

Nuclear single-copy orthologous genes as phylogenomic markers for resolving the closely related firefly genera Pteroptyx, Medeopteryx, and Trisinuata (Coleoptera: Lampyridae: Luciolinae).

Fireflies (Lampyridae) are bioluminescent beetles with broad ecological roles across temperate and tropical ecosystems, occupying diverse habitats including forests, wetlands, grasslands, mangroves, and riverine systems. The subfamily Luciolinae is primarily distributed across Asia and the Indo-Pacific. Phylogenetic relationships among three closely related Luciolinae genera - Medeopteryx, Pteroptyx, and Trisinuata - remain unresolved using mitochondrial genome data alone. This study used nuclear genome data to resolve relationships among these genera and identify a lighter-weight nuclear marker panel for expanding taxon sampling. Draft genomes were reconstructed for fifteen firefly species, eight from the focal genera, and analyzed with five published firefly genomes. Using BUSCO and OrthoFinder, 1,011 nuclear single-copy orthologs (SCOs) were identified for phylogenomic inference. Discordance between concatenation- and coalescence-based phylogenies indicated incomplete lineage sorting (ILS). The coalescence-based phylogeny recoveredPteroptyxas monophyletic and sister to a (Medeopteryx,Trisinuata) clade, with Trisinuata nested within a non-monophyletic Medeopteryx; however, quartet support at the base of Pteroptyx, particularly at Pt. valida, was low.Filtering for compositional homogeneity, clock-likeness, and species-tree concordance yielded 103 SCOs with a significantly higher proportion of parsimony-informative sites than non-selected loci, retaining the backbone topology with higher gene concordance support at scored clades, while ILS-driven discordance at Pt. valida persists - confirming that the reduced panel retains phylogenetic resolving power for future taxon sampling. These findings demonstrate a practical framework for using nuclear SCOs to resolve close phylogenetic relationships within Luciolinae. Future work should expand taxon sampling - especially forTrisinuata - alongside long-read assemblies, for a more robust phylogenomic framework.

Fireflies↗