Integrative overview of the herpetofauna from Serra da Mocidade, a granitic mountain range in northern Brazil

Abstract The Brazilian mountain ranges from the Guiana Shield highlands are largely unexplored, with an understudied herpetofauna. Here the amphibian and reptile species diversity of the remote Serra da Mocidade mountain range, located in extreme northern Brazil, is reported upon, and biogeographical affinities and taxonomic highlights are discussed. A 22-days expedition to this mountain range was undertaken during which specimens were sampled at four distinct altitudinal levels (600, 960, 1,060 and 1,365 m above sea level) using six complementary methods. Specimens were identified through an integrated approach that considered morphological, bioacoustical, and molecular analyses. Fifty-one species (23 amphibians and 28 reptiles) were found, a comparable richness to other mountain ranges in the region. The recorded assemblage showed a mixed compositional influence from assemblages typical of other mountain ranges and lowland forest habitats in the region. Most of the taxa occupying the Serra da Mocidade mountain range are typical of the Guiana Shield or widely distributed in the Amazon. Extensions of known distribution ranges and candidate undescribed taxa are also recorded. This is the first herpetofaunal expedition that accessed the higher altitudinal levels of this mountain range, contributing to the basic knowledge of these groups in remote areas.


Introduction
Scientific interest in mountain ranges arises primarily because they are characterized by a geographical isolation associated with differential availability of topographical, climatic, and edaphic conditions along the altitudinal gradients (Haslett 1997, Martinelli 2007, Körner et al. 2017. These characteristics provide ideal conditions for the development of unique evolutionary lineages and occurrence of a variety of endemic taxa (Lomolino 2001, Rull 2005, Hoorn et al. 2013, Nogué et al. 2013). This speciation hotspot pattern occurs in the mountain ranges of northern South America (Guiana Shield highlands), which lies on the ancient terrain of the Cratonic Guiana Shield (Hershkovitz 1969, Désamoré et al. 2014, Salerno et al. 2012, Bonaccorso and Guayasamin 2013. Despite a history of different concepts and geographic boundaries (Mayr and Phelps 1967, Huber 1988a, Kok 2013a, these high altitude areas are considered as a distinct biogeographic region (Morrone 2014), presenting biotic affinities with the megadiverse Amazon and Andean regions (Duellman 1979, Salerno et al. 2012, Mannion et al. 2014.
Initial discoveries concerning the amphibian and reptile diversity from Guiana Shield highlands were reported by localized expeditions at the transition between the nineteenth and twentieth centuries (Boulenger 1900, Roze 1958a. Knowledge increased exponentially when helicopters facilitated access to remote mountains (Aubrecht et al. 2012), leading to several expeditions focused on describing the assemblages of particular localities. Such studies brought to prominence several endemic taxa from the highlands (e.g., Gorzula 1992, Myers 1997, Myers and Donnelly 1992, MacCulloch et al. 2007, Barrio-Amorós and Brewer-Carias 2008, Kok 2008, 2009a, 2009b, 2010, 2013b, Kok and Rivas 2011, and some resulted in broad outlines of the main biogeographical patterns of these taxa (e.g., Hoogmoed 1979, Duellman 1999, Gorzula and Señaris 1999, McDiarmid and Donnelly 2005. More recently, studies have shown that diversification and evolutionary patterns of distinct species were associated with the landscape history of the region , 2017, Salerno et al. 2012, Vacher et al. 2017. While several intriguing patterns have been found for some highland lineages, such as recent diversification and low genetic divergence among mountains (Salerno et al. 2012, 2017, overall knowledge on the geographical and altitudinal distribution patterns of amphibians and reptiles from Guiana Shield highlands and their drivers are far from being fully understood. One of the main reasons is the occurrence of huge sampling gaps, mostly due to the short-term nature of inventories and to the difficult and costly access to highland areas, some of which remain unexplored (Aubrecht et al. 2012). The sampling deficiency in the Brazilian region of the Guiana Shield highlands highlights the impor-tance of exploring these areas for biodiversity and biogeography assessments. Sampling in novel mountain ranges will most likely result in the discovery of new taxa and unique lineages, and geographical range extensions, all of which can contribute to the conservation of these threatened regions (Rull et al. 2016).
Furthermore, most specimens currently collected during biological inventories are identified using morphological characteristics. However, given the pervasive occurrence of cryptic diversity (Vences et al. 2005, Vences andWake 2007), particularly in Neotropical amphibians and reptiles (Fouquet et al. 2007a, c, Geurgas and Rodrigues 2010, Oliveira et al. 2016, Kok et al. 2016, other methods are being used to reveal the hidden diversity in remote areas with difficult access. For example, molecular techniques may contribute to indicate the presence of undescribed species and detect cryptic speciation through divergence in DNA sequences, and this short-term result might take longer to be achieved using a single taxonomic data source (Vences et al. 2005, Vences and Wake 2007, Fouquet et al. 2007a, Paz and Crawford 2012. Recently, a multidisciplinary initiative conducted an expedition ("Biodiversity of the Serra da Mocidade") to inventory the biological diversity of distinct taxonomic groups in the poorly known region of the Serra da Mocidade, a remote granitic mountain range located in northern Brazil (INPA, 2016), highly isolated from other mountains and with difficult access. Here we present and discuss the diversity of amphibians and reptiles of the area and their biotic affinities, using an integrative approach combining morphological, bioacoustical and molecular analyses to identify specimens. We found remarkable records, and make observations concerning species taxonomy, ecology and distribution patterns in the Guiana Shield highlands region.

Study area
The Guiana Shield highlands region is located in northern South America, within the limits of Venezuela, Guyana, Suriname, and Brazil ( Fig. 1), and is composed of mountain ranges covered by dense forests, shrubby vegetation or moss forests, surrounded by a lowland matrix of either tropical forests or savanna ecosystems (Mayr and Phelps 1967, Huber 1988b, McDiarmid and Donnelly 2005. These ancient mountains date to the Precambrian period (1.8-2.5 billion years) (Santos et al. 2003, Nogué et al. 2009, Kok 2013a, and have two main geological origins: sedimentary rocks, which are currently exposed as abrupt vertical elevations with tabular tops (called "tepuis") and igneous-metamorphic rocks, which lie beneath the sedimentary rocks and may also be exposed as granitic uplifts with a mountain-like appearance (Steyermark 1986, Schubert and Briceño 1987, Huber 1995, Hoorn and Wesselingh 2010. This study was conducted at Serra da Mocidade (Figs 1, 2), a complex of granitic mountains located in extreme northern Brazil, within the limits of Caracaraí municipality, Roraima state, with a mean altitudinal level of 1,000 m above sea level (hereafter, asl) and peaks over 1,900 m asl. The Serra da Mocidade area is protected by a federal conservation unit (Parque Nacional da Serra da Mocidade), an indigenous Yanomami territory and a military area, property of the Brazilian Army (Ferreira et al. 2014, Ministério do Meio Ambiente 2016. This mountain range is isolated in relation to other Guiana Shield highlands mountains, with the closest mountain range (Serra do Aracá) ca. 100 km distant to the west. The lowland regional climate has low annual temperature variation, ranging from 24° C to 27° C (average 26° C), a rainy season from April to September and a dry season from October to March (Alvares et al. 2013). Annual rainfall exceeds 2,500 mm (Sombroek 2001, Alvares et al. 2013. Daily thermal oscillation increases at higher altitudes, where temperatures are lower and the influence of wind is stronger, preventing establishment of taller forests (Haslett 1997). The base and lower slopes of the mountain range are covered by dense tropical submontane and montane rainforests, while increased humidity at higher altitudes fosters the development of extensive epiphyte and moss coverage on tree trunks (Ministério do Meio Ambiente 2016).

Sampling areas and species survey
Logistic support from the Brazilian Army allowed aerial access at altitudes only accessible by helicopters, and the installation of two base camps from which it was possible to reach different altitudinal levels ( Fig. 1). At base camp #1 at 600 m asl (01°36'N, 61°54'W), we sampled areas in the foothills of the mountains, while base camp #2 at 1,060 m (01°42'N, 61°47'W) allowed access to altitudinal levels of 960 m asl (01°42'N, 61°48'W) and 1,365 m asl (01°43'N, 61°45'W). Altitudinal levels closer to the camp (900-1,100 m asl) were sampled more intensively due to difficult access in the areas located at 1,365 m asl, which we sampled during a single operation. Although herpetofaunal rapid inventories already have been conducted in the lowlands of this region (Ministério do Meio Ambiente 2016), this is the first expedition that reached the higher altitudes of this mountain range.
Surveys were conducted in two teams, each with three trained herpetologists, one from each base camp (first team at base camp #1 for 15 to 23 January 2016, and second team at base camp #2 for 25 January to 06 February 2016), totaling 22 days of field work and 66 man-days of sampling effort. To maximize characterization of the regional herpetofaunal diversity we used six complementary survey methods to detect individuals: (1) active surveys  conducted during day and night, where individuals were visually detected or heard in as many microhabitat as possible and manually captured; (2) pitfall traps Christman 1982, Heyer et al. 1994), installed in two sampling lines, each of ten 60-liters buckets spaced every 10 meters (total sampling line length 100 m), and buried in the ground with the opening at the surface level, interleaved with a pole-supported plastic fence, the lower part of which was buried in the ground; (3) trammel nets, which were installed in streams near the base camps, in order to catch turtles, aquatic amphibians and squamates (Campbell and Christman 1982); (4) hook with meat bait, specifically for turtles; (5) glue traps, installed in fallen logs and tree trunks aiming capture of arboreal squamates and (6) shotgun. All traps were visited every 24 hours or less throughout the sampling period at each base camp. Specimens obtained using methods other than the ones cited above were considered as occasional encounters (Martins andOliveira 1998, Heyer et al. 1994).

Collection and species identification
Specimens were killed with an injection of Thiopental® or Lidocaine, fixed with 10% formalin, and preserved in 70% ethanol. We removed muscle or liver tissue from specimens before the fixation process and stored it in absolute ethanol. Voucher specimens and tissue samples were deposited in the Collections of Amphibians and Reptiles (IN-PA-H) and of Genetics Resources (INPA-HT) of the Instituto Nacional de Pesquisas da Amazônia (INPA), Manaus, Amazonas, Brazil, respectively.

Bioacoustics
The calls from some amphibian species were recorded in uncompressed wav format, with a Zoom H1 Handy Recorder (Zoom Corporations, Tokyo, Japan) equipped with an internal microphone, distant about 1-2 m to the emitter. The calls were digitised and analysed using Raven Pro. 1.5 (Cornell Laboratory of Ornithology) at a sampling frequency of 44 KHz and 16-bit resolution. Call structures were visually analysed in the spectrograms, and we measured the following quantitative parameters, considered in amphibian taxonomy (Köhler et al. 2017): call duration (s), inter-call interval (s), pulse duration (s), pulse rate (pulses/s) and dominant frequency (kHz), summarized in mean values ± standard deviation. The data were compared to calls described in the literature for each analysed species (Cardoso and Haddad 1984, Hoogmoed 1990, Donnelly and Myers 1991, Fouquet et al. 2007b, Morais et al. 2012, Fouquet et al. 2015.

Molecular analyses for uncertain species
When specimens were part of groups already considered as a complex of multiple species or when we considered the possibility of an undescribed taxon, we explored their taxonomic status using DNA sequences of mitochondrial gene 16S, a standard marker for amphibians and reptiles (Vences et al. 2012). Genomic DNA was isolated from collected tissues using a commercial kit (Wizard, Promega Corp., Madison, WI), and target region amplified via the Polymerase Chain Reaction (PCR) using primers 16Sar and 16Sbr (Palumbi et al. 1991). Purified PCR products were sequenced using the Big Dye Terminator sequencing kit (Applied Biosystems, Waltham, USA) in automated sequencer ABI 3130 XL (Applied Biosystems, Waltham, USA) at Thematic Laboratory of Molecular Biology at INPA.
DNA sequences obtained were compared with those available for closely related taxa in GenBank . Sequences were manually edited and aligned with CLUSTAL X algorithm run on MEGA 6.06 software (Tamura et al. 2013). Using the same software, we generated maximum likelihood (ML) phylogenetic trees for each taxonomic group inferred by 5,000 bootstrap replicates, and estimated genetic distances between main lineages of each taxon using uncorrected pairwise distances. Although the threshold of the genetic distance percentage that represents an interspecific variation is arbitrary and varies according to the group diversification, we followed Fouquet et al. (2007a, c) and considered the possibility of new taxa when genetic distances were above 3%. In these cases, the final definition of the taxonomic status of a given specimen was thus the result of an integrated interpretation of morphological, bioacoustical, and molecular results.
To identify biogeographical patterns and evaluate the contribution of adjacent regions to the composition of the local herpetofauna at Serra da Mocidade, we classified species according to their known geographical distribution, delimiting boundaries of biogeographic regions based on the main geological compartments for Northern South America (Gibbs and Barron 1993, Aleixo and Rossetti 2007, Hoorn and Wesselingh 2010, already known to influence the biotic distribution (Ávila-Pires et al. 1995, Aleixo andRossetti 2007). Therefore, the species were classified in the following categories, which decrease in geographical scale: Widely distributed in Amazonia (WD), for species with wide geographical ranges throughout this region; Western Amazonia (WA) and Eastern Amazonia (EA) for species typical of these macro-regions, influenced by sedimentary basin of the Amazonas River and crystalline shields, respectively; Andes (AN) for species typical of this mountain range of western South America; Guiana Shield (GS), for species occurring in the lowlands of crystalline basement north of the Amazonas River; Guiana Shield highlands (GH), for species occurring mainly in the uplands of this mountainous complex; Potentially endemic to the Serra da Mocidade region (PE), for species potentially restricted to the studied mountain range. Finally, we also considered species that occur at punctual restricted localities outside of their main geographical range (PR) and at transition zones between Amazonia and other biomes (TZ). To evaluate the contribution of these regions in shaping the local herpetofauna, we calculated a relative percentage, dividing the number of species from each biogeographical region in relation to the total number of recorded species, for both amphibians and reptiles.
To verify the contribution of altitudinal generalists and specialists to the species composition of the Serra da Mocidade herpetofauna, we also classified species according to their known altitudinal range as reported in the literature, as lowland (occuring mainly below 500 m asl) and upland (occuring mainly above 500 m asl) species (Hoogmoed 1979, McDiarmid andDonnely 2005). The relative contribution of these assemblages was also evaluated by dividing the number of species from each group in relation to the total number of recorded species, for both amphibians and reptiles, considering the total sampling and samples from distinct altitudinal levels.

Results
A total of 305 specimens (232 amphibians, 58 squamates, 13 chelonians, and two crocodilians) was recorded, belonging to 51 species (23 amphibians, 24 squamates, three chelonians, and one crocodilian) from 25 families ( Table 1). The most diverse Table 1.   families were Hylidae for amphibians (six species), Dactyloidae for lizards (three species), Colubridae for snakes (six species), and Chelidae for chelonians (two species). Some examples of this diversity and the sampled habitats appear in Figs 3-8. As expected, the species rarefaction curves show that the Serra da Mocidade mountain range still has potential to harbour a greater amphibian and reptile diversity than we record in this short-term sampling (Fig. 9a). Furthermore, the species richness differs between sampled altitudinal levels, tending to decrease in the higher altitudes (above 900 m asl) (Fig. 9b).
The number of specimens recorded during sampling was similar in both camps (149 at base camp #1 vs. 158 at base camp #2), with some species exclusively recorded at distinct altitudinal levels, as with frogs of the genus Leptodactylus Fitzinger, 1826 only recorded at 600 m asl or the snakes Micrurus remotus Roze, 1987 and Chironius septentrionalis Dixon, Wiest & Cei, 1993 only recorded at 1,060 m asl (Fig. 10). Regarding sampling methods, the colubrid Atractus riveroi Roze, 1961 and the gymnophthalmid Cercosaura ocellata Wagler, 1830 were exclusively recorded using pitfall traps, a method that recorded a total of 15 specimens. Active surveys recorded 252 specimens, eight specimens were collected using shotguns, three using glue traps, ten using trammel nets and hook with bait, and 17 by occasional encounters, including all three caecilians.

Diversity and biogeographical affinities
The species richness of Serra da Mocidade is comparable to that recorded for other mountain ranges in the Guiana Shield highlands region, which harbour a low number of species compared to Guiana Shield lowland forest habitats (Table 2). Regarding the species composition, the amphibians and reptile assemblages of the Serra da Mocidade were evidenced in a intermediate position between the Guiana Shield highlands and lowland forest assemblages, which are distinct from the cluster generated by the lowland open habitat assemblages (R² = 0.74, stress = 0.14). We note a high divergence in species com-      position of the mountain ranges included in the analyses, because many of these mountains harbour several endemic species. Among these mountain ranges, the herpetofauna of Serra da Mocidade is more similar to that occurring at the geographically close Apiaú and Neblina, as well as to that occurs at the uplands of Nassau and Lely plateaus, which reach lower altitudes than the others mountain ranges included in this analyses (Fig. 9c).
This same mixed influence is corroborated regarding biogeographical regions, as the composition of the amphibian assemblage was most strongly influenced by widely distributed Amazonian forest lineages (47%), but also included species restricted to, or   typical for, the Guiana Shield (22%), Guiana Shield highlands (8%) and Eastern Amazonia (4%). The reptile assemblage composition showed a similar pattern, with most species widely distributed in Amazonia (75%) and smaller proportions of exclusive lineages from the Guiana Shield (7%), Guiana Shield highlands (7%), Eastern Amazonia (7%) and Western Amazonia (3.5%). Some of the taxa we collected are potentially endemic to the Serra da Mocidade mountain range, such as Stefania sp., Epicrionops sp., Brasilotyphlus sp. and Hyalinobatrachium aff. taylori (Table 1).
Most of the amphibian species we recorded on Serra da Mocidade occur across wide altitudinal ranges along their geographical distribution (43%). Several species are typical from the Amazonian lowlands (35%), and some are typical of the uplands (22%). The reptilian assemblage showed a similar altitudinal pattern: most of the Serra da Mocidade species are generalists in terms of altitudinal range (53.5%) or lowland forms (43%), and few are typical of the uplands (3.5%). When we analysed the distinct assemblages recorded at 600 m asl and above 900 m asl, we found different distribution patterns at a finer scale: at 600 m asl, assemblages were mainly composed of lowland species (47% of amphibians and 52% of reptiles) and altitudinal generalists (42% of amphibians and 48% of reptiles). Above 900 m asl, no elements from lowland forests were recorded for either taxonomic group and assemblages were composed of upland species (33% of amphibians and 11% of reptiles) and altitudinal generalists (67% of amphibians and 89% of reptiles) (Fig. 10).

Remarkable records
Accounts of the remarkable species are provided below, with comments on ecology, taxonomy, morphology, evolutionary distinctiveness, and biogeography.
Epicrionops sp. -The genus Epicrionops Boulenger, 1883 is distributed in mountainous habitats of the Andes (seven species) and the Guiana Shield highlands [Epicrionops niger (Dunn, 1942)] (Frost 2017). However, the generic position of E. niger is uncertain, as recent studies have shown this species to be genetically closer to species in the genus Rhinatrema (distributed mainly throughout the Guiana Shield) than to Andean species of Epicrionops Wiens 2011, San Mauro et al. 2014). The juvenile specimen recorded at Serra da Mocidade (Fig. 4h) is the first of the genus recorded from Brazil and has a high genetic distance to E. niger from Guyana (>15% on the 16S fragment used, Fig. 11). Considering such high molecular divergence and degree of geographical isolation, this taxon represents a new candidate species, needing further studies and samples for its formal description.
Brasilotyphlus sp. -The genera Microcaecilia Taylor, 1968 andBrasilotyphlus Taylor, 1968 (possibly synonymous by lack of diagnosability-see Maciel and Hoogmoed 2011) are distributed throughout eastern Amazonia, with highest diversity in the Guiana Shield. We found two caeciliid specimens at Serra da Mocidade, both at 1,060 m asl, possessing a diastema between palatine and vomerine teeths, the only known morphological characteristic diagnostic for the genus Brasilotyphlus (Maciel et al. 2009, Wilkinson et al. 2013). However, some morphological characteristic distinguish those specimens for other Brasilotyphlus, e.g. the lower number of primary annulus (less than 140 vs. more than 141 in the two known species). Although no Brasilotyphlus or Microcaecilia is known from the Brazilian side of the Guiana Shield highlands region, there is evidence that at least two undescribed species occur there, and one of them is in description process (Maciel and Hoogmoed 2011, Pedro Nunes pers. comm.). Our molecular analyses showed this taxon to be highly divergent from other Microcaecilia with available sequences in GenBank, with a minimum genetic distance of 17% on the 16S fragment used (Fig. 11). Based on these morphological and molecular results, we consider the Brasilotyphlus collected at Serra da Mocidade as a new candidate species.
Anomaloglossus apiau -The high endemism levels of the small cryptically coloured terrestrial frogs of the genus Anomaloglossus recorded from Guiana Shield highlands mountains , Vacher et al. 2017, together with the geographical isolation of Serra da Mocidade, led us to the immediately suspect that the population found in this mountain range was a new taxon. However, detailed analyses showed that their morphological (Fig. 3a,b), acoustic (Fig. 12, Table 3) and molecular variation (Fig. 13) are consistent with the description of A. apiau (Fouquet et al. 2015a), a species previously known only from its type locality, at Serra do Apiaú, a mountain range distant ca. 80 km northeast of Serra da Mocidade. Here we extend the known geographical range of this species. We found the species in all altitudinal levels sampled at Serra da Mocidade, near rocky streams with rapids and waterfalls (Fig. 8d, e). The population of Serra da Mocidade had adult males with 15-20 mm SVL, and females slightly larger than originally described for A. apiau (19-23 mm SVL) (see Fouquet et al. 2015a). As in the population from Serra do Apiaú, specimens from Serra da Mocidade had a high intraspecific polymorphism in colour pattern, and strong sexual dichromatism (Fig. 3a, b). The advertisement call has the same temporal and spectral structure as reported in the species description (long series of paired notes, followed by intervals of silence) (Fig. 12, Table 3). Additionally, the population from Serra da Mocidade occurs within the altitudinal range cited in the original description for the type locality (500-1,400 m asl). Molecular analyses confirmed the similarity between the two populations, with the sample of A. apiau from the type locality nested within the Serra da Mocidade clade (Fig. 13). Two subclades from Serra da Mocidade are separated by a low genetic distance (maximum 2% on the 16S fragment used) and occur at different altitudes and drainages. However, this difference is most likely due to natural intraspecific variation, since the genetic distance between populations of A. apiau from the type locality and Serra da Mocidade is also below 2%. This record of A. apiau from the Serra da Mocidade (first outside the type locality) is remarkable, since the species was not recorded from the mountain range of the Serra da Maroquinha, located at less than 5 km southeast from the Serra do Apiaú. New studies on the cryptic diversity, evolution and biogeography of the genus Anomaloglossus (such as Vacher et al. 2017) should reveal the species limits, geographical patterns and drivers of the diversification in low-and uplands of the Guiana Shield.
Rhinella martyi -Rhinella martyi has a confusing taxonomic history, as do many of the small-medium terrestrial forest bufonids in the Rhinella margaritifera group (dos Santos et al. 2015). Until recently, this taxon was allocated in the problematic taxon Rhinella margaritifera (Laurenti, 1768), but Fouquet et al. (2007a) recognized the presence of distinct lineages occurring in the Guiana Shield. Two of these lineages were later described as the new species R. lescurei and R. martyi (Fouquet et al. 2007b), but the specific status of the latter was rejected by Ávila-Pires et al.
. However, awaiting further integrative revisions of the group, the name R. martyi was retained (Frost 2017). Both the neotype designated to Rhinella margaritifera (Lavilla et al. 2013) and the holotype recently rediscovered (Lavilla et al. 2017), have a different morphology from populations named as R. martyi (e.g., by having larger cranial crests and body size). Apparently, R. martyi is widely distributed along the northern Guiana Shield, occurring along a broad altitudinal range (Fouquet et al. 2007b). Individuals of R. martyi were found on leaf-litter inside primary forests, and calling males on the banks of streams, mostly at 1,060 m asl (Fig. 3d). Morphological variation among specimens from the Serra da Mocidade mostly fits with the variation proposed in the original species' description, except in a lower variation in adult body size (at Serra da Mocidade, males 50-60 mm SVL, females 62-70 mm SVL), smaller height of supratympanic crests and the bony knob at angle of jaws less developed. The advertisement call recorded at Serra da Mocidade also differs slightly from the species description, with longer pulses and higher dominant frequency (Fig. 12, Table 3). However, the overall morphological and acoustic variation in R. martyi is still being clarified (A. Fouquet pers. comm.). Elucidation of the taxonomic status of populations of Rhinella margaritifera group from the Guiana Shield still depends on broader integrative revisions.
Pristimantis aff. vilarsi -Pristimantis is one of the most speciose genus of vertebrates (Padial et al. 2014, Frost 2017. Such high diversity is accompanied by a problematic taxonomy and difficulties in species delimitation (Padial et al. 2014). In the Guiana Shield highlands region, several lineages of Pristimantis diversified into lowlands and highland forms ). We found a species from the Pristimantis conspicillatus group at the Serra da Mocidade (Fig. 4g). It may represent a new taxon, but its taxonomic status is being evaluated in a broader sense in relation to other Pristimantis. Our molecular analyses showed that samples from the Serra da Mocidade are more similar to Pristimantis vilarsi (Melin, 1941, redescribed in Barrio-Amorós andMolina 2006) and Pristimantis zeuctotylus (Lynch & Hoogmoed, 1977), both species of the P. conspicillatus group from the lowlands and uplands of Guiana Shield. The population from Serra da Mocidade has morphological details that distinguish it from these two genetically-related species, as for example the adult body size (SVL): 22-33 mm in males and 31-49 mm in females of P. vilarsi, 20-29 mm in males and 30-43 mm in females of P. zeuctotylus (Barrio-Amorós and Molina 2006) and 36-46 mm in males and 48-57 mm in females from Serra da Mocidade. Pristimantis aff. vilarsi occurs across the altitudinal range sampled at Serra da Mocidade and was one of the most abundant anurans in our sampling (along with A. apiau). It was recorded on leaf-litter of primary dense forests, but also inhabits rocky outcrops within streams with fast flowing water (Fig. 8d, e).
Hyalinobatrachium aff. taylori -The species of Hyalinobatrachium from Serra da Mocidade (Fig. 3e) (Señaris and Ayarzagüena 2005, Guayasamin et al. 2009, Castroviejo-Fisher et al. 2011. It is possible that this difference represents an intraspecific variation, since the genetic distance on the 16S fragment used was less than 4% between both taxa (Fig. 14), but bone colouration is a strong character in the genus taxonomy, and green bones is a rare character in Hyalinobatrachium (two species), but common in other centrolenid genera, such as Vitreorana and Centrolene (Guayasamin et al. 2009). In addition, Serra da Mocidade is located far from the known geographical range of H. taylori, and the occurrence of this species at this locality would be the first record of the species in Brazil. The diversity of Hyalinobatrachium in the Guiana Shield highlands region appears underestimated, since some lineages were recently described as new species after an integrative revision (Castroviejo-Fisher et al. 2011), and the evolution of lineages in the H. taylori complex seem strongly influenced and structured by the altitudinal levels in which they occur, with clades from lowlands and highlands separated by low genetic distances (Castroviejo-Fisher et al. 2011). Considering the uncertainties, and waiting for the results of future detailed studies, we opted to keep this taxon as H. aff. taylori.
Vitreorana ritae -The small glassfrog Vitreorana ritae is apparently widely distributed in the Guiana Shield (Guayasamin et al. 2009), occupying a wide altitudinal range (see Cisneros-Heredia 2013 for taxonomical accounts and synonymization of V. oyampiensis with V. ritae). However, single-site records in western Amazonia and south of the Amazon River cloud the delineation of the total species' range and possible connections between populations. The species was found at Serra da Mocidade (Fig. 3f ) in riparian vegetation at three altitudinal levels, but with greatest abundance at 1,060 m asl. Specimens from Serra da Mocidade are differentiated from Vitreorana helenae (Ayarzagüena, 1992), a morphologically similar taxon from the Guiana Shield highlands region that has yellow eyes and a lime-green dorsum (eyes predominantly gray and a darker green dorsum in V. ritae) (Guayasamin et al. 2008). Molecular data also grouped sequences from the Serra da Mocidade with sequences of V. ritae from other Guiana Shield localities with less than 2% of genetic distance on the 16S fragment used, while the distinction between V. ritae and V. helenae lies between 4 and 5% (Fig. 15).
Stefania sp. -The marsupial frog genus Stefania have an evolutionary history intrinsically linked to the evolution of the Guiana Shield highlands landscape, as the  genus has a geographical range restricted to this region and high endemism levels in distinct mountain ranges (Duellman and Hoogmoed 1984, Salerno et al. 2012, Duellman 2015, Kok et al. 2016, 2017. We found two adult Stefania at 960 m asl on rocky outcrops in rapidly flowing streams (Fig. 8d, e). The specimens from Serra da Mocidade (Fig. 4c) have a large body size (52-54 mm SVL), head wider than long (head length / width = 0.93-0.94), snout rounded in dorsal and lateral views, shagreneed dorsum, basal toe webbing, colour pattern in life dorsally yellowishbrown with large rectangular blotches, ventrally gray (lighter in the abdomen) with large whitish patches, and large t-shaped finger and toe discs. The combination of these morphological characters distinguishes these specimens from other congeners (see variation in Rivero 1970, Myers and Donnelly 1997, MacCulloch and Lathrop 2002, Carvalho et al. 2010). The most similar species considering the external morphology is also the one with the closest known geographical distribution (Stefania tamacuarina Myers & Donnelly, 1997, which occurs ca. 300 km distant to the west of the Serra da Mocidade) (Caramaschi and Niemeyer 2005), however these taxa slightly differ in length and shape of snout, dorsal and ventral colouration pattern, as well as size of eye and hands (Myers and Donnelly 1997). The morphological divergence compared to the known congeners is sustained by a high genetic distance (>10% on the 16S gene fragment used), even considering other undescribed taxa detected in a recent phylogenetic analyses of the genus   (Fig. 16). The population of Stefania sp. from Serra da Mocidade is apparently isolated from known ranges of both described and undescribed species, and when confirmed as a new species, will be the third known Stefania endemic to a granitic mountain, while other taxa in this genus mainly occur at summits of sandstone mountains (Myers and Donnelly 1997). While new specimen samples (especially adults of S. tamacuarina) and DNA sequences are not available, we opted to keep this taxon as a candidate species.
Dendropsophus minutus -The species complex of small arboreal Dendropsophus minutus hylids is widely distributed in South America east of Andes, occupying a broad altitudinal range (Frost 2017). The evolutionary history of this complex has been recently revised, showing multiple distinct lineages throughout its wide distribution (Gehara et al. 2014). Several distinct cryptic lineages occur in the Guiana Shield (already noted by Hawkins et al. 2007), which suggests the elevation of the synonymized taxon Dendropsophus goughi [species described by Boulenger (1911) and synonymized to D. minutus by Lutz (1973)] to specific level for at least one of these lineages. While new studies do not clarify this taxonomic issue, the available name for the population recorded in the Serra da Mocidade is D. minutus (Fig. 3g). The specimens recorded share the presence of distinct light stripes in the cloacal region and in some specimens in the heels, and emitted mainly territorial calls (type B sensu Cardoso and Haddad 1984), with acoustic parameters within the known variation for the species (Fig. 12, Table 3).
Dendropsophus parviceps -Recent molecular studies on the evolution of the genus Dendropsophus found a polyphyletic Dendropsophus parviceps species group sensu Faivovich et al. 2005(Fouquet et al. 2011, 2015b, Motta et al. 2012, which harbours small cryptically coloured treefrogs. Recently, a new species of this group (Dendropsophus counani Fouquet, Orrico, Ernst, Blanc, Martinez, Vacher, Rodrigues, Ouboter, Jairam & Ron 2015) historically misidentified as D. parviceps, was described from the Guiana Shield (Fouquet et al. 2015b). We compared the specimens collected in Serra da Mocidade (Fig. 3h) with known taxa of the D. parviceps group, and both morphologic variation and molecular similarity (Fig. 17) reveals their identity as D. parviceps. We found a genetic distance of 3% on the 16S fragment used between samples from Serra da Mocidade and the D. parviceps sequence available in GenBank, from southwestern Amazonia (Acre state, Brazil), and this phylogenetic divergence can most likely be attributed to intraspecific variation due to wide geographical distance. This species has an ex-  tensive distribution throughout Amazonia, including historical records in other mountain ranges of Guayan Highlands Paolillo 1988, Schlüter and, and present high levels of genetic divergence reported between some populations (Fouquet et al. 2015b). It is likely that future studies will reveal other independent lineages hidden under this name, as was the case of D. counani.
Boana boans, B. multifasciata, B. xerophylla and Osteocephalus taurinus -These medium to large sized arboreal hylids have wide geographical ranges in Amazonia: B. boans (Fig. 3i) and O. taurinus (Fig. 4b) occur throughout almost the entire basin, B. xerophylla (Fig. 3j) occurs on the northern Guiana Shield and B. multifasciata (Fig. 4a) in eastern Amazonia and in the transition zone between Amazonia and Cerrado, an open savanna biome with its core distribution in central Brazil (Frost 2017). All these taxa are considered as lowland species, but they may occur in lower densities at higher altitudinal levels, having been recorded on other Guiana Shield highlands mountain ranges (Aubrecht et al. 2012 Fig. 12, Table 3). Previous studies investigating the intraspecific variation of these taxa have been conducted (e.g., Jungfer et al. 2013, Guarnizo et al. 2016, Orrico et al. 2017), but a clear definition of specific limits still needs further definition.
Adenomera andreae and Leptodactylus petersii -These species are both terrestrial frogs, typical of Amazonian lowlands (Sá et al. 2014, Frost 2017. The intraspecific molecular variation for A. andreae, known to be widely distributed in forest habitats of Amazonia, shows divergent lineages endemic to the Guiana Shield (Fouquet et al. 2014). As for L. petersii, several studies have investigated the morphological and bioacoustics variation of the species and their close relatives (Leptodactylus melanonotus species group) (Heyer 1970, 1994, de Sá et al. 2014), but the high polymorphism and intraspecific variation of species from this group still hinder definitive taxonomic assignment. Further studies with broader datasets may reveal restricted lineages within L. petersii.
Lithobates palmipes -This large ranid inhabit forest habitats near to slow-flowing water bodies (lakes, ponds and rivers), is widely distributed in Amazonia, Atlantic Forest and transitional habitats and may represent a species complex. The presence of several gaps in the known species distribution hampers a precise geographical determination (La . Although the expected occurrence of the species reaches Roraima state and it has already been recorded in nearby areas, no state record has appeared in the literature until this study (nearest documented record at 400 km west of the Serra da Mocidade in Rodrigues et al. 2013). Lithobates palmipes abundance also varies considerably across its distributional range (Acosta-Galvis 1999), but the environmental determinants of such variation are still unknown (Ramalho et al. 2011). In the foothills of Serra da Mocidade, the species was abundant in the Pacú River floodplain (Fig. 8b), a tributary of the Branco River, with more than 11 specimens (adults and larvae) recorded in three sampling nights.
Pseudogonatodes guianensis -Pseudogonatodes guianensis is a small leaf-litter lizard widely distributed within Amazonian lowlands (Ávila-Pires 1995, Ribeiro-Júnior 2015b). However, it is a species apparently rare in the Guiana Shield highlands region, where Coleodactylus septentrionalis Vanzolini, 1980 is the more abundant sphaerodactylid (Ribeiro-Júnior 2015b). No record of this species from the Roraima state appears in the previous literature, and our finding filling this distribution gap.
Plica plica -This widely distributed arboreal lizard (Ávila-Pires 1995, Rivas et al. 2012 has several distinct lineages with restricted geographical ranges (Murphy andJowers 2013, Oliveira et al. 2016). However, the genetic divergence between lineages is not clearly reflected in the morphologic variation. A recent study on the morphological variation of the Plica plica complex led to the description of several new taxa with restricted distributions and kept P. plica as the species with broadest distribution within the Guiana Shield (Murphy and Jowers 2013). Also, a study on the intraspecific genetic variation of this species showed at least two distinct lineages from the region of Serra da Mocidade (Oliveira et al. 2016). We recorded P. plica (Fig. 5h) at three altitudinal ranges, and considering such current state of knowledge for the species complex, we retain the specimens from the Serra da Mocidade under this epithet, based in their external morphology.
Mabuya nigropunctata -Specimens of Mabuya from Serra da Mocidade (Fig. 5g) were assigned to this name based in morphological characters within the known variation of the species, as for example, the paired prefrontal scales, two pairs of nuchal scales, five supraciliar scales, dorsals scales tricarinate, dark ventral surfaces of hand and feet, which are covered by small tubercles, and a dark lateral band, not limited by dorsal and ventral light stripe. This morphological variation promptly differ this specimens from other Mabuya with a geographical distribution known for this region: Mabuya carvalhoi Rebouças-Spieker and Vanzolini, 1990 (with fused prefrontals, large granules in ventral surface of hands and feets, three to five pairs of nuchals, five longitudinal light stripes along the body and a blue tail) and Mabuya bistriata (Spix, 1825) (with four supraciliars, dorsals smooth, one pair of nuchals, distinct light stripes limiting the dark lateral band and at dorsum, and light ventral surfaces of hands, which are covered by moderately large granules). However, the recorded specimens have an interesting characteristic that differs from the species known morphology (Ávila-Pires 1995, Hedges andConn 2012): the fusion of frontoparietal scales in a single butterfly-shaped scale. Morphological comparisons with other M. nigropunctata specimens collected in the Brazilian Guiana Shield highlands region (on Pico da Neblina, deposited at INPA-H collection) reveal the same pattern on the frontoparietal scales for some specimens, while other specimens from the same locality have the typical pattern of two frontoparietals.
Intraspecific analyses of molecular variation indicated a strong genetic structure and multiple lineages within this widely distributed Amazonian taxon (Miralles et al. 2005, Miralles andCarranza 2010). One of those lineages (which have Mabuya surinamensis Hallowell, 1857 as available name) occurs on the Guiana and Brazilian Shields, including the Guiana Shield highlands region (Miralles and Carranza 2010). The dissimilarity in the number of frontoparietals found in the specimens from Northern Brazil may represent typical variation within this lineage, but confirmation awaits further taxonomic studies.
Tretioscincus oriximinensis -The small cryptic specimens of Tretioscincus from Serra da Mocidade (Fig. 5d) were collected on leaf-litter of primary dense forest and have morphological characters within the known diagnosis of T. oriximinensis: high number of dorsal scales (>30 rows), polygonal scales on tail, in 12 rows and keeled, prefrontal in contact, six gular scales anteriorly, dorsolateral light stripe become paler at middle of the dorsum (Ávila-Pires 1995). In addition, molecular analyses showed low genetic distance on the 16S fragment used (<0.1%) between Serra da Mocidade specimens and a sample of T. oriximinensis from Venezuela. Although the species is apparently more abundant in open habitats, the original description also cites some individuals found in forested areas (Ávila-Pires 1995).
There is a subtle morphological divergence between T. oriximinensis populations from northern (including Serra da Mocidade) and eastern Amazonia (Ávila-Pires 1995), with differences in number of ventral scales and extension of keels in scales from tail to posterior dorsals. This morphological divergence is reflected at the molecular level (Fig. 18), as we found ca. 4% of genetic distance between samples from this populations. Dissimilarities between the T. oriximinensis populations indicate that their taxonomic status deserves to be further investigated.
Atractus riveroi -This groundsnake typical from Guiana Shield highlands was previously known by only two specimens (Roze 1961, Passos et al. 2013. We found three specimens in pitfall traps at the 1,060 m asl altitudinal level: two with a brown dorsal background with dark markings (Fig. 6d) and one with a black dorsal background with small white spots (Fig. 6e) (Passos et al. 2013). For detailed information on morphologic variation and geographical distribution of this species based in these records and new specimens from nearby mountain ranges, see Fraga et al. (2017).
Chironius fuscus and C. septentrionalis -Two species of diurnal Chironius snakes were recorded at Serra da Mocidade: one is a typical lowland species widely distributed in Amazonia and other ecosystems in South America (Chironius fuscus - Fig. 6b), while the other is a upland inhabitant (Chironius septentrionalis - Fig. 6a, adult male recorded only at 1,060 m asl, with 1,480 mm SVL, 350 mm caudal length -tail damaged, dorsals 12/12/8, ventrals 179, subcaudals 60, anal plate divided, apical pits on neck scales) (Dixon et al. 1993). The latter species was described as a subspecies of the widely distributed Amazonian species Chironius multiventris Schmidt andWalker, 1943 (Dixon et al. 1993), but a morphologic taxonomic revision elevated it to specific level (Dixon et al. 1993). This species occurs at high altitudinal levels in adjacent Venezuela, but may occur in lowlands on the island of Trinidad (Dixon et al. 1993). This is the first record from Brazil, extending the known distribution by more than 500 km south.
Drymobius rhombifer -Despite being widely distributed in Amazonia (Rivas et al. 2012), this diurnal snake (Fig. 6c) is rarely recorded, apparently due to low densities throughout its range (Stafford and Castro 2010 Micrurus remotus -The holotype of this small monadal coral snake is from the Guiana Shield highlands region [Cerro de la Neblina (Roze 1987)]. Additional specimens have been recorded in this region and in southern Amazonia (Roze 1987, Bernarde et al. 2011, Bernarde and Gomes 2012. Given the great geographical distance and environmental dissimilarity between these two localities, and the overall difficulty in species delimitation due to high intraspecific polymorphism in the species group (Feitosa et al. 2015), a taxonomic revision is required to clarify their status. The specimen from Serra da Mocidade (Fig. 7b) is one of the few documented for the Roraima state (see plate 172 in Campbell and Lamar 2004) and differs slightly in the body colour from patterns described in the literature: it has two light spots on the rostrum, greater extent of white on the head and the white ring after the black nuchal collar is almost imperceptible.
Platemys platycephala melanonota -Platemys platycephala (Schneider, 1792) is a solitary and nocturnal chelid, which inhabits shallow temporary pools within Amazonian lowland rainforests (Vogt 2008). Several years after their original description, a melanistic colour form was described as a distinct subspecies, possessing a large amount of black pigmentation in the carapace and head and some differences in head scalation in comparison to nominal subspecies (Mendes-Pinto et al. 2012).
As this subspecies is more rarely recorded, the total geographical distribution is unknown, but there are records from Peru, Ecuador and Brazil (Mendes-Pinto et al. 2012, van Dijk et al. 2014, which suggest a wide distribution and low abundances throughout the range. The specimen from Serra da Mocidade (Fig. 7g) is the first record of the subspecies in the Roraima state, more than 900 km from the nearest known location (Mendes-Pinto et al. 2012).

Discussion
Several remarkable herpetofaunal records were found during the first large biological expedition conducted at the Serra da Mocidade mountain range. Based on our results, it is evident that the herpetofauna inhabitant of this mountain range has a greater biotic affinity with lowlands from the Amazon region, but some elements typical of the uplands from Guiana Shield highlands region occurs above 900 m asl. However, the local species diversity is certainly underestimated in our results, especially at higher altitudinal zones, as this mountain range has an extensive unexplored area of upland forests, which may harbour populations of undescribed amphibians and reptiles typical of the Guiana Shield highlands region. Additionally, the sampling period covered the dry season in the region. This, linked to the fact that the strong El Niño event of 2015-2016 produced higher temperatures throughout the year (Varotsos et al. 2016), may have potentially limited encounter rates of some herpetological groups.
Integrative approaches are increasingly being used in biological inventories (e.g., Vieites et al. 2009, Jansen et al. 2011, Paz and Crawford 2012, Moraes et al. 2016, and are based on several proposed methods (e.g., Padial et al. 2010). The integrative identification approach used here to reveal the species diversity at Serra da Mocidade illustrates the relevance of using different evolutionary data sources to identify taxa from remote and unexplored Amazonian areas. Replication of this method for future Amazonian biodiversity inventories will certainly contribute to a more precise evaluation of species diversity and distribution, as well as the origin, diversification drivers and conservation status of such species.
With the molecular approach, based on reciprocal monophyly, high nodal support and genetic distances of mtDNA, we detected additional samples and extended the distribution of known lineages (A. apiau, R. martyi, and T. oriximinensis, with genetic distances less than 2%), discovered putative divergent lineages of known species (H. aff. taylori, P. aff. vilarsi, V. ritae, and D. parviceps, with genetic distances between 3-6%), as well as candidate new species (Stefania sp., Epicrionops sp. and Brasilotyphlus sp., with genetic distances between 10-17%). The taxonomic status of the putative new species detected in this study needs to be confirmed using a broader and more detailed analyses of data sources, which may increase the known species diversity at Serra da Mocidade.
Several attempts have been made to define biogeographical sub-regions within Guiana Shield highlands, mainly based on assemblage similarities among isolated mountain ranges (e.g., Mayr and Phelps 1967, Huber 1988a, McDiarmid and Donnelly 2005. Despite differences in their overall boundaries, all such studies seem to agree on the existence of two main sub-biogeographical regions with high similarity in species composition: the western and eastern Guiana Shield highlands (Mayr and Phelps 1967). This distinction is also observed in the geological origin of regions, as each sub-region is distinct in age and landscape history (Santos et al. 2003). Western Guiana Shield highlands may have experienced a greater faunal exchange during Pleistocene glacial events because the mountains involved lie on a higher basement matrix than Eastern Guiana Shield highlands (Mayr and Phelps 1967, Cook 1974, Cracraft 1985. At Serra da Mocidade, the herpetofauna composition at highest altitudes was predominantly composed of widespread and altitudinal generalists, with greater similarity to the western sub-region. Occurrence of upland species known to occur in other Guiana Shield highlands mountain ranges (e.g., A. apiau, A. riveroi) and species closely related to other Guiana Shield highlands endemics (e.g., Stefania sp., Epicrionops sp.) reinforce the postulated recent high connectivity between highland assemblages of distinct mountain ranges in this region , Salerno et al. 2012. Increased gene flow may have occurred during Pleistocene glacial events (Noonan and Gaucher 2005, Lötters et al. 2010) which made the lowland matrix more permeable, leading to recent events of colonization and speciation. Phylogeographical studies and more extensive sampling are required to address these possibilities more explicitly.
Some species that are often altitudinal-generalists were only found on the Serra da Mocidade at altitudes higher than 900 m (e.g., D. minutus, D. parviceps, B. xerophylla, A. punctatus, B. b. bilineatus). The Serra da Mocidade mountain range has a complex hydrologic mosaic, with streams that vary in their basin origins and amounts of dissolved sediments (Barbosa 2005, Ministério do Meio Ambiente 2016), which generates high habitat heterogeneity. The absence of some altitudinal generalists at lower altitudes (<900 m asl) may be related to a balance between physiological tolerance and habitat conditions, or to sampling bias. Long-term sampling would be necessary to fully elucidate species distribution patterns along these altitudinal gradients.
The results from this short-term multidisciplinary expedition (see other scientific results in Bastos et al. 2016, Neto et al. 2016, Dantas and Hamada 2017, Fraga et al. 2017, Lourenço 2017, Raimundi et al. 2017, and the documentary associated with the expedition at http://www.grifafilmes.com/en/new-species) highlight our poor knowledge of the Brazilian Guiana Shield highlands region. Our study can be used as primary data source for future biodiversity, biogeography, and conservation assessments that consider Guiana Shield highlands on a larger scale. We encourage additional initiatives to enable costly and logistically difficult expeditions to remote Amazonian areas, to fill-in basic knowledge gaps of biological diversity in remote areas, and investigate the processes that led to the currently observed biodiversity patterns.