Skip to main content
Download PDF

Patterns of enrichment and acceleration in evolutionary rates of promoters suggest a role of regulatory regions in cetacean gigantism

BMC Ecology and Evolution volume 23, Article number: 62 (2023) Cite this article

Abstract

Background

Cetaceans (whales, porpoises, and dolphins) are a lineage of aquatic mammals from which some species became giants. Only recently, gigantism has been investigated from the molecular point of view. Studies focused mainly on coding regions, and no data on the influence of regulatory regions on gigantism in this group was available. Accordingly, we investigated the molecular evolution of non-coding regulatory regions of genes already described in the literature for association with size in mammals, focusing mainly on the promoter regions. For this, we used Ciiider and phyloP tools. Ciiider identifies significantly enriched transcription factor binding sites, and phyloP estimates the molecular evolution rate of the promoter.

Results

We found evidence of enrichment of transcription binding factors related to large body size, with distinct patterns between giant and non-giant cetaceans in the IGFBP7 and NCAPG promoters, in which repressive agents are present in small cetaceans and those that stimulate transcription, in giant cetaceans. In addition, we found evidence of acceleration in the IGF2, IGFBP2, IGFBP7, and ZFAT promoters.

Conclusion

Our results indicate that regulatory regions may also influence cetaceans’ body size, providing candidate genes for future research to understand the molecular basis of the largest living animals.

Peer Review reports

Background

Cetaceans (whales, porpoises, and dolphins) compose a lineage exclusively of aquatic mammals, classified into two groups: the odontocetes—animals with teeth—and mysticetes—animals with baleen plates that allow for food filtration [1]. Cetaceans evolved from small-sized terrestrial ancestors approximately 50 Myr ago, during the Eocene [2]. By then, cetaceans started recolonizing the aquatic environment, a process followed by extensive morphological and physiological modifications such as reducing olfactory and gustatory systems, loss of hind limbs, and modifications toward a hydrodynamic body [3]. Some cetacean species have become gigantic, with colossal measures that are not achieved by living animals. Gigantism results from species evolving enormous body sizes compared with their small-sized ancestors. This feature affects critical life-history traits, such as fecundity, due to the consequent lower reproductive rate and an overall reduction in effective population genetic size (Ne) due to lower population densities [4, 5]. Despite this, some cetacean species reach large body sizes that are unique among living animals, ranging from the impressive gray whale (Eschrichtius robustus) with 15 m to the colossal blue whale (Balaenoptera musculus) that reaches up to 30 m [6, 7].

Some ecological hypotheses have been proposed to explain the large body proportions in cetaceans, such as thermoregulation [8], a wider space available in the aquatic environment to explore new niches [9], and also food acquisition, which in mysticetes is associated with filtration of small prey [10], and in sperm whale, the largest odontocete at 20 m in length [11], with the ability to dive to extraordinary depths to capture their prey [12]. In addition to these ecological causes, the genetics behind body size has been recently investigated, taking advantage of the sequenced cetacean genomes. For example, evolutionary analyses have shown signatures of positive selection on size-related genes in cetaceans. Sun et al. 2019 [13] found evidence of selection in genes related to small size in cetaceans, such as ACAN, OBSL1, and GRB10 genes; whereas, in giant cetaceans, genes possibly evolving under positive selection were those with known roles in promoting growth and large sizes, such as CBS, EIF2AK3, and PLOD1 genes. Still, these studies focused only on coding regions, and information on the influence of regulatory regions on gigantism in this group is scarce.

Non-coding sequences with regulatory functions (e.g., promoters and enhancers) coordinate the spatial-temporal expression of genes [14]. Although regulatory regions are not under the same constraints as coding sequences, highly conserved sequence blocks in different species indicate evolutionarily conserved functions [15, 16]. On the other hand, modifications of gene regulatory elements have been associated with phenotypic changes in animal evolution, such as pigmentation changes in dogs [17], bristle patterns in flies [18], and skeletal differences in fish [19]. This approach to studying transcription factors is currently facilitated by computational methods that can identify potential candidate gene regulatory elements by detecting regions of the genome that exhibit evolutionary conservation or acceleration [16].

Comparative genome-wide regulatory sequence approaches can provide insights into the evolutionary history of large body size in cetaceans. Specifically, our study focuses on species that are at least 10 m long and classified as giants. In Fig. 1 we present the cetacean species included in our investigation, highlighted in blue for giants and red for non-giants, based on their average size. Accordingly, we investigated the molecular evolution of non-coding regulatory regions of genes previously described in the literature as being associated with size in mammals, such as EGF, GHSR, IGF2, IGFBP2, IGFBP7, LCORL, NCAPG, PLAG1, and ZFAT, focusing on cetaceans. Our analyses were performed within a phylogenetic framework, where each promoter contained a consistent set of 52 species, including 39 from different orders of mammals and 13 species of cetaceans, of which eight were classified as giants with a minimum length of 10 m. The objective was to investigate differences in the enrichment of Transcription Factor Binding Sites (TFBS) between giant and non-giant cetaceans, as well as to identify potential evolutionary acceleration in these animals with large body sizes.

Results

Phylogenetic reconstruction

To identify potential candidate genes that may contribute to gigantism in cetaceans, we examined the evolution of their regulatory regions. To this end, we generated phylogenetic trees for each gene selected in this study using a maximum likelihood approach. These trees were constructed to visually explore the evolutionary relationships among the promoter sequences of the species included in our dataset, focusing on potential convergence among giant species. Specifically, we generated phylogenetic trees for the promoter region (-1500 bp to + 500 bp from TSS) of each gene. The promoter regions were defined as the sequences upstream of the transcription start site, as this is where regulatory elements, such as transcription factor binding sites, are typically located. The correct grouping of mysticetes and odontocetes was observed in most of the phylogenetic trees, along with other groups of mammals, such as Artiodactyla, Carnivora, Primate, Cingulata, and Chiroptera (Additional file 1: Supplementary Figures S1-S8 show the phylogenetic trees for each promoter).

One exception was observed in the NCAPG promoter, in which the odontocete sperm whale (Physeter catodon) was grouped within the mysticete clade, and the mysticete minke whale (Balaenoptera acutorostrata) was grouped with the odontocetes (Fig. 2). Thus, the NCAPG tree had a clade formed by the gigantic animals included in this dataset, Balaenoptera musculus, Physeter catodon, Eschrichtius robustus, Megaptera novaeangliae, Balaenoptera physalus, Eubalaena australis, Eubalaena glacialis, and Eubalaena japonica. To confirm this scenario, we performed a Bayesian approach, which returned the same grouping by size presented previously. This may be due to factors such as evolutionary convergence or rapid evolution of regulatory elements in this particular gene. To explore this further, we performed additional analyses to investigate the rate of evolution and the dynamics of changes in the TFBS of the promoters.

Fig. 1
figure 1

The adult average size, in meters, of all cetacean species included in this study. Blue values indicate giant cetaceans, and red values non-giant cetaceans. Gigantism in this group is defined by body measurements with an average length of 10 m. Physeter catodon, Eschrichtius robustus, Eubalaena japonica, Eubalaena australis, Eubalaena glacialis, Megaptera novaeangliae, Balaenoptera physalus, and Balaenoptera musculus are classified as giants. Size values are from the “Encyclopedia of Marine Mammals” and the phylogeny from McGowen et al., 2020 [20]

Full size image
Fig. 2
figure 2

The formation of a clade with only giant cetaceans (the sperm whale odontocete and the other large mysticetes) and another with smaller cetaceans (the mysticete minke whale [Balaenoptera acutorostrata] alongside the remaining odontocetes) in both Maximum likelihood tree generated by IQ − TREE and Bayesian tree generated by Mr Bayes v3.2.6 constructed from the promoter region of the NCAPG gene. Numbers under nodes represent bootstrap support (right) and Bayesian posterior probability (left)

Full size image

Regulatory regions analyses

To gain insights into the molecular evolution of non-coding regulatory regions of genes associated with body size in cetaceans, we employed a scanning approach using Ciiider to identify transcription factor binding sites (TFBS) for the nine promoters of interest across all species in our dataset. The identification of TFBS is crucial to understanding the regulatory mechanisms that control gene expression, and that may contribute to the evolution of morphological traits such as body size. By analyzing the presence and distribution of TFBS in these promoters, we aimed to identify potential regulatory modifications that may have contributed to gigantism in cetaceans and to gain a better understanding of the molecular basis of body size evolution in this group of mammals. The scanning approach performed in Ciider identified TFBS for all nine promoters in all species. The results of the enrichment analyses showed conservation patterns across different mammalian groups, suggesting that some transcription factors are evolutionarily conserved, as all mammals have the same patterns of transcription factor enrichment in the same approximate location of the promoter (Table 1, also Additional file 1: Supplementary Figures S9-S14 shows GHSR, IGF2, IGFBP2, LCORL, PLAG1, and ZFAT). For example, in the EGF promoter, TCF7 and CDX1 were found to be spatially conserved across phylogeny at the − 500 bp position within the promoter region, as shown in Fig. 3, which includes all species analysed to highlight this conservation. On the other hand, TFBS exclusive to certain groups was observed in the NCAPG and IGFBP7 promoters, as demonstrated in Figs. 4 and 5, respectively. Figures 4 and 5 only show cetaceans to highlight the specificity of the TFBS in the promoters related to giant and non-giant cetaceans. The complete figures, including other species, can be found in the supplementary material.

Table 1 The top ten over-represented transcription factors for each promoter of interest. These results are obtained when comparing the studied promoters against a background of human genes
Full size table
Fig. 3
figure 3

Enrichment pattern for the EGF promoter implemented in the Ciiider program. The result shows transcription factors in bars, such as TCF7 and CDX1, conserved in the phylogeny. Mammals are cetaceans, artiodactyls, carnivores, primates, bats, and cingulate (marked with an asterisk). The giant cetaceans are: Balaenopteramusculus, Balaenoptera physalus, Eschrichtius robustus, Eubalaenaaustralis, Eubalaenaglacialis, Eubalaena japonica, Megaptera novaeangliae, and Physeter catodon

Full size image
Fig. 4
figure 4

Closer look at the enrichment pattern for the NCAPG promoter implemented in the Ciiider program. The result shows transcription factors in bars. TEF and PBX1, highlighted in blue, are present only in Physeter catodon, Eschrichtius robustus, Eubalaena japonica, Eubalaena australis, Eubalaena glacialis, Megaptera novaeangliae, Balaenoptera physalus, and Balaenoptera musculus that are giant cetaceans, and FOXP3 and ZBTB33, highlighted in red, only in non-giant cetaceans such Tursiops truncatus, Orcinus orca, Lipotes vexillifer, Phocoena sinus, and Balaenoptera acutorostrata scammoni

Full size image
Fig. 5
figure 5

Closer look at the enrichment pattern for the IGFBP7 promoter implemented in the Ciiider program. The result shows transcription factors in bars. PAX2, highlighted in blue, is present in a triple pattern only in Physeter catodon, Eschrichtius robustus, Eubalaena japonica, Eubalaena australis, Eubalaena glacialis, Megaptera novaeangliae, Balaenoptera physalus, and Balaenoptera musculus, which are classified as giants

Full size image

In the NCAPG promoter, we identified a pattern of enrichment that split cetaceans into two groups: giants and non-giants. The giant mysticetes had an enrichment pattern with the transcription factors TEF and PBX1 in the region between − 1300 and − 1200 bp positions, shared only with the sperm whale (Physeter catodon), a species of odontocete that can exceed 20 m in length. In contrast, all cetaceans not classified as giants showed the enrichment of the transcription factor FOXP3 in the − 800 bp position and the transcription factor ZBTB33 in the − 200 bp position.

Similar patterns were also found in the IGFBP7 gene, with giant mysticetes presenting a unique triple enrichment pattern at the − 1100 bp position with the transcription factor PAX2 shared only with the sperm whale.

Additionally, we used phyloP from the PHAST package to estimate the molecular evolution rate of the promoters and identify signals of evolutionary acceleration in specific branches. Specifically, we aimed to identify whether promoters of giant cetaceans underwent accelerated evolution compared to non-giants. To achieve this, phyloP calculated the conservation and acceleration scores in a partitioned tree through a set of named branches, the giant cetaceans, and all remaining species .Thus, the tests for conservation/acceleration occur in the set of named branches relative to the others. Positive scores indicate conservation and negative scores indicate acceleration. A substitution model, against which all subtrees were compared, was derived from the phyloFit program, from the same PHAST package. Our analysis revealed possible evidence of accelerated evolution in the promoters of gigantic cetaceans, as evidenced by negative scores in the IGF2, IGFBP2, IGFBP7, and ZFAT promoters(Table 2).

Table 2 Conservation or acceleration in promoter sequences of the nine genes studied in this work, estimated based on the likelihood ratio test of phyloP for the subtree comparing an alternative model (alt_subscale) with a free scale parameter (alt_scale) within the given REV substitution model (null_scale). Positive scores indicate evolutionary conservation, and negative scores denote evolutionary acceleration, as observed in the IGF2, IGFBP2, IGFBP7, and ZFAT promoters
Full size table

Discussion

This study investigates the molecular evolution of regulatory regions of genes potentially linked to cetaceans’ gigantism, focusing on the promoter region. We found evidence of enrichment of transcription factors binding sites potentially related to large body size, with distinct patterns between giant and non-giant cetaceans in the IGFBP7 and NCAPG promoters. We also found evidence of acceleration in the IGF2, IGFBP2, IGFBP7, and ZFAT promoters. We will focus our discussion on these 5 promoters, as the other four (EGF, GHSR, LCORL, and PLAG1) did not yield relevant results for our research question.

Despite being non-coding regions, which are often known to be difficult to align and contain many neutrally evolving sites and potentially a few constrained ones, we obtained high-quality alignments from our promoters, with ultimately recovered phylogenetic trees consistent with known relationships among species, except for the NCAPG. In this case, the phylogenetic signal was strong enough for the sperm whale (Physeter catodon) odontocete to be grouped with the other mysticetes, excluding the minke whale, grouped with the other odontocetes. In this way, two clades of cetaceans are redeemed: one that contains only those classified as giants and the other with non-giant cetaceans. The use of the Bayesian approach also resulted in the formation of the same clades divided by size. The phylogenetic incongruence between the highly reliable species tree and the promoter tree is a common phenomenon across the Tree of Life, as different regions can have different evolutionary histories [21] due to mechanisms such as incomplete lineage sorting (ILS), introgression, or convergent evolution [22,23,24]. The last one, convergent evolution, could fit the scenario of this work since we have species from two evolutionarily distinct groups (odontocetes and mysticetes) with similar gigantism-related mechanisms. Moreover, as discussed in the following paragraph, the enrichment analysis provides evidence that the convergent evolution of this region is a plausible explanation for this case. Regarding the other eight promoters, the recovery of trees consistent with the most accepted phylogenetic hypotheses for the groups included in the study gives us more confidence that we are indeed using a fundamental regulatory region of the genes in our dataset. Additionally, as discussed further, the identification of evolutionarily conserved TFBS across different mammalian groups in our study further supports the functional importance and conservation of these regulatory regions.

The analyses implemented in the Ciiider program identified the transcription factors binding sites in promoters. Subsequently, the enrichment test revealed some patterns in the promoters of our dataset. First, the same transcription factors were found in the same approximate position in different mammalian lineages, demonstrating evolutionary conservation, and this is the case for the EGF promoter (Fig. 3). Regulatory elements spatially conserved among different lineages suggest an important biological role, as observed between humans and mice for the Cd247, a gene with functional consequences in systemic autoimmunity [25], and in transcription factors related to growth and development in monocot and dicot lineages [26]. Second, NCAPG and IGFBP7 promoters presented different patterns for giant and non-giant cetaceans. In the NCAPG promoter, this transcription factors distribution pattern is likely responsible for the phylogenetic signal in the promoter tree discussed before. The sperm whale (Physeter catodon) has the transcription factors TEF and PBX1 in the region between − 1300 and − 1200 bp position like other giant mysticetes. In contrast, the minke whale has more similarities with smaller odontocetes than its giant mysticetes relatives. Our results suggest that these regions have undergone different selective pressures and that some of the TFBS may have evolved more rapidly in certain lineages. These findings provide further evidence that the NCAPG promoter has experienced unique evolutionary processes that could contribute to the observed incongruence in the phylogenetic tree.

TEF (Thyrotroph embryonic factor) is a protein that belongs to the proline- and acidic amino acid-rich (PAR) bZIP family and is expressed initially in the embryonic anterior pituitary, whereas in adults, it is involved in controlling the cell cycle and the death of hematopoietic cells [27, 28]. These features make TEF a possible tumor suppressor, as demonstrated in bladder cancer (BC). The upregulation significantly retarded BC cell growth by inhibiting the G1/S transition via regulating AKT/FOXOs signaling [28]. In the same way, PBX1 (Pre-B-cell leukemia homeobox 1) is a member of the Three Amino acid Loop Extension (TALE)-class homeodomain family. It is responsible for diverse developmental processes, including skeleton patterning, hematopoiesis, pancreas, and urogenital systems organogenesis [29,30,31,32,33]. It is also involved in fetal growth in activity with decidual natural killer (dNK) cells, driving transcription of pleiotrophin and osteoglycin in dNK cells. On the other hand, the PBX1 inactivation in mouse dNK cells impairs fetal development by decreasing growth-promoting factors that result in fetal growth restriction [34].

Together, both TEF and PBX1 factors are related to general growth processes, such as the control of cell proliferation or directly linked to embryonic growth like PBX1, highlighting the biological meaning of their enrichment pattern only in giant cetaceans, mainly when this enrichment occurs in the promoter of a gene strongly associated with increased body size, such as NCAPG.

The NCAPG (Non-SMC Condensin I Complex Subunit G) gene was previously associated with increased body size and weight gain in horses, donkeys, pigs, humans, and chickens [35,36,37,38,39,40,41]. In bovine species, evolutionarily close to the cetaceans, NCAPG is associated with many essential features such as birth weight, wither height, feeding efficiency, and pubertal growth [42,43,44]. In previous work from our group—focusing on coding regions—evolutionary analyses showed that the NCAPG gene has evidence of positive selection in giant cetaceans [45]. Our promoter and coding regions results imply this gene’s essential role in cetacean gigantism.

The IGFBP7 promoter also showed a specific triple pattern transcription factor only shared by giant cetaceans: the PAX2 (Paired Box Gene 2), which is critical during the embryonic development of systems such as the central nervous system (brain and spinal cord), kidney, eye, ear, and urogenital tract [46, 47]. PAX2 deficiency has been associated with various growth defects, such as kidney hypoplasia, optic coloboma, and vesicoureteral reflux [48]. Furthermore, PAX2 role in embryo development and oncogenesis suggests that it works as a regulatory factor in cell growth [49, 50]. This feature is similar to the IGFBP7 gene, a member of the IGFBP superfamily responsible for the viability of insulin-like growth factors (IGFs) that are molecules involved in promoting cell growth and division [51]. This gene also acts as an oncosuppressor in prostate, breast, lung, and colorectal cancer due to its regulatory action related to cell proliferation, cell adhesion, cell senescence, and angiogenesis [52,53,54]. One of the main challenges of gigantism is the suppression of tumors due to a large number of cells. Therefore, mechanisms that manage to mitigate cancerous processes were crucial during the evolutionary history of the giants.

The cetaceans not classified as giants in this work comprise Tursiops truncatus, Orcinus orca, Lipotes vexillifer, Phocoena sinus, and Balaenoptera acutorostrata scammoni. In the enrichment analyses performed in Ciiider, only these cetacean species share the transcription factors FOXP3 and ZBTB33 in the NCAPG promoter. The first, FOXP3 (Forkhead box protein P3), is a transcription factor belonging to the forkhead box protein family and may act as a transcriptional activator or repressor [55]. It is also associated with the differentiation and function of regulatory T (Treg) cells, which are responsible for suppressing the activation of other leukocytes and thus contribute to immune homeostasis [56,57,58].

The ZBTB33 (Zinc finger and BTB domain-containing 33, also known as Kaiso), exhibits bimodal DNA recognition and acts as a transcriptional repressor and activator depending on the sequence context and cellular phenotype [59]. As a repressor, it recruits other repressors, forming further complexes and aiding in dampening the transcription of the target gene by blocking the binding of transcriptional activators [60]. One of the targets of the transcriptional repressor action of ZBTB33 is the Wnt signaling pathway, associated with critical physiological activities such as growth, differentiation, and migration during development [61]. Focusing on growth, Wnt signaling shapes growing tissues while inducing cells to proliferate, acting as growth factors, and directly affecting cellular organization by the cytoskeleton and mitotic spindle [62]. In summary, the presence of transcription factors that can act as repressors in the promoter of the NCAPG gene related to body growth only in small cetaceans may indicate how these animals did not develop giant sizes.

We found evidence of accelerated evolution in IGF2, IGFBP2, IGFBP7, and ZFAT promoters. The first three (IGF2, IGFBP2, and IGFBP7) are a group of genes that work together to promote growth. The insulin-like growth factors (IGFs), such as IGF2, are important in somatic growth and cell proliferation and responsible for fetal and post-natal growth [63]. This action is only completed by the modulation of insulin-like growth factor binding proteins (IGFBPs), a group that serves as transport proteins for insulin-like growth factors, regulating the bioavailability and function of IGFs [64]. For this direct growth-promoting action, the evidence of evolution acceleration on the promoters found by phyloP in the giant cetaceans follows the knowledge about these genes and reinforces their coordinated performance. Furthermore, the IGFBP7 coding sequence was also associated with positive selection in investigating gigantism in cetaceans [45]. Likewise, the ZFAT gene has been associated with height in multiple human populations in horse body size and has been reported to have crucial roles in the maintenance and differentiation of the adipocytes, the number of T cells, and embryonic development [65,66,67,68]. Therefore, they are likely associated with growth due to controlling various aspects of body enlargement and acting as tumor suppressors. The remaining promoters (EGF, GHSR, LCORL, NCAPG, and PLAG1) exhibit conservation, as identified by CiiiDER, which found highly conserved patterns in most of the genes of interest, with NCAPG and IGFBP7 showing conservation specifically in cetacean groups. Notably, IGFBP7, which also underwent evolutionary acceleration as detected by phyloP, may be associated with multiple gene functions, including body growth and tumor suppression.

Recent studies in other lineages have also highlighted the importance of regulatory regions in controlling body size. For instance, a deletion in the promoter region of IGF2BP1 has been associated with larger body sizes in chickens [69], and variation in the STAT3 promoter has been shown to contribute to larger body size traits in cattle [70]. Additionally, the control of growth hormone IGF1 protein levels by long non-coding RNA has been implicated in the size of large dogs [71]. These findings, along with our own, underscore the critical role that regulatory regions play in determining size characteristics across diverse taxa. Further studies on the molecular evolution of these regions are needed, and future experimental testing will provide further insights into the regulatory mechanisms underlying body size variation.

Although with some limitations, such as the number of genes used, our study provides the first steps toward what other works can reach, especially those related to experimental validation. It is far from the definitive answer to a complex question. Still, this start could be useful in future research, indicating which genes are possibly related to gigantism in cetaceans and that this phenomenon must be understood in an integrated way.

Conclusions

We investigated the promoter regions of genes possibly associated with increased body size in giant cetaceans. In summary, we found evolutionary conservation and evidence of differential transcription factors enrichment, with distinct patterns between giants and non-giants cetaceans for IGFBP7 and the NCAPG promoters. In NCAPG, observing the presence of repressive transcription factors only in cetaceans of small body-size was also possible. Furthermore, evolutionarily acceleration was detected in the promoters of the IGF2, IGFBP2, IGFBP7, and ZFAT genes. In conclusion, our study provides evidence of the evolution of cetacean gigantism from a regulatory approach.

Materials and methods

Sample data

The promoters of nine genes were chosen because they have been described in the scientific literature as associated with changes in body size. The EGF (Epidermal Growth Factor), GHSR (Growth Hormone Secretagogue Receptor), IGF2 (Insulin-Like Growth Factor 2), IGFBP2 (Insulin-Like Growth Factor Binding Protein 2), and IGFBP7 (Insulin-Like Growth Factor Binding Protein 7) are part of the growth hormone/insulin-like growth factor (GH-IGF) axis, which plays a central role in regulating growth in vertebrates [72, 73]. The LCORL (Ligand Dependent Nuclear Receptor Corepressor Like), NCAPG (Non-SMC Condensin I Complex Subunit G), PLAG1 (Pleomorphic Adenoma Gene 1), and ZFAT (Zinc Finger And AT-Hook Domain Containing) are associated with the body enlargement of species such as cows, pigs, sheep, and goats, which are artiodactyls, evolutionarily close to cetaceans [74, 75]. The sequences of these promoters were retrieved in the Eukaryotic Promoter Database (EPD) from the Swiss Institute of Bioinformatics. Firstly, we located the transcription start site (TSS) for the human species and selected a 1500 bp region upstream of the TSS. Then, the promoter sequences of cetacean and other mammalian species were searched in public databases, such as Ensembl and GenBank (NCBI), using BLAST (Basic Local Alignment Search Tool), which compares nucleotide or protein sequences and calculates the statistical significance, finding similarity regions among sequences of interest.

For cetaceans, we used sequences from 13 species, five odontocetes (Tursiops truncatus, Orcinus orca, Lipotes vexillifer, Physeter catodon, and Phocoena sinus), and eight mysticetes (Balaenoptera acutorostrata scammoni, Eschrichtius robustus, Megaptera novaeangliae, Balaenoptera physalus, Balaenoptera musculus, Eubalaena australis, Eubalaena glacialis, and Eubalaena japonica). The sequences for Eubalaena australis and Eubalaena glacialis were retrieved from genomes available on the public platform DNA Zoo. All other cetacean sequences were retrieved from GenBank, and the Additional file 1: Supplementary Table 1 shows the accession numbers.

Following Lambert et al. 2010 [76], gigantism is attributed to species larger than 10 m. In our dataset, the following species fit this definition: blue whale (Balaenoptera musculus), sperm whale (Physeter catodon), gray whale (Eschrichtius robustus), humpback whale (Megaptera novaeangliae), fin whale (Balaenoptera physalus), South Atlantic right whale (Eubalaena australis), North Atlantic right whale (Eubalaena glacialis), and Pacific right whale (Eubalaena japonica).

In addition to cetaceans, we included 39 other species to represent the major mammalian groups, such as the order Artiodactyla (Bos taurus, Capra hircus, Bison bison, Odocoileus virginianus, Ovis aries, Sus scrofa, Camelus dromedarius, Camelus ferus, Camelus bactrianus), Carnivora (Panthera leo, Panthera onca, Panthera tigris altaica, Panthera pardus, Felis catus, Prionailurus bengalensis, Canis lupus familiaris, Canis lupus dingo, Vulpes lagopus, Vulpes vulpes, Ursus arctos horribilis, Ursus thibetanus), Primate (Homo sapiens, Pan paniscus, Pan troglodytes, Gorilla gorilla gorilla, Pongo abelii, Callithrix jacchus, Rhinopithecus roxellana, Papio anubis, Nomascus leucogenys, Macaca fascicularis, Macaca mulatta and Chlorocebus sabaeus), Cingulata (Dasypus novemcinctus), and Chiroptera (Artibeus jamaicensis, Hipposideros armiger, Phyllostomus discolor, Pipistrellus pipistrellus, Rhinolophus ferrumequinum). Thus, there were the same 52 species in each promoter studied.

Phylogenetic reconstructions

The sequences were aligned using the MUSCLE program [77] and visualized in AliView [78]. After this, phylogenetic trees were constructed for each promoter using the IQ-TREE program’s maximum likelihood strategy, 1,000 bootstrap replicates to estimate branch confidence, 1,000 maximum iterations, 1,000 number of bootstrap alignments, 0.5 perturbation strength, 100 IQ-TREE stopping rule, 0.99 minimum correlation coefficient, and “auto” in substitution model. This entire process was done directly on the IQ-TREE Web Server portal [79]. For Bayesian analysis, we determine the optimal number of partitions and evolutionary models for each promoter using PARTITION FINDER software v2.1.1 [80], which employed the Bayesian Information Criterion (BIC). Subsequently, Bayesian phylogenetic trees were constructed using MrBayes v3.2.6 [81]. The Markov chain Monte Carlo (MCMC) algorithm was run for 5,000,000 generations with four chains, and trees were sampled every 100 generations, utilizing the molecular evolution model selected by PARTITION FINDER v2.1.1. The resulting trees were visualized using FigTree v1.3.1. Finally, we visualized the results in the program FigTree v1.3.1.

Regulatory regions analyses

Promoter analyses were performed using Ciiider and phyloP tools. Ciiider was used to predict and to analyze transcription factor binding sites within a sequence and identify significantly enriched ones [82]. This is important since over-represented transcription factors are more likely to regulate gene expression that ultimately alters the phenotype [83]. We used scanning and enrichment approaches in Ciiider.

Given a sequence, the scanning predicts potential transcription factors in the region of interest. The MATCH algorithm searches for transcription factor binding sites in DNA sequences [84] using a Position Frequency Matrix (PFM). A set of non-redundant profiles derived from experimentally defined transcription factor binding sites for eukaryotes is used in this work, derived from the JASPAR database containing position matrixes of these elements [85]. Since PFMs generally have a highly conserved core-binding region flanked by areas of higher variability, a core PFM is created for the five most conserved consecutive bases. To search for transcription factor binding sites, sequences are divided into regions of five overlapping bases compared to the core PFM. If the similarity score between a five-base sequence and the core PFM meets a defined threshold, then the sequence window is increased to the full length of the transcription factors, and the similarity score to the full PFM is calculated. The default deficit is 0.15, meaning the scan will accept any transcription factors with MATCH scores of 0.85 or above [84].

The enrichment approach allow us to identify those transcription factor binding sites that are significantly over- or under-represented in the regions of interest when compared to the background regions used in the analysis. To reduce the possibility of chance findings, we used a comparative background consisting of several other genes provided by the Ciiider program. Thus, we reduce the chances of the results being stochastic. In short, Ciiider scans these background sequences using the same criteria for the sequences of interest and for the background sequences. The program determines the over- and under-representation of transcription factors by comparing the number of sequences containing these factors to those without them, followed by a statistical test such as Fisher’s exact test [84].

We used phyloP tool from PHAST (Phylogenetic Analysis with Space/Time) package to estimate the molecular evolution rate of the promoter and detect signals of evolutionary acceleration in specific branches [86, 87]. First, we generated a substitution model using the phyloFit program, which fits one or more tree models to multiple alignments of DNA sequences using maximum likelihood, and the substitution model used was REV (Reversible Evolutionary Model), the default of phyloFit, which is more realistic and flexible than simpler neutral models and can capture variations in nucleotide substitution rates at different positions in the alignment [88, 89]. Using REVl, we calculated conservation and acceleration scores with the “branch” option, which partitions the tree into named branches and tests for conservation/acceleration in the named branches relative to the others. We compared the set of named branches containing giant cetaceans against the remaining species. We selected the LRT option, which compares an alternative model having a free scale parameter with the substitution model, and the CONACC mode, which allows for acceleration as well as conservation, assigning positive values (scores) to indicate conservation and negative values to indicate acceleration. Thus, CONACC mode summarizes conservation and acceleration.

Data Availability

All the data supporting our findings are contained within the manuscript and in the supplemental file.

Abbreviations

ACAN:

Aggrecan

BLAST:

Basic Local Alignment Search Tool

CBS :

Cystathionine β – synthase

EGF :

Epidermal Growth Factor

EIF2AK3 :

Translation initiation factor2-α kinase 3

EPD:

Eukaryotic Promoter Database

GHSR :

Growth Hormone Secretagogue Receptor

GRB10 :

Growth factor receptor-bound protein 10

IGF :

Insulin-Like Growth Factor 2

IGFBP2 :

Insulin-Like Growth Factor Binding Protein 2

IGFBP7 :

Insulin-Like Growth Factor Binding Protein 7

LCORL :

Ligand Dependent Nuclear Receptor Corepressor Like

NCAPG :

Non-SMC Condensin I Complex Subunit G

OBSL1 :

Obscurin-like protein 1

PFM:

Position Frequency Matrix

PLAG1 :

Pleomorphic Adenoma Gene 1

PLOD1:

lysyl hydroxylase 1 gene

TSS:

Transcription Start Site

ZFAT :

Zinc Finger And AT-Hook Domain Containing

References

  1. Mead JG, Brownell R. L. Order Cetacea in Mammal Species of the World: A Taxonomic and Geographic Reference (eds. Wilson, D. E. & Reeder, D. M.). United States of America: University Press, Cambridge; 2005;723–743.

  2. Thewissen J. Whales originated from aquatic artiodactyls in the Eocene epoch of India. Nature. 2007;450:1190–4.

    CAS  PubMed  Google Scholar 

  3. Berta A. etal.Marine mammals: evolutionary biology. United States of America. Cambridge: Academic Press; 2005;178–94.

    Google Scholar 

  4. Leffler EM. Revisiting an old riddle: what determines genetic diversity levels within species? PLoS Biol. 2012. https://doi.org/10.1371/journal.pbio.1001388.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Damuth J. Population density and body size in mammals. Nature. 1981;290:699–700.

    Google Scholar 

  6. Jones ML, Swartz L. Gray Whale: Eschrichtius robustus. Encyclopedia of Marine Mammals 2nd edition. United States of America: Academic Press, Cambridge; 2009;503–511.

  7. Sears R. Perrin WF. Blue Whale: Balaenoptera musculus. Encyclopedia of Marine Mammals 2nd edition. United States of America: Academic Press, Cambridge; 2009;120–124.

    Google Scholar 

  8. Downhower JF, Bulmer LS. Calculating just how small a whale can be. Nature. 1988;335:675.

    Google Scholar 

  9. Smith FA, Lyons SK. How big should a mammal be? A macroecological look at mammalian body size over space and time. Phil Trans R Soc B. 2011;366:2364–78.

    PubMed  PubMed Central  Google Scholar 

  10. Goldbogen JA, Madsen PT. The evolution of foraging capacity and gigantism in cetaceans. Exp Biol. 2018. https://doi.org/10.1242/jeb.166033.

    Article  Google Scholar 

  11. Whitehead H. Sperm Whale: Physeter macrocephalus. Encyclopedia of Marine Mammals 3th edition. United States of America: Academic Press, Cambridge; 2018;919–925.

  12. Goldbogen JA. Why whales are big but not bigger: physiological drivers and ecological limits in the age of ocean giants. Science. 2019;366:1367–72.

    CAS  PubMed  Google Scholar 

  13. Sun Y. Insights into body size variation in cetaceans from the evolution of body-size related genes. BMC Evol Bio. 2019. https://doi.org/10.1186/s12862-019-1461-9.

    Article  Google Scholar 

  14. Carroll SB. Evo-Devo and an expanding evolutionary synthesis: a genetic theory of morphological evolution. Cell. 2008;134:25–36.

    CAS  PubMed  Google Scholar 

  15. Dermitzakis ET, Clark AG. Evolution of transcription factor binding sites in mammalian gene regulatory regions: conservation and turnover. Mol Biol Evol. 2002;19(7):1114–21.

    CAS  PubMed  Google Scholar 

  16. Lowe CB. Three periods of regulatory innovation during vertebrate evolution. Science. 2011;333:1019–24.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Karlsson EK. Efficient mapping of mendelian traits in dogs through genome-wide association. Nat Genet. 2007;39(11):1321–8.

    CAS  PubMed  Google Scholar 

  18. McGregor AP. Morphological evolution through multiple cis-regulatory mutations at a single gene. Nature. 2007;448:587–91.

    CAS  PubMed  Google Scholar 

  19. Chan YF. Adaptive evolution of pelvic reduction in sticklebacks by recurrent deletion of a Pitx1 enhancer. Science. 2010;327:302–5.

    CAS  PubMed  Google Scholar 

  20. McGowen MR. Phylogenomic resolution of the Cetacean Tree of Life using Target sequence capture. Syst Biol. 2020;69:479–501.

    CAS  PubMed  Google Scholar 

  21. Degnan JH, Rosenberg NA. Discordance of Species Trees with their most likely gene trees. PLoS Genet. 2006. https://doi.org/10.1371/journal.pgen.0020068.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Maddison WP. Gene trees in species trees. Syst Biol. 1997;46:523–36.

    Google Scholar 

  23. Satta Y. DNA archives and our nearest relative: the trichotomy problem revisited. Mol Phylogenet Evol. 2000;14:259–75.

    CAS  PubMed  Google Scholar 

  24. Parker J. Genome-wide signatures of convergent evolution in echolocating mammals. Nature. 2013;502:228–31.

    CAS  PubMed  Google Scholar 

  25. Pundhir S. Spatially conserved regulatory elements identified within human and mouse Cd247 gene using high-throughput sequencing data from the ENCODE project. Gene. 2014;545:80–7.

    CAS  PubMed  Google Scholar 

  26. Xu F. Cis-regulatory signatures of orthologous stress-associated bZIP transcription factors from rice, sorghum and Arabidopsis based on phylogenetic footprints. BMC Genom. 2012;13:497–512.

    CAS  Google Scholar 

  27. Inukai T. TEF, an antiapoptotic bZIP transcription factor related to the oncogenic E2A-HLF chimera, inhibits cell growth by down-regulating expression of the common β chain of cytokine receptors. Blood. 2005;105:4437–44.

    CAS  PubMed  Google Scholar 

  28. Yang J. Thyrotroph embryonic factor is downregulated in Bladder cancer and suppresses proliferation and tumorigenesis via the AKT/FOXOs signalling pathway. Cell Prolif. 2019. https://doi.org/10.1111/cpr.12560.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Selleri L. Requirement for Pbx1 in skeletal patterning and programming chondrocyte proliferation and differentiation. Development. 2001;128:3543–57.

    CAS  PubMed  Google Scholar 

  30. Specchia G. Extramedullary involvement at relapse in acute promyelocytic Leukemia patients treated or not with all-trans retinoic acid: a report by the Gruppo Italiano Malattie Ematologiche dell’Adulto. J Clin Oncol. 2001;19:4023–8.

    CAS  PubMed  Google Scholar 

  31. Kim SK. Pbx1 inactivation disrupts pancreas development and in Ipf1-deficient mice promotes Diabetes Mellitus. Nat Genet. 2002;30:430–5.

    CAS  PubMed  Google Scholar 

  32. Schnabel CA. Pbx1 is essential for adrenal development and urogenital differentiation. Genesis. 2003;37:123–30.

    CAS  PubMed  Google Scholar 

  33. Magnani L. PBX1 genomic pioneer function drives ERα Signaling underlying progression in Breast Cancer. PLoS Genet. 2011. https://doi.org/10.1371/journal.pgen.1002368.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Zhou Y. PBX1 expression in uterine natural killer cells drives fetal growth. Sci Transl Med. 2020. https://doi.org/10.1126/scitranslmed.aax1798.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Tetens J. A genome-wide association study indicates LCORL/NCAPG as a candidate locus for withers height in German Warmblood horses. Anim Genet. 2013;44:467–71.

    CAS  PubMed  Google Scholar 

  36. Shen J. Genomic Analyses Reveal Distinct Genetic Architectures and selective pressures in Chinese donkeys. J Genet Genom. 2021;48:737–45.

    Google Scholar 

  37. Rubin CJ, et al. Strong signatures of selection in the domestic pig genome. Proc. Natl. Acad. Sci. USA 2012;109:19529–19536.

  38. Gudbjartsson DF. Many sequence variants affecting diversity of adult human heigh. Nat Genet. 2008;40:609–15.

    CAS  PubMed  Google Scholar 

  39. Lettre G. Identification of ten loci associated with height highlights new biological pathways in human growth. Nat Genet. 2008;40:584–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Weedon MN. Genome-wide association analysis identifies 20 loci that influence adult height. Nat Genet. 2008;40:575–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Sasaki S. Genetic mapping of quantitative trait loci affecting body weight, egg character and egg production in F2intercross chickens. Anim Genet. 2004;35:188–94.

    CAS  PubMed  Google Scholar 

  42. Eberlein A. Dissection of genetic factors modulating fetal growth in cattle indicates a substantial role of the Non-SMC condensin I Complex, Subunit G (NCAPG) Gene. Genetics. 2009;183:951–64.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. Weikard R. Metabolomic profiles indicate distinct physiological pathways affected by two loci with major divergent effect on Bos taurus growth and lipid deposition. Physiol Genom. 2010;42:79–88.

    Google Scholar 

  44. Setoguchi K. The SNP c.1326T > G in the non-SMC condensin I complex, subunit G (NCAPG) gene encoding a p.Ile442Met variant is associated with an increase in body frame size at puberty in cattle. Anim Genet. 2011;42:650–5.

    CAS  PubMed  Google Scholar 

  45. Silva FA. The molecular evolution of genes previously associated with large sizes reveals possible pathways to cetacean gigantism. Sci Rep. 2022. https://doi.org/10.1038/s41598-022-24529-3.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Eccles MR. PAX genes in development and Disease: the role of PAX2 in urogenital tract development. Int J Dev Biol. 2002;46:535–44.

    CAS  PubMed  Google Scholar 

  47. Porteous S. Primary renal hypoplasia in humans and mice with PAX2 mutations: evidence of increased apoptosis in fetal kidneys of Pax21Neu +/– mutant mice. Hum Mol Genet. 2000;9:1–11.

    CAS  PubMed  Google Scholar 

  48. Weber S. SIX2 and BMP4 mutations associate with anomalous kidney development. J Am Soc Nephrol. 2008;19:891–903.

    CAS  PubMed  PubMed Central  Google Scholar 

  49. Jahangiri R. PAX2 expression is correlated with better survival in tamoxifen-treated breast carcinoma patients. Tissue Cell. 2018;52:135–42.

    CAS  PubMed  Google Scholar 

  50. Song H. PAX2 expression in Ovarian Cancer. Int J Mol Sci. 2013;14:6090–105.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Burger AM. Essential roles of IGFBP-3 and IGFBP-rP1 in Breast cancer. Eur J Cancer. 2005;41:1515–27.

    CAS  PubMed  Google Scholar 

  52. Akaogi K, et al. Specific accumulation of tumor-derived adhesion factor in tumor blood vessels and in capillary tube-like structures of cultured vascular endothelial cells. Proc. Natl. Acad. Sci. USA. 1996;93:8384–8389.

  53. Sprenger CC. Insulin-like growth factor binding protein-related protein 1 (IGFBP-rP1) is a potential Tumor suppressor protein for Prostate cancer. Cancer Res. 1999;59:2370–5.

    CAS  PubMed  Google Scholar 

  54. Wilson HM. Insulin-like growth factor binding protein-related protein 1 inhibits proliferation of MCF-7 Breast cancer cells via a senescence-like mechanism. Cell Growth Differ. 2002;13:205–13.

    CAS  PubMed  Google Scholar 

  55. Jia H. The expression of FOXP3 and its role in human cancers. Biochim Biophys Acta Rev Cancer. 2019;187:170–8.

    Google Scholar 

  56. Hori S. Control of regulatory T cell development by the transcription factor Foxp3. Science. 2003;299:1057–61.

    CAS  PubMed  Google Scholar 

  57. Fontenot JD. Foxp3 programs the development and function of CD4 + CD25 + regulatory T cells. Nat Immunol. 2003;4:330–6.

    CAS  PubMed  Google Scholar 

  58. Lu L. The regulation of immune tolerance by FOXP3. Nature. 2017;17:703–17.

    CAS  Google Scholar 

  59. Pozner A. Cell specific Kaiso (ZBTB33) regulation of cell cycle through cyclin D1 and cyclin E1. J Biol Chem. 2016;291:24538–50.

    CAS  PubMed  PubMed Central  Google Scholar 

  60. Yoon HG. N-CoR mediates DNA methylation-dependent repression through a methyl CpG binding protein Kaiso. Mol Cell. 2003;12:723–34.

    CAS  PubMed  Google Scholar 

  61. Azbazdar Y 1. Regulation of wnt signaling pathways at the plasma membrane and their Misregulation in Cancer. Front Cell Dev Biol. 2021. https://doi.org/10.3389/fcell.2021.631623.

  62. Nusse R, Clevers H. Wnt/β-Catenin signaling, Disease, and emerging therapeutic modalities. Cell. 2017;169:985–99.

    CAS  PubMed  Google Scholar 

  63. Monzavi R, Cohen P. IGFs and IGFBPs: role in health and Disease. Best Pract Res Clin Endocrinol Metab. 2002;16:433–47.

    CAS  PubMed  Google Scholar 

  64. Ding H, Wu T. Insulin-like growth factor binding proteins in Autoimmune Diseases. Front Endocrinol. 2018. https://doi.org/10.3389/fendo.2018.00499.

    Article  Google Scholar 

  65. Takeuchi F. Evaluation of genetic loci influencing adult height in the Japanese population. J Hum Genet. 2009;54:749–52.

    PubMed  Google Scholar 

  66. Allen HL. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010;467:832–8.

    Google Scholar 

  67. Signer-Hasler H. A genome-wide association study reveals loci influencing height and other conformation traits in horses. PLoS ONE. 2012. https://doi.org/10.1371/journal.pone.0037282.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Makvandi-Nejad S. Four loci explain 83% of size variation in the horse. PLoS ONE. 2012. https://doi.org/10.1371/journal.pone.0039929.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Wang K. The Chicken Pan-genome reveals Gene Content Variation and a promoter region deletion in IGF2BP1 affecting body size. Mol Biol Evol. 2021;38:5066–81.

    CAS  PubMed  PubMed Central  Google Scholar 

  70. Wu S. Genetic variants in STAT3 promoter regions and their application in molecular breeding for body size traits in Qinchuan Cattle. Int J Mol Sci. 2018. https://doi.org/10.3390/ijms19041035.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Plassais J. Natural and human-driven selection of a single non-coding body size variant in ancient and modern canids. Curr Biol. 2022;32:889–97.

    CAS  PubMed  PubMed Central  Google Scholar 

  72. Canosa LF. Neuroendocrine control of growth hormone in fish. Gen Comp Endocrinol. 2007;151:1–26.

    CAS  PubMed  Google Scholar 

  73. Bergan-Roller HE, Sheridan MA. The growth hormone signaling system: insights into coordinating the anabolic and catabolic actions of growth hormone. Gen Comp Endocrinol. 2018;258:119–33.

    CAS  PubMed  Google Scholar 

  74. Takasuga A. PLAG1 and NCAPG-LCORL in livestock. Anim Sci J. 2016;87:159–67.

    CAS  PubMed  Google Scholar 

  75. Wang W. Molecular characterization and expression of SPP1, LAP3 and LCORL and their association with growth traits in sheep. Genes. 2019. https://doi.org/10.3390/genes10080616.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Lambert O. The giant bite of a new raptorial sperm whale from the Miocene epoch of Peru. Nature. 2010;466:105–8.

    CAS  PubMed  Google Scholar 

  77. Edgar RC, Muscle. Multiple sequence alignment with high accuracy and high-throughput. Nucleic Acids Res. 2004;32:1792–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  78. Larsson A, Aliview. A fast and lightweight alignment viewer and editor for large data sets. Bioinfo. 2014;22:3276–8.

    Google Scholar 

  79. Nguyen LT. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.

    CAS  PubMed  Google Scholar 

  80. Lanfear R, et al. PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2016;34:772–3.

    Google Scholar 

  81. Ronquist F, Huelsenbeck JP. MrBayes 3: bayesian phylogenetic inference under mixed models. J Bioinform. 2003;19:1572–4.

    CAS  Google Scholar 

  82. Gearing LJ. CiiiDER: a tool for predicting and analysing transcription factor binding sites. PLoS ONE. 2019. https://doi.org/10.1371/journal.pone.0215495.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Boeva V. Analysis of genomic sequence motifs for deciphering transcription factor binding and transcriptional regulation in eukaryotic cells. Front Genet. 2016. https://doi.org/10.3389/fgene.2016.00024.

    Article  PubMed  PubMed Central  Google Scholar 

  84. Kel AE. MATCH: a tool for searching transcription factor binding sites in DNA sequences. Nucleic Acids Res. 2003;31:3576–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  85. Castro-Mondragon, JA, et al. JASPAR 2020: Update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. 2020;48:87–92.

  86. Pollard KS. Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 2010;20:110–21.

    CAS  PubMed  PubMed Central  Google Scholar 

  87. Hubisz MJ. PHAST and RPHAST: phylogenetic analysis with space/time models. Brief Bioinform. 2011;12:41–51.

    CAS  PubMed  Google Scholar 

  88. Siepel A, Pollard KS, Haussler D. New methods for detecting lineage-specific selection. Comput Mol Biol. 2006. https://doi.org/10.1007/11732990_17.

    Article  Google Scholar 

  89. Siepel A, Haussler D. Phylogenetic estimation of context-dependent substitution rates by maximum likelihood. Mol Biol Evol. 2004;21:468–88.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was funded by Coordination for the Improvement of Higher Education Personnel—Brazil (CAPES)—Finance Code 001 and FAPESP (2015/18269-1).

Author information

Authors and Affiliations

  1. Dept of Genetics, Evolution, Microbiology & Immunology, Institute of Biology, University of Campinas, Rua Monteiro Lobato, 255, Campinas, 13083-862, SP, Brazil

    Felipe A. Silva, Agnello C. R. Picorelli, Giovanna S. Veiga & Mariana F. Nery

Authors
  1. Felipe A. Silva
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Agnello C. R. Picorelli
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Giovanna S. Veiga
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Mariana F. Nery
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

M.F.N. conceived the research hypothesis; F.A.S., A.C.R.P., and G.S.V. analyzed the data and drafted the manuscript; all authors contributed equally to the final version. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Mariana F. Nery.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

12862_2023_2171_MOESM1_ESM.docx

Supplementary Material 1: Additional file 1: Table S1: Species used in this study and the respective accession numbers from NCBI; Figure S1: Maximum likelihood tree generated by IQ − TREE constructed from the promoter region of the EGF gene. Numbers under nodes represent bootstrap support; Figure S2: Maximum likelihood tree generated by IQ − TREE constructed from the promoter region of the GHSR gene. Numbers under nodes represent bootstrap support; Figure S3: Maximum likelihood tree generated by IQ − TREE constructed from the promoter region of the IGF2 gene. Numbers under nodes represent bootstrap support; Figure S4: Maximum likelihood tree generated by IQ − TREE constructed from the promoter region of the IGFBP2 gene. Numbers under nodes represent bootstrap support; Figure S5: Maximum likelihood tree generated by IQ − TREE constructed from the promoter region of the IGFBP7 gene. Numbers under nodes represent bootstrap support; Figure S6: Maximum likelihood tree generated by IQ − TREE constructed from the promoter region of the LCORL gene. Numbers under nodes represent bootstrap support; Figure S7: Maximum likelihood tree generated by IQ − TREE constructed from the promoter region of the PLAG1 gene. Numbers under nodes represent bootstrap support; Figure S8: Maximum likelihood tree generated by IQ − TREE constructed from the promoter region of the ZFAT gene. Numbers under nodes represent bootstrap support; Figure S9: Enrichment pattern for the GHSR promoter implemented in the Ciiider program. The result shows transcription factors in bars. Mammals are cetaceans, artiodactyls, carnivores, primates, bats, and cingulates; Figure S10: Enrichment pattern for the IGF2 promoter implemented in the Ciiider program. The result shows transcription factors in bars. Mammals are cetaceans, artiodactyls, carnivores, primates, bats, and cingulates; Figure S11: Enrichment pattern for the IGFBP2 promoter implemented in the Ciiider program. The result shows transcription factors in bars. Mammals are cetaceans, artiodactyls, carnivores, primates, bats, and cingulates; Figure S12: Enrichment pattern for the LCORL promoter implemented in the Ciiider program. The result shows transcription factors in bars. Mammals are cetaceans, artiodactyls, carnivores, primates, bats, and cingulates; Figure S13: Enrichment pattern for the PLAG1 promoter implemented in the Ciiider program. The result shows transcription factors in bars. Mammals are cetaceans, artiodactyls, carnivores, primates, bats, and cingulates; Figure S14: Enrichment pattern for the ZFAT promoter implemented in the Ciiider program. The result shows transcription factors in bars. Mammals are cetaceans, artiodactyls, carnivores, primates, bats, and cingulates.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Silva, F.A., Picorelli, A.C.R., Veiga, G.S. et al. Patterns of enrichment and acceleration in evolutionary rates of promoters suggest a role of regulatory regions in cetacean gigantism. BMC Ecol Evo 23, 62 (2023). https://doi.org/10.1186/s12862-023-02171-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-023-02171-5

Keywords

BMC Ecology and Evolution

ISSN: 2730-7182

Contact us