Dataset Information

Mirage: estimation of ancestral gene-copy numbers by considering different evolutionary patterns among gene families.

ABSTRACT:

Motivation

Reconstruction of gene copy number evolution is an essential approach for understanding how complex biological systems have been organized. Although various models have been proposed for gene copy number evolution, existing evolutionary models have not appropriately addressed the fact that different gene families can have very different gene gain/loss rates.

Results

In this study, we developed Mirage (MIxtuRe model for Ancestral Genome Estimation), which allows different gene families to have flexible gene gain/loss rates. Mirage can use three models for formulating heterogeneous evolution among gene families: the discretized Γ model, probability distribution-free model and pattern mixture (PM) model. Simulation analysis showed that Mirage can accurately estimate heterogeneous gene gain/loss rates and reconstruct gene-content evolutionary history. Application to empirical datasets demonstrated that the PM model fits genome data from various taxonomic groups better than the other heterogeneous models. Using Mirage, we revealed that metabolic function-related gene families displayed frequent gene gains and losses in all taxa investigated.

Availability and implementation

The source code of Mirage is freely available at https://github.com/fukunagatsu/Mirage.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

SUBMITTER: Fukunaga T

PROVIDER: S-EPMC9710636 | biostudies-literature | 2021

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Mirage: estimation of ancestral gene-copy numbers by considering different evolutionary patterns among gene families.

Fukunaga Tsukasa T Iwasaki Wataru W

Bioinformatics advances 20210730 1

<h4>Motivation</h4>Reconstruction of gene copy number evolution is an essential approach for understanding how complex biological systems have been organized. Although various models have been proposed for gene copy number evolution, existing evolutionary models have not appropriately addressed the fact that different gene families can have very different gene gain/loss rates.<h4>Results</h4>In this study, we developed Mirage (MIxtuRe model for Ancestral Genome Estimation), which allows differen ...[more]

PMID: 36700099

Similar Datasets

Project description:Carbapenemase production is one of the leading mechanisms of carbapenem resistance in Gram-negative bacteria. An increase in carbapenemase gene (blaCarb) copies is an important mechanism of carbapenem resistance. No currently available bioinformatics tools allow for reliable detection and reporting of carbapenemase gene copy numbers. Here, we describe the carbapenemase-encoding gene copy number estimator (CCNE), a ready-to-use bioinformatics tool that was developed to estimate blaCarb copy numbers from whole-genome sequencing data. Its performance on Klebsiella pneumoniae carbapenemase gene (blaKPC) copy number estimation was evaluated by simulation and quantitative PCR (qPCR), and the results were compared with available algorithms. CCNE has two components, CCNE-acc and CCNE-fast. CCNE-acc detects blaCarb copy number in a comprehensive and high-accuracy way, while CCNE-fast rapidly screens blaCarb copy numbers. CCNE-acc achieved the best accuracy (100%) and the lowest root mean squared error (RMSE; 0.07) in simulated noise data sets, compared to the assembly-based method (23.4% accuracy, 1.697 RMSE) and the OrthologsBased method (78.9% accuracy, 0.395 RMSE). In the qPCR validation, a high consistency was observed between the blaKPC copy number determined by qPCR and that determined with CCNE. Reverse transcription-qPCR transcriptional analysis of 40 isolates showed that blaKPC expression was positively correlated with the blaKPC copy numbers detected by CCNE (P < 0.001). An association study of 357 KPC-producing K. pneumoniae isolates and their antimicrobial susceptibility identified a significant association between the estimated blaKPC copy number and MICs of imipenem (P < 0.001) and ceftazidime-avibactam (P < 0.001). Overall, CCNE is a useful genomic tool for the analysis of antimicrobial resistance genes copy number; it is available at https://github.com/biojiang/ccne. IMPORTANCE Globally disseminated carbapenem-resistant Enterobacterales is an urgent threat to public health. The most common carbapenem resistance mechanism is the production of carbapenemases. Carbapenemase-producing isolates often exhibit a wide range of carbapenem MICs. Higher carbapenem MICs have been associated with treatment failure. The increase of carbapenemase gene (blaCarb) copy numbers contributes to increased carbapenem MICs. However, blaCarb gene copy number detection is not routinely conducted during a genomic analysis, in part due to the lack of optimal bioinformatics tools. In this study, we describe a ready-to-use tool we developed and designated the carbapenemase-encoding gene copy number estimator (CCNE) that can be used to estimate the blaCarb copy number directly from whole-genome sequencing data, and we extended the data to support the analysis of all known blaCarb genes and some other antimicrobial resistance genes. Furthermore, CCNE can be used to interrogate the correlations between genotypes and susceptibility phenotypes and to improve our understanding of antimicrobial resistance mechanisms.

Project description:BackgroundFamily Rhinolophidae (horseshoe bats), Hipposideridae (leaf-nosed bats) and Rhinonycteridae (trident bats) are exclusively distributed in the Old-World, and their biogeography reflects the complex historic geological events throughout the Cenozoic. Here we investigated the origin of these families and unravel the conflicting family origin theories using a high resolution tree covering taxa from each zoogeographic realm from Africa to Australia. Ancestral range estimations were performed using a probabilistic approach implemented in BioGeoBEARS with subset analysis per biogeographic range [Old-World as whole, Australia-Oriental-Oceania (AOO) and Afrotropical-Madagascar-Palearctic (AMP)].ResultOur result supports an Oriental origin for Rhinolophidae, whereas Hipposideridae originated from the Oriental and African regions in concordance with fossil evidence of both families. The fossil evidence indicates that Hipposideridae has diversified across Eurasia and the Afro-Arabian region since the Middle Eocene. Meanwhile, Rhinonycteridae (the sister family of Hipposideridae) appears to have originated from the Africa region splitting from the common ancestor with Hipposideridae in Africa. Indomalaya is the center of origin of Rhinolophidae AOO lineages, and Indomalayan + Philippines appears to be center of origin of Hipposideridae AOO lineage indicating allopatric speciation and may have involved jump-dispersal (founder-event) speciation within AOO lineage. Wallacea and the Philippines may have been used as stepping stones for dispersal towards Oceania and Australia from the Oriental region. Multiple colonization events via different routes may have occurred in the Philippines (i.e., Palawan and Wallacea) since the Late Miocene. The colonization of Rhinolophidae towards Africa from Asia coincided with the estimated time of Tethys Ocean closure around the Oligocene to Miocene (around 27 Ma), allowing species to disperse via the Arabian Peninsula. Additionally, the number of potential cryptic species in Rhinolophidae in Southeast Asia may have increased since Plio-Pleistocene and late Miocene.ConclusionOverall, we conclude an Oriental origin for Rhinolophidae, and Oriental + African for Hipposideridae. The result demonstrates that complex historical events, in addition to species specific ecomorphology and specialization of ecological niches may shape current distributions.

Dataset Information

Mirage: estimation of ancestral gene-copy numbers by considering different evolutionary patterns among gene families.

Motivation

Results

Availability and implementation

Supplementary information

Publications

Mirage: estimation of ancestral gene-copy numbers by considering different evolutionary patterns among gene families.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets