Quantifying similarity in animal vocal sequences: Which metric performs best?

1. Many animals communicate using sequences of discrete acoustic elements which can be complex, vary in their degree of stereotypy, and are potentially open-ended. Variation in sequences can provide important ecological, behavioural, or evolutionary information about the structure and connectivity of populations, mechanisms for vocal cultural evolution, and the underlying drivers responsible for these processes. Various mathematical techniques have been used to form a realistic approximation of sequence similarity for such tasks. 2. Here, we use both simulated and empirical datasets from animal vocal sequences (rock hyrax, Procavia capensis; humpback whale, Megaptera novaeangliae; bottlenose dolphin, Tursiops truncatus; and Carolina chickadee, Poecile carolinensis) to test which of eight sequence analysis metrics are more likely to reconstruct the information encoded in the sequences, and to test the fidelity of estimation of model parameters, when the sequences are assumed to conform to particular statistical models. 3. Results from the simulated data indicated that multiple metrics were equally successful in reconstructing the information encoded in the sequences of simulated individuals (Markov chains, n-gram models, repeat distribution, and edit distance), and data generated by different stochastic processes (entropy rate and n-grams). However, the string edit (Levenshtein) distance performed consistently and significantly better than all other tested metrics (including entropy, Markov chains, n-grams, mutual information) for all empirical datasets, despite being less commonly used in the field of animal acoustic communication. 4. The Levenshtein distance metric provides a robust analytical approach that should be considered in the comparison of animal acoustic sequences in preference to other commonly employed techniques (such as Markov chains, hidden Markov models, or Shannon entropy). The recent discovery that non-Markovian vocal sequences may be more common in animal communication than previously thought, provides a rich area for future research that requires non-Markovian based analysis techniques to investigate animal grammars and potentially the origin of human language.We thank Melinda Rekdahl, Todd Freeberg and his graduate students, Amiyaal Ilany, Elizabeth Hobson, and Jessica Crance for providing comments of on a previous version of this manuscript. We thank Mike Noad, Melinda Rekdahl, and Claire Garrigue for assistance with humpback whale song collection and initial categorisation of the song, Vincent Janik and Laela Sayigh for assistance with signature whistle collection, Todd Freeberg with chickadee recordings, and Eli Geffen and Amiyaal Ilany for assistance with hyrax song collection and analysis. E.C.G is supported by a Newton International Fellowship. Part of this work was conducted while E.C.G. was supported by a National Research Council (National Academy of Sciences) Postdoctoral Fellowship at the National Marine Mammal Laboratory, AFSC, NMFS, NOAA. The findings and conclusions in this paper are those of the authors and do not necessarily represent the views of the National Marine Fisheries Service. We would also like to thank Randall Wells and the Sarasota Dolphin Research Program for the opportunity to record the Sarasota dolphins, where data were collected under a series of National Marine Fisheries Service Scientific Research Permits issued to Randall Wells. A.K. is supported by the Herchel Smith Postdoctoral Fellowship Fund. Part of this work was conducted while A.K. was a Postdoctoral Fellow at the National Institute for Mathematical and Biological Synthesis, an Institute sponsored by the National Science Foundation through NSF Award #DBI-1300426, with additional support from The University of Tennessee, Knoxville.This is the author accepted manuscript. The final version is available from Wiley via http://dx.doi.org/10.1111/2041-210X.1243

Niche partitioning by three Pterodroma petrel species during non-breeding in the equatorial Pacific Ocean

Niche divergence is expected for species that compete for shared resources, including migrants that occupy similar regions during the non-breeding season. Studies of temperate seabirds indicate that both spatial and behavioural segregation can be important mechanisms for reducing competition, but there have been few investigations of resource partitioning by closely related taxa in low productivity, tropical environments. We investigated niche partitioning in 3 gadfly petrel taxa, Pterodroma leucoptera leucoptera (n = 22), P. leucoptera caledonica (n = 7) and P. pycrofti (n = 12), during their non-breeding season in the eastern tropical Pacific Ocean by combining tracking data from geolocator-immersion loggers with remotely sensed environmental data in species distribution models (SDMs), and by comparing feather stable isotope ratios. The 3 taxa showed spatial partitioning: two foraged in the North Equatorial Counter Current and one in the South Equatorial Current. This reflected differences in their realised habitat niches, with significant taxon-specific responses to thermocline depth, sea surface temperature and bathymetry. There were also differences among taxa in activity patterns, and all birds spent a much larger proportion of time in flight at night than during the day, suggesting predominance of nocturnal foraging behaviour. Comparison of stable isotope ratios in feathers suggests that P. l. leucoptera and P. pycrofti mainly consume vertically migrating mesopelagic fishes, whereas the diet of P. l. caledonica also includes some lower trophic levels including crustaceans and squid. Unique insights can be gained from studies of the foraging ecology of tropical pelagic seabirds, in comparison with temperate and polar waters, and are urgently required for understanding and protecting tropical avifauna in key marine habitats

Population genomics of a predatory mammal reveals patterns of decline and impacts of exposure to toxic toads

Mammal declines across northern Australia are one of the major biodiversity loss events occurring globally. There has been no regional assessment of the implications of these species declines for genomic diversity. To address this, we conducted a species-wide assessment of genomic diversity in the northern quoll (Dasyurus hallucatus), an Endangered marsupial carnivore. We used next generation sequencing methods to genotype 10,191 single nucleotide polymorphisms (SNPs) in 352 individuals from across a 3220-km length of the continent, investigating patterns of population genomic structure and diversity, and identifying loci showing signals of putative selection. We found strong heterogeneity in the distribution of genomic diversity across the continent, characterized by (i) biogeographical barriers driving hierarchical population structure through long-term isolation, and (ii) severe reductions in diversity resulting from population declines, exacerbated by the spread of introduced toxic cane toads (Rhinella marina). These results warn of a large ongoing loss of genomic diversity and associated adaptive capacity as mammals decline across northern Australia. Encouragingly, populations of the northern quoll established on toad-free islands by translocations appear to have maintained most of the initial genomic diversity after 16 years. By mapping patterns of genomic diversity within and among populations, and investigating these patterns in the context of population declines, we can provide conservation managers with data critical to informed decision-making. This includes the identification of populations that are candidates for genetic management, the importance of remnant island and insurance/translocated populations for the conservation of genetic diversity, and the characterization of putative evolutionarily significant units

GEANT4 : a simulation toolkit

Abstract Geant4 is a toolkit for simulating the passage of particles through matter. It includes a complete range of functionality including tracking, geometry, physics models and hits. The physics processes offered cover a comprehensive range, including electromagnetic, hadronic and optical processes, a large set of long-lived particles, materials and elements, over a wide energy range starting, in some cases, from 250 eV and extending in others to the TeV energy range. It has been designed and constructed to expose the physics models utilised, to handle complex geometries, and to enable its easy adaptation for optimal use in different sets of applications. The toolkit is the result of a worldwide collaboration of physicists and software engineers. It has been created exploiting software engineering and object-oriented technology and implemented in the C++ programming language. It has been used in applications in particle physics, nuclear physics, accelerator design, space engineering and medical physics. PACS: 07.05.Tp; 13; 2

Search for CP Violation in the Decay Z -> b (b bar) g

About three million hadronic decays of the Z collected by ALEPH in the years 1991-1994 are used to search for anomalous CP violation beyond the Standard Model in the decay Z -> b \bar{b} g. The study is performed by analyzing angular correlations between the two quarks and the gluon in three-jet events and by measuring the differential two-jet rate. No signal of CP violation is found. For the combinations of anomalous CP violating couplings, ${\hat{h}}_b = {\hat{h}}_{Ab}g_{Vb}-{\hat{h}}_{Vb}g_{Ab}$ and $h^{\ast}_b = \sqrt{\hat{h}_{Vb}^{2}+\hat{h}_{Ab}^{2}}$, limits of \hat{h}_b < 0.59$and$h^{\ast}_{b} < 3.02$ are given at 95\% CL.Comment: 8 pages, 1 postscript figure, uses here.sty, epsfig.st

Search for supersymmetry with a dominant R-parity violating LQDbar couplings in e+e- collisions at centre-of-mass energies of 130GeV to 172 GeV

A search for pair-production of supersymmetric particles under the assumption that R-parity is violated via a dominant LQDbar coupling has been performed using the data collected by ALEPH at centre-of-mass energies of 130-172 GeV. The observed candidate events in the data are in agreement with the Standard Model expectation. This result is translated into lower limits on the masses of charginos, neutralinos, sleptons, sneutrinos and squarks. For instance, for m_0=500 GeV/c^2 and tan(beta)=sqrt(2) charginos with masses smaller than 81 GeV/c^2 and neutralinos with masses smaller than 29 GeV/c^2 are excluded at the 95% confidence level for any generation structure of the LQDbar coupling.Comment: 32 pages, 30 figure

Search for the glueball candidates f0(1500) and fJ(1710) in gamma gamma collisions

Data taken with the ALEPH detector at LEP1 have been used to search for gamma gamma production of the glueball candidates f0(1500) and fJ(1710) via their decay to pi+pi-. No signal is observed and upper limits to the product of gamma gamma width and pi+pi- branching ratio of the f0(1500) and the fJ(1710) have been measured to be Gamma_(gamma gamma -> f0(1500)). BR(f0(1500)->pi+pi-) < 0.31 keV and Gamma_(gamma gamma -> fJ(1710)). BR(fJ(1710)->pi+pi-) < 0.55 keV at 95% confidence level.Comment: 10 pages, 3 figure