Toadlets of the genus Brachycephalus are endemic to the Atlantic rainforests of southeastern and southern Brazil. The 14 species currently described have snout-vent lengths less than 18 mm and are thought to have evolved through miniaturization: an evolutionary process leading to an extremely small adult body size. Here, we present the first comprehensive phylogenetic analysis for Brachycephalus, using a multilocus approach based on two nuclear (Rag-1 and Tyr) and three mitochondrial (Cyt b, 12S, and 16S rRNA) gene regions. Phylogenetic relationships were inferred using a partitioned Bayesian analysis of concatenated sequences and the hierarchical Bayesian method (BEST) that estimates species trees based on the multispecies coalescent model. Individual gene trees showed conflict and also varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. The concatenated gene tree was completely resolved and is identical in topology and degree of statistical support to the individual mtDNA gene tree. On the other hand, the BEST species tree showed reduced significant node support relative to the concatenate tree and recovered a basal trichotomy, although some bipartitions were significantly supported at the tips of the species tree. Comparison of the log likelihoods for the concatenated and BEST trees suggests that the method implemented in BEST explains the multilocus data for Brachycephalus better than the Bayesian analysis of concatenated data. Landmark-based geometric morphometrics revealed marked variation in cranial shape between the species of Brachycephalus. In addition, a statistically significant association was demonstrated between variation in cranial shape and genetic distances estimated from the mtDNA and nuclear loci. Notably, B. ephippium and B. garbeana that are predicted to be sister-species in the individual and concatenated gene trees and the BEST species tree share an evolutionary novelty, the hyperossified dorsal plate.

No abstract

Species of the genus Brachycephalus, have a snout-vent length of less than 18 mm and are believed to have evolved through miniaturization. Brachycephalus ephippium, is particularly interesting; because its entire skull is hyperossified, and the presacral vertebrae and transverse processes are covered by a dorsal shield. We demonstrate in this paper that, at the macroscopic level, a completely hyperossified skull and dorsal shield occur only in B. ephippium, but not in B. ferruginus, B. izechsohni, B. pernix, B. pombali, B. brunneus, B. didactylus, and B. hermogenesi. An intermediate condition, in which the skull is hyperossified but a dorsal shield is absent, occurs in B. vertebralis, B. nodoterga, B. pitanga, and B. alipioi. The microscopic structure of hyperossification was examined in skulls of B. ephippium and B. pitanga, revealing a complex organization involving the presence of Sharpey fibers, which in humans are characteristic of periodontal connections.

Genetic monitoring using noninvasive samples provides a complement or alternative to traditional population monitoring methods. However, next-generation sequencing approaches to monitoring typically require high quality DNA and the use of noninvasive samples (e.g., scat) is often challenged by poor DNA quality and contamination by nontarget species. One promising solution is a highly multiplexed sequencing approach called genotyping-in-thousands by sequencing (GT-seq), which can enable cost-efficient genomics-based monitoring for populations based on noninvasively collected samples. Here, we develop and validate a GT-seq panel of 324 single nucleotide polymorphisms (SNPs) optimized for genotyping of polar bears based on DNA from noninvasively collected faecal samples. We demonstrate (1) successful GT-seq genotyping of DNA from a range of sample sources, including successful genotyping (>50% loci) of 62.9% of noninvasively collected faecal samples determined to contain polar bear DNA; and (2) that we can reliably differentiate individuals, ascertain sex, assess relatedness, and resolve population structure of Canadian polar bear subpopulations based on a GT-seq panel of 324 SNPs. Our GT-seq data reveal spatial-genetic patterns similar to previous polar bear studies but at lesser cost per sample and through use of noninvasively collected samples, indicating the potential of this approach for population monitoring. This GT-seq panel provides the foundation for a noninvasive toolkit for polar bear monitoring and can contribute to community-based programmes -a framework which may serve as a model for wildlife conservation and management for species worldwide.

No abstract