Disentangling the drivers of continental mammalian endemism

Endemic species and species with small ranges are ecologically and evolutionarily distinct and are vulnerable to extinction. Determining which abiotic and biotic factors structure patterns of endemism on continents can advance our understanding of global biogeographic processes, but spatial patterns of mammalian endemism have not yet been effectively predicted and reconstructed. Using novel null model techniques, we reconstruct trends in mammalian endemism and describe the isolated and combined effects of physiographic, ecological, and evolutionary factors on endemism. We calculated weighted endemism for global continental ecoregions and compared the spatial distribution of endemism to niche-based, geographic null models of endemism. These null models distribute species randomly across continents, simulating their range sizes from their degree of climatic specialization. They isolate the effects of physiography (topography and climate) and species richness on endemism. We then ran linear and structural models to determine how topography and historical climate stability influence endemism. The highest rates of mammalian endemism were found in topographically rough, climatically stable ecoregions with many species. The null model that isolated physiography did not closely approximate the observed distribution of endemism (r² = .09), whereas the null model that incorporated both physiography and species richness did (r² = .59). The linear models demonstrate that topography and climatic stability both influenced endemism values, but that average climatic niche breadth was not highly correlated with endemism. Climate stability and topography both influence weighted endemism in mammals, but the spatial distribution of mammalian endemism is driven by a combination of physiography and species richness. Despite its relationship to individual range size, average climate niche breadth has only a weak influence on endemism. The results highlight the importance of historical biogeographic processes (e.g. centers of speciation) and geography in driving endemism patterns, and disentangle the mechanisms structuring species ranges worldwide.     

The relation between the neutrality index for mitochondrial genes and the distribution of mutational effects on fitness

We explore factors affecting patterns of polymorphism and divergence (as captured by the neutrality index) at mammalian mitochondrial loci. To do this, we develop a population genetic model that incorporates a fraction of neutral amino acid sites, mutational bias, and a probability distribution of selection coefficients against new nonsynonymous mutations. We confirm, by reanalyzing publicly available datasets, that the mitochondrial cyt-b gene shows a broad range of neutrality indices across mammalian taxa, and explore the biological factors that can explain this observation. We find that observed patterns of differences in the neutrality index, polymorphism, and divergence are not caused by differences in mutational bias. They can, however, be explained by a combination of a small fraction of neutral amino acid sites, weak selection acting on most amino acid mutations, and differences in effective population size among taxa.     

Spatial requirements of jaguars and pumas in Southern Mexico

Understanding how large felids use space is essential for the design of conservation plans that are required for their survival. Jaguars (Panthera onca) and pumas (Puma concolor) are the largest felids in the Neotropics, and they are sympatric throughout the entire range of the jaguar. However, there is very little information about the spatial requirements of these two species in the tropical rainforests of Central America. Using satellite GPS collars, we compared the spatial ecology of jaguars and pumas in a tropical rainforest in southern Mexico. We found that jaguars had home ranges that were 2–6 times larger than those of pumas. The mean annual home range was 181.4 ± 4.0 km² for female jaguars and 431.6 ± 152.6 km² for males. Annual home range for the only female puma tracked was 34.3 km², and 72.0 ± 85.2 km² for males. Jaguars and pumas with overlapping home ranges showed little overlap of core areas and avoided using the same sites at the same time, which suggested that territoriality and avoidance were in play. Further evidence of avoidance was derived from our observation that pumas exhibited greater movement during the lightest periods of the day and jaguars moved most during the darkest. This temporal separation probably facilitates coexistence. Our data suggest that habitat destruction and fragmentation has more severe effects on jaguars than on pumas. According to our estimates, patches of at least 180 km² of primary forest are required to meet the annual spatial requirements of female jaguars. Thus, to conserve jaguars in this region, large tracts of primary forest should be preserved. Importantly, this population of jaguars depends on the adequate preservation of connectivity between natural reserves in Mexico and Guatemala.     

Levels of endemism are not necessarily biased by the co-presence of species with different range sizes: a case study of Vilenkin and Chikatunov's models

Aim  To illustrate problems in the methods proposed by B. Vilenkin and V. Chikatunov to study levels of endemism and species–area relationships.
Location  The study used data on the distribution of tenebrionid beetles (Coleoptera, Tenebrionidae) on the Aegean Islands (Greece).
Methods  A total of 32 islands and 170 taxa (species and subspecies) were included in this study. Levels of endemism were evaluated both as the proportion of endemic taxa, and according to the methods proposed by Vilenkin and Chikatunov, which are based on the number of non-endemic taxa and various relationships with area. A model of the species–area relationship proposed by these authors was also analysed.
Results  The number of endemic taxa was positively correlated with the number of taxa with different distribution types, but this positive correlation did not influence the estimation of the level of endemism. In fact, the commonly used estimate of endemicity as a percentage was strongly correlated with the endemism values calculated according to the method of Vilenkin and Chikatunov. The usual power function fitted the species–area relationship as well as the most complicated method of Vilenkin and Chikatunov.
Main conclusions  As hypothesized by Vilenkin and Chikatunov, the number of endemic taxa was influenced both by the number of taxa of other biogeographical ranks, and by an island's area. However, explanations for the positive relationship between the number of endemic taxa and taxa of different biogeographical ranks are equivocal. Importantly, this relationship did not necessarily influence the level of endemism, which could be expressed adequately by percentages. The method proposed by Vilenkin and Chikatunov to estimate the species–area relationship cannot be clearly justified on theoretical grounds and is of questionable practical utility.     

Evaluation of the enteric microflora of captive whooping cranes (Grus americana) and sandhill cranes (Grus canadensis)

The enteric flora of captive whooping cranes (Grus americana) and sandhill cranes (Grus canadensis) has not been well described, despite its potential importance in the understanding of both the normal condition of the intestinal physiology of these animals and the altered colonization within disease states in these birds. Nineteen whooping cranes and 23 sandhill cranes housed currently at the Calgary Zoo or its affiliated Devonian Wildlife Conservation Centre (DWCC) in Calgary, Alberta were sampled from October 2004–February 2005 by collecting aerobic and anaerobic cloacal swabs from each bird. There were seven major groupings of bacteria isolated from both species of crane. Gram‐positive cocci, coliforms, and gram‐negative bacilli were the most prevalent types of bacteria isolated for both crane species, with Escherichia coli, Enterococcus faecalis, and Streptococcus Group D, not Enterococcus the bacterial species isolated most commonly. There was a significant difference in the average number of isolates per individual between the two crane species but no differences between age or gender categories within crane species. Campylobacter sp. were isolated from five whooping cranes. The potential zoonotic pathogen Campylobacter jejuni was isolated from one whooping crane and C. upsaliensis was isolated from a second. Three other isolates were unspeciated members of the Campylobacter genus and likely belong to a species undescribed previously. The evaluation of the enteric cloacal flora of whooping cranes and sandhill cranes illustrates that differences exist between these two closely related crane species, and highlights the potential implications these differences may have for current practices involving captive wildlife. Zoo Biol 0:1–13, 2007. © 2007 Wiley‐Liss, Inc.     

Quantitative delineation of how breathing motions open ligand migration channels in myoglobin and its mutants

Ligand migration and binding are central to the biological functions of many proteins such as myoglobin (Mb) and it is widely thought that protein breathing motions open up ligand channels dynamically. However, how a protein exerts its control over the opening and closing of these channels through its intrinsic dynamics is not fully understood. Specifically, a quantitative delineation of the breathing motions that are needed to open ligand channels is lacking. In this work, we present and apply a novel normal mode‐based method to quantitatively delineate what and how breathing motions open ligand migration channels in Mb and its mutants. The motivation behind this work springs from the observation that normal mode motions are closely linked to the breathing motions that are thought to open ligand migration channels. In addition, the method provides a direct and detailed depiction of the motions of each and every residue that lines a channel and can identify key residues that play a dominating role in regulating the channel. The all‐atom model and the full force‐field employed in the method provide a realistic energetics on the work cost required to open a channel, and as a result, the method can be used to efficiently study the effects of mutations on ligand migration channels and on ligand entry rates. Our results on Mb and its mutants are in excellent agreement with MD simulation results and experimentally determined ligand entry rates. Proteins 2015; 83:757–770. © 2015 Wiley Periodicals, Inc.     

The performance of fine‐grained and coarse‐grained elastic network models and its dependence on various factors

In a recent work we developed a method for deriving accurate simplified models that capture the essentials of conventional all‐atom NMA and identified two best simplified models: ssNMA and eANM, both of which have a significantly higher correlation with NMA in mean square fluctuation calculations than existing elastic network models such as ANM and ANMr2, a variant of ANM that uses the inverse of the squared separation distances as spring constants. Here, we examine closely how the performance of these elastic network models depends on various factors, namely, the presence of hydrogen atoms in the model, the quality of input structures, and the effect of crystal packing. The study reveals the strengths and limitations of these models. Our results indicate that ssNMA and eANM are the best fine‐grained elastic network models but their performance is sensitive to the quality of input structures. When the quality of input structures is poor, ANMr2 is a good alternative for computing mean‐square fluctuations while ANM model is a good alternative for obtaining normal modes. Proteins 2015; 83:1273–1283. © 2015 Wiley Periodicals, Inc.     

Vibrational entropy differences between mesophile and thermophile proteins and their use in protein engineering

We recently introduced ENCoM, an elastic network atomic contact model, as the first coarse-grained normal mode analysis method that accounts for the nature of amino acids and can predict the effect of mutations on thermostability based on changes vibrational entropy. In this proof-of-concept article, we use pairs of mesophile and thermophile homolog proteins with identical structures to determine if a measure of vibrational entropy based on normal mode analysis can discriminate thermophile from mesophile proteins. We observe that in around 60% of cases, thermophile proteins are more rigid at equivalent temperatures than their mesophile counterpart and this difference can guide the design of proteins to increase their thermostability through series of mutations. We observe that mutations separating thermophile proteins from their mesophile orthologs contribute independently to a decrease in vibrational entropy and discuss the application and implications of this methodology to protein engineering.