1
|
Yao HT, Marchand B, Berkemer SJ, Ponty Y, Will S. Infrared: a declarative tree decomposition-powered framework for bioinformatics. Algorithms Mol Biol 2024; 19:13. [PMID: 38493130 PMCID: PMC10943887 DOI: 10.1186/s13015-024-00258-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 02/13/2024] [Indexed: 03/18/2024] Open
Abstract
MOTIVATION Many bioinformatics problems can be approached as optimization or controlled sampling tasks, and solved exactly and efficiently using Dynamic Programming (DP). However, such exact methods are typically tailored towards specific settings, complex to develop, and hard to implement and adapt to problem variations. METHODS We introduce the Infrared framework to overcome such hindrances for a large class of problems. Its underlying paradigm is tailored toward problems that can be declaratively formalized as sparse feature networks, a generalization of constraint networks. Classic Boolean constraints specify a search space, consisting of putative solutions whose evaluation is performed through a combination of features. Problems are then solved using generic cluster tree elimination algorithms over a tree decomposition of the feature network. Their overall complexities are linear on the number of variables, and only exponential in the treewidth of the feature network. For sparse feature networks, associated with low to moderate treewidths, these algorithms allow to find optimal solutions, or generate controlled samples, with practical empirical efficiency. RESULTS Implementing these methods, the Infrared software allows Python programmers to rapidly develop exact optimization and sampling applications based on a tree decomposition-based efficient processing. Instead of directly coding specialized algorithms, problems are declaratively modeled as sets of variables over finite domains, whose dependencies are captured by constraints and functions. Such models are then automatically solved by generic DP algorithms. To illustrate the applicability of Infrared in bioinformatics and guide new users, we model and discuss variants of bioinformatics applications. We provide reimplementations and extensions of methods for RNA design, RNA sequence-structure alignment, parsimony-driven inference of ancestral traits in phylogenetic trees/networks, and design of coding sequences. Moreover, we demonstrate multidimensional Boltzmann sampling. These applications of the framework-together with our novel results-underline the practical relevance of Infrared. Remarkably, the achieved complexities are typically equivalent to the ones of specialized algorithms and implementations. AVAILABILITY Infrared is available at https://amibio.gitlabpages.inria.fr/Infrared with extensive documentation, including various usage examples and API reference; it can be installed using Conda or from source.
Collapse
Affiliation(s)
- Hua-Ting Yao
- LIX, CNRS UMR 7161, Ecole Polytechnique, Institut Polytechnique de Paris, Palaiseau, France.
- Department of Theoretical Chemistry, University of Vienna, Vienna, Austria.
- School of Computer Science, McGill University, Montreal, Canada.
| | - Bertrand Marchand
- LIX, CNRS UMR 7161, Ecole Polytechnique, Institut Polytechnique de Paris, Palaiseau, France
| | - Sarah J Berkemer
- LIX, CNRS UMR 7161, Ecole Polytechnique, Institut Polytechnique de Paris, Palaiseau, France
- Earth-Life Science Institute, Tokyo Institute of Technology, Tokyo, Japan
| | - Yann Ponty
- LIX, CNRS UMR 7161, Ecole Polytechnique, Institut Polytechnique de Paris, Palaiseau, France
| | - Sebastian Will
- LIX, CNRS UMR 7161, Ecole Polytechnique, Institut Polytechnique de Paris, Palaiseau, France.
| |
Collapse
|
2
|
Eggenhofer F, Höner Zu Siederdissen C. Evolutionary Structure Conservation and Covariance Scores. Methods Mol Biol 2024; 2726:255-284. [PMID: 38780735 DOI: 10.1007/978-1-0716-3519-3_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
Effective homology search for non-coding RNAs is frequently not possible via sequence similarity alone. Current methods leverage evolutionary information like structure conservation or covariance scores to identify homologs in organisms that are phylogenetically more distant. In this chapter, we introduce the theoretical background of evolutionary structure conservation and covariance score, and we show hands-on how current methods in the field are applied on example datasets.
Collapse
Affiliation(s)
- Florian Eggenhofer
- Bioinformatics Group, Department of Computer Science University of Freiburg, Freiburg, Germany
| | - Christian Höner Zu Siederdissen
- Bioinformatics Group, Department of Computer Science, University of Leipzig, Leipzig, Germany.
- Interdisciplinary Center for Bioinformatics, University of Leipzig, Leipzig, Germany.
- Bioinformatics/High-Throughput Analysis, Faculty of Mathematics and Computer Science, Friedrich Schiller University Jena, Jena, Germany.
| |
Collapse
|
3
|
Kinateder T, Drexler L, Straub K, Merkl R, Sterner R. Experimental and computational analysis of the ancestry of an evolutionary young enzyme from histidine biosynthesis. Protein Sci 2023; 32:e4536. [PMID: 36502290 PMCID: PMC9798254 DOI: 10.1002/pro.4536] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 12/01/2022] [Accepted: 12/03/2022] [Indexed: 12/14/2022]
Abstract
The conservation of fold and chemistry of the enzymes associated with histidine biosynthesis suggests that this pathway evolved prior to the diversification of Bacteria, Archaea, and Eukaryotes. The only exception is the histidinol phosphate phosphatase (HolPase). So far, non-homologous HolPases that possess distinct folds and belong to three different protein superfamilies have been identified in various phylogenetic clades. However, their evolution has remained unknown to date. Here, we analyzed the evolutionary history of the HolPase from γ-Proteobacteria (HisB-N). It has been argued that HisB-N and its closest homologue d-glycero-d-manno-heptose-1,7-bisphosphate 7-phosphatase (GmhB) have emerged from the same promiscuous ancestral phosphatase. GmhB variants catalyze the hydrolysis of the anomeric d-glycero-d-manno-heptose-1,7-bisphosphate (αHBP or βHBP) with a strong preference for one anomer (αGmhB or βGmhB). We found that HisB-N from Escherichia coli shows promiscuous activity for βHBP but not αHBP, while βGmhB from Crassaminicella sp. shows promiscuous activity for HolP. Accordingly, a combined phylogenetic tree of αGmhBs, βGmhBs, and HisB-N sequences revealed that HisB-Ns form a compact subcluster derived from βGmhBs. Ancestral sequence reconstruction and in vitro analysis revealed a promiscuous HolPase activity in the resurrected enzymes prior to functional divergence of the successors. The following increase in catalytic efficiency of the HolP turnover is reflected in the shape and electrostatics of the active site predicted by AlphaFold. An analysis of the phylogenetic tree led to a revised evolutionary model that proposes the horizontal gene transfer of a promiscuous βGmhB from δ- to γ-Proteobacteria where it evolved to the modern HisB-N.
Collapse
Affiliation(s)
- Thomas Kinateder
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| | - Lukas Drexler
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| | - Kristina Straub
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| | - Rainer Merkl
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| | - Reinhard Sterner
- Institute of Biophysics and Physical Biochemistry and Regensburg Center for Biochemistry, University of RegensburgRegensburgGermany
| |
Collapse
|
4
|
Kuang M, Zhang Y, Lam TW, Ting HF. MLProbs: A Data-Centric Pipeline for Better Multiple Sequence Alignment. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2023; 20:524-533. [PMID: 35120007 DOI: 10.1109/tcbb.2022.3148382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
In this paper, we explore using the data-centric approach to tackle the Multiple Sequence Alignment (MSA) construction problem. Unlike the algorithm-centric approach, which reduces the construction problem to a combinatorial optimization problem based on an abstract mathematical model, the data-centric approach explores using classification models trained from existing benchmark data to guide the construction. We identified two simple classifications to help us choose a better alignment tool and determine whether and how much to carry out realignment. We show that shallow machine-learning algorithms suffice to train sensitive models for these classifications. Based on these models, we implemented a new multiple sequence alignment pipeline, called MLProbs. Compared with 10 other popular alignment tools over four benchmark databases (namely, BAliBASE, OXBench, OXBench-X and SABMark), MLProbs consistently gives the highest TC score. More importantly, MLProbs shows non-trivial improvement for protein families with low similarity; in particular, when evaluated against the 1,356 protein families with similarity ≤ 50%, MLProbs achieves a TC score of 56.93, while the next best three tools are in the range of [55.41, 55.91] (increased by more than 1.8%). We also compared the performance of MLProbs and other MSA tools in two real-life applications - Phylogenetic Tree Construction Analysis and Protein Secondary Structure Prediction - and MLProbs also had the best performance. In our study, we used only shallow machine-learning algorithms to train our models. It would be interesting to study whether deep-learning methods can help make further improvements, so we suggest some possible research directions in the conclusion section.
Collapse
|
5
|
Molecular phylogenetics and systematics of two enteric helminth parasites (Baylisascaris laevis and Diandrya vancouverensis) in the Vancouver Island marmot (Marmota vancouverensis). INTERNATIONAL JOURNAL FOR PARASITOLOGY: PARASITES AND WILDLIFE 2022; 19:301-310. [DOI: 10.1016/j.ijppaw.2022.11.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2022] [Revised: 11/15/2022] [Accepted: 11/15/2022] [Indexed: 11/18/2022]
|
6
|
Marchant DB, Chen G, Cai S, Chen F, Schafran P, Jenkins J, Shu S, Plott C, Webber J, Lovell JT, He G, Sandor L, Williams M, Rajasekar S, Healey A, Barry K, Zhang Y, Sessa E, Dhakal RR, Wolf PG, Harkess A, Li FW, Rössner C, Becker A, Gramzow L, Xue D, Wu Y, Tong T, Wang Y, Dai F, Hua S, Wang H, Xu S, Xu F, Duan H, Theißen G, McKain MR, Li Z, McKibben MTW, Barker MS, Schmitz RJ, Stevenson DW, Zumajo-Cardona C, Ambrose BA, Leebens-Mack JH, Grimwood J, Schmutz J, Soltis PS, Soltis DE, Chen ZH. Dynamic genome evolution in a model fern. NATURE PLANTS 2022; 8:1038-1051. [PMID: 36050461 PMCID: PMC9477723 DOI: 10.1038/s41477-022-01226-7] [Citation(s) in RCA: 43] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Accepted: 07/15/2022] [Indexed: 05/31/2023]
Abstract
The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.
Collapse
Affiliation(s)
| | - Guang Chen
- Central Laboratory, State Key Laboratory for Managing Biotic and Chemical Threats to the Quality and Safety of Agro-products, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
- College of Agriculture, Yangtze University, Jingzhou, China
| | - Shengguan Cai
- College of Agriculture and Biotechnology, Zhejiang University, Hangzhou, China
- School of Science, Western Sydney University, Penrith, New South Wales, Australia
| | - Fei Chen
- College of Life and Environmental Sciences, Hangzhou Normal University, Hangzhou, China
| | | | - Jerry Jenkins
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Shengqiang Shu
- United States Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Chris Plott
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Jenell Webber
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - John T Lovell
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
- United States Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Guifen He
- United States Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Laura Sandor
- United States Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Melissa Williams
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Shanmugam Rajasekar
- Arizona Genomics Institute, School of Plant Sciences, University of Arizona, Tucson, AZ, USA
| | - Adam Healey
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Kerrie Barry
- United States Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Yinwen Zhang
- Institute of Bioinformatics, University of Georgia, Athens, GA, USA
| | - Emily Sessa
- Department of Biology, University of Florida, Gainesville, FL, USA
| | - Rijan R Dhakal
- Department of Biological Sciences, University of Alabama in Huntsville, Huntsville, AL, USA
| | - Paul G Wolf
- Department of Biological Sciences, University of Alabama in Huntsville, Huntsville, AL, USA
| | - Alex Harkess
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
- Department of Crop, Soil, and Environmental Sciences, Auburn University, Auburn, AL, USA
| | - Fay-Wei Li
- Boyce Thompson Institute, Ithaca, NY, USA
- Plant Biology Section, Cornell University, Ithaca, NY, USA
| | - Clemens Rössner
- Justus-Liebig-University, Department of Biology and Chemistry, Institute of Botany, Gießen, Germany
| | - Annette Becker
- Justus-Liebig-University, Department of Biology and Chemistry, Institute of Botany, Gießen, Germany
| | - Lydia Gramzow
- Matthias Schleiden Institute/Genetics, Friedrich Schiller University Jena, Jena, Germany
| | - Dawei Xue
- College of Life and Environmental Sciences, Hangzhou Normal University, Hangzhou, China
| | - Yuhuan Wu
- College of Life and Environmental Sciences, Hangzhou Normal University, Hangzhou, China
| | - Tao Tong
- College of Agriculture, Yangtze University, Jingzhou, China
| | - Yuanyuan Wang
- Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan, China
| | - Fei Dai
- College of Agriculture and Biotechnology, Zhejiang University, Hangzhou, China
| | - Shuijin Hua
- Institute of Crops and Nuclear Technology Utilization, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
| | - Hua Wang
- State Key Laboratory for Managing Biotic and Chemical Threats to the Quality and Safety of Agro-products, Institute of Virology and Biotechnology, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
| | - Shengchun Xu
- College of Agriculture and Biotechnology, Zhejiang University, Hangzhou, China
| | - Fei Xu
- College of Agriculture and Biotechnology, Zhejiang University, Hangzhou, China
| | - Honglang Duan
- Institute for Forest Resources & Environment of Guizhou, Key Laboratory of Forest Cultivation in Plateau Mountain of Guizhou Province, College of Forestry, Guizhou University, Guiyang, China
| | - Günter Theißen
- Matthias Schleiden Institute/Genetics, Friedrich Schiller University Jena, Jena, Germany
| | - Michael R McKain
- Department of Biological Sciences, University of Alabama, Tuscaloosa, AL, USA
| | - Zheng Li
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, USA
| | - Michael T W McKibben
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ, USA
| | - Michael S Barker
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ, USA
| | | | | | | | | | | | - Jane Grimwood
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - Jeremy Schmutz
- Genome Sequencing Center, HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
- United States Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Pamela S Soltis
- Florida Museum of Natural History, University of Florida, Gainesville, FL, USA.
| | - Douglas E Soltis
- Department of Biology, University of Florida, Gainesville, FL, USA.
- Florida Museum of Natural History, University of Florida, Gainesville, FL, USA.
| | - Zhong-Hua Chen
- School of Science, Western Sydney University, Penrith, New South Wales, Australia.
- Hawkesbury Institute for the Environment, Western Sydney University, Penrith, New South Wales, Australia.
| |
Collapse
|
7
|
Hassler HB, Probert B, Moore C, Lawson E, Jackson RW, Russell BT, Richards VP. Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies. MICROBIOME 2022; 10:104. [PMID: 35799218 PMCID: PMC9264627 DOI: 10.1186/s40168-022-01295-y] [Citation(s) in RCA: 44] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Accepted: 05/23/2022] [Indexed: 05/02/2023]
Abstract
BACKGROUND The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82]. RESULTS At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%). CONCLUSIONS Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.
Collapse
Affiliation(s)
- Hayley B. Hassler
- Department of Biological Sciences, College of Science, Clemson University, Clemson, SC 29634 USA
| | - Brett Probert
- Department of Biological Sciences, College of Science, Clemson University, Clemson, SC 29634 USA
| | - Carson Moore
- Department of Biological Sciences, College of Science, Clemson University, Clemson, SC 29634 USA
| | - Elizabeth Lawson
- Department of Biological Sciences, College of Science, Clemson University, Clemson, SC 29634 USA
| | | | - Brook T. Russell
- School of Mathematical and Statistical Sciences, Clemson University, Clemson, SC 29634 USA
| | - Vincent P. Richards
- Department of Biological Sciences, College of Science, Clemson University, Clemson, SC 29634 USA
| |
Collapse
|
8
|
Nayeem MA, Bayzid MS, Rahman AH, Shahriyar R, Rahman MS. Multiobjective Formulation of Multiple Sequence Alignment for Phylogeny Inference. IEEE TRANSACTIONS ON CYBERNETICS 2022; 52:2775-2786. [PMID: 33044939 DOI: 10.1109/tcyb.2020.3020308] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Multiple sequence alignment (MSA) is a preliminary task for estimating phylogenies. It is used for homology inference among the sequences of a set of species. Generally, the MSA task is handled as a single-objective optimization process. The alignments computed under one criterion may be different from the alignments generated by other criteria, inferring discordant homologies and thus leading to different hypothesized evolutionary histories relating the sequences. The multiobjective (MO) formulation of MSA has recently been advocated by several researchers, to address this issue. An MO approach independently optimizes multiple (often conflicting) objective functions at the same time and outputs a set of competitive alignments. However, no conceptual or experimental rational from a real-world application perspective has been reported so far for any MO formulation of MSA. This article work investigates the impact of MO formulation in the context of an important scientific problem, namely, phylogeny estimation. Employing popular evolutionary MO algorithms, we show that: 1) trees inferred based on alignments produced by the existing MSA methods used in practice are substantially worse in quality than the trees inferred based on the alignment's output by an MO algorithm and 2) even high-quality alignments (according to popular measures available in the literature) may fail to achieve acceptable accuracy in generating phylogenetic trees. Thus, we essentially ask the following natural question: "can a phylogeny-aware (i.e., application-aware) metric guide in selecting appropriate MO formulations to ensure better phylogeny estimation?" Here, we report a carefully designed extensive experimental study that positively answers this question.
Collapse
|
9
|
Chao J, Tang F, Xu L. Developments in Algorithms for Sequence Alignment: A Review. Biomolecules 2022; 12:biom12040546. [PMID: 35454135 PMCID: PMC9024764 DOI: 10.3390/biom12040546] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2022] [Revised: 03/29/2022] [Accepted: 03/31/2022] [Indexed: 01/27/2023] Open
Abstract
The continuous development of sequencing technologies has enabled researchers to obtain large amounts of biological sequence data, and this has resulted in increasing demands for software that can perform sequence alignment fast and accurately. A number of algorithms and tools for sequence alignment have been designed to meet the various needs of biologists. Here, the ideas that prevail in the research of sequence alignment and some quality estimation methods for multiple sequence alignment tools are summarized.
Collapse
Affiliation(s)
- Jiannan Chao
- Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, Chengdu 610054, China;
| | - Furong Tang
- Yangtze Delta Region Institute (Quzhou), University of Electronic Science and Technology of China, Quzhou 324003, China;
- School of Electronic and Communication Engineering, Shenzhen Polytechnic, Shenzhen 518055, China
| | - Lei Xu
- School of Electronic and Communication Engineering, Shenzhen Polytechnic, Shenzhen 518055, China
- Correspondence:
| |
Collapse
|
10
|
Zhang Y, Zhang Q, Zhou J, Zou Q. A survey on the algorithm and development of multiple sequence alignment. Brief Bioinform 2022; 23:6546258. [PMID: 35272347 DOI: 10.1093/bib/bbac069] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2021] [Revised: 01/30/2022] [Accepted: 02/09/2022] [Indexed: 12/21/2022] Open
Abstract
Multiple sequence alignment (MSA) is an essential cornerstone in bioinformatics, which can reveal the potential information in biological sequences, such as function, evolution and structure. MSA is widely used in many bioinformatics scenarios, such as phylogenetic analysis, protein analysis and genomic analysis. However, MSA faces new challenges with the gradual increase in sequence scale and the increasing demand for alignment accuracy. Therefore, developing an efficient and accurate strategy for MSA has become one of the research hotspots in bioinformatics. In this work, we mainly summarize the algorithms for MSA and its applications in bioinformatics. To provide a structured and clear perspective, we systematically introduce MSA's knowledge, including background, database, metric and benchmark. Besides, we list the most common applications of MSA in the field of bioinformatics, including database searching, phylogenetic analysis, genomic analysis, metagenomic analysis and protein analysis. Furthermore, we categorize and analyze classical and state-of-the-art algorithms, divided into progressive alignment, iterative algorithm, heuristics, machine learning and divide-and-conquer. Moreover, we also discuss the challenges and opportunities of MSA in bioinformatics. Our work provides a comprehensive survey of MSA applications and their relevant algorithms. It could bring valuable insights for researchers to contribute their knowledge to MSA and relevant studies.
Collapse
Affiliation(s)
- Yongqing Zhang
- School of Computer Science, Chengdu University of Information Technology, 610225, Chengdu, China.,School of Computer Science and Engineering, University of Electronic Science and Technology of China, 611731, Chengdu, China
| | - Qiang Zhang
- School of Computer Science, Chengdu University of Information Technology, 610225, Chengdu, China
| | - Jiliu Zhou
- School of Computer Science, Chengdu University of Information Technology, 610225, Chengdu, China
| | - Quan Zou
- Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, 610054, Chengdu, China
| |
Collapse
|
11
|
Root ZD, Allen C, Gould C, Brewer M, Jandzik D, Medeiros DM. A Comprehensive Analysis of Fibrillar Collagens in Lamprey Suggests a Conserved Role in Vertebrate Musculoskeletal Evolution. Front Cell Dev Biol 2022; 10:809979. [PMID: 35242758 PMCID: PMC8887668 DOI: 10.3389/fcell.2022.809979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2021] [Accepted: 01/18/2022] [Indexed: 12/03/2022] Open
Abstract
Vertebrates have distinct tissues which are not present in invertebrate chordates nor other metazoans. The rise of these tissues also coincided with at least one round of whole-genome duplication as well as a suite of lineage-specific segmental duplications. Understanding whether novel genes lead to the origin and diversification of novel cell types, therefore, is of great importance in vertebrate evolution. Here we were particularly interested in the evolution of the vertebrate musculoskeletal system, the muscles and connective tissues that support a diversity of body plans. A major component of the musculoskeletal extracellular matrix (ECM) is fibrillar collagens, a gene family which has been greatly expanded upon in vertebrates. We thus asked whether the repertoire of fibrillar collagens in vertebrates reflects differences in the musculoskeletal system. To test this, we explored the diversity of fibrillar collagens in lamprey, a jawless vertebrate which diverged from jawed vertebrates (gnathostomes) more than five hundred million years ago and has undergone its own gene duplications. Some of the principal components of vertebrate hyaline cartilage are the fibrillar collagens type II and XI, but their presence in cartilage development across all vertebrate taxa has been disputed. We particularly emphasized the characterization of genes in the lamprey hyaline cartilage, testing if its collagen repertoire was similar to that in gnathostomes. Overall, we discovered thirteen fibrillar collagens from all known gene subfamilies in lamprey and were able to identify several lineage-specific duplications. We found that, while the collagen loci have undergone rearrangement, the Clade A genes have remained linked with the hox clusters, a phenomenon also seen in gnathostomes. While the lamprey muscular tissue was largely similar to that seen in gnathostomes, we saw considerable differences in the larval lamprey skeletal tissue, with distinct collagen combinations pertaining to different cartilage types. Our gene expression analyses were unable to identify type II collagen in the sea lamprey hyaline cartilage nor any other fibrillar collagen during chondrogenesis at the stages observed, meaning that sea lamprey likely no longer require these genes during early cartilage development. Our findings suggest that fibrillar collagens were multifunctional across the musculoskeletal system in the last common ancestor of vertebrates and have been largely conserved, but these genes alone cannot explain the origin of novel cell types.
Collapse
Affiliation(s)
- Zachary D Root
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, United States
| | - Cara Allen
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, United States
| | - Claire Gould
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, United States
| | - Margaux Brewer
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, United States
| | - David Jandzik
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, United States.,Department of Zoology, Comenius University in Bratislava, Bratislava, Slovakia
| | - Daniel M Medeiros
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, United States
| |
Collapse
|
12
|
Lladós J, Cores F, Guirado F, Lérida JL. Accurate consistency-based MSA reducing the memory footprint. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2021; 208:106237. [PMID: 34198017 DOI: 10.1016/j.cmpb.2021.106237] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Accepted: 06/08/2021] [Indexed: 06/13/2023]
Abstract
BACKGROUND AND OBJECTIVE The emergence of Next-Generation sequencing has created a push for faster and more accurate multiple sequence alignment tools. The growing number of sequences and their longer sizes, which require the use of increased system resources and produce less accurate results, are heavily challenging to these applications. Consistency-based methods have the most intensive CPU and memory usage requirements. We hypothesize that library reductions can enhance the scalability and performance of consistency-based multiple sequence alignment tools; however, we have previously shown a noticeable impact on the accuracy when extreme reductions were performed. METHODS In this study, we propose Matrix-Based T-Coffee, a consistency-based method that uses library reductions in conjunction with a complementary objective function. The proposed method, implemented in T-Coffee, can mitigate the accuracy loss caused by low memory resources. RESULTS The use of a complementary objective function with a library reduction of ≥30% improved the accuracy of T-Coffee. Interestingly, ≥50% library reduction achieved lower execution times and better overall scalability. CONCLUSIONS Matrix-Based T-Coffee benefits from accurate alignments while achieving better scalability. This leads to a reduction in memory footprint and execution time. In addition, these enhancements could be applied to other aligners based on consistency.
Collapse
Affiliation(s)
- Jordi Lladós
- INSPIRES Research Center, Universitat de Lleida. Jaume II. 69, 25001 Lleida, Spain.
| | - Fernando Cores
- INSPIRES Research Center, Universitat de Lleida. Jaume II. 69, 25001 Lleida, Spain.
| | - Fernando Guirado
- INSPIRES Research Center, Universitat de Lleida. Jaume II. 69, 25001 Lleida, Spain.
| | - Josep L Lérida
- INSPIRES Research Center, Universitat de Lleida. Jaume II. 69, 25001 Lleida, Spain.
| |
Collapse
|
13
|
Metaheuristics for multiple sequence alignment: A systematic review. Comput Biol Chem 2021; 94:107563. [PMID: 34425495 DOI: 10.1016/j.compbiolchem.2021.107563] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2020] [Revised: 08/04/2021] [Accepted: 08/09/2021] [Indexed: 11/21/2022]
Abstract
The Multiple Sequence Alignment (MSA) is a key task in bioinformatics, because it is used in different important biological analysis, such as function and structure prediction of unknown proteins. There are several approaches to perform MSA and the use of metaheuristics stands out because of the search ability of these methods, which generally leads to good results in a reasonable amount of time. This paper presents a Systematic Literature Review (SLR) on metaheuristics for MSA, compiling relevant works published between 2014 and 2019. The results of our SLR show the constant interest in this subject, due to the several recent publications that use different metaheuristics to obtain more accurate alignments. Moreover, the final results of our SLR show a multi-objective and hybrid approaches trends, which generally leads these methods to achieve even better results. Thus, we show in this work how the use of metaheuristics to perform MSA still remains an important and promising open research field.
Collapse
|
14
|
Rehman HA, Zafar K, Khan A, Imtiaz A. Multiple sequence alignment using enhanced bird swarm align algorithm. JOURNAL OF INTELLIGENT & FUZZY SYSTEMS 2021. [DOI: 10.3233/jifs-210055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Abstract
Discovering structural, functional and evolutionary information in biological sequences have been considered as a core research area in Bioinformatics. Multiple Sequence Alignment (MSA) tries to align all sequences in a given query set to provide us ease in annotation of new sequences. Traditional methods to find the optimal alignment are computationally expensive in real time. This research presents an enhanced version of Bird Swarm Algorithm (BSA), based on bio inspired optimization. Enhanced Bird Swarm Align Algorithm (EBSAA) is proposed for multiple sequence alignment problem to determine the optimal alignment among different sequences. Twenty-one different datasets have been used in order to compare performance of EBSAA with Genetic Algorithm (GA) and Particle Swarm Align Algorithm (PSAA). The proposed technique results in better alignment as compared to GA and PSAA in most of the cases.
Collapse
Affiliation(s)
- Hafiz Asadul Rehman
- Department of Computer Science, NationalUniversity of Computer and Emerging Science Lahore, Pakistan
| | - Kashif Zafar
- Department of Computer Science, NationalUniversity of Computer and Emerging Science Lahore, Pakistan
| | - Ayesha Khan
- University of Management & Technology, Lahore, Pakistan
| | | |
Collapse
|
15
|
Jia H, Aadland K, Kolaczkowski O, Kolaczkowski B. Direct molecular evidence for an ancient, conserved developmental toolkit controlling post-transcriptional gene regulation in land plants. Mol Biol Evol 2021; 38:4765-4777. [PMID: 34196710 PMCID: PMC8557471 DOI: 10.1093/molbev/msab201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
In plants, miRNA production is orchestrated by a suite of proteins that control transcription of the pri-miRNA gene, post-transcriptional processing and nuclear export of the mature miRNA. Post-transcriptional processing of miRNAs is controlled by a pair of physically interacting proteins, hyponastic leaves 1 (HYL1) and Dicer-like 1 (DCL1). However, the evolutionary history and structural basis of the HYL1–DCL1 interaction is unknown. Here we use ancestral sequence reconstruction and functional characterization of ancestral HYL1 in vitro and in Arabidopsis thaliana to better understand the origin and evolution of the HYL1–DCL1 interaction and its impact on miRNA production and plant development. We found the ancestral plant HYL1 evolved high affinity for both double-stranded RNA (dsRNA) and its DCL1 partner before the divergence of mosses from seed plants (∼500 Ma), and these high-affinity interactions remained largely conserved throughout plant evolutionary history. Structural modeling and molecular binding experiments suggest that the second of two dsRNA-binding motifs (DSRMs) in HYL1 may interact tightly with the first of two C-terminal DCL1 DSRMs to mediate the HYL1–DCL1 physical interaction necessary for efficient miRNA production. Transgenic expression of the nearly 200 Ma-old ancestral flowering-plant HYL1 in A. thaliana was sufficient to rescue many key aspects of plant development disrupted by HYL1− knockout and restored near-native miRNA production, suggesting that the functional partnership of HYL1–DCL1 originated very early in and was strongly conserved throughout the evolutionary history of terrestrial plants. Overall, our results are consistent with a model in which miRNA-based gene regulation evolved as part of a conserved plant “developmental toolkit.”
Collapse
Affiliation(s)
- Haiyan Jia
- Department of Biology, University of North Carolina, Chapel Hill, NC
| | - Kelsey Aadland
- Department of Microbiology and Plant Pathology, University of California, Riverside, CA
| | - Oralia Kolaczkowski
- Department of Microbiology and Cell Science, Institute for Food and Agricultural Sciences, University of Florida, Gainesville, FL
| | - Bryan Kolaczkowski
- Department of Microbiology and Cell Science, Institute for Food and Agricultural Sciences, University of Florida, Gainesville, FL
| |
Collapse
|
16
|
Aadland K, Kolaczkowski B. Alignment-Integrated Reconstruction of Ancestral Sequences Improves Accuracy. Genome Biol Evol 2021; 12:1549-1565. [PMID: 32785673 PMCID: PMC7523730 DOI: 10.1093/gbe/evaa164] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/03/2020] [Indexed: 12/31/2022] Open
Abstract
Ancestral sequence reconstruction (ASR) uses an alignment of extant protein sequences, a phylogeny describing the history of the protein family and a model of the molecular-evolutionary process to infer the sequences of ancient proteins, allowing researchers to directly investigate the impact of sequence evolution on protein structure and function. Like all statistical inferences, ASR can be sensitive to violations of its underlying assumptions. Previous studies have shown that, whereas phylogenetic uncertainty has only a very weak impact on ASR accuracy, uncertainty in the protein sequence alignment can more strongly affect inferred ancestral sequences. Here, we show that errors in sequence alignment can produce errors in ASR across a range of realistic and simplified evolutionary scenarios. Importantly, sequence reconstruction errors can lead to errors in estimates of structural and functional properties of ancestral proteins, potentially undermining the reliability of analyses relying on ASR. We introduce an alignment-integrated ASR approach that combines information from many different sequence alignments. We show that integrating alignment uncertainty improves ASR accuracy and the accuracy of downstream structural and functional inferences, often performing as well as highly accurate structure-guided alignment. Given the growing evidence that sequence alignment errors can impact the reliability of ASR studies, we recommend that future studies incorporate approaches to mitigate the impact of alignment uncertainty. Probabilistic modeling of insertion and deletion events has the potential to radically improve ASR accuracy when the model reflects the true underlying evolutionary history, but further studies are required to thoroughly evaluate the reliability of these approaches under realistic conditions.
Collapse
Affiliation(s)
- Kelsey Aadland
- Department of Microbiology and Cell Science, Institute of Food and Agricultural Sciences, University of Florida
| | - Bryan Kolaczkowski
- Department of Microbiology and Cell Science, Institute of Food and Agricultural Sciences, University of Florida
| |
Collapse
|
17
|
Lamprey lecticans link new vertebrate genes to the origin and elaboration of vertebrate tissues. Dev Biol 2021; 476:282-293. [PMID: 33887266 DOI: 10.1016/j.ydbio.2021.03.020] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2020] [Revised: 03/19/2021] [Accepted: 03/22/2021] [Indexed: 11/23/2022]
Abstract
The evolution of vertebrates from an invertebrate chordate ancestor involved the evolution of new organs, tissues, and cell types. It was also marked by the origin and duplication of new gene families. If, and how, these morphological and genetic innovations are related is an unresolved question in vertebrate evolution. Hyaluronan is an extracellular matrix (ECM) polysaccharide important for water homeostasis and tissue structure. Vertebrates possess a novel family of hyaluronan binding proteins called Lecticans, and studies in jawed vertebrates (gnathostomes) have shown they function in many of the cells and tissues that are unique to vertebrates. This raises the possibility that the origin and/or expansion of this gene family helped drive the evolution of these vertebrate novelties. In order to better understand the evolution of the lectican gene family, and its role in the evolution of vertebrate morphological novelties, we investigated the phylogeny, genomic arrangement, and expression patterns of all lecticans in the sea lamprey (Petromyzon marinus), a jawless vertebrate. Though both P. marinus and gnathostomes each have four lecticans, our phylogenetic and syntenic analyses are most consistent with the independent duplication of one of more lecticans in the lamprey lineage. Despite the likely independent expansion of the lamprey and gnathostome lectican families, we find highly conserved expression of lecticans in vertebrate-specific and mesenchyme-derived tissues. We also find that, unlike gnathostomes, lamprey expresses its lectican paralogs in distinct subpopulations of head skeleton precursors, potentially reflecting an ancestral diversity of skeletal tissue types. Together, these observations suggest that the ancestral pre-duplication lectican had a complex expression pattern, functioned to support mesenchymal histology, and likely played a role in the evolution of vertebrate-specific cell and tissue types.
Collapse
|
18
|
Abstract
Multiple sequence alignment is a core first step in many bioinformatics analyses, and errors in these alignments can have negative consequences for scientific studies. In this article, we review some of the recent literature evaluating multiple sequence alignment methods and identify specific challenges that arise when performing these evaluations. In particular, we discuss the different trends observed in simulation studies and when using biological benchmarks. Overall, we find that multiple sequence alignment, far from being a "solved problem," would benefit from new attention.
Collapse
Affiliation(s)
- Tandy Warnow
- Department of Computer Science, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
19
|
Fernandez‐Pozo N, Metz T, Chandler JO, Gramzow L, Mérai Z, Maumus F, Mittelsten Scheid O, Theißen G, Schranz ME, Leubner‐Metzger G, Rensing SA. Aethionema arabicum genome annotation using PacBio full-length transcripts provides a valuable resource for seed dormancy and Brassicaceae evolution research. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2021; 106:275-293. [PMID: 33453123 PMCID: PMC8641386 DOI: 10.1111/tpj.15161] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Revised: 12/31/2020] [Accepted: 01/08/2021] [Indexed: 05/06/2023]
Abstract
Aethionema arabicum is an important model plant for Brassicaceae trait evolution, particularly of seed (development, regulation, germination, dormancy) and fruit (development, dehiscence mechanisms) characters. Its genome assembly was recently improved but the gene annotation was not updated. Here, we improved the Ae. arabicum gene annotation using 294 RNA-seq libraries and 136 307 full-length PacBio Iso-seq transcripts, increasing BUSCO completeness by 11.6% and featuring 5606 additional genes. Analysis of orthologs showed a lower number of genes in Ae. arabicum than in other Brassicaceae, which could be partially explained by loss of homeologs derived from the At-α polyploidization event and by a lower occurrence of tandem duplications after divergence of Aethionema from the other Brassicaceae. Benchmarking of MADS-box genes identified orthologs of FUL and AGL79 not found in previous versions. Analysis of full-length transcripts related to ABA-mediated seed dormancy discovered a conserved isoform of PIF6-β and antisense transcripts in ABI3, ABI4 and DOG1, among other cases found of different alternative splicing between Turkey and Cyprus ecotypes. The presented data allow alternative splicing mining and proposition of numerous hypotheses to research evolution and functional genomics. Annotation data and sequences are available at the Ae. arabicum DB (https://plantcode.online.uni-marburg.de/aetar_db).
Collapse
Affiliation(s)
- Noe Fernandez‐Pozo
- Plant Cell BiologyDepartment of BiologyUniversity of MarburgMarburgGermany
| | - Timo Metz
- Plant Cell BiologyDepartment of BiologyUniversity of MarburgMarburgGermany
| | - Jake O. Chandler
- School of Biological SciencesRoyal Holloway University of LondonEghamSurreyUK
| | - Lydia Gramzow
- Matthias Schleiden Institute/GeneticsFriedrich Schiller University JenaJenaGermany
| | - Zsuzsanna Mérai
- Gregor Mendel Institute of Molecular Plant BiologyAustrian Academy of SciencesVienna BioCenter (VBC)ViennaAustria
| | | | - Ortrun Mittelsten Scheid
- Gregor Mendel Institute of Molecular Plant BiologyAustrian Academy of SciencesVienna BioCenter (VBC)ViennaAustria
| | - Günter Theißen
- Matthias Schleiden Institute/GeneticsFriedrich Schiller University JenaJenaGermany
| | - M. Eric Schranz
- Biosystematics GroupWageningen UniversityWageningenThe Netherlands
| | - Gerhard Leubner‐Metzger
- School of Biological SciencesRoyal Holloway University of LondonEghamSurreyUK
- Laboratory of Growth RegulatorsCentre of the Region Haná for Biotechnological and Agricultural ResearchPalacký University and Institute of Experimental BotanyAcademy of Sciences of the Czech RepublicOlomoucCzech Republic
| | - Stefan A. Rensing
- Plant Cell BiologyDepartment of BiologyUniversity of MarburgMarburgGermany
- BIOSS Centre for Biological Signaling StudiesUniversity of FreiburgFreiburgGermany
- LOEWE Center for Synthetic Microbiology (SYNMIKRO)University of MarburgMarburgGermany
| |
Collapse
|
20
|
Abstract
Multiple sequence alignment (MSA) is a central step in many bioinformatics and computational biology analyses. Although there exist many methods to perform MSA, most of them fail when dealing with large datasets due to their high computational cost. MSAProbs-MPI is a publicly available tool ( http://msaprobs.sourceforge.net ) that provides highly accurate results in relatively short runtime thanks to exploiting the hardware resources of multicore clusters. In this chapter, I explain the statistical and biological concepts employed in MSAProbs-MPI to complete the alignments, as well as the high-performance computing techniques used to accelerate it. Moreover, I provide some hints about the configuration parameters that should be used to guarantee high-performance executions.
Collapse
|
21
|
Proposal for a subdivision of the family Psathyrellaceae based on a taxon-rich phylogenetic analysis with iterative multigene guide tree. Mycol Prog 2020. [DOI: 10.1007/s11557-020-01606-3] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
22
|
Zhan Q, Fu Y, Jiang Q, Liu B, Peng J, Wang Y. SpliVert: A Protein Multiple Sequence Alignment Refinement Method Based on Splitting-Splicing Vertically. Protein Pept Lett 2020; 27:295-302. [PMID: 31385760 DOI: 10.2174/0929866526666190806143959] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2019] [Revised: 04/26/2019] [Accepted: 06/14/2019] [Indexed: 11/22/2022]
Abstract
BACKGROUND Multiple Sequence Alignment (MSA) is a fundamental task in bioinformatics and is required for many biological analysis tasks. The more accurate the alignments are, the more credible the downstream analyses. Most protein MSA algorithms realign an alignment to refine it by dividing it into two groups horizontally and then realign the two groups. However, this strategy does not consider that different regions of the sequences have different conservation; this property may lead to incorrect residue-residue or residue-gap pairs, which cannot be corrected by this strategy. OBJECTIVE In this article, our motivation is to develop a novel refinement method based on splitting- splicing vertically. METHODS Here, we present a novel refinement method based on splitting-splicing vertically, called SpliVert. For an alignment, we split it vertically into 3 parts, remove the gap characters in the middle, realign the middle part alone, and splice the realigned middle parts with the other two initial pieces to obtain a refined alignment. In the realign procedure of our method, the aligner will only focus on a certain part, ignoring the disturbance of the other parts, which could help fix the incorrect pairs. RESULTS We tested our refinement strategy for 2 leading MSA tools on 3 standard benchmarks, according to the commonly used average SP (and TC) score. The results show that given appropriate proportions to split the initial alignment, the average scores are increased comparably or slightly after using our method. We also compared the alignments refined by our method with alignments directly refined by the original alignment tools. The results suggest that using our SpliVert method to refine alignments can also outperform direct use of the original alignment tools. CONCLUSION The results reveal that splitting vertically and realigning part of the alignment is a good strategy for the refinement of protein multiple sequence alignments.
Collapse
Affiliation(s)
- Qing Zhan
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China
| | - Yilei Fu
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China
| | - Qinghua Jiang
- School of Life Science and Technology, Harbin Institute of Technology, Harbin, China
| | - Bo Liu
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China
| | - Jiajie Peng
- School of Computer Science, Northwestern Polytechnical University, Xi'an, China
| | - Yadong Wang
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China
| |
Collapse
|
23
|
A bi-objective function optimization approach for multiple sequence alignment using genetic algorithm. Soft comput 2020. [DOI: 10.1007/s00500-020-04917-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
|
24
|
Doll NM, Royek S, Fujita S, Okuda S, Chamot S, Stintzi A, Widiez T, Hothorn M, Schaller A, Geldner N, Ingram G. A two-way molecular dialogue between embryo and endosperm is required for seed development. Science 2020; 367:431-435. [PMID: 31974252 DOI: 10.1126/science.aaz4131] [Citation(s) in RCA: 69] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2019] [Accepted: 12/18/2019] [Indexed: 12/17/2023]
Abstract
The plant embryonic cuticle is a hydrophobic barrier deposited de novo by the embryo during seed development. At germination, it protects the seedling from water loss and is, thus, critical for survival. Embryonic cuticle formation is controlled by a signaling pathway involving the ABNORMAL LEAF SHAPE1 subtilase and the two GASSHO receptor-like kinases. We show that a sulfated peptide, TWISTED SEED1 (TWS1), acts as a GASSHO ligand. Cuticle surveillance depends on the action of the subtilase, which, unlike the TWS1 precursor and the GASSHO receptors, is not produced in the embryo but in the neighboring endosperm. Subtilase-mediated processing of the embryo-derived TWS1 precursor releases the active peptide, triggering GASSHO-dependent cuticle reinforcement in the embryo. Thus, a bidirectional molecular dialogue between embryo and endosperm safeguards cuticle integrity before germination.
Collapse
Affiliation(s)
- N M Doll
- Laboratoire Reproduction et Développement des Plantes, University of Lyon, ENS de Lyon, UCB Lyon 1, CNRS, INRAE, F-69342, Lyon, France
| | - S Royek
- Department of Plant Physiology and Biochemistry, University of Hohenheim, 70599 Stuttgart, Germany
| | - S Fujita
- Department of Plant Molecular Biology, University of Lausanne, 1015 Lausanne, Switzerland
| | - S Okuda
- Structural Plant Biology Laboratory, Department of Botany and Plant Biology, University of Geneva, 1211 Geneva, Switzerland
| | - S Chamot
- Laboratoire Reproduction et Développement des Plantes, University of Lyon, ENS de Lyon, UCB Lyon 1, CNRS, INRAE, F-69342, Lyon, France
| | - A Stintzi
- Department of Plant Physiology and Biochemistry, University of Hohenheim, 70599 Stuttgart, Germany
| | - T Widiez
- Laboratoire Reproduction et Développement des Plantes, University of Lyon, ENS de Lyon, UCB Lyon 1, CNRS, INRAE, F-69342, Lyon, France
| | - M Hothorn
- Structural Plant Biology Laboratory, Department of Botany and Plant Biology, University of Geneva, 1211 Geneva, Switzerland
| | - A Schaller
- Department of Plant Physiology and Biochemistry, University of Hohenheim, 70599 Stuttgart, Germany
| | - N Geldner
- Department of Plant Molecular Biology, University of Lausanne, 1015 Lausanne, Switzerland
| | - G Ingram
- Laboratoire Reproduction et Développement des Plantes, University of Lyon, ENS de Lyon, UCB Lyon 1, CNRS, INRAE, F-69342, Lyon, France.
| |
Collapse
|
25
|
Gramzow L, Lobbes D, Innard N, Theißen G. Independent origin of MIRNA genes controlling homologous target genes by partial inverted duplication of antisense-transcribed sequences. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2020; 101:401-419. [PMID: 31571291 DOI: 10.1111/tpj.14550] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2019] [Revised: 08/21/2019] [Accepted: 09/18/2019] [Indexed: 06/10/2023]
Abstract
Some microRNAs (miRNAs) are key regulators of developmental processes, mainly by controlling the accumulation of transcripts encoding transcription factors that are important for morphogenesis. MADS-box genes encode a family of transcription factors which control diverse developmental processes in flowering plants. Here we study the convergent evolution of two MIRNA (MIR) gene families, named MIR444 and MIR824, targeting members of the same clade of MIKCC -group MADS-box genes. We show that these two MIR genes most likely originated independently in monocots (MIR444) and in Brassicales (eudicots, MIR824). We provide evidence that, in both cases, the future target gene was transcribed in antisense prior to the evolution of the MIR genes. Both MIR genes then likely originated by a partial inverted duplication of their target genes, resulting in natural antisense organization of the newly evolved MIR gene and its target gene at birth. We thus propose a model for the origin of MIR genes, MEPIDAS (MicroRNA Evolution by Partial Inverted Duplication of Antisense-transcribed Sequences). MEPIDAS is a refinement of the inverted duplication hypothesis. According to MEPIDAS, a MIR gene evolves at a genomic locus at which the future target gene is also transcribed in the antisense direction. A partial inverted duplication at this locus causes the antisense transcript to fold into a stem-loop structure that is recognized by the miRNA biogenesis machinery to produce a miRNA that regulates the gene at this locus. Our analyses exemplify how to elucidate the origin of conserved miRNAs by comparative genomics and will guide future studies. OPEN RESEARCH BADGE: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/genbank/.
Collapse
Affiliation(s)
- Lydia Gramzow
- Department of Genetics, Friedrich Schiller University Jena, Philosophenweg 12, 07743, Jena, Germany
| | - Dajana Lobbes
- Department of Genetics, Friedrich Schiller University Jena, Philosophenweg 12, 07743, Jena, Germany
| | - Nathan Innard
- Department of Genetics, Friedrich Schiller University Jena, Philosophenweg 12, 07743, Jena, Germany
| | - Günter Theißen
- Department of Genetics, Friedrich Schiller University Jena, Philosophenweg 12, 07743, Jena, Germany
| |
Collapse
|
26
|
Maldonado E, Antunes A. LMAP_S: Lightweight Multigene Alignment and Phylogeny eStimation. BMC Bioinformatics 2019; 20:739. [PMID: 31888452 PMCID: PMC6937843 DOI: 10.1186/s12859-019-3292-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2019] [Accepted: 11/26/2019] [Indexed: 01/22/2023] Open
Abstract
Background Recent advances in genome sequencing technologies and the cost drop in high-throughput sequencing continue to give rise to a deluge of data available for downstream analyses. Among others, evolutionary biologists often make use of genomic data to uncover phenotypic diversity and adaptive evolution in protein-coding genes. Therefore, multiple sequence alignments (MSA) and phylogenetic trees (PT) need to be estimated with optimal results. However, the preparation of an initial dataset of multiple sequence file(s) (MSF) and the steps involved can be challenging when considering extensive amount of data. Thus, it becomes necessary the development of a tool that removes the potential source of error and automates the time-consuming steps of a typical workflow with high-throughput and optimal MSA and PT estimations. Results We introduce LMAP_S (Lightweight Multigene Alignment and Phylogeny eStimation), a user-friendly command-line and interactive package, designed to handle an improved alignment and phylogeny estimation workflow: MSF preparation, MSA estimation, outlier detection, refinement, consensus, phylogeny estimation, comparison and editing, among which file and directory organization, execution, manipulation of information are automated, with minimal manual user intervention. LMAP_S was developed for the workstation multi-core environment and provides a unique advantage for processing multiple datasets. Our software, proved to be efficient throughout the workflow, including, the (unlimited) handling of more than 20 datasets. Conclusions We have developed a simple and versatile LMAP_S package enabling researchers to effectively estimate multiple datasets MSAs and PTs in a high-throughput fashion. LMAP_S integrates more than 25 software providing overall more than 65 algorithm choices distributed in five stages. At minimum, one FASTA file is required within a single input directory. To our knowledge, no other software combines MSA and phylogeny estimation with as many alternatives and provides means to find optimal MSAs and phylogenies. Moreover, we used a case study comparing methodologies that highlighted the usefulness of our software. LMAP_S has been developed as an open-source package, allowing its integration into more complex open-source bioinformatics pipelines. LMAP_S package is released under GPLv3 license and is freely available at https://lmap-s.sourceforge.io/.
Collapse
Affiliation(s)
- Emanuel Maldonado
- CIIMAR/CIMAR - Interdisciplinary Centre of Marine and Environmental Research, University of Porto, Terminal de Cruzeiros do Porto de Leixões, Av. General Norton de Matos, s/n, 4450-208, Porto, Portugal
| | - Agostinho Antunes
- CIIMAR/CIMAR - Interdisciplinary Centre of Marine and Environmental Research, University of Porto, Terminal de Cruzeiros do Porto de Leixões, Av. General Norton de Matos, s/n, 4450-208, Porto, Portugal. .,Department of Biology, Faculty of Sciences, University of Porto, Rua do Campo Alegre, 4169-007, Porto, Portugal.
| |
Collapse
|
27
|
Zhan Q, Wang N, Jin S, Tan R, Jiang Q, Wang Y. ProbPFP: a multiple sequence alignment algorithm combining hidden Markov model optimized by particle swarm optimization with partition function. BMC Bioinformatics 2019; 20:573. [PMID: 31760933 PMCID: PMC6876095 DOI: 10.1186/s12859-019-3132-7] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND During procedures for conducting multiple sequence alignment, that is so essential to use the substitution score of pairwise alignment. To compute adaptive scores for alignment, researchers usually use Hidden Markov Model or probabilistic consistency methods such as partition function. Recent studies show that optimizing the parameters for hidden Markov model, as well as integrating hidden Markov model with partition function can raise the accuracy of alignment. The combination of partition function and optimized HMM, which could further improve the alignment's accuracy, however, was ignored by these researches. RESULTS A novel algorithm for MSA called ProbPFP is presented in this paper. It intergrate optimized HMM by particle swarm with partition function. The algorithm of PSO was applied to optimize HMM's parameters. After that, the posterior probability obtained by the HMM was combined with the one obtained by partition function, and thus to calculate an integrated substitution score for alignment. In order to evaluate the effectiveness of ProbPFP, we compared it with 13 outstanding or classic MSA methods. The results demonstrate that the alignments obtained by ProbPFP got the maximum mean TC scores and mean SP scores on these two benchmark datasets: SABmark and OXBench, and it got the second highest mean TC scores and mean SP scores on the benchmark dataset BAliBASE. ProbPFP is also compared with 4 other outstanding methods, by reconstructing the phylogenetic trees for six protein families extracted from the database TreeFam, based on the alignments obtained by these 5 methods. The result indicates that the reference trees are closer to the phylogenetic trees reconstructed from the alignments obtained by ProbPFP than the other methods. CONCLUSIONS We propose a new multiple sequence alignment method combining optimized HMM and partition function in this paper. The performance validates this method could make a great improvement of the alignment's accuracy.
Collapse
Affiliation(s)
- Qing Zhan
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, 150001, China
| | - Nan Wang
- Department of Mathematics, Harbin Institute of Technology, Harbin, 150001, China
| | - Shuilin Jin
- Department of Mathematics, Harbin Institute of Technology, Harbin, 150001, China
| | - Renjie Tan
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, 150001, China
| | - Qinghua Jiang
- School of Life Science and Technology, Harbin Institute of Technology, Harbin, 150001, China
| | - Yadong Wang
- School of Computer Science and Technology, Harbin Institute of Technology, Harbin, 150001, China.
| |
Collapse
|
28
|
Hoffmeier A, Gramzow L, Bhide AS, Kottenhagen N, Greifenstein A, Schubert O, Mummenhoff K, Becker A, Theißen G. A Dead Gene Walking: Convergent Degeneration of a Clade of MADS-Box Genes in Crucifers. Mol Biol Evol 2019; 35:2618-2638. [PMID: 30053121 DOI: 10.1093/molbev/msy142] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Genes are "born," and eventually they "die." These processes shape the phenotypic evolution of organisms and are hence of great biological interest. If genes die in plants, they generally do so quite rapidly. Here, we describe the fate of GOA-like genes that evolve in a dramatically different manner. GOA-like genes belong to the subfamily of Bsister genes of MIKC-type MADS-box genes. Typical MIKC-type genes encode conserved transcription factors controlling plant development. We show that ABS-like genes, a clade of Bsister genes, are indeed highly conserved in crucifers (Brassicaceae) maintaining the ancestral function of Bsister genes in ovule and seed development. In contrast, their closest paralogs, the GOA-like genes, have been undergoing convergent gene death in Brassicaceae. Intriguingly, erosion of GOA-like genes occurred after millions of years of coexistence with ABS-like genes. We thus describe Delayed Convergent Asymmetric Degeneration, a so far neglected but possibly frequent pattern of duplicate gene evolution that does not fit classical scenarios. Delayed Convergent Asymmetric Degeneration of GOA-like genes may have been initiated by a reduction in the expression of an ancestral GOA-like gene in the stem group of Brassicaceae and driven by dosage subfunctionalization. Our findings have profound implications for gene annotations in genomics, interpreting patterns of gene evolution and using genes in phylogeny reconstructions of species.
Collapse
Affiliation(s)
- Andrea Hoffmeier
- Genetics, Matthias Schleiden Institute, Friedrich-Schiller-University Jena, Jena, Germany
| | - Lydia Gramzow
- Genetics, Matthias Schleiden Institute, Friedrich-Schiller-University Jena, Jena, Germany
| | - Amey S Bhide
- Plant Developmental Biology Group, Institute of Botany, Justus-Liebig-University Giessen, Giessen, Germany
| | - Nina Kottenhagen
- Genetics, Matthias Schleiden Institute, Friedrich-Schiller-University Jena, Jena, Germany
| | - Andreas Greifenstein
- Genetics, Matthias Schleiden Institute, Friedrich-Schiller-University Jena, Jena, Germany
| | - Olesia Schubert
- Plant Developmental Biology Group, Institute of Botany, Justus-Liebig-University Giessen, Giessen, Germany
| | - Klaus Mummenhoff
- Department of Biology/Botany, University of Osnabrück, Osnabrück, Germany
| | - Annette Becker
- Plant Developmental Biology Group, Institute of Botany, Justus-Liebig-University Giessen, Giessen, Germany
| | - Günter Theißen
- Genetics, Matthias Schleiden Institute, Friedrich-Schiller-University Jena, Jena, Germany
| |
Collapse
|
29
|
Nute M, Saleh E, Warnow T. Evaluating Statistical Multiple Sequence Alignment in Comparison to Other Alignment Methods on Protein Data Sets. Syst Biol 2019; 68:396-411. [PMID: 30329135 PMCID: PMC6472439 DOI: 10.1093/sysbio/syy068] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 09/27/2018] [Accepted: 10/11/2018] [Indexed: 01/15/2023] Open
Abstract
The estimation of multiple sequence alignments of protein sequences is a basic step in many bioinformatics pipelines, including protein structure prediction, protein family identification, and phylogeny estimation. Statistical coestimation of alignments and trees under stochastic models of sequence evolution has long been considered the most rigorous technique for estimating alignments and trees, but little is known about the accuracy of such methods on biological benchmarks. We report the results of an extensive study evaluating the most popular protein alignment methods as well as the statistical coestimation method BAli-Phy on 1192 protein data sets from established benchmarks as well as on 120 simulated data sets. Our study (which used more than 230 CPU years for the BAli-Phy analyses alone) shows that BAli-Phy has better precision and recall (with respect to the true alignments) than the other alignment methods on the simulated data sets but has consistently lower recall on the biological benchmarks (with respect to the reference alignments) than many of the other methods. In other words, we find that BAli-Phy systematically underaligns when operating on biological sequence data but shows no sign of this on simulated data. There are several potential causes for this change in performance, including model misspecification, errors in the reference alignments, and conflicts between structural alignment and evolutionary alignments, and future research is needed to determine the most likely explanation. We conclude with a discussion of the potential ramifications for each of these possibilities. [BAli-Phy; homology; multiple sequence alignment; protein sequences; structural alignment.]
Collapse
Affiliation(s)
- Michael Nute
- Department of Statistics, University of Illinois at Urbana-Champaign, 725 S Wright St #101, Champaign, IL 61820, USA
| | - Ehsan Saleh
- Department of Computer Science, University of Illinois at Urbana-Champaign, 201 N. Goodwin Ave, Urbana, IL 61801, USA
| | - Tandy Warnow
- Department of Computer Science, University of Illinois at Urbana-Champaign, 201 N. Goodwin Ave, Urbana, IL 61801, USA.,Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign, 1205 W. Clark St., Urbana, IL 61801, USA.,National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
30
|
Velsko IM, Perez MS, Richards VP. Resolving Phylogenetic Relationships for Streptococcus mitis and Streptococcus oralis through Core- and Pan-Genome Analyses. Genome Biol Evol 2019; 11:1077-1087. [PMID: 30847473 PMCID: PMC6461889 DOI: 10.1093/gbe/evz049] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/05/2019] [Indexed: 02/06/2023] Open
Abstract
Taxonomic and phylogenetic relationships of Streptococcus mitis and Streptococcus oralis have been difficult to establish biochemically and genetically. We used core-genome analyses of S. mitis and S. oralis, as well as the closely related species Streptococcus pneumoniae and Streptococcus parasanguinis, to clarify the phylogenetic relationships between S. mitis and S. oralis, as well as within subclades of S. oralis. All S. mitis (n = 67), S. oralis (n = 89), S. parasanguinis (n = 27), and 27 S. pneumoniae genome assemblies were downloaded from NCBI and reannotated. All genes were delineated into homologous clusters and maximum-likelihood phylogenies built from putatively nonrecombinant core gene sets. Population structure was determined using Bayesian genome clustering, and patristic distance was calculated between populations. Population-specific gene content was assessed using a phylogenetic-based genome-wide association approach. Streptococcus mitis and S. oralis formed distinct clades, but species mixing suggests taxonomic misassignment. Patristic distance between populations suggests that S. oralis subsp. dentisani is a distinct species, whereas S. oralis subsp. tigurinus and subsp. oralis are supported as subspecies, and that S. mitis comprises two subspecies. None of the genes within the pan-genomes of S. mitis and S. oralis could be statistically correlated with either, and the dispensable genomes showed extensive variation among isolates. These are likely important factors contributing to established overlap in biochemical characteristics for these taxa. Based on core-genome analysis, the substructure of S. oralis and S. mitis should be redefined, and species assignments within S. oralis and S. mitis should be made based on whole-genome analysis to be robust to misassignment.
Collapse
Affiliation(s)
| | - Megan S Perez
- Department of Biological Sciences, Clemson University
- Department of Arts and Sciences, LeTourneau University
| | | |
Collapse
|
31
|
Sng NJ, Kolaczkowski B, Ferl RJ, Paul AL. A member of the CONSTANS-Like protein family is a putative regulator of reactive oxygen species homeostasis and spaceflight physiological adaptation. AOB PLANTS 2019; 11:ply075. [PMID: 30705745 PMCID: PMC6348315 DOI: 10.1093/aobpla/ply075] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2018] [Accepted: 12/12/2018] [Indexed: 05/20/2023]
Abstract
A feature of the physiological adaptation to spaceflight in Arabidopsis thaliana (Arabidopsis) is the induction of reactive oxygen species (ROS)-associated gene expression. The patterns of ROS-associated gene expression vary among Arabidopsis ecotypes, and the role of ROS signalling in spaceflight acclimation is unknown. What could differences in ROS gene regulation between ecotypes on orbit reveal about physiological adaptation to novel environments? Analyses of ecotype-dependent responses to spaceflight resulted in the elucidation of a previously uncharacterized gene (OMG1) as being ROS-associated. The OMG1 5' flanking region is an active promoter in cells where ROS activity is commonly observed, such as in pollen tubes, root hairs, and in other tissues upon wounding. qRT-PCR analyses revealed that upon wounding on Earth, OMG1 is an apparent transcriptional regulator of MYB77 and GRX480, which are associated with the ROS pathway. Fluorescence-based ROS assays show that OMG1 affects ROS production. Phylogenetic analysis of OMG1 and closely related homologs suggests that OMG1 is a distant, unrecognized member of the CONSTANS-Like protein family, a member that arose via gene duplication early in the angiosperm lineage and subsequently lost its first DNA-binding B-box1 domain. These data illustrate that members of the rapidly evolving COL protein family play a role in regulating ROS pathway functions, and their differential regulation on orbit suggests a role for ROS signalling in spaceflight physiological adaptation.
Collapse
Affiliation(s)
- Natasha J Sng
- Plant Molecular and Cellular Biology, University of Florida, Gainesville, FL, USA
| | - Bryan Kolaczkowski
- Microbiology and Cell Science, University of Florida, Gainesville, FL, USA
| | - Robert J Ferl
- Plant Molecular and Cellular Biology, University of Florida, Gainesville, FL, USA
- Horticultural Science Department, University of Florida, Gainesville, FL, USA
- Interdisciplinary Center for Biotechnology Research (ICBR), University of Florida, Gainesville, FL, USA
| | - Anna-Lisa Paul
- Plant Molecular and Cellular Biology, University of Florida, Gainesville, FL, USA
- Horticultural Science Department, University of Florida, Gainesville, FL, USA
- Corresponding author’s e-mail address:
| |
Collapse
|
32
|
High-Throughput Reconstruction of Ancestral Protein Sequence, Structure, and Molecular Function. Methods Mol Biol 2019; 1851:135-170. [PMID: 30298396 DOI: 10.1007/978-1-4939-8736-8_8] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Ancestral protein sequence reconstruction is a powerful technique for explicitly testing hypotheses about the evolution of molecular function, allowing researchers to meticulously dissect how historical changes in protein sequence impacted functional repertoire by altering the protein's 3D structure. These techniques have provided concrete, experimentally validated insights into ancient evolutionary processes and help illuminate the complex relationship between protein sequence, structure, and function. Inferring the protein family phylogenies on which ancestral sequence reconstruction depends and reconstructing the sequences, themselves, are amenable to high-throughput computational analysis. However, determining the structures of ancestral-reconstructed proteins and characterizing their functions typically rely on time-consuming and expensive laboratory analyses, limiting most current studies to examining a relatively small number of specific hypotheses. For this reason, we have little detailed, unbiased information about how molecular function evolves across large protein family phylogenies. Here we describe a generalized protocol that integrates ancestral sequence reconstruction with structural homology modeling and structure-based molecular affinity prediction to characterize historical changes in protein function across families with thousands of individual sequences. We highlight key steps in the analysis protocol requiring particularly careful attention to avoid introducing potential errors as well as steps for which computationally efficient subroutines can be substituted for more intensive approaches, allowing researchers to scale the analysis up or down, depending on available resources and requirements for reproducibility and scientific rigor. In our view, this approach provides a compelling compliment to more laboratory-intensive procedures, generating important contextual information that can help guide detailed experiments.
Collapse
|
33
|
Chaabane L. A hybrid solver for protein multiple sequence alignment problem. J Bioinform Comput Biol 2018; 16:1850015. [DOI: 10.1142/s0219720018500154] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In this work, a novel hybrid model called PSOSA for solving multiple sequence alignment (MSA) problem is proposed. The developed approach is a combination between particle swarm optimization (PSO) algorithm and simulated annealing (SA) technique. In our PSOSA approach, PSO is exploited in global search, but it is easily trapped into local optimum and may lead to premature convergence. SA is incorporated as local improvement approach to overcome local optimum problem and intensify the search in local regions to improve solution quality. Numerical results on BAliBASE benchmark have shown the effectiveness of the proposed method and its ability to achieve good quality solutions when compared with those given by other existing methods.
Collapse
Affiliation(s)
- Lamiche Chaabane
- Department of Computer Science, Mohamed Boudiaf University, BP. 166 M’sila 28000, Algeria
| |
Collapse
|
34
|
Orlando G, Raimondi D, Khan T, Lenaerts T, Vranken WF. SVM-dependent pairwise HMM: an application to protein pairwise alignments. Bioinformatics 2018; 33:3902-3908. [PMID: 28666322 DOI: 10.1093/bioinformatics/btx391] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2016] [Accepted: 06/12/2017] [Indexed: 12/27/2022] Open
Abstract
Motivation Methods able to provide reliable protein alignments are crucial for many bioinformatics applications. In the last years many different algorithms have been developed and various kinds of information, from sequence conservation to secondary structure, have been used to improve the alignment performances. This is especially relevant for proteins with highly divergent sequences. However, recent works suggest that different features may have different importance in diverse protein classes and it would be an advantage to have more customizable approaches, capable to deal with different alignment definitions. Results Here we present Rigapollo, a highly flexible pairwise alignment method based on a pairwise HMM-SVM that can use any type of information to build alignments. Rigapollo lets the user decide the optimal features to align their protein class of interest. It outperforms current state of the art methods on two well-known benchmark datasets when aligning highly divergent sequences. Availability and implementation A Python implementation of the algorithm is available at http://ibsquare.be/rigapollo. Contact wim.vranken@vub.be. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Gabriele Orlando
- Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, La Plaine Campus, Triomflaan.,Structural Biology Brussels, Vrije Universiteit Brussel, Pleinlaan 2.,Structural Biology Research Center, VIB.,Structural Machine Learning Group, Université Libre de Bruxelles
| | - Daniele Raimondi
- Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, La Plaine Campus, Triomflaan.,Structural Biology Brussels, Vrije Universiteit Brussel, Pleinlaan 2.,Structural Biology Research Center, VIB.,Structural Machine Learning Group, Université Libre de Bruxelles
| | - Taushif Khan
- Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, La Plaine Campus, Triomflaan.,Structural Biology Brussels, Vrije Universiteit Brussel, Pleinlaan 2
| | - Tom Lenaerts
- Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, La Plaine Campus, Triomflaan.,Structural Machine Learning Group, Université Libre de Bruxelles.,Artificial Intelligence Lab, Vrije Universiteit Brussel, 1050 Brussels, Belgium
| | - Wim F Vranken
- Interuniversity Institute of Bioinformatics in Brussels, ULB-VUB, La Plaine Campus, Triomflaan.,Structural Biology Brussels, Vrije Universiteit Brussel, Pleinlaan 2.,Structural Biology Research Center, VIB
| |
Collapse
|
35
|
Rubio-Largo Á, Vanneschi L, Castelli M, Vega-Rodríguez MA. Multiobjective characteristic-based framework for very-large multiple sequence alignment. Appl Soft Comput 2018. [DOI: 10.1016/j.asoc.2017.06.022] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
36
|
Lladós J, Cores F, Guirado F. Scalable Consistency in T-Coffee Through Apache Spark and Cassandra Database. J Comput Biol 2018; 25:894-906. [PMID: 30004242 DOI: 10.1089/cmb.2018.0084] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Next-generation sequencing, also known as high-throughput sequencing, has increased the volume of genetic data processed by sequencers. In the bioinformatic scientific area, highly rated multiple sequence alignment tools, such as MAFFT, ProbCons, and T-Coffee (TC), use the probabilistic consistency as a prior step to the progressive alignment stage to improve the final accuracy. However, such methods are severely limited by the memory required to store the consistency information. Big data processing and persistence techniques are used to manage and store the huge amount of information that is generated. Although these techniques have significant advantages, few biological applications have adopted them. In this article, a novel approach named big data tree-based consistency objective function for alignment evaluation (BDT-Coffee) is presented. BDT-Coffee is based on the integration of consistency information through Cassandra database in TC, previously generated by the MapReduce processing paradigm, to enable large data sets to be processed with the aim of improving the performance and scalability of the original algorithm.
Collapse
Affiliation(s)
- Jordi Lladós
- INSPIRES Research Center, Universitat de Lleida , Lleida, Spain
| | - Fernando Cores
- INSPIRES Research Center, Universitat de Lleida , Lleida, Spain
| | | |
Collapse
|
37
|
Sun J, Chen K, Hao Z. Pairwise alignment for very long nucleic acid sequences. Biochem Biophys Res Commun 2018; 502:313-317. [PMID: 29800571 DOI: 10.1016/j.bbrc.2018.05.134] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2018] [Accepted: 05/18/2018] [Indexed: 02/03/2023]
Abstract
Sequence alignment is one of the fundamental problems in computational biology and has numerous applications. The Smith-Waterman algorithm generates optimal local alignment for pairwise alignment task and has become a standard algorithm in its field. However, the current version of the Smith-Waterman algorithm demands a significant amount of memory and is not suitable for alignment of very long sequences. On the hand, the recent DNA sequencing technologies have produced vast amounts of biological sequences. Some nucleic acid sequences are very long and cannot employ the Smith-Waterman algorithm. To this end, this study proposes the PAAVLS algorithm that follows the dynamic programming technique employed by the Smith-Waterman algorithm and largely reduces the demand of memory. The proposed PAAVLS algorithm can be employed for alignment of very long sequences, i.e., sequences contain more than 100,000,000 nucleotides, on a personal computer. Additionally, the running time of the proposed PAAVLS algorithm is comparable with the running time of the standard Smith-Waterman algorithm.
Collapse
Affiliation(s)
- Jie Sun
- School of Computer Science and Software Engineering, Tianjin Polytechnic University, Tianjin, 300387, China
| | - Ke Chen
- School of Computer Science and Software Engineering, Tianjin Polytechnic University, Tianjin, 300387, China.
| | - Zhixiang Hao
- School of Computer Science and Software Engineering, Tianjin Polytechnic University, Tianjin, 300387, China
| |
Collapse
|
38
|
Alves RJE, Minh BQ, Urich T, von Haeseler A, Schleper C. Unifying the global phylogeny and environmental distribution of ammonia-oxidising archaea based on amoA genes. Nat Commun 2018; 9:1517. [PMID: 29666365 PMCID: PMC5904100 DOI: 10.1038/s41467-018-03861-1] [Citation(s) in RCA: 157] [Impact Index Per Article: 26.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2017] [Accepted: 03/19/2018] [Indexed: 12/30/2022] Open
Abstract
Ammonia-oxidising archaea (AOA) are ubiquitous and abundant in nature and play a major role in nitrogen cycling. AOA have been studied intensively based on the amoA gene (encoding ammonia monooxygenase subunit A), making it the most sequenced functional marker gene. Here, based on extensive phylogenetic and meta-data analyses of 33,378 curated archaeal amoA sequences, we define a highly resolved taxonomy and uncover global environmental patterns that challenge many earlier generalisations. Particularly, we show: (i) the global frequency of AOA is extremely uneven, with few clades dominating AOA diversity in most ecosystems; (ii) characterised AOA do not represent most predominant clades in nature, including soils and oceans; (iii) the functional role of the most prevalent environmental AOA clade remains unclear; and (iv) AOA harbour molecular signatures that possibly reflect phenotypic traits. Our work synthesises information from a decade of research and provides the first integrative framework to study AOA in a global context.
Collapse
Affiliation(s)
- Ricardo J Eloy Alves
- Archaea Biology and Ecogenomics Division, Department of Ecogenomics and Systems Biology, University of Vienna, Althanstrasse 14, 1090, Vienna, Austria
| | - Bui Quang Minh
- Center for Integrative Bioinformatics Vienna, Max F. Perutz Laboratories, University of Vienna, Medical University of Vienna, Campus Vienna Biocenter 5, Dr. Bohr Gasse 9, 1030, Vienna, Austria
- Ecology and Evolution, Research School of Biology, Australian National University, 2601, Canberra, ACT, Australia
| | - Tim Urich
- Institute of Microbiology, Ernst-Moritz-Arndt University, Felix-Hausdorff-Strasse 8, 17487, Greifswald, Germany
| | - Arndt von Haeseler
- Center for Integrative Bioinformatics Vienna, Max F. Perutz Laboratories, University of Vienna, Medical University of Vienna, Campus Vienna Biocenter 5, Dr. Bohr Gasse 9, 1030, Vienna, Austria
| | - Christa Schleper
- Archaea Biology and Ecogenomics Division, Department of Ecogenomics and Systems Biology, University of Vienna, Althanstrasse 14, 1090, Vienna, Austria.
| |
Collapse
|
39
|
Rubio-Largo A, Vanneschi L, Castelli M, Vega-Rodriguez MA. A Characteristic-Based Framework for Multiple Sequence Aligners. IEEE TRANSACTIONS ON CYBERNETICS 2018; 48:41-51. [PMID: 27831898 DOI: 10.1109/tcyb.2016.2621129] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The multiple sequence alignment is a well-known bioinformatics problem that consists in the alignment of three or more biological sequences (protein or nucleic acid). In the literature, a number of tools have been proposed for dealing with this biological sequence alignment problem, such as progressive methods, consistency-based methods, or iterative methods; among others. These aligners often use a default parameter configuration for all the input sequences to align. However, the default configuration is not always the best choice, the alignment accuracy of the tool may be highly boosted if specific parameter configurations are used, depending on the biological characteristics of the input sequences. In this paper, we propose a characteristic-based framework for multiple sequence aligners. The idea of the framework is, given an input set of unaligned sequences, extract its characteristics and run the aligner with the best parameter configuration found for another set of unaligned sequences with similar characteristics. In order to test the framework, we have used the well-known multiple sequence comparison by log-expectation (MUSCLE) v3.8 aligner with different benchmarks, such as benchmark alignments database v3.0, protein reference alignment benchmark v4.0, and sequence alignment benchmark v1.65. The results shown that the alignment accuracy and conservation of MUSCLE might be greatly improved with the proposed framework, specially in those scenarios with a low percentage of identity. The characteristic-based framework for multiple sequence aligners is freely available for downloading at http://arco.unex.es/arl/fwk-msa/cbf-msa.zip.
Collapse
|
40
|
Smith MA, Seemann SE, Quek XC, Mattick JS. DotAligner: identification and clustering of RNA structure motifs. Genome Biol 2017; 18:244. [PMID: 29284541 PMCID: PMC5747123 DOI: 10.1186/s13059-017-1371-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2017] [Accepted: 12/05/2017] [Indexed: 01/01/2023] Open
Abstract
The diversity of processed transcripts in eukaryotic genomes poses a challenge for the classification of their biological functions. Sparse sequence conservation in non-coding sequences and the unreliable nature of RNA structure predictions further exacerbate this conundrum. Here, we describe a computational method, DotAligner, for the unsupervised discovery and classification of homologous RNA structure motifs from a set of sequences of interest. Our approach outperforms comparable algorithms at clustering known RNA structure families, both in speed and accuracy. It identifies clusters of known and novel structure motifs from ENCODE immunoprecipitation data for 44 RNA-binding proteins.
Collapse
Affiliation(s)
- Martin A Smith
- RNA Biology and Plasticity Group, Garvan Institute of Medical Research, 384 Victoria Street, Sydney, NSW 2010, Australia. .,St Vincent's Clinical School, Faculty of Medicine, UNSW Australia, Sydney, NSW 2010, Australia.
| | - Stefan E Seemann
- Center for non-coding RNA in Technology and Health (RTH), University of Copenhagen, Groennegaardsvej 3, Frederiksberg, 1870, Denmark.,Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, DK-1870, Frederiksberg, Denmark
| | - Xiu Cheng Quek
- RNA Biology and Plasticity Group, Garvan Institute of Medical Research, 384 Victoria Street, Sydney, NSW 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW Australia, Sydney, NSW 2010, Australia
| | - John S Mattick
- RNA Biology and Plasticity Group, Garvan Institute of Medical Research, 384 Victoria Street, Sydney, NSW 2010, Australia.,St Vincent's Clinical School, Faculty of Medicine, UNSW Australia, Sydney, NSW 2010, Australia
| |
Collapse
|
41
|
Ascunce MS, Huguet-Tapia JC, Ortiz-Urquiza A, Keyhani NO, Braun EL, Goss EM. Phylogenomic analysis supports multiple instances of polyphyly in the oomycete peronosporalean lineage. Mol Phylogenet Evol 2017. [DOI: 10.1016/j.ympev.2017.06.013] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
42
|
Bacterial Foraging Optimization -Genetic Algorithm for Multiple Sequence Alignment with Multi-Objectives. Sci Rep 2017; 7:8833. [PMID: 28821841 PMCID: PMC5562892 DOI: 10.1038/s41598-017-09499-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2017] [Accepted: 07/27/2017] [Indexed: 01/06/2023] Open
Abstract
This research work focus on the multiple sequence alignment, as developing an exact multiple sequence alignment for different protein sequences is a difficult computational task. In this research, a hybrid algorithm named Bacterial Foraging Optimization-Genetic Algorithm (BFO-GA) algorithm is aimed to improve the multi-objectives and carrying out measures of multiple sequence alignment. The proposed algorithm employs multi-objectives such as variable gap penalty minimization, maximization of similarity and non-gap percentage. The proposed BFO-GA algorithm is measured with various MSA methods such as T-Coffee, Clustal Omega, Muscle, K-Align, MAFFT, GA, ACO, ABC and PSO. The experiments were taken on four benchmark datasets such as BAliBASE 3.0, Prefab 4.0, SABmark 1.65 and Oxbench 1.3 databases and the outcomes prove that the proposed BFO-GA algorithm obtains better statistical significance results as compared with the other well-known methods. This research study also evaluates the practicability of the alignments of BFO-GA by applying the optimal sequence to predict the phylogenetic tree by using ClustalW2 Phylogeny tool and compare with the existing algorithms by using the Robinson-Foulds (RF) distance performance metric. Lastly, the statistical implication of the proposed algorithm is computed by using the Wilcoxon Matched-Pair Signed- Rank test and also it infers better results.
Collapse
|
43
|
Gudyś A, Deorowicz S. QuickProbs 2: Towards rapid construction of high-quality alignments of large protein families. Sci Rep 2017; 7:41553. [PMID: 28139687 PMCID: PMC5282490 DOI: 10.1038/srep41553] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2016] [Accepted: 12/21/2016] [Indexed: 01/05/2023] Open
Abstract
The ever-increasing size of sequence databases caused by the development of high throughput sequencing, poses to multiple alignment algorithms one of the greatest challenges yet. As we show, well-established techniques employed for increasing alignment quality, i.e., refinement and consistency, are ineffective when large protein families are investigated. We present QuickProbs 2, an algorithm for multiple sequence alignment. Based on probabilistic models, equipped with novel column-oriented refinement and selective consistency, it offers outstanding accuracy. When analysing hundreds of sequences, Quick-Probs 2 is noticeably better than ClustalΩ and MAFFT, the previous leaders for processing numerous protein families. In the case of smaller sets, for which consistency-based methods are the best performing, QuickProbs 2 is also superior to the competitors. Due to low computational requirements of selective consistency and utilization of massively parallel architectures, presented algorithm has similar execution times to ClustalΩ, and is orders of magnitude faster than full consistency approaches, like MSAProbs or PicXAA. All these make QuickProbs 2 an excellent tool for aligning families ranging from few, to hundreds of proteins.
Collapse
Affiliation(s)
- Adam Gudyś
- Institute of Informatics, Silesian University of Technology, Akademicka 16, 44-100 Gliwice, Poland
| | - Sebastian Deorowicz
- Institute of Informatics, Silesian University of Technology, Akademicka 16, 44-100 Gliwice, Poland
| |
Collapse
|
44
|
Marra NJ, Richards VP, Early A, Bogdanowicz SM, Pavinski Bitar PD, Stanhope MJ, Shivji MS. Comparative transcriptomics of elasmobranchs and teleosts highlight important processes in adaptive immunity and regional endothermy. BMC Genomics 2017; 18:87. [PMID: 28132643 PMCID: PMC5278576 DOI: 10.1186/s12864-016-3411-x] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2016] [Accepted: 12/12/2016] [Indexed: 02/08/2023] Open
Abstract
BACKGROUND Comparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy. RESULTS Across seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species' heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as "antigen processing and presentation of exogenous peptide antigen via MHC class II", and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism. CONCLUSIONS This collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.
Collapse
Affiliation(s)
- Nicholas J Marra
- Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY, 14853, USA.,Save Our Seas Shark Research Center and Guy Harvey Research Institute, Nova Southeastern University, 8000 North Ocean Drive, Dania Beach, FL, 33004, USA
| | - Vincent P Richards
- Department of Biological Sciences, Clemson University, Clemson, SC, 29634, USA
| | - Angela Early
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, 14853, USA
| | - Steve M Bogdanowicz
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, 14853, USA
| | - Paulina D Pavinski Bitar
- Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY, 14853, USA
| | - Michael J Stanhope
- Department of Population Medicine and Diagnostic Sciences, College of Veterinary Medicine, Cornell University, Ithaca, NY, 14853, USA.
| | - Mahmood S Shivji
- Save Our Seas Shark Research Center and Guy Harvey Research Institute, Nova Southeastern University, 8000 North Ocean Drive, Dania Beach, FL, 33004, USA.
| |
Collapse
|
45
|
Pugh C, Kolaczkowski O, Manny A, Korithoski B, Kolaczkowski B. Resurrecting ancestral structural dynamics of an antiviral immune receptor: adaptive binding pocket reorganization repeatedly shifts RNA preference. BMC Evol Biol 2016; 16:241. [PMID: 27825296 PMCID: PMC5101713 DOI: 10.1186/s12862-016-0818-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2016] [Accepted: 10/28/2016] [Indexed: 02/07/2023] Open
Abstract
Background Although resurrecting ancestral proteins is a powerful tool for understanding the molecular-functional evolution of gene families, nearly all studies have examined proteins functioning in relatively stable biological processes. The extent to which more dynamic systems obey the same ‘rules’ governing stable processes is unclear. Here we present the first detailed investigation of the functional evolution of the RIG-like receptors (RLRs), a family of innate immune receptors that detect viral RNA in the cytoplasm. Results Using kinetic binding assays and molecular dynamics simulations of ancestral proteins, we demonstrate how a small number of adaptive protein-coding changes repeatedly shifted the RNA preference of RLRs throughout animal evolution by reorganizing the shape and electrostatic distribution across the RNA binding pocket, altering the hydrogen bond network between the RLR and its RNA target. In contrast to observations of proteins involved in metabolism and development, we find that RLR-RNA preference ‘flip flopped’ between two functional states, and shifts in RNA preference were not always coupled to gene duplications or speciation events. We demonstrate at least one reversion of RLR-RNA preference from a derived to an ancestral function through a novel structural mechanism, indicating multiple structural implementations of similar functions. Conclusions Our results suggest a model in which frequent shifts in selection pressures imposed by an evolutionary arms race preclude the long-term functional optimization observed in stable biological systems. As a result, the evolutionary dynamics of immune receptors may be less constrained by structural epistasis and historical contingency. Electronic supplementary material The online version of this article (doi:10.1186/s12862-016-0818-6) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Charles Pugh
- Department of Microbiology & Cell Science and Institute for Food and Agricultural Sciences, University of Florida, Gainesville, USA
| | - Oralia Kolaczkowski
- Department of Microbiology & Cell Science and Institute for Food and Agricultural Sciences, University of Florida, Gainesville, USA
| | - Austin Manny
- Department of Microbiology & Cell Science and Institute for Food and Agricultural Sciences, University of Florida, Gainesville, USA
| | - Bryan Korithoski
- Department of Microbiology & Cell Science and Institute for Food and Agricultural Sciences, University of Florida, Gainesville, USA
| | - Bryan Kolaczkowski
- Department of Microbiology & Cell Science and Institute for Food and Agricultural Sciences, University of Florida, Gainesville, USA. .,Genetics Institute, University of Florida, Gainesville, USA.
| |
Collapse
|
46
|
Yadav RK, Banka H. IBBOMSA: An Improved Biogeography-based Approach for Multiple Sequence Alignment. Evol Bioinform Online 2016; 12:237-246. [PMID: 27812276 PMCID: PMC5084829 DOI: 10.4137/ebo.s40457] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2016] [Revised: 08/10/2016] [Accepted: 08/16/2016] [Indexed: 11/05/2022] Open
Abstract
In bioinformatics, multiple sequence alignment (MSA) is an NP-hard problem. Hence, nature-inspired techniques can better approximate the solution. In the current study, a novel biogeography-based optimization (NBBO) is proposed to solve an MSA problem. The biogeography-based optimization (BBO) is a new paradigm for optimization. But, there exists some deficiencies in solving complicated problems such as low population diversity and slow convergence rate. NBBO is an enhanced version of BBO, in which, a new migration operation is proposed to overcome the limitations of BBO. The new migration adopts more information from other habitats, maintains population diversity, and preserves exploitation ability. In the performance analysis, the proposed and existing techniques such as VDGA, MOMSA, and GAPAM are tested on publicly available benchmark datasets (ie, Bali base). It has been observed that the proposed method shows the superiority/competitiveness with the existing techniques.
Collapse
Affiliation(s)
- Rohit Kumar Yadav
- Research Scholar, Department of Computer Science and Engineering, Indian Institute of Technology (ISM) Dhanbad, Jharkhand, India
| | - Haider Banka
- Associate Professor, Department of Computer Science and Engineering, Indian Institute of Technology (ISM) Dhanbad, Jharkhand, India
| |
Collapse
|
47
|
Li N, Wang K, Williams HN, Sun J, Ding C, Leng X, Dong K. Analysis of gene gain and loss in the evolution of predatory bacteria. Gene 2016; 598:63-70. [PMID: 27825775 DOI: 10.1016/j.gene.2016.10.039] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2016] [Revised: 10/11/2016] [Accepted: 10/27/2016] [Indexed: 01/18/2023]
Abstract
Predatory bacteria are ubiquitously distributed in nature in including in aquatic environments, sewage, intestinal tracts of animals and humans, rhizophere and, soils. However, our understanding of their evolutionary history is limited. Results of recent studies have shown that acquiring novel genes is a major force driving bacterial evolution. Therefore, to gain a better understanding of the impact of gene gain and loss in the evolution of bacterial predators, this study employed comparative genomic approaches to identify core-set gene families and species-specific gene families, and model gene gain and loss events among 11 genomes that represented diverse lineages. In total, 1977 gene families were classified. Of these 509 (pattern 11111111111) were present all of the 11 species. Among the non-core set gene families, 52 were present only in saltwater bacteria predators and had no ortholog in the other genomes. Similarly 109 and 44 were present only in the genomes of Micavibrio spp. and Bdellovibrio spp., respectively. In this study, the gain loss mapping engine GLOOME was selected to analyze and estimate the expectations and probabilities of both gain and loss events in the predatory bacteria. In total, 354 gene families were involved in significant gene gain events, and 407 gene families were classified into gene loss events with high supported value. Moreover, 18 families from the core set gene family were identified as putative genes under positive selection. The results of this study suggest that acquisition of particular genes that encode functional proteins in metabolism and cellular processes and signaling, especially ABC systems, may help bacterial predators adapt to surrounding environmental changes and present different predation strategies for survival in their habitats.
Collapse
Affiliation(s)
- Nan Li
- College of Marine Science and Environmental Sciences, Tianjin University of Science and Technology, Tianjin 300457, China; School of the Environment, Florida A&M University, Tallahassee, FL, USA.
| | - Kai Wang
- Institute of Biochemistry and Cell Biology, Shanghai Institute of Biological Science, Chinese Academy of Science, Shanghai, China
| | - Henry N Williams
- School of the Environment, Florida A&M University, Tallahassee, FL, USA
| | - Jun Sun
- College of Marine Science and Environmental Sciences, Tianjin University of Science and Technology, Tianjin 300457, China
| | - Changling Ding
- College of Marine Science and Environmental Sciences, Tianjin University of Science and Technology, Tianjin 300457, China
| | - Xiaoyun Leng
- College of Marine Science and Environmental Sciences, Tianjin University of Science and Technology, Tianjin 300457, China
| | - Ke Dong
- Department of Biological Sciences, Seoul National University, Seoul 151-742, South Korea
| |
Collapse
|
48
|
Palmer A, Painter J, Hassler H, Richards VP, Bruce T, Morrison S, Brown E, Kozak-Muiznieks NA, Lucas C, McNealy TL. Legionella clemsonensis
sp. nov.: a green fluorescing Legionella
strain from a patient with pneumonia. Microbiol Immunol 2016; 60:694-701. [DOI: 10.1111/1348-0421.12439] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2016] [Revised: 08/23/2016] [Accepted: 09/05/2016] [Indexed: 11/29/2022]
Affiliation(s)
- Allison Palmer
- Department of Biological Sciences; Life Sciences Building
| | - Joseph Painter
- Department of Biological Sciences; Life Sciences Building
| | - Hayley Hassler
- Department of Biological Sciences; Life Sciences Building
| | | | - Terri Bruce
- Clemson Light Imaging Facility; Clemson University, 190 Collings Street; Clemson SC 29631
| | - Shatavia Morrison
- Centers for Disease Control, 1600 Clifton Road; Atlanta GA 30333 USA
| | - Ellen Brown
- Centers for Disease Control, 1600 Clifton Road; Atlanta GA 30333 USA
| | | | - Claressa Lucas
- Centers for Disease Control, 1600 Clifton Road; Atlanta GA 30333 USA
| | | |
Collapse
|
49
|
Rani RR, Ramyachitra D. Multiple sequence alignment using multi-objective based bacterial foraging optimization algorithm. Biosystems 2016; 150:177-189. [PMID: 27784624 DOI: 10.1016/j.biosystems.2016.10.005] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2016] [Revised: 10/18/2016] [Accepted: 10/18/2016] [Indexed: 10/20/2022]
Abstract
Multiple sequence alignment (MSA) is a widespread approach in computational biology and bioinformatics. MSA deals with how the sequences of nucleotides and amino acids are sequenced with possible alignment and minimum number of gaps between them, which directs to the functional, evolutionary and structural relationships among the sequences. Still the computation of MSA is a challenging task to provide an efficient accuracy and statistically significant results of alignments. In this work, the Bacterial Foraging Optimization Algorithm was employed to align the biological sequences which resulted in a non-dominated optimal solution. It employs Multi-objective, such as: Maximization of Similarity, Non-gap percentage, Conserved blocks and Minimization of gap penalty. BAliBASE 3.0 benchmark database was utilized to examine the proposed algorithm against other methods In this paper, two algorithms have been proposed: Hybrid Genetic Algorithm with Artificial Bee Colony (GA-ABC) and Bacterial Foraging Optimization Algorithm. It was found that Hybrid Genetic Algorithm with Artificial Bee Colony performed better than the existing optimization algorithms. But still the conserved blocks were not obtained using GA-ABC. Then BFO was used for the alignment and the conserved blocks were obtained. The proposed Multi-Objective Bacterial Foraging Optimization Algorithm (MO-BFO) was compared with widely used MSA methods Clustal Omega, Kalign, MUSCLE, MAFFT, Genetic Algorithm (GA), Ant Colony Optimization (ACO), Artificial Bee Colony (ABC), Particle Swarm Optimization (PSO) and Hybrid Genetic Algorithm with Artificial Bee Colony (GA-ABC). The final results show that the proposed MO-BFO algorithm yields better alignment than most widely used methods.
Collapse
Affiliation(s)
- R Ranjani Rani
- Department of Computer Science, Bharathiar University, Coimbatore, Tamilnadu, India.
| | - D Ramyachitra
- Department of Computer Science, Bharathiar University, Coimbatore, Tamilnadu, India.
| |
Collapse
|
50
|
Jeong H, Qian X, Yoon BJ. Effective comparative analysis of protein-protein interaction networks by measuring the steady-state network flow using a Markov model. BMC Bioinformatics 2016; 17:395. [PMID: 27766938 PMCID: PMC5073945 DOI: 10.1186/s12859-016-1215-2] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Background Comparative analysis of protein-protein interaction (PPI) networks provides an effective means of detecting conserved functional network modules across different species. Such modules typically consist of orthologous proteins with conserved interactions, which can be exploited to computationally predict the modules through network comparison. Results In this work, we propose a novel probabilistic framework for comparing PPI networks and effectively predicting the correspondence between proteins, represented as network nodes, that belong to conserved functional modules across the given PPI networks. The basic idea is to estimate the steady-state network flow between nodes that belong to different PPI networks based on a Markov random walk model. The random walker is designed to make random moves to adjacent nodes within a PPI network as well as cross-network moves between potential orthologous nodes with high sequence similarity. Based on this Markov random walk model, we estimate the steady-state network flow – or the long-term relative frequency of the transitions that the random walker makes – between nodes in different PPI networks, which can be used as a probabilistic score measuring their potential correspondence. Subsequently, the estimated scores can be used for detecting orthologous proteins in conserved functional modules through network alignment. Conclusions Through evaluations based on multiple real PPI networks, we demonstrate that the proposed scheme leads to improved alignment results that are biologically more meaningful at reduced computational cost, outperforming the current state-of-the-art algorithms. The source code and datasets can be downloaded from http://www.ece.tamu.edu/~bjyoon/CUFID.
Collapse
Affiliation(s)
- Hyundoo Jeong
- Department of Electrical and Computer Engineering, Texas A&M University, College Station, USA
| | - Xiaoning Qian
- Department of Electrical and Computer Engineering, Texas A&M University, College Station, USA
| | - Byung-Jun Yoon
- Department of Electrical and Computer Engineering, Texas A&M University, College Station, USA.
| |
Collapse
|