Convergence of gut microbiota in myrmecophagous amphibians

The gut microbiome composition of terrestrial vertebrates is known to converge in response to common specialized dietary strategies, like leaf-eating (folivory) or ant- and termite-eating (myrmecophagy). To date, such convergence has been studied in mammals and birds, but has been neglected in amphibians. Here, we analysed 15 anuran species (frogs and toads) representing five Neotropical families and demonstrated the compositional convergence of the gut microbiomes of distantly related myrmecophagous species. Specifically, we found that the gut microbial communities of bufonids and microhylids, which have independently evolved myrmecophagy, were significantly more similar than expected based on their hosts' evolutionary divergence. Conversely, we found that gut microbiome composition was significantly associated with host evolutionary history in some cases. For instance, the microbiome composition of Xenohyla truncata, one of the few known amphibians that eat fruits, was not different from those of closely related tree frogs with an arthropod generalist diet. Bacterial taxa overrepresented in myrmecophagous species relative to other host families include Paludibacter, Treponema, and Rikenellaceae, suggesting diet-mediated selection and prey-to-predator transmission likely driving the observed compositional convergence. This study provides a basis for examining the roles of the gut microbiome in host tolerance and sequestration of toxic alkaloids from ants and termites.

Host macrocyclic acylcarnitines mediate symbiotic interactions between frogs and their skin microbiome

The host-microbiome associations occurring on the skin of vertebrates significantly influence hosts' health. However, the factors mediating their interactions remain largely unknown. Herein, we used integrated technical and ecological frameworks to investigate the skin metabolites sustaining a beneficial symbiosis between tree frogs and bacteria. We characterize macrocyclic acylcarnitines as the major metabolites secreted by the frogs' skin and trace their origin to an enzymatic unbalance of carnitine palmitoyltransferases. We found that these compounds colocalize with bacteria on the skin surface and are mostly represented by members of the Pseudomonas community. We showed that Pseudomonas sp. MPFS isolated from frogs' skin can exploit acylcarnitines as its sole carbon and nitrogen source, and this metabolic capability is widespread in Pseudomonas. We summarize frogs' multiple mechanisms to filter environmental bacteria and highlight that acylcarnitines likely evolved for another function but were co-opted to provide nutritional benefits to the symbionts.

Habitat split as a driver of disease in amphibians

Anthropogenic habitat disturbance is fundamentally altering patterns of disease transmission and immunity across the vertebrate tree of life. Most studies linking anthropogenic habitat change and disease focus on habitat loss and fragmentation, but these processes often lead to a third process that is equally important: habitat split. Defined as spatial separation between the multiple classes of natural habitat that many vertebrate species require to complete their life cycles, habitat split has been linked to population declines in vertebrates, e.g. amphibians breeding in lowland aquatic habitats and overwintering in fragments of upland terrestrial vegetation. Here, we link habitat split to enhanced disease risk in amphibians (i) by reviewing the biotic and abiotic forces shaping elements of immunity and (ii) through a spatially oriented field study focused on tropical frogs. We propose a framework to investigate mechanisms by which habitat split influences disease risk in amphibians, focusing on three broad host factors linked to immunity: (i) composition of symbiotic microbial communities, (ii) immunogenetic variation, and (iii) stress hormone levels. Our review highlights the potential for habitat split to contribute to host-associated microbiome dysbiosis, reductions in immunogenetic repertoire, and chronic stress, that often facilitate pathogenic infections and disease in amphibians and other classes of vertebrates. We highlight that targeted habitat-restoration strategies aiming to connect multiple classes of natural habitats (e.g. terrestrial-freshwater, terrestrial-marine, marine-freshwater) could enhance priming of the vertebrate immune system through repeated low-load exposure to enzootic pathogens and reduced stress-induced immunosuppression.

Low microbiome diversity in threatened amphibians from two biodiversity hotspots

Microbial diversity positively influences community resilience of the host microbiome. However, extinction risk factors such as habitat specialization, narrow environmental tolerances, and exposure to anthropogenic disturbance may homogenize host-associated microbial communities critical for stress responses including disease defense. In a dataset containing 43 threatened and 90 non-threatened amphibian species across two biodiversity hotspots (Brazil's Atlantic Forest and Madagascar), we found that threatened host species carried lower skin bacterial diversity, after accounting for key environmental and host factors. The consistency of our findings across continents suggests the broad scale at which low bacteriome diversity may compromise pathogen defenses in species already burdened with the threat of extinction.

Draft genome and multi-tissue transcriptome assemblies of the Neotropical leaf-frog Phyllomedusa bahiana

Amphibians are increasingly threatened worldwide, but the availability of genomic resources that could be crucial for implementing informed conservation practices lags well behind that for other vertebrate groups. Here, we describe draft de novo genome, mitogenome, and transcriptome assemblies for the Neotropical leaf-frog Phyllomedusa bahiana native to the Brazilian Atlantic Forest and Caatinga. We used a combination of PacBio long reads and Illumina sequencing to produce a 4.74-Gbp contig-level genome assembly, which has a contiguity comparable to other recent nonchromosome level assemblies. The assembled mitogenome comprises 16,239 bp and the gene content and arrangement are similar to other Neobratrachia. RNA-sequencing from 8 tissues resulted in a highly complete (86.3%) reference transcriptome. We further use whole-genome resequencing data from P. bahiana and from its sister species Phyllomedusa burmeisteri, to demonstrate how our assembly can be used as a backbone for population genomics studies within the P. burmeisteri species group. Our assemblies thus represent important additions to the catalog of genomic resources available from amphibians.

A tale of two bellies: systematics of the oval frogs (Anura: Microhylidae: Elachistocleis Parker, 1927)

Oval frogs (Elachistocleis) have a broad geographic distribution covering nearly all of South America and parts of Central America. They also have a large inter- and intraspecific variation of the few morphological characters commonly used as diagnostic traits among species of the genus. Based on molecular data, we provide the most complete phylogeny of Elachistocleis to date, and explore its genetic diversity using distance-based and tree-based methods for putative species delimitation. Our results show that at least two of the most relevant traditional characters used in the taxonomy of this group (belly pattern and dorsal median white line) carry less phylogenetic information than previously thought. Based on our results, we propose some synonymizations and some candidate new species. This study is a first major step in disentangling the current systematics of Elachistocleis. Yet, a comprehensive review of morphological data is needed before any new species descriptions can be properly made.

Habitat and Host Species Drive the Structure of Bacterial Communities of Two Neotropical Trap-Jaw Odontomachus Ants : Habitat and Host Species Drive the Structure of Bacterial Communities of Two Neotropical Trap-Jaw Odontomachus Ants

Ants have long been known for their associations with other taxa, including macroscopic fungi and symbiotic bacteria. Recently, many ant species have had the composition and function of their bacterial communities investigated. Due to its behavioral and ecological diversity, the subfamily Ponerinae deserves more attention regarding its associated microbiota. Here, we used the V4 region of the 16S rRNA gene to characterize the bacterial communities of Odontomachus chelifer (ground-nesting) and Odontomachus hastatus (arboreal), two ponerine trap-jaw species commonly found in the Brazilian savanna ("Cerrado") and Atlantic rainforest. We investigated habitat effects (O. chelifer in the Cerrado and the Atlantic rainforest) and species-specific effects (both species in the Atlantic rainforest) on the bacterial communities' structure (composition and abundance) in two different body parts: cuticle and gaster. Bacterial communities differed in all populations studied. Cuticular communities were more diverse, while gaster communities presented variants common to other ants, including Wolbachia and Candidatus Tokpelaia hoelldoblerii. Odontomachus chelifer populations presented different communities in both body parts, highlighting the influence of habitat type. In the Atlantic rainforest, the outcome depended on the body part targeted. Cuticular communities were similar between species, reinforcing the habitat effect on bacterial communities, which are mainly composed of environmentally acquired taxa. Gaster communities, however, differed between the two Odontomachus species, suggesting species-specific effects and selective filters. Unclassified Firmicutes and uncultured Rhizobiales variants are the main components accounting for the observed differences. Our study indicates that both host species and habitat act synergistically, but to different degrees, to shape the bacterial communities in these Odontomachus species.

Signatures of functional bacteriome structure in a tropical direct-developing amphibian species

Background

Host microbiomes may differ under the same environmental conditions and these differences may influence susceptibility to infection. Amphibians are ideal for comparing microbiomes in the context of disease defense because hundreds of species face infection with the skin-invading microbe Batrachochytrium dendrobatidis (Bd), and species richness of host communities, including their skin bacteria (bacteriome), may be exceptionally high. We conducted a landscape-scale Bd survey of six co-occurring amphibian species in Brazil’s Atlantic Forest. To test the bacteriome as a driver of differential Bd prevalence, we compared bacteriome composition and co-occurrence network structure among the six focal host species.

Results

Intensive sampling yielded divergent Bd prevalence in two ecologically similar terrestrial-breeding species, a group with historically low Bd resistance. Specifically, we detected the highest Bd prevalence in Ischnocnema henselii but no Bd detections in Haddadus binotatus. Haddadus binotatus carried the highest bacteriome alpha and common core diversity, and a modular network partitioned by negative co-occurrences, characteristics associated with community stability and competitive interactions that could inhibit Bd colonization.

Conclusions

Our findings suggest that community structure of the bacteriome might drive Bd resistance in H. binotatus, which could guide microbiome manipulation as a conservation strategy to protect diverse radiations of direct-developing species from Bd-induced population collapses.

Formation of topological defects in nematic shells with a dumbbell-like shape

The present study investigates dumbbell-shaped nematic liquid crystal shells. Using molecular dynamics (MD) simulations, we consider the effects of an external electric field on nematic ordering by computing the average molecular alignment's time evolution and equilibrium configuration. We show that the number and location of topological defects are strongly affected by the external field, with the orientational ordering's equilibrium configuration depending on field direction about the shell's long axis. For a transverse external field, it is verified that the defect rearrangement presents a non-linear dynamics, with a field independent characteristic time scale delimiting the short and long time regimes. Effects associated with varying the shell's Gaussian curvature are also analyzed.

Molecular basis of a bacterial-amphibian symbiosis revealed by comparative genomics, modeling, and functional testing

The molecular bases for the symbiosis of the amphibian skin microbiome with its host are poorly understood. Here, we used the odor-producer Pseudomonas sp. MPFS and the treefrog Boana prasina as a model to explore bacterial genome determinants and the resulting mechanisms facilitating symbiosis. Pseudomonas sp. MPFS and its closest relatives, within a new clade of the P. fluoresens Group, have large genomes and were isolated from fishes and plants, suggesting environmental plasticity. We annotated 16 biosynthetic gene clusters from the complete genome sequence of this strain, including those encoding the synthesis of compounds with known antifungal activity and of odorous methoxypyrazines that likely mediate sexual interactions in Boana prasina. Comparative genomics of Pseudomonas also revealed that Pseudomonas sp. MPFS and its closest relatives have acquired specific resistance mechanisms against host antimicrobial peptides (AMPs), specifically two extra copies of a multidrug efflux pump and the same two-component regulatory systems known to trigger adaptive resistance to AMPs in P. aeruginosa. Subsequent molecular modeling indicated that these regulatory systems interact with an AMP identified in Boana prasina through the highly acidic surfaces of the proteins comprising their sensory domains. In agreement with a symbiotic relationship and a highly selective antibacterial function, this AMP did not inhibit the growth of Pseudomonas sp. MPFS but inhibited the growth of another Pseudomonas species and Escherichia coli in laboratory tests. This study provides deeper insights into the molecular interaction of the bacteria-amphibian symbiosis and highlights the role of specific adaptive resistance toward AMPs of the hosts.

Ecological divergence and synchronous Pleistocene diversification in the widespread South American butter frog complex

Phylogeographic studies primarily focus on the major role of landscape topography in driving lineage diversification. However, populational phylogeographic breaks may also occur as a result of either niche conservatism or divergence, in the absence of geographic barriers to gene flow. Furthermore, these two factors are not mutually exclusive and can act in concert, making it challenging to evaluate their relative importance on explaining genetic variation in nature. Herein, we use sequences of two mitochondrial and four nuclear genes to investigate the timing and diversification patterns of species pertaining to the Leptodactylus latrans complex, which harbors four morphologically cryptic species with broad distributions across environmental gradients in eastern South America. The origin of this species complex dates back to the late Miocene (ca. 5.5 Mya), but most diversification events occurred synchronically during the late Pleistocene likely as the result of ecological divergence driven by Quaternary climatic oscillations. Further, significant patterns of environmental niche divergences among species in the L. latrans complex imply that ecological isolation is the primary mode of genetic diversification, mostly because phylogenetic breaks are associated with environmental transitions rather than topographic barriers at both species and populational scale. We provided new insights about diversification patterns and processes within a species complex of broadly and continuously distributed group of frogs along South America.

From genes to ecosystems: a synthesis of amphibian biodiversity research in Brazil

Abstract Here, we summarize examples of significant advances in amphibian research supported by the São Paulo Research Foundation (FAPESP), focusing on recent discoveries in the fields of community ecology, habitat change, infection diseases, and multipurpose DNA sequencing. We demonstrated that FAPESP has been fundamental not only by directly funding research projects and scholarships, but also through its science training policy, fostering international collaborations with world-class research institutions, improving and consolidating new lines of research that often depended on a synergetic combination of different knowledge and complex tools. We emphasized that future studies will continue to focus on basic questions, such as description of new species, as well as taxonomic and systematic corrections. Furthermore, we also expect that there will be a strong integration among different disciplines using novel bioinformatics tools and modeling approaches, such as machine learning. These new approaches will be critical to further develop our understanding of foundational questions of amphibian life-history trait variation, disease transmission, community assembly, biogeography, and population forecasts under different global change scenarios such as agricultural expansion, agrochemical use, habitat loss, and climate change.

Universal dynamical scaling laws in three-state quantum walks

We perform a finite-time scaling analysis over the detrapping point of a three-state quantum walk on the line. The coin operator is parametrized by ρ that controls the wave packet spreading velocity. The input state prepared at the origin is set as a symmetric linear combination of two eigenstates of the coin operator with a characteristic mixing angle θ, one of them being the component that results in full spreading occurring at θ_{c}(ρ) for which no fraction of the wave packet remains trapped near the initial position. We show that relevant quantities, such as the survival probability and the participation ratio assume single parameter scaling forms at the vicinity of the detrapping angle θ_{c}. In particular, we show that the participation ratio grows linearly in time with a logarithmic correction, thus, shedding light on previous reports of sublinear behavior.

Initial Phylotranscriptomic Confirmation of Homoplastic Evolution of the Conspicuous Coloration and Bufoniform Morphology of Pumpkin-Toadlets in the Genus Brachycephalus

The genus Brachycephalus is a fascinating group of miniaturized anurans from the Brazilian Atlantic Forest, comprising the conspicuous, brightly colored pumpkin-toadlets and the cryptic flea-toads. Pumpkin-toadlets are known to contain tetrodotoxins and therefore, their bright colors may perform an aposematic function. Previous studies based on a limited number of mitochondrial and nuclear-encoded markers supported the existence of two clades containing species of pumpkin-toadlet phenotype, but deep nodes remained largely unresolved or conflicting between data sets. We use new RNAseq data of 17 individuals from nine Brachycephalus species to infer their evolutionary relationships from a phylogenomic perspective. Analyses of almost 5300 nuclear-encoded ortholog protein-coding genes and full mitochondrial genomes confirmed the existence of two separate pumpkin-toadlet clades, suggesting the convergent evolution (or multiple reversals) of the bufoniform morphology, conspicuous coloration, and probably toxicity. In addition, the study of the mitochondrial gene order revealed that three species (B. hermogenesi, B. pitanga, and B. rotenbergae) display translocations of different tRNAs (NCY and CYA) from the WANCY tRNA cluster to a position between the genes ATP6 and COIII, showing a new mitochondrial gene order arrangement for vertebrates. The newly clarified phylogeny suggests that Brachycephalus has the potential to become a promising model taxon to understand the evolution of coloration, body plan and toxicity. Given that toxicity information is available for only few species of Brachycephalus, without data for any flea-toad species, we also emphasize the need for a wider screening of toxicity across species, together with more in-depth functional and ecological study of their phenotypes.

A comprehensive phylogeny of dwarf geckos of the genus Lygodactylus, with insights into their systematics and morphological variation

The 71 currently known species of dwarf geckos of the genus Lygodactylus are a clade of biogeographic interest due to their occurrence in continental Africa, Madagascar, and South America. Furthermore, because many species are morphologically cryptic, our knowledge of species-level diversity within this genus is incomplete, as indicated by numerous unnamed genetic lineages revealed in previous molecular studies. Here we provide an extensive multigene phylogeny covering 56 of the named Lygodactylus species, four named subspecies, and 34 candidate species of which 19 are newly identified in this study. Phylogenetic analyses, based on ∼10.1 kbp concatenated sequences of eight nuclear-encoded and five mitochondrial gene fragments, confirm the monophyly of 14 Lygodactylus species groups, arranged in four major clades. We recover two clades splitting from basal nodes, one comprising exclusively Malagasy species groups, and the other containing three clades. In the latter, there is a clade with only Madagascar species, which is followed by a clade containing three African and one South American species groups, and its sister clade containing six African and two Malagasy species groups. Relationships among species groups within these latter clades remain weakly supported. We reconstruct a Lygodactylus timetree based on a novel fossil-dated phylotranscriptomic tree of squamates, in which we included data from two newly sequenced Lygodactylus transcriptomes. We estimate the crown diversification of Lygodactylus started at 46 mya, and the dispersal of Lygodactylus among the main landmasses in the Oligocene and Miocene, 35-22 mya, but emphasize the wide confidence intervals of these estimates. The phylogeny suggests an initial out-of-Madagascar dispersal as most parsimonious, but accounting for poorly resolved nodes, an out-of-Africa scenario may only require one extra dispersal step. More accurate inferences into the biogeographic history of these geckos will likely require broader sampling of related genera and phylogenomic approaches to provide better topological support. A survey of morphological characters revealed that most of the major clades and species groups within Lygodactylus cannot be unambiguously characterized by external morphology alone, neither by unique character states nor by a diagnostic combination of character states. Thus, any future taxonomic work will likely benefit from integrative, phylogenomic approaches.

The phylogeny of the Casque-headed Treefrogs (Hylidae: Hylinae: Lophyohylini)

The South American and West Indian Casque-headed Treefrogs (Hylidae: Hylinae: Lophyohylini) include 85 species. These are notably diverse in morphology (e.g. disparate levels of cranial hyperossification) and life history (e.g. different reproductive modes, chemical defences), have a wide distribution, and occupy habitats from the tropical rainforests to semiarid scrubland. In this paper, we present a phylogenetic analysis of this hylid tribe based on sequence fragments of up to five mitochondrial (12S, 16S, ND1, COI, Cytb) and six nuclear genes (POMC, RAG-1, RHOD, SIAH, TNS3, TYR). We included most of its species (> 96%), in addition to a number of new species. Our results indicate: (i) the paraphyly of Trachycephalus with respect to Aparasphenodon venezolanus; (ii) the nonmonophyly of Aparasphenodon, with Argenteohyla siemersi, Corythomantis galeata and Nyctimantis rugiceps nested within it, and Ap. venezolanus nested within Trachycephalus; (iii) the polyphyly of Corythomantis; (iv) the nonmonophyly of the recognized species groups of Phyllodytes; and (v) a pervasive low support for the deep relationships among the major clades of Lophyohylini, including C. greeningi and the monotypic genera Itapotihyla and Phytotriades. To remedy the nonmonophyly of Aparasphenodon, Corythomantis, and Trachycephalus, we redefined Nyctimantis to include Aparasphenodon (with the exception of Ap. venezolanus, which we transferred to Trachycephalus), Argenteohyla, and C. galeata. Additionally, our results indicate the need for taxonomic work in the following clades: (i) Trachycephalus dibernardoi and Tr. imitatrix; (ii) Tr. atlas, Tr. mambaiensis and Tr. nigromaculatus; and (iii) Phyllodytes. On the basis of our phylogenetic results, we analyzed the evolution of skull hyperossification and reproductive biology, with emphasis on the multiple independent origins of phytotelm breeding, in the context of Anura. We also analyzed the inter-related aspects of chemical defences, venom delivery, phragmotic behaviour, co-ossification, and prevention of evaporative water loss.

The phylogeny of Dendropsophini (Anura: Hylidae: Hylinae)

The relationships of the hyline tribe Dendropsophini remain poorly studied, with most published analyses dealing with few of the species groups of Dendropsophus. In order to test the monophyly of Dendropsophini, its genera, and the species groups currently recognized in Dendropsophus, we performed a total evidence phylogenetic analysis. The molecular dataset included sequences of three mitochondrial and five nuclear genes from 210 terminals, including 12 outgroup species, the two species of Xenohyla, and 93 of the 108 recognized species of Dendropsophus. The phenomic dataset includes 46 terminals, one per species (34 Dendropsophus, one Xenohyla, and 11 outgroup species). Our results corroborate the monophyly of Dendropsophini and the reciprocal monophyly of Dendropsophus and Xenohyla. Some species groups of Dendropsophus are paraphyletic (the D. microcephalus, D. minimus, and D. parviceps groups, and the D. rubicundulus clade). On the basis of our results, we recognize nine species groups; for three of them (D. leucophyllatus, D. microcephalus, and D. parviceps groups) we recognize some nominal clades to highlight specific morphology or relationships and facilitate species taxonomy. We further discuss the evolution of oviposition site selection, where our results show multiple instances of independent evolution of terrestrial egg clutches during the evolutionary history of Dendropsophus.

The first mitochondrial genome of a South America parthenogenetic lizard (Squamata: Gymnophthalmidae)

The mitogenome of the South American parthenogenetic lizard Loxopholis percarinatum Müller, 1923 (Squamata: Gymnophthalmidae), a uni-bisexual species complex, was recovered for three individuals from Rio Negro region, Amazonas, Brazil. The content and order of genes are typical for vertebrate mitochondrial genomes, and we recovered 13 protein-coding genes, 22 tRNA, and two rRNA (12S and 16S), in addition to partial fragments of the Control Region. A maximum likelihood phylogenetic analysis with mitogenomes of selected lizard families recovered L. percarinatum with Iphisa elegans Gray, 1851, the only other Gymnophthalmidae species available in GenBank.

Geography, Host Genetics, and Cross-Domain Microbial Networks Structure the Skin Microbiota of Fragmented Brazilian Atlantic Forest Frog Populations

The host-associated microbiome plays a significant role in health. However, the roles of factors such as host genetics and microbial interactions in determining microbiome diversity remain unclear. We examined these factors using amplicon-based sequencing of 175 Thoropa taophora frog skin swabs collected from a naturally fragmented landscape in southeastern Brazil. Specifically, we examined (1) the effects of geography and host genetics on microbiome diversity and structure; (2) the structure of microbial eukaryotic and bacterial co-occurrence networks; and (3) co-occurrence between microeukaryotes with bacterial OTUs known to affect growth of the fungal pathogen Batrachochytrium dendrobatidis (Bd). While bacterial alpha diversity varied by both site type and host MHC IIB genotype, microeukaryotic alpha diversity varied only by site type. However, bacteria and microeukaryote composition showed variation according to both site type and host MHC IIB genotype. Our network analysis showed the highest connectivity when both eukaryotes and bacteria were included, implying that ecological interactions may occur among domains. Lastly, anti-Bd bacteria were not broadly negatively co-associated with the fungal microbiome and were positively associated with potential amphibian parasites. Our findings emphasize the importance of considering both domains in microbiome research and suggest that for effective probiotic strategies for amphibian disease management, considering potential interactions among all members of the microbiome is crucial.

Asymmetric acoustic wave scattering by a nonreciprocal and position-dependent mass defect

We investigate the asymmetric wave scattering in a phononic one-dimensional lattice with a nonreciprocal defect and position dependent masses coupled by the defect spring. The nonreciprocal interaction is characterized by a single parameter Δ while the nonlinear contribution due to position-dependent masses are controlled by a parameterχ. The transmission and reflection coefficients are analytically computed and the effects of the nonreciprocity and nonlinearity are detailed. We show that, in opposite with the linear case, the rectification factor has a frequency dependence, which leads to a more efficient diode-like action at large wavevectors. Further, the nonlinearity leads to an asymmetry of the reflected component, absent in the linear regime. We extend our analysis to a system with frictional forces which suppresses the multistability window promoted by the nonlinear mass contribution without compromising the rectification action.

Successful application of ancient DNA extraction and library construction protocols to museum wet collection specimens

Millions of scientific specimens are housed in museum collections, a large part of which are fluid preserved. The use of formaldehyde as fixative and subsequent storage in ethanol is especially common in ichthyology and herpetology. This type of preservation damages DNA and reduces the chance of successful retrieval of genetic data. We applied ancient DNA extraction and single stranded library construction protocols to a variety of vertebrate samples obtained from wet collections and of different ages. Our results show that almost all samples tested yielded endogenous DNA. Archival DNA extraction was successful across different tissue types as well as using small amounts of tissue. Conversion of archival DNA fragments into single-stranded libraries resulted in usable data even for samples with initially undetectable DNA amounts. Subsequent target capture approaches for mitochondrial DNA using homemade baits on a subset of 30 samples resulted in almost complete mitochondrial genome sequences in several instances. Thus, application of ancient DNA methodology makes wet collection specimens, including type material as well as rare, old or extinct species, accessible for genetic and genomic analyses. Our results, accompanied by detailed step-by-step protocols, are a large step forward to open the DNA archive of museum wet collections for scientific studies.

Low-temperature pseudo-phase-transition in an extended Hubbard diamond chain

We consider the extended Hubbard diamond chain with an arbitrary number of particles driven by chemical potential. The interaction between dimer diamond chain and nodal couplings is considered in the atomic limit (no hopping), whereas the dimer interaction includes the hopping term. We demonstrate that this model exhibits a pseudo-transition effect in the low-temperature regime. Here, we explore the pseudo-transition rigorously by analyzing several physical quantities. The internal energy and entropy depict sudden, although continuous, jumps which closely resembles discontinuous or first-order phase-transition. At the same time, the correlation length and specific heat exhibit astonishing strong sharp peaks quite similar to a second-order phase-transition. We associate the ascending and descending parts of the peak with power-law "pseudo-critical" exponents. We determine the pseudo-critical exponents in the temperature range where these peaks are developed, namely, ν=1 for the correlation length and α=3 for the specific heat. We also study the behavior of the electron density and isothermal compressibility around the pseudo-critical temperature.

Mitogenome analyses elucidate the evolutionary relationships of a probable Eocene wet tropics relic in the xerophile lizard genus Acanthodactylus

Climate has a large impact on diversity and evolution of the world's biota. The Eocene-Oligocene transition from tropical climate to cooler, drier environments was accompanied by global species turnover. A large number of Old World lacertid lizard lineages have diversified after the Eocene-Oligocene boundary. One of the most speciose reptile genera in the arid Palearctic, Acanthodactylus, contains two sub-Saharan species with unresolved phylogenetic relationship and unknown climatic preferences. We here aim to understand how and when adaptation to arid conditions occurred in Acanthodactylus and when tropical habitats where entered. Using whole mitogenomes from fresh and archival DNA and published sequences we recovered a well-supported Acanthodactylus phylogeny and underpinned the timing of diversification with environmental niche analyses of the sub-Saharan species A. guineensis and A. boueti in comparison to all arid Acanthodactylus. We found that A. guineensis represents an old lineage that splits from a basal node in the Western clade, and A. boueti is a derived lineage and probably not its sister. Their long branches characterize them-and especially A. guineensis-as lineages that may have persisted for a long time without further diversification or have undergone multiple extinctions. Environmental niche models verified the occurrence of A. guineensis and A. boueti in hot humid environments different from the other 42 arid Acanthodactylus species. While A. guineensis probably remained in tropical habitat from periods prior to the Eocene-Oligocene boundary, A. boueti entered tropical environments independently at a later period. Our results provide an important baseline for studying adaptation and the transition from humid to arid environments in Lacertidae.

Skin microbiome correlates with bioclimate and Batrachochytrium dendrobatidis infection intensity in Brazil's Atlantic Forest treefrogs

In Brazil’s Atlantic Forest (AF) biodiversity conservation is of key importance since the fungal pathogen Batrachochytrium dendrobatidis (Bd) has led to the rapid loss of amphibian populations here and worldwide. The impact of Bd on amphibians is determined by the host's immune system, of which the skin microbiome is a critical component. The richness and diversity of such cutaneous bacterial communities are known to be shaped by abiotic factors which thus may indirectly modulate host susceptibility to Bd. This study aimed to contribute to understanding the environment-host–pathogen interaction determining skin bacterial communities in 819 treefrogs (Anura: Hylidae and Phyllomedusidae) from 71 species sampled across the AF. We investigated whether abiotic factors influence the bacterial community richness and structure on the amphibian skin. We further tested for an association between skin bacterial community structure and Bd co-occurrence. Our data revealed that temperature, precipitation, and elevation consistently correlate with richness and diversity of the skin microbiome and also predict Bd infection status. Surprisingly, our data suggest a weak but significant positive correlation of Bd infection intensity and bacterial richness. We highlight the prospect of future experimental studies on the impact of changing environmental conditions associated with global change on environment-host–pathogen interactions in the AF.

Critical properties of the SIS model on the clustered homophilic network

The spreading of epidemics in complex networks has been a subject of renewed interest of several scientific branches. In this regard, we have focused our attention on the study of the susceptible-infected-susceptible (SIS) model, within a Monte Carlo numerical simulation approach, representing the spreading of epidemics in a clustered homophilic network. The competition between infection and recovery that drives the system either to an absorbing or to an active phase is analyzed. We estimate the static critical exponents β ∕ ν , 1 ∕ ν and γ ∕ ν , through finite-size scaling (FSS) analysis of the order parameter ρ and its fluctuations, showing that they differ from those associated with the contact process on a scale-free network, as well as those predicted by the heterogeneous mean-field theory.

Multiple origins of green coloration in frogs mediated by a novel biliverdin-binding serpin

Many vertebrates have distinctive blue-green bones and other tissues due to unusually high biliverdin concentrations-a phenomenon called chlorosis. Despite its prevalence, the biochemical basis, biology, and evolution of chlorosis are poorly understood. In this study, we show that the occurrence of high biliverdin in anurans (frogs and toads) has evolved multiple times during their evolutionary history, and relies on the same mechanism-the presence of a class of serpin family proteins that bind biliverdin. Using a diverse combination of techniques, we purified these serpins from several species of nonmodel treefrogs and developed a pipeline that allowed us to assemble their complete amino acid and nucleotide sequences. The described proteins, hereafter named biliverdin-binding serpins (BBS), have absorption spectra that mimic those of phytochromes and bacteriophytochromes. Our models showed that physiological concentration of BBSs fine-tune the color of the animals, providing the physiological basis for crypsis in green foliage even under near-infrared light. Additionally, we found that these BBSs are most similar to human glycoprotein alpha-1-antitrypsin, but with a remarkable functional diversification. Our results present molecular and functional evidence of recurrent evolution of chlorosis, describe a biliverdin-binding protein in vertebrates, and introduce a function for a member of the serpin superfamily, the largest and most ubiquitous group of protease inhibitors.

Emerging extreme value and Fermi-Dirac distributions in the Lévy branching and annihilating process

We study the dynamics of the branching and annihilating process with long-range interactions. Static particles generate an offspring and annihilate upon contact. The branching distance is supposed to follow a Lévy-like power-law distribution with P(r)∝1/r^{α}. We analyze the long term behavior of the mean particles number and its fluctuations as a function of the parameter α that controls the range of the branching process. We show that the dynamic exponent associated with the particle number fluctuations varies continuously for α<4 while the particle number exponent only changes for α<3. A crossover from extreme value Frechet (at α=3) and Gumbell (for 2<α<3) distributions is developed, similar to the one reported in recent experiments with cw-pumped random fiber lasers presenting underlying gain and Lévy processes. We report the dependence of the relevant dynamical power-law exponents on α showing that explosive growth takes place for α≤2. Further, the average occupation number distribution is shown to evolve from the standard Fermi-Dirac form to the generalized one within the context of nonextensive statistics.

Symbiotic bacterial communities in rainforest fungus-farming ants: evidence for species and colony specificity

Animals may host diverse bacterial communities that can markedly affect their behavioral physiology, ecology, and vulnerability to disease. Fungus-farming ants represent a classical example of mutualism that depends on symbiotic microorganisms. Unraveling the bacterial communities associated with fungus-farming ants is essential to understand the role of these microorganisms in the ant-fungus symbiosis. The bacterial community structure of five species of fungus-farmers (non-leaf-cutters; genera Mycocepurus, Mycetarotes, Mycetophylax, and Sericomyrmex) from three different environments in the Brazilian Atlantic rainforest (lowland forest, restinga forest, and sand dunes) was characterized with amplicon-based Illumina sequencing of 16 S ribosomal RNA gene. Possible differences in bacterial communities between ants internal to the nest (on the fungus garden) and external foragers were also investigated. Our results on the richness and diversity of associated bacteria provide novel evidence that these communities are host- and colony-specific in fungus-farming ants. Indeed, the bacterial communities associated with external foragers differ among the five species, and among colonies of the same species. Furthermore, bacterial communities from internal ants vs. foragers do not differ or differ only slightly within each ant species. This study highlights the importance of describing ant-associated bacterial communities to better understand this host-bacterial interaction in the social environment of insect colonies and provides the foundation for future studies on the ecological and evolutionary processes that drive the success of fungus-farming ants.

High-throughput DNA sequencing of museum specimens sheds light on the long-missing species of the Bokermannohyla claresignata group (Anura: Hylidae: Cophomantini)

The two species of the Bokermannohyla claresignata species group (Anura: Hylidae) have not been collected for the last four decades. It is the only species group of the hyline tribe Cophomantini that has not yet been analysed genetically. Its phylogenetic position is thus uncertain, and it has a combination of adult and larval character states that make this group a crucial missing piece that hinders our understanding of Cophomantini phylogenetics and character evolution. We obtained DNA sequences from a museum larval specimen of Bok. claresignata, using specialized extraction methods and high-throughput DNA sequencing, and combined the molecular phylogenetic results with available phenotypic information to provide new insights into the taxonomy and phylogenetic relationships of its species group. Our phylogenetic results place Bok. claresignata as sister to the Boana pulchella group, supporting its inclusion in Boana, together with Bokermannohyla clepsydra. In light of this new finding, we recognize a newly defined Boana claresignata group to accommodate these species, thus resolving both the polyphyly of Bokermannohyla and the paraphyly of Boana. Considering the phylogenetic relationships of the Boana claresignata group, we also discuss the evolution of suctorial tadpoles and mature oocyte/egg pigmentation in Cophomantini.

Rectification of acoustic phonons in harmonic chains with nonreciprocal spring defects

The scattering of acoustic phonons by nonreciprocal spring defects inserted in an harmonic chain is investigated. The degree of nonreciprocity of the forces mediated by the defect springs is parameterized by a single quantity Δ that effectively takes into account the interaction of the coupled masses with hidden degrees of freedom of an underlying nonequilibrium system. We demonstrate a pronounced rectification effect with transmission having a preferential direction. Nonreciprocity also allows energy exchange between the system and the medium. Further, we show a cooperative action between defects mediated by resonant cavity modes. The influence of damping forces is also explored and shown to promote the rectification of the reflected vibrational wave component.

Anurans of the Caparaó National Park and surroundings, southeast Brazil

Abstract The Atlantic Forest (AF) is one of the biodiversity hotspots of the world, and the most fragmented biome of Brazil. This biome includes different phytophysiognomies, as riparian, slope, cloudy forests, and grasslands. Such complexity, allied to huge latitudinal and high elevational range, provides diverse habitats and conditions for amphibian speciation. As a result, there are over 600 amphibian species known to occur in the AF. Within this biome the Caparaó National Park (CNP) is relevant, as it includes the highest peak of the biome, the Pico da Bandeira at almost 3,000 m above sea level, as well as different phytophysiognomies as rocky fields and humid forests. In spite of that, its amphibian fauna is still poorly described. We inventoried amphibians at the CNP and surrounding areas from 2016 to 2018 and recorded 47 anuran species, of which two are locally endemic and at least six have not been described yet. Additionally, we compiled data from previous surveys (2004 to 2008) and secondary data from scientific collections. All together, we registered a total of 61 anuran species from 12 families for the CNP and surroundings, placing this area among the 10 amphibian richest sites in the AF. Some of these species are represented by only one or two collected specimens and have not been registered in the CNP since the 1980’s, such as Thoropa lutzi and Hylodes vanzolinii. These species could be examples of population declines or even past local extinctions, highlighting the need of further sampling efforts in that highly biodiverse site.

A review of the elusive bicolored iris Snouted Treefrogs (Anura: Hylidae:Scinax uruguayus group)

The genus Scinax currently includes more than 120 species, recovered in two major clades, the S. catharinae and the S. ruber clades. The latter comprises 75 species, most of which remain unassigned to any species groups, while 12 are included in the S. rostratus and S. uruguayus groups. In this paper we present a taxonomic review of the two species currently included in the S. uruguayus group, discussing some putative phenotypic synapomorphies of this group. Although S. pinima and S. uruguayus have been considered as distinct species, this has been based on scant evidence, and several authors doubted of their distinctiveness. Our study of available specimens of S. pinima and S. uruguayus corroborates that both are valid and diagnosable species based on phenotypic evidence. Furthermore, our results show that S. pinima previously known only from its type locality, has a much widespread distribution than previously thought (including the Brazilian states of Paraná, Santa Catarina, and Rio Grande do Sul), which, added to the biological information presented here allows to suggest the removal of this species from the "Data Deficient" IUCN Red List category to "Least Concern". Also, we describe a new species formerly reported as S. aff. pinima and S. uruguayus from NE Argentina and some localities from the Brazilian State of Rio Grande do Sul. All species are diagnosed and characterized using adult and larval morphology, osteology, vocalizations, cytogenetics, and natural history.

Arthropod-bacteria interactions influence assembly of aquatic host microbiome and pathogen defense

The host-associated microbiome is vital to host immunity and pathogen defense. In aquatic ecosystems, organisms may interact with environmental bacteria to influence the pool of potential symbionts, but the effects of these interactions on host microbiome assembly and pathogen resistance are unresolved. We used replicated bromeliad microecosystems to test for indirect effects of arthropod-bacteria interactions on host microbiome assembly and pathogen burden, using tadpoles and the fungal amphibian pathogen Batrachochytrium dendrobatidis as a model host-pathogen system. Arthropods influenced host microbiome assembly by altering the pool of environmental bacteria, with arthropod-bacteria interactions specifically reducing host colonization by transient bacteria and promoting antimicrobial components of aquatic bacterial communities. Arthropods also reduced fungal zoospores in the environment, but fungal infection burdens in tadpoles corresponded most closely with arthropod-mediated patterns in microbiome assembly. This result indicates that the cascading effects of arthropods on the maintenance of a protective host microbiome may be more strongly linked to host health than negative effects of arthropods on pools of pathogenic zoospores. Our work reveals tight links between healthy ecosystem dynamics and the functioning of host microbiomes, suggesting that ecosystem disturbances such as loss of arthropods may have downstream effects on host-associated microbial pathogen defenses and host fitness.

Localization-delocalization transition in discrete-time quantum walks with long-range correlated disorder

We study the effects of spatially long-range correlated phase disorder on the Hadamard quantum walk on a line. The shift operator is built to exhibit an intrinsic disorder distribution featuring long-range correlations. To impose such, we resort to fractional Brownian motion with power-law spectrum 1/k^{2α} with α≥0 being the exponent that controls the degree of correlations. We discuss the scaling behavior of the walker's wave packet and report a localization-delocalization transition controlled by α. We unveil two intermediate dynamical regimes between exponential localization and full delocalization.

High Genetic Diversity and No Population Structure of the New World Screwworm Fly Cochliomyia hominivorax (Diptera: Calliphoridae) on a Microgeographic Scale: Implications for Management Units

The New World screwworm fly Cochliomyia hominivorax (Coquerel, 1858) (Diptera: Calliphoridae) is an important livestock pest endemic to the Americas that has been eradicated from North and continental Central America with a control program based on the Sterile Insect Technique (SIT). The establishment of target management units is a strategic step in the implementation of new control programs, which can be achieved using genetic studies of natural populations. Previous studies of New World screwworm fly populations were conducted on the continental scale and identified four main groups: two in South America and two in the Caribbean. However, studies within these groups are needed to determine which smaller geographic areas can be treated as management units. Here, we analyze the genetic variability distribution and the population demographic signals of the New World screwworm fly in a 6,000 km2 area located along the border of Brazil and Uruguay. This area was the subject of the first control pilot program conducted in South America. We studied eight microsatellite loci and sequences from two mitochondrial DNA regions in individuals sampled at 20-25 livestock breeding farms. We observed no population structure and found high genetic variability on the geographical scale sampled for both molecular markers. Our microsatellite data suggest that these populations are not in equilibrium, and demographic analyses based on mitochondrial data indicate population expansion. These results suggest that this geographic scale is not adequate for future New World screwworm fly management in South America.

The Intestinal Microbiota of Tadpoles Differs from Those of Syntopic Aquatic Invertebrates

Bacterial communities associated to eukaryotes play important roles in the physiology, development, and health of their hosts. Here, we examine the intestinal microbiota in tadpoles and aquatic invertebrates (insects and gastropods) to better understand the degree of specialization in the tadpole microbiotas. Samples were collected at the same time in one pond, and the V4 region of the bacterial 16S rRNA gene was sequenced with Illumina amplicon sequencing. We found that bacterial richness and diversity were highest in two studied snail individuals, intermediate in tadpoles, and lowest in the four groups of aquatic insects. All groups had substantial numbers of exclusive bacterial operational taxonomic units (OTUs) in their guts, but also shared a high proportion of OTUs, probably corresponding to transient environmental bacteria. Significant differences were found for all pairwise comparisons of tadpoles and snails with the major groups of insects, but not among insect groups or between snails and tadpoles. The similarity between tadpoles and snails may be related to similar feeding mode as both snails and tadpoles scratch biofilms and algae from surfaces; however, this requires confirmation due to low sample sizes. Overall, the gut microbiota differences found among syntopic aquatic animals are likely shaped by both food preferences and host identity.

Phase diagram and re-entrant fermionic entanglement in a hybrid Ising-Hubbard ladder

The degree of fermionic entanglement is examined in an exactly solvable Ising-Hubbard ladder, which involves interacting electrons on the ladder's rungs described by Hubbard dimers at half-filling on each rung, accounting for intrarung hopping and Coulomb terms. The coupling between neighboring Hubbard dimers is assumed to have an Ising-like nature. The ground-state phase diagram consists of four distinct regions corresponding to the saturated paramagnetic, the classical antiferromagnetic, the quantum antiferromagnetic, and the mixed classical-quantum phase. We have exactly computed the fermionic concurrence, which measures the degree of quantum entanglement between the pair of electrons on the ladder rungs. The effects of the hopping amplitude, the Coulomb term, temperature, and magnetic fields on the fermionic entanglement are explored in detail. It is shown that the fermionic concurrence displays a re-entrant behavior when quantum entanglement is being generated at moderate temperatures above the classical saturated paramagnetic ground state.

Ground-state phase diagram, fermionic entanglement and kinetically-induced frustration in a hybrid ladder with localized spins and mobile electrons

We introduce an exactly solvable hybrid spin-ladder model containing localized nodal Ising spins and interstitial mobile electrons, which are allowed to perform a quantum-mechanical hopping between the ladder's legs. The quantum-mechanical hopping process induces an antiferromagnetic coupling between the ladder's legs that competes with a direct exchange coupling of the nodal spins. The model is exactly mapped onto the Ising spin ladder with temperature-dependent two- and four-spin interactions, which is subsequently solved using the transfer-matrix technique. We report the ground-state phase diagram and compute the fermionic concurrence to characterize the quantum entanglement between the pair of interstitial mobile electrons. We further provide a detailed analysis of the local spin ordering including the pair and four-spin correlation functions around an elementary plaquette, as well as, the local ordering diagrams. It is shown that a complex sequence of distinct local orderings and frustrated correlations takes place when the model parameters drive the investigated system close to a zero-temperature triple coexistence point.

Bose-Einstein condensation in chains with power-law hoppings: Exact mapping on the critical behavior in d-dimensional regular lattices

Recent experimental progress on the realization of quantum systems with highly controllable long-range interactions has impelled the study of quantum phase transitions in low-dimensional systems with power-law couplings. Long-range couplings mimic higher-dimensional effects in several physical contexts. Here, we provide the exact relation between the spectral dimension d at the band bottom and the exponent α that tunes the range of power-law hoppings of a one-dimensional ideal lattice Bose gas. We also develop a finite-size scaling analysis to obtain some relevant critical exponents and the critical temperature of the BEC transition. In particular, an irrelevant dangerous scaling field has to be taken into account when the hopping range is sufficiently large to make the effective dimensionality d>4.

The mitochondrial genomes of three species of poison frogs (Anura: Dendrobates)

We reconstructed nearly complete mitogenomes for three species of poison frogs, Dendrobates auratus, D. leucomelas, and D. tinctorius, from RNAseq data. We recovered the 13 protein-coding genes, 22 tRNA genes (except tRNA-Val for D. leucomelas), and two rRNA genes for all three species, plus partial sequences of the control region. The order of genes agrees with that known from a previously sequenced D. auratus, being the most commonly found for neobatrachian frogs. Based on full-sibling comparisons we estimate the probable error rate of Illumina-RNAseq reconstructed mitogenomes of up to 0.15%.

Meeting the challenge of DNA barcoding Neotropical amphibians: polymerase chain reaction optimization and new COI primers

Amphibians are one of the most threatened vertebrate classes, yet at the same time new species are being described every year, demonstrating that the number of existing species is grossly underestimated. In groups such as amphibians, with high extinction rates and poorly known species boundaries, DNA barcoding is a tool that can rapidly assess genetic diversity and estimate species richness for prioritizing conservation decisions. However, reliable recovery of the 5' region of the cytochrome c oxidase subunit 1 (COI) gene is critical for the ongoing effort to gather DNA barcodes for all amphibian species. Here, we provide new PCR conditions and tested new primers that increase the efficiency of barcode recovery in amphibians. We found that a low extension temperature for PCR cycles significantly improves the efficiency of amplification for all combinations of primers. Combining low PCR extension temperature and primers AnF1 + AnR1, we were able to recover COI sequences for 100% of the species analysed (N = 161), encompassing ~15% of the species known from Brazil (representing 77 genera and 23 families), which is an important improvement over previous studies. The preliminary assessment of species diversity suggested that number of species might be underestimated by about 25%. We conclude that DNA barcoding is an efficient, simple, and standardized protocol for identifying cryptic diversity in amphibians and advocate for its use in biodiversity inventories and across widespread populations within known species.