Larval fitness, microsatellite diversity and MHC class II diversity in common frog (Rana temporaria) populations

City divided: Unveiling family ties and genetic structuring of coyotes in Seattle

Linear barriers pose significant challenges for wildlife gene flow, impacting species persistence, adaptation, and evolution. While numerous studies have examined the effects of linear barriers (e.g., fences and roadways) on partitioning urban and non-urban areas, understanding their influence on gene flow within cities remains limited. Here, we investigated the impact of linear barriers on coyote (Canis latrans) population structure in Seattle, Washington, where major barriers (i.e., interstate highways and bodies of water) divide the city into distinct quadrants. Just under 1000 scats were collected to obtain genetic data between January 2021 and December 2022, allowing us to identify 73 individual coyotes. Notably, private allele analysis underscored limited interbreeding among quadrants. When comparing one quadrant to each other, there were up to 16 private alleles within a single quadrant, representing nearly 22% of the population allelic diversity. Our analysis revealed weak isolation by distance, and despite being a highly mobile species, genetic structuring was apparent between quadrants even with extremely short geographic distance between individual coyotes, implying that Interstate 5 and the Ship Canal act as major barriers. This study uses coyotes as a model species for understanding urban gene flow and its consequences in cities, a crucial component for bolstering conservation of rarer species and developing wildlife friendly cities.

Importance of Genetic-Fitness Correlations for the Conservation of Amphibians

Endangered animals suffer from isolation of their habitats. Isolation leads to a reduction in population size as well as a decrease in genetic diversity and a concomitant increase in the risk of extinction. Amphibians are the most endangered vertebrate class. Besides habitat loss, fragmentation and isolation, amphibians are threatened by emerging diseases e.g., chytrid fungus or Ranavirus. By employing experiments, researchers investigate whether changes in genetic diversity within or among isolated populations affect amphibian fitness. While genetic diversity estimates are based on molecular markers, typically microsatellites, fitness is mostly measured as tadpole performance in rearing experiments often under varying environmental conditions. Tadpole performances (e.g., body mass, growth rate and survival) have been found to be negatively affected by low genetic diversity, as several studies have found a positive association between genetic diversity and these fitness traits. Moreover, infection with pathogens also seems to be more likely in individuals or populations with lower genetic diversity. Overall, these genetic-fitness correlations seem to be more pronounced or detectable in smaller, declining populations but not in larger populations. Genomic studies, which sample a larger fraction of the genome, are still scarce in the conservation genetic literature on amphibians. These are likely to increase in upcoming years and may reveal adaptive variants that protect against dangerous pathogens or environmental changes. Altogether, genetic-fitness correlation studies should be a priority in order to develop effective management plans for the genetic rescue of isolated, imperilled amphibian populations.

The effects of urban life on animal immunity: Adaptations and constraints

Land transformation, including urbanization, is a dominant form of anthropogenic change to the global environment at the dawn of the Anthropocene epoch. More and more species are brought into direct contact with humans, being either required to develop broad-scale adaptations to urban environment or filtered out from urbanized areas. While behavioural or physiological adaptations are at the forefront of urban biology research, there is accumulating evidence for divergent pathogen pressure across urbanization gradients, requiring adjustments in host immune function. At the same time, host immunity may be constrained by unfavourable components of an urban environment, such as poor-quality food resources, disturbance, or pollution. Here, I reviewed existing evidence for adaptations and constrains in the immune system of urban animals, focusing on the recent implementation of metabarcoding, genomic, transcriptomic, and epigenomic approaches in urban biology research. I show that spatial variation in pathogen pressure across urban and non-urban landscapes is highly complex and may be context-dependent, but there is solid evidence for pathogen-driven immunostimulation in urban-dwelling animals. I also show that genes coding for molecules directly involved in interactions with pathogens are the prime candidates for immunogenetic adaptations to urban life. Evidence emerging from landscape genomics and transcriptomics show that immune adaptations to urban life may have a polygenic nature, but immune traits may not be among the key biological functions experiencing broad-scale microevolutionary changes in response to urbanization. Finally, I provided recommendations for future research, including i) a better integration of different 'omic' approaches to obtain a more complete picture of immune adaptations to urban life in non-model animal taxa, ii) quantification of fitness landscapes for immune phenotypes and genotypes across urbanization gradient, and iii) much broader taxonomic coverage (including invertebrates) necessary to draw more robust conclusions on how general (or taxa-specific) are immune responses of animals to urbanization.

Evolution of life in urban environments

Our planet is an increasingly urbanized landscape, with over half of the human population residing in cities. Despite advances in urban ecology, we do not adequately understand how urbanization affects the evolution of organisms, nor how this evolution may affect ecosystems and human health. Here, we review evidence for the effects of urbanization on the evolution of microbes, plants, and animals that inhabit cities. Urbanization affects adaptive and nonadaptive evolutionary processes that shape the genetic diversity within and between populations. Rapid adaptation has facilitated the success of some native species in urban areas, but it has also allowed human pests and disease to spread more rapidly. The nascent field of urban evolution brings together efforts to understand evolution in response to environmental change while developing new hypotheses concerning adaptation to urban infrastructure and human socioeconomic activity. The next generation of research on urban evolution will provide critical insight into the importance of evolution for sustainable interactions between humans and our city environments.

Limited polymorphism of the functional MHC class II B gene in the black-spotted frog (Pelophylax nigromaculatus) identified by locus-specific genotyping

Amphibians can be more vulnerable to environmental changes and diseases than other species because of their complex life cycle and physiological requirements. Therefore, understanding the adaptation of amphibians to environmental changes is crucial for their conservation. Major histocompatibility complex (MHC) presents an excellent tool for the investigation of adaptive variations and the assessment of adaptive potential because it can be under strong diversifying selection. Here, we isolated the MHC class II B (MHCIIB) gene from cDNA sequences of the black-spotted frog (Pelophylax nigromaculatus), a widespread amphibian species in China, and designed locus-specific primers to characterize adaptive variability of this amphibian. Ten alleles were identified from 67 individual frogs of three populations and no more than two alleles were present in each individual animal. Furthermore, none of the sequences had indels or/and stop codons, which is in good agreement with locus-specific amplification of a functional gene. However, we found low polymorphism at both nucleotide and amino acid levels, even in the antigen-binding region. Purifying selection acting at this locus was supported by the findings that the d_N/d_S ratio across all alleles was below 1 and that negatively selected sites were detected by different methods. Allele frequency distributions were significantly different among geographic populations, indicating that physiographic factors may have strong effect on the genetic structure of the black-spotted frog. This study revealed limited polymorphism of three adjacent black-spotted frog populations at the functional MHCIIB locus, which may be attributed to region-specific differences. The locus-specific genotyping technique developed in this study would provide a foundation for future studies on adaptive divergence among different frog populations.

Anthropogenic and ecological drivers of amphibian disease (ranavirosis)

Ranaviruses are causing mass amphibian die-offs in North America, Europe and Asia, and have been implicated in the decline of common frog (Rana temporaria) populations in the UK. Despite this, we have very little understanding of the environmental drivers of disease occurrence and prevalence. Using a long term (1992-2000) dataset of public reports of amphibian mortalities, we assess a set of potential predictors of the occurrence and prevalence of Ranavirus-consistent common frog mortality events in Britain. We reveal the influence of biotic and abiotic drivers of this disease, with many of these abiotic characteristics being anthropogenic. Whilst controlling for the geographic distribution of mortality events, disease prevalence increases with increasing frog population density, presence of fish and wild newts, increasing pond depth and the use of garden chemicals. The presence of an alternative host reduces prevalence, potentially indicating a dilution effect. Ranavirosis occurrence is associated with the presence of toads, an urban setting and the use of fish care products, providing insight into the causes of emergence of disease. Links between occurrence, prevalence, pond characteristics and garden management practices provides useful management implications for reducing the impacts of Ranavirus in the wild.

Drift rather than selection dominates MHC class II allelic diversity patterns at the biogeographical range scale in natterjack toads Bufo calamita

Study of major histocompatibility complex (MHC) loci has gained great popularity in recent years, partly due to their function in protecting vertebrates from infections. This is of particular interest in amphibians on account of major threats many species face from emergent diseases such as chytridiomycosis. In this study we compare levels of diversity in an expressed MHC class II locus with neutral genetic diversity at microsatellite loci in natterjack toad (Bufo (Epidalea) calamita) populations across the whole of the species' biogeographical range. Variation at both classes of loci was high in the glacial refugium areas (REF) and much lower in postglacial expansion areas (PGE), especially in range edge populations. Although there was clear evidence that the MHC locus was influenced by positive selection in the past, congruence with the neutral markers suggested that historical demographic events were the main force shaping MHC variation in the PGE area. Both neutral and adaptive genetic variation declined with distance from glacial refugia. Nevertheless, there were also some indications from differential isolation by distance and allele abundance patterns that weak effects of selection have been superimposed on the main drift effect in the PGE zone.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

Comparison of 454 pyrosequencing methods for characterizing the major histocompatibility complex of nonmodel species and the advantages of ultra deep coverage

Characterization and population genetic analysis of multilocus genes, such as those found in the major histocompatibility complex (MHC) is challenging in nonmodel vertebrates. The traditional method of extensive cloning and Sanger sequencing is costly and time-intensive and indirect methods of assessment often underestimate total variation. Here, we explored the suitability of 454 pyrosequencing for characterizing multilocus genes for use in population genetic studies. We compared two sample tagging protocols and two bioinformatic procedures for 454 sequencing through characterization of a 185-bp fragment of MHC DRB exon 2 in wolverines (Gulo gulo) and further compared the results with those from cloning and Sanger sequencing. We found 10 putative DRB alleles in the 88 individuals screened with between two and four alleles per individual, suggesting amplification of a duplicated DRB gene. In addition to the putative alleles, all individuals possessed an easily identifiable pseudogene. In our system, sequence variants with a frequency below 6% in an individual sample were usually artefacts. However, we found that sample preparation and data processing procedures can greatly affect variant frequencies in addition to the complexity of the multilocus system. Therefore, we recommend determining a per-amplicon-variant frequency threshold for each unique system. The extremely deep coverage obtained in our study (approximately 5000×) coupled with the semi-quantitative nature of pyrosequencing enabled us to assign all putative alleles to the two DRB loci, which is generally not possible using traditional methods. Our method of obtaining locus-specific MHC genotypes will enhance population genetic analyses and studies on disease susceptibility in nonmodel wildlife species.