No link between population isolation and speciation rate in squamate reptiles.

Rates of species formation vary widely across the tree of life and contribute to massive disparities in species richness among clades. This variation can emerge from differences in metapopulation-level processes that affect the rates at which lineages diverge, persist, and evolve reproductive barriers and ecological differentiation. For example, populations that evolve reproductive barriers quickly should form new species at faster rates than populations that acquire reproductive barriers more slowly. This expectation implicitly links microevolutionary processes (the evolution of populations) and macroevolutionary patterns (the profound disparity in speciation rate across taxa). Here, leveraging extensive field sampling from the Neotropical Cerrado biome in a biogeographically controlled natural experiment, we test the role of an important microevolutionary process-the propensity for population isolation-as a control on speciation rate in lizards and snakes. By quantifying population genomic structure across a set of codistributed taxa with extensive and phylogenetically independent variation in speciation rate, we show that broad-scale patterns of species formation are decoupled from demographic and genetic processes that promote the formation of population isolates. Population isolation is likely a critical stage of speciation for many taxa, but our results suggest that interspecific variability in the propensity for isolation has little influence on speciation rates. These results suggest that other stages of speciation-including the rate at which reproductive barriers evolve and the extent to which newly formed populations persist-are likely to play a larger role than population isolation in controlling speciation rate variation in squamates.

Biochemical characterization of a novel sphingomyelinase-like protein from the Rhipicephalus microplus tick.

Sphingomyelinase D is a toxin present in venomous spiders and bacteria and is associated with infection symptoms in patients affected by spider bites. It was observed that in Ixodes scapularis ticks, sphingomyelinase-like protein secreted in saliva can modulate the host immune response, affecting the transmission of flavivirus to the host via exosomes. In this work, a sphingomyelinase D-like protein (RmSMase) from R. microplus, a tick responsible for economic losses and a vector of pathogens for cattle, was investigated. The amino acid sequence revealed the lack of important residues for enzymatic activity, but the recombinant protein showed sphingomyelinase D activity. RmSMase shows Ca2+ and Mg2+ dependence in acidic pH, differing from IsSMase, which has Mg2+ dependence in neutral pH. Due to the difference between RmSMase and other SMases described, the data suggest that RmSMase belongs to SMase D class IIc. RmSMase mRNA transcription levels are upregulated during tick feeding, and the recombinant protein was recognized by host antibodies elicited after heavy tick infestation, indicating that RmSMase is present in tick saliva and may play a role in the tick feeding process.

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective.

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

Amphibian Speciation Rates Support a General Role of Mountains as Biodiversity Pumps.

AbstractContinental mountain areas cover <15% of global land surface, yet these regions concentrate >80% of global terrestrial diversity. One prominent hypothesis to explain this pattern proposes that high mountain diversities could be explained by higher diversification rates in regions of high topographic complexity (HTC). While high speciation in mountains has been detected for particular clades and regions, the global extent to which lineages experience faster speciation in mountains remains unknown. Here we addressed this issue using amphibians as a model system (>7,000 species), and we found that families showing high speciation rates contain a high proportion of species distributed in mountains. Moreover, we found that lineages inhabiting areas of HTC speciate faster than lineages occupying areas that are topographically less complex. When comparing across regions, we identified the same pattern in five biogeographical realms where higher speciation rates are associated with higher levels of complex topography. Low-magnitude differences in speciation rates between some low and high complex topographies suggest that high mountain diversity is also affected by low extinction and/or high colonization rates. Nevertheless, our results bolster the importance of mountains as engines of speciation at different geographical scales and highlight their importance for the conservation of global biodiversity.

Idiosyncratic responses to drivers of genetic differentiation in the complex landscapes of Isthmian Central America.

Isthmian Central America (ICA) is one of the most biodiverse regions in the world, hosting an exceptionally high number of species per unit area. ICA was formed <25 million years ago and, consequently, its biotic assemblage is relatively young and derived from both colonization and in situ diversification. Despite intensive taxonomic work on the local fauna, the potential forces driving genetic divergences and ultimately speciation in ICA remain poorly studied. Here, we used a landscape genetics approach to test whether isolation by distance, topography, habitat suitability, or environment drive the genetic diversity of the regional frog assemblage. To this end, we combined data on landscape features and mitochondrial DNA sequence variation for nine codistributed amphibian species with disparate life histories. In five species, we found that at least one of the factors tested explained patterns of genetic divergence. However, rather than finding a general pattern, our results revealed idiosyncratic responses to historical and ecological processes, indicating that intrinsic life-history characteristics may determine the effect of different drivers of isolation on genetic divergence in ICA. Our work also suggests that the convergence of several factors promoting isolation among populations over a heterogeneous landscape might maximize genetic differentiation, despite short geographical distances. In conclusion, abiotic factors and geographical features have differentially affected the genetic diversity across the regional frog assemblage. Much more complex models (i.e., considering multiple drivers), beyond simple vicariance of Caribbean and Pacific lineages, are needed to better understand the evolutionary history of ICA's diverse biotas.

Global priorities for conservation across multiple dimensions of mammalian diversity.

Conservation priorities that are based on species distribution, endemism, and vulnerability may underrepresent biologically unique species as well as their functional roles and evolutionary histories. To ensure that priorities are biologically comprehensive, multiple dimensions of diversity must be considered. Further, understanding how the different dimensions relate to one another spatially is important for conservation prioritization, but the relationship remains poorly understood. Here, we use spatial conservation planning to (i) identify and compare priority regions for global mammal conservation across three key dimensions of biodiversity-taxonomic, phylogenetic, and traits-and (ii) determine the overlap of these regions with the locations of threatened species and existing protected areas. We show that priority areas for mammal conservation exhibit low overlap across the three dimensions, highlighting the need for an integrative approach for biodiversity conservation. Additionally, currently protected areas poorly represent the three dimensions of mammalian biodiversity. We identify areas of high conservation priority among and across the dimensions that should receive special attention for expanding the global protected area network. These high-priority areas, combined with areas of high priority for other taxonomic groups and with social, economic, and political considerations, provide a biological foundation for future conservation planning efforts.

Cetartiodactyla: Updating a time-calibrated molecular phylogeny.

Cetartiodactyla comprises one of the most diverse mammal radiations. Currently, 23 families, 131 genera and more than 330 species are recognized. Several studies have been trying to resolve its phylogenetic relationships. The most comprehensive dated phylogenetic hypothesis available includes only 55% of the extant species, precluding a clear understanding of ecological and evolutionary patterns in Cetartiodactyla. Here, we gathered all mitochondrial genetic data available in GenBank to build a robust Cetartiodactyla calibrated phylogenetic tree using 21 fossil calibration points. We found mitogenomic data for 225 species and included other 93 species from which there was at least one mitochondrial gene available. Using a Bayesian approach, we generated a dated tree comprising 90% of the extant Cetartiodactyla species (n = 318). The major lineages showed robust support and families divergence times are congruent with the available fossil evidence and with previously published phylogenetic hypotheses. By making available a dated phylogeny with extensively sampled clades, we expect to foster future studies on the origin, tempo and mode of Cetartiodactyla diversification.

Evolutionary time drives global tetrapod diversity.

Global variation in species richness is widely recognized, but the explanation for what drives it continues to be debated. Previous efforts have focused on a subset of potential drivers, including evolutionary rate, evolutionary time (maximum clade age of species restricted to a region), dispersal (migration from one region to another), ecological factors and climatic stability. However, no study has evaluated these competing hypotheses simultaneously at a broad spatial scale. Here, we examine their relative contribution in determining the richness of the most comprehensive dataset of tetrapods to our knowledge (84% of the described species), distinguishing between the direct influences of evolutionary rate, evolutionary time and dispersal, and the indirect influences of ecological factors and climatic stability through their effect on direct factors. We found that evolutionary time exerted a primary influence on species richness, with evolutionary rate being of secondary importance. By contrast, dispersal did not significantly affect richness patterns. Ecological and climatic stability factors influenced species richness indirectly by modifying evolutionary time (i.e. persistence time) and rate. Overall, our findings suggest that global heterogeneity in tetrapod richness is explained primarily by the length of time species have had to diversify.

Estimating synchronous demographic changes across populations using hABC and its application for a herpetological community from northeastern Brazil.

Many studies propose that Quaternary climatic cycles contracted and/or expanded the ranges of species and biomes. Strong expansion-contraction dynamics of biomes presume concerted demographic changes of associated fauna. The analysis of temporal concordance of demographic changes can be used to test the influence of Quaternary climate on diversification processes. Hierarchical approximate Bayesian computation (hABC) is a powerful and flexible approach that models genetic data from multiple species, and can be used to estimate the temporal concordance of demographic processes. Using available single-locus data, we can now perform large-scale analyses, both in terms of number of species and geographic scope. Here, we first compared the power of four alternative hABC models for a collection of single-locus data. We found that the model incorporating an a priori hypothesis about the timing of simultaneous demographic change had the best performance. Second, we applied the hABC models to a data set of seven squamate and four amphibian species occurring in the Seasonally Dry Tropical Forests (Caatinga) in northeastern Brazil, which, according to paleoclimatic evidence, experienced an increase in aridity during the Pleistocene. If this increase was important for the diversification of associated xeric-adapted species, simultaneous population expansions should be evident at the community level. We found a strong signal of synchronous population expansion in the Late Pleistocene, supporting the increase of the Caatinga during this time. This expansion likely enhanced the formation of communities adapted to high aridity and seasonality and caused regional extirpation of taxa adapted to wet forest.

Life-History Patterns of Lizards of the World.

Identification of mechanisms that promote variation in life-history traits is critical to understand the evolution of divergent reproductive strategies. Here we compiled a large life-history data set (674 lizard populations, representing 297 species from 263 sites globally) to test a number of hypotheses regarding the evolution of life-history traits in lizards. We found significant phylogenetic signal in most life-history traits, although phylogenetic signal was not particularly high. Climatic variables influenced the evolution of many traits, with clutch frequency being positively related to precipitation and clutches of tropical lizards being smaller than those of temperate species. This result supports the hypothesis that in tropical and less seasonal climates, many lizards tend to reproduce repeatedly throughout the season, producing smaller clutches during each reproductive episode. Our analysis also supported the hypothesis that viviparity has evolved in lizards as a response to cooler climates. Finally, we also found that variation in trait values explained by clade membership is unevenly distributed among lizard clades, with basal clades and a few younger clades showing the most variation. Our global analyses are largely consistent with life-history theory and previous results based on smaller and scattered data sets, suggesting that these patterns are remarkably consistent across geographic and taxonomic scales.

Global mammal beta diversity shows parallel assemblage structure in similar but isolated environments.

The taxonomic, phylogenetic and trait dimensions of beta diversity each provide us unique insights into the importance of historical isolation and environmental conditions in shaping global diversity. These three dimensions should, in general, be positively correlated. However, if similar environmental conditions filter species with similar trait values, then assemblages located in similar environmental conditions, but separated by large dispersal barriers, may show high taxonomic, high phylogenetic, but low trait beta diversity. Conversely, we expect lower phylogenetic diversity, but higher trait biodiversity among assemblages that are connected but are in differing environmental conditions. We calculated all pairwise comparisons of approximately 110 × 110 km grid cells across the globe for more than 5000 mammal species (approx. 70 million comparisons). We considered realms as units representing geographical distance and historical isolation and biomes as units with similar environmental conditions. While beta diversity dimensions were generally correlated, we highlight geographical regions of decoupling among beta diversity dimensions. Our analysis shows that assemblages from tropical forests in different realms had low trait dissimilarity while phylogenetic beta diversity was significantly higher than expected, suggesting potential convergent evolution. Low trait beta diversity was surprisingly not found between isolated deserts, despite harsh environmental conditions. Overall, our results provide evidence for parallel assemblage structure of mammal assemblages driven by environmental conditions at a global scale.

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome.

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.

Evaluation of the biological function of ribosomal protein S18 from cattle tick Rhipicephalus microplus.

Rhipicephalus (Boophilus) microplus, also known as the cattle tick, causes severe parasitism and transmits different pathogens to vertebrate hosts, leading to massive economic losses. In the present study, we performed a functional characterization of a ribosomal protein from R. microplus to investigate its importance in blood feeding, egg production and viability. Ribosomal protein S18 (RPS18) is part of the 40S subunit, associated with 18S rRNA, and has been previously pointed to have a secondary role in different organisms. Rhipicephalus microplus RPS18 (RmRPS18) gene expression levels were modulated in female salivary glands during blood feeding. Moreover, mRNA levels in this tissue were 10 times higher than those in the midgut of fully engorged female ticks. Additionally, recombinant RmRPS18 was recognized by IgG antibodies from sera of cattle naturally or experimentally infested with ticks. RNAi-mediated knockdown of the RmRPS18 gene was performed in fully engorged females, leading to a significant (29 %) decrease in egg production. Additionally, egg hatching was completely impaired, suggesting that no viable eggs were produced by the RmRPS18-silenced group. Furthermore, antimicrobial assays revealed inhibitory activities against gram-negative Escherichia coli and gram-positive Staphylococcus aureus bacteria, affecting bacterial growth. Data presented here show the important role of RmRPS18 in tick physiology and suggest that RmRPS18 can be a potential target for the development of novel strategies for tick control.

Gut membrane proteins as candidate antigens for immunization of mice against the tick Amblyomma sculptum.

Amblyomma sculptum is widely distributed in Brazil and is the main vector of Rickettsia rickettsii, the causative agent of the Brazilian spotted fever (BSF). Tick gut proteins play an essential role in blood feeding, digestion, and protection of gut epithelium. Therefore, many of these were investigated as potential vaccine targets for tick-control strategies. The present study aimed to select transcripts corresponding to putative immunogenic proteins in the A. sculptum gut epithelial membrane, produce recombinant proteins and evaluate them as antigens against A. sculptum infestations. Three gut proteins - AsMucin, AsAPP, and AsLAMP - and a chimeric protein (rAsChimera) based on 22 peptides containing putative B cell epitopes from seven different gut proteins were evaluated as anti-A. sculptum antigens. Mice immunizations revealed that all recombinant targets elicited humoral response with significantly increased IgG levels compared to controls. For rAsChimera, IgG levels remained significantly higher than controls up to 75 days after the end of the immunization. Challenge trials revealed that vaccination with the chimeric protein was the most effective against A. sculptum, inducing 100 % nymph mortality and reaching 80.8 % efficacy against females. The other three proteins did not induce relevant protection, as AsAPP had only 26.6 % efficacy, whereas AsMucin and AsLAMP induced no protection. These data indicate that targeting gut protein immunogenic regions may be an effective strategy for a vaccine formulation againstA. sculptum.

The macroevolutionary singularity of snakes.

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Predator size, prey size, and dietary niche breadth relationships in marine predators.

Based on geographical and home range sizes, physiology, and gape limitation, a positive relationship between predator size and diet breadth is expected. Alternatively, larger predators might avoid smaller prey; in this case no relationship would be found. Here, I used a large data set on the diets of marine predators to describe and identify mechanisms responsible for the relationships among predator body size, diet breadth, and the mean, minimum, maximum, and variance of prey size. I found no relationship between predator size and diet breadth. Mean, minimum, maximum, and variance of prey size were all positively associated with predator size. I found that larger predators increase their minimum and maximum prey size with similar slopes, which explains the lack of relationship between predator size and diet breadth. The results support predictions of the hypothesis that optimal foraging is the main factor constraining the shape of the relationships among predator size, prey size, and diet breadth. Future research should focus on examining the relationship between body size and the breadth of different niche axis across different groups of organisms to assess whether a positive relationship between body size and niche breadth is a general rule in macroecology.

Occurrence of the Parasitic Fly Philornis torquans on Fledglings of the Rufous-Fronted Thornbird ( Phacellodomus rufifrons) in Southeast Brazil.

Philornis is a neotropical genus of muscid fly that interacts with birds and may affect the development and survival of the birds' offspring. Although Philornis is a relatively common parasite, there is a lack of information about Philornis hosts in several parts of the Americas. In this study, two nests of the Rufousfronted Thornbird ( Phacellodomus rufifrons) were collected in Pedro Leopoldo, southeast Brazil. The first contained four nestlings of advanced age (about 20 d old) and a recently emerged Philornis torquans female adult fly. The second nest contained three nestlings (less than 7 d old) and several Philornis torquans subcutaneous larvae. One of the nestlings was infested by 53 larvae, which had attacked several parts of its body and caused individual wounds containing 1 to more than 15 larvae. The length of the larvae ranged from 3 to 18 mm and only one was a second instar; the remaining 69 were third instars. The pupal period lasted 9-13 d. In total, 71 larvae were collected from the nest, with nestling parasitism varying from 7 to 53 larvae (mean- 23.7±25.5 larvae/nestling).

Species diversity as a surrogate for conservation of phylogenetic and functional diversity in terrestrial vertebrates across the Americas.

Preserving the evolutionary history and ecological functions that different species embody, in addition to species themselves, is a growing concern for conservation. Recent studies warn that conservation priority regions identified using species diversity differ from those based on phylogenetic or functional diversity. However, spatial mismatches in conservation priority regions need not indicate low surrogacy among these dimensions in conservation planning. Here, we use data for 10,213 terrestrial vertebrate species across the Americas to evaluate surrogacy; that is, the proportion of phylogenetic or functional diversity represented in conservation plans targeting species. We find that most conservation plans targeting species diversity also represent phylogenetic and functional diversity well, despite spatial mismatches in the priority regions identified by each plan. However, not all phylogenetic and functional diversity is represented within species-based plans, with the highest-surrogacy conservation strategy depending on the proportion of land area included in plans. Our results indicate that targeting species diversity could be sufficient to preserve much of the phylogenetic and functional dimensions of biodiversity in terrestrial vertebrates of the Americas. Incorporating phylogenetic and functional data in broad-scale conservation planning may not always be necessary, especially when the cost of doing so is high.

RNA-seq analysis of the salivary glands and midgut of the Argasid tick Ornithodoros rostratus.

Ornithodoros rostratus is a South American argasid tick which importance relies on its itchy bite and potential as disease vector. They feed on a wide variety of hosts and secrete different molecules in their saliva and intestinal content that counteract host defences and help to accommodate and metabolize the relatively large quantity of blood upon feeding. The present work describes the transcriptome profile of salivary gland (SG) and midgut (MG) of O. rostratus using Illumina sequencing. A total of 8,031 contigs were assembled and assigned to different functional classes. Secreted proteins were the most abundant in the SG and accounted for ~67% of all expressed transcripts with contigs with identity to lipocalins and acid tail proteins being the most representative. On the other hand, immunity genes were upregulated in MG with a predominance of defensins and lysozymes. Only 10 transcripts in SG and 8 in MG represented ~30% of all RNA expressed in each tissue and one single contig (the acid tail protein ORN-9707) represented ~7% of all expressed contigs in SG. Results highlight the functional difference of each organ and identified the most expressed classes and contigs of O. rostratus SG and MG.

Niche expansion and the niche variation hypothesis: does the degree of individual variation increase in depauperate assemblages?

The niche expansion and niche variation hypotheses predict that release from interspecific competition will promote niche expansion in depauperate assemblages. Niche expansion can occur by different mechanisms, including an increase in within-individual, among-individual, or bimodal variation (sexual dimorphism). Here we explore whether populations with larger niche breadth have a higher degree of diet variation. We also test whether populations from depauperate lizard assemblages differ in dietary resource use with respect to variation within and/or among individuals and sexual dimorphism. We found support for the niche expansion and niche variation hypotheses. Populations in assemblages with low phylogenetic diversity had a higher degree of individual variation, suggesting a tendency for niche expansion. We also found evidence suggesting that the mechanism causing niche expansion is an increase in variation among individuals rather than an increase in within-individual variation or an increase in bimodal variation due to sexual dimorphism.