The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation.

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomics resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomics resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti- predator strategies, and resilience and adaptive responses. They also serve as critical models for understanding widespread genomic characteristics, including evolutionary genome expansions and contractions given they have the largest range in genome sizes of any animal taxon and multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The advent of long-read sequencing technologies, along with computational techniques that enhance scaffolding capabilities and streamline computational workload is now enabling the ability to overcome some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC) in early 2023. This burgeoning community already has more than 282 members from 41 countries (6 in Africa, 131 in the Americas, 27 in Asia, 29 in Australasia, and 89 in Europe). The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and outline how the AGC can enable amphibian genomics research to "leap" to the next level.

Effects of overwintering on the transcriptome and fitness traits in a damselfly with variable voltinism across two latitudes.

Winter diapause consists of cessation of development that allows individuals to survive unfavourable conditions. Winter diapause may bear various costs and questions have been raised about the evolutionary mechanisms maintaining facultative diapause. Here, we explored to what extent a facultative winter diapause affects life-history traits and the transcriptome in the damselfly Ischnura elegans, and whether these effects were latitude-specific. We collected adult females at central and high latitudes and raised their larvae in growth chambers. Larvae were split into a non-diapausing and post-winter (diapausing) cohort, were phenotyped and collected for a gene expression analysis. At the phenotypic level, we found no difference in survival between the two cohorts, and the post-winter cohort was larger and heavier than the non-winter cohort. These effects were mostly independent of the latitude of origin. At the transcriptomic level, wintering affected gene expression with a small fraction of genes significantly overlapping across latitudes, especially those related to morphogenesis. In conclusion, we found clear effects of diapause on the phenotype but little evidence for latitudinal-specific effects of diapause. Our results showed a shared transcriptomic basis underpinning diapause demonstrated, here, at the intraspecific level and supported the idea of evolutionary convergence of the response to diapause across organisms.

Lizard host abundances and climatic factors explain phylogenetic diversity and prevalence of blood parasites on an oceanic island.

Host abundance might favour the maintenance of a high phylogenetic diversity of some parasites via rapid transmission rates. Blood parasites of insular lizards represent a good model to test this hypothesis because these parasites can be particularly prevalent in islands and host lizards highly abundant. We applied deep amplicon sequencing and analysed environmental predictors of blood parasite prevalence and phylogenetic diversity in the endemic lizard Gallotia galloti across 24 localities on Tenerife, an island in the Canary archipelago that has experienced increasing warming and drought in recent years. Parasite prevalence assessed by microscopy was over 94%, and a higher proportion of infected lizards was found in warmer and drier locations. A total of 33 different 18s rRNA parasite haplotypes were identified, and the phylogenetic analyses indicated that they belong to two genera of Adeleorina (Apicomplexa: Coccidia), with Karyolysus as the dominant genus. The most important predictor of between-locality variation in parasite phylogenetic diversity was the abundance of lizard hosts. We conclude that a combination of climatic and host demographic factors associated with an insular syndrome may be favouring a rapid transmission of blood parasites among lizards on Tenerife, which may favour the maintenance of a high phylogenetic diversity of parasites.

Genomic analysis reveals complex population structure within the smooth newt, Lissotriton vulgaris, in Central Europe.

Species with wide-range distributions usually display high genetic variation. This variation can be partly explained by historical lineages that were temporally isolated from each other and are back into secondary reproductive contact, and partly by local adaptations. The smooth newt (Lissotriton vulgaris) is one of the most widely distributed amphibians species across Eurasia and forms a species complex with a partially overlapping distribution and morphology. In the present study, we explored the population genomic structure of smooth newt lineages in the Carpathian Basin (CB) relying on single-nucleotide polymorphisms. Our dataset included new and previously published data to study the secondary contact zone between lineages in the CB and also tested for the barrier effect of rivers to gene flow between these lineages. We confirmed the presence of the South L. v. vulgaris Lineage distributed in Transdanubia and we provided new distribution records of L. v. ampelensis inhabiting the eastern territories of the CB. High genetic diversity of smooth newts was observed, especially in the North Hungarian Mountains and at the interfluves of the main rivers in the South with four distinct lineages of L. v. vulgaris and one lineage of L. v. ampelensis showing a low level of admixture with the spatially closest L. v. vulgaris lineage. Moreover, admixture detected at the interfluve of the main rivers (i.e. Danube and Tisza) suggested a secondary contact zone in the area. Finally, we found that the river Danube has a very weak effect on population divergence, while the river Tisza is a geographical barrier limiting gene flow between smooth newt lineages. As the range boundaries of L. v. vulgaris and L. v. ampelensis in the CB coincide with the river Tisza, our study underpins the influence of rivers in lineage diversification.

Latitude-specific urbanization effects on life history traits in the damselfly Ischnura elegans.

Many species are currently adapting to cities at different latitudes. Adaptation to urbanization may require eco-evolutionary changes in response to temperature and invasive species that may differ between latitudes. Here, we studied single and combined effects of increased temperatures and an invasive alien predator on the phenotypic response of replicated urban and rural populations of the damselfly Ischnura elegans and contrasted these between central and high latitudes. Adult females were collected in rural and urban ponds at central and high latitudes. Their larvae were exposed to temperature treatments (current [20°C], mild warming [24°C], and heat wave [28°C; for high latitude only]) crossed with the presence or absence of chemical cues released by the spiny-cheek crayfish (Faxonius limosus), only present at the central latitude. We measured treatment effects on larval development time, mass, and growth rate. Urbanization type affected all life history traits, yet these responses were often dependent on latitude, temperature, and sex. Mild warming decreased mass in rural and increased growth rate in urban populations. The effects of urbanization type on mass were latitude-dependent, with central-latitude populations having a greater phenotypic difference. Urbanization type effects were sex-specific with urban males being lighter and having a lower growth rate than rural males. At the current temperature and mild warming, the predator cue reduced the growth rate, and this independently of urbanization type and latitude of origin. This pattern was reversed during a heat wave in high-latitude damselflies. Our results highlight the context-dependency of evolutionary and plastic responses to urbanization, and caution for generalizing how populations respond to cities based on populations at a single latitude.

Understanding the evolution of immune genes in jawed vertebrates.

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

The effect of temperature and invasive alien predator on genetic and phenotypic variation in the damselfly Ischnura elegans: cross-latitude comparison.

BACKGROUND: Understanding and predicting how organisms respond to human-caused environmental changes has become a major concern in conservation biology. Here, we linked gene expression and phenotypic data to identify candidate genes underlying existing phenotypic trait differentiation under individual and combined environmental variables. For this purpose, we used the damselfly Ischnura elegans. Egg clutches from replicated high- (southern Sweden) and central-latitude (southern Poland) populations facing different degrees of seasonal time constraints were collected. Damselfly larvae were exposed to experimental treatments: current and mild warming temperatures crossed with the presence or absence of an invasive alien predator cue released by the spiny-cheek crayfish, Faxonius limosus, which is only present in Poland to date. We measured the following traits: larval development time, body size, mass and growth rate, and used the larvae for gene expression analysis by RNA-seq. Data were analysed using a multivariate approach. RESULTS: We showed latitudinal differences in coping with mild warming and predator cues. When exposed to an increased temperature and a predator cue, central-latitude individuals had the shortest development and the fastest growth compared to high-latitude individuals. There was a general effect of predator cues regarding mass and growth rate reduction independent of latitude. Transcriptome analysis revealed that metabolic pathways related to larval anatomy and development tended to be upregulated in response to mild warming but only in fast-growing central-latitude individuals. Metabolic pathways linked to oxidative stress tended to be downregulated in response to a predator cue, especially in central-latitude individuals. CONCLUSION: Different phenotypic and transcriptomic responses to environmental factors might be attributed to the variability in I. elegans life history strategies between the two latitudes caused by seasonal time constraints and to its coexistence with the invasive alien predator in nature. By providing insights into how organisms may respond to future anthropogenic changes, our results may be of particular interest in conservation biology.

Salamanders reveal novel trajectories of amphibian MHC evolution.

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I) in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution.

Emerging infectious diseases of amphibians in Poland: distribution and environmental drivers.

Emerging infectious diseases are a threat to biodiversity and have taken a large toll on amphibian populations worldwide. The chytrid fungi Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal), and the iridovirus Ranavirus (Rv), are of concern as all have contributed to amphibian declines. In central and eastern Europe, their geographical and host distributions and main environmental drivers determining prevalence are poorly known. We screened over 1000 amphibians from natural and captive populations in Poland for the presence of Bd, Bsal and Rv. In wild amphibian populations, we found that Bd is widespread, present in 46 out of 115 sampled localities as well as 2 captive colonies, and relatively common with overall prevalence at 14.4% in 9 species. We found lower prevalence of Rv at 2.4%, present in 11 out of 92 sampling sites, with a taxonomic breadth of 8 different amphibian species. Bsal infection was not detected in any individuals. In natural populations, Pelophylax esculentus and Bombina variegata accounted for 75% of all Bd infections, suggesting a major role for these 2 species as pathogen reservoirs in Central European freshwater habitats. General linear models showed that climatic as well as landscape features are associated with Bd infection in Poland. We found that higher average annual temperature constrains Bd infection, while landscapes with numerous water bodies or artificial elements (a surrogate for urbanization) increase the chances of infection. Our results show that a combination of climatic and landscape variables may drive regional and local pathogen emergence.

Coevolution between MHC Class I and Antigen-Processing Genes in Salamanders.

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.

Molecular Evolution of Antigen-Processing Genes in Salamanders: Do They Coevolve with MHC Class I Genes?

Proteins encoded by antigen-processing genes (APGs) prepare antigens for presentation by the major histocompatibility complex class I (MHC I) molecules. Coevolution between APGs and MHC I genes has been proposed as the ancestral gnathostome condition. The hypothesis predicts a single highly expressed MHC I gene and tight linkage between APGs and MHC I. In addition, APGs should evolve under positive selection, a consequence of the adaptive evolution in MHC I. The presence of multiple highly expressed MHC I genes in some teleosts, birds, and urodeles appears incompatible with the coevolution hypothesis. Here, we use urodele amphibians to test two key expectations derived from the coevolution hypothesis: 1) the linkage between APGs and MHC I was studied in Lissotriton newts and 2) the evidence for adaptive evolution in APGs was assessed using 42 urodele species comprising 21 genera from seven families. We demonstrated that five APGs (PSMB8, PSMB9, TAP1, TAP2, and TAPBP) are tightly linked (<0.5 cM) to MHC I. Although all APGs showed some codons under episodic positive selection, we did not find a pervasive signal of positive selection expected under the coevolution hypothesis. Gene duplications, putative gene losses, and divergent allelic lineages detected in some APGs demonstrate considerable evolutionary dynamics of APGs in salamanders. Overall, our results indicate that if coevolution between APGs and MHC I occurred in urodeles, it would be more complex than envisaged in the original formulation of the hypothesis.

Evolutionary principles guiding amphibian conservation.

The Anthropocene has witnessed catastrophic amphibian declines across the globe. A multitude of new, primarily human-induced drivers of decline may lead to extinction, but can also push species onto novel evolutionary trajectories. If these are recognized by amphibian biologists, they can be engaged in conservation actions. Here, we summarize how principles stemming from evolutionary concepts have been applied for conservation purposes, and address emerging ideas at the vanguard of amphibian conservation science. In particular, we examine the consequences of increased drift and inbreeding in small populations and their implications for practical conservation. We then review studies of connectivity between populations at the landscape level, which have emphasized the limiting influence of anthropogenic structures and degraded habitat on genetic cohesion. The rapid pace of environmental changes leads to the central question of whether amphibian populations can cope either by adapting to new conditions or by shifting their ranges. We gloomily conclude that extinction seems far more likely than adaptation or range shifts for most species. That said, conservation strategies employing evolutionary principles, such as selective breeding, introduction of adaptive variants through translocations, ecosystem interventions aimed at decreasing phenotype-environment mismatch, or genetic engineering, may effectively counter amphibian decline in some areas or for some species. The spread of invasive species and infectious diseases has often had disastrous consequences, but has also provided some premier examples of rapid evolution with conservation implications. Much can be done in terms of setting aside valuable amphibian habitat that should encompass both natural and agricultural areas, as well as designing protected areas to maximize the phylogenetic and functional diversity of the amphibian community. We conclude that an explicit consideration and application of evolutionary principles, although certainly not a silver bullet, should increase effectiveness of amphibian conservation in both the short and long term.

Differential introgression across newt hybrid zones: Evidence from replicated transects.

Genomic heterogeneity of divergence between hybridizing species may reflect heterogeneity of introgression, but also processes unrelated to hybridization. Heterogeneous introgression and its repeatability can be directly tested in natural hybrid zones by examining multiple transects. Here, we studied hybrid zones between the European newts Lissotriton montandoni and two lineages of Lissotriton vulgaris, with replicate transects within each zone. Over 1,000 nuclear genes located on a linkage map and mitochondrial DNA were investigated using geographical and genomic clines. Overall, the five transects were all similar, showing hallmarks of strong reproductive isolation: bimodal distribution of genotypes in central populations and narrow allele frequency clines. However, the extent of introgression differed between the zones, possibly as a consequence of their different ages, as suggested by the analysis of heterozygosity runs in diagnostic markers. In three transects genomic signatures of small-scale (~2 km) zone movements were detected. We found limited overlap of cline outliers between transects, and only weak evidence of stronger differentiation of introgression between zones than between transects within zones. Introgression was heterogeneous across linkage groups, with patterns of heterogeneity similar between transects and zones. Predefined candidates for increased or reduced introgression exhibited only a subtle tendency in the expected direction, suggesting that interspecific differentiation is not a reliable indicator for the strength of introgression. These hierarchically sampled hybrid zones of apparently different ages show how introgression unfolds with time and offer an excellent opportunity to dissect the dynamics of hybridization and architecture of reproductive isolation at advanced stages of speciation.

Mapping of quantitative trait loci for life history traits segregating within common frog populations.

The evolution of complex traits is often shaped by adaptive divergence. However, very little is known about the number, effect size, and location of the genomic regions influencing the variation of these traits in natural populations. Based on a dense linkage map of the common frog, Rana temporaria, we have localized, for the first time in amphibians, three significant and nine suggestive quantitative trait loci (QTLs) for metabolic rate, growth rate, development time, and weight at metamorphosis, explaining 5.6-18.9% of the overall phenotypic variation in each trait. We also found a potential pleiotropic QTL between development time and size at metamorphosis that, if confirmed, might underlie the previously reported genetic correlation between these traits. Furthermore, we demonstrate that the genetic variation linked to fitness-related larval traits segregates within Rana temporaria populations. This study provides the first insight into the genomic regions that affect larval life history traits in anurans, providing a valuable resource to delve further into the genomic basis of evolutionary change in amphibians.

Comparative High-Density Linkage Mapping Reveals Conserved Genome Structure but Variation in Levels of Heterochiasmy and Location of Recombination Cold Spots in the Common Frog.

By combining 7077 SNPs and 61 microsatellites, we present the first linkage map for some of the early diverged lineages of the common frog, Rana temporaria, and the densest linkage map to date for this species. We found high homology with the published linkage maps of the Eastern and Western lineages but with differences in the order of some markers. Homology was also strong with the genome of the Tibetan frog Nanorana parkeri and we found high synteny with the clawed frog Xenopus tropicalis We confirmed marked heterochiasmy between sexes and detected nonrecombining regions in several groups of the male linkage map. Contrary to the expectations set by the male heterogamety of the common frog, we did not find male heterozygosity excess in the chromosome previously shown to be linked to sex determination. Finally, we found blocks of loci showing strong transmission ratio distortion. These distorted genomic regions might be related to genetic incompatibilities between the parental populations, and are promising candidates for further investigation into the genetic basis of speciation and adaptation in the common frog.

Impact of an intervention on the prescription of aliskiren after new evidence on safety reported.

PURPOSE: The purpose of the study is to analyze the impact of an intervention to disseminate safety alerts on the utilization of Aliskiren added to angiotensin converting enzyme inhibitor (ACEI) or an angiotensin-receptor blocker (ARB). METHODS: Quasi-experimental design (non-randomized intervention) comparing the utilization of Aliskiren + ACEI or ARB in a primary care area-intervention (PCA-I) with a primary care area-control (PCA-C) following a safety alert. All physicians were provided with a list of diabetic patients (DM) on Aliskiren + ACEI or ARB. Physicians in the PCA-I received also a non-DM patients list, a report with recommendations and information on the utilization of Aliskiren + ACEI or ARB in their area. Information was obtained from electronic medical records, period from May 2010 to December 2012. Interrupted time series analysis were used to assess the effect of the intervention on the number of patients on Aliskiren + ACEI or ARB. RESULTS: The number of DM receiving Aliskiren + ACEI or ARB at the time of the alert (23 December 2011) was 106 in the PCA-I (91 non-DM) and 45 in the PCA-C (25 non-DM). After the alert, a decreased in the number of patients on Aliskiren + ACEI or ARB was noted at both PCAs, although the average of daily treatments ended was significantly higher in the PCA-I, both in the DM group (slope after alert: -0.81, 95%CI -0.91 to -0.71 vs. -0.30, 95%CI -0.37 to -0.22) as well as in the non-DM group (-0.56, 95%CI -0.67 to -0.45 vs. -0.10 95%CI -0.17 to -0.04). CONCLUSIONS: The prescription of Aliskiren + ACEI or ARB decreased at both PCAs, albeit such decreased was more significant at the PCA-I. The intervention led to a more expeditious implementation of the safety alert recommendations. Copyright © 2016 John Wiley & Sons, Ltd.

Heritability of Batrachochytrium dendrobatidis burden and its genetic correlation with development time in a population of Common toad (Bufo spinosus).

Despite the important threat that emerging pathogens pose for the conservation of biodiversity as well as human health, very little is known about the adaptive potential of host species to withstand infections. We studied the quantitative genetic architecture responsible for the burden of the fungal pathogen Batrachochytrium dendrobatidis in a population of common toads in conjunction with other life-history traits (i.e., body size and development rate) that may be affected by common selective pressures. We found a significant heritable component that is associated with fungal burden, which may allow for local adaptation to this pathogen to proceed. In addition, the high genetic correlation found between fungal burden and development time suggests that both traits have to be taken into account in order to assess the adaptive response of host populations to this emerging pathogen.