1
|
Prillo S, Ravoor A, Yosef N, Song YS. ConvexML: Scalable and accurate inference of single-cell chronograms from CRISPR/Cas9 lineage tracing data. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.12.03.569785. [PMID: 38076815 PMCID: PMC10705529 DOI: 10.1101/2023.12.03.569785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
CRISPR/Cas9 gene editing technology has enabled lineage tracing for thousands of cells in vivo. However, most of the analysis of CRISPR/Cas9 lineage tracing data has so far been limited to the reconstruction of single-cell tree topologies, which depict lineage relationships between cells, but not the amount of time that has passed between ancestral cell states and the present. Time-resolved trees, known as chronograms, would allow one to study the evolutionary dynamics of cell populations at an unprecedented level of resolution. Indeed, time-resolved trees would reveal the timing of events on the tree, the relative fitness of subclones, and the dynamics underlying phenotypic changes in the cell population - among other important applications. In this work, we introduce the first scalable and accurate method to refine any given single-cell tree topology into a single-cell chronogram by estimating its branch lengths. To do this, we leverage a statistical model of CRISPR/Cas9 cutting with missing data, paired with a conservative version of maximum parsimony that reconstructs only the ancestral states that we are confident about. As part of our method, we propose a novel approach to represent and handle missing data - specifically, double-resection events - which greatly simplifies and speeds up branch length estimation without compromising quality. All this leads to a convex maximum likelihood estimation (MLE) problem that can be readily solved in seconds with off-the-shelf convex optimization solvers. To stabilize estimates in low-information regimes, we propose a simple penalized version of MLE using a minimum branch length and pseudocounts. We benchmark our method using simulations and show that it performs well on several tasks, outperforming more naive baselines. Our method, which we name 'ConvexML', is available through the cassiopeia open source Python package.
Collapse
Affiliation(s)
| | - Akshay Ravoor
- Computer Science Division, University of California, Berkeley
| | - Nir Yosef
- Computer Science Division, University of California, Berkeley
- Department of Systems Immunology, Weizmann Institute of Science
| | - Yun S. Song
- Computer Science Division, University of California, Berkeley
- Department of Statistics, University of California, Berkeley
| |
Collapse
|
2
|
Lu B, Curtius K, Graham TA, Yang Z, Barnes CP. CNETML: maximum likelihood inference of phylogeny from copy number profiles of multiple samples. Genome Biol 2023; 24:144. [PMID: 37340508 DOI: 10.1186/s13059-023-02983-0] [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: 03/24/2022] [Accepted: 06/08/2023] [Indexed: 06/22/2023] Open
Abstract
Phylogenetic trees based on copy number profiles from multiple samples of a patient are helpful to understand cancer evolution. Here, we develop a new maximum likelihood method, CNETML, to infer phylogenies from such data. CNETML is the first program to jointly infer the tree topology, node ages, and mutation rates from total copy numbers of longitudinal samples. Our extensive simulations suggest CNETML performs well on copy numbers relative to ploidy and under slight violation of model assumptions. The application of CNETML to real data generates results consistent with previous discoveries and provides novel early copy number events for further investigation.
Collapse
Affiliation(s)
- Bingxin Lu
- Department of Cell and Developmental Biology, University College London, London, UK.
- UCL Genetics Institute, University College London, London, UK.
| | - Kit Curtius
- Barts Cancer Institute, Barts and the London School of Medicine and Dentistry, Queen Mary University of London, London, UK
- Division of Biomedical Informatics, Department of Medicine, University of California San Diego, La Jolla, CA, USA
| | - Trevor A Graham
- Barts Cancer Institute, Barts and the London School of Medicine and Dentistry, Queen Mary University of London, London, UK
- Centre for Evolution and Cancer, Institute of Cancer Research, London, UK
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London, UK
| | - Chris P Barnes
- Department of Cell and Developmental Biology, University College London, London, UK.
- UCL Genetics Institute, University College London, London, UK.
| |
Collapse
|
3
|
Lasley RM, Evans N, Paulay G, Michonneau F, Windsor A, Irwansyah, Ng PKL. Allopatric mosaics in the Indo-West Pacific crab subfamily Chlorodiellinae reveal correlated patterns of sympatry, genetic divergence, and genitalic disparity. Mol Phylogenet Evol 2023; 181:107710. [PMID: 36707010 DOI: 10.1016/j.ympev.2023.107710] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2022] [Revised: 12/31/2022] [Accepted: 01/18/2023] [Indexed: 01/26/2023]
Abstract
Molecular studies have revealed that many species once thought to be wide-ranging in the Indo-West Pacific contain allopatric mosaics of endemic lineages. These lineages provide compelling evidence that substantial time is needed to evolve isolating mechanisms sufficient to permit successful secondary sympatry, and that divergence is initiated in allopatry. In this context, questions arise regarding the nature, timing, and origin of isolating mechanisms that permit secondary sympatry. We present a phylogeny of the crab subfamily Chlorodiellinae which displays allopatric mosaics within species. These allopatric lineages typically do not have divergent male genitalia, while older sympatric lineages do. We tested the relationship between genetic distance (proxy for time), sympatry, and the divergence of male genitalic morphology. Our results suggest that male genitalic divergence is not involved in the initiation of speciation in chlorodielline crabs, having likely occurred only after isolation began in allopatry. However, morphological evolution of genitalia seemingly does play an important role in completing the process of speciation in these crabs.
Collapse
Affiliation(s)
- Robert M Lasley
- Florida Museum of Natural History, University of Florida, Gainesville, FL 32611-7800, USA; University of Guam, EPSCoR, Mangilao 96926, Guam; National Museum of Natural History, Smithsonian Institution, Suitland, MD 20746, USA.
| | - Nathaniel Evans
- Florida Museum of Natural History, University of Florida, Gainesville, FL 32611-7800, USA; National Museum of Natural History, Smithsonian Institution, Suitland, MD 20746, USA; Tropical Biosphere Research Center, University of the Ryukyus, Nishihara, Okinawa 903-0213, Japan
| | - Gustav Paulay
- Florida Museum of Natural History, University of Florida, Gainesville, FL 32611-7800, USA
| | - Francois Michonneau
- Florida Museum of Natural History, University of Florida, Gainesville, FL 32611-7800, USA
| | - Amanda Windsor
- National Museum of Natural History, Smithsonian Institution, Suitland, MD 20746, USA
| | - Irwansyah
- Syiah Kuala University, Aceh 23111, Indonesia
| | - Peter K L Ng
- Lee Kong Chian Natural History Museum, National University of Singapore, 117377, Singapore
| |
Collapse
|
4
|
Bergeron LA, Besenbacher S, Zheng J, Li P, Bertelsen MF, Quintard B, Hoffman JI, Li Z, St Leger J, Shao C, Stiller J, Gilbert MTP, Schierup MH, Zhang G. Evolution of the germline mutation rate across vertebrates. Nature 2023; 615:285-291. [PMID: 36859541 PMCID: PMC9995274 DOI: 10.1038/s41586-023-05752-y] [Citation(s) in RCA: 73] [Impact Index Per Article: 73.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Accepted: 01/23/2023] [Indexed: 03/03/2023]
Abstract
The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself1. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies2. Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis3. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.
Collapse
Affiliation(s)
- Lucie A Bergeron
- Villum Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Søren Besenbacher
- Department of Molecular Medicine, Aarhus University, Aarhus, Denmark
| | - Jiao Zheng
- BGI-Shenzhen, Shenzhen, China
- BGI Education Center, University of Chinese Academy of Sciences, Shenzhen, China
| | | | | | | | - Joseph I Hoffman
- Department of Animal Behaviour, Bielefeld University, Bielefeld, Germany
- British Antarctic Survey, High Cross, Cambridge, UK
| | - Zhipeng Li
- College of Animal Science and Technology, Jilin Agricultural University, Changchun, China
| | - Judy St Leger
- Department of Biomedical Sciences, Cornell University, Ithaca, NY, USA
| | - Changwei Shao
- Key Lab of Sustainable Development of Marine Fisheries, Ministry of Agriculture and Rural Affairs, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao, China
| | - Josefin Stiller
- Villum Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - M Thomas P Gilbert
- Center for Evolutionary Hologenomics, The GLOBE Institute, University of Copenhagen, Copenhagen, Denmark
- University Museum, NTNU, Trondheim, Norway
| | | | - Guojie Zhang
- Villum Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
- Centre for Evolutionary & Organismal Biology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China.
- Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, China.
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China.
| |
Collapse
|
5
|
Valderrama E, Landis JB, Skinner D, Maas PJM, Maas-van de Kramer H, André T, Grunder N, Sass C, Pinilla-Vargas M, Guan CJ, Phillips HR, de Almeida AMR, Specht CD. The genetic mechanisms underlying the convergent evolution of pollination syndromes in the Neotropical radiation of Costus L. FRONTIERS IN PLANT SCIENCE 2022; 13:874322. [PMID: 36161003 PMCID: PMC9493542 DOI: 10.3389/fpls.2022.874322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/12/2022] [Accepted: 06/27/2022] [Indexed: 06/16/2023]
Abstract
Selection together with variation in floral traits can act to mold floral form, often driven by a plant's predominant or most effective pollinators. To investigate the evolution of traits associated with pollination, we developed a phylogenetic framework for evaluating tempo and mode of pollination shifts across the genus Costus L., known for its evolutionary toggle between traits related to bee and bird pollination. Using a target enrichment approach, we obtained 957 loci for 171 accessions to expand the phylogenetic sampling of Neotropical Costus. In addition, we performed whole genome resequencing for a subset of 20 closely related species with contrasting pollination syndromes. For each of these 20 genomes, a high-quality assembled transcriptome was used as reference for consensus calling of candidate loci hypothesized to be associated with pollination-related traits of interest. To test for the role these candidate genes may play in evolutionary shifts in pollinators, signatures of selection were estimated as dN/dS across the identified candidate loci. We obtained a well-resolved phylogeny for Neotropical Costus despite conflict among gene trees that provide evidence of incomplete lineage sorting and/or reticulation. The overall topology and the network of genome-wide single nucleotide polymorphisms (SNPs) indicate that multiple shifts in pollination strategy have occurred across Costus, while also suggesting the presence of previously undetected signatures of hybridization between distantly related taxa. Traits related to pollination syndromes are strongly correlated and have been gained and lost in concert several times throughout the evolution of the genus. The presence of bract appendages is correlated with two traits associated with defenses against herbivory. Although labellum shape is strongly correlated with overall pollination syndrome, we found no significant impact of labellum shape on diversification rates. Evidence suggests an interplay of pollination success with other selective pressures shaping the evolution of the Costus inflorescence. Although most of the loci used for phylogenetic inference appear to be under purifying selection, many candidate genes associated with functional traits show evidence of being under positive selection. Together these results indicate an interplay of phylogenetic history with adaptive evolution leading to the diversification of pollination-associated traits in Neotropical Costus.
Collapse
Affiliation(s)
- Eugenio Valderrama
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | - Jacob B. Landis
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
- BTI Computational Biology Center, Boyce Thompson Institute, Ithaca, NY, United States
| | - Dave Skinner
- Le Jardin Ombragé, Tallahassee, FL, United States
| | - Paul J. M. Maas
- Section Botany, Naturalis Biodiversity Center, Leiden, Netherlands
| | | | - Thiago André
- Departamento de Botânica, Instituto de Ciências Biológicas, Universidade de Brasília, Brasília, DF, Brazil
| | - Nikolaus Grunder
- Department of Biological Sciences, California State University, East Bay, Hayward, CA, United States
| | - Chodon Sass
- University and Jepson Herbaria, University of California, Berkeley, Berkeley, CA, United States
| | - Maria Pinilla-Vargas
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | - Clarice J. Guan
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | - Heather R. Phillips
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | | | - Chelsea D. Specht
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| |
Collapse
|
6
|
Hewitt TL, Haponski AE, Foighil DÓ. Evolution of diverse host infection mechanisms delineates an adaptive radiation of lampsiline freshwater mussels centered on their larval ecology. PeerJ 2021; 9:e12287. [PMID: 34820162 PMCID: PMC8603817 DOI: 10.7717/peerj.12287] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Accepted: 09/21/2021] [Indexed: 11/25/2022] Open
Abstract
North American watersheds contain a high diversity of freshwater mussels (Unionoida). During the long-lived, benthic phase of their life cycle, up to 40 species can co-occur in a single riffle and there is typically little evidence for major differences in their feeding ecology or microhabitat partitioning. In contrast, their brief parasitic larval phase involves the infection of a wide diversity of fish hosts and female mussels have evolved a spectrum of adaptations for infecting host fish with their offspring. Many species use a passive broadcast strategy: placing high numbers of larvae in the water column and relying on chance encounters with potential hosts. Many other species, including most members of the Lampsilini, have a proactive strategy that entails the use of prey-mimetic lures to change the behavior of the hosts, i.e., eliciting a feeding response through which they become infected. Two main lure types are collectively produced: mantle tissue lures (on the female’s body) and brood lures, containing infective larvae, that are released into the external environment. In this study, we used a phylogenomic approach (ddRAD-seq) to place the diversity of infection strategies used by 54 North American lampsiline mussels into an evolutionary context. Ancestral state reconstruction recovered evidence for the early evolution of mantle lures in this clade, with brood lures and broadcast infection strategies both being independently derived twice. The most common infection strategy, occurring in our largest ingroup clade, is a mixed one in which mimetic mantle lures are apparently the predominant infection mechanism, but gravid females also release simple, non-mimetic brood lures at the end of the season. This mixed infection strategy clade shows some evidence of an increase in diversification rate and most members use centrarchids (Micropterus & Lepomis spp.) as their predominant fish hosts. Broad linkage between infection strategies and predominant fish host genera is also seen in other lampsiline clades: worm-like mantle lures of Toxolasma spp. with sunfish (Lepomis spp.); insect larvae-like brood lures (Ptychobranchus spp.), or mantle lures (Medionidus spp., Obovaria spp.), or mantle lures combined with host capture (Epioblasma spp.) with a spectrum of darter (Etheostoma & Percina spp.) and sculpin (Cottus spp.) hosts, and tethered brood lures (Hamiota spp.) with bass (Micropterus spp.). Our phylogenetic results confirm that discrete lampsiline mussel clades exhibit considerable specialization in the primary fish host clades their larvae parasitize, and in the host infection strategies they employ to do so. They are also consistent with the hypothesis that larval resource partitioning of fish hosts is an important factor in maintaining species diversity in mussel assemblages. We conclude that, taking their larval ecology and host-infection mechanisms into account, lampsiline mussels may be legitimately viewed as an adaptive radiation.
Collapse
Affiliation(s)
- Trevor L Hewitt
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan - Ann Arbor, Ann Arbor, MI, United States of America
| | - Amanda E Haponski
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan - Ann Arbor, Ann Arbor, MI, United States of America
| | - Diarmaid Ó Foighil
- Department of Ecology and Evolutionary Biology and Museum of Zoology, University of Michigan - Ann Arbor, Ann Arbor, MI, United States of America
| |
Collapse
|
7
|
FENG JEAN, DEWITT WILLIAMS, MCKENNA AARON, SIMON NOAH, WILLIS AMYD, MATSEN FREDERICKA. ESTIMATION OF CELL LINEAGE TREES BY MAXIMUM-LIKELIHOOD PHYLOGENETICS. Ann Appl Stat 2021; 15:343-362. [PMID: 35990087 PMCID: PMC9387344 DOI: 10.1214/20-aoas1400] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2023]
Abstract
CRISPR technology has enabled cell lineage tracing for complex multicellular organisms through insertion-deletion mutations of synthetic genomic barcodes during organismal development. To reconstruct the cell lineage tree from the mutated barcodes, current approaches apply general-purpose computational tools that are agnostic to the mutation process and are unable to take full advantage of the data's structure. We propose a statistical model for the CRISPR mutation process and develop a procedure to estimate the resulting tree topology, branch lengths, and mutation parameters by iteratively applying penalized maximum likelihood estimation. By assuming the barcode evolves according to a molecular clock, our method infers relative ordering across parallel lineages, whereas existing techniques only infer ordering for nodes along the same lineage. When analyzing transgenic zebrafish data from McKenna, Findlay and Gagnon et al. (2016), we find that our method recapitulates known aspects of zebrafish development and the results are consistent across samples.
Collapse
Affiliation(s)
- JEAN FENG
- Department of Epidemiology and Biostatistics, University of California, San Francisco
| | | | - AARON MCKENNA
- Department of Molecular and Systems Biology, Dartmouth College
| | - NOAH SIMON
- Department of Biostatistics, University of Washington
| | - AMY D. WILLIS
- Department of Biostatistics, University of Washington
| | | |
Collapse
|
8
|
Jandzik D, Stock DW. Differences in developmental potential predict the contrasting patterns of dental diversification in characiform and cypriniform fishes. Proc Biol Sci 2021; 288:20202205. [PMID: 33563123 PMCID: PMC7893225 DOI: 10.1098/rspb.2020.2205] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 01/12/2021] [Indexed: 11/12/2022] Open
Abstract
Morphological diversification during adaptive radiation may depend on factors external or internal to the lineage. We provide evidence for the latter in characiform fishes (tetras and piranhas), which exhibit extensive dental diversity. Phylogenetic character mapping supported regain of lost teeth as contributing to this diversity. To test for latent potential for dentition that would facilitate its evolutionary expansion, we overexpressed a tooth initiation signal, the tumour necrosis factor pathway ligand ectodysplasin, in a model characiform, the Mexican tetra (Astyanax mexicanus). This manipulation resulted in extensive ectopic dentition, in contrast with its previously reported limited effect in the zebrafish (Danio rerio). Tooth location in the order Cypriniformes, to which the zebrafish belongs, is much more restricted than in characiforms, a pattern that may be explained by differences in the retention of ancestral developmental potential. Our results suggest that differences in evolvability between lineages may lead to contrasting patterns of diversification.
Collapse
Affiliation(s)
- David Jandzik
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO 80309, USA
- Department of Zoology, Comenius University in Bratislava, Bratislava 84215, Slovakia
| | - David W. Stock
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO 80309, USA
| |
Collapse
|
9
|
Meerow AW, Gardner EM, Nakamura K. Phylogenomics of the Andean Tetraploid Clade of the American Amaryllidaceae (Subfamily Amaryllidoideae): Unlocking a Polyploid Generic Radiation Abetted by Continental Geodynamics. FRONTIERS IN PLANT SCIENCE 2020; 11:582422. [PMID: 33250911 PMCID: PMC7674842 DOI: 10.3389/fpls.2020.582422] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2020] [Accepted: 10/12/2020] [Indexed: 05/27/2023]
Abstract
One of the two major clades of the endemic American Amaryllidaceae subfam. Amaryllidoideae constitutes the tetraploid-derived (n = 23) Andean-centered tribes, most of which have 46 chromosomes. Despite progress in resolving phylogenetic relationships of the group with plastid and nrDNA, certain subclades were poorly resolved or weakly supported in those previous studies. Sequence capture using anchored hybrid enrichment was employed across 95 species of the clade along with five outgroups and generated sequences of 524 nuclear genes and a partial plastome. Maximum likelihood phylogenetic analyses were conducted on concatenated supermatrices, and coalescent-based species tree analyses were run on the gene trees, followed by hybridization network, age diversification and biogeographic analyses. The four tribes Clinantheae, Eucharideae, Eustephieae, and Hymenocallideae (sister to Clinantheae) are resolved in all analyses with > 90 and mostly 100% support, as are almost all genera within them. Nuclear gene supermatrix and species tree results were largely in concordance; however, some instances of cytonuclear discordance were evident. Hybridization network analysis identified significant reticulation in Clinanthus, Hymenocallis, Stenomesson and the subclade of Eucharideae comprising Eucharis, Caliphruria, and Urceolina. Our data support a previous treatment of the latter as a single genus, Urceolina, with the addition of Eucrosia dodsonii. Biogeographic analysis and penalized likelihood age estimation suggests an origin in the Cauca, Desert and Puna Neotropical bioprovinces for the complex in the mid-Oligocene, with more dispersals than vicariances in its history, but no extinctions. Hymenocallis represents the only instance of long-distance vicariance from the tropical Andean origin of its tribe Hymenocallideae. The absence of extinctions correlates with the lack of diversification rate shifts within the clade. The Eucharideae experienced a sudden lineage radiation ca. 10 Mya. We tie much of the divergences in the Andean-centered lineages to the rise of the Andes, and suggest that the Amotape-Huancabamba Zone functioned as both a corridor (dispersal) and a barrier to migration (vicariance). Several taxonomic changes are made. This is the largest DNA sequence data set to be applied within Amaryllidaceae to date.
Collapse
Affiliation(s)
- Alan W. Meerow
- USDA-ARS-SHRS, National Clonal Germplasm Repository, Miami, FL, United States
| | - Elliot M. Gardner
- Singapore Botanic Gardens, National Parks Board, Singapore, Singapore
- Institute of Environment, Florida International University, Miami, FL, United States
| | - Kyoko Nakamura
- USDA-ARS-SHRS, National Clonal Germplasm Repository, Miami, FL, United States
| |
Collapse
|
10
|
Valderrama E, Sass C, Pinilla-Vargas M, Skinner D, Maas PJM, Maas-van de Kamer H, Landis JB, Guan CJ, Specht CD. Unraveling the Spiraling Radiation: A Phylogenomic Analysis of Neotropical Costus L. FRONTIERS IN PLANT SCIENCE 2020; 11:1195. [PMID: 32922414 PMCID: PMC7456938 DOI: 10.3389/fpls.2020.01195] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Accepted: 07/23/2020] [Indexed: 06/01/2023]
Abstract
The family of pantropical spiral gingers (Costaceae Nakai; c. 125 spp.) can be used as a model to enhance our understanding of the mechanisms underlying Neotropical diversity. Costaceae has higher taxonomic diversity in South and Central America (c. 72 Neotropical species, c. 30 African, c. 23 Southeast Asian), particularly due to a radiation of Neotropical species of the genus Costus L. (c. 57 spp.). However, a well-supported phylogeny of the Neotropical spiral gingers including thorough sampling of proposed species encompassing their full morphologic and geographic variation is lacking, partly due to poor resolution recovered in previous analyses using a small sampling of loci. Here we use a phylogenomic approach to estimate the phylogeny of a sample of Neotropical Costus species using a targeted enrichment approach. Baits were designed to capture conserved elements' variable at the species level using available genomic sequences of Costus species and relatives. We obtained 832 loci (generating 791,954 aligned base pairs and 31,142 parsimony informative sites) for samples that encompassed the geographical and/or morphological diversity of some recognized species. Higher support values that improve the results of previous studies were obtained when including all the available loci, even those producing unresolved gene trees and having a low proportion of variable sites. Concatenation and coalescent-based species trees methods converge in almost the same topology suggesting a robust estimation of the relationships, even under the high levels of gene tree conflict presented here. The bait set design here presented made inferring a robust phylogeny to test taxonomic hypotheses possible and will improve our understanding of the origins of the charismatic diversity of the Neotropical spiral gingers.
Collapse
Affiliation(s)
- Eugenio Valderrama
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | - Chodon Sass
- The University and Jepson Herbaria, University of California, Berkeley, Berkeley, CA, United States
| | - Maria Pinilla-Vargas
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | | | - Paul J. M. Maas
- Section Botany, Naturalis Biodiversity Center, Leiden, Netherlands
| | | | - Jacob B. Landis
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | - Clarice J. Guan
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| | - Chelsea D. Specht
- School of Integrative Plant Science, Section of Plant Biology and the L.H. Bailey Hortorium, Cornell University, Ithaca, NY, United States
| |
Collapse
|
11
|
Zhang C, Dinh VU, Matsen FA. Nonbifurcating Phylogenetic Tree Inference via the Adaptive LASSO. J Am Stat Assoc 2020; 116:858-873. [PMID: 34305211 DOI: 10.1080/01621459.2020.1778481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
Phylogenetic tree inference using deep DNA sequencing is reshaping our understanding of rapidly evolving systems, such as the within-host battle between viruses and the immune system. Densely sampled phylogenetic trees can contain special features, including sampled ancestors in which we sequence a genotype along with its direct descendants, and polytomies in which multiple descendants arise simultaneously. These features are apparent after identifying zero-length branches in the tree. However, current maximum-likelihood based approaches are not capable of revealing such zero-length branches. In this paper, we find these zero-length branches by introducing adaptive-LASSO-type regularization estimators for the branch lengths of phylogenetic trees, deriving their properties, and showing regularization to be a practically useful approach for phylogenetics.
Collapse
Affiliation(s)
- Cheng Zhang
- School of Mathematical Sciences and Center for Statistical Science, Peking University
| | - V U Dinh
- Department of Mathematical Sciences, University of Delaware
| | | |
Collapse
|
12
|
Tiny Plants with Enormous Potential: Phylogeny and Evolution of Duckweeds. THE DUCKWEED GENOMES 2020. [DOI: 10.1007/978-3-030-11045-1_2] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
13
|
Grover CE, Arick MA, Thrash A, Conover JL, Sanders WS, Peterson DG, Frelichowski JE, Scheffler JA, Scheffler BE, Wendel JF. Insights into the Evolution of the New World Diploid Cottons (Gossypium, Subgenus Houzingenia) Based on Genome Sequencing. Genome Biol Evol 2019; 11:53-71. [PMID: 30476109 PMCID: PMC6320677 DOI: 10.1093/gbe/evy256] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/20/2018] [Indexed: 12/24/2022] Open
Abstract
We employed phylogenomic methods to study molecular evolutionary processes and phylogeny in the geographically widely dispersed New World diploid cottons (Gossypium, subg. Houzingenia). Whole genome resequencing data (average of 33× genomic coverage) were generated to reassess the phylogenetic history of the subgenus and provide a temporal framework for its diversification. Phylogenetic analyses indicate that the subgenus likely originated following transoceanic dispersal from Africa about 6.6 Ma, but that nearly all of the biodiversity evolved following rapid diversification in the mid-Pleistocene (0.5-2.0 Ma), with multiple long-distance dispersals required to account for range expansion to Arizona, the Galapagos Islands, and Peru. Comparative analyses of cpDNAversus nuclear data indicate that this history was accompanied by several clear cases of interspecific introgression. Repetitive DNAs contribute roughly half of the total 880 Mb genome, but most transposable element families are relatively old and stable among species. In the genic fraction, pairwise synonymous mutation rates average 1% per Myr, with nonsynonymous changes being about seven times less frequent. Over 1.1 million indels were detected and phylogenetically polarized, revealing a 2-fold bias toward deletions over small insertions. We suggest that this genome down-sizing bias counteracts genome size growth by TE amplification and insertions, and helps explain the relatively small genomes that are restricted to this subgenus. Compared with the rate of nucleotide substitution, the rate of indel occurrence is much lower averaging about 17 nucleotide substitutions per indel event.
Collapse
Affiliation(s)
- Corrinne E Grover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University
| | - Mark A Arick
- Institute for Genomics, Biocomputing, and Biotechnology, Mississippi State University
| | - Adam Thrash
- Institute for Genomics, Biocomputing, and Biotechnology, Mississippi State University
| | - Justin L Conover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University
| | - William S Sanders
- Institute for Genomics, Biocomputing, and Biotechnology, Mississippi State University
- Department of Computer Science & Engineering, Mississippi State University
- The Jackson Laboratory, Connecticut
| | - Daniel G Peterson
- Institute for Genomics, Biocomputing, and Biotechnology, Mississippi State University
| | | | | | - Brian E Scheffler
- USDA, Genomics and Bioinformatics Research Unit, Stoneville, Mississippi
| | - Jonathan F Wendel
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University
| |
Collapse
|
14
|
Conover JL, Karimi N, Stenz N, Ané C, Grover CE, Skema C, Tate JA, Wolff K, Logan SA, Wendel JF, Baum DA. A Malvaceae mystery: A mallow maelstrom of genome multiplications and maybe misleading methods? JOURNAL OF INTEGRATIVE PLANT BIOLOGY 2019; 61:12-31. [PMID: 30474311 DOI: 10.1111/jipb.12746] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/23/2018] [Accepted: 11/18/2018] [Indexed: 06/09/2023]
Abstract
Previous research suggests that Gossypium has undergone a 5- to 6-fold multiplication following its divergence from Theobroma. However, the number of events, or where they occurred in the Malvaceae phylogeny remains unknown. We analyzed transcriptomic and genomic data from representatives of eight of the nine Malvaceae subfamilies. Phylogenetic analysis of nuclear data placed Dombeya (Dombeyoideae) as sister to the rest of Malvadendrina clade, but the plastid DNA tree strongly supported Durio (Helicteroideae) in this position. Intraspecific Ks plots indicated that all sampled taxa, except Theobroma (Byttnerioideae), Corchorus (Grewioideae), and Dombeya (Dombeyoideae), have experienced whole genome multiplications (WGMs). Quartet analysis suggested WGMs were shared by Malvoideae-Bombacoideae and Sterculioideae-Tilioideae, but did not resolve whether these are shared with each other or Helicteroideae (Durio). Gene tree reconciliation and Bayesian concordance analysis suggested a complex history. Alternative hypotheses are suggested, each involving two independent autotetraploid and one allopolyploid event. They differ in that one entails an allopolyploid origin for the Durio lineage, whereas the other invokes an allopolyploid origin for Malvoideae-Bombacoideae. We highlight the need for more genomic information in the Malvaceae and improved methods to resolve complex evolutionary histories that may include allopolyploidy, incomplete lineage sorting, and variable rates of gene and genome evolution.
Collapse
Affiliation(s)
- Justin L Conover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, 50011, USA
| | - Nisa Karimi
- Department of Botany, University of Wisconsin Madison, WI, 53706, USA
| | - Noah Stenz
- Department of Botany, University of Wisconsin Madison, WI, 53706, USA
| | - Cécile Ané
- Department of Botany, University of Wisconsin Madison, WI, 53706, USA
- Department of Statistics, University of Wisconsin Madison, WI, 53706, USA
| | - Corrinne E Grover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, 50011, USA
| | - Cynthia Skema
- Morris Arboretum of the University of Pennsylvania, 100 E. Northwestern Avenue, Philadelphia, PA, 19118, USA
| | - Jennifer A Tate
- Institute of Fundamental Sciences, Massey University, Palmerston North 4442, New Zealand
| | - Kirsten Wolff
- School of Natural and Environmental Sciences, Newcastle University, NE1 7RU, UK
| | - Samuel A Logan
- School of Natural and Environmental Sciences, Newcastle University, NE1 7RU, UK
| | - Jonathan F Wendel
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, 50011, USA
| | - David A Baum
- Department of Botany, University of Wisconsin Madison, WI, 53706, USA
- Wisconsin Institute for Discovery, 330 N Orchard St, Madison, WI 53715, USA
| |
Collapse
|
15
|
Clavel J, Aristide L, Morlon H. A Penalized Likelihood Framework for High-Dimensional Phylogenetic Comparative Methods and an Application to New-World Monkeys Brain Evolution. Syst Biol 2018; 68:93-116. [PMID: 29931145 DOI: 10.1093/sysbio/syy045] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2017] [Accepted: 06/13/2018] [Indexed: 01/03/2023] Open
Abstract
Working with high-dimensional phylogenetic comparative data sets is challenging because likelihood-based multivariate methods suffer from low statistical performances as the number of traits $p $ approaches the number of species $n $ and because some computational complications occur when $p $ exceeds $n$. Alternative phylogenetic comparative methods have recently been proposed to deal with the large $p $ small $n $ scenario but their use and performances are limited. Herein, we develop a penalized likelihood (PL) framework to deal with high-dimensional comparative data sets. We propose various penalizations and methods for selecting the intensity of the penalties. We apply this general framework to the estimation of parameters (the evolutionary trait covariance matrix and parameters of the evolutionary model) and model comparison for the high-dimensional multivariate Brownian motion (BM), Early-burst (EB), Ornstein-Uhlenbeck (OU), and Pagel's lambda models. We show using simulations that our PL approach dramatically improves the estimation of evolutionary trait covariance matrices and model parameters when $p$ approaches $n$, and allows for their accurate estimation when $p$ equals or exceeds $n$. In addition, we show that PL models can be efficiently compared using generalized information criterion (GIC). We implement these methods, as well as the related estimation of ancestral states and the computation of phylogenetic principal component analysis in the R package RPANDA and mvMORPH. Finally, we illustrate the utility of the new proposed framework by evaluating evolutionary models fit, analyzing integration patterns, and reconstructing evolutionary trajectories for a high-dimensional 3D data set of brain shape in the New World monkeys. We find a clear support for an EB model suggesting an early diversification of brain morphology during the ecological radiation of the clade. PL offers an efficient way to deal with high-dimensional multivariate comparative data.
Collapse
Affiliation(s)
- Julien Clavel
- École Normale Supérieure, Paris Sciences et Lettres (PSL) Research University, Institut de Biologie de l'École Normale Supérieure (IBENS), CNRS UMR 8197, INSERM U1024, 46 rue d'Ulm, F-75005 Paris, France
| | - Leandro Aristide
- École Normale Supérieure, Paris Sciences et Lettres (PSL) Research University, Institut de Biologie de l'École Normale Supérieure (IBENS), CNRS UMR 8197, INSERM U1024, 46 rue d'Ulm, F-75005 Paris, France
| | - Hélène Morlon
- École Normale Supérieure, Paris Sciences et Lettres (PSL) Research University, Institut de Biologie de l'École Normale Supérieure (IBENS), CNRS UMR 8197, INSERM U1024, 46 rue d'Ulm, F-75005 Paris, France
| |
Collapse
|
16
|
Diversity begets diversity in competition for space. Nat Ecol Evol 2017; 1:156. [DOI: 10.1038/s41559-017-0156] [Citation(s) in RCA: 65] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2016] [Accepted: 04/04/2017] [Indexed: 11/09/2022]
|
17
|
Naumenko SA, Logacheva MD, Popova NV, Klepikova AV, Penin AA, Bazykin GA, Etingova AE, Mugue NS, Kondrashov AS, Yampolsky LY. Transcriptome‐based phylogeny of endemic Lake Baikal amphipod species flock: fast speciation accompanied by frequent episodes of positive selection. Mol Ecol 2017; 26:536-553. [DOI: 10.1111/mec.13927] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2016] [Revised: 10/11/2016] [Accepted: 10/12/2016] [Indexed: 12/15/2022]
Affiliation(s)
- Sergey A. Naumenko
- Belozersky Institute of Physico‐Chemical Biology Lomonosov Moscow State University Moscow Russia
- Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences Moscow Russia
- Genetics and Genome Biology Program The Hospital For Sick Children Toronto ON Canada
| | - Maria D. Logacheva
- Belozersky Institute of Physico‐Chemical Biology Lomonosov Moscow State University Moscow Russia
- Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences Moscow Russia
- Pirogov Russian National Research Medical University Moscow Russia
| | - Nina V. Popova
- Belozersky Institute of Physico‐Chemical Biology Lomonosov Moscow State University Moscow Russia
| | - Anna V. Klepikova
- Belozersky Institute of Physico‐Chemical Biology Lomonosov Moscow State University Moscow Russia
- Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences Moscow Russia
| | - Aleksey A. Penin
- Belozersky Institute of Physico‐Chemical Biology Lomonosov Moscow State University Moscow Russia
- Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences Moscow Russia
| | - Georgii A. Bazykin
- Belozersky Institute of Physico‐Chemical Biology Lomonosov Moscow State University Moscow Russia
- Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences Moscow Russia
- Pirogov Russian National Research Medical University Moscow Russia
- Skolkovo Institute of Science and Technology Skolkovo Russia
| | - Anna E. Etingova
- Baikal Museum Irkutsk Research Center Russian Academy of Sciences Listvyanka, Irkutsk region Russia
| | - Nikolai S. Mugue
- Laboratory of Molecular Genetics Russian Institute for Fisheries and Oceanography (VNIRO) Moscow Russia
- Laboratory of Experimental Embryology Koltsov Institute of Developmental Biology Moscow Russia
| | - Alexey S. Kondrashov
- Belozersky Institute of Physico‐Chemical Biology Lomonosov Moscow State University Moscow Russia
- Department of Ecology and Evolution University of Michigan Ann Arbor MI USA
| | - Lev Y. Yampolsky
- Department of Biological Sciences East Tennessee State University Johnson City TN USA
| |
Collapse
|
18
|
Potter GE, Hens N. A penalized likelihood approach to estimate within-household contact networks from egocentric data. J R Stat Soc Ser C Appl Stat 2013; 62:629-648. [PMID: 23935218 DOI: 10.1111/rssc.12011] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Acute infectious diseases are transmitted over networks of social contacts. Epidemic models are used to predict the spread of emergent pathogens and compare intervention strategies. Many of these models assume equal probability of contact within mixing groups (homes, schools, etc.), but little work has inferred the actual contact network, which may influence epidemic estimates. We develop a penalized likelihood method to infer contact networks within households, a key area for disease transmission. Using egocentric surveys of contact behavior in Belgium, we estimate within-household contact networks for six different age compositions. Our estimates show dependency in contact behavior and vary substantively by age composition, with fewer contacts occurring in older households. Our results are relevant for epidemic models used to make policy recommendations.
Collapse
Affiliation(s)
- Gail E Potter
- California Polytechnic State University, San Luis Obispo, CA, U.S.A., Center for Statistics and Quantitative Infectious Diseases, Fred Hutchinson Cancer Research Center, Seattle, WA, U.S.A
| | | |
Collapse
|
19
|
Paradis E. Molecular dating of phylogenies by likelihood methods: a comparison of models and a new information criterion. Mol Phylogenet Evol 2013; 67:436-44. [PMID: 23454091 DOI: 10.1016/j.ympev.2013.02.008] [Citation(s) in RCA: 96] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2012] [Revised: 01/18/2013] [Accepted: 02/08/2013] [Indexed: 10/27/2022]
Abstract
Dating the divergence in a phylogenetic tree is a fundamental step in evolutionary analysis. Some extensions and improvements of the penalised likelihood method originally presented by Sanderson are introduced. The improvements are the introduction of alternative models, including one with non-correlated rates of molecular substitution ("relaxed" model), a completely reworked fitting algorithm that considers the high-dimensionality of the optimisation problem, and the development of a new information criterion for model selection in the presence of a penalised term. It is also shown that the strict clock model is a special case of the present approach. An extensive simulation study was conducted to assess the statistical performance of these improvements. Overall, the different estimators studied here appeared as unbiased though their variance varied depending on the fitted and the simulated models and on the number of calibration points. The strict clock model gave good estimates of branch lengths even in the presence of heterogeneous substitution rates. The correlated model gave the best estimates of substitution rates whatever the model used to simulate the data. These results, which are certainly the first from an extensive simulation study of a molecular dating method, call for more comparison with alternative methods, as well as further work on the developments introduced here.
Collapse
Affiliation(s)
- Emmanuel Paradis
- Institut de Recherche pour le Développement, ISEM UMR 226/5554 - UM2/CNRS/IRD, Jl. Taman Kemang 32B, Jakarta 12730, Indonesia.
| |
Collapse
|
20
|
Huelsenbeck JP, Alfaro ME, Suchard MA. Biologically inspired phylogenetic models strongly outperform the no common mechanism model. Syst Biol 2011; 60:225-32. [PMID: 21252385 PMCID: PMC3038349 DOI: 10.1093/sysbio/syq089] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2009] [Revised: 06/29/2009] [Accepted: 09/22/2010] [Indexed: 11/13/2022] Open
Abstract
But Tuffley and Steel (1997) introduced a model called No Common Mechanism (NCM), in which characters may-but are not required to-vary their relative rates independently, both within and between branches. Because the independent variation is taken only as a possibility, not as a requirement, NCM would apply to almost any situation, and so may be accepted as realistic. This is useful because Tuffley and Steel also showed that maximum likelihood under NCM selects the same trees as does parsimony. With the realistic NCM in the background, then, most parsimonious trees have greatest power to explain available observations. -Farris (2008).
Collapse
Affiliation(s)
- John P Huelsenbeck
- Department of Integrative Biology, University of California, Berkeley, CA 94720-3140, USA.
| | | | | |
Collapse
|
21
|
Randle CP, Pickett KM. The conflation of ignorance and knowledge in the inference of clade posteriors. Cladistics 2010; 26:550-559. [DOI: 10.1111/j.1096-0031.2009.00301.x] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
|
22
|
Holder MT, Lewis PO, Swofford DL. The akaike information criterion will not choose the no common mechanism model. Syst Biol 2010; 59:477-85. [PMID: 20547783 DOI: 10.1093/sysbio/syq028] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mark T Holder
- Department of Ecology and Evolutionary Biology, University of Kansas, 1200 Sunnyside Avenue, Lawrence, KS 66045, USA.
| | | | | |
Collapse
|
23
|
|
24
|
Abstract
Heterotachy is a general term to describe positions in a sequence that evolve at different rates in different lineages. Kolaczkowski and Thornton (2004. Performance of maximum parsimony and likelihood phylogenetics when evolution is heterogeneous. Nature 431:980-984.) recently described an intriguing heterotachy model that leads to topological bias for likelihood-based methods and parsimony methods. In this article, we show that heterotachy can generally be viewed as multivariate rates-across-sites variation, which can be described as randomly drawing rates (or branch lengths) from a multivariate distribution for each branch at each site. Motivated by this idea, we propose a pairwise alpha heterotachy adjustment model, which gives us much improved topological estimation in the settings by Kolaczkowski and Thornton (2004).
Collapse
Affiliation(s)
- Jihua Wu
- Department of Mathematics and Statistics, Dalhousie University, Halifax, Nova Scotia, Canada.
| | | |
Collapse
|