Applications of Evolutionary Bioinformatics in Basic and Biomedical Research

With the revolutionary progress in sequencing technologies, computational biology emerged as a game-changing field which is applied in understanding molecular events of life for not only complementary but also exploratory purposes. Bioinformatics resources and tools significantly help in data generation, organization and analysis. However, there is still a need for developing new approaches built based on a biologist’s point of view. In protein bioinformatics, there are several fundamental problems such as (i) determining protein function; (ii) identifying protein-protein interactions; (iii) predicting the effect of amino acid variants. Here, I present three chapters addressing these problems from an evolutionary perspective. Firstly, I describe a novel search pipeline for protein domain identification. The algorithm chain provides sensitive domain assignments with the highest possible specificity. Secondly, I present a tool enabling large-scale visualization of presences and absences of proteins in hierarchically clustered genomes. This tool visualizes multi-layer information of any kind of genome-linked data with a special focus on domain architectures, enabling identification of coevolving domains/proteins, which can eventually help in identifying functionally interacting proteins. And finally, I propose an approach for distinguishing between benign and damaging missense mutations in a human disease by establishing the precise evolutionary history of the associated gene. This part introduces new criteria on how to determine functional orthologs via phylogenetic analysis. All three parts use comparative genomics and/or sequence analyses. Taken together, this study addresses important problems in protein bioinformatics and as a whole it can be utilized to describe proteins by their domains, coevolving partners and functionally important residues

Dynamic maps of UV damage formation and repair for the human genome

Nucleotide excision repair removes DNA damage caused by carcinogens, such as UV and anticancer drugs, such as cisplatin. We have developed two methods, high-sensitivity damage sequencing and excision repair sequencing that map the formation and repair of damage in the human genome at single-nucleotide resolution. The combination of dynamic damage and repair maps provides a holistic perspective of UV damage and repair of the human genome and has potential applications in cancer prevention and chemotherapy

Aquerium: A web application for comparative exploration of domain-based protein occurrences on the taxonomically clustered genome tree: Architecture Querying Podium

Gene duplication and loss are major driving forces in evolution. While many important genomic resources provide information on gene presence, there is a lack of tools giving equal importance to presence and absence information as well as web platforms enabling easy visual comparison of multiple domain-based protein occurences at once. Here, we present Aquerium, a platform for visualizing genomic presence and absence of biomolecules with a focus on protein domain architectures. The web server offers advanced domain organization querying against the database of pre-computed domains for ~26000 organisms and it can be utilized for identification of evolutionary events, such as fusion, disassociation, duplication and shuffling of protein domains. The tool also allows alternative inputs of custom entries or BLASTP results for visualization. Aquerium will be a useful tool for biologists who perform comparative genomic and evolutionary analyses. The web server is freely accessible at http://aquerium.utk.edu

Differential damage and repair of DNA-adducts induced by anti-cancer drug cisplatin across mouse organs

The platinum-based drug cisplatin is a widely used first-line therapy for several cancers. Cisplatin interacts with DNA mainly in the form of Pt-d(GpG) di-adduct, which stalls cell proliferation and activates DNA damage response. Although cisplatin shows a broad spectrum of anticancer activity, its utility is limited due to acquired drug resistance and toxicity to non-targeted tissues. Here, by integrating genome-wide high-throughput Damage-seq, XR-seq, and RNA-seq approaches, along with publicly available epigenomic data, we systematically study the genome-wide profiles of cisplatin damage formation and excision repair in mouse kidney, liver, lung and spleen. We find different DNA damage and repair spectra across mouse organs, which are associated with tissue-specific transcriptomic and epigenomic profiles. The framework and the multi-omics data we present here constitute an unbiased foundation for understanding the mechanisms of cellular response to cisplatin. Our approach should be applicable for studying drug resistance and for tailoring cancer chemotherapy regimens

Genome-wide transcription-coupled repair in Escherichia coli is mediated by the Mfd translocase

In transcription-coupled repair (TCR), nucleotide excision repair occurs most rapidly in the template strand of actively transcribed genes. TCR has been observed in a limited set of genes directly assayed in Escherichia coli cells. In vitro, Mfd translocase performs reactions necessary to mediate TCR: It removes RNA polymerase blocked by a template strand lesion and rapidly delivers repair enzymes to the lesion. This study applied excision repair sequencing methodology to map the location of repair sites in different E. coli strains. Results showed that Mfd-dependent TCR is widespread in the E. coli genome. Results with UvrD helicase demonstrated its role in basal repair, but no overall role in TCR

MiST 3.0: an updated microbial signal transduction database with an emphasis on chemosensory systems

Bacteria and archaea employ dedicated signal transduction systems that modulate gene expression, second-messenger turnover, quorum sensing, biofilm formation, motility, host-pathogen and beneficial interactions. The updated MiST database provides a comprehensive classification of microbial signal transduction systems. This update is a result of a substantial scaling to accommodate constantly growing microbial genomic data. More than 125 000 genomes, 516 million genes and almost 100 million unique protein sequences are currently stored in the database. For each bacterial and archaeal genome, MiST 3.0 provides a complete signal transduction profile, thus facilitating theoretical and experimental studies on signal transduction and gene regulation. New software infrastructure and distributed pipeline implemented in MiST 3.0 enable regular genome updates based on the NCBI RefSeq database. A novel MiST feature is the integration of unique profile HMMs to link complex chemosensory systems with corresponding chemoreceptors in bacterial and archaeal genomes. The data can be explored online or via RESTful API (freely available at https://mistdb.com)

Effects of replication domains on genome-wide UV-induced DNA damage and repair

Nucleotide excision repair is the primary repair mechanism that removes UV-induced DNA lesions in placentals. Unrepaired UV-induced lesions could result in mutations during DNA replication. Although the mutagenesis of pyrimidine dimers is reasonably well understood, the direct effects of replication fork progression on nucleotide excision repair are yet to be clarified. Here, we applied Damage-seq and XR-seq techniques and generated replication maps in synchronized UV-treated HeLa cells. The results suggest that ongoing replication stimulates local repair in both early and late replication domains. Additionally, it was revealed that lesions on lagging strand templates are repaired slower in late replication domains, which is probably due to the imbalanced sequence context. Asymmetric relative repair is in line with the strand bias of melanoma mutations, suggesting a role of exogenous damage, repair, and replication in mutational strand asymmetry

Establishing the precise evolutionary history of a gene improves prediction of disease-causing missense mutations

PURPOSE: Predicting the phenotypic effects of mutations has become an important application in clinical genetic diagnostics. Computational tools evaluate the behavior of the variant over evolutionary time and assume that variations seen during the course of evolution are probably benign in humans. However, current tools do not take into account orthologous/paralogous relationships. Paralogs have dramatically different roles in Mendelian diseases. For example, whereas inactivating mutations in the NPC1 gene cause the neurodegenerative disorder Niemann-Pick C, inactivating mutations in its paralog NPC1L1 are not disease-causing and, moreover, are implicated in protection from coronary heart disease. METHODS: We identified major events in NPC1 evolution and revealed and compared orthologs and paralogs of the human NPC1 gene through phylogenetic and protein sequence analyses. We predicted whether an amino acid substitution affects protein function by reducing the organism’s fitness. RESULTS: Removing the paralogs and distant homologs improved the overall performance of categorizing disease-causing and benign amino acid substitutions. CONCLUSION: The results show that a thorough evolutionary analysis followed by identification of orthologs improves the accuracy in predicting disease-causing missense mutations. We anticipate that this approach will be used as a reference in the interpretation of variants in other genetic diseases as well. Genet Med 18 10, 1029–1036

bioliners/projects: MiST v4.0

This release includes: incorporation of metagenome-assembled genomes into MiST their integration with genomes within a single website biosample information scaled representation of proteins taxonomy updates due to the recent inclusion of the rank "phylum" in the International Code of Nomenclature for Prokaryote, which effected 42 phyl