Advancing Understanding of Amphibian Evolution, Ecology, Behavior, and Conservation with Massively Parallel Sequencing

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

Genetic insight into a polygenic trait using a novel genome-wide association approach in a wild amphibian population

Body size variation is central in the evolution of life-history traits in amphibians, but the underlying genetic architecture of this complex trait is still largely unknown. Herein, we studied the genetic basis of body size and fecundity of the alternative morphotypes in a wild population of the Greek smooth newt (Lissotriton graecus). By combining a genome-wide association approach with linkage disequilibrium network analysis, we were able to identify clusters of highly correlated loci thus maximizing sequence data for downstream analysis. The putatively associated variants explained 12.8% to 44.5% of the total phenotypic variation in body size and were mapped to genes with functional roles in the regulation of gene expression and cell cycle processes. Our study is the first to provide insights into the genetic basis of complex traits in newts and provides a useful tool to identify loci potentially involved in fitness-related traits in small data sets from natural populations in non-model species.

Genomic footprints of (pre) colonialism: Population declines in urban and forest túngara frogs coincident with historical human activity

Urbanisation is rapidly altering ecosystems, leading to profound biodiversity loss. To mitigate these effects, we need a better understanding of how urbanisation impacts dispersal and reproduction. Two contrasting population demographic models have been proposed that predict that urbanisation either promotes (facilitation model) or constrains (fragmentation model) gene flow and genetic diversity. Which of these models prevails likely depends on the strength of selection on specific phenotypic traits that influence dispersal, survival, or reproduction. Here, we a priori examined the genomic impact of urbanisation on the Neotropical túngara frog (Engystomops pustulosus), a species known to adapt its reproductive traits to urban selective pressures. Using whole-genome resequencing for multiple urban and forest populations we examined genomic diversity, population connectivity and demographic history. Contrary to both the fragmentation and facilitation models, urban populations did not exhibit substantial changes in genomic diversity or differentiation compared with forest populations, and genomic variation was best explained by geographic distance rather than environmental factors. Adopting an a posteriori approach, we additionally found both urban and forest populations to have undergone population declines. The timing of these declines appears to coincide with extensive human activity around the Panama Canal during the last few centuries rather than recent urbanisation. Our study highlights the long-lasting legacy of past anthropogenic disturbances in the genome and the importance of considering the historical context in urban evolution studies as anthropogenic effects may be extensive and impact nonurban areas on both recent and older timescales.

Chromosome-level assembly, annotation and phylome of Pelobates cultripes, the western spadefoot toad

Genomic resources for amphibians are still hugely under-represented in vertebrate genomic research, despite being a group of major interest for ecology, evolution and conservation. Amphibians constitute a highly threatened group of vertebrates, present a vast diversity in reproductive modes, are extremely diverse in morphology, occupy most ecoregions of the world, and present the widest range in genome sizes of any major group of vertebrates. We combined Illumina, Nanopore and Hi-C sequencing technologies to assemble a chromosome-level genome sequence for an anuran with a moderate genome size (assembly span 3.09 Gb); Pelobates cultripes, the western spadefoot toad. The genome has an N50 length of 330 Mb with 98.6% of the total sequence length assembled into 14 super scaffolds, and 87.7% complete BUSCO genes. We use published transcriptomic data to provide annotations, identifying 32,684 protein-coding genes. We also reconstruct the P. cultripes phylome and identify 2,527 gene expansions. We contribute the first draft of the genome of the western spadefoot toad, P. cultripes. This species represents a relatively basal lineage in the anuran tree with an interesting ecology and a high degree of developmental plasticity, and thus is an important resource for amphibian genomic research.

Molecular mechanisms of local adaptation for salt-tolerance in a treefrog

Salinization is a global phenomenon affecting ecosystems and forcing freshwater organisms to deal with increasing levels of ionic stress. However, our understanding of mechanisms that permit salt tolerance in amphibians is limited. This study investigates mechanisms of salt tolerance in locally adapted, coastal populations of a treefrog, Hyla cinerea. Using a common garden experiment, we (i) determine the extent that environment (i.e., embryonic and larval saltwater exposure) or genotype (i.e., coastal vs. inland) affects developmental benchmarks and transcriptome expression, and (ii) identify genes that may underpin differences in saltwater tolerance. Differences in gene expression, survival, and plasma osmolality were most strongly associated with genotype. Population genetic analyses on expressed genes also delineated coastal and inland groups based on genetic similarity. Coastal populations differentially expressed osmoregulatory genes including ion transporters (atp1b1, atp6V1g2, slc26a), cellular adhesion components (cdh26, cldn1, gjb3, ocln), and cytoskeletal components (odc1-a, tgm3). Several of these genes are the same genes expressed by euryhaline fish after exposure to freshwater, which is a novel finding for North American amphibians and suggests that these genes may be associated with local salinity adaptation. Coastal populations also highly expressed glycerol-3-phosphate dehydrogenase 1 (gpd1), which indicates they use glycerol as a compatible osmolyte to reduce water loss - another mechanism of saltwater tolerance previously unknown in frogs. These data signify that Hyla cinerea inhabiting coastal, brackish wetlands have evolved a salt-tolerant ecotype, and highlights novel candidate pathways that can lead to salt tolerance in freshwater organisms facing habitat salinization.

Analysis of muntjac deer genome and chromatin architecture reveals rapid karyotype evolution

Closely related muntjac deer show striking karyotype differences. Here we describe chromosome-scale genome assemblies for Chinese and Indian muntjacs, Muntiacus reevesi (2n = 46) and Muntiacus muntjak vaginalis (2n = 6/7), and analyze their evolution and architecture. The genomes show extensive collinearity with each other and with other deer and cattle. We identified numerous fusion events unique to and shared by muntjacs relative to the cervid ancestor, confirming many cytogenetic observations with genome sequence. One of these M. muntjak fusions reversed an earlier fission in the cervid lineage. Comparative Hi-C analysis showed that the chromosome fusions on the M. muntjak lineage altered long-range, three-dimensional chromosome organization relative to M. reevesi in interphase nuclei including A/B compartment structure. This reshaping of multi-megabase contacts occurred without notable change in local chromatin compaction, even near fusion sites. A few genes involved in chromosome maintenance show evidence for rapid evolution, possibly associated with the dramatic changes in karyotype.