Ancient and Recent Hybridization in the Oreochromis Cichlid Fishes.

Cichlid fishes of the genus Oreochromis (tilapia) are among the most important fish for inland capture fisheries and global aquaculture. Deliberate introductions of non-native species for fisheries improvement and accidental escapees from farms have resulted in admixture with indigenous species. Such hybridization may be detrimental to native biodiversity, potentially leading to genomic homogenization of populations and the loss of important genetic material associated with local adaptation. By contrast, introgression may fuel diversification when combined with ecological opportunity, by supplying novel genetic combinations. To date, the role of introgression in the evolutionary history of tilapia has not been explored. Here we studied both ancient and recent hybridization in tilapia, using whole genome resequencing of 575 individuals from 23 species. We focused on Tanzania, a natural hotspot of tilapia diversity, and a country where hybridization between exotic and native species in the natural environment has been previously reported. We reconstruct the first genome-scale phylogeny of the genus and reveal prevalent ancient gene flow across the Oreochromis phylogeny. This has likely resulted in the hybrid speciation of one species, O. chungruruensis. We identify multiple cases of recent hybridization between native and introduced species in the wild, linked to the use of non-native species in both capture fisheries improvement and aquaculture. This has potential implications for both conservation of wild populations and the development of the global tilapia aquaculture industry.

EBP-Colombia and the bioeconomy: Genomics in the service of biodiversity conservation and sustainable development.

The 2016 Peace Agreement has increased access to Colombia's unique ecosystems, which remain understudied and increasingly under threat. The Colombian government has recently announced its National Bioeconomic Strategy (NBS), founded on the sustainable characterization, management, and conservation of the nation's biodiversity as a means to achieve sustainability and peace. Molecular tools will accelerate such endeavors, but capacity remains limited in Colombia. The Earth Biogenome Project's (EBP) objective is to characterize the genomes of all eukaryotic life on Earth through networks of partner institutions focused on sequencing either specific taxa or eukaryotic communities at regional or national scales. Colombia's immense biodiversity and emerging network of stakeholders have inspired the creation of the national partnership "EBP-Colombia." Here, we discuss how this Colombian-driven collaboration between government, academia, and the private sector is integrating research with sustainable, environmentally focused strategies to develop Colombia's postconflict bioeconomy and conserve biological and cultural diversity. EBP-Colombia will accelerate the uptake of technology and promote partnership and exchange of knowledge among Colombian stakeholders and the EBP's global network of experts; assist with conservation strategies to preserve Colombia's vast biological wealth; and promote innovative approaches among public and private institutions in sectors such as agriculture, tourism, recycling, and medicine. EBP-Colombia can thus support Colombia's NBS with the objective of sustainable and inclusive development to address the many social, environmental, and economic challenges, including conflict, inequality, poverty, and low agricultural productivity, and so, offer an alternative model for economic development that similarly placed countries can adopt.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Standards recommendations for the Earth BioGenome Project.

A global international initiative, such as the Earth BioGenome Project (EBP), requires both agreement and coordination on standards to ensure that the collective effort generates rapid progress toward its goals. To this end, the EBP initiated five technical standards committees comprising volunteer members from the global genomics scientific community: Sample Collection and Processing, Sequencing and Assembly, Annotation, Analysis, and IT and Informatics. The current versions of the resulting standards documents are available on the EBP website, with the recognition that opportunities, technologies, and challenges may improve or change in the future, requiring flexibility for the EBP to meet its goals. Here, we describe some highlights from the proposed standards, and areas where additional challenges will need to be met.

Characterizing subgenome recombination and chromosomal imbalances in banana varietal lineages.

BACKGROUND: Bananas and plantains (Musa spp.) are among the most important crops worldwide. The cultivated varieties are vegetatively propagated, so their genetic diversity is essentially fixed over time. Musa acuminata, M. balbisiana and M. schizocarpa have provided the named A, B and S subgenomes that predominantly constitute these varieties. Here we aimed to characterize intergenetic recombination and chromosomal imbalances between these A/B/S subgenomes, which often result in copy-number variants (CNVs) leading to changes in gene dosage and phenotype, in a diverse panel of bananas and plantains. This will allow us to characterize varietal lineages better and identify sources of genetic variation. METHODS: We delimited population structure and clonal lineages in a diverse panel of 188 banana and plantain accessions from the most common cultivars using admixture, principal component and phylogenetic analyses. We used new scalable alignment-based methods, Relative Averaged Alignment (RAA) and Relative Coverage, to infer subgenome composition (AA, AAB, etc.) and interspecific recombination. RESULTS: In our panel, we identified ten varietal lineages composed of somatic clones, plus three groups of tetraploid accessions. We identified chromosomal exchanges resulting in gains/losses in chromosomal segments (CNVs), particularly in AAB and ABB varieties. CONCLUSIONS: We demonstrated alignment-based RAA and Relative Coverage can identify subgenome composition and introgressions with similar results to more complex approaches based on single nucleotide polymorphism (SNP) databases. These ab initio species-agnostic methods can be used without sequencing a panel of wild ancestors to find private SNPs, or in recently diverged pools where private SNPs are uncommon. The extensive A/B/S exchanges and the variation in the length of some introgressions between lineages further support multiple foundational events of hybridization and residual backcrossing. Imbalances between A/B/S may have resulted in CNVs and gene dosage variation. Since most edible banana genomes are fixed on time, these CNVs are stable genetic variations probably associated with phenotypic variation for future genetic studies.

Bacteraemia associated with multiple septic localizations caused by Klebsiella pneumoniae sequence type ST660.

We report a case of Klebsiella pneumoniae bacteraemia in an 80-year-old man in France with no history of travel to Asia, complicated by endogenous endophthalmitis, multiple cerebral microbleeds and hepatic microabscesses, associated with a Bentall endocarditis. Hypervirulence pathotype was suggested based on clinical picture, bacterial isolate genomic sequence and hypermucoidy. Interestingly, the isolate had the non-K1/K2-capsular serotype locus KL113-like, carried a KpVP-1-like virulence plasmid, and belonged to the emerging sublineage SL660 (comprising the sequence type ST660).

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia.

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

Evolution of miRNA-Binding Sites and Regulatory Networks in Cichlids.

The divergence of regulatory regions and gene regulatory network (GRN) rewiring is a key driver of cichlid phenotypic diversity. However, the contribution of miRNA-binding site turnover has yet to be linked to GRN evolution across cichlids. Here, we extend our previous studies by analyzing the selective constraints driving evolution of miRNA and transcription factor (TF)-binding sites of target genes, to infer instances of cichlid GRN rewiring associated with regulatory binding site turnover. Comparative analyses identified increased species-specific networks that are functionally associated to traits of cichlid phenotypic diversity. The evolutionary rewiring is associated with differential models of miRNA- and TF-binding site turnover, driven by a high proportion of fast-evolving polymorphic sites in adaptive trait genes compared with subsets of random genes. Positive selection acting upon discrete mutations in these regulatory regions is likely to be an important mechanism in rewiring GRNs in rapidly radiating cichlids. Regulatory variants of functionally associated miRNA- and TF-binding sites of visual opsin genes differentially segregate according to phylogeny and ecology of Lake Malawi species, identifying both rewired, for example, clade-specific and conserved network motifs of adaptive trait associated GRNs. Our approach revealed several novel candidate regulators, regulatory regions, and three-node motifs across cichlid genomes with previously reported associations to known adaptive evolutionary traits.

In vitro and in silico parameters for precise cgMLST typing of Listeria monocytogenes.

BACKGROUND: Whole genome sequencing analyzed by core genome multi-locus sequence typing (cgMLST) is widely used in surveillance of the pathogenic bacteria Listeria monocytogenes. Given the heterogeneity of available bioinformatics tools to define cgMLST alleles, our aim was to identify parameters influencing the precision of cgMLST profiles. METHODS: We used three L. monocytogenes reference genomes from different phylogenetic lineages and assessed the impact of in vitro (i.e. tested genomes, successive platings, replicates of DNA extraction and sequencing) and in silico parameters (i.e. targeted depth of coverage, depth of coverage, breadth of coverage, assembly metrics, cgMLST workflows, cgMLST completeness) on cgMLST precision made of 1748 core loci. Six cgMLST workflows were tested, comprising assembly-based (BIGSdb, INNUENDO, GENPAT, SeqSphere and BioNumerics) and assembly-free (i.e. kmer-based MentaLiST) allele callers. Principal component analyses and generalized linear models were used to identify the most impactful parameters on cgMLST precision. RESULTS: The isolate's genetic background, cgMLST workflows, cgMLST completeness, as well as depth and breadth of coverage were the parameters that impacted most on cgMLST precision (i.e. identical alleles against reference circular genomes). All workflows performed well at ≥40X of depth of coverage, with high loci detection (> 99.54% for all, except for BioNumerics with 97.78%) and showed consistent cluster definitions using the reference cut-off of ≤7 allele differences. CONCLUSIONS: This highlights that bioinformatics workflows dedicated to cgMLST allele calling are largely robust when paired-end reads are of high quality and when the sequencing depth is ≥40X.

A quantitative framework for characterizing the evolutionary history of mammalian gene expression.

The evolutionary history of a gene helps predict its function and relationship to phenotypic traits. Although sequence conservation is commonly used to decipher gene function and assess medical relevance, methods for functional inference from comparative expression data are lacking. Here, we use RNA-seq across seven tissues from 17 mammalian species to show that expression evolution across mammals is accurately modeled by the Ornstein-Uhlenbeck process, a commonly proposed model of continuous trait evolution. We apply this model to identify expression pathways under neutral, stabilizing, and directional selection. We further demonstrate novel applications of this model to quantify the extent of stabilizing selection on a gene's expression, parameterize the distribution of each gene's optimal expression level, and detect deleterious expression levels in expression data from individual patients. Our work provides a statistical framework for interpreting expression data across species and in disease.

Extensive genome introgression between domestic ferret and European polecat during population recovery in Great Britain.

The European polecat (Mustela putorius) is a mammalian predator which occurs across much of Europe east to the Ural Mountains. In Great Britain, following years of persecution the range of the European polecat contracted and by the early 1900s was restricted to unmanaged forests of central Wales. The European polecat has recently undergone a population increase due to legal protection and its range now overlaps that of feral domestic ferrets (Mustela putorius furo). During this range expansion, European polecats hybridized with feral domestic ferrets producing viable offspring. Here, we carry out population-level whole-genome sequencing on 8 domestic ferrets, 19 British European polecats, and 15 European polecats from the European mainland. We used a range of population genomics methods to examine the data, including phylogenetics, phylogenetic graphs, model-based clustering, phylogenetic invariants, ABBA-BABA tests, topology weighting, and Fst. We found high degrees of genome introgression in British polecats outside their previous stronghold, even in those individuals phenotyped as "pure" polecats. These polecats ranged from presumed F1 hybrids (gamma = 0.53) to individuals that were much less introgressed (gamma = 0.2). We quantify this introgression and find introgressed genes containing Fst outliers associated with cognitive function and sight.

Genomic elements located in the accessory repertoire drive the adaptation to biocides in Listeria monocytogenes strains from different ecological niches.

In response to the massive use of biocides for controlling Listeria monocytogenes (hereafter Lm) contaminations along the food chain, strains showing biocide tolerance emerged. Here, accessory genomic elements were associated with biocide tolerance through pangenome-wide associations performed on 197 Lm strains from different lineages, ecological, geographical and temporal origins. Mobile elements, including prophage-related loci, the Tn6188_qacH transposon and pLMST6_emrC plasmid, were widespread across lineage I and II food strains and associated with tolerance to benzalkonium-chloride (BC), a quaternary ammonium compound (QAC) widely used in food processing. The pLMST6_emrC was also associated with tolerance to another QAC, the didecyldimethylammonium-chloride, displaying a pleiotropic effect. While no associations were detected for chemically reactive biocides (alcohols and chlorines), genes encoding for cell-surface proteins were associated with BC or polymeric biguanide tolerance. The latter was restricted to lineage I strains from animal and the environment. In conclusion, different genetic markers, with polygenic nature or not, appear to have driven the Lm adaptation to biocide, especially in food strains but also from animal and the environment. These markers could aid to monitor and predict the spread of biocide tolerant Lm genotypes across different ecological niches, finally reducing the risk of such strains in food industrial settings.

ROS-mediated PI3K activation drives mitochondrial transfer from stromal cells to hematopoietic stem cells in response to infection.

Hematopoietic stem cells (HSCs) undergo rapid expansion in response to stress stimuli. Here we investigate the bioenergetic processes which facilitate the HSC expansion in response to infection. We find that infection by Gram-negative bacteria drives an increase in mitochondrial mass in mammalian HSCs, which results in a metabolic transition from glycolysis toward oxidative phosphorylation. The initial increase in mitochondrial mass occurs as a result of mitochondrial transfer from the bone marrow stromal cells (BMSCs) to HSCs through a reactive oxygen species (ROS)-dependent mechanism. Mechanistically, ROS-induced oxidative stress regulates the opening of connexin channels in a system mediated by phosphoinositide 3-kinase (PI3K) activation, which allows the mitochondria to transfer from BMSCs into HSCs. Moreover, mitochondria transfer from BMSCs into HSCs, in the response to bacterial infection, occurs before the HSCs activate their own transcriptional program for mitochondrial biogenesis. Our discovery demonstrates that mitochondrial transfer from the bone marrow microenvironment to HSCs is an early physiologic event in the mammalian response to acute bacterial infection and results in bioenergetic changes which underpin emergency granulopoiesis.

Whole genome resequencing data enables a targeted SNP panel for conservation and aquaculture of Oreochromis cichlid fishes.

Cichlid fish of the genus Oreochromis form the basis of the global tilapia aquaculture and fisheries industries. Broodstocks for aquaculture are often collected from wild populations, which in Africa may be from locations containing multiple Oreochromis species. However, many species are difficult to distinguish morphologically, hampering efforts to maintain good quality farmed strains. Additionally, non-native farmed tilapia populations are known to be widely distributed across Africa and to hybridize with native Oreochromis species, which themselves are important for capture fisheries. The morphological identification of these hybrids is particularly unreliable. Here, we describe the development of a single nucleotide polymorphism (SNP) genotyping panel from whole-genome resequencing data that enables targeted species identification in Tanzania. We demonstrate that an optimized panel of 96 genome-wide SNPs based on FST outliers performs comparably to whole genome resequencing in distinguishing species and identifying hybrids. We also show this panel outperforms microsatellite-based and phenotype-based classification methods. Case studies indicate several locations where introduced aquaculture species have become established in the wild, threatening native Oreochromis species. The novel SNP markers identified here represent an important resource for assessing broodstock purity in hatcheries and helping to conserve unique endemic biodiversity.

Newly discovered cichlid fish biodiversity threatened by hybridization with non-native species.

Invasive freshwater fishes are known to readily hybridize with indigenous congeneric species, driving loss of unique and irreplaceable genetic resources. Here we reveal that newly discovered (2013-2016) evolutionarily significant populations of Korogwe tilapia (Oreochromis korogwe) from southern Tanzania are threatened by hybridization with the larger invasive Nile tilapia (Oreochromis niloticus). We use a combination of morphology, microsatellite allele frequencies and whole genome sequences to show that O. korogwe from southern lakes (Nambawala, Rutamba and Mitupa) are distinct from geographically disjunct populations in northern Tanzania (Zigi River and Mlingano Dam). We also provide genetic evidence of O. korogwe × niloticus hybrids in three southern lakes and demonstrate heterogeneity in the extent of admixture across the genome. Finally, using the least admixed genomic regions we estimate that the northern and southern O. korogwe populations most plausibly diverged ~140,000 years ago, suggesting that the geographical separation of the northern and southern groups is not a result of a recent translocation, and instead these populations represent independent evolutionarily significant units. We conclude that these newly discovered and phenotypically unique cichlid populations are already threatened by hybridization with an invasive species, and propose that these irreplaceable genetic resources would benefit from conservation interventions.

Acute myeloid leukemia induces protumoral p16INK4a-driven senescence in the bone marrow microenvironment.

Acute myeloid leukemia (AML) is an age-related disease that is highly dependent on the bone marrow (BM) microenvironment. With increasing age, tissues accumulate senescent cells, characterized by an irreversible arrest of cell proliferation and the secretion of a set of proinflammatory cytokines, chemokines, and growth factors, collectively known as the senescence-associated secretory phenotype (SASP). Here, we report that AML blasts induce a senescent phenotype in the stromal cells within the BM microenvironment and that the BM stromal cell senescence is driven by p16INK4a expression. The p16INK4a-expressing senescent stromal cells then feed back to promote AML blast survival and proliferation via the SASP. Importantly, selective elimination of p16INK4a+ senescent BM stromal cells in vivo improved the survival of mice with leukemia. Next, we find that the leukemia-driven senescent tumor microenvironment is caused by AML-induced NOX2-derived superoxide. Finally, using the p16-3MR mouse model, we show that by targeting NOX2 we reduced BM stromal cell senescence and consequently reduced AML proliferation. Together, these data identify leukemia-generated NOX2-derived superoxide as a driver of protumoral p16INK4a-dependent senescence in BM stromal cells. Our findings reveal the importance of a senescent microenvironment for the pathophysiology of leukemia. These data now open the door to investigate drugs that specifically target the "benign" senescent cells that surround and support AML.

A new genome allows the identification of genes associated with natural variation in aluminium tolerance in Brachiaria grasses.

Toxic concentrations of aluminium cations and low phosphorus availability are the main yield-limiting factors in acidic soils, which represent half of the potentially available arable land. Brachiaria grasses, which are commonly sown as forage in the tropics because of their resilience and low demand for nutrients, show greater tolerance to high concentrations of aluminium cations (Al3+) than most other grass crops. In this work, we explored the natural variation in tolerance to Al3+ between high and low tolerant Brachiaria species and characterized their transcriptional differences during stress. We identified three QTLs (quantitative trait loci) associated with root vigour during Al3+ stress in their hybrid progeny. By integrating these results with a new Brachiaria reference genome, we identified 30 genes putatively responsible for Al3+ tolerance in Brachiaria. We observed differential expression during stress of genes involved in RNA translation, response signalling, cell wall composition, and vesicle location homologous to aluminium-induced proteins involved in limiting uptake or localizing the toxin. However, there was limited regulation of malate transporters in Brachiaria, which suggests that exudation of organic acids and other external tolerance mechanisms, common in other grasses, might not be relevant in Brachiaria. The contrasting regulation of RNA translation and response signalling suggests that response timing is critical in high Al3+-tolerant Brachiaria.

A comprehensive genomic history of extinct and living elephants.

Elephantids are the world's most iconic megafaunal family, yet there is no comprehensive genomic assessment of their relationships. We report a total of 14 genomes, including 2 from the American mastodon, which is an extinct elephantid relative, and 12 spanning all three extant and three extinct elephantid species including an ∼120,000-y-old straight-tusked elephant, a Columbian mammoth, and woolly mammoths. Earlier genetic studies modeled elephantid evolution via simple bifurcating trees, but here we show that interspecies hybridization has been a recurrent feature of elephantid evolution. We found that the genetic makeup of the straight-tusked elephant, previously placed as a sister group to African forest elephants based on lower coverage data, in fact comprises three major components. Most of the straight-tusked elephant's ancestry derives from a lineage related to the ancestor of African elephants while its remaining ancestry consists of a large contribution from a lineage related to forest elephants and another related to mammoths. Columbian and woolly mammoths also showed evidence of interbreeding, likely following a latitudinal cline across North America. While hybridization events have shaped elephantid history in profound ways, isolation also appears to have played an important role. Our data reveal nearly complete isolation between the ancestors of the African forest and savanna elephants for ∼500,000 y, providing compelling justification for the conservation of forest and savanna elephants as separate species.

Analysis of structural variants in four African cichlids highlights an association with developmental and immune related genes.

BACKGROUND: East African lake cichlids are one of the most impressive examples of an adaptive radiation. Independently in Lake Victoria, Tanganyika, and Malawi, several hundreds of species arose within the last 10 million to 100,000 years. Whereas most analyses in cichlids focused on nucleotide substitutions across species to investigate the genetic bases of this explosive radiation, to date, no study has investigated the contribution of structural variants (SVs) in the evolution of adaptive traits across the three Great Lakes of East Africa. RESULTS: Here, we annotate and characterize the repertoires and evolutionary potential of different SV classes (deletion, duplication, inversion, insertions and translocations) in four cichlid species: Haplochromis burtoni, Metriaclima zebra, Neolamprologus brichardi and Pundamilia nyererei. We investigate the patterns of gain and loss evolution for each SV type, enabling the identification of lineage specific events. Both deletions and inversions show a significant overlap with SINE elements, while inversions additionally show a limited, but significant association with DNA transposons. Inverted regions are enriched for genes regulating behaviour, or involved in skeletal and visual system development. We also find that duplicated regions show enrichment for genes associated with "antigen processing and presentation" and other immune related categories. Our pipeline and results were further tested by PCR validation of selected deletions and inversions, which confirmed respectively 7 out of 10 and 6 out of 9 events. CONCLUSIONS: Altogether, we provide the first comprehensive overview of rearrangement evolution in East African cichlids, and some important insights into their likely contribution to adaptation.

microRNA profiling in the Weddell seal suggests novel regulatory mechanisms contributing to diving adaptation.

BACKGROUND: The Weddell Seal (Leptonychotes weddelli) represents a remarkable example of adaptation to diving among marine mammals. This species is capable of diving > 900 m deep and remaining underwater for more than 60 min. A number of key physiological specializations have been identified, including the low levels of aerobic, lipid-based metabolism under hypoxia, significant increase in oxygen storage in blood and muscle; high blood volume and extreme cardiovascular control. These adaptations have been linked to increased abundance of key proteins, suggesting an important, yet still understudied role for gene reprogramming. In this study, we investigate the possibility that post-transcriptional gene regulation by microRNAs (miRNAs) has contributed to the adaptive evolution of diving capacities in the Weddell Seal. RESULTS: Using small RNA data across 4 tissues (brain, heart, muscle and plasma), in 3 biological replicates, we generate the first miRNA annotation in this species, consisting of 559 high confidence, manually curated miRNA loci. Evolutionary analyses of miRNA gain and loss highlight a high number of Weddell seal specific miRNAs. Four hundred sixteen miRNAs were differentially expressed (DE) among tissues, whereas 80 miRNAs were differentially expressed (DE) across all tissues between pups and adults and age differences for specific tissues were detected in 188 miRNAs. mRNA targets of these altered miRNAs identify possible protective mechanisms in individual tissues, particularly relevant to hypoxia tolerance, anti-apoptotic pathways, and nitric oxide signal transduction. Novel, lineage-specific miRNAs associated with developmental changes target genes with roles in angiogenesis and vasoregulatory signaling. CONCLUSIONS: Altogether, we provide an overview of miRNA composition and evolution in the Weddell seal, and the first insights into their possible role in the specialization to diving.