skip to main content


Title: Phylogenomics reveals the origin of mammal lice out of Afrotheria
Mammals host a wide diversity of parasites. Lice, comprising more than 5,000 species, are one group of ectoparasites whose major lineages have a somewhat patchwork distribution across the major groups of mammals. Here we explored patterns in the diversification of mammalian lice by reconstructing a higher-level phylogeny of these lice, leveraging whole genome sequence reads to assemble single-copy orthologue genes across the genome. The evolutionary tree of lice indicated that three of the major lineages of placental mammal lice had a single common ancestor. Comparisons of this parasite phylogeny with that for their mammalian hosts indicated that the common ancestor of elephants, elephant shrews and hyraxes (that is, Afrotheria) was the ancestral host of this group of lice. Other groups of placental mammals obtained their lice via host-switching out of these Afrotherian ancestors. In addition, reconstructions of the ancestral host group (bird versus mammal) for all parasitic lice supported an avian ancestral host, indicating that the ancestor of Afrotheria acquired these parasites via host-switching from an ancient avian host. These results shed new light on the long-standing question of why the major groups of parasitic lice are not uniformly distributed across mammals and reveal the origins of mammalian lice.  more » « less
Award ID(s):
1925487 1926919
NSF-PAR ID:
10345580
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Nature Ecology & Evolution
Volume:
6
ISSN:
2397-334X
Page Range / eLocation ID:
1205-1210
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Lice are considered a model system for studying the process of cospeciation because they are obligate and permanent parasites and are often highly host‐specific. Among lice, species in the family Echinophthiriidae Enderlein (Anoplura) are unique in that they infest mammalian hosts with an amphibious lifestyle, i.e. pinnipeds and the river otter. There is evidence that the ancestor of this group infested the terrestrial ancestor of pinnipeds, which suggests these parasites coevolved with their hosts during the transition to marine environments. However, there has been no previous study investigating the phylogenetic relationships among sucking lice parasitizing seals and sea lions. To uncover the evolutionary history of these parasites, we obtained genomic data forAntarctophthirus microchirTrouessart and Neumann (from two hosts),Antarctophthirus carliniiLeonardiet al.,Antarctophthirus lobodontisEnderlein,Antarctophthirus ogmorhiniEnderlein,Lepidophthirus macrorhiniEnderlein, andProechinophthirus fluctusFerris. From genomic sequence reads, we assembled > 1000 nuclear genes and used these data to infer a phylogenetic tree for these lice. We also used the assembled genes in combination with read‐mapping to estimate heterozygosity and effective population size from individual lice. Our analysis supports the monophyly of lice from pinnipeds and uncovers phylogenetic relationships within the group. Surprisingly, we found thatA. carlinii,A. lobodontis, andA. ogmorhinihave very little genetic divergence among them, whereas the divergence between different geographic representatives ofA. microchirindicate that they are possibly different species. Nevertheless, our phylogeny of Echinophthiriidae suggests that these lice have consistently codiverged with their hosts with minimal host switching. Population genomic metrics indicate that louse effective population size is linked to host demographics, which further highlights the close association between pinnipeds and their lice.

     
    more » « less
  2. INTRODUCTION Resolving the role that different environmental forces may have played in the apparent explosive diversification of modern placental mammals is crucial to understanding the evolutionary context of their living and extinct morphological and genomic diversity. RATIONALE Limited access to whole-genome sequence alignments that sample living mammalian biodiversity has hampered phylogenomic inference, which until now has been limited to relatively small, highly constrained sequence matrices often representing <2% of a typical mammalian genome. To eliminate this sampling bias, we used an alignment of 241 whole genomes to comprehensively identify and rigorously analyze noncoding, neutrally evolving sequence variation in coalescent and concatenation-based phylogenetic frameworks. These analyses were followed by validation with multiple classes of phylogenetically informative structural variation. This approach enabled the generation of a robust time tree for placental mammals that evaluated age variation across hundreds of genomic loci that are not restricted by protein coding annotations. RESULTS Coalescent and concatenation phylogenies inferred from multiple treatments of the data were highly congruent, including support for higher-level taxonomic groupings that unite primates+colugos with treeshrews (Euarchonta), bats+cetartiodactyls+perissodactyls+carnivorans+pangolins (Scrotifera), all scrotiferans excluding bats (Fereuungulata), and carnivorans+pangolins with perissodactyls (Zooamata). However, because these approaches infer a single best tree, they mask signatures of phylogenetic conflict that result from incomplete lineage sorting and historical hybridization. Accordingly, we also inferred phylogenies from thousands of noncoding loci distributed across chromosomes with historically contrasting recombination rates. Throughout the radiation of modern orders (such as rodents, primates, bats, and carnivores), we observed notable differences between locus trees inferred from the autosomes and the X chromosome, a pattern typical of speciation with gene flow. We show that in many cases, previously controversial phylogenetic relationships can be reconciled by examining the distribution of conflicting phylogenetic signals along chromosomes with variable historical recombination rates. Lineage divergence time estimates were notably uniform across genomic loci and robust to extensive sensitivity analyses in which the underlying data, fossil constraints, and clock models were varied. The earliest branching events in the placental phylogeny coincide with the breakup of continental landmasses and rising sea levels in the Late Cretaceous. This signature of allopatric speciation is congruent with the low genomic conflict inferred for most superordinal relationships. By contrast, we observed a second pulse of diversification immediately after the Cretaceous-Paleogene (K-Pg) extinction event superimposed on an episode of rapid land emergence. Greater geographic continuity coupled with tumultuous climatic changes and increased ecological landscape at this time provided enhanced opportunities for mammalian diversification, as depicted in the fossil record. These observations dovetail with increased phylogenetic conflict observed within clades that diversified in the Cenozoic. CONCLUSION Our genome-wide analysis of multiple classes of sequence variation provides the most comprehensive assessment of placental mammal phylogeny, resolves controversial relationships, and clarifies the timing of mammalian diversification. We propose that the combination of Cretaceous continental fragmentation and lineage isolation, followed by the direct and indirect effects of the K-Pg extinction at a time of rapid land emergence, synergistically contributed to the accelerated diversification rate of placental mammals during the early Cenozoic. The timing of placental mammal evolution. Superordinal mammalian diversification took place in the Cretaceous during periods of continental fragmentation and sea level rise with little phylogenomic discordance (pie charts: left, autosomes; right, X chromosome), which is consistent with allopatric speciation. By contrast, the Paleogene hosted intraordinal diversification in the aftermath of the K-Pg mass extinction event, when clades exhibited higher phylogenomic discordance consistent with speciation with gene flow and incomplete lineage sorting. 
    more » « less
  3. Abstract

    The Cyathocotylidae is a globally distributed family of digeneans parasitic as adults in fish, reptiles, birds and mammals in both freshwater and marine environments. Molecular phylogenetic analysis of interrelationships among cyathocotylids is lacking with only a few species included in previous studies. We used sequences of the nuclear 28S rRNA gene to examine phylogenetic affinities of 11 newly sequenced taxa of cyathocotylids and the closely related family Brauninidae collected from fish, reptiles, birds and dolphins from Australia, Southeast Asia, Europe, North America and South America. This is the first study to provide sequence data from adult cyathocotylids parasitic in fish and reptiles. Our analyses demonstrated that the members of the genusBraunina(family Brauninidae) belong to the Cyathocotylidae, placing the Brauninidae into synonymy with the Cyathocotylidae. In addition, our DNA sequences supported the presence of a second species in the currently monotypicBraunina. Our phylogeny revealed thatCyathocotylespp. from crocodilians belong to a separate genus (Suchocyathocotyle, previously proposed as a subgenus) and subfamily (Suchocyathocotylinae subfam. n.). Morphological study ofGogatea serpentum indicumsupported its elevation to species asGogatea mehri. The phylogeny did not supportHolostephanoideswithin the subfamily Cyathocotylinae; instead,Holostephanoidesformed a strongly supported clade with members of the subfamily Szidatiinae (GogateaandNeogogatea). Therefore, we transferHolostephanoidesinto the Szidatiinae. DNA sequence data revealed the potential presence of cryptic species reported under the nameMesostephanus microbursa. Our phylogeny indicated at least two major host switching events in the evolutionary history of the subfamily Szidatiinae which likely resulted in the transition of these parasites from birds to fish and snakes. Likewise, the transition to dolphins byBrauninarepresents another major host switching event among the Cyathocotylidae. In addition, our phylogeny revealed more than a single transition between freshwater and marine environments demonstrated in our dataset byBrauninaand someMesostephanus.

     
    more » « less
  4. null (Ed.)
    Most animals have a conserved mitochondrial genome structure composed of a single chromosome. However, some organisms have their mitochondrial genes separated on several smaller circular or linear chromosomes. Highly fragmented circular chromosomes (“minicircles”) are especially prevalent in parasitic lice (Insecta: Phthiraptera), with 16 species known to have between nine and 20 mitochondrial minicircles per genome. All of these species belong to the same clade (mammalian lice), suggesting a single origin of drastic fragmentation. Nevertheless, other work indicates a lesser degree of fragmentation (2–3 chromosomes/genome) is present in some avian feather lice (Ischnocera: Philopteridae). In this study, we tested for minicircles in four species of the feather louse genus Columbicola (Philopteridae). Using whole genome shotgun sequence data, we applied three different bioinformatic approaches for assembling the Columbicola mitochondrial genome. We further confirmed these approaches by assembling the mitochondrial genome of Pediculus humanus from shotgun sequencing reads, a species known to have minicircles. Columbicola spp. genomes are highly fragmented into 15–17 minicircles between ∼1,100 and ∼3,100 bp in length, with 1–4 genes per minicircle. Subsequent annotation of the minicircles indicated that tRNA arrangements of minicircles varied substantially between species. These mitochondrial minicircles for species of Columbicola represent the first feather lice (Philopteridae) for which minicircles have been found in a full mitochondrial genome assembly. Combined with recent phylogenetic studies of parasitic lice, our results provide strong evidence that highly fragmented mitochondrial genomes, which are otherwise rare across the Tree of Life, evolved multiple times within parasitic lice. 
    more » « less
  5. INTRODUCTION Thousands of genetic variants have been associated with human diseases and traits through genome-wide association studies (GWASs). Translating these discoveries into improved therapeutics requires discerning which variants among hundreds of candidates are causally related to disease risk. To date, only a handful of causal variants have been confirmed. Here, we leverage 100 million years of mammalian evolution to address this major challenge. RATIONALE We compared genomes from hundreds of mammals and identified bases with unusually few variants (evolutionarily constrained). Constraint is a measure of functional importance that is agnostic to cell type or developmental stage. It can be applied to investigate any heritable disease or trait and is complementary to resources using cell type– and time point–specific functional assays like Encyclopedia of DNA Elements (ENCODE) and Genotype-Tissue Expression (GTEx). RESULTS Using constraint calculated across placental mammals, 3.3% of bases in the human genome are significantly constrained, including 57.6% of coding bases. Most constrained bases (80.7%) are noncoding. Common variants (allele frequency ≥ 5%) and low-frequency variants (0.5% ≤ allele frequency < 5%) are depleted for constrained bases (1.85 versus 3.26% expected by chance, P < 2.2 × 10 −308 ). Pathogenic ClinVar variants are more constrained than benign variants ( P < 2.2 × 10 −16 ). The most constrained common variants are more enriched for disease single-nucleotide polymorphism (SNP)–heritability in 63 independent GWASs. The enrichment of SNP-heritability in constrained regions is greater (7.8-fold) than previously reported in mammals and is even higher in primates (11.1-fold). It exceeds the enrichment of SNP-heritability in nonsynonymous coding variants (7.2-fold) and fine-mapped expression quantitative trait loci (eQTL)–SNPs (4.8-fold). The enrichment peaks near constrained bases, with a log-linear decrease of SNP-heritability enrichment as a function of the distance to a constrained base. Zoonomia constraint scores improve functionally informed fine-mapping. Variants at sites constrained in mammals and primates have greater posterior inclusion probabilities and higher per-SNP contributions. In addition, using both constraint and functional annotations improves polygenic risk score accuracy across a range of traits. Finally, incorporating constraint information into the analysis of noncoding somatic variants in medulloblastomas identifies new candidate driver genes. CONCLUSION Genome-wide measures of evolutionary constraint can help discern which variants are functionally important. This information may accelerate the translation of genomic discoveries into the biological, clinical, and therapeutic knowledge that is required to understand and treat human disease. Using evolutionary constraint in genomic studies of human diseases. ( A ) Constraint was calculated across 240 mammal species, including 43 primates (teal line). ( B ) Pathogenic ClinVar variants ( N = 73,885) are more constrained across mammals than benign variants ( N = 231,642; P < 2.2 × 10 −16 ). ( C ) More-constrained bases are more enriched for trait-associated variants (63 GWASs). ( D ) Enrichment of heritability is higher in constrained regions than in functional annotations (left), even in a joint model with 106 annotations (right). ( E ) Fine-mapping (PolyFun) using a model that includes constraint scores identifies an experimentally validated association at rs1421085. Error bars represent 95% confidence intervals. BMI, body mass index; LF, low frequency; PIP, posterior inclusion probability. 
    more » « less