Click here to close
Hello! We notice that you are using Internet Explorer, which is not supported by Xenbase and may cause the site to display incorrectly.
We suggest using a current version of Chrome,
FireFox, or Safari.
Nat Commun
2017 May 15;8:15213. doi: 10.1038/ncomms15213.
Show Gene links
Show Anatomy links
Transcriptomic and macroevolutionary evidence for phenotypic uncoupling between frog life history phases.
Wollenberg Valero KC
,
Garcia-Porta J
,
Rodríguez A
,
Arias M
,
Shah A
,
Randrianiaina RD
,
Brown JL
,
Glaw F
,
Amat F
,
Künzel S
,
Metzler D
,
Isokpehi RD
,
Vences M
.
???displayArticle.abstract???
Anuran amphibians undergo major morphological transitions during development, but the contribution of their markedly different life-history phases to macroevolution has rarely been analysed. Here we generate testable predictions for coupling versus uncoupling of phenotypic evolution of tadpole and adult life-history phases, and for the underlying expression of genes related to morphological feature formation. We test these predictions by combining evidence from gene expression in two distantly related frogs, Xenopus laevis and Mantidactylus betsileanus, with patterns of morphological evolution in the entire radiation of Madagascan mantellid frogs. Genes linked to morphological structure formation are expressed in a highly phase-specific pattern, suggesting uncoupling of phenotypic evolution across life-history phases. This gene expression pattern agrees with uncoupled rates of trait evolution among life-history phases in the mantellids, which we show to have undergone an adaptive radiation. Our results validate a prevalence of uncoupling in the evolution of tadpole and adult phenotypes of frogs.
Figure 1. Predictions for coupled and uncoupled trait evolution in life-history phases.
(a) ‘Phase-Pleiotropy' under coupled evolution. The same genes contribute to tadpole and adult morphological features (here: body size). Non-synonymous DNA substitutions are rare due to phase-specific selective regimes. The same genes are equally expressed across phases. Phenotypic evolutionary rates are correlated and rate changes are concentrated on the same parts of the phylogeny (t represents evolutionary time). (b) ‘Phase-Polygeny' model 1: uncoupled evolution with different genes coding for tadpole and adult phenotype. Non-synonymous DNA substitutions are common due to phase-specific selective regimes. Genes are unequally expressed across phases. Parallel genetic changes are rare, only if favored by selection across phases, and facilitated by linkage. Phenotypic evolutionary rates are not correlated and rate changes are unequally distributed on the phylogeny. (c) ‘Phase-Polygeny' model 2: uncoupled evolution with regulatory (gene expression) or epistatic changes among tadpole and adult phenotype. The same genes contribute to tadpole and adult phenotype, but are regulated/expressed differently (mRNA—bars, epistatic gene products—tear drops, crosses, grey circles). Independent changes in expression/regulation are common due to phase-specific selective regimes. Transcripts are unequally expressed across phases. Parallel changes in expression/regulation are rare, occurring only if favoured by selection across phases and facilitated by functional interactions between regulators. Phenotypic evolutionary rates are not correlated and rate changes are unequally distributed on the phylogeny.
Figure 2. Life-history specific gene expression and contribution to phenotype in frogs.
(a,b) Venn diagram showing summary statistics of genes differentially expressed between life-history phases tadpole and adult, and combinations of these in (a) Xenopus laevis (from Xenbase) and (b) Mantidactylus betsileanus (from two newly sequenced transcriptomes). Numbers in Venn diagrams: top row, number of annotated genes expressed per phase; middle row, number of genes associated with morphological feature formation expressed per phase; bottom row, percentage of morphological-feature-formation-associated genes. (c) Morphology-associated gene expression in adult and tadpole M. betsileanus shown as scatterplot of log tadpole FPKM (fragments per kilobase of exon per million reads mapped) versus log adult FPKM. Symbols denote different categories of morphology-associated gene ontology terms; if one gene has more than one term the symbols are partly overlapping. Symbol size represents combined adult and tadpole FPKM value, as shown in the stacked FPKM scale (right). Red labelled symbols denote significantly differentially expressed transcripts (protein symbols). (d) STRING subnetwork (confidence view) of functional interactions between some genes expressed in adult life phase in Xenopus (blue gene symbols) and M. betsileanus (red gene symbols). Expressed genes of both species are linked within the same functional network with ‘osteoclast differentiation' as the significantly overexpressed function. Line thickness represents strength of confidence for interactions.
Figure 3. Phenotypic optima and uncoupling of character evolution in mantellid frogs.
(a) Phylogeny showing high (red) and low (blue) relative rates of phenotypic evolution in mantellids. Black branches show no rate acceleration or deceleration. Triangles mark shifts between selective optima (black, shared in adult and tadpole phenotypes; grey, unique to either of the two phases). A, Aglyptodactylus; Bl, Blommersia; Ge, Gephyromantis; Gu, Guibemantis; Mt, Mantella; Sp, Spinomantis. The focal species Mantidactylus betsileanus is marked with an asterisk. (b) Rates of phenotypic evolution through time, showing uncoupled rates at ca. 24 and 7.5 Myr ago (indicated by a star and highlighted by grey bars). The continuous line represents adults, the hatched line tadpoles. (c) Plot of the number of shifts between selective optima through time, showing two peaks at ca. 25 and 12.5 Myr ago. The second adaptive peak is characterized by uncoupling of shifts between adults and tadpoles (star). (d) Variation of the diversification rate (in orange) and the predicted increase of species richness (in blue) according to the best-supported diversification model found in our study (model 3). This model of diversification assumes an exponential variation of speciation rate with no extinction. (e) Disparity through time (DTT) plots visualizing the dynamics of phenotypic diversification in adults and tadpoles of the family Mantellidae, for the first MDS variable. See Supplementary Fig. 5 for full plots including comparisons with simulated curves. Inset photos show representative mantellids.
Ashburner,
Gene ontology: tool for the unification of biology. The Gene Ontology Consortium.
2000, Pubmed
Ashburner,
Gene ontology: tool for the unification of biology. The Gene Ontology Consortium.
2000,
Pubmed
Berry,
The expression pattern of thyroid hormone response genes in remodeling tadpole tissues defines distinct growth and resorption gene expression programs.
1998,
Pubmed
,
Xenbase
Bossuyt,
Convergent adaptive radiations in Madagascan and Asian ranid frogs reveal covariation between larval and adult traits.
2000,
Pubmed
Christidis,
Insights into the evolution of a cryptic radiation of bats: dispersal and ecological radiation of Malagasy Miniopterus (Chiroptera: Miniopteridae).
2014,
Pubmed
Conesa,
Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research.
2005,
Pubmed
Eastman,
A novel comparative method for identifying shifts in the rate of character evolution on trees.
2011,
Pubmed
Flatt,
Hormonal pleiotropy and the juvenile hormone regulation of Drosophila development and life history.
2005,
Pubmed
Flicek,
Ensembl 2011.
2011,
Pubmed
Gentleman,
Bioconductor: open software development for computational biology and bioinformatics.
2004,
Pubmed
Haas,
Phylogeny of frogs as inferred from primarily larval characters (Amphibia:Anura).
2003,
Pubmed
Haas,
De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis.
2013,
Pubmed
Hansen,
STABILIZING SELECTION AND THE COMPARATIVE ANALYSIS OF ADAPTATION.
1997,
Pubmed
Harmon,
GEIGER: investigating evolutionary radiations.
2008,
Pubmed
Harmon,
Tempo and mode of evolutionary radiation in iguanian lizards.
2003,
Pubmed
Hunter,
InterPro: the integrative protein signature database.
2009,
Pubmed
Jønsson,
Ecological and evolutionary determinants for the adaptive radiation of the Madagascan vangas.
2012,
Pubmed
Kanehisa,
KEGG as a reference resource for gene and protein annotation.
2016,
Pubmed
Karpinka,
Xenbase, the Xenopus model organism database; new virtualized system, data types and genomes.
2015,
Pubmed
,
Xenbase
Martin,
Adaptive radiation and behaviour of the Malagasy lemurs.
1972,
Pubmed
Moen,
Testing Convergence Versus History: Convergence Dominates Phenotypic Evolution for over 150 Million Years in Frogs.
2016,
Pubmed
Moore,
Explosive adaptive radiation and extreme phenotypic diversity within ant-nest beetles.
2014,
Pubmed
Morlon,
Inferring the dynamics of diversification: a coalescent approach.
2010,
Pubmed
Pagel,
Bayesian estimation of ancestral character states on phylogenies.
2004,
Pubmed
Paradis,
APE: Analyses of Phylogenetics and Evolution in R language.
2004,
Pubmed
Poux,
Molecular phylogeny and divergence times of Malagasy tenrecs: influence of data partitioning and taxon sampling on dating analyses.
2008,
Pubmed
Pryor,
Tadpole nutritional ecology and digestive physiology: Implications for captive rearing of larval anurans.
2014,
Pubmed
Reddy,
Diversification and the adaptive radiation of the vangas of Madagascar.
2012,
Pubmed
Reimand,
g:Profiler--a web-based toolset for functional profiling of gene lists from large-scale experiments.
2007,
Pubmed
Reimand,
g:Profiler--a web server for functional interpretation of gene lists (2011 update).
2011,
Pubmed
Roelants,
Anuran radiations and the evolution of tadpole morphospace.
2011,
Pubmed
Samonds,
Imperfect isolation: factors and filters shaping Madagascar's extant vertebrate fauna.
2013,
Pubmed
Scantlebury,
Diversification rates have declined in the Malagasy herpetofauna.
2013,
Pubmed
Shubin,
Fossils, genes and the evolution of animal limbs.
1997,
Pubmed
Strauß,
Opposing Patterns of Seasonal Change in Functional and Phylogenetic Diversity of Tadpole Assemblages.
2016,
Pubmed
Szklarczyk,
STRING v10: protein-protein interaction networks, integrated over the tree of life.
2015,
Pubmed
Tarazona,
Data quality aware analysis of differential expression in RNA-seq with NOISeq R/Bioc package.
2015,
Pubmed
Wirta,
An old adaptive radiation of forest dung beetles in Madagascar.
2008,
Pubmed
Wollenberg,
Speciation in little: the role of range and body size in the diversification of Malagasy mantellid frogs.
2011,
Pubmed