Phylogenetic Clustering of Genes Reveals Shared Evolutionary Trajectories and Putative Gene Functions.
Ontology highlight
ABSTRACT:
SUBMITTER: Liu C
PROVIDER: S-EPMC6130602 | biostudies-literature | 2018 Sep
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Phylogenetic inference can potentially result in a more accurate tree using data from multiple loci. However, if the loci are incongruent-due to events such as incomplete lineage sorting or horizontal gene transfer-it can be misleading to infer a single tree. To address this, many previous contributions have taken a mechanistic approach, by modeling specific processes. Alternatively, one can cluster loci without assuming how these incongruencies might arise. Such "process-agnostic" approaches typically infer a tree for each locus and cluster these. There are, however, many possible combinations of tree distance and clustering methods; their comparative performance in the context of tree incongruence is largely unknown. Furthermore, because standard model selection criteria such as AIC cannot be applied to problems with a variable number of topologies, the issue of inferring the optimal number of clusters is poorly understood. Here, we perform a large-scale simulation study of phylogenetic distances and clustering methods to infer loci of common evolutionary history. We observe that the best-performing combinations are distances accounting for branch lengths followed by spectral clustering or Ward's method. We also introduce two statistical tests to infer the optimal number of clusters and show that they strongly outperform the silhouette criterion, a general-purpose heuristic. We illustrate the usefulness of the approach by 1) identifying errors in a previous phylogenetic analysis of yeast species and 2) identifying topological incongruence among newly sequenced loci of the globeflower fly genus Chiastocheta We release treeCl, a new program to cluster genes of common evolutionary history (http://git.io/treeCl).
| S-EPMC4868114 | biostudies-literature
Project description:Escherichia coli (E. coli) ST1193 is an emerging fluoroquinolones-resistant and virulent lineage. Large gaps remain in our understanding of the evolutionary processes and differences of this lineage. Therefore, we used 76 E. coli ST1193 genomes to detect strain-level genetic diversity and phylogeny of this lineage globally. All E. coli ST1193 possessed fimH64, filCH5, and fumC14. There was 94.7% of isolates classified as O-type O75. There was 9.33% of E. coli ST1193 that possessed K5 capsular, while 90.67% of isolates possessed K1 capsular. The core genome analysis revealed that all isolates were divided into two phylogenetic clades (clade A and B). Clade A included 25 non-Chinese E. coli ST1193, and clade B contained all isolates collected from Fuzhou, China, respectively. The results of comparative genomics indicated Indels were identified in 150 clade-specific genes, which were enriched into the biological process and molecular function. Accessory genome phylogenetic tree showed a high degree of correlation between accessory genome clusters and core genome clades. There was significant difference in antibiotic resistance genes (ARGs) [bla CTX-M-55 , bla TEM-1 , sul2, tet(B), tet(R), APH(6)-Id, and AAC(3)-IId], virulence factors (cia, neuC, gad, and traT), and plasmid replicon types (IncQ1, Col156, and IncB/O/K/Z) between clade A (non-Chinese isolates) and clade B (Chinese isolates) (p < 0.05). Further analysis of the genetic environments of bla CTX-M-55 demonstrated that the flanking contexts of bla CTX-M-55 were diverse. In conclusion, our results reveal the distinct evolutionary trajectories of the spread of E. coli ST1193 in Fuzhou, China and non-China regions. This supports both global transmission and localized lineage expansion of this lineage following specific introductions into a geographic locality.
| S-EPMC8602892 | biostudies-literature
Project description:Major depressive disorder (MDD) affects around 15% of the population at some stage in their lifetime. It can be gravely disabling and it is associated with increased risk of suicide. Genetics play an important role; however, there are additional environmental contributions to the pathogenesis. A number of possible risk genes that increase liability for developing symptoms of MDD have been identified in genome-wide association studies (GWAS). The goal of this study was to characterize the MDD risk genes with respect to the degree of evolutionary conservation in simpler model organisms such as Caenorhabditis elegans and zebrafish, the phenotypes associated with variation in these genes and the extent of network connectivity. The MDD risk genes showed higher conservation in C. elegans and zebrafish than genome-to-genome comparisons. In addition, there were recurring themes among the phenotypes associated with variation of these risk genes in C. elegans. The phenotype analysis revealed enrichment for essential genes with pleiotropic effects. Moreover, the MDD risk genes participated in more interactions with each other than did randomly-selected genes from similar-sized gene sets. Syntenic blocks of risk genes with common functional activities were also identified. By characterizing evolutionarily-conserved counterparts to the MDD risk genes, we have gained new insights into pathogenetic processes relevant to the emergence of depressive symptoms in man.
| S-EPMC8319724 | biostudies-literature
Project description:Resistance to the broad-spectrum antibiotic ciprofloxacin is detected at high rates for a wide range of bacterial pathogens. To investigate the dynamics of ciprofloxacin resistance development, we applied a comparative resistomics workflow for three clinically relevant species of Gram-negative bacteria: Escherichia coli, Acinetobacter baumannii, and Pseudomonas aeruginosa. We combined experimental evolution in a morbidostat with deep sequencing of evolving bacterial populations in time series to reveal both shared and unique aspects of evolutionary trajectories. Representative clone characterization by sequencing and MIC measurements enabled direct assessment of the impact of mutations on the extent of acquired drug resistance. In all three species, we observed a two-stage evolution: (i) early ciprofloxacin resistance reaching 4- to 16-fold the MIC for the wild type, commonly as a result of single mutations in DNA gyrase target genes (gyrA or gyrB), and (ii) additional genetic alterations affecting the transcriptional control of the drug efflux machinery or secondary target genes (DNA topoisomerase parC or parE). IMPORTANCE The challenge of spreading antibiotic resistance calls for systematic efforts to develop more "irresistible" drugs based on a deeper understanding of dynamics and mechanisms of antibiotic resistance acquisition. To address this challenge, we have established a comparative resistomics approach which combines experimental evolution in a continuous-culturing device, the morbidostat, with ultradeep sequencing of evolving microbial populations to identify evolutionary trajectories (mutations and genome rearrangements) leading to antibiotic resistance over a range of target pathogens. Here, we report the comparative resistomics study of three Gram-negative bacteria (Escherichia coli, Acinetobacter baumannii, and Pseudomonas aeruginosa), which revealed shared and species-specific aspects of the evolutionary landscape leading to robust resistance against the clinically important antibiotic ciprofloxacin. Despite some differences between morbidostat-deduced mutation profiles and those observed in clinical isolates of individual species, a cross-species comparative resistomics approach allowed us to recapitulate all types of clinically relevant ciprofloxacin resistance mechanisms. This observation supports the anticipated utility of this approach in guiding rational optimization of treatment regimens for current antibiotics and the development of novel antibiotics with minimized resistance propensities.
| S-EPMC8262867 | biostudies-literature
Project description:Dosage compensation (DC) on the X Chromosome counteracts the deleterious effects of gene loss on the Y Chromosome. However, DC is not efficient if the X Chromosome also degenerates. This indeed occurs in Drosophila miranda, in which both the neo-Y and the neo-X are under accelerated pseudogenization. To examine the generality of this pattern, we investigated the evolution of two additional neo-sex chromosomes that emerged independently in D. albomicans and D. americana and reanalyzed neo-sex chromosome evolution in D. miranda Comparative genomic and transcriptomic analyses revealed that the pseudogenization rate on the neo-X is also accelerated in D. albomicans and D. americana although to a lesser extent than in D. miranda In males, neo-X-linked genes whose neo-Y-linked homologs are pseudogenized tended to be up-regulated more than those whose neo-Y-linked homologs remain functional. Moreover, genes under strong functional constraint and genes highly expressed in the testis tended to remain functional on the neo-X and neo-Y, respectively. Focusing on the D. miranda and D. albomicans neo-sex chromosomes that emerged independently from the same autosome, we further found that the same genes tend to become pseudogenized in parallel on the neo-Y. These genes include Idgf6 and JhI-26, which may be unnecessary or even harmful in males. Our results indicate that neo-sex chromosomes in Drosophila share a common evolutionary trajectory after their emergence, which may prevent sex chromosomes from being an evolutionary dead end.
| S-EPMC8559708 | biostudies-literature
Project description:BackgroundClustering is a fundamental operation in the analysis of biological sequence data. New DNA sequencing technologies have dramatically increased the rate at which we can generate data, resulting in datasets that cannot be efficiently analyzed by traditional clustering methods.This is particularly true in the context of taxonomic profiling of microbial communities through direct sequencing of phylogenetic markers (e.g. 16S rRNA) - the domain that motivated the work described in this paper. Many analysis approaches rely on an initial clustering step aimed at identifying sequences that belong to the same operational taxonomic unit (OTU). When defining OTUs (which have no universally accepted definition), scientists must balance a trade-off between computational efficiency and biological accuracy, as accurately estimating an environment's phylogenetic composition requires computationally-intensive analyses. We propose that efficient and mathematically well defined clustering methods can benefit existing taxonomic profiling approaches in two ways: (i) the resulting clusters can be substituted for OTUs in certain applications; and (ii) the clustering effectively reduces the size of the data-sets that need to be analyzed by complex phylogenetic pipelines (e.g., only one sequence per cluster needs to be provided to downstream analyses).ResultsTo address the challenges outlined above, we developed DNACLUST, a fast clustering tool specifically designed for clustering highly-similar DNA sequences.Given a set of sequences and a sequence similarity threshold, DNACLUST creates clusters whose radius is guaranteed not to exceed the specified threshold. Underlying DNACLUST is a greedy clustering strategy that owes its performance to novel sequence alignment and k-mer based filtering algorithms.DNACLUST can also produce multiple sequence alignments for every cluster, allowing users to manually inspect clustering results, and enabling more detailed analyses of the clustered data.ConclusionsWe compare DNACLUST to two popular clustering tools: CD-HIT and UCLUST. We show that DNACLUST is about an order of magnitude faster than CD-HIT and UCLUST (exact mode) and comparable in speed to UCLUST (approximate mode). The performance of DNACLUST improves as the similarity threshold is increased (tight clusters) making it well suited for rapidly removing duplicates and near-duplicates from a dataset, thereby reducing the size of the data being analyzed through more elaborate approaches.
| S-EPMC3213679 | biostudies-literature
Project description:The characterization of functional redundancy and divergence between duplicate genes is an important step in understanding the evolution of genetic systems. Large-scale genetic network analysis in Saccharomyces cerevisiae provides a powerful perspective for addressing these questions through quantitative measurements of genetic interactions between pairs of duplicated genes, and more generally, through the study of genome-wide genetic interaction profiles associated with duplicated genes. We show that duplicate genes exhibit fewer genetic interactions than other genes because they tend to buffer one another functionally, whereas observed interactions are non-overlapping and reflect their divergent roles. We also show that duplicate gene pairs are highly imbalanced in their number of genetic interactions with other genes, a pattern that appears to result from asymmetric evolution, such that one duplicate evolves or degrades faster than the other and often becomes functionally or conditionally specialized. The differences in genetic interactions are predictive of differences in several other evolutionary and physiological properties of duplicate pairs.
| S-EPMC3010121 | biostudies-literature
Project description:Phylogenic association between bacteria living under harsh conditions can provide important information on adaptive mechanism, survival strategy and their potential application. Indigenous microorganisms isolated from toxic refinery oily sludge with ability to degrade a diverse range of hydrocarbons were identified and characterized. The strains including Pseudomonas aeruginosa RS1, Microbacterium sp. RS2, Bacillus sp. RS3, Acinetobacter baumannii RS4 and Stenotrophomonas sp. RS5 could utilize n-alkanes, cycloalkanes, polynuclear aromatic hydrocarbons (PAHs) with 2-4 rings and also substituted PAHs as sole substrate. The phylogenetic position of Bacillus sp. RS3 and Pseudomonas sp. RS1 was tested by applying the maximum likelihood (ML) method to the aligned 16S rRNA nucleotide sequences of PAH and aliphatic hydrocarbon degrading strains belonging to the corresponding genus. The base substitution matrix created with each set of organisms capable of degrading aromatic and aliphatic hydrocarbons showed significant transitional event with high values of transition: transversion ratio (R) under all conditions. The guanine-cytosine (GC) content of the hydrocarbon degrading test strains was also found to be highest for the clade which harbored them. The test strains consistently occupied a distinct terminal end within the phylogenetic tree constructed by ML analysis. This study reveals that the refinery sludge imposed environmental stress on the bacterial strains which possibly caused significant genetic alteration and phenotypic adaptation. Due to the divergent evolution of the Pseudomonas and Bacillus strains in the sludge, they appeared distinctly different from other hydrocarbon degrading strains of the same genus.
| S-EPMC5971019 | biostudies-literature
Project description:Human isolates of serotype III Streptococcus agalactiae (group B streptococcus [GBS]) can be divided into three separate phylogenetic lineages based on analysis of the restriction digest patterns (RDPs) of chromosomal DNA. Nine DNA sequences that are present in all isolates of the RDP III-3 phylogenetic lineage, but not in the other lineages, were identified by genomic subtractive hybridization. A complete physical map of a III-3 chromosome was constructed. Six of the nine III-3-specific sequences mapped to a 340-kb Sse8387I fragment which contains or is located close to known GBS virulence genes. One of the III-3-specific probes, AW-10, encodes part of GBSi1, a group II intron that is inserted at two sites within the GBS genome. The second chromosomal site for GBSi1 was isolated, sequenced, and mapped to a location near the locus responsible for hemolysin production. These findings suggest that the genetic variation that distinguishes the RDP type III-3 strains from other serotype III strains occurs largely within localized areas of the genome containing known or putative virulence genes.
| S-EPMC127635 | biostudies-literature
Project description:BackgroundWnt5 genes belong to the large Wnt family, encoding proteins implicated into several tumorigenic and developmental processes. Phylogenetic analyses showed that Wnt5 gene has been duplicated at the divergence time of gnathostomata from agnatha. Interestingly, experimental data for some species indicated that only one of the two Wnt5 paralogs participates in the development of the endocrine pancreas. The purpose of this paper is to reexamine the phylogenetic history of the Wnt5 developmental regulators and investigate the functional shift between paralogs through comparative genomics.ResultsIn this study, the phylogeny of Wnt5 genes was investigated in species belonging to protostomia and deuterostomia. Furthermore, an in silico regulatory region analysis of Wnt5 paralogs was conducted, limited to those species with insulin producing cells and pancreas, covering the evolutionary distance from agnatha to gnathostomata. Our results confirmed the Wnt5 gene duplication and additionally revealed that this duplication event included also the upstream region. Moreover, within this latter region, a conserved module was detected to which a complex of transcription factors, known to be implicated in embryonic pancreas formation, bind.ConclusionsResults and observations presented in this study, allow us to conclude that during evolution, the Wnt5 gene has been duplicated in early vertebrates, and that some paralogs conserved a module within their regulatory region, functionally related to embryonic development of pancreas. Interestingly, our results allowed advancing a possible explanation on why the Wnt5 orthologs do not share the same function during pancreas development. As a final remark, we suggest that an in silico comparative analysis of regulatory regions, especially when associated to published experimental data, represents a powerful approach for explaining shift of roles among paralogs.
| S-EPMC2922100 | biostudies-literature