1
|
Pratomo A, Bengen DG, Zamani NP, Lane C, Humphries AT, Borbee E, Subhan B, Madduppa H. Diversity and distribution of Symbiodiniaceae detected on coral reefs of Lombok, Indonesia using environmental DNA metabarcoding. PeerJ 2022; 10:e14006. [PMID: 36312748 PMCID: PMC9610659 DOI: 10.7717/peerj.14006] [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: 10/06/2021] [Accepted: 08/14/2022] [Indexed: 01/19/2023] Open
Abstract
Background Dinoflagellates of family Symbiodiniaceae are important to coral reef ecosystems because of their contribution to coral health and growth; however, only a few studies have investigated the function and distribution of Symbiodiniaceae in Indonesia. Understanding the distribution of different kinds of Symbiodiniaceae can improve forecasting of future responses of various coral reef systems to climate change. This study aimed to determine the diversity of Symbiodiniaceae around Lombok using environmental DNA (eDNA). Methods Seawater and sediment samples were collected from 18 locations and filtered to obtain fractions of 0.4-12 and >12 µm. After extraction, molecular barcoding polymerase chain reaction was conducted to amplify the primary V9-SSU 18S rRNA gene, followed by sequencing (Illumina MiSeq). BLAST, Naïve-fit-Bayes, and maximum likelihood routines were used for classification and phylogenetic reconstruction. We compared results across sampling sites, sample types (seawater/sediment), and filter pore sizes (fraction). Results Phylogenetic analyses resolved the amplicon sequence variants into 16 subclades comprising six Symbiodiniaceae genera (or genera-equivalent clades) as follows: Symbiodinium, Breviolum, Cladocopium, Durusdinium, Foraminifera Clade G, and Halluxium. Comparative analyses showed that the three distinct lineages within Cladocopium, Durusdinium, and Foraminifera Clade G were the most common. Most of the recovered sequences appeared to be distinctive of different sampling locations, supporting the possibility that eDNA may resolve regional and local differences among Symbiodiniaceae genera and species. Conclusions eDNA surveys offer a rapid proxy for evaluating Symbiodiniaceae species on coral reefs and are a potentially useful approach to revealing diversity and relative ecological dominance of certain Symbiodiniaceae organisms. Moreover, Symbiodiniaceae eDNA analysis shows potential in monitoring the local and regional stability of coral-algal mutualisms.
Collapse
Affiliation(s)
- Arief Pratomo
- Raja Ali Haji Maritime University, Tanjungpinang, Indonesia,Department of Marine Science and Technology, Institut Pertanian Bogor, Bogor, Indonesia
| | - Dietriech G. Bengen
- Department of Marine Science and Technology, Institut Pertanian Bogor, Bogor, Indonesia
| | - Neviaty P. Zamani
- Department of Marine Science and Technology, Institut Pertanian Bogor, Bogor, Indonesia
| | - Christopher Lane
- Department of Biological Sciences, University of Rhode Island, Rhode Island, United States of America
| | - Austin T. Humphries
- Department of Fisheries, Animal and Veterinary Sciences, University of Rhode Island, Rhode Island, United States of America
| | - Erin Borbee
- Department of Biological Sciences, University of Rhode Island, Rhode Island, United States of America
| | - Beginer Subhan
- Department of Marine Science and Technology, Institut Pertanian Bogor, Bogor, Indonesia
| | - Hawis Madduppa
- Department of Marine Science and Technology, Institut Pertanian Bogor, Bogor, Indonesia,Oceanogen Research Center, Bogor, Indonesia
| |
Collapse
|
2
|
Zhang C, Zhao Y, Braun EL, Mirarab S. TAPER: Pinpointing errors in multiple sequence alignments despite varying rates of evolution. Methods Ecol Evol 2021. [DOI: 10.1111/2041-210x.13696] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Chao Zhang
- Bioinformatics and Systems Biology Program University of California San Diego CA USA
| | - Yiming Zhao
- Electrical and Computer Engineering Department University of California San Diego CA USA
| | - Edward L. Braun
- Department of Biology and Genetics Institute University of Florida Gainesville FL USA
| | - Siavash Mirarab
- Electrical and Computer Engineering Department University of California San Diego CA USA
| |
Collapse
|
3
|
Wang W, Hejasebazzi A, Zheng J, Liu KJ. Build a better bootstrap and the RAWR shall beat a random path to your door: phylogenetic support estimation revisited. Bioinformatics 2021; 37:i111-i119. [PMID: 34252944 PMCID: PMC8336443 DOI: 10.1093/bioinformatics/btab263] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
Motivation The standard bootstrap method is used throughout science and engineering to perform general-purpose non-parametric resampling and re-estimation. Among the most widely cited and widely used such applications is the phylogenetic bootstrap method, which Felsenstein proposed in 1985 as a means to place statistical confidence intervals on an estimated phylogeny (or estimate ‘phylogenetic support’). A key simplifying assumption of the bootstrap method is that input data are independent and identically distributed (i.i.d.). However, the i.i.d. assumption is an over-simplification for biomolecular sequence analysis, as Felsenstein noted. Results In this study, we introduce a new sequence-aware non-parametric resampling technique, which we refer to as RAWR (‘RAndom Walk Resampling’). RAWR consists of random walks that synthesize and extend the standard bootstrap method and the ‘mirrored inputs’ idea of Landan and Graur. We apply RAWR to the task of phylogenetic support estimation. RAWR’s performance is compared to the state-of-the-art using synthetic and empirical data that span a range of dataset sizes and evolutionary divergence. We show that RAWR support estimates offer comparable or typically superior type I and type II error compared to phylogenetic bootstrap support. We also conduct a re-analysis of large-scale genomic sequence data from a recent study of Darwin’s finches. Our findings clarify phylogenetic uncertainty in a charismatic clade that serves as an important model for complex adaptive evolution. Availability and implementation Data and software are publicly available under open-source software and open data licenses at: https://gitlab.msu.edu/liulab/RAWR-study-datasets-and-scripts.
Collapse
Affiliation(s)
- Wei Wang
- Department of Computer Science and Engineering, Michigan State University, East Lansing, MI 48824, USA
| | - Ahmad Hejasebazzi
- Department of Computer Science and Engineering, Michigan State University, East Lansing, MI 48824, USA
| | - Julia Zheng
- Department of Computer Science and Engineering, Michigan State University, East Lansing, MI 48824, USA.,Ecology, Evolution, and Behavior Program, Michigan State University, East Lansing, MI 48824, USA
| | - Kevin J Liu
- Department of Computer Science and Engineering, Michigan State University, East Lansing, MI 48824, USA.,Ecology, Evolution, and Behavior Program, Michigan State University, East Lansing, MI 48824, USA.,Genetics and Genome Sciences Program, Michigan State University, East Lansing, MI 48824, USA
| |
Collapse
|
4
|
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
|
5
|
Esparza M, Jedlicki E, González C, Dopson M, Holmes DS. Effect of CO 2 Concentration on Uptake and Assimilation of Inorganic Carbon in the Extreme Acidophile Acidithiobacillus ferrooxidans. Front Microbiol 2019; 10:603. [PMID: 31019493 PMCID: PMC6458275 DOI: 10.3389/fmicb.2019.00603] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Accepted: 03/11/2019] [Indexed: 02/01/2023] Open
Abstract
This study was motivated by surprising gaps in the current knowledge of microbial inorganic carbon (Ci) uptake and assimilation at acidic pH values (pH < 3). Particularly striking is the limited understanding of the differences between Ci uptake mechanisms in acidic versus circumneutral environments where the Ci predominantly occurs either as a dissolved gas (CO2) or as bicarbonate (HCO3 -), respectively. In order to gain initial traction on the problem, the relative abundance of transcripts encoding proteins involved in Ci uptake and assimilation was studied in the autotrophic, polyextreme acidophile Acidithiobacillus ferrooxidans whose optimum pH for growth is 2.5 using ferrous iron as an energy source, although they are able to grow at pH 5 when using sulfur as an energy source. The relative abundance of transcripts of five operons (cbb1-5) and one gene cluster (can-sulP) was monitored by RT-qPCR and, in selected cases, at the protein level by Western blotting, when cells were grown under different regimens of CO2 concentration in elemental sulfur. Of particular note was the absence of a classical bicarbonate uptake system in A. ferrooxidans. However, bioinformatic approaches predict that sulP, previously annotated as a sulfate transporter, is a novel type of bicarbonate transporter. A conceptual model of CO2 fixation was constructed from combined bioinformatic and experimental approaches that suggests strategies for providing ecological flexibility under changing concentrations of CO2 and provides a portal to elucidating Ci uptake and regulation in acidic conditions. The results could advance the understanding of industrial bioleaching processes to recover metals such as copper at acidic pH. In addition, they may also shed light on how chemolithoautotrophic acidophiles influence the nutrient and energy balance in naturally occurring low pH environments.
Collapse
Affiliation(s)
- Mario Esparza
- Laboratorio de Biominería, Departamento de Biotecnología, Facultad de Ciencias del Mar y Recursos Biológicos, Universidad de Antofagasta, Antofagasta, Chile
| | - Eugenia Jedlicki
- Center for Bioinformatics and Genome Biology, Fundación Ciencia & Vida, Santiago, Chile
| | - Carolina González
- Center for Bioinformatics and Genome Biology, Fundación Ciencia & Vida, Santiago, Chile
| | - Mark Dopson
- Centre for Ecology and Evolution in Microbial Model Systems (EEMiS), Linnaeus University, Kalmar, Sweden
| | - David S. Holmes
- Center for Bioinformatics and Genome Biology, Fundación Ciencia & Vida, Santiago, Chile
- Centro de Genómica y Bioinformática, Facultad de Ciencias, Universidad Mayor, Santiago, Chile
| |
Collapse
|
6
|
Ashkenazy H, Sela I, Levy Karin E, Landan G, Pupko T. Multiple Sequence Alignment Averaging Improves Phylogeny Reconstruction. Syst Biol 2018; 68:117-130. [PMID: 29771363 DOI: 10.1093/sysbio/syy036] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2017] [Accepted: 05/09/2018] [Indexed: 01/11/2023] Open
Abstract
The classic methodology of inferring a phylogenetic tree from sequence data is composed of two steps. First, a multiple sequence alignment (MSA) is computed. Then, a tree is reconstructed assuming the MSA is correct. Yet, inferred MSAs were shown to be inaccurate and alignment errors reduce tree inference accuracy. It was previously proposed that filtering unreliable alignment regions can increase the accuracy of tree inference. However, it was also demonstrated that the benefit of this filtering is often obscured by the resulting loss of phylogenetic signal. In this work we explore an approach, in which instead of relying on a single MSA, we generate a large set of alternative MSAs and concatenate them into a single SuperMSA. By doing so, we account for phylogenetic signals contained in columns that are not present in the single MSA computed by alignment algorithms. Using simulations, we demonstrate that this approach results, on average, in more accurate trees compared to 1) using an unfiltered MSA and 2) using a single MSA with weights assigned to columns according to their reliability. Next, we explore in which regions of the MSA space our approach is expected to be beneficial. Finally, we provide a simple criterion for deciding whether or not the extra effort of computing a SuperMSA and inferring a tree from it is beneficial. Based on these assessments, we expect our methodology to be useful for many cases in which diverged sequences are analyzed. The option to generate such a SuperMSA is available at http://guidance.tau.ac.il.
Collapse
Affiliation(s)
- Haim Ashkenazy
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel Aviv University, Ramat Aviv 69978, Tel Aviv, Israel
| | - Itamar Sela
- National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894, USA
| | - Eli Levy Karin
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel Aviv University, Ramat Aviv 69978, Tel Aviv, Israel.,Department of Molecular Biology & Ecology of Plants, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Giddy Landan
- Institute of Microbiology, Christian-Albrechts-University of Kiel, 24118 Kiel, Germany
| | - Tal Pupko
- Department of Cell Research and Immunology, George S. Wise Faculty of Life Sciences, Tel Aviv University, Ramat Aviv 69978, Tel Aviv, Israel
| |
Collapse
|
7
|
Gates DJ, Pilson D, Smith SD. Filtering of target sequence capture individuals facilitates species tree construction in the plant subtribe Iochrominae (Solanaceae). Mol Phylogenet Evol 2018; 123:26-34. [DOI: 10.1016/j.ympev.2018.02.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2017] [Revised: 01/30/2018] [Accepted: 02/01/2018] [Indexed: 10/18/2022]
|
8
|
Sun M, Soltis DE, Soltis PS, Zhu X, Burleigh JG, Chen Z. Deep phylogenetic incongruence in the angiosperm clade Rosidae. Mol Phylogenet Evol 2015; 83:156-66. [DOI: 10.1016/j.ympev.2014.11.003] [Citation(s) in RCA: 82] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2014] [Revised: 11/01/2014] [Accepted: 11/05/2014] [Indexed: 10/24/2022]
|
9
|
Drew BT, Ruhfel BR, Smith SA, Moore MJ, Briggs BG, Gitzendanner MA, Soltis PS, Soltis DE. Another Look at the Root of the Angiosperms Reveals a Familiar Tale. Syst Biol 2014; 63:368-82. [DOI: 10.1093/sysbio/syt108] [Citation(s) in RCA: 64] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
|