Whale phylogeny and rapid radiation events revealed using novel retroposed elements and their flanking sequences

Improved mammalian family phylogeny using gap-rare multiple sequence alignment: A timetree of extant placentals and marsupials

The timing of mammalian diversification in relation to the Cretaceous-Paleogene (KPg) mass extinction continues to be a subject of substantial debate. Previous studies have either focused on limited taxonomic samples with available whole-genome data or relied on short sequence alignments coupled with extensive species samples. In the present study, we improved an existing dataset from the landmark study of Meredith et al. (2011) by filling in missing fragments and further generated another dataset containing 120 taxa and 98 exonic markers. Using these two datasets, we then constructed phylogenies for extant mammalian families, providing improved resolution of many conflicting relationships. Moreover, the timetrees generated, which were calibrated using appropriate molecular clock models and multiple fossil records, indicated that the interordinal diversification of placental mammals initiated before the Late Cretaceous period. Additionally, intraordinal diversification of both extant placental and marsupial lineages accelerated after the KPg boundary, supporting the hypothesis that the availability of numerous vacant ecological niches subsequent to the mass extinction event facilitated rapid diversification. Thus, our results support a scenario of placental radiation characterized by both basal cladogenesis and active interordinal divergences spanning from the Late Cretaceous into the Paleogene.

Homoplasy of Retrotransposon Insertions in Toothed Whales

Retrotransposon insertion patterns facilitate a virtually homoplasy-free picture of phylogenetic history. Still, a few most likely random parallel insertions or deletions result in rare cases of homoplasy in primates. The following question arises: how frequent is retrotransposon homoplasy in other phylogenetic clades? Here, we derived genome insertion data of toothed whales to evaluate the extension of homoplasy in a representative laurasiatherian group. Among more than a thousand extracted and aligned retrotransposon loci, we detected 37 cases of precise parallel insertions in species that are separated by over more than 10 million years, a time frame which minimizes the effects of incomplete lineage sorting. We compared the phylogenetic signal of insertions with the flanking sequences of these loci to further exclude potential polymorphic loci derived by incomplete lineage sorting. We found that the phylogenetic signals of retrotransposon insertion patterns exhibiting true homoplasy differ from the signals of their flanking sequences. In toothed whales, precise parallel insertions account for around 0.18-0.29% of insertion cases, which is about 12.5 times the frequency of such insertions among Alus in primates. We also detected five specific deletions of retrotransposons on various lineages of toothed whale evolution, a frequency of 0.003%, which is slightly higher than such occurrences in primates. Overall, the level of retrotransposon homoplasy in toothed whales is still marginal compared to the phylogenetic diagnostic retrotransposon presence/absence signal.

Vertebral trabecular bone mechanical properties vary among functional groups of cetaceans

Since their appearance in the fossil record 34 million years ago, modern cetaceans (dolphins, whales, and porpoises) have radiated into diverse habitats circumglobally, developing vast phenotypic variations among species. Traits such as skeletal morphology and ecologically linked behaviors denote swimming activity; trade-offs in flexibility and rigidity along the vertebral column determine patterns of caudal oscillation. Here, we categorized 10 species of cetaceans (families Delphinidae and Kogiidae; N = 21 animals) into functional groups based on vertebral centra morphology, swimming speeds, diving behavior, and inferred swimming patterns. We quantified trabecular bone mechanical properties (yield strength, apparent stiffness, and resilience) among functional groups and regions of the vertebral column (thoracic, lumbar, and caudal). We extracted 6 mm³ samples from vertebral bodies and tested them in compression in 3 orientations (rostrocaudal, dorsoventral, and mediolateral) at 2 mm min⁻¹. Overall, bone from the pre-fluke/fluke boundary had the greatest yield strength and resilience, indicating that the greatest forces are translated to the tail during caudal oscillatory swimming. Group 1, composed of 5 shallow-diving delphinid species, had the greatest vertebral trabecular bone yield strength, apparent stiffness, and resilience of all functional groups. Conversely, Group 3, composed of 2 deep-diving kogiid species, had the least strong, stiff, and resilient bone, while Group 2 (3 deep-diving delphinid species) exhibited intermediate values. These data suggest that species that incorporate prolonged glides during deep descents in the water column actively swim less, and place relatively smaller loads on their vertebral columns, compared with species that execute shallower dives. We found that cetacean vertebral trabecular bone properties differed from the properties of terrestrial mammals; for every given bone strength, cetacean bone was less stiff by comparison. This relative lack of material rigidity within vertebral bone may be attributed to the non-weight-bearing locomotor modes of fully aquatic mammals.

Extensive Interspecific Gene Flow Shaped Complex Evolutionary History and Underestimated Species Diversity in Rapidly Radiated Dolphins

Recently diverged taxa are often characterized by high rates of hybridization, which can complicate phylogenetic reconstruction. For this reason, the phylogenetic relationships and evolutionary history of dolphins are still not very well resolved; the question of whether the genera Tursiops and Stenella are monophyletic is especially controversial. Here, we performed re-sequencing of six dolphin genomes and combined them with eight previously published dolphin SRA datasets and six whole-genome datasets to investigate the phylogenetic relationships of dolphins and test the monophyly hypothesis of Tursiops and Stenella. Phylogenetic reconstruction with the maximum likelihood and Bayesian methods of concatenated loci, as well as with coalescence analyses of sliding window trees, produced a concordant and well-supported tree. Our studies support the non-monophyletic status of Tursiops and Stenella because the species referred these genera do not form exclusive monophyletic clades. This suggests that the current taxonomy of both genera might not reflect their evolutionary history and may underestimate their diversity. A four-taxon D-statistic (ABBA-BABA) test, five-taxon DFOIL test, and tree-based PhyloNet analyses all showed extensive gene flow across dolphin species, which could explain the instability in resolving phylogenetic relationship of oceanic dolphins with different and limited markers. This study could be a good case to demonstrate how genomic data can reveal complex speciation and phylogeny in rapidly radiating animal groups.

Comparative genomics provides insights into the aquatic adaptations of mammals

Divergent lineages can respond to common environmental factors through convergent processes involving shared genomic components or pathways, but the molecular mechanisms are poorly understood. Here, we provide genomic resources and insights into the evolution of mammalian lineages adapting to aquatic life. Our data suggest convergent evolution, for example, in association with thermoregulation through genes associated with a surface heat barrier (NFIA) and internal heat exchange (SEMA3E). Combined with the support of previous reports showing that the UCP1 locus has been lost in many marine mammals independently, our results suggest that the thermostatic strategy of marine mammals shifted from enhancing heat production to limiting heat loss.

The ancestors of marine mammals once roamed the land and independently committed to an aquatic lifestyle. These macroevolutionary transitions have intrigued scientists for centuries. Here, we generated high-quality genome assemblies of 17 marine mammals (11 cetaceans and six pinnipeds), including eight assemblies at the chromosome level. Incorporating previously published data, we reconstructed the marine mammal phylogeny and population histories and identified numerous idiosyncratic and convergent genomic variations that possibly contributed to the transition from land to water in marine mammal lineages. Genes associated with the formation of blubber (NFIA), vascular development (SEMA3E), and heat production by brown adipose tissue (UCP1) had unique changes that may contribute to marine mammal thermoregulation. We also observed many lineage-specific changes in the marine mammals, including genes associated with deep diving and navigation. Our study advances understanding of the timing, pattern, and molecular changes associated with the evolution of mammalian lineages adapting to aquatic life.

Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture

The evolution of cetaceans, from their early transition to an aquatic lifestyle to their subsequent diversification, has been the subject of numerous studies. However, although the higher-level relationships among cetacean families have been largely settled, several aspects of the systematics within these groups remain unresolved. Problematic clades include the oceanic dolphins (37 spp.), which have experienced a recent rapid radiation, and the beaked whales (22 spp.), which have not been investigated in detail using nuclear loci. The combined application of high-throughput sequencing with techniques that target specific genomic sequences provide a powerful means of rapidly generating large volumes of orthologous sequence data for use in phylogenomic studies. To elucidate the phylogenetic relationships within the Cetacea, we combined sequence capture with Illumina sequencing to generate data for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}3200 protein-coding genes for 68 cetacean species and their close relatives including the pygmy hippopotamus. By combining data from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$>$\end{document}38,000 exons with existing sequences from 11 cetaceans and seven outgroup taxa, we produced the first comprehensive comparative genomic data set for cetaceans, spanning 6,527,596 aligned base pairs (bp) and 89 taxa. Phylogenetic trees reconstructed with maximum likelihood and Bayesian inference of concatenated loci, as well as with coalescence analyses of individual gene trees, produced mostly concordant and well-supported trees. Our results completely resolve the relationships among beaked whales as well as the contentious relationships among oceanic dolphins, especially the problematic subfamily Delphinidae. We carried out Bayesian estimation of species divergence times using MCMCTree and compared our complete data set to a subset of clocklike genes. Analyses using the complete data set consistently showed less variance in divergence times than the reduced data set. In addition, integration of new fossils (e.g., Mystacodon selenensis) indicates that the diversification of Crown Cetacea began before the Late Eocene and the divergence of Crown Delphinidae as early as the Middle Miocene. [Cetaceans; phylogenomics; Delphinidae; Ziphiidae; dolphins; whales.]

Widespread positive selection on cetacean TLR extracellular domain

Toll like receptors (TLRs), key members of innate immune system, can recognize a wide diversity of pathogens and initiate both innate and adaptive immune responses in vertebrate. Cetaceans must have faced new challenges of pathogens when their terrestrial relatives transitioned from the terrestrial to aquatic environment. Here, we sequenced the extracellular domain (ECD) of 10 TLRs in cetacean lineages because this region involved in the recognition of pathogens. A total of 148 sites ranging between 5-26 codons (0.01%-4.83%) were identified to be robust candidates of positive selection at the ECD of 10 TLRs. In addition, the majority (90.54%) of these positively selected codons were found to have radical amino acid changes, which strengthen the evidence of positive selection. Importantly, more radical amino acid changes in selected sites were enriched in the period of early evolutionary transition from land to semi-aquatic and from semi-aquatic to full-aquatic habitat, which might endow cetaceans with a faster adaptation to new pathogens as they transitioned into novel habitat. Interestingly, similar selective intensity was detected in both viral and non-viral TLRs in cetaceans, which was not in line with previous studies on primates and birds that reported stronger positive selection in non-viral TLRs than in viral TLRs. This result may be explained by the fact that cetaceans might have faced diversity of bacteria and viruses during its transitions from terrestrial to aquatic environment whereas both primates and birds probably being affected by only a restricted number of related viruses due to their homogeneous habitat.

Incomplete lineage sorting impacts the inference of macroevolutionary regimes from molecular phylogenies when concatenation is employed: An analysis based on Cetacea

Interest in methods that estimate speciation and extinction rates from molecular phylogenies has increased over the last decade. The application of such methods requires reliable estimates of tree topology and node ages, which are frequently obtained using standard phylogenetic inference combining concatenated loci and molecular dating. However, this practice disregards population-level processes that generate gene tree/species tree discordance. We evaluated the impact of employing concatenation and coalescent-based phylogeny inference in recovering the correct macroevolutionary regime using simulated data based on the well-established diversification rate shift of delphinids in Cetacea. We found that under scenarios of strong incomplete lineage sorting, macroevolutionary analysis of phylogenies inferred by concatenating loci failed to recover the delphinid diversification shift, while the coalescent-based tree consistently retrieved the correct rate regime. We suggest that ignoring microevolutionary processes reduces the power of methods that estimate macroevolutionary regimes from molecular data.

Back to Water: Signature of Adaptive Evolution in Cetacean Mitochondrial tRNAs

The mitochondrion is the power plant of the eukaryotic cell, and tRNAs are the fundamental components of its translational machinery. In the present paper, the evolution of mitochondrial tRNAs was investigated in the Cetacea, a clade of Cetartiodactyla that retuned to water and thus had to adapt its metabolism to a different medium than that of its mainland ancestors. Our analysis focussed on identifying the factors that influenced the evolution of Cetacea tRNA double-helix elements, which play a pivotal role in the formation of the secondary and tertiary structures of each tRNA and consequently manipulate the whole translation machinery of the mitochondrion. Our analyses showed that the substitution pathways in the stems of different tRNAs were influenced by various factors, determining a molecular evolution that was unique to each of the 22 tRNAs. Our data suggested that the composition, AT-skew, and GC-skew of the tRNA stems were the main factors influencing the substitution process. In particular, the range of variation and the fluctuation of these parameters affected the fate of single tRNAs. Strong heterogeneity was observed among the different species of Cetacea. Finally, it appears that the evolution of mitochondrial tRNAs was also shaped by the environments in which the Cetacean taxa differentiated. This latter effect was particularly evident in toothed whales that either live in freshwater or are deep divers.

Comparative genomics reveals conservation of filaggrin and loss of caspase-14 in dolphins

The expression of filaggrin and its stepwise proteolytic degradation are critical events in the terminal differentiation of epidermal keratinocytes and in the formation of the skin barrier to the environment. Here, we investigated whether the evolutionary transition from a terrestrial to a fully aquatic lifestyle of cetaceans, that is dolphins and whales, has been associated with changes in genes encoding filaggrin and proteins involved in the processing of filaggrin. We used comparative genomics, PCRs and re-sequencing of gene segments to screen for the presence and integrity of genes coding for filaggrin and proteases implicated in the maturation of (pro)filaggrin. Filaggrin has been conserved in dolphins (bottlenose dolphin, orca and baiji) but has been lost in whales (sperm whale and minke whale). All other S100 fused-type genes have been lost in cetaceans. Among filaggrin-processing proteases, aspartic peptidase retroviral-like 1 (ASPRV1), also known as saspase, has been conserved, whereas caspase-14 has been lost in all cetaceans investigated. In conclusion, our results suggest that filaggrin is dispensable for the acquisition of fully aquatic lifestyles of whales, whereas it appears to confer an evolutionary advantage to dolphins. The discordant evolution of filaggrin, saspase and caspase-14 in cetaceans indicates that the biological roles of these proteins are not strictly interdependent.

Exploring Massive Incomplete Lineage Sorting in Arctoids (Laurasiatheria, Carnivora)

Freed from the competition of large raptors, Paleocene carnivores could expand their newly acquired habitats in search of prey. Such changing conditions might have led to their successful distribution and rapid radiation. Today, molecular evolutionary biologists are faced, however, with the consequences of such accelerated adaptive radiations, because they led to sequential speciation more rapidly than phylogenetic markers could be fixed. The repercussions being that current genealogies based on such markers are incongruent with species trees.Our aim was to explore such conflicting phylogenetic zones of evolution during the early arctoid radiation, especially to distinguish diagnostic from misleading phylogenetic signals, and to examine other carnivore-related speciation events. We applied a combination of high-throughput computational strategies to screen carnivore and related genomes in silico for randomly inserted retroposed elements that we then used to identify inconsistent phylogenetic patterns in the Arctoidea group, which is well known for phylogenetic discordances.Our combined retrophylogenomic and in vitro wet lab approach detected hundreds of carnivore-specific insertions, many of them confirming well-established splits or identifying and solving conflicting species distributions. Our systematic genome-wide screens for Long INterspersed Elements detected homoplasy-free markers with insertion-specific truncation points that we used to distinguish phylogenetically informative markers from conflicting signals. The results were independently confirmed by phylogenetic diagnostic Short INterspersed Elements. As statistical analysis ruled out ancestral hybridization, these doubly verified but still conflicting patterns were statistically determined to be genomic remnants from a time of ancestral incomplete lineage sorting that especially accompanied large parts of Arctoidea evolution.

The population history of endogenous retroviruses in mule deer (Odocoileus hemionus)

Mobile elements are powerful agents of genomic evolution and can be exceptionally informative markers for investigating species and population-level evolutionary history. While several studies have utilized retrotransposon-based insertional polymorphisms to resolve phylogenies, few population studies exist outside of humans. Endogenous retroviruses are LTR-retrotransposons derived from retroviruses that have become stably integrated in the host genome during past infections and transmitted vertically to subsequent generations. They offer valuable insight into host-virus co-evolution and a unique perspective on host evolutionary history because they integrate into the genome at a discrete point in time. We examined the evolutionary history of a cervid endogenous gammaretrovirus (CrERVγ) in mule deer (Odocoileus hemionus). We sequenced 14 CrERV proviruses (CrERV-in1 to -in14), and examined the prevalence and distribution of 13 proviruses in 262 deer among 15 populations from Montana, Wyoming, and Utah. CrERV absence in white-tailed deer (O. virginianus), identical 5' and 3' long terminal repeat (LTR) sequences, insertional polymorphism, and CrERV divergence time estimates indicated that most endogenization events occurred within the last 200000 years. Population structure inferred from CrERVs (F ST = 0.008) and microsatellites (θ = 0.01) was low, but significant, with Utah, northwestern Montana, and a Helena herd being particularly differentiated. Clustering analyses indicated regional structuring, and non-contiguous clustering could often be explained by known translocations. Cluster ensemble results indicated spatial localization of viruses, specifically in deer from northeastern and western Montana. This study demonstrates the utility of endogenous retroviruses to elucidate and provide novel insight into both ERV evolutionary history and the history of contemporary host populations.

Characterization of hairless (Hr) and FGF5 genes provides insights into the molecular basis of hair loss in cetaceans

Background

Hair is one of the main distinguishing characteristics of mammals and it has many important biological functions. Cetaceans originated from terrestrial mammals and they have evolved a series of adaptations to aquatic environments, which are of evolutionary significance. However, the molecular mechanisms underlying their aquatic adaptations have not been well explored. This study provided insights into the evolution of hair loss during the transition from land to water by investigating and comparing two essential regulators of hair follicle development and hair follicle cycling, i.e., the Hairless (Hr) and FGF5 genes, in representative cetaceans and their terrestrial relatives.

Results

The full open reading frame sequences of the Hr and FGF5 genes were characterized in seven cetaceans. The sequence characteristics and evolutionary analyses suggested the functional loss of the Hr gene in cetaceans, which supports the loss of hair during their full adaptation to aquatic habitats. By contrast, positive selection for the FGF5 gene was found in cetaceans where a series of positively selected amino acid residues were identified.

Conclusions

This is the first study to investigate the molecular basis of the hair loss in cetaceans. Our investigation of Hr and FGF5, two indispensable regulators of the hair cycle, provide some new insights into the molecular basis of hair loss in cetaceans. The results suggest that positive selection for the FGF5 gene might have promoted the termination of hair growth and early entry into the catagen stage of hair follicle cycling. Consequently, the hair follicle cycle was disrupted and the hair was lost completely due to the loss of the Hr gene function in cetaceans. This suggests that cetaceans have evolved an effective and complex mechanism for hair loss.

A phylogenetic blueprint for a modern whale

The emergence of Cetacea in the Paleogene represents one of the most profound macroevolutionary transitions within Mammalia. The move from a terrestrial habitat to a committed aquatic lifestyle engendered wholesale changes in anatomy, physiology, and behavior. The results of this remarkable transformation are extant whales that include the largest, biggest brained, fastest swimming, loudest, deepest diving mammals, some of which can detect prey with a sophisticated echolocation system (Odontoceti - toothed whales), and others that batch feed using racks of baleen (Mysticeti - baleen whales). A broad-scale reconstruction of the evolutionary remodeling that culminated in extant cetaceans has not yet been based on integration of genomic and paleontological information. Here, we first place Cetacea relative to extant mammalian diversity, and assess the distribution of support among molecular datasets for relationships within Artiodactyla (even-toed ungulates, including Cetacea). We then merge trees derived from three large concatenations of molecular and fossil data to yield a composite hypothesis that encompasses many critical events in the evolutionary history of Cetacea. By combining diverse evidence, we infer a phylogenetic blueprint that outlines the stepwise evolutionary development of modern whales. This hypothesis represents a starting point for more detailed, comprehensive phylogenetic reconstructions in the future, and also highlights the synergistic interaction between modern (genomic) and traditional (morphological+paleontological) approaches that ultimately must be exploited to provide a rich understanding of evolutionary history across the entire tree of Life.