Ontogenetic development of intestinal length and relationships to diet in an Australasian fish family (Terapontidae)

Abstract

Background

One of the most widely accepted ecomorphological relationships in vertebrates is the negative correlation between intestinal length and proportion of animal prey in diet. While many fish groups exhibit this general pattern, other clades demonstrate minimal, and in some cases contrasting, associations between diet and intestinal length. Moreover, this relationship and its evolutionary derivation have received little attention from a phylogenetic perspective. This study documents the phylogenetic development of intestinal length variability, and resultant correlation with dietary habits, within a molecular phylogeny of 28 species of terapontid fishes. The Terapontidae (grunters), an ancestrally euryhaline-marine group, is the most trophically diverse of Australia’s freshwater fish families, with widespread shifts away from animal-prey-dominated diets occurring since their invasion of fresh waters.

Results

Description of ontogenetic development of intestinal complexity of terapontid fishes, in combination with ancestral character state reconstruction, demonstrated that complex intestinal looping (convolution) has evolved independently on multiple occasions within the family. This modification of ontogenetic development drives much of the associated interspecific variability in intestinal length evident in terapontids. Phylogenetically informed comparative analyses (phylogenetic independent contrasts) showed that the interspecific differences in intestinal length resulting from these ontogenetic developmental mechanisms explained ~65% of the variability in the proportion of animal material in terapontid diets.

Conclusions

The ontogenetic development of intestinal complexity appears to represent an important functional innovation underlying the extensive trophic differentiation seen in Australia’s freshwater terapontids, specifically facilitating the pronounced shifts away from carnivorous (including invertebrates and vertebrates) diets evident across the family. The capacity to modify intestinal morphology and physiology may also be an important facilitator of trophic diversification during other phyletic radiations.

Background

Morphological divergence associated with dietary shifts has played a major role in the phyletic radiation of many vertebrates [1–5]. These evolutionary changes in diet and trophic morphology can occur rapidly [6, 7], even within ecological timescales [8]. However, the frequency with which particular dietary modes have evolved varies considerably across different vertebrate lineages. While plant-based diets have a broad taxonomic distribution among mammals (>25%) [9], the occurrence of herbivory is much more restricted (2–5% of species) amongst other vertebrate groups [6, 10]. Despite the wide array of feeding modes amongst fishes and the biomass dominance of herbivorous and detritivorous fishes in many assemblages [11, 12], the development of these non-animal prey based trophic habits has been an infrequent evolutionary phenomenon, largely confined to a few families of teleosts [10, 13–16]. The morphological and physiological specializations that facilitate these trophic shifts have accordingly attracted considerable interest from ecologists and evolutionary biologists [10, 16–19].

One of the most widely identified ecomorphological relationships between vertebrate morphology and ecology, and one particularly relevant to dietary radiations involving shifts from carnivory to plant-detrital diets, is intestinal length. The vertebrate digestive tract represents a functional link between foraging (energy intake) and energy management and allocation, but is energetically costly to maintain, and may account for 20–25% of an animal’s metabolic rate [20]. A core prediction of digestive theory [sensu [18, 21] is that the consumption of food with a high content of indigestible material results in an increase in gut dimensions. Numerous studies have shown that digestive tracts across all vertebrate classes tend to be shortest in carnivores, intermediate in omnivores and longest in herbivores and detritivorous species, [20, 22, 23]. The functional significance of this association lies in the need for species on diets that are low in protein and high in roughage to have longer guts in order to ingest large volumes of low-quality food, increase absorptive surface area and maximise digestive efficiency [19]. While a range of fish families display this diet-morphology relationship [17, 24–29], other fish groups demonstrate minimal, and in some cases contrasting relationships between intestinal length and diet [30, 31].

Much of the literature on diet-intestinal length relationships makes little acknowledgment of the evolutionary history of the studied species [18]. Species sharing a common ancestor are not evolutionarily independent, and phylogenetic proximity voids the assumption of sample independence underpinning many conventional statistical tests, thereby creating difficulties in attributing morphological-ecological relationships to adaptive causes rather than phylogenetic artefacts [32]. Applying caution to inferences drawn from phylogenetically naive diet-intestinal length studies is being increasingly advocated [18, 24, 25, 33]. While an abundance of comparative ecomorphological studies of oral kinematics, food procurement and dietary habits in vertebrates has recently emerged [34–37], the association between diet and intestinal length has received surprisingly little phylogenetically informed attention; although recent exceptions have occurred [25, 26].

While developmental plasticity has long been posited as a driver of the origin and diversification of novel traits [38], study of the evolutionary and developmental processes underpinning interspecific differences in intestinal length has been largely neglected. Interspecific variations in intestinal length between closely related species are largely driven by variations in allometric intestinal growth during ontogeny [39, 40]. Substantive allometric increases in intestinal length typically involve additional intestinal looping or convolution that must be accommodated in the body cavity [41]. Previous research has suggested looping patterns are not random, with an underlying phylogenetic component, so that patterns of development of intestinal looping have been used to reconstruct the phylogenetic systematics of a number of fish lineages [41–43]. Yamaoka [43] noted that use of intestinal complexity as a tool in systematic research involves a ‘two-storey’ structure, with the first storey comprising a qualitative aspect (coiling pattern), and the second storey composed of the quantitative (functional) component of intestinal length. To our knowledge, integration of ontogenetic development patterns with molecular phylogenetic reconstruction and comparative approaches has not been attempted. Concurrent appraisal of the ontogenetic processes producing variation in intestinal length, and the functional significance of these processes (i.e., associations with diet) in a phylogenetic context is similarly lacking.

Northern Australia’s Terapontidae (grunters) offer a promising candidate for examining the relationship between intestinal length and diet, and the phylogenetic context for ontogenetic development of intestinal length. The Terapontidae is one of the most speciose and trophically diverse of Australia’s freshwater fish families, exhibiting carnivorous, omnivorous, herbivorous and detritivorous feeding habits [44]. A genus-level phylogeny of the family by Vari [45] relied heavily on differences in ontogenetic development of intestinal configuration as a diagnostic character. Vari’s morphological character analysis suggested a sequence of four intestinal patterns of increasing complexity, beginning at the plesiomorphic condition of a simple two-loop intestine throughout life history in the genera Leiopotherapon, Amniataba, Hannia, Variichthys, Lagusia, Terapon, Pelates, Pelsartia, Rhyncopelates and Mesopristes (Figure 1). The genera Hephaestus, Bidyanus and Scortum have an intermediate “six-loop” pattern. Juveniles of these genera initially possess the “two-loop” pattern before undergoing an ontogenetic elongation and folding to produce more complex patterns as adults. Vari noted that this pattern appears to have been secondarily lost in a subunit of Hephaestus referred to as Hephaestus “genus b”. The adult life stages of the genera Pingalla and Syncomistes purportedly undergo a further ontogenetic shift to produce a highly convoluted and elaborate intestinal pattern, with the final and most complex intestinal pattern unique (autoapomorphic) to Syncomistes. Juveniles of the species in Pingalla and Syncomistes possess a similar intestinal convolution to adults of genera exhibiting the adult “six-loop” pattern, with Vari presuming these species pass through the simple “two-loop” pattern earlier in ontogeny.

Figure 1

Cladogram depicting terapontid generic relationships derived from comparative morphology, adapted from Vari [42]. Cladogram depicting terapontid generic relationships derived from comparative morphology, adapted from [42], showing intestinal convolution and dentition characters used to differentiate genera. Note that Amniataba, Hannia and Variichthys form an unresolved trichotomy. Vari [42] also identified two distinct sub-clades within the genus Hephaestus – “genus a” which develops a “6-loop” intestinal pattern, and “genus b” which retains the plesiomorphic “2-loop” intestine.

A recent molecular-based phylogeny suggests a different topology for this phylogeny, as well as substantial lineage and dietary diversification, particularly the adoption of plant and detritus-based diets, upon a single historical invasion of fresh water environments by euryhaline-marine terapontids [45]. Here we utilise a suite of phylogenetic comparative methods to address two study aims: firstly we re-examine the process of ontogenetic development of intestinal length in Terapontidae within the context of a molecular phylogeny. Patterns of ontogenetic intestinal configuration are described and then combined with ancestral character state reconstruction to examine the evolutionary history of intestinal complexity within terapontids, including the number of gains/losses of particular intestinal patterns within the family. Secondly, in line with predictions of digestive theory, we predict shorter intestinal length in species that consume higher proportions of animal prey than those consuming greater amounts of plant and/or detrital material. If this hypothesized relationship exists, it will provide evidence for dietary ecomorphological diversification, based on modification of intestinal length, which is likely to be a significant driver of the phyletic and trophic radiation evident in Australia’s freshwater terapontids.

Methods

Taxon sampling, molecular markers and phylogeny reconstruction

A phylogenetic analysis of 28 terapontid species was performed based on nuclear and mitochondrial DNA (mtDNA) sequences from Davis and others [45]. The ingroup consisted of 28 species, including nine Australian marine-euryhaline species, all genera and 18 of 24 species of Australian freshwater terapontids, and one species present only in New Guinea. Two representative sequences of one species (Hannia greenwayi) were included due to their different placement in the topology. Distribution of species used in this study in relation to Davis and others [45] are presented in supporting information (see Additional file 1: Figure S1). On the basis of earlier stomach-content based classifications of diet, these selected species exhibit all of the major trophic habits displayed by Australia’s freshwater, euryhaline and marine terapontids: invertivory, generalised carnivory, omnivory, herbivory and detritivory-algivory [45].

Sequence data consisted of an 1141 base-pair (bp) fragment of the mtDNA gene cytochrome b (cytb) and a 3896 and 905 bp fragment of the nuclear recombination activation genes RAG1 and RAG2 (hereafter referred to as RAG) respectively for a total of 5942 bp for each individual included in our study. We used our previous dataset [45]; Dryad Digital Repository doi:http://10.5061/dryad.4r7b7hg1, trimmed out taxa for which we lacked ontogenetic data and realigned the dataset. Cytb was aligned by eye while RAG sequences were aligned using the online version of MAFFT 6.822 [46] using the accurate G-INS-i algorithm with the scoring matrix for nucleotide sequences set to 1PAM/K=2, a gap opening penalty of 1.53 and an offset value of 0.1. Combined partitioned phylogenetic analyses were performed with maximum likelihood (ML) using GARLI 2.0 [47]. We identified the best-fitting model of molecular evolution using the Akaike Information Criterion (AIC) in Modeltest 3.7 [48] using PAUP* 4.0b10 [49]. For cytb Modeltest identified TrN+I+G as the best model and for RAG GTR+I+G was the best model. We ran GARLI with 10 search replicates using the default settings with two partitions representing cytb and RAG with their respective models. For bootstrapping we ran 1000 replicates with the previous settings except that the options genthreshfortopoterm was reduced to 10,000 and treerejectionthreshold was reduced to 20 as suggested in the GARLI manual to speed up bootstrapping. The concatenated sequence data file and tree files were deposited in Dryad, doi:http://10.5061/dryad.h30t5. Trees were rooted with representatives from several related families based on the work of Yagishita and others [45, 50].

Specimen collection

Fish for dietary and morphometric quantification were collected from a number of fish survey studies conducted across fresh water and marine habitats across Australia [44] and Papua New Guinea, as well as being sourced from museum collections. Fish were preserved in either buffered formalin or ethanol. Larger specimens had incisions in the body wall or fixative injected via hypodermic syringe into the body cavity to aid fixation of internal organs.

Intestinal coiling pattern description and intestinal length measurement

After weighing fish and measuring standard length (SL, in mm), specimens were dissected and the entire digestive system and viscera were removed from the body cavity. All terapontids possess a Y-shaped stomach with a straight descending limb from the oesophagus, followed by a blind sac at the bend of the stomach, which leads anteriorly to the pyloric limb on the left side of the body [42]. Intestinal convolution patterns posterior to the pyloric outlet were observed using a dissecting microscope and sketched and photographed from dorsal, ventral, left and right aspects. While Vari [42] described intestinal patterns from the left side of the body, we followed Yamaoka [43] by defining intestinal patterns from the ventral aspect, which facilitates definition of the bilaterally symmetrical body structure of fishes. After description of intestinal coiling structure, the intestine was carefully uncoiled to avoid stretching and intestinal length (IL) was measured as the distance from the pyloric outlet to the rectum. Species’ means for standard length and intestinal length were log₁₀ transformed to homogenise variance prior to analysis and to increase data independence.

Reconstructing the evolutionary history of terapontid intestinal length development

The historical patterns of terapontid intestinal development were hypothesized utilising ancestral character reconstruction techniques in Mesquite 2.75 [51]. We used the “Trace Over Trees” function in Mesquite, which reconstructs ancestral history on multiple phylogenies, to incorporate phylogenetic uncertainty in ancestral reconstructions of character states. In order to generate a collection of trees we used the Bayesian method BEAST 1.7.1 [52] and generated input files using BEAUti 1.7.1. The analysis used an uncorrelated lognormal relaxed molecular clock with rate variation following a tree prior using the speciation birth-death process, and the same models of sequence evolution for the nuclear and mtDNA partitions as per our ML analysis above. BEAST analyses were run for 50 million generations, with parameters logged every 100,000 generations. Multiple runs were conducted to check for stationarity and that independent runs were converging on a similar result. The treefile was summarized using TreeAnnotator 1.7.1 with the mean values placed on the maximum clade credibility tree. The first 10% of trees were removed as burn-in, providing 450 trees for reconstructing ancestral states, with ancestral states summarized onto the maximum clade credibility tree. States were summarized for each node by counting all trees with uniquely best states. If no state was more parsimonious than the other, the reconstruction at that node was classed as equivocal. The frequency of each state was reported for all trees containing that ancestral node, with the variability of inferred states among trees providing a measure of the degree to which ancestral state reconstructions for the node concerned are affected by uncertainty in tree topology and branch lengths. Adult intestinal configurations were coded as discrete (categorical) character states and optimised onto the molecular phylogeny. Because alternative methods of character state reconstruction can produce conflicting results, both maximum parsimony (MP) and maximum likelihood ML methods of ancestral state reconstruction were employed [53, 54]. Parsimony ancestral state reconstruction, which minimizes the amount of character change given a tree topology and character state distribution, has been widely utilised but may over-represent confidence in ancestral character states [53]. For the MP analysis, character transitions were considered to be unordered (changes between any character state are equally costly). A character was assigned to a node if it created fewer steps, otherwise the node was considered equivocal.

ML ancestral character state reconstruction finds the ancestral states that maximize the probability that the observed states would evolve under a stochastic model of evolution [53, 54]. A symmetrical Mk1 model [55], which assumes equal forward and backward character transition rates (i.e., all changes equally probable), was used as the evolutionary model. A major advantage of ML is that the analysis takes branch lengths into account, allows the uncertainty associated with each reconstructed ancestral state to be quantified, and is preferable for medium-sized trees [54, 56]. Likelihood ratios at internal nodes were compared by pairs, and were reported as proportional likelihoods. While likelihoods do not necessarily translate into levels of statistical significance, a difference of 2 log units for a character (i.e., ~7.4 times more probable than any other alternative state) was employed to assign states at a node, otherwise the node was considered equivocal (defined as ‘the rule of thumb’) [54].

Dietary data

Pronounced ontogenetic diet shifts in association with significant allometric growth in many diet-ecomorphological characters are a prominent feature of terapontid ontogenetic biology [40, 44]. To limit any confounding effects of ontogeny on comparative analyses in the present study, assessment focused on the morphologies and dietary habits of the largest size classes only (i.e., when intestinal length was most fully developed). Although the full range of items contributing to the diet of the examined terapontids have been quantified (22 different food classes [40]), in this study, gut contents were simply categorised as the percent contribution of animal material to species’ diet (i.e., the combined contribution of fish, insects and crustaceans). Arcsine transformations of dietary percentages were conducted prior to further analysis to improve normality [57].

Body size-intestinal length correction

Appropriately correcting for body size effects and allometric scaling of morphological traits, while simultaneously taking phylogeny into account, poses an ongoing challenge for comparative studies [58]. To remove effects of body size and allometric scaling of intestinal length between terapontid species, the “phyl_resid” function outlined by Revell [58] was used to regress mean species’ intestinal lengths against mean standard lengths to produce phylogenetically size-corrected residuals in the R package “phytools” [59, 60]. Hereafter, reference to intestinal length refers to the phylogenetically size-corrected estimate.

Testing for phylogenetic signal

To test whether the traits considered in this study (intestinal length and volumetric plant-detrital proportions in diet) individually showed evidence of phylogenetic signal two metrics were utilised – the K statistic [61] and Pagel’s λ[62]. These statistics compare the observed fit of the data to the phylogeny with the analytical expectation based on the topology and branch lengths of the phylogeny, assuming a Brownian (random walk) model of character evolution. Blomberg’s K quantifies the amount of phylogenetic signal in the tip data relative to the expectation (K = 1) for a trait that evolved by Brownian motion along the specified topology and branch lengths [61]. Values of K close to 0 indicate random evolution of traits, values close to 1 correspond to a Brownian-motion-type evolution, and values < 1 indicate strong phylogenetic signal and trait conservatism. Following Blomberg [61], K's significance was assessed using a data randomization test conducted by randomly permutating the tips of the phylogeny 1000 times. A significant phylogenetic signal was indicated if the observed K value was greater than across 95% of the randomizations.

Pagel’s λ provides the best fit of the Brownian motion model to the tip data by means of a maximum likelihood approach [63]. Thus, if λ = 1, the trait evolved according to the Brownian motion, and λ can take any value from 0 (i.e., a star phylogeny, where the trait shows no phylogenetic signal) to >1 (more phylogenetic signal than expected under the Brownian motion). The significance of λ can be assessed by a likelihood ratio comparison of nested models with particular values (i.e., 0 or 1).Tests for phylogenetic signal were implemented using the “phylosignal” and “Kcalc” functions in “phytools” [59]. Both statistics were calculated for traits based on the maximum clade credibility tree.

Phylogenetic comparative analyses

Correlations between intestinal length and dietary composition were examined both with and without phylogenetic correction. To remove the possible correlation associated with phylogenetic relatedness, we calculated phylogenetically independent contrasts (PIC; [32]) of intestinal length and proportion of animal material in species’ diets. For PIC analysis, the molecular topology with branch lengths was imported into Mesquite 2.75 [51]. The PDAP (Phenotypic Diversity Analysis Package) module [64, 65] implemented in Mesquite was used to calculate standardised independent contrasts for the correlation between size-corrected intestinal length and arcsine-transformed proportion of animal material in diet at 28 internal nodes on the terapontid phylogeny. The Pearson product-moment correlation coefficient r (computed through the origin) and its associated P value are reported. The relationship between the phylogenetically independent contrasts was then determined by using a reduced major axis regression (RMA) as there is considerable variation in calculation of both morphological and dietary variables.

Initial diagnostic plots of the absolute values of the standardized phylogenetically independent contrasts versus their standard deviations revealed that branch lengths of the phylogenetic tree adequately fitted the tip data, indicating that estimated branch lengths were adequate for the assumptions of independent contrasts [64]. While PIC is reasonably robust to violations of branch length assumptions [66], additional PICs were calculated using topologies with several arbitrary branch lengths as a sensitivity analysis for any potential uncertainty associated with branch lengths derived in the molecular phylogeny: branch lengths set to unity (1.0 – similar to a speciation model of character evolution), contemporaneous tips with internodes set to one [67], contemporaneous tips with internodes set to one less that the number of descendant tip species [68], and contemporaneous tips with internodes set to the log of number of descendant tip species [68]. All tree manipulation was done using Mesquite (version 2.75).

To assess the effects of failing to control for phylogenetic relatedness, a phylogenetically naive RMA regression (i.e., assuming a star phylogeny) was conducted to investigate the relationship between intestinal length residuals (calculated from an ordinary least squares regression of standard length versus intestinal length) versus arcsine-transformed proportion of animal material in diet.

Results

Phylogenetic analysis

Maximum likelihood recovered one tree with a likelihood score of -34413.698284 (Figure 2). Most nodes within the fresh water radiation in the tree were well resolved with strong support [69] evidenced by bootstrap values mostly >80. Marine-euryhaline species relationships mostly had no bootstrap support. A “highest clade credibility tree” generated from BEAST analyses also had a very similar structure to the ML approach, especially when considering well supported nodes (see Additional file 2: Figure S2), highlighting the general similarities in tree structure regardless of phylogenetic construction approach used.

Figure 2

Maximum likelihood phylogeny for 28 terapontid species based on analysis of combined nuclear and mitochondrial DNA. All bootstrap values are based on 1000 pseudoreplicates. Outgroup species were pruned from the tree. Images are identified by initials of genus and species nearby in the tree. Taxon names are colour-coded according to macrohabitat associations identified in [44]: red = marine, green = euryhaline, and black = freshwater. The node signifying invasion of Australasian freshwater habitat is indicated.

Dietary and morphological quantification

Data on diet, length and intestinal length of each species and trophic classifications from the literature are presented in Table 1. There is broad variability in terapontid diets, from those comprising only animal prey to those consuming almost no animal material. Presenting the often flexible dietary habits of fishes as averages can mask important spatio-temporal variability. Previous studies on the trophic ecology of terapontids (including specimens used in this study) identified low spatial and temporal dietary variability in many species, particularly herbivore-detritivores [70]. Dietary data derived in this study also agreed closely with often more seasonally and spatially comprehensive data from other sources and geographic areas [45, 71, 72]. This suggests available data should provide a robust approximation of typical diet for most species, particularly at the broad level of comparative amounts of animal prey in diet.

Table 1 Summary data on terapontid morphology and diet used in study

Relative intestinal length (IL/SL) is the most commonly used descriptor in diet-morphology assessments [17], so RIL ranges are provided for comparison with published data (Table 1). Reduced major axis regressions of log₁₀–transformed standard length versus log₁₀–transformed intestinal length for each species over the available studied size are also outlined in supporting information (see Additional file 3: Table S1).

Ontogenetic development of intestinal morphology

Fish digestive tracts were examined for elaborations such as hindgut chambers, caecal pouches and valves. The only external differences in intestinal structure between terapontid species appear to be length and coiling patterns. The simplest intestinal layout consisted of two loops and was evident immediately after post-larval metamorphosis in all species examined (Figure 3, configuration 1A-1B). The first loop occurred posterior of the pylorus near the rear of the body cavity, after a slight dextral curve immediately posterior to the pylorus. The intestine continued anteriorly until a second loop occurred ventral to the stomach, after which the intestine continued posteriorly to the anus, producing an “s-shaped or two-loop” layout. This simple configuration was evident throughout the life history of Leiopotherapon unicolor, Amniataba percoides, A. caudavittatus, Hannia greenwayi, Hephaestus carbo, Hep. epirrhinos, Hep. transmontanus, Pelates sexlineatus, P. quadrilineatus, Terapon theraps, T. puta, T. jarbua and Varichthys lacustris (see Additional file 2: Figure S3); however, significant allometric increases in intestinal length were achieved in several species by increasing the length of each loop in both anterior and posterior directions (Figure 3, configuration 1B). RILs of <1.2 typify terapontid species that retained the two-loop intestinal configuration throughout their life history (Table 1). The “two-loop” intestinal configuration was the juvenile morphology of all remaining species.

Figure 3

Patterns of ontogenetic development of intestinal layout in terapontids. Intestinal tracts are viewed ventrally, the anterior-most portion of the intestine (outlet of the pylorus) is always located to the top of each figure.

In all other species except L. aheneus and Hel. sexlineatus, a transverse folding and elongation of the middle portion of the two-loop intestinal pattern occurred, directing the elongated section to the left of the body cavity (Figure 3, configuration 1C and 1D) and ventrally beneath the posteriorly directed section of the two-loop pattern. This produced the “six-loop” configuration described by Vari [42], which was the layout throughout the remaining life history of Bidyanus welchi, Hep. fuliginosus, Hep. jenkinsi and Hep. tulliensis (see Additional file 2: Figure S4). Adult RILs of ~2 to 2.5 characterized these species (Table 1). In Pingalla and Scortum species, the loops on the right-hand side of the body cavity continued dorso-anteriorly before turning to lengthen in a posterior direction (Figure 3, configuration 1E-1F). In Pingalla species this remained the intestinal layout of adults. A further increase in intestinal complexity occurred in Scortum species and was characterised by additional convolution in a spiral configuration (Figure 3, configuration 1G-1H). In all of these species the majority of convolution occurred on the left-hand side of the body cavity. The RILs of Scortum species averaged ~4.5, and reached over 7 in some specimens (Table 1; Additional file 2: Figure S7).

A different development of intestinal configuration was evident in Syncomistes species. Like Hephaestus, Pingalla and Scortum species, Syncomistes species developed the “six-loop” pattern, but the subsequent looping in Syncomistes proceeded anteriorly before folding and lengthening to the right-hand side of the body cavity. At the same time, posterior looping from the “six-loop” configuration proceeded to the left-hand side of the body cavity behind the stomach (Figure 3, configuration 2C-2E). This was followed by a reversal of looping directions in both the anterior and posterior sections, looping back to the left- and right-hand side of the body cavity respectively (Figure 3, configuration 2F- 2I). This complex intestinal configuration resulted in RILs of Syncomistes reaching over 6 in some specimens (Table 1; Additional file 2: Figure S10).

Another distinct pattern of ontogenetic intestinal looping was evident in Leiopotherapon aheneus. From the initial two-loop pattern the anterior loop lengthened anteriorly along the ventral surface of the stomach close to the pyloric outlet (Figure 3, configuration 3A-3B). This was followed by a folding in the middle section of the intestine (Figure 3, configuration 3C-3E). This folding initially proceeded anteriorly along the dorso-ventral plane of the body before turning to the right-hand side of the body cavity (Figure 3, configuration 3F-3G). The majority of folding in this pattern occurred on the right-hand side of the body. Intestinal lengths of L. aheneus typically reached between 2-3 times standard length in larger specimens (Table 1; Additional file 2: Figure S8).

A final distinct pattern of ontogenetic intestinal looping was evident in Hel. sexlineatus. From the initial two-loop pattern, the posterior and anterior loops extended in both directions during ontogeny. The anterior loop then extended past the pyloric outlet, before looping around the anterior aspect of the stomach, crossing the dorso-ventral plane to lengthen into the anterior, right-hand side of the body cavity (Figure 3, configuration 4D-4E). While only a comparatively minor increase in complexity, this configuration produced higher RILs compared to the standard “two-loop” intestinal layout (Table 1; Additional file 2: Figure S11).

Character optimisations and reconstruction of ancestral character states

Optimising adult intestinal configuration patterns onto the maximum clade credibility phylogeny indicated that the ontogenetic development of increased intestinal complexity has evolved independently on three occasions in terapontid fishes. While a range of patterns of ontogenetic increase in intestinal complexity have evolved in the clade containing Hephaestus, Scortum, Bidyanus, Syncomistes and Pingalla species, ontogenetic increases in intestinal convolution were limited to just a single species (L. aheneus) in the other major freshwater clade, as well as on a single occasion in the euryhaline/marine clade (Hel. sexlineatus). An examination of ancestral state reconstructions across the 450 trees from the BEAST analysis yielded very similar predictions between parsimony and likelihood analyses (Figure 4) and the inferred ancestral states for terapontid intestinal length configuration were not substantially affected by uncertainty in tree topology, branch lengths, or character state reconstruction methods. Both MP and ML analysis indicated that the “two-loop” pattern is unequivocally plesiomorphic within Terapontidae, and that the “two-loop” intestinal pattern was exhibited by the most recent common ancestor of all freshwater species (i.e., at the time of fresh water invasion). Both reconstruction approaches also indicated that the evolution of adult intestinal complexity followed a complex pattern of multiple independent gains and one loss within both major freshwater clades. Both approaches indicated that the “six-loop” intestinal configuration was a precursor to the range of more complex intestinal patterns evident in Pingalla, Scortum and Syncomistes species. Character state reconstruction also suggested that the two similar patterns of increase evident in Pingalla and Scortum species evolved independently. An apparent reversion to the plesiomorphic state of an adult “two-loop” intestinal pattern was also evident in Hep. epirrhinos, the only species within this clade to retain this intestinal configuration as an adult.

Figure 4

Ancestral character state reconstruction of terapontid intestinal morphology. Summary of maximum likelihood (left graph) and maximum parsimony (right) ancestral character reconstruction of adult intestinal configuration for 450 terapontid trees displayed on the maximum clade credibility tree. Circles at terminal nodes represent the observed character state for extant species. Pie charts for ancestral nodes show estimated proportions for reconstructed character states at that internal node.

Phylogenetic signal

Blomberg’s K and Pagel’s λ for proportion of animal prey in diet and intestinal length both demonstrated significant levels of phylogenetic signal, indicating that neither variable was independent and, therefore, phylogenetic comparative methods were justified in further analyses. While the estimates of phylogenetic signal for the two variables were both significant, the patterns of phylogenetic signal were not convergent. Phylogeny was a significant predictor of variation in animal material in terapontid diet (K = 0.73, observed PIC variance = 1.01, P < 0.001, Pagel’s λ = 0.88, P < 0.001). However, both K and λ were estimated to be considerably less than 1, suggesting a phylogenetic signal lower than the one expected under Brownian motion and, accordingly, substantial evolutionary lability in terapontid diet, even between closely related species. Phylogeny accounted for a larger component of variability in intestinal length in the terapontids (K = 1.05, observed PIC variance = 0.294, P < 0.001, Pagel’s λ = 0.94, P < 0.001), suggesting a phylogenetic signal close to what would be expected under Brownian motion in both statistics.

Comparative analyses

After correcting for phylogenetic proximity, the independent contrasts of intestinal length versus diet were significantly, and negatively, correlated with the percentage of animal material in terapontid diet, explaining 65% of variation in diet composition (r²= 0.65, RMA slope = -1.48, P < 0.001). Twenty-two of the 28 independent contrasts were negative, and occurred across both deep and shallow nodes of the phylogeny (Figure 5). Several of the most negative contrasts occurred at nodes within the phylogeny (nodes 38, 55, 24, 29, 7, 18 and 19) that were precursors to gains/losses in intestinal complexity identified in the character mapping and ancestral character reconstruction (Figure 4). This highlights the importance of gains in intestinal complexity in facilitating dietary radiation. PICs with branch lengths set to unity (r² = 0.46, RMA slope = -1.94, P < 0.001), Nee in Purvis [68] branch lengths (r² = 0.52, RMA slope = -1.80, P < 0.001), Grafen [73] branch lengths (r² = 0.54, RMA slope = -1.71, P < 0.001) and Pagel’s [67] branch lengths (r² = 0.49, RMA slope = -1.82, P < 0.001) all produced similar results to the molecular phylogeny. A phylogenetically naive RMA regression also identified a significant negative relationship between intestinal length residuals and arcsine-transformed proportion of animal material in diet, and, furthermore, this analysis explained a greater proportion of data variation than any of the phylogenetic comparative analyses (r² = 0.72, RMA slope = -1.50, P < 0.001).

Figure 5

Phylogenetic independent contrasts of terapontid diet versus intestinal length. (a) Relationship between phylogenetically independent contrasts of intestinal length residuals and contrasts of arcsine transformed proportion of animal material in diet. Numbers represent the nodes (contrasts) indicated in the phylogeny in (b).

Discussion

Evolution of intestinal length and dietary radiation in terapontids

Several patterns of ontogenetic development of increased intestinal length were evident in the terapontid species examined. Like previous studies [41–43], results highlighted an underlying phylogenetic component to these developmental patterns. The interspecific differences in intestinal length resulting from these ontogenetic developmental mechanisms explained a significant amount of the variability in the volume of animal prey in terapontid diets. Results indicate that the widely held ecomorphological maxim of longer digestive tracts equating with increasing consumption of non-animal prey, holds true for terapontids, even when accounting for phylogenetic relationships between species. Study outcomes align with a growing number of studies, where if phylogeny is taken into account, carnivores have shorter intestines than related species consuming larger amounts on non-animal prey [25, 26].

This study produced a number of commonalities as well as contrasts to the previous work on the family outlined in Vari [42]. Both studies identified the “two-loop” intestinal configuration as being the plesiomorphic adult pattern within Terapontidae. This study suggested a number of different contrasts to the patterns of intestinal development across the family, at both species and family levels. The secondary loss of the “six-loop” intestinal layout Vari [42] proposed in “Hephaestus genus b” instead appears due to the polyphyly of Hephaestus and phylogenetic location of this “Hephaestus genus b” in a separate clade of species with a “two-loop” intestinal layout. Vari [42] suggested that Scortum species shared the same adult “six-loop” intestinal pattern as Bidyanus and Hephaestus species (Figure 1). The current studyinstead highlighted Scortum and Syncomistes species as developing the most complex intestinal patterns of any terapontid species. This study also identified previously undescribed pattern of ontogenetic intestinal length increase in L. aheneus and Hel. octolineatus. The different topology emerging from molecular relationships compared to Vari’s [42] phylogeny also suggested a different sequence of intestinal length complexity across the family. Rather than the progressive increase in complexity as genera become more derived, proposed by Vari [42] (Figure 1), a more complex historical process of development was predicted from molecular relationships. Character-state reconstruction inferred that the relatively complex intestinal configurations of adult Pingalla, Scortum and Syncomistes species all evolved from the six-loop pattern on three separate occasions. The novel ontogenetic development documented in both L. aheneus and Hel. octolineatus also demonstrated that the capacity for significantly increasing intestinal length during ontogeny has evolved independently in both major clades of freshwater terapontids as well as euryhaline-marine species. These multiple independent origins of increased intestinal complexity across several clades suggests convergent evolution toward increased intestinal length in terapontids having diets with lower proportions of animal prey.

Although the role of ontogenetic phenomena in phyletic evolution remains strongly debated [74–76], modification of ontogenetic development is proposed as one of the most common mechanisms through which morphological change and novelties originate during phyletic evolution [74, 75]. The development of intestinal complexity in terapontids exhibits several elements of heterochronic processes [74, 75], where ontogeny is modified to produce morphological novelty. Several possible peramorphic (recapitulatory) processes, for example, could explain the apparent repetition of adult intestinal layouts (two-loop and six-loop patterns) of ancestral forms during the ontogeny of many descendent terapontid taxa, before additional intestinal complexity is added to ancestral configurations. A range of associated heterochronic processes (acceleration, hypermorphosis and pre-displacement) can all produce descendent phenotypes that transcend the ancestral form [74, 75]. Similarly, paedomorphic phenomena, where adults retain the juvenile morphology of putative ancestral taxa, could explain the apparent retention of two-loop intestinal layout throughout the life history of Hep. epirrhinos, within a clade of closely related Hephaestus species that have a six-loop configuration (Figure 4). Without a range of additional size/age and possibly shape-based data on terapontid ontogenetic trajectories [71, 77], the exact role of heterochronic processes can only be speculated upon. However, recapitulation does appear to be a recurrent theme in the development of intestinal length complexity in a number of fish lineages [43]. With additional genetic and morphological data, terapontids may provide a valuable model lineage for elucidating the role of modification of ontogeny as a driver of evolutionary diversification.

The utility of intestinal length as a predictor of diet

While standard regression and PIC approaches both highlighted significant relationships between intestinal length and animal material in the diet, the amount of variability explained was lower in the PIC analysis. This difference underlines the importance of comparative methods in not overstating the strength of the association between morphology and ecology [27]. Although intestinal length emerged from the phylogenetically informed analysis as a useful predictor of diet, a substantial amount of unexplained variability was also evident in the relationship. Behavioural, ecological, physiological and historical factors can interact to influence the strength of the congruence between morphological and ecological characters [78]. Issues associated with age, phenotypic plasticity, antecedent food availability (i.e., periods of starvation) as well as the relative levels of different dietary substrates have emerged from both field and controlled laboratory studies as possibly inducing changes in intestinal length [17, 33, 79, 80]. While intestinal length is clearly a somewhat plastic character, ontogenetic and phylogenetic factors appear more influential than diet on gut dimensions in some fish clades [33], suggesting a precedence of genetic adaptation over phenotypic plasticity as the major force acting on the digestive system. Intestinal looping patterns identified in this study were largely consistent with previous research (see Vari; [42]), and seemed species/genus-specific, but strength of any underlying genetic component to their expression needs to be tested with controlled feeding experiments e.g., [33, 79]. The capacity for at least some phenotypic plasticity in intestinal length in response to different trophic opportunities could promote initial divergence in dietary habits, and potentially provide scope for natural selection to extend and consolidate the phenotypic response.

While intestinal length may be a useful predictor of broad dietary habits, it may have a variable capacity to predict finer scale dietary divisions among omnivores [28]. Many of the terapontids examined here consume varying proportions of both animal and plant or detrital material, and would require more robust dietary and morphological data to adequately test ecomorphological relationships at a finer scale. Omnivory has been interpreted as a compromise strategy in which protein from scarce animal prey is complemented by energy from abundant primary foods [81]. Omnivory and generalist diets are also regarded as an adaptive response to seasonal variations in water level and trophic resources that characterise hydrologically variable tropical river systems [82]. With the wet-dry tropical catchments that harbour the majority of terapontid diversity ranking among some of the most hydrologically variable globally, versatility in feeding habits is, not unexpectedly, a common feature of many terapontid diets [44, 45].

We also used intestinal length as a dietary predictor in relation to stomach content data. Classifying diets on the basis of stomach content analysis can be problematic for fishes, particularly nominal herbivores and detritivores, dietary habits expressed by several species in this study (at least on the basis of stomach content data). Conventional macroscopic dietary quantification can be prone to inadequately identifying the actual nutritional targets of ingestion, and often require integration with microscopic, histological or stable isotopic approaches to accurately define dietary ecology. Many marine ‘herbivores’ once commonly perceived to be algivores have been revealed by detailed dietary analyses to be highly dependent on amorphic detritus scraped from epilithic algal complexes [16, 83, 84]. Similarly, recent studies have indicated that freshwater ‘detritivorous’ fishes assimilate carbon from biofilm and seston, and nitrogen from intermediate microbial decomposers in the environment, and are not capable of direct assimilation of vascular plant carbon [16]. In contrast to the abundant research on terrestrial vertebrate ‘nutritional ecology’, the nutritional targets, food composition and associated digestive functioning of herbivorous-detritivorous fish are poorly defined [10, 18, 84]. While these gaps are being addressed in the marine environment (incrementally in some areas; [84]), they are equally, if not more pronounced in freshwater species, and pose a considerable impediment to understanding the trophic ecology and food web function of herbivorous-detritivorous freshwater fishes [16, 85, 86].

Alimentary anatomy is frequently an unreliable indicator of functional capacity of herbivorous fishes, particularly if the digestive tract is considered in isolation [84]. Morphological and functional changes to the biomechanics and musculoskeletal functional morphology related to food procurement and handling are considered critical components in the impressive evolutionary diversification and ecological success of teleosts, including many herbivorous and detritivorous fishes [87, 88]. There are marked changes in oral anatomy (flattened, depressible dentition, dentary rotation etc.) across several of the freshwater genera within the Terapontidae, such as Scortum, Pingalla and Syncomistes species, that have adopted diets volumetrically dominated by plant and/or detrital material Figure 1; [42]. Interestingly, the marine herbivore Hel. sexlineatus, recently separated from the Pelates genus [89], also appears to have evolved flattened, tricuspidite dentition similar to that of freshwater herbivores [90]. Assessment of these changes to oral anatomy and feeding kinematics in relation to terapontid trophic diversification would be a valuable complement to the role of intestinal modification documented in this study.

Intestinal length considered in isolation is also in many ways a simplistic indicator of the functional morphology of fish intestinal tracts. Other aspects of digestive morphology and physiology such as intestinal diameter, digesta passage rates, ultra-structural surface area and digestive enzyme profiles can also have significant associations with diet [24, 25, 91–94]. Beyond the significant correlation between increasing intestinal length and decreasing animal prey in diet, interpretation of the evolution of specific dietary habits in the Terapontidae must be treated with caution as we currently have limited insights into the physiological processes associated with extracting and utilizing nutrients from consumed foods. A nutritional ecological approach, however, incorporating knowledge of diet, functional morphology, intake, digestive physiology and dietary assimilation [sensu 18,86] would provide a more robust foundation with which to resolve the trophic habits of terapontids. Regardless of the underlying nutritional targets and associated digestive mechanisms, however, the pronounced shifts toward non-animal prey evident in many terapontids are clearly associated with significant modification of intestinal length.

Terapontids as a model system for studying dietary diversification

The capacity to increase intestinal length, and associated shifts away from carnivory, have evolved independently across multiple marine-euryhaline and freshwater genera within Terapontidae, but are especially pronounced in freshwater species. Shifts away from carnivory and evolution of herbivory and plant-detrital diets are prominent in many of the more speciose and ecologically diverse marine and freshwater fish lineages, often marking a profound shift in the phylogenetic trajectories, species diversity and ecological impact of certain clades [87, 95]. The significant diet-intestinal length relationship evident in approximately 55% of extant terapontid species suggests that the capacity to develop long intestines during ontogeny has facilitated the widespread shifts away from carnivorous diets across the family.

Studies of trophically diverse lineages using cladistics and assessment of digestive tract characters could be useful in elucidating the process of evolution of herbivorous-detritivorous trophic habits [17]. Terapontids could provide such a model to demonstrate the process of evolution of non-animal prey-based diets from an ancestrally carnivorous lineage. With carnivory the likely ancestral habit of the euryhaline-marine ancestors of Australia’s freshwater terapontids, the invasion of fresh waters saw adoption of a variety of omnivorous, herbivorous and detritivorous dietary habits during the terapontid fresh water radiation [45]. Due to its biogeographic isolation the Australian freshwater fish fauna is particularly unusual for its prevalence of acanthopterygian percomorph fishes (which typically dominate marine habitats), and for an almost complete lack of ostariophysan fishes which dominate fresh water environments on other continents. The timing of Australia’s break-up from Gondwana precluded the presence of cichlids, characiforms, cypriniformes and most siluriformes, which represent the dominant proportion of herbivores and detritivores in other continents’ freshwater fish faunas [95, 96]. The majority of Australia’s freshwater fishes are ‘secondary’ freshwater teleostean species (i.e., freshwater species derived from marine ancestors), often with strong affinities to tropical Indo-Pacific marine taxa [97, 98].

Fossil evidence suggests that the Terapontidae has had a long evolutionary history (≥ 40-45 Ma) in Australian fresh waters [99]. Paleoecological conditions that may have facilitated the dietary diversification of early fresh water-invading terapontids, particularly shifts away from carnivory, probably include a range of vacant niches due to a lack of an incumbent herbivorous-detritivorous fish fauna [45]. Similar processes and timescales relating to ecological opportunity and release from competitive constraints have been proposed to explain the significant morphological disparification and lineage diversification evident in Australasian ariid catfishes following a similar fresh water invasion [100]. Following invasion of a new habitat, species may show a rapid burst of cladogenesis and associated ecomorphological (often diet-related) diversification [3, 101–103]. The majority of morphological divergence in characters like intestinal convolution and dentition appear to have occurred independently on several occasions in freshwater terapontids [43]; this study. The significant relationship between intestinal length and shifts away from animal prey in the diet of terapontids suggests that the evolution of longer intestines, in particular, facilitated much of the dietary diversification evident in Australian fresh water environments.

Conclusions

Intestinal length is a significant correlate to interspecific dietary variation in terapontids. The ontogenetic development of intestinal complexity appears to represent an important functional innovation driving much of the ecological (trophic) radiation evident within Terapontidae. The significant negative correlation between trophic morphology (intestinal length) and proportion of animal material in terapontid diet suggests resource-based divergent selection as an important diversifying force in the adaptive radiation of Australia’s freshwater terapontids, particularly the pronounced shifts away from ancestral carnivorous dietary habits evident across the family. Much previous research has suggested that modifications of oral anatomy and functional associations with initial food procurement are one of the primary drivers of fish lineage diversification [36, 37, 104, 105]. The capacity to modify intestinal morphology-physiology in light of new digestive challenges may also be an important facilitator of trophic diversification during phyletic radiations see also [8, 26, 106]. Moreover, the ontogenetic development of a range of intestinal convolutions being limited to freshwater terapontids is suggestive of ecomorphological character release within the family following invasion of fresh waters by ancestral euryhaline-marine species. Assessment of the relative patterns of lineage diversification between freshwater and euryhaline-marine terapontids in other aspects of trophic morphology sensu [100] and ecology would be fruitful avenues for research on the phylogenetic effects of adaptive zone shifts.