Antimicrobial peptide and sequence variation along a latitudinal gradient in two anurans

BACKGROUND

While there is evidence of both purifying and balancing selection in immune defense genes, large-scale genetic diversity in antimicrobial peptides (AMPs), an important part of the innate immune system released from dermal glands in the skin, has remained uninvestigated. Here we describe genetic diversity at three AMP loci (Temporin, Brevinin and Palustrin) in two ranid frogs (Rana arvalis and R. temporaria) along a 2000 km latitudinal gradient. We amplified and sequenced part of the Acidic Propiece domain and the hypervariable Mature Peptide domain (~ 150-200 bp) in the three genes using Illumina Miseq and expected to find decreased AMP genetic variation towards the northern distribution limit of the species similarly to studies on MHC genetic patterns.

RESULTS

We found multiple loci for each AMP and relatively high gene diversity, but no clear pattern of geographic genetic structure along the latitudinal gradient. We found evidence of trans-specific polymorphism in the two species, indicating a common evolutionary origin of the alleles. Temporin and Brevinin did not form monophyletic clades suggesting that they belong to the same gene family. By implementing codon evolution models we found evidence of strong positive selection acting on the Mature Peptide. We also found evidence of diversifying selection as indicated by divergent allele frequencies among populations and high Theta k values.

CONCLUSION

Our results suggest that AMPs are an important source of adaptive diversity, minimizing the chance of microorganisms developing resistance to individual peptides.

Genome Assembly of the Common Pheasant Phasianus colchicus: A Model for Speciation and Ecological Genomics

The common pheasant (Phasianus colchicus) in the order Galliformes and the family Phasianidae, has 30 subspecies distributed across its native range in the Palearctic realm and has been introduced to Europe, North America, and Australia. It is an important game bird often subjected to wildlife management as well as a model species to study speciation, biogeography, and local adaptation. However, the genomic resources for the common pheasant are generally lacking. We sequenced a male individual of the subspecies torquatus of the common pheasant with the Illumina HiSeq platform. We obtained 94.88 Gb of usable sequences by filtering out low-quality reads of the raw data generated. This resulted in a 1.02 Gb final assembly, which equals the estimated genome size. BUSCO analysis using chicken as a model showed that 93.3% of genes were complete. The contig N50 and scaffold N50 sizes were 178 kb and 10.2 Mb, respectively. All these indicate that we obtained a high-quality genome assembly. We annotated 16,485 protein-coding genes and 123.3 Mb (12.05% of the genome) of repetitive sequences by ab initio and homology-based prediction. Furthermore, we applied a RAD-sequencing approach for another 45 individuals of seven representative subspecies in China and identified 4,376,351 novel single nucleotide polymorphism (SNPs) markers. Using this unprecedented data set, we uncovered the geographic population structure and genetic introgression among common pheasants in China. Our results provide the first high-quality reference genome for the common pheasant and a valuable genome-wide SNP database for studying population genomics and demographic history.

Endless forms of sexual selection

In recent years, the field of sexual selection has exploded, with advances in theoretical and empirical research complementing each other in exciting ways. This perspective piece is the product of a “stock-taking” workshop on sexual selection and sexual conflict. Our aim is to identify and deliberate on outstanding questions and to stimulate discussion rather than provide a comprehensive overview of the entire field. These questions are organized into four thematic sections we deem essential to the field. First we focus on the evolution of mate choice and mating systems. Variation in mate quality can generate both competition and choice in the opposite sex, with implications for the evolution of mating systems. Limitations on mate choice may dictate the importance of direct vs. indirect benefits in mating decisions and consequently, mating systems, especially with regard to polyandry. Second, we focus on how sender and receiver mechanisms shape signal design. Mediation of honest signal content likely depends on integration of temporally variable social and physiological costs that are challenging to measure. We view the neuroethology of sensory and cognitive receiver biases as the main key to signal form and the ‘aesthetic sense’ proposed by Darwin. Since a receiver bias is sufficient to both initiate and drive ornament or armament exaggeration, without a genetically correlated or even coevolving receiver, this may be the appropriate ‘null model’ of sexual selection. Thirdly, we focus on the genetic architecture of sexually selected traits. Despite advances in modern molecular techniques, the number and identity of genes underlying performance, display and secondary sexual traits remains largely unknown. In-depth investigations into the genetic basis of sexual dimorphism in the context of long-term field studies will reveal constraints and trajectories of sexually selected trait evolution. Finally, we focus on sexual selection and conflict as drivers of speciation. Population divergence and speciation are often influenced by an interplay between sexual and natural selection. The extent to which sexual selection promotes or counteracts population divergence may vary depending on the genetic architecture of traits as well as the covariance between mating competition and local adaptation. Additionally, post-copulatory processes, such as selection against heterospecific sperm, may influence the importance of sexual selection in speciation. We propose that efforts to resolve these four themes can catalyze conceptual progress in the field of sexual selection, and we offer potential avenues of research to advance this progress.

Genotype by environment interaction for female fertility traits under conventional and organic production systems in Danish Holsteins

Conventional and organic production systems mainly differ in feeding strategies, outdoor and pasture access, and the use of antibiotic treatments. These environmental differences could lead to a genotype by environment interaction (G × E) and a requirement for including G × E in breeding decisions. The objectives of this study were to estimate variance components and heritabilities for conventional and organic production systems and investigate G × E under these 2 production systems for female fertility traits in Danish Holsteins. The analyzed traits included the interval from calving to first insemination (ICF), the interval from first to last insemination, number of inseminations per conception (NINS), and non-return rate within 56 d after the first insemination. Records of female fertility in heifers and the first 3 lactations in cows as well as grass ratio of feed at herd level were collected during the period from 2011 to 2016. The performances of a trait in heifers and cows (lactation 1 to 3) were considered as different traits. The (co)variance components and the resulting heritabilities and genetic correlations were estimated using 2 models. One was a bivariate model treating performances of a trait under organic and conventional production systems as 2 different traits using a reduced data set, and the other was a reaction norm model with random regression on the production system and the grass ratio of feed using a full data set. The full data set comprised records of 37,836 females from 112 organic herds and 513,599 females from 1,224 conventional herds, whereas the reduced data set comprised records from all these 112 organic herds and 92,696 females from 185 convention herds extracted from the full data set with grass ratio of feed lower than 0.20. All female fertility performances of the organic production system were superior to those of the conventional production system. Besides, heterogeneities in additive genetic variances and heritabilities were observed between conventional and organic production systems for all traits. Furthermore, genetic correlations between these 2 production systems ranged from 0.607 to 1.000 estimated from bivariate models and from 0.848 to 0.999 estimated from reaction norm models. Statistically significant G × E were observed for NINS in heifers, non-return rate within 56 d after the first insemination in heifers, and ICF from the bivariate model, and for ICF and NINS in cows from the reaction norm model.

Urban environment and reservoir host species are associated with Batrachochytrium dendrobatidis infection prevalence in the common toad

Human-induced changes of the environment, including landscape alteration and habitat loss, may affect wildlife disease dynamics and have important ramifications for wildlife conservation. Amphibians are among the vertebrate taxa most threatened by anthropogenic habitat change. The emerging fungal pathogen Batrachochytrium dendrobatidis (Bd) has caused extinctions and population declines in hundreds of anuran species globally. We studied how the urban landscape is associated with the prevalence of Bd infections by sampling 655 anurans of 3 species (mainly the common toad Bufo bufo) in 42 ponds surrounded by different amounts of urban habitat (defined as towns, cities or villages). We also examined the association between Bd infections and a potential reservoir host species (the moor frog Rana arvalis). We found that 38% of the sites were positive for Bd with an infection prevalence of 4.4%. The extent of urban landscape was negatively correlated with Bd infection prevalence. However, the positive association of Bd with the presence of the possible reservoir species was substantially stronger than the urban effects. The body condition index of B. bufo was negatively associated with Bd infection. This Bd effect was stronger than the negative effect of urban landscape on body condition. Our results suggest that urban environments in Sweden have a negative impact on Bd infections, while the presence of the reservoir species has a positive impact on Bd prevalence. Our study also highlights the potential importance of Bd infection on host fitness, especially in rural landscapes.

Genomic regions of speciation and adaptation among three species of grouse

Understanding the molecular basis of adaption is one of the central goals in evolutionary biology and when investigated across sister species it can provide detailed insight into the mechanisms of speciation. Here, we sequence the genomes of 34 individuals from three closely related grouse species in order to uncover the genomic architecture of speciation and the genes involved in adaptation. We identify 6 regions, containing 7 genes that show lineage specific signs of differential selection across the species. These genes are involved in a variety of cell processes ranging from stress response to neural, gut, olfactory and limb development. Genome wide neutrality test statistics reveal a strong signal of population expansion acting across the genomes. Additionally, we uncover a 3.5 Mb region on chromosome 20 that shows considerably lower levels of differentiation across the three grouse lineages, indicating possible action of uniform selection in this region.

Effects of host species and environmental factors on the prevalence of Batrachochytrium dendrobatidis in northern Europe

The fungal pathogen Batrachochytrium dendrobatidis (Bd) poses a major threat to amphibian populations. To assist efforts to address such threats, we examined differences in Bd host infection prevalence among amphibian species and its relations to both local environmental factors in breeding habitats and landscape variables measured at three scales (500, 2000 and 5000 m radii) around breeding sites in southernmost Sweden. We sampled 947 anurans of six species in 31 ponds and assessed their infection status. We then examined correlations of infection prevalence with canopy cover, pond perimeter and pH (treated as local-scale pond characteristics), and the number of ponds, area of arable land, area of mature forest, number of resident people and presence of sea within the three radii (treated as landscape variables). The Bd infection prevalence was very low, 0.5–1.0%, in two of the six anuran species (Bufo bufo and Rana temporaria), and substantially higher (13–64%) in the other four (Bombina bombina, Bufotes variabilis, Epidalea calamita, Rana arvalis). In the latter four species Bd infection prevalence was positively associated with ponds’ pH (site range: 5.3–8.1), and negatively associated with areas of mature forest and/or wetlands in the surroundings. Our results show that the infection dynamics of Bd are complex and associated with host species, local pond characteristics and several landscape variables at larger spatial scales. Knowledge of environmental factors associated with Bd infections and differences in species’ susceptibility may help to counter further spread of the disease and guide conservation action plans, especially for the most threatened species.

Past and potential future population dynamics of three grouse species using ecological and whole genome coalescent modeling

Studying demographic history of species provides insight into how the past has shaped the current levels of overall biodiversity and genetic composition of species, but also how these species may react to future perturbations. Here we investigated the demographic history of the willow grouse (Lagopus lagopus), rock ptarmigan (Lagopus muta), and black grouse (Tetrao tetrix) through the Late Pleistocene using two complementary methods and whole genome data. Species distribution modeling (SDM) allowed us to estimate the total range size during the Last Interglacial (LIG) and Last Glacial Maximum (LGM) as well as to indicate potential population subdivisions. Pairwise Sequentially Markovian Coalescent (PSMC) allowed us to assess fluctuations in effective population size across the same period. Additionally, we used SDM to forecast the effect of future climate change on the three species over the next 50 years. We found that SDM predicts the largest range size for the cold-adapted willow grouse and rock ptarmigan during the LGM. PSMC captured intraspecific population dynamics within the last glacial period, such that the willow grouse and rock ptarmigan showed multiple bottlenecks signifying recolonization events following the termination of the LGM. We also see signals of population subdivision during the last glacial period in the black grouse, but more data are needed to strengthen this hypothesis. All three species are likely to experience range contractions under future warming, with the strongest effect on willow grouse and rock ptarmigan due to their limited potential for northward expansion. Overall, by combining these two modeling approaches, we have provided a multifaceted examination of the biogeography of these species and how they have responded to climate change in the past. These results help us understand how cold-adapted species may respond to future climate changes.

Development and worldwide use of non-lethal, and minimal population-level impact, protocols for the isolation of amphibian chytrid fungi

Parasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to researchers as part of the BiodivERsA project RACE; here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been applied across 5 continents, 23 countries and in 62 amphibian species. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this group of emerging pathogens.

Recent Asian origin of chytrid fungi causing global amphibian declines

Globalized infectious diseases are causing species declines worldwide, but their source often remains elusive. We used whole-genome sequencing to solve the spatiotemporal origins of the most devastating panzootic to date, caused by the fungus Batrachochytrium dendrobatidis, a proximate driver of global amphibian declines. We traced the source of B. dendrobatidis to the Korean peninsula, where one lineage, BdASIA-1, exhibits the genetic hallmarks of an ancestral population that seeded the panzootic. We date the emergence of this pathogen to the early 20th century, coinciding with the global expansion of commercial trade in amphibians, and we show that intercontinental transmission is ongoing. Our findings point to East Asia as a geographic hotspot for B. dendrobatidis biodiversity and the original source of these lineages that now parasitize amphibians worldwide.

Signatures of historical selection on MHC reveal different selection patterns in the moor frog (Rana arvalis)

MHC genes are key components in disease resistance and an excellent system for studying selection acting on genetic variation in natural populations. Current patterns of variation in MHC genes are likely to be influenced by past and ongoing selection as well as demographic fluctuations in population size such as those imposed by post-glacial recolonization processes. Here, we investigated signatures of historical selection and demography on an MHC class II gene in 12 moor frog populations along a 1700-km latitudinal gradient. Sequences were obtained from 207 individuals and consecutively assigned into two different clusters (northern and southern clusters, respectively) in concordance with a previously described dual post-glacial colonization route. Selection analyses comparing the relative rates of non-synonymous to synonymous substitutions (dN/dS) suggested evidence of different selection patterns in the northern and the southern clusters, with divergent selection prevailing in the south but uniform positive selection predominating in the north. Also, models of codon evolution revealed considerable differences in the strength of selection: The southern cluster appeared to be under strong selection while the northern cluster showed moderate signs of selection. Our results indicate that the MHC alleles in the north diverged from southern MHC alleles as a result of differential selection patterns.

Mind the gaps in research on the control of gastrointestinal nematodes of farmed ruminants and pigs

Gastrointestinal (GI) nematode control has an important role to play in increasing livestock production from a limited natural resource base and to improve animal health and welfare. In this synthetic review, we identify key research priorities for GI nematode control in farmed ruminants and pigs, to support the development of roadmaps and strategic research agendas by governments, industry and policymakers. These priorities were derived from the DISCONTOOLS gap analysis for nematodes and follow-up discussions within the recently formed Livestock Helminth Research Alliance (LiHRA). In the face of ongoing spread of anthelmintic resistance (AR), we are increasingly faced with a failure of existing control methods against GI nematodes. Effective vaccines against GI nematodes are generally not available, and anthelmintic treatment will therefore remain a cornerstone for their effective control. At the same time, consumers and producers are increasingly concerned with environmental issues associated with chemical parasite control. To address current challenges in GI nematode control, it is crucial to deepen our insights into diverse aspects of epidemiology, AR, host immune mechanisms and the socio-psychological aspects of nematode control. This will enhance the development, and subsequent uptake, of the new diagnostics, vaccines, pharma-/nutraceuticals, control methods and decision support tools required to respond to the spread of AR and the shifting epidemiology of GI nematodes in response to climatic, land-use and farm husbandry changes. More emphasis needs to be placed on the upfront evaluation of the economic value of these innovations as well as the socio-psychological aspects to prioritize research and facilitate uptake of innovations in practice. Finally, targeted regulatory guidance is needed to create an innovation-supportive environment for industries and to accelerate the access to market of new control tools.

Phylogeography and larval spine length of the dragonfly Leucorhinia dubia in Europe

Presence or absence of predators selects for different kind of morphologies. Hence, we expect variation in traits that protect against predators to vary over geographical areas where predators vary in past and present abundance. Abdominal larval spines in dragonfly larvae provide protection against fish predators. We studied geographical variation in larval spine length of the dragonfly Leucorrhinia dubia across Western Europe using a phylogenetic approach. Larvae were raised in a common garden laboratory experiment in the absence of fish predators. Results show that larvae from northern Europe (Sweden and Finland) had significantly longer larval spines compared to larvae from western and central Europe. A phylogeny based on SNP data suggests that short larval spines is the ancestral stage in the localities sampled in this study, and that long spines have evolved in the Fenno-Scandian clade. The role of predators in shaping the morphological differences among the sampled localities is discussed.

Drift, selection, or migration? Processes affecting genetic differentiation and variation along a latitudinal gradient in an amphibian

BACKGROUND

Past events like fluctuations in population size and post-glacial colonization processes may influence the relative importance of genetic drift, migration and selection when determining the present day patterns of genetic variation. We disentangle how drift, selection and migration shape neutral and adaptive genetic variation in 12 moor frog populations along a 1700 km latitudinal gradient. We studied genetic differentiation and variation at a MHC exon II locus and a set of 18 microsatellites.

RESULTS

Using outlier analyses, we identified the MHC II exon 2 (corresponding to the β-2 domain) locus and one microsatellite locus (RCO8640) to be subject to diversifying selection, while five microsatellite loci showed signals of stabilizing selection among populations. STRUCTURE and DAPC analyses on the neutral microsatellites assigned populations to a northern and a southern cluster, reflecting two different post-glacial colonization routes found in previous studies. Genetic variation overall was lower in the northern cluster. The signature of selection on MHC exon II was weaker in the northern cluster, possibly as a consequence of smaller and more fragmented populations.

CONCLUSION

Our results show that historical demographic processes combined with selection and drift have led to a complex pattern of differentiation along the gradient where some loci are more divergent among populations than predicted from drift expectations due to diversifying selection, while other loci are more uniform among populations due to stabilizing selection. Importantly, both overall and MHC genetic variation are lower at northern latitudes. Due to lower evolutionary potential, the low genetic variation in northern populations may increase the risk of extinction when confronted with emerging pathogens and climate change.

Blood transcriptomes and de novo identification of candidate loci for mating success in lekking great snipe (Gallinago media)

We assembled the great snipe blood transcriptome using data from fourteen lekking males, in order to de novo identify candidate genes related to sexual selection, and determined the expression profiles in relation to mating success. The three most highly transcribed genes were encoding different haemoglobin subunits. All tended to be overexpressed in males with high mating success. We also called single nucleotide polymorphisms (SNPs) from the transcriptome data and found considerable genetic variation for many genes expressed during lekking. Among these, we identified 14 polymorphic candidate SNPs that had a significant genotypic association with mating success (number of females mated with) and/or mating status (mated or not). Four of the candidate SNPs were found in HBAA (encoding the haemoglobin α-chain). Heterozygotes for one of these and one SNP in the gene PABPC1 appeared to enjoy higher mating success compared to males homozygous for either of the alleles. In a larger data set of individuals, we genotyped 38 of the identified SNPs but found low support for consistent selection as only one of the zygosities of previously identified candidate SNPs and none of their genotypes were associated with mating status. However, candidate SNPs generally showed lower levels of spatial genetic structure compared to noncandidate markers. We also scored the prevalence of avian malaria in a subsample of birds. Males infected with avian malaria parasites had lower mating success in the year of sampling than noninfected males. Parasite infection and its interaction with specific genes may thus affect performance on the lek.

Range shifts or extinction? Ancient DNA and distribution modelling reveal past and future responses to climate warming in cold-adapted birds

Global warming is predicted to cause substantial habitat rearrangements, with the most severe effects expected to occur in high-latitude biomes. However, one major uncertainty is whether species will be able to shift their ranges to keep pace with climate-driven environmental changes. Many recent studies on mammals have shown that past range contractions have been associated with local extinctions rather than survival by habitat tracking. Here, we have used an interdisciplinary approach that combines ancient DNA techniques, coalescent simulations and species distribution modelling, to investigate how two common cold-adapted bird species, willow and rock ptarmigan (Lagopus lagopus and Lagopus muta), respond to long-term climate warming. Contrary to previous findings in mammals, we demonstrate a genetic continuity in Europe over the last 20 millennia. Results from back-casted species distribution models suggest that this continuity may have been facilitated by uninterrupted habitat availability and potentially also the greater dispersal ability of birds. However, our predictions show that in the near future, some isolated regions will have little suitable habitat left, implying a future decrease in local populations at a scale unprecedented since the last glacial maximum.

Looking into the past - the reaction of three grouse species to climate change over the last million years using whole genome sequences

Tracking past population fluctuations can give insight into current levels of genetic variation present within species. Analysing population dynamics over larger timescales can be aligned to known climatic changes to determine the response of species to varying environments. Here, we applied the Pairwise Sequentially Markovian Coalescent (psmc) model to infer past population dynamics of three widespread grouse species; black grouse, willow grouse and rock ptarmigan. This allowed the tracking of the effective population size (Ne ) of all three species beyond 1 Mya, revealing that (i) early Pleistocene cooling (~2.5 Mya) caused an increase in the willow grouse and rock ptarmigan populations, (ii) the mid-Brunhes event (~430 kya) and following climatic oscillations decreased the Ne of willow grouse and rock ptarmigan, but increased the Ne of black grouse and (iii) all three species reacted differently to the last glacial maximum (LGM) - black grouse increased prior to it, rock ptarmigan experienced a severe bottleneck and willow grouse was maintained at large population size. We postulate that the varying psmc signal throughout the LGM depicts only the local history of the species. Nevertheless, the large population fluctuations in willow grouse and rock ptarmigan indicate that both species are opportunistic breeders while black grouse tracks the climatic changes more slowly and is maintained at lower Ne . Our results highlight the usefulness of the psmc approach in investigating species' reaction to climate change in the deep past, but also that caution should be taken in drawing general conclusions about the recent past.