Benchmarks in antimicrobial peptide prediction are biased due to the selection of negative data.

Antimicrobial peptides (AMPs) are a heterogeneous group of short polypeptides that target not only microorganisms but also viruses and cancer cells. Due to their lower selection for resistance compared with traditional antibiotics, AMPs have been attracting the ever-growing attention from researchers, including bioinformaticians. Machine learning represents the most cost-effective method for novel AMP discovery and consequently many computational tools for AMP prediction have been recently developed. In this article, we investigate the impact of negative data sampling on model performance and benchmarking. We generated 660 predictive models using 12 machine learning architectures, a single positive data set and 11 negative data sampling methods; the architectures and methods were defined on the basis of published AMP prediction software. Our results clearly indicate that similar training and benchmark data set, i.e. produced by the same or a similar negative data sampling method, positively affect model performance. Consequently, all the benchmark analyses that have been performed for AMP prediction models are significantly biased and, moreover, we do not know which model is the most accurate. To provide researchers with reliable information about the performance of AMP predictors, we also created a web server AMPBenchmark for fair model benchmarking. AMPBenchmark is available at http://BioGenies.info/AMPBenchmark.

Evolutionary responses of a reef-building coral to climate change at the end of the last glacial maximum.

Climate change threatens the survival of coral reefs on a global scale, primarily through mass bleaching and mortality as a result of marine heatwaves. While these short-term effects are clear, predicting the fate of coral reefs over the coming century is a major challenge. One way to understand the longer-term effects of rapid climate change is to examine the response of coral populations to past climate shifts. Coastal and shallow-water marine ecosystems such as coral reefs have been reshaped many times by sea-level changes during the Pleistocene, yet, few studies have directly linked this with its consequences on population demographics, dispersal, and adaptation. Here we use powerful analytical techniques, afforded by haplotype phased whole-genomes, to establish such links for the reef-building coral, Acropora digitifera. We show that three genetically distinct populations are present in northwestern Australia, and that their rapid divergence since the last glacial maximum (LGM) can be explained by a combination of founder-effects and restricted gene flow. Signatures of selective sweeps, too strong to be explained by demographic history, are present in all three populations and overlap with genes that show different patterns of functional enrichment between inshore and offshore habitats. In contrast to rapid divergence in the host, we find that photosymbiont communities are largely undifferentiated between corals from all three locations, spanning almost 1000 km, indicating that selection on host genes and not acquisition of novel symbionts, has been the primary driver of adaptation for this species in northwestern Australia.

SNP data reveals the complex and diverse evolutionary history of the blue-ringed octopus genus (Octopodidae: Hapalochlaena) in the Asia-Pacific.

The blue-ringed octopus species complex (Hapalochlaena spp.), known to occur from Southern Australia to Japan, currently contains four formally described species (Hapalochlaena maculosa, Hapalochlaena fasciata, Hapalochlaena lunulata and Hapalochlaena nierstraszi). These species are distinguished based on morphological characters (iridescent blue rings and/or lines) along with reproductive strategies. However, the observation of greater morphological diversity than previously captured by the current taxonomic framework indicates that a revision is required. To examine species boundaries within the genus we used mitochondrial (12S rRNA, 16S rRNA, cytochrome c oxidase subunit 1 [COI], cytochrome c oxidase subunit 3 [COIII] and cytochrome b [Cytb]) and genome-wide SNP data (DaRT seq) from specimens collected across its geographic range including variations in depth from 3 m to >100 m. This investigation indicates substantially greater species diversity present within the genus Hapalochlaena than is currently described. We identified 10,346 SNPs across all locations, which when analysed support a minimum of 11 distinct clades. Bayesian phylogenetic analysis of the mitochondrial COI gene on a more limited sample set dates the diversification of the genus to âˆ¼30 mya and corroborates eight of the lineages indicated by the SNP analyses. Furthermore, we demonstrate that the diagnostic lined patterning of H. fasciata found in North Pacific waters and NSW, Australia is polyphyletic and therefore likely the result of convergent evolution. Several "deep water" (>100 m) lineages were also identified in this study with genetic convergence likely to be driven by external selective pressures. Examination of morphological traits, currently being undertaken in a parallel morphological study, is required to describe additional species within the complex.

ampir: an R package for fast genome-wide prediction of antimicrobial peptides.

SUMMARY: Antimicrobial peptides (AMPs) are the key components of the innate immune system that protect against pathogens, regulate the microbiome and are promising targets for pharmaceutical research. Computational tools based on machine learning have the potential to aid discovery of genes encoding novel AMPs but existing approaches are not designed for genome-wide scans. To facilitate such genome-wide discovery of AMPs we developed a fast and accurate AMP classification framework, ampir. ampir is designed for high throughput, integrates well with existing bioinformatics pipelines, and has much higher classification accuracy than existing methods when applied to whole genome data. AVAILABILITY AND IMPLEMENTATION: ampir is implemented primarily in R with core feature calculation methods written in C++. Release versions are available via CRAN and work on all major operating systems. The development version is maintained at https://github.com/legana/ampir. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Newly Discovered Peptides from the Coral Heliofungia actiniformis Show Structural and Functional Diversity.

Scleractinian corals are crucially important to the health of some of the world's most biodiverse, productive, and economically important marine habitats. Despite this importance, analysis of coral peptidomes is still in its infancy. Here we show that the tentacle extract from the stony coral Heliofungia actiniformis is rich in peptides with diverse and novel structures. We have characterized the sequences and three-dimensional structures of four new peptides, three of which have no known homologues. We show that a 2 kDa peptide, Hact-2, promotes significant cell proliferation on human cells and speculate this peptide may be involved in the remarkable regenerative capacity of corals. We found a 3 kDa peptide, Hact-3, encoded within a fascin-like domain, and homologues of Hact-3 are present in the genomes of other coral species. Two additional peptides, Hact-4 and Hact-SCRiP1, with limited sequence similarity, both contain a beta-defensin-like fold and highlight a structural link with the small cysteine-rich proteins (SCRiP) family of proteins found predominantly in corals. Our results provide a first glimpse into the remarkable and unexplored structural diversity of coral peptides, providing insight into their diversity and putative functions and, given the ancient lineage of corals, potential insight into the evolution of structural motifs.

Establishing a distributed national research infrastructure providing bioinformatics support to life science researchers in Australia.

EMBL Australia Bioinformatics Resource (EMBL-ABR) is a developing national research infrastructure, providing bioinformatics resources and support to life science and biomedical researchers in Australia. EMBL-ABR comprises 10 geographically distributed national nodes with one coordinating hub, with current funding provided through Bioplatforms Australia and the University of Melbourne for its initial 2-year development phase. The EMBL-ABR mission is to: (1) increase Australia's capacity in bioinformatics and data sciences; (2) contribute to the development of training in bioinformatics skills; (3) showcase Australian data sets at an international level and (4) enable engagement in international programs. The activities of EMBL-ABR are focussed in six key areas, aligning with comparable international initiatives such as ELIXIR, CyVerse and NIH Commons. These key areas-Tools, Data, Standards, Platforms, Compute and Training-are described in this article.

Biogeography, reproductive biology and phylogenetic divergence within the Fungiidae (mushroom corals).

While the escalating impacts of climate change and other anthropogenic pressures on coral reefs are well documented at the coral community level, studies of species-specific trends are less common, owing mostly to the difficulties and uncertainties in delineating coral species. It has also become clear that traditional coral taxonomy based largely on skeletal macromorphology has underestimated the diversity of many coral families. Here, we use targeted enrichment methods to sequence 2476 ultraconserved elements (UCEs) and exonic loci to investigate the relationship between populations of Fungia fungites from Okinawa, Japan, where this species reproduces by brooding (i.e., internal fertilization), and Papua New Guinea and Australia, where it reproduces by broadcast-spawning (i.e., external fertilization). Moreover, we analyzed the relationships between populations of additional fungiid species (Herpolitha limax and Ctenactis spp.) that reproduce only by broadcast-spawning. Our phylogenetic and species delimitation analyses reveal strong biogeographic structuring in both F. fungites and Herpolitha limax, consistent with cryptic speciation in Okinawa in both species and additionally for H. limax in the Red Sea. By combining UCE/exon data and mitochondrial sequences captured in off-target reads, we reinforce earlier findings that Ctenactis, a genus consisting of three nominal morphospecies, is not a natural group. Our results highlight the need for taxonomic and systematic re-evaluations of some species and genera within the family Fungiidae. This work demonstrates that sequence data generated by the application of targeted capture methods can provide objective criteria by which we can test phylogenetic hypotheses based on morphological and/or life history traits.

Coupled Genomic Evolutionary Histories as Signatures of Organismal Innovations in Cephalopods: Co-evolutionary Signatures Across Levels of Genome Organization May Shed Light on Functional Linkage and Origin of Cephalopod Novelties.

How genomic innovation translates into organismal organization remains largely unanswered. Possessing the largest invertebrate nervous system, in conjunction with many species-specific organs, coleoid cephalopods (octopuses, squids, cuttlefishes) provide exciting model systems to investigate how organismal novelties evolve. However, dissecting these processes requires novel approaches that enable deeper interrogation of genome evolution. Here, the existence of specific sets of genomic co-evolutionary signatures between expanded gene families, genome reorganization, and novel genes is posited. It is reasoned that their co-evolution has contributed to the complex organization of cephalopod nervous systems and the emergence of ecologically unique organs. In the course of reviewing this field, how the first cephalopod genomic studies have begun to shed light on the molecular underpinnings of morphological novelty is illustrated and their impact on directing future research is described. It is argued that the application and evolutionary profiling of evolutionary signatures from these studies will help identify and dissect the organismal principles of cephalopod innovations. By providing specific examples, the implications of this approach both within and beyond cephalopod biology are discussed.

Quantitative Proteomic Analysis of the Slime and Ventral Mantle Glands of the Striped Pyjama Squid (Sepioloidea lineolata).

Cephalopods are known to produce an extensive range of secretions including ink, mucus, and venom. Sepiadariidae, a family of small, benthic bobtail squids, are notable for the high volume of viscous slime they emit when stressed. One species, Sepioloidea lineolata (striped pyjama squid), is covered with glands along the perimeter of the ventral mantle, and these structures are hypothesized to be the source of its slime. Using label-free quantitative proteomics, we analyzed five tissue types (dorsal and ventral mantle muscle, dorsal and ventral epithelium, and ventral mantle glands) and the slime from four individuals. In doing so, we were able to determine the relationship between the slime and the tissues as well as highlight proteins that were specifically identified within the slime and ventral mantle glands. A total of 28 proteins were identified to be highly enriched in slime, and these were composed of peptidases and protease inhibitors. Seven of these proteins contained predicted signal peptides, indicating classical secretion, with four proteins having no identifiable domains or similarity to any known proteins. The ventral mantle glands also appear to be the tissue with the closest overall proteomic composition to the slime; therefore, it is likely that the slime originates, at least in part, from these glands.

Identification and Characterization of a Peptide from the Stony Coral Heliofungia actiniformis.

Marine organisms produce a diverse range of toxins and bioactive peptides to support predation, competition, and defense. The peptide repertoires of stony corals (order Scleractinia) remain relatively understudied despite the presence of tentacles used for predation and defense that are likely to contain a range of bioactive compounds. Here, we show that a tentacle extract from the mushroom coral, Heliofungia actiniformis, contains numerous peptides with a range of molecular weights analogous to venom profiles from species such as cone snails. Using NMR spectroscopy and mass spectrometry we characterized a 12-residue peptide (Hact-1) with a new sequence (GCHYTPFGLICF) and well-defined ß-hairpin structure stabilized by a single disulfide bond. The sequence is encoded within the genome of the coral and expressed in the polyp body tissue. The structure present is common among toxins and venom peptides, but Hact-1 does not show activity against select examples of Gram-positive and Gram-negative bacteria or a range of ion channels, common properties of such peptides. Instead, it appears to have a limited effect on human peripheral blood mononuclear cells, but the ecological function of the peptide remains unknown. The discovery of this peptide from H. actiniformis is likely to be the first of many from this and related species.

Comparative Proteomic Analysis of Slime from the Striped Pyjama Squid, Sepioloidea lineolata, and the Southern Bottletail Squid, Sepiadarium austrinum (Cephalopoda: Sepiadariidae).

Sepioloidea lineolata, the striped pyjama squid (family Sepiadariidae), is a small species of benthic bobtail squid distributed along the Southern Indo-Pacific coast of Australia. Like other sepiadariid squids, it is known to secrete large volumes of viscous slime when stressed. In order to identify key proteins involved in the function of sepiadariid slimes, we compared the slime proteome of Sepioloidea lineolata with that of a closely related species, Sepiadarium austrinum. Of the 550 protein groups identified in Sepioloidea lineolata slime, 321 had orthologs in Sepiadarium austrinum, and the abundance of these (iBAQ) was highly correlated between species. Both slimes were dominated by a small number of abundant proteins, and several of these were short secreted proteins with no homologues outside the class Cephalopoda. No mucins were identified within either species' slime, suggesting that it is structurally distinct from mucin polymer-based gels found in many vertebrate and echinoderm secretions. The extent of N-glycosylation in the slime of Sepioloidea lineolata was also studied via glycan cleavage with Peptide: N-glycosidase F (PNGase-F). Although very few (four) proteins showed strong evidence of N-glycosylation, we found that treatment with PNGase-F led to a slight increase in peptide identification rates compared with controls.

Shotgun Proteomics Analysis of Saliva and Salivary Gland Tissue from the Common Octopus Octopus vulgaris.

The salivary apparatus of the common octopus ( Octopus vulgaris) has been the subject of biochemical study for over a century. A combination of bioassays, behavioral studies and molecular analysis on O. vulgaris and related species suggests that its proteome should contain a mixture of highly potent neurotoxins and degradative proteins. However, a lack of genomic and transcriptomic data has meant that the amino acid sequences of these proteins remain almost entirely unknown. To address this, we assembled the posterior salivary gland transcriptome of O. vulgaris and combined it with high resolution mass spectrometry data from the posterior and anterior salivary glands of two adults, the posterior salivary glands of six paralarvae and the saliva from a single adult. We identified a total of 2810 protein groups from across this range of salivary tissues and age classes, including 84 with homology to known venom protein families. Additionally, we found 21 short secreted cysteine rich protein groups of which 12 were specific to cephalopods. By combining protein expression data with phylogenetic analysis we demonstrate that serine proteases expanded dramatically within the cephalopod lineage and that cephalopod specific proteins are strongly associated with the salivary apparatus.

Liver proteome response of pre-harvest Atlantic salmon following exposure to elevated temperature.

BACKGROUND: Atlantic salmon production in Tasmania (Southern Australia) occurs near the upper limits of the species thermal tolerance. Summer water temperatures can average over 19 °C over several weeks and have negative effects on performance and health. Liver tissue exerts important metabolic functions in thermal adaptation. With the aim of identifying mechanisms underlying liver plasticity in response to chronic elevated temperature in Atlantic salmon, label-free shotgun proteomics was used to explore quantitative protein changes after 43 days of exposure to elevated temperature. RESULTS: A total of 276 proteins were differentially (adjusted p-value < 0.05) expressed between the control (15 °C) and elevated (21 °C) temperature treatments. As identified by Ingenuity Pathway Analysis (IPA), transcription and translation mechanisms, protein degradation via the proteasome, and cytoskeletal components were down-regulated at elevated temperature. In contrast, an up-regulated response was identified for NRF2-mediated oxidative stress, endoplasmic reticulum stress, and amino acid degradation. The proteome response was paralleled by reduced fish condition factor and hepato-somatic index at elevated temperature. CONCLUSIONS: The present study provides new evidence of the interplay among different cellular machineries in a scenario of heat-induced energy deficit and oxidative stress, and refines present understanding of how Atlantic salmon cope with chronic exposure to temperature near the upper limits of thermal tolerance.

AlignStat: a web-tool and R package for statistical comparison of alternative multiple sequence alignments.

BACKGROUND: Alternative sequence alignment algorithms yield different results. It is therefore useful to quantify the similarities and differences between alternative alignments of the same sequences. These measurements can identify regions of consensus that are likely to be most informative in downstream analysis. They can also highlight systematic differences between alignments that relate to differences in the alignment algorithms themselves. RESULTS: Here we present a simple method for aligning two alternative multiple sequence alignments to one another and assessing their similarity. Differences are categorised into merges, splits or shifts in one alignment relative to the other. A set of graphical visualisations allow for intuitive interpretation of the data. CONCLUSIONS: AlignStat enables the easy one-off online use of MSA similarity comparisons or into R pipelines. The web-tool is available at AlignStat.Science.LaTrobe.edu.au. The R package, readme and example data are available on CRAN and GitHub.com/TS404/AlignStat.

Combined Transcriptomic and Proteomic Analysis of the Posterior Salivary Gland from the Southern Blue-Ringed Octopus and the Southern Sand Octopus.

This study provides comprehensive proteomic profiles from the venom producing posterior salivary glands of octopus (superorder Octopodiformes) species. A combined transcriptomic and proteomic approach was used to identify 1703 proteins from the posterior salivary gland of the southern blue-ringed octopus, Hapalochlaena maculosa and 1300 proteins from the posterior salivary gland of the southern sand octopus, Octopus kaurna. The two proteomes were broadly similar; clustering of proteins into orthogroups revealed 937 that were shared between species. Serine proteases were particularly diverse and abundant in both species. Other abundant proteins included a large number of secreted proteins, many of which had no known conserved domains, or homology to proteins with known function. On the basis of homology to known venom proteins, 23 putative toxins were identified in H. maculosa and 24 in O. kaurna. These toxins span nine protein families: CAP (cysteine rich secretory proteins, antigen 5, parthenogenesis related), chitinase, carboxylesterase, DNase, hyaluronidase, metalloprotease, phospholipase, serine protease and tachykinin. Serine proteases were responsible for 70.9% and 86.3% of putative toxin expression in H. maculosa and O. kaurna, respectively, as determined using intensity based absolute quantification (iBAQ) measurements. Phylogenetic analysis of the putative toxin serine proteases revealed a similar suite of diverse proteins present in both species. Posterior salivary gland composition of H. maculosa and O. kaurna differ in several key aspects. While O. kaurna expressed the proteinaceous neurotoxin, tachykinin, this was absent from H. maculosa, perhaps reflecting the acquisition of a potent nonproteinaceous neurotoxin, tetrodotoxin (TTX) produced by bacteria in the salivary glands of that species. The dispersal factor, hyaluronidase was particularly abundant in H. maculosa. Chitinase was abundant in both species and is believed to facilitate envenomation in chitinous prey such as crustaceans. Cephalopods represent a largely unexplored source of novel proteins distinct from all other venomous taxa and are of interest for further inquiry, as novel proteinaceous toxins derived from venoms may contribute to pharmaceutical design.

Proteogenomic analysis of the Venturia pirina (Pear Scab Fungus) secretome reveals potential effectors.

A proteogenomic analysis is presented for Venturia pirina, a fungus that causes scab disease on European pear (Pyrus communis). V. pirina is host-specific, and the infection is thought to be mediated by secreted effector proteins. Currently, only 36 V. pirina proteins are catalogued in GenBank, and the genome sequence is not publicly available. To identify putative effectors, V. pirina was grown in vitro on and in cellophane sheets mimicking its growth in infected leaves. Secreted extracts were analyzed by tandem mass spectrometry, and the data (ProteomeXchange identifier PXD000710) was queried against a protein database generated by combining in silico predicted transcripts with six frame translations of a whole genome sequence of V. pirina (GenBank Accession JEMP00000000 ). We identified 1088 distinct V. pirina protein groups (FDR 1%) including 1085 detected for the first time. Thirty novel (not in silico predicted) proteins were found, of which 14 were identified as potential effectors based on characteristic features of fungal effector protein sequences. We also used evidence from semitryptic peptides at the protein N-terminus to corroborate in silico signal peptide predictions for 22 proteins, including several potential effectors. The analysis highlights the utility of proteogenomics in the study of secreted effectors.

Microbiome manipulation by corals and other Cnidaria via quorum quenching.

A dynamic mucous layer containing numerous micro-organisms covers the surface of corals and has multiple functions including both removal of sediment and "food gathering."1 It is likely to also act as the primary barrier to infection; various proteins and compounds with antimicrobial activity have been identified in coral mucus, though these are thought to be largely or exclusively of microbial origin. As in Hydra,2 anti-microbial peptides (AMPs) are likely to play major roles in regulating the microbiomes of corals.3,4 Some eukaryotes employ a complementary but less obvious approach to manipulate their associated microbiome by interfering with quorum signaling, effectively preventing bacteria from coordinating gene expression across a population. Our investigation of immunity in the reef-building coral Acropora millepora,5 however, led to the discovery of a coral gene referred to here as AmNtNH1 that can inactivate a range of acyl homoserine lactones (AHLs), common bacterial quorum signaling molecules, and is induced on immune challenge of adult corals and expressed during the larval settlement process. Closely related proteins are widely distributed within the Scleractinia (hard corals) and some other cnidarians, with multiple paralogs in Acropora, but their closest relatives are bacterial, implying that these are products of one or more lateral gene transfer events post-dating the cnidarian-bilaterian divergence. The deployment by corals of genes used by bacteria to compete with other bacteria reflects a mechanism of microbiome manipulation previously unknown in Metazoa but that may apply more generally.

Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef.

Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.

Genomic evidence for West Antarctic Ice Sheet collapse during the Last Interglacial.

The marine-based West Antarctic Ice Sheet (WAIS) is considered vulnerable to irreversible collapse under future climate trajectories, and its tipping point may lie within the mitigated warming scenarios of 1.5° to 2°C of the United Nations Paris Agreement. Knowledge of ice loss during similarly warm past climates could resolve this uncertainty, including the Last Interglacial when global sea levels were 5 to 10 meters higher than today and global average temperatures were 0.5° to 1.5°C warmer than preindustrial levels. Using a panel of genome-wide, single-nucleotide polymorphisms of a circum-Antarctic octopus, we show persistent, historic signals of gene flow only possible with complete WAIS collapse. Our results provide the first empirical evidence that the tipping point of WAIS loss could be reached even under stringent climate mitigation scenarios.

Spatially varying selection between habitats drives physiological shifts and local adaptation in a broadcast spawning coral on a remote atoll in Western Australia.

At the Rowley Shoals in Western Australia, the prominent reef flat becomes exposed on low tide and the stagnant water in the shallow atoll lagoons heats up, creating a natural laboratory for characterizing the mechanisms of coral resilience to climate change. To explore these mechanisms in the reef coral Acropora tenuis, we collected samples from lagoon and reef slope habitats and combined whole-genome sequencing, ITS2 metabarcoding, experimental heat stress, and transcriptomics. Despite high gene flow across the atoll, we identified clear shifts in allele frequencies between habitats at relatively small linked genomic islands. Common garden heat stress assays showed corals from the lagoon to be more resistant to bleaching, and RNA sequencing revealed marked differences in baseline levels of gene expression between habitats. Our results provide new insight into the complex mechanisms of coral resilience to climate change and highlight the potential for spatially varying selection across complex coral reef seascapes to drive pronounced ecological divergence in climate-related traits.