A New 13 Million Year Old Gavialoid Crocodylian from Proto-Amazonian Mega-Wetlands Reveals Parallel Evolutionary Trends in Skull Shape Linked to Longirostry

Gavialoid crocodylians are the archetypal longirostrine archosaurs and, as such, understanding their patterns of evolution is fundamental to recognizing cranial rearrangements and reconstructing adaptive pathways associated with elongation of the rostrum (longirostry). The living Indian gharial Gavialis gangeticus is the sole survivor of the group, thus providing unique evidence on the distinctive biology of its fossil kin. Yet phylogenetic relationships and evolutionary ecology spanning ~70 million-years of longirostrine crocodylian diversification remain unclear. Analysis of cranial anatomy of a new proto-Amazonian gavialoid, Gryposuchus pachakamue sp. nov., from the Miocene lakes and swamps of the Pebas Mega-Wetland System reveals that acquisition of both widely separated and protruding eyes (telescoped orbits) and riverine ecology within South American and Indian gavialoids is the result of parallel evolution. Phylogenetic and morphometric analyses show that, in association with longirostry, circumorbital bone configuration can evolve rapidly for coping with trends in environmental conditions and may reflect shifts in feeding strategy. Our results support a long-term radiation of the South American forms, with taxa occupying either extreme of the gavialoid morphospace showing preferences for coastal marine versus fluvial environments. The early biogeographic history of South American gavialoids was strongly linked to the northward drainage system connecting proto-Amazonian wetlands to the Caribbean region.


Introduction
Numerous unresolved issues hinder understanding of the origin, time of divergence, and patterns of adaptive radiation of gavialoid crocodylians. Whereas molecular data sets favor a close relationship and an Eocene [1] or even Neogene [2][3][4] divergence between the Indian gharial Gavialis gangeticus and the Indonesian false gharial Tomistoma schlegelii, morphological phylogenies suggest that: 1) these two extant longirostrine species are much more distantly related, 2) their elongated skulls evolved convergently, and 3) the oldest fossil gavialoid dates back to the Cretaceous Period in North America, Europe, and Africa [5,6]. In fact, analyses of extant Gavialis and its nearest fossil relatives do not provide strong support for their phylogenetic affinities with any other crocodylian clade, probably as part of what Clark [7] called the "longirostrine problem". Clark [7] identified several cranial features related to longirostry in crocodyliforms that might have appeared independently in different groups, including possessing widely separated and protruding eyes, namely "telescoped" orbits. Potential convergent evolution of this distinctive cranial morphotype in gavialoids not only obscured the ancestral anatomical condition for the entire clade, but also created a challenge for deciphering ingroup affinities. For example, Indian Gavialis and South American Gryposuchus species are extremely similar in cranial morphology and both have fully "telescoped" orbits [8,9]. Whether this distinctive pattern results from common or independent origin has been uncertain and debated [6,10]. Here, we describe a gavialoid with "non-telescoped" orbits from the Middle Miocene of the Pebas Formation of northeastern Peru, that provides evidence of early morphological stages of the evolution of the Gryposuchus lineage in the Amazonian Neotropics relevant to resolving shared or independent origin of "telescoped" orbits. This new gavialoid is the only longirostrine species in the hyperdiverse crocodylian community that inhabited lakes, swamps, and deltas of the Pebas Mega-Wetland System [11], a huge proto-Amazonian biome that appears to have played a crucial role for marine to freshwater transitions in many vertebrate groups occurring today in various river drainage systems of tropical South America [12]. Its lineage survived within the paleo-Orinoco drainage throughout the Late Miocene, providing further evidence for the persistence of Pebas-like mega-wetland conditions in aquatic environments of northernmost South America [11]. These records highlight the biogeographic role of the long-lasting drainage linking western proto-Amazonia with the Caribbean region prior to onset of the eastward flowing, transcontinental Amazonian River system approximately 10.5 million years ago. We analyze the phylogenetic relationships of this new Pebasian species to test whether parallel evolution occurred within an adaptive radiation of Indian and South American gavialoids. Mapping our phylogenetic hypothesis onto a morphometric space, the pattern of the remarkable taxonomic and anatomical diversification of South American forms provides novel insights into the ecological significance underlying circumorbital skull configurations throughout gavialoid history.

Materials and Methods
The Natural History Museum of National University of San Marcos (MUSM), Lima, Peru, provided all permits for this study and houses permanently the specimens from the Pebas Formation described in this publication (Credential Number MUSM 2008-8). All data has been collected, prepared and deposited in this public national Peruvian repository by the authors of the article. The current data is based only in fossil specimens, to which it has been assigned a specific catalog number. Specimen

Phylogenetic analysis
To determine the phylogenetic relationships of the new Pebasian gavialoid species, we included it in a data matrix of morphological characters revised and updated from Salas-Gismondi et al. [11], that built on characters initially developed by Brochu [13,14] and Jouve [6], as well as additional characters compiled from other contributions (S1 Appendix). The present analysis focuses on testing the phylogenetic relationships among gavialoid crocodylians rather than more broadly across all Crocodylia as in the former analyses. The complete data matrix consists of 206 morphological characters for 42 eusuchian taxa, with Bernissartia fagesii as an outgroup and including most members of the gavialoid clade but only representative taxa of Brevirostres (S1 Appendix). This matrix was analyzed with maximum parsimony methods using TNT 1.1 [15]. All characters had equal weighting and were treated as unordered and non-additive. To assess nodal support, branches with a minimum length of 0 were collapsed and Bremer support values (decay indices) were calculated and shown on the strict consensus phylogeny.

Cranial circumorbital morphospace analysis and phylogenetic mapping
In order to understand morphological evolution linked to longirostry, we applied a geometric morphometric approach to analysis of the circumorbital anatomy of 22 species of crocodylians, including extant and extinct species (S1 Appendix). As the circumorbital region has shown only minor intraspecific variation within extant adult crocodylians [16], all taxa are represented by one fully adult specimen, each preserving in dorsal view: (1) rostral sutures, (2) orbital shape pattern, and (3) no significant postmortem distortion. Morphometric analyses require complete datasets for all landmarks, thus MUSM 1981 was partially reconstructed from the morphology of the complete left side via geometric reflection [17], and included in the analysis. We selected twenty discrete landmark loci based on their availability, relative coplanarity, and clear demarcation within images (S1 Fig). Two-dimensional landmark digitization, Procrustes superimposition, and principal component analyses (PCA) were performed with the Geomorph package in R software [18,19]. The 22 landmark configurations were superimposed and scaled to unit centroid size following the generalized Procrustes method [20]. The coordinates of the superimposed configurations were later imported into the Euclidean tangent shape space for subsequent statistical analysis [20]. In order to determine components of shape variation, the 22 principal components from that analysis were obtained by applying a principal component analysis (PCA) on the variance covariance of the projected coordinates. Principal component axes (PCs) 1 and 2 (estimated cumulative variance = 70.6%) were plotted in principal components space. A simplified version of the parsimony-based phylogenetic hypothesis unweighted for branch length and including only the taxa sampled for the twenty landmarks, was mapped in the circumorbital morphospace with MorphoJ software [21].

Nomenclatural acts
The electronic edition of this article conforms to the requirements of the amended International Code of Zoological Nomenclature, and hence the new names contained herein are available under that Code from the electronic edition of this article. This published work and the nomenclatural acts it contains have been registered in ZooBank, the online registration system for the ICZN. The ZooBank LSIDs (Life Science Identifiers) can be resolved and the associated information viewed through any standard web browser by appending the LSID to the prefix "http://zoobank.org/". The LSID for this publication is: urn:lsid:zoobank.org: pub:31C1EBEF-A7B1-44BA-ACB4-1552AD2BC650. The electronic edition of this work was published in a journal with an ISSN, and has been archived and is available from the following digital repositories: PubMed Central, LOCKSS.

Geological Context, Paleoenvironment, and Age
In the study area, the Pebas Formation is comprised of transgressive and regressive bay-margin deposits. Parasequences are formed by fossiliferous blue to gray clays interbedded with unconsolidated sands, which are typically capped by lignite layers and mollusc shell beds [22][23][24]. Lacustrine brackish-water, sometimes tidally influenced, mega-wetland deposits are indicated by most of the studies in the lower and middle parts of the Pebas Formation. The time span of the Pebas Formation ranges approximately from the early Miocene (ca. 23 Ma) to the early Late Miocene (ca. 10.5 Ma) [25]. The biostratigraphic framework of the Pebas Formation is based on pollen, ostracod, and particularly on the abundant and diverse molluscan faunas preserved throughout the unit [26][27][28]. Fossiliferous deposits are located along the Amazon River banks and tributaries of the so-called Iquitos Arch, which corresponds to the modern forebulge of the northwestern Amazonia foreland basin [24]. Fossil vertebrates are found mainly in lignitic bonebeds [11]. These localities (IQ) are mapped within the Molluscan Zones (MZ) proposed by Wesselingh et al., [28] for the study area (Fig 1). Most of the gavialoid remains come from localities IQ26 and IQ114, and correspond to MZ8 (late Middle Miocene; ca. 13 Ma). These localities have documented the highest diversity within any single known extant or fossil crocodylian community, including the shoveling mollusc-crushing caiman Gnatusuchus pebasensis, as well as Kuttanacaiman iquitosensis, Caiman wannlangstoni, Purussaurus neivensis, Mourasuchus atopus, Paleosuchus sp. [11], and the new gavialoid taxon described here. Other fossil vertebrate remains, such as fishes, aquatic turtles, and mammals, also are abundant. IQ116 is a less productive locality than IQ26 and IQ114, but might correspond to the same time interval. This later outcrop is restricted to a 1-2 meter thick capping lignite level overlying a grey mud with shell beds [11]. Two specimens of the new Pebasian gavialoid, MUSM 1681 and MUSM 1682, were found in locality IQ136, probably corresponding to MZ5 (Middle Miocene; ca. . These are the oldest gavialoids known from the Amazon Basin. An additional gavialoid specimen, MUSM 987 was unearthed at IQ101, in the Momón River banks, and might correspond in age to MZ6 (ca. . Pebas Formation deposits represent a highly dynamic mosaic of shallow aquatic settings (e.g., lakes, swamps, and rivers) of varying salinity and marine influence [29][30][31][32], and dysoxic muddy bottoms within lakes and swamps were common in the Pebas System. Several gavialoid bones, including a diagnostic postorbital bone (MUSM 2430) belonging to Gryposuchus colombianus or Gr. croizati, were found at locality IQ125 in the Nueva Unión area. This locality belongs to the "Uppermost Pebas Formation" (Fig 1) [11,33] and might correspond to MZ9 or younger intervals (Middle-Late Miocene boundary; ca. 12-11 Ma). Outcrops there consist of fine-grained fluvial sandstones, floodplain clays and silts, and paleosols. These levels are lignite-poor and lack mollusc remains [33].
Among the copious gavialoid material reported from the Late Miocene Urumaco Formation of Venezuela, we identified a taxon consistent in morphology with the new Pebasian gavialoid [34]. The Urumaco Formation consists of complex intercalation of sandstone, organic-rich mudstone, coal, shale, and thick-bedded coquinoidal limestones with abundant mollusc fragments [35]. The Urumaco specimen comprises a partial skull (AMU CURS 12; Fig 2D) collected within the upper member of the Urumaco Formation, at the Domo de Agua Blanca Locality [34]. This member is characterized by organic-rich, dark-gray laminated mudstone and shale, and abundant vertebrate fragments [35]. The upper member of the Urumaco Formation was deposited in a delta plain [35]. 3M); MUSM 1727, snout of a juvenile, Locality IQ26 (MZ8; Fig 3K and 3L); MUSM 1439 symphyseal mandible of a juvenile, Locality IQ26 (MZ8; Fig 3N).
in the referred specimens MUSM 1681 (Table 1) and MUSM 2032, and partially distorted in MUSM 900 ( Table 1). The preserved length (from tip of the snout to the posterior angle of left supratemporal fenestra) of the type skull is 623.2 mm (Fig 1A-1C Table 2), including cranial and mandibular remains, represent individuals of equivalent size to the holotype with the exception of MUSM 900, which is a larger skull lacking the anterior half of the snout (Fig 1E). In contrast, MUSM 1439, MUSM  Fig 1D). This specimen only lacks the coronoid bone and posterior end of the retroarticular process. Parts of the posterior lamina of the splenial that are in contact with the surangular are missing or collapsed. No postcranial material has been not definitively identified as pertaining to this new species. Skull. Gryposuchus pachakamue has a long and slender skull. The snout is parallel-sided and tubular in cross section. The anterior maxillary snout has slightly sinuate margins. Proportions of the rostrum correspond to those of Gavialis gangeticus: the rostral length/skull length index is 0.75 (0.76 in G. gangeticus) and the rostral width/postorbital width index is 0.27 (0.25 in G. gangeticus). In front of the orbits, the initial lateral expansion of the skull occurs at the level of the sixteenth maxillary alveolus and continues to widen gradually posteriorly. In dorsal view, the post-rostral skull outline is roughly triangular and the skull table is wide, trapezoidal in shape, and perforated by large supratemporal fenestrae. The supratemporal fenestrae are irregular in outline and widest posteriorly. The orbits are slightly wider than long and markedly smaller than the supratemporal fenestrae. The infratemporal fenestra is roughly ellipsoid in shape. The narial opening is heart-shaped, and its narial rim surface is interrupted anteriorly and anterolaterally by distinctive grooves, but no plateau-like shelf typical of adult males of Gavialis is observed in any specimen. The incisive foramen forms a slender, elongated isosceles triangle. The occipital plate is posteriorly inclined but appears to be at a lesser degree than in other South American gavialoids. Suborbital fenestrae are proportionally longer than those of Gryposuchus colombianus and Gavialis gangeticus, with the anterior end acute and the posterior margin broadly rounded. Pterygoid bullae are located laterodorsal to the posterior half of the palatine bridge.
The premaxillae are expanded at the level of the anterior end of the external naris but not as much as in Gavialis or Gryposuchus colombianus (Fig 4). Long and slender posterior processes of the premaxillae reach the level of the fifth maxillary alveolus (Fig 4A). Thin anterior processes of the nasals are in extensive contact medially with the posterior processes of the premaxillae, a condition also observed in Eogavialis africanus ( Fig 4C). Although the extension of these processes is not precisely symmetrical, their contact length roughly equals the narial opening length. As in all South American gavialoids, each premaxilla bears four alveoli and not five as in Gavialis, Eogavialis, and Eothoracosaurus, most probably by the loss of the second tooth loci of these latter taxa [5] (Fig 4H and 4I). In Gr. pachakamue and other South American gavialoids, where known, the first premaxillary alveoli are close one to another and separated by large gaps from the second alveoli. First and second tooth loci bear alveolar collars projecting ventrally relative to the palatal plane. Second alveoli are the biggest in the premaxilla whereas the fourth premaxillary alveoli are the smallest, as in Gryposuchus species [42] ( Fig  4F). Relatively big foramina are located medial to the third and fourth premaxillary alveoli. The foramina are connected by shallow and bowed grooves. As in other Gryposuchus species, posterior ventral processes of the premaxillae are relatively short, reaching the level of the Evidence for Parallel Evolutionary Ecology in Proto-Amazonian Gavialoid Crocodylians second maxillary alveoli. In Piscogavialis, they surpass the level of the fifth maxillary alveoli ( Fig 4G). Ventrally, the premaxilla bears two anterior medial processes that project into the incisive foramen, and the premaxillary-maxillary suture is jagged and not linear as in all other gavialoids. The maxillae are long and tubular resembling those of Eogavialis, and not dorsoventrally flattened as in Gryposuchus colombianus and Piscogavialis. Dorsally, the maxillae are not in contact due to the forward projecting nasals that contact the slender posterior processes of the premaxillae. As in other Gryposuchus species, modest alveolar salients are observed along most of the maxillary margins. Each maxilla bears 22 tooth positions, all of them subequal in size, exclusive of the last two maxillary alveoli that are the only significantly smaller alveoli ( Fig 2B). Ventrally, the intermaxillary suture extends posteriorly until the level of the fifteenth tooth position. The edge of the maxillary tooth alveoli is higher than the palatal space between both tooth rows. The palatal surface is relatively convex. The anterior interalveolar length is notably larger than the diameter of adjacent alveoli. This length decreases progressively posteriorly, to become equal to alveoli diameter at the fifteenth tooth position and even smaller posterior to it. Dorsally, at the level of the thirteenth tooth position, the maxillae contact the anterior extension of the lacrimals. There is no posterior process of the maxilla between the lacrimal and jugal, in contrast to the condition in Gryposuchus colombianus. The maxillae have a long edentulous posterior process that reaches the level of the postorbital bar (e.g., MUSM 900). Maxillae have no contact with the prefrontals.
The nasals are extremely long, slender bones with close and extensive anterolateral contact with the premaxillae. Additionally, the nasals contact the maxillae and lacrimals laterally, the prefrontals posteriorly, and the frontal posteromedially. The nasals are in contact with each other along the sagittal axis for most of their length. Posteriorly, they are separated by the pointed anterior process of the frontal, almost at the level of the anterior extension of the jugals. The lateral margins of the nasals slightly diverge posteriorly, reaching their greatest transverse diameter just ahead of the anterior processes of the prefrontals. The posterior end of the nasals reach the level of the anterior margins of the orbits, resembling in this aspect Eogavialis and Gryposuchus colombianus, and in distinction from Piscogavialis and Ikanogavialis in which the posterior end of the nasals is located far anterior to this position. The condition in Gavialis is intermediate in this aspect. Whereas the posterior process of the nasal is pointed in all South American gavialoids where this area is known, including Gryposuchus pachakamue, in Gavialis it is strongly denticulate.
The lacrimals are large and roughly triangular in shape. They contact the maxillae and jugals laterally and the nasals and prefrontals medially. Posteriorly, they form the anterior margin of the orbits. The anterior process of the lacrimal exceeds in length that of both the frontal and prefrontal. A small, discrete knob occurs at the orbital margin, lateral to the prefrontal-lacrimal suture.
The jugals form the ventral margins of the orbits and infratemporal fenestrae. In Gryposuchus pachakamue this bone differs from other Gryposuchus species and most other gavialoids, and it more closely resembles the general pattern observed in non-gavialoid crocodylians. Gavialoids, such as Gryposuchus colombianus, Gr. croizati, Ikanogavialis, Gavialis, and to a lesser degree Siquisiquesuchus, present a deep notch immediately anterior to the postorbital pillar. However, in non-gavialoid crocodylians as well as some gavialoids such as Gryposuchus pachakamue, Eothoracosaurus, Eosuchus, and Piscogavialis, the jugal orbital rim progressively descends lateral to the postorbital pillar, thus the ventral margin of the orbit is gently circular. Additionally, in these latter gavialoid taxa and in Siquisiquesuchus, the postorbital pillar reaches the horizontal bar of the jugal medially and a longitudinal sulcus is present between these structures. This latter condition also differs from the distinctive feature observed in gavialoids such as Gavialis and Gryposuchus species other than Gr. pachakamue, in which the postorbital bar lies flush with the lateral surface of the jugal.
The prefrontals are short and rhomboid in shape. They are separated from each other by the frontal and nasals and form part of the anterodorsal orbital margin. From the orbit, the prefrontal-frontal suture follows a gentle semicircular path. In contrast, in Gavialis and Gryposuchus colombianus this suture is sharply angulated, a condition associated with their possession of fully "telescoped orbits" (see below).
The frontal bears a long anterior process that largely exceeds the anterior end of the prefrontal and that of the orbit too. Posteriorly, the main surface of the frontal is slightly concave and weakly sculptured (Fig 2A). Laterally, along the fronto-postorbital suture the skull table is markedly higher than the surrounding areas. Frontal participation in the orbital rim is limited to the posteromedial corner. The frontal in the skull table is anteroposteriorly short due to the large size of the supratemporal fenestrae. Frontal sutures on the skull table are well preserved in MUSM 1681 and MUSM 2032 (Fig 3A and 3B). The fronto-parietal suture is markedly concavoconvex (i.e., M-shaped) between the supratemporal fenestrae and is linear laterally. It runs along the anterior border of, and only briefly enters into, the supratemporal fenestra as in Gryposuchus colombianus. In MUSM 900 the suture runs along the dorsal surface of the skull table, just grazing the margin of the fenestra. The frontal bone and skull table are not elevated relative to the rostrum. The skull table presents an uneven surface, being depressed between the orbits but higher along the postorbital-frontal suture and parietal-supraoccipital region.
The postorbital bones on the skull table are relatively flat and lightly built. They form the anterior corner of the skull table and lack the anterolateral postorbital process characteristic of Gavialis, Gryposuchus colombianus, and Gr. neogaeus. The development of this anterolateral process is variable in Gryposuchus croizati and is small in Siquisiquesuchus and Piscogavialis. The postorbital medially contacts the frontal and slightly contacts the parietal, at the rim of the supratemporal fenestra. The postorbital bone forms part of the dorsal portion of the postorbital bar. As in other gavialoids, the postorbital bar is robust and longer than wide in cross section. On the lateral side of the postorbital bar, lying entirely on the postorbital bone, there is an anteroposteriorly-oriented, crest-like bump, similar to that of Gryposuchus colombianus and Piscogavialis. In MUSM 900, the largest specimen, the bump presents an anterolateral spine similar to that of Gavialis [43]. The suture of the postorbital bone with the jugal is located lower on the postorbital bar, although its precise pattern is not clearly preserved. The lowermost descending process of the postorbital is located posteriorly; it reaches the level of the horizontal bar of the jugal, and contacts the ectopterygoid medially. The postorbital bar is inset from the anterolateral margin of the skull table. Beneath this margin, on the excavated lateral surface, most specimens possess two big foraminae although their size, number, and positions are variable across specimens. It is impossible to determine if the postorbital contacts the quadratojugal medially.
The squamosals are incomplete in the holotype but well preserved in MUSM 1681, MUSM 2032, and particularly in MUSM 900; they occupy the posterolateral corner of the skull table in each. Gryposuchus pachakamue exhibits prong-like posterior extensions of the squamosal, but they are comparatively shorter than those of most South American gavialoids and Argochampsa ( Fig 1E) [6,9,10]. In Gryposuchus colombianus, this region is not well preserved, but it seems that long "prongs" were present as well [8]. On the skull table, the bar that posteriorly limits the supratemporal fenestrae is relatively thin (MUSM 1681), similar to that of Gavialis and Gryposuchus colombianus, but to a lesser degree than in Piscogavialis and Gr. neogaeus [44,45]; in Eogavialis, this bar is comparatively wider. The anterior process of the squamosal lateral to the skull table runs anteroventrally, reaching the posterodorsal end of the postorbital pillar.
The parietal forms the medial and posteromedial margins of the supratemporal fenestrae. Anteriorly, the parietal surface between the supratemporal fenestrae is horizontal. The posterior surface of the parietal gently ascends, thus its main body slopes anteriorly and also laterally, as is typical in gavialoids. Posteriorly, it contacts the supraoccipital along the sagittal axis. The parietal interfenestral bar is comparatively thinner than that in Gavialis and Siquisiquesuchus [10].
The supraoccipital forms an inverse triangle in the occipital plate, excluded from the foramen magnum. It also bears a small, rhomboid-shaped, dorsomedial projection on the cranial roof that points posteriorly. Together with the parietal, this posterior extension becomes a prominent vertical medial crest on the occipital plate (observed in MUSM 1681). The postemporal fenestrae are hardly discernable in any of the specimens, and they seem to be small. Available specimens show asymmetrical development of the postoccipital processes (processus postoccipitales of Kälin) [46] and surrounding area of the supraoccipital, although the so-called "nuchal crest" is not hypertrophied as it is in Gryposuchus colombianus [8].
The quadratojugal forms the posterior margin and corner of the infratemporal fenestra. The anterior process of the quadratojugal comprises a robust spine, as is usually observed in gavialoids. Dorsomedially, the ascending process that bounds most of the posterior margins of the infratemporal fenestra is very thin and long; its dorsal end is not discernible. Sutural contacts of the quadratojugal with the jugal and the quadrate are not parallel, in contrast with the condition in Piscogavialis.
The quadrates are robust and short. The quadrate ventrally bounds the otic aperture and, although incomplete, it seems to occupy part of its posterior margin; thus the quadrate-squamosal suture might have reached the otic aperture along the posterior wall. The small foramen aëreum is located on the mediodorsal surface of the quadrate. The axis of the mandibular condyles is oblique, as it generally is in gavialoids. Medial and lateral condyles are clearly discernable, with the medial one deflected ventrally as in South American gharials, where known, such as in Piscogavialis, Siquisiquesuchus, and Gryposuchus croizati.
The laterosphenoids are partially preserved in the holotype, as well as in MUSM 2032 and MUSM 1681. The anterior dorsal margin of the laterosphenoid, as well as the capitate process, is oriented lateromedially. The dorsal capitate process is massive, and its whole anterior portion is relatively flat.
The palatine bones cover the entire ventral surface of the bridge between the suborbital fenestrae. They are roughly parallel-sided and transversely convex along this bridge. From the anterior limit of the suborbital fenestrae, the lateral margins of the palatine bones converge to a point approximately at the level of the fifteenth maxillary alveolus. The anterior limit of the suborbital fenestra is located at the level of the nineteenth maxillary alveolus (Fig 5A). The palatine-pterygoid suture lies anterior to the rear margin of the suborbital fenestra.
The ectopterygoids are incomplete and badly damaged, but preserved regions are informative. The anterior process of the ectopterygoid is thin and long, resembling that of Piscogavialis, and comprises the posterior half of the margin of the suborbital fenestra. The anterior tip runs medially to the maxillary tooth row until the level of the anterior limit of the penultimate alveolus, as in Piscogavialis. The distance between the tooth row and the anterior process of the ectopterygoid is less than in Gavialis. In ventral view, Gryposuchus pachakamue and Piscogavialis display a substantial contact between the ectopterygoid and jugal, whereas this contact is smaller in the extant gharial Gavialis gangeticus. Available evidence from the same anatomical area is not clear in other South American fossil gharials. Eogavialis specimens do not preserve details of this region, other than in exhibiting a large ectopterygoid-jugal contact. Differing from all other gavialoids, in Argochampsa the ectopterygoid stops far behind the tooth row [47]. The posterior process of the ectopterygoid is short in Gryposuchus pachakamue; thus contact with the pterygoid flanges is reduced relative to the condition in other gavialoids.
The pterygoid widely separates the ectopterygoids from the palatines. In the holotype of Gryposuchus pachakamue, pterygoid bullae in intimate contact with the palatines are observed dorsal to the posterior portion of the palatal bridge. The bullae are neither ovoid in shape nor smooth-surfaced, as in other gavialoids; instead they are kidney-like in shape with irregular prominences and depressions. They are relatively small, flat, and project slightly laterally into the suborbital fenestrae, as in the type specimen of Gryposuchus colombianus. The choana is distorted in MUSM 900 but is not preserved at all in any other specimen.
The exoccipitals are best preserved in MUSM 1681. As in most gavialoids, they are visible in dorsal view, but not as much as in Piscogavialis, Gryposuchus colombianus, and Gr. croizati. The paraoccipital processes are long and encompass a ventral plate-like expansion that completely covers the cranioquadrate foramen. Dorsal to the foramen magnum, the exoccipitals and the supraoccipital are collapsed by deformation, thus positional relationships of these bones are obscured. The ventral processes of the exoccipitals laterally embrace the basioccipital tubera. These processes are robust and anteroposteriorly extended.
The basioccipital is typically gavialoid in being low and laterally expanded ventral to the condyle, producing two pendulous tuberae. The bassioccipital plate is smooth, exclusive of the lateral and ventral margins, in which tuberosities are well developed. Medially, the ventral margin of the plate is excavated. The basisphenoid is anteroposteriorly broad anterior to the basioccipital. The medial eustachian foramen is wider than long.
The braincase region is badly damaged around the prootic in all of the specimens.
Mandible. The mandibular rami (MUSM 987; Fig 1D) are sutured anteriorly, through a long rostral symphysis, to form a Y-shaped structure. Although long, the symphyseal region is proportionally shorter than in Gavialis, Ikanogavialis, Siquisiquesuchus, and Piscogavialis. The mandible is low along the symphysis; and its height progressively increases posteriorly until the level of the glenoid fossa. The external mandibular fenestra is small, eye-shaped and occurs comparatively closer to the articular region than in Gavialis.
doi:10.1371/journal.pone.0152453.g005 external mandibular fenestra (EMF). The first 15 alveoli are implanted within salients along the lateral border of the dentary, providing a sinuous profile to this margin. The dentary (MUSM 987) bears 22 alveoli, all of them sub-equal in size, exclusive of the first and fourth, which are of greater diameter: the second and third are the smallest of the whole series. Gryposuchus colombianus and Gr. croizati present similar alveolar counts (i.e., [22][23] and the same general dental pattern described above. Differences in Gr. pachakamue relative to these other Gryposuchus species include: proportionally higher and more tubular dentary, stronger alveolar salients, and no constriction between the fourth and fifth alveoli (this last character only pertains to comparisons with Gr. colombianus). There are at least four tooth positions located posterior to the level of the mandibular symphysis, as in Piscogavialis, whereas other gavialoids generally possess no more than three. These alveoli are located along the medial limit of the dentary and their lingual wall is formed either by the splenial (i.e., nineteenth and twentieth alveoli) or the surangular (i.e., twenty first and twenty second alveoli).
The splenials wedge out between the dentaries at the level of the twelfth dentary alveolus. As also present in South American gavialoids for which the splenial is known, the anterior process is long, slender, and constricted between the dentaries along the symphysis. Its lateral margin is bowed medially. In MUSM 1428, a right mandible, the medial surface of the splenial, situated within the symphysis, shows no perforation of the foramen intermandibularis oralis. Absent in Tomistoma schlegeli, this foramen occurs within the symphysis of Gavialis gangeticus at the level of the twenty-second alveolus [43]. We identified a distinctive foramen in Gryposuchus pachakamue (MUSM 987 and MUSM1428) and Eogavialis africanus (AMNH 5069)-probably homologous to that of Gavialis-situated just behind the symphysis, within the medially divergent walls of the splenial, at the level of the twentieth alveolus (Fig 5B). Posterior to the symphysis, the splenial occupies the internal half of the mandibular rami. The splenial is excluded from the margins of the foramen mandibularis caudalis, whereas its ventral posterior process extends beyond the rear limit of this foramen (Fig 6B). Other details of the splenial anatomy, particularly the positional relationship to the angular or coronoid, are not preserved in any specimen.
The surangular extends from the medial anterior border of the twentieth alveolus to the lateral surface of the retroarticular process. Within this process, the surangular is partially preserved in MUSM 987, but fully preserved in MUSM1440 (Fig 5D). This latter specimen shows that the surangular fails to reach the posterior tip of the retroarticular process. In lateral view the surangular is low. It is restricted to the posterior angle of the EMF due to the posterior expansion of the dentary and remarkable depth of the angular that embraces this fenestra posteriorly. This condition probably also pertains to Piscogavialis (Fig 6C) and Gryposuchus colombianus [8]. Externally, at the level of the glenoid fossa, there is a large anterodorsally facing foramen (Fig 2D). Among gavialoids, a similar foramen is observed only in Eogavialis africanus (SMNS 11785). Just behind this foramen, the surangular covers the postglenoid process of the articular; therefore this process is hardly seen in lateral view as in Gavialis (Fig 6A and  6C) but in contrast to Piscogavialis in which the surangular ascending lamina is relatively small (Fig 6E). Medially, the surangular almost reaches the lower margin of the EMF as observed in most crocodylians and South American gharials [8,44], except for Piscogavialis (Fig 6E). Posteriorly, the surangular-angular suture reaches the articular at its ventral tip. Proportions of the foramen mandibularis caudalis are similar in Gryposuchus pachakamue and Piscogavialis, with this foramen being relatively longer than in Gavialis.
The angular is excluded from the ventral margin of the EMF being only restricted to its posterior margin, a feature also observed in the newly recovered material of Piscogavialis (MUSM 449; Fig 6E). Other gavialoids such as Gavialis, Gryposuchus neogaeus, and Gr. colombianus [8,44] typically possess a small posterior process of the dentary; therefore the angular borders most of the EMF ventrally. The angular dorsally approaches the posterior process of the dentary, but fails to contact it. An angular-dentary contact in this position is described for Gr. colombianus, but to our knowledge this area is damaged in all referred specimens of that species and this condition is best considered to be unknown in Gr. colombianus. The angular forms the ventral section of the retroarticular process in Gryposuchus pachakamue. The posterior tip of the angular is not preserved.
The articular bone overlies the angular. It is posteriorly and medially inclined relative to the condition in Gavialis (Fig 6). The articular also is inclined in Piscogavialis, a condition probably correlated with the ventromedial projection of the medial hemicondyle of the quadrate. As a consequence, the medial fossa of the retroarticular dorsal surface can be observed in medial view. A prominent longitudinal crest on the dorsal surface of the retroarticular process limits this fossa laterally. This crest is also well developed in Piscogavialis, and other South American gharials, such as Gryposuchus colombianus, and Siquisiquesuchus [10]. Within the glenoid fossa, the articular-surangular suture runs diagonally from its concave anterior margin to the lateral limit of the postglenoid crest. Evidence for Parallel Evolutionary Ecology in Proto-Amazonian Gavialoid Crocodylians Dentition. The adult upper dentition formula is 4 premaxillary + 22 maxillary teeth. This count is similar to other Gryposuchus species [8,9] but much less than that of Piscogavialis [45] and Ikanogavialis [48]. The number of tooth loci in the lower jaw is 22.
Tooth crowns are conical in shape and inclined posterolingually in the plane formed by the fore and aft carinae. Beyond this general morphological similarity, proportions vary substantially relative to tooth locus position in both the upper and lower jaws. Although anterior teeth are long and slender, they do not show the sigmoid-shaped crown of other gavialoids. Posterior teeth are short, robust, and slightly blunt. Although weak longitudinal striae are observed, the surface can be described as virtually smooth.
Juvenile specimens. Iquitos gavialoid specimens of relatively smaller sizes are referred here as juvenile individuals of Gryposuchus pachakamue. MUSM 1988 is an incomplete skull lacking the snout (Fig 3C, 3G and 3H). The width of the skull across the postorbital bars is 90.5 mm whereas in the holotype this diameter is around 178.0 mm. This specimen preserves the skull table as well as the orbital and occipital regions. The medial portion of the posterior region of the skull table has been compressed, thus the foramen magnum is collapsed and the occipital plate is partially distorted. Compared with specimens described above as possessing the adult morphology of Gryposuchus pachakamue, the juvenile MUSM 1988 has a slender postrostral skull, due to the lesser degree of posterior divergence of the jugal lateral margins. The orbits are elongated and equivalent in size to the supratemporal fenestrae. In contrast to the juvenile condition, the adult morphology shows short, wide orbits, and large supratemporal fenestrae, suggesting that a marked allometric development involving the orbital and postorbital regions occurs during ontogeny [16]. Juvenile fenestral shape and proportions resemble those of the holotype of Piscogavialis jugaliperforatus, the latter clearly representing an adult individual. The fronto-parietal suture lies entirely on the skull table without contacting the margin of the supratemporal fenestra. The posterior process of the nasals surpasses the anterior margin of the orbits. The interorbital bridge width is equivalent to the orbit transverse diameter, as in adult individuals. The basioccipital plate and tuberae are well preserved in MUSM 1988 (Fig 3H). The basioccipital plate is wider and dorsoventrally shorter than in the adult condition represented by MUSM 1681 (Fig 1E). The foramen carotideum posterior is well exposed in posterior view, next to the lateral margin of the exoccipital ventral processes. These processes are robust and reach the tuberae. The tuberae are separated medially by a depression posterior to the medial Eustachian foramen. The lateral margins of the basioccipital plate are parallel and composed entirely by the exoccipitals. In adult specimens (Fig 1E), lateral margins become divergent ventrally and a significant portion of the tuberae are extended ventral to the exoccipital, yielding a pendulous shape for the basioccipital-exoccipital structure.
MUSM 1727 comprises a partial snout preserved until the level of the ninth maxillary alveolus (Fig 3K and 3L). As in MUSM 1988, this juvenile snout is comparatively more slender than in specimens representing the adult morphology. Major differences between juveniles and adults lie in the relatively smaller size of the alveoli and, consequently the larger diastemata between adjacent tooth loci, within the juvenile specimens in both upper and lower dental quadrants (Fig 3L, 3N and 3O). Partial mandibles comprising the symphyseal region have been recovered from Localities IQ 136 (MUSM 1682: MZ5; Fig 6O) and IQ 26 (MUSM 1439: MZ8; Fig 6N). We estimate the time span separating these localities in around 3 million years. These represent animals of equivalent size, although smaller alveoli in the specimen recovered from the younger deposits indicate that it probably represents an earlier ontogenetic stage and would have had a larger adult size. Other features of the juvenile specimens, such as the extension of the splenial symphysis and the postsymphyseal tooth loci, seem to be consistent with those of adult individuals.
Gryposuchus cf. pachakamue from Urumaco. Other than differences in size, the anatomical traits of AMU CURS 12, a partial skull from the Late Miocene Urumaco Formation of Venezuela [34] are essentially identical to those of the new species Gryposuchus pachakamue ( Fig  3D). As in the Peruvian specimens, AMU CURS 12 bears a trapezoidal skull table that is extensively perforated by supratemporal fenestrae, and circular and only moderately "telescoped" orbits. Additionally, the Venezuelan specimen resembles Gr. pachakamue and also differs from other Gryposuchus species in having its interorbital bridge width equivalent to the transverse diameter of the orbit, and a postorbital pillar contacting the horizontal bar of the jugal medially. Bone sutures are not discernable on the Venezuelan specimen.

Anatomy of the orbital region in gavialoids
Circumorbital bones include the frontal, prefrontals, lacrimals, jugals, and postorbitals. The configuration of these bones determines the orbital shape and the general position of the eyeballs in crocodylians. In the adult individuals of Gryposuchus pachakamue, the orbits are slightly wider than long, closer in shape to those of Eogavialis africanus than to the circular orbits of Gavialis and other Gryposuchus species. In contrast to the adult condition, the orbital shape and proportions of a juvenile specimen (MUSM 1988; Fig 3C) resemble those of most crocodylians, including gavialoids such as Piscogavialis and Eothoracosaurus, in which the orbits are comparatively longer anteroposteriorly. The relative diameter of the interorbital bridge (i.e., frontal bone width/orbital width ratio) defines the separation between the orbits. Brevirostrine crocodylians, Eogavialis, Piscogavialis, and Ikanogavialis generally bear a narrow interorbital bridge diameter (i.e., ratio <1.0) compared to Gryposuchus pachakamue, in which the width of the frontal bone is essentially equivalent to the width of the orbit and this ratio is consistently higher (more slender interorbital bridge) than that of adult Gavialis and other species of Gryposuchus (i.e., ratio >1; Fig 3 and Fig 7C). The dorsal orbital margins are upturned in many alligatoroids and crocodyloids [13,16]. Among gavialoids, Eogavialis, Gryposuchus pachakamue, Piscogavialis, and Siquisiquesuchus have also upturned dorsal orbital margins, but to a lesser degree than in Ikanogavialis, Gavialis, and other Gryposuchus species. In these latter taxa, dorsal and posterior margins are sharp and projected towards the orbits. Similarly, the anterior orbital margin (i.e., lacrimal bones) of Eogavialis, Gr. pachakamue, and Piscogavialis lies flush with the rostral surface contrasting with the raised anterior margins of other gavialoids. The jugals of Gr. pachakamue differ from those of other Gryposuchus species and most other gavialoids, more closely resembling the general pattern of non-gavialoid crocodylians. Gavialoids such as Gryposuchus colombianus, Gr. croizati, Ikanogavialis, Gavialis, and to a lesser degree Siquisiquesuchus, present a deep notch immediately anterior to the postorbital pillar. However, in non-gavialoid crocodylians and some gavialoids such as Gryposuchus pachakamue, Eothoracosaurus, Eosuchus, and Piscogavialis, the jugal orbital rim progressively descends lateral to the postorbital pillar, thus in those taxa the ventral margin of the orbit is gently circular. Additionally, in those taxa as well as in Siquisiquesuchus, the postorbital pillar reaches the horizontal bar of the jugal medially and a longitudinal sulcus is present between these structures. This latter condition also differs from a distinctive feature observed in many other gavialoids, such as Gavialis and other Gryposuchus species, in which the postorbital bar lies flush with the lateral surface of the jugal (Fig 7C). Postorbital bones on the skull table are relatively flat in Gr. pachakamue. They form the anterior corner of the skull table and lack the anterolateral postorbital process characteristic of Gavialis, Gryposuchus colombianus, and Gr. neogaeus.
In summary, in contrast to the relatively conservative morphology of the orbital region of most brevirostrine crocodylians, all gavialoids present conspicuous modifications within the circumorbital anatomy associated with possessing protruding eyes or "telescoped" orbits. The degree of development of the "telescoped" orbit condition varies among gavialoid taxa and usually includes at least some of the following character states: upturned anterior, dorsal, and posterior orbital margins; anterolateral postorbital process; ventral orbital margin with a prominent notch; postorbital pillar laterally displaced; and orbits widely separated. Adult Gavialis, Gryposuchus colombianus, and Gr. croizati possess all of these traits and therefore the fully "telescoped" orbit condition. Juvenile individuals of Gavialis lack widely separated orbits; other characters associated with the "telescoped" orbits condition are recognized but they are less marked than in adult skulls. The general configuration of this area in the new gavialoid Gryposuchus pachakamue depicts weak development of "telescoped" orbits. This is interpreted as an incipient condition for more extensive telescoping of the orbits in later diverging relatives, based on ancestral state transformations determined from the maximum parsimony phylogeny.

Phylogenetic analysis
Our first parsimony analysis retained 45 equally optimal trees with a minimum length of 538 steps. The strict consensus phylogeny (Fig 7A) calculated from those trees yielded the following statistics: length = 553; consistency index (CI) = 0.461; retention index (RI) = 0.743. Our strict consensus tree shows general coincidence with previous morphological and molecular analyses for major relationships within crocodylian clades [5,6,13], but as in other morphological analyses, it differs markedly from results based only on molecular data in the hypothesized affinities between Gavialis and Tomistoma (see Results in [3,4]). Thus, we found strong Bremer support for the monophyly of Gavialoidea, with Gavialis having a much closer relationship to Cretaceous gavialoid taxa such as Eothoracosaurus than to extant Tomistoma, the latter being most closely allied with Crocodylus within the Crocodyloidea. Within the Alligatoroidea, Culebrasuchus mesoamericanus is closer to Alligator mississippiensis than to Caiman crocodilus [11].
The new Pebasian species, Gryposuchus pachakamue, is recovered within gavialoids as the sister taxon of the Gr. colombianus + Gr. croizati clade. Thus, our results suggest that all known Amazonian gavialoids belong to a single monophyletic taxon (the inclusive species of Gryposuchus) characterized by association with riverine and lacustrine-tidal paleoenvironments. On the other hand, relationships among South American taxa usually associated with coastal marine paleoenvironments show no resolution: Siquisiquesuchus, Piscogavialis, and Ikanogavialis lie in a polytomy together with Amazonian Gryposuchus and Indo-Asian Gavialis. As a consequence, this first analysis finds no support for the monophyly of a clade comprising South American gavialoids, as had been suggested by previous studies [10,39,49]. This clade, namely Gryposuchinae as proposed by Vélez-Juarbe et al. [39], was usually linked only with weak support or collapsed when the analysis included Argochampsa [9]. It is noteworthy that we found support, although low, for a novel association between African Argochampsa and the Caribbean taxon Aktiogavialis. This clade is supported by the presence of long supratemporal fenestrae (character 191-0), a character state that implies a reversal within these gavialoids. This might be considered problematic since Argochampsa, as an early representative of the clade, could have retained the ancestral condition for Crocodylia if character polarization based on the phylogenetic position of "thoracosaurs" is erroneous. Nevertheless, further comparisons between Argochampsa and Aktiogavialis allows us to recognize other similarities not yet included in phylogenetic analyses, but potentially supporting their close relationship, such as the relative proportions of the skull table in dorsal view and the presence of a distinctive shallow fossa in the anterior margin of the supratemporal fenestra. This fossa has not been identified in any other crocodylian species at any ontogenetic stage [39], but it is clearly visible in the holotype of both Argochampsa and Aktiogavialis, although barely discernable in other individuals of the African taxon Argochampsa (Stéphane Jouve, pers. comm., 2015). The Argochampsa + Aktiogavialis group is found as most closely related to the clade of South American gavialoids plus Indo-Asian Gavialis, and with Eogavialis as the nearest outgroup to all of these taxa. Therefore, African Argochampsa is more closely related to Indo-Asian Gavialis as was previously suggested by others [6,9], than it is to African Eogavialis. In fact, both Paleogene African taxa, Argochampsa and Eogavialis, share key characters with more anatomicallyderived gavialoids like Gavialis. The clade encompassing this subset of late-diverging gavialoids is here termed "gharials" (Fig 7) and it is characterized by a posteriorly pointing supraoccipital (character 160-1), anteroposteriorly wide basisphenoid (character 172-1), robust exoccipital ventral process to basioccipital tubera (character 176-1), and basioccipital plate with ventrally divergent sides (character 196-1). Other conspicuous features of the gavialoid temporal region present in Gryposuchus pachakamue, such as ovoid infratemporal fenestra (character 204-1) and anteriorly flaring squamosal groove (character 147-1), also are observed in the tomistomine Thecachampsa americana and might represent independent, convergent acquisitions. As in most crocodylians, the infratemporal fenestra is triangular in Gavialis (character 147-0), and is most parsimoniously regarded as a reversal in this latter taxon. Paleocene Eosuchus and Cretaceous "thoracosaurs" from the northern hemisphere are identified in this analysis as the first and second basalmost diverging branches among the gavialoids, respectively.
To better understand morphological transformations, we performed a second analysis, this time excluding individual characters revealed in the first analysis as highly homoplastic, since results suggested that they were acquired independently up to three times among later-diverging gavialoids, including Gryposuchus and Gavialis. In addition, this analysis collapsed states of character 138 (as they may not be independent) and deleted state 2 from character 137, both of which describe morphology of the orbital margin associated with "telescoped" orbits, which the all-data parsimony analysis indicate clearly evolved multiple times independently within gavialoids (see Discussion). Regarding character 138, two states were combined in this second analysis as all taxa with state 2 also preserve state 1 (non-independent states), and all taxa previously coded as having state 2 (i.e., 138-2: dorsal and posterior orbital edges upturned) were recoded as 1 (i.e., 138-1: dorsal edges of orbits upturned). This search retained 24 equally optimal trees, with a length of 533 steps. The corresponding strict consensus tree is 542 steps long (CI = 0.467; RI = 0.750; Fig 7B and S2 Fig).
Compared to the first analysis, this approach substantially increased resolution of South American gavialoid relationships. The monophyly of gryposuchines, exclusive of the Caribbean taxon Aktiogavialis, is recovered and Gavialis is identified as the sister clade of gryposuchines. Character support for the gryposuchines includes having four premaxillary alveoli (character 87-1), lack of exposure of the prootic on external braincase wall (character 164-1), a quadrate with ventromedially projected medial hemicondyle (character 181-4), and a retroarticular longitudinal crest (character 203-1). This last character is present in all South American gharials preserving this region [9,10], including an unnamed Late Oligocene gharial from Pirabas, Brazil [50]. The distinctive long posterior or posterolateral projections of the squamosal, (i.e., squamosal prongs) no longer diagnose Gryposuchinae as had been hypothesized previously (see Results in [10]), because Argochampsa also clearly possesses this feature [6]. Within gryposuchines, there is low support for a clade comprising Ikanogavialis, Siquisiquesuchus, and Piscogavialis. Other gavialoid relationships were unaffected relative to the initial analysis. Discussions below are based on this second phylogenetic analysis, unless otherwise noted.

Discussion
The evolutionary ecology of gavialoids: Evidence from Amazonia Gryposuchus pachakamue inhabited the heart of the proto-Amazonian mega-wetlands ecosystem during the Middle Miocene, from around 16 to 13 Ma, and it represents the oldest known record of a gavialoid from this area. Remains belonging to this new taxon were consistently recovered from deposits depicting shallow lacustrine paleoenvironments. As the basalmost species of the Gryposuchus clade, it provides essential evidence for accurately reconstructing the ancestral anatomy and ecology of this clade and provides evidence of parallel evolution of distinctive rostral and orbital anatomy in gavialoids. Our phylogenetic analyses reconstruct the acquisition of widely separated orbits as independent evolutionary events in Asian Gavialis and later-diverging Gryposuchus species (Gr. colombianus + Gr. croizati) in South America. As a consequence, the comparatively slender interorbital bridge of Gr. pachakamue is primitive for all gavialoids (Fig 7C). A wide interorbital bridge is associated with possessing "telescoped" orbits. Traits associated with fully "telescoped" orbits, as is observed in Gavialis and advanced Gryposuchus species (but absent in Gr. pachakamue) include: postorbital bar flush with lateral jugal surface (character 135-1), upturned dorsal and posterior orbital margins (character 137-2), and ventral orbital margin with a prominent notch (character 138-1). Parsimony analyses also suggest parallel development for two of these character states, indicating that the "telescoped" orbit condition is homoplastic in gavialoids and occurred independently in advanced South American Gryposuchus and Asian Gavialis species.
By mapping the phylogenetic tree onto the first two axes of the PCA axes, we examined ecological associations of circumorbital bone arrangements throughout gavialoid evolution (Fig 8). Higher PC1 and PC2 scores define a morphospace of gavialoids found in coastal marine deposits, such as Piscogavialis and Argochampsa, whereas lower scores along these PC axes correspond to the morphospace including taxa with fully "telescoped" orbits. Lacustrine Gryposuchus pachakamue and Eogavialis africanus, the latter lacking definitive data on its paleoenvironmental provenance, exemplify the intermediate morphospace. The fully "telescoped" orbit morphospace is represented by Gavialis and Gryposuchus colombianus; since they are distantly related, this morphospace depicts the parallel evolution of this distinctive cranial anatomy, apparently associated with convergent specialization in a freshwater habitat and similar visually enhanced feeding strategies. Distinct features of this morphotype include: upturned anterior, dorsal, and posterior orbital margins; ventral orbital margin with a prominent notch; postorbital pillar laterally displaced; and orbits widely separated. Gavialoids display different degrees of the "telescoped" condition. extreme of PC1 present slender skull tables and interorbital bridge, long orbits and prefrontals, and laterally oriented anterior processes of the jugals, whereas those on the negative extreme bear broad skull tables, wide posterior portion of the interorbital bridge and orbits, short orbits, short prefrontals, and medially oriented anterior processes of the jugals. PC2 correlates with the relative length of the pre-orbital bones, involving mostly the frontal and lacrimals and the width of the prefrontals. Species with higher PC2 scores have comparatively longer and more slender frontals and a narrow interorbital bridge. Taxa with lower scores present short frontal and lacrimal bones and short and wide anterior portion of the interorbital bridge. The phylogenetic morphospace of the orbital and circumorbital region in Miocene South American gavialoids covers most of the variation of the whole clade. Taxon abbreviations: Bor., Borealosuchus; Pal., Paleosuchus; The., Thecachampsa; Tho., Thoracosaurus (S1 Appendix).
The feeding strategy of the extant gharial Gavialis gangeticus involve active use of the telescoped eyes and integumentary sense organs in capturing fishes in streams [51], as its habitat is confined to freshwater settings of the Indian subcontinent, notably restricted to riverine environments [52]. Fossil gharials with well-developed "telescoped" orbits are usually found in depositional settings documenting fluvial-dominated paleoenvironments. This similar habitat association in the extant and extinct convergent taxa offers further support for the adaptive value of possessing "telescoped" orbits and suggests that this morphotype in fossil gavialoids typically is correlated with taxa living in riverine ecosystems. Plio-Pleistocene fossil Gavialis species inhabited, and might have dispersed geographically via, fluvial systems then occurring between Indo-Pakistan and Southeast Asia [53]. In South America, late Middle Miocene Gryposuchus colombianus is restricted to fluvial-influenced settings of the Pebas Mega-Wetland System at La Venta (Colombia) [8,54] and Fitzcarrald (Peru) [11,55,56], close to the rapidly rising Andes, whereas this species is absent from coeval deposits representing brackish, dysoxic lakes, swamps and deltas of the same biome at Iquitos (Peru) where Gryposuchus pachakamue consistently occurs. This latter species not only is the sole gavialoid in the Iquitos fauna, but the sole longirostrine crocodylian within a highly diversified community dominated by small caimans with crushing dentitions reflecting a specialized malacophagous diet [11]. In younger beds at Iquitos (Molluscan Zone MZ9 or younger), we only recovered a "telescoped" orbit gavialoid from one locality (IQ125) in the "Uppermost Pebas Formation" [33], at a time corresponding to development of new regional fluvial-dominated conditions, attributable to a peak in Andean uplift and the onset of the Amazon River System [11].
Within a wider environmental context, the demise of the Pebas Mega-Wetland System and subsequent establishment of the fluvio-tidal Acre Phase in the Amazonian Basin at around 10.5 Ma seems to have promoted diversification, size increase, and specialization in gharials (Fig 9) compared with the presence of only one species per Middle Miocene-aged crocodylian fauna (i.e., Iquitos, La Venta, and Fitzcarrald). In contrast, Late Miocene assemblages of Acre (Brazil) contain four gavialoid taxa [57], including the conspicuous record of forms with "telescoped" orbits (RS-G, pers. obs). Although the same number of taxa is observed in the Late Miocene Urumaco Formation of Venezuela, matches between the dominant environment and specific morphotypes is far more complicated to assess there. The Urumaco Formation preserves several aquatic environments that were potential habitats for crocodylians within the Paleo-Orinoco Basin, such as delta plains, swamps, marginal marine embayments, and rivers, all with a marine influence from the Caribbean [35]. Consistent with this patchwork of environments, the Urumaco gavialoid species exhibit a highly disparate array of cranial circumorbital configurations [9,34,48]. In Urumaco a telescoped orbit specialization was fully attained by the giant Gryposuchus croizati, most likely a dweller in fluvial settings. Based on its paleoenvironmental preference in the Iquitos fauna, Gryposuchus cf. pachakamue from Urumaco instead might have occupied relicts of the brackish lacustrine ecosystems once widely distributed in proto-Amazonia [32]. Therefore, the record of this gavialoid at Urumaco adds evidence for the persistence of Pebasian aquatic conditions during the Late Miocene in northernmost South America, where the lower course of the Pebasian proto-Amazonian System drainage was formerly situated [11,58].

Caribbean and proto-Amazonian biogeographic relationships and the origin of Neotropical gharials
Inasmuch as gharials were inhabitants in the Caribbean region at least since the Late Oligocene [39], the Caribbean Portal at the mouth of the proto-Amazonian aquatic system might have played an important role for the invasion of gavialoids into the inland areas. Siquisiquesuchus was an early Miocene Caribbean inhabitant of present-day Venezuelan coasts [10]. Other marine and fresh water units in the region document numerous remains of gharials throughout the Miocene [9,34,62,63], showing persistent faunas living close to the Caribbean Portal. Either by means of marine transgressions or via riverine drainage systems, a continuous aquatic corridor united the Caribbean Sea with western proto-Amazonia for most of the Paleogene and early Neogene [19,32,61]. The heart of proto-Amazonia during the extensive and persistent Pebas Mega-Wetland System was continuously connected with the Caribbean Portal by a northward-flowing trunk drainage [32,61]. Although not coeval, evidence of the prevailing aquatic connections that linked both areas through the Llanos (Colombia) until the Late-Middle Miocene boundary is founded on the record of Gryposuchus pachakamue from both Urumaco and Iquitos, as well as other biotic indicators [11,[64][65][66][67][68]. Specific environmental conditions fostering high diversity of Late Miocene Urumaco gavialoid diversity remain elusive.
Our phylogenetic analyses propose sister-group relationships of Caribbean Aktiogavialis and African Argochampsa (Fig 7 and Fig 9). In the context of our time-calibrated phylogenetic tree, this association would suggest an African origin for the Caribbean gavialoid, probably by western transatlantic dispersal, as suggested by early authors [42,69] and more recently by Vélez-Juarbe et al. [39]. Argochampsa, Aktiogavialis, and the oldest records of South American Evidence for Parallel Evolutionary Ecology in Proto-Amazonian Gavialoid Crocodylians gharials all were found in deposits from coastal marine settings [10,39,47,50,70]. Whether fossil gharials were strictly marine or not, distinct lines of evidence indicate that gavialoids flourished during high sea surface levels and temperatures (SST) of the Paleocene and Eocene epochs [59]. During this time interval (~60-45 Ma), paleogeographic reconstructions depict a tropical marine connection from India and Africa to the Caribbean and northern South America through the Tethys Ocean and Caribbean Sea [60]. Although the fossil record is far from complete in these areas, our time-calibrated phylogenetic tree and the occurrence of advanced gavialoids along both the Tethys and Caribbean coasts and islands, suggest that this marine realm could have served as a preferred habitat and dispersal system for gavialoids during the Paleogene (Fig 9). The presence of gharials in South America most likely resulted from a single transoceanic colonization event and subsequent diversification in South America, distinct from the one that gave rise to the Caribbean taxon Aktiogavialis.

Conclusions
Gavialoid history exhibits independent acquisitions of the "telescoped" orbits condition. Analyses of the new Pebasian species Gryposuchus pachakamue and other South American fossil gavialoids document high plasticity in orbital anatomy, which appears to have been strongly correlated with a visually enhanced feeding strategy and environmental circumstances. Morphospaces occupied by fluvial and coastal marine specialists are identified by quantitative analysis of orbital and circumorbital shape variation. In light of the phylogenetic history, a fluvial habitus in South American gharials is derived from ancestral lacustrine-deltaic forms with incipient development of protruding eyes or telescoped orbits. The circumorbital region of coastal marine gavialoids is closer in morphology to that of brevirostrine crocodylians. Identifying morphological steps of parallel evolution and ancestral ecological habitus in gavialoids provides models for reconstructing puzzling phylogenetic histories and adaptive radiations within extinct crocodylomorphs clades with elongated rostrums, such as thalattosuchians, dyrosaurids, and pholidosaurids [71]. Proto-Amazonian connections with the Caribbean Sea to the north, and the subsequent onset of the transcontinental Amazon River System draining eastward, provided multiple habitats and conditions for gavialoid colonizations of new areas and extensive morphological diversification in South America throughout the mid-late Cenozoic.
Supporting Information S1 Appendix. Phylogenetic and morphometric data. (AMNH); E. "Dino" Frey (SMNK), and R. Schoch (SMNS) for access to comparative collections; A. Valdés-Velasquez, M. Orliac, G. Billet, J. Martin, R. Lebrun, and F. Pujos for invaluable discussions; V. DeLeon for providing substantial insights on morphometric analyses and interpretation, and M. Burga for assistance in gathering data for this analysis; and P. Bona, C. Brochu, S. Jouve, T. Scheyer, and J. Moreno-Bernal for providing visual material of crocodylians. We are much indebted to C. Brochu and J. Martin for their constructive critical reviews. This is ISEM publication 2016-023 S.