Genomic signatures of hybridization between Ixodes ricinus and Ixodes persulcatus in natural populations.

Identifying hybridization between common pathogen vectors is essential due to the major public health implications through risks associated with hybrid's enhanced pathogen transmission potential. The hard-ticks Ixodes ricinus and Ixodes persulcatus are the two most common vectors of tick-borne pathogens that affect human and animal health in Europe. Ixodes ricinus is a known native species in Finland with a well-known distribution, whereas I. persulcatus has expanded in range and abundance over the past 60 years, and currently it appears the most common tick species in certain areas in Finland. Here we used double-digest restriction site-associated DNA (ddRAD) sequencing on 186 ticks (morphologically identified as 92 I. ricinus, and 94 I. persulcatus) collected across Finland to investigate whether RAD generated single nucleotide polymorphisms (SNPs) can discriminate tick species and identify potential hybridization events. Two different clustering methods were used to assign specific species based on how they clustered and identified hybrids among them. We were able to discriminate between the two tick species and identified 11 putative hybrids with admixed genomic proportions ranging from approximately 24 to 76 percent. Four of these hybrids were morphologically identified as I. ricinus while the remaining seven were identified as I. persulcatus. Our results thus indicate that RAD SNPs are robust in identifying both species of the ticks as well as putative hybrids. These results further suggest ongoing hybridization between I. ricinus and I. persulcatus in their natural populations in Finland. The unique ability of RAD markers to discriminate between tick species and hybrids adds a useful aspect to tick evolutionary studies. Our findings align with previous studies and suggest a shared evolutionary history between the species, with instances of individuals possessing a considerable proportion of the other species' genome. This study is a significant step in understanding the formation of hybridization zones due to range expansion potentially associated with climate change.

First evidence of sex-specific responses to chemical cues in tardigrade mate searching behaviour.

Chemical cues are widely used in intraspecific and interspecific communication, either as substances deposited in the substrate or as molecules diffused in water or air. In tardigrades, an emerging microscopic study system, chemical communication and its role in reproduction are poorly known. Here, we assessed sex differences in the detection of (a) short-range diffusing signals and (b) deposited cue trails during the mate-searching behaviour of freely moving virgin male and female Macrobiotus polonicus. We tracked individual behaviour (a) in simultaneous double-choice chambers, where live conspecifics of each sex were presented in water and (b) of freely moving pairs on agar without water. We found that males, but not females, preferentially associated with opposite-sex individuals in trials conducted in water. In contrast, neither sex detected nor followed cues deposited on agar. In conclusion, our study suggests that mate discrimination and approach are male-specific traits and are limited to waterborne chemical cues. These results support the existence of Darwinian sex roles in pre-mating behaviour in an animal group with virtually non-existing sex differences in morphology or ecology.

Coinfection patterns of two marine apicomplexans are not associated with genetic diversity of their polychaete host.

Coinfections of two or more parasites within one host are more of a rule than an exception in nature. Interactions between coinfecting parasites can greatly affect their abundance and prevalence. Characteristics of the host, such as genetic diversity, can also affect the infection dynamics of coinfecting parasites. Here, we investigate for the first time the association of coinfection patterns of two marine apicomplexans, Rhytidocystis sp. and Selenidium pygospionis, with the genetic diversity of their host, the polychaete Pygospio elegans, from natural populations. Host genetic diversity was determined with seven microsatellite loci and summarized as allelic richness, inbreeding coefficient, and individual heterozygosity. We detected nonsignificant correlations between infection loads and both individual host heterozygosity and population genetic diversity. Prevalence and infection load of Rhytidocystis sp. were higher than those of S. pygospionis, and both varied spatially. Coinfections were common, and almost all hosts infected by S. pygospionis were also infected by Rhytidocystis sp. Rhytidocystis sp. infection load was significantly higher in dual infections. Our results suggest that factors other than host genetic diversity might be more important in marine apicomplexan infection patterns and experimental approaches would be needed to further determine how interactions between the apicomplexans and their host influence infection.

Species and genetic diversity relationships in benthic macroinvertebrate communities along a salinity gradient.

BACKGROUND: Species- and genetic diversity can change in parallel, resulting in a species-genetic diversity correlation (SGDC) and raising the question if the same drivers influence both biological levels of diversity. The SGDC can be either positive or negative, depending on whether the species diversity and the genetic diversity of the measured species respond in the same or opposite way to drivers. Using a traditional species diversity approach together with ultra-conserved elements and high throughput sequencing, we evaluated the SGDCs in benthic macrofauna communities in the Baltic Sea, a geologically young brackish water sea characterised by its steep salinity gradient and low species richness. Assessing SGDCs from six focal marine invertebrate species from different taxonomic groups and with differing life histories and ecological functions on both a spatial and temporal scale gives a more comprehensive insight into the community dynamics of this young ecosystem and the extrinsic factors that might drive the SGDCs. RESULTS: No significant correlations between species diversity and genetic diversity were found for any of the focal species. However, both negative and positive trends of SGDCs for the individual focal species were observed. When examining the environmental drivers, no common trends between the species were found, even when restricting the analysis to specific taxonomic classes. Additionally, there were no common environmental factors driving the diversity relationships for species sharing the same SGDC trend (positive or negative). Local population dynamics, together with the invasion history of the individual species and their unique adaptation to the distinctive environment of the Baltic Sea, are expected to be of major influence on the outcome of the SGDCs. CONCLUSIONS: The present results highlight the importance of assessing SGDCs using multiple species, not just a single indicator species. This emphasises a need to pay attention to the ecology and life history of the focal species. This study also provides insight into the large differences in both patterns and drivers of genetic diversity, which is important when including genetic biodiversity in conservation plans. We conclude that the effects of environmental and biological factors and processes that affects diversity patterns at both the community and genetic levels are likely species dependent, even in an environment such as the Baltic Sea with strong environmental gradients.

De novo transcriptome assembly of a facultative parasitic nematode Pelodera (syn. Rhabditis) strongyloides.

Pelodera strongyloides is a generally free-living gonochoristic facultative nematode. The whole genomic sequence of P. strongyloides remains unknown but 4 small subunit ribosomal RNA (ssrRNA) gene sequences are available. This project launched a de novo transcriptome assembly with 100 bp paired-end RNA-seq reads from normal, starved and wet-plate cultured animals. Trinity assembly tool generated 104,634 transcript contigs with N50 contig being 2195 bp and average contig length at 1103 bp. Transcriptome BLASTX matching results of five nematodes (C. elegans, Strongyloides stercoralis, Necator americanus, Trichuris trichiura, and Pristionchus pacificus) were consistent with their evolutionary relationships. Sixteen genes were identified to be homologous to key elements of the C. elegans RNA interference system, such as Dicer, Argonaute, RNA-dependent RNA polymerase and double strand RNA transport proteins. In starved samples, we observed up-regulation of cuticle related genes and 3 dauer formation genes. Dauer morphology was captured with enlarged phasmid under light microscopy, and dauer and normal larvae counts in clumps had a Pearson's product-moment correlation of 0.805 with P-value = 0.0088. Our results demonstrate that P. strongyloides could be used for studying nematode-related human or pet parasitic diseases. The sequenced assembled transcriptome reported here may be useful to understand the evolution of parasitism in Nematoda.

miRToolsGallery: a tag-based and rankable microRNA bioinformatics resources database portal.

Database URL: http://www.mirtoolsgallery.org.

Seasonal genetic variation associated with population dynamics of a poecilogonous polychaete worm.

Poecilogonous species show variation in developmental mode, with larvae that differ both morphologically and ecologically. The spionid polychaete Pygospio elegans shows variation in developmental mode not only between populations, but also seasonally within populations. We investigated the consequences of this developmental polymorphism on the spatial and seasonal genetic structure of P. elegans at four sites in the Danish Isefjord-Roskilde-Fjord estuary at six time points, from March 2014 until February 2015. We found genetic differentiation between our sampling sites as well as seasonal differentiation at two of the sites. The seasonal genetic shift correlated with the appearance of new size cohorts in the populations. Additionally, we found that the genetic composition of reproductive individuals did not always reflect the genetic composition of the entire sample, indicating that variance in reproductive success among individuals is a likely explanation for the patterns of chaotic genetic patchiness observed during this and previous studies. The heterogeneous, unpredictable character of the estuary might maintain poecilogony in P. elegans as a bet-hedging strategy in the Isefjord-Roskilde-Fjord complex in comparison with other sites where P. elegans are expected to be fixed to a certain mode of development.

De novo transcriptome assembly and developmental mode specific gene expression of Pygospio elegans.

Species with multiple different larval developmental modes are interesting models for the study of mechanisms underlying developmental mode transitions and life history evolution. Pygospio elegans, a small, tube-dwelling polychaete worm commonly found in estuarine and marine habitats around the northern hemisphere, is one species with variable developmental modes. To provide new genomic resources for studying P. elegans and to address the differences in gene expression between individuals producing offspring with different larval developmental modes, we performed whole transcriptome Illumina RNA sequencing of adult worms from two populations and prepared a de novo assembly of the P. elegans transcriptome. The transcriptome comprises 66,233 unigenes, of which 33,807 contain predicted coding sequences, 26,448 have at least one functional annotation, and 3,076 are classified as putative long non-coding RNAs. We found more than 8,000 unigenes significantly differentially expressed between adult worms from populations producing either planktonic or benthic larvae. This comprehensive transcriptome resource for P. elegans adds to the available genomic data for annelids and can be used to uncover mechanisms allowing developmental variation in this and potentially other marine invertebrate species.

Variation in gene expression within clones of the earthworm Dendrobaena octaedra.

Gene expression is highly plastic, which can help organisms to both acclimate and adapt to changing environments. Possible variation in gene expression among individuals with the same genotype (among clones) is not widely considered, even though it could impact the results of studies that focus on gene expression phenotypes, for example studies using clonal lines. We examined the extent of within and between clone variation in gene expression in the earthworm Dendrobaena octaedra, which reproduces through apomictic parthenogenesis. Five microsatellite markers were developed and used to confirm that offspring are genetic clones of their parent. After that, expression of 12 genes was measured from five individuals each from six clonal lines after exposure to copper contaminated soil. Variation in gene expression was higher over all genotypes than within genotypes, as initially assumed. A subset of the genes was also examined in the offspring of exposed individuals in two of the clonal lines. In this case, variation in gene expression within genotypes was as high as that observed over all genotypes. One gene in particular (chymotrypsin inhibitor) also showed significant differences in the expression levels among genetically identical individuals. Gene expression can vary considerably, and the extent of variation may depend on the genotypes and genes studied. Ensuring a large sample, with many different genotypes, is critical in studies comparing gene expression phenotypes. Researchers should be especially cautious inferring gene expression phenotypes when using only a single clonal or inbred line, since the results might be specific to only certain genotypes.

DNA Methylation and Potential for Epigenetic Regulation in Pygospio elegans.

Transitions in developmental mode are common evolutionarily, but how and why they occur is not understood. Developmental mode describes larval phenotypes, including morphology, ecology and behavior of larvae, which typically are generalized across different species. The polychaete worm Pygospio elegans is one of few species polymorphic in developmental mode, with multiple larval phenotypes, providing a possibility to examine the potential mechanisms allowing transitions in developmental mode. We investigated the presence of DNA methylation in P. elegans, and, since maternal provisioning is a key factor determining eventual larval phenotype, we compared patterns of DNA methylation in females during oogenesis in this species. We demonstrate that intragenic CpG site DNA methylation and many relevant genes necessary for DNA methylation occur in P. elegans. Methylation-sensitive AFLP analysis showed that gravid females with offspring differing in larval developmental mode have significantly different methylation profiles and that the females with benthic larvae and non-reproductive females from the same location also differ in their epigenetic profiles. Analysis of CpG sites in transcriptome data supported our findings of DNA methylation in this species and showed that CpG observed/expected ratios differ among females gravid with embryos destined to different developmental modes. The differences in CpG site DNA methylation patterns seen among the samples suggest a potential for epigenetic regulation of gene expression (through DNA methylation) in this species.

Low levels of mitochondrial DNA and symbiont diversity in the worldwide agricultural pest, the greenhouse whitefly Trialeurodes vaporariorum (Hemiptera: Aleyrodidae).

Trialeurodes vaporariorum, the greenhouse whitefly, is a cosmopolitan agricultural pest. Little is known about the genetic diversity of T. vaporariorum and the bacterial symbionts associated with this species. Here, we undertook a large phylogeographic study by investigating both the mitochondrial (mt) diversity and the infection status of 38 T. vaporariorum collections from 18 countries around the world. Genetic diversity of T. vaporariorum was studied by analyzing sequence data from the mt cytochrome oxidase I, cytochrome b, and NADH dehydrogenase subunit 5 genes. Maximum-likelihood (ML) phylogeny reconstruction delineated 2 clades characterized by limited sequence divergence: one clade comprised samples only from the Northern hemisphere whereas the other comprised samples from a broader geographical range. The presence of secondary symbionts was determined by PCR using primers specific for Hamiltonella, Rickettsia, Arsenophonus, Cardinium, Wolbachia, and Fritschea. Most individuals examined harbored at least one secondary endosymbiont, and Arsenophonus was detected in almost all male and female individuals. Wolbachia was present at a much lower frequency, and Cardinium was detected in only a few individuals from Greece. Rickettsia, Hamiltonella, and Fritschea were not found. Additionally, we set out to further analyze Arsenophonus diversity by multilocus sequence typing analysis; however, the Arsenophonus sequences did not exhibit any polymorphism. Our results revealed remarkably low diversity in both mtDNA and symbionts in this worldwide agricultural pest, contrasting sharply with that of the ecologically similar Bemisia tabaci.

Interspecific interactions influence contrasting spatial genetic structures in two closely related damselfly species.

Spatial genetic structure (SGS) is largely determined by colonization history, landscape and ecological characteristics of the species. Therefore, sympatric and ecologically similar species are expected to exhibit similar SGSs, potentially enabling prediction of the SGS of one species from that of another. On the other hand, due to interspecific interactions, ecologically similar species could have different SGSs. We explored the SGSs of the closely related Calopteryx splendens and Calopteryx virgo within Finland and related the genetic patterns to characteristics of the sampling localities. We observed different SGSs for the two species. Genetic differentiation even within short distances in C. splendens suggests genetic drift as an important driver. However, we also observed indication of previous gene flow (revealed by a negative relationship between genetic differentiation and increasing potential connectivity of the landscape). Interestingly, genetic diversity of C. splendens was negatively related to density of C. virgo, suggesting that interspecific interactions influence the SGS of C. splendens. In contrast, genetic differentiation between C. virgo subpopulations was low and only exhibited relationships with latitude, pointing to high gene flow, colonization history and range margin effects as the drivers of SGS. The different SGSs of the two ecologically similar species caution indirect inferences of SGS based on ecologically similar surrogate species.

Temporal genetic structure in a poecilogonous polychaete: the interplay of developmental mode and environmental stochasticity.

BACKGROUND: Temporal variation in the genetic structure of populations can be caused by multiple factors, including natural selection, stochastic environmental variation, migration, or genetic drift. In benthic marine species, the developmental mode of larvae may indicate a possibility for temporal genetic variation: species with dispersive planktonic larvae are expected to be more likely to show temporal genetic variation than species with benthic or brooded non-dispersive larvae, due to differences in larval mortality and dispersal ability. We examined temporal genetic structure in populations of Pygospio elegans, a poecilogonous polychaete with within-species variation in developmental mode. P. elegans produces either planktonic, benthic, or intermediate larvae, varying both among and within populations, providing a within-species test of the generality of a relationship between temporal genetic variation and larval developmental mode. RESULTS: In contrast to our expectations, our microsatellite analyses of P. elegans revealed temporal genetic stability in the UK population with planktonic larvae, whereas there was variation indicative of drift in temporal samples of the populations from the Baltic Sea, which have predominantly benthic and intermediate larvae. We also detected temporal variation in relatedness within these populations. A large temporal shift in genetic structure was detected in a population from the Netherlands, having multiple developmental modes. This shift could have been caused by local extiction due to extreme environmental conditions and (re)colonization by planktonic larvae from neighboring populations. CONCLUSIONS: In our study of P. elegans, temporal genetic variation appears to be due to not only larval developmental mode, but also the stochastic environment of adults. Large temporal genetic shifts may be more likely in marine intertidal habitats (e.g. North Sea and Wadden Sea) which are more prone to environmental stochasticity than the sub-tidal Baltic habitats. Sub-tidal and/or brackish (less saline) habitats may support smaller P. elegans populations and these may be more susceptible to the effects of random genetic drift. Moreover, higher frequencies of asexual reproduction and the benthic larval developmental mode in these populations leads to higher relatedness and contributes to drift. Our results indicate that a general relationship between larval developmental mode and temporal genetic variation may not exist.

The effect of inbreeding rate on fitness, inbreeding depression and heterosis over a range of inbreeding coefficients.

Understanding the effects of inbreeding and genetic drift within populations and hybridization between genetically differentiated populations is important for many basic and applied questions in ecology and evolutionary biology. The magnitudes and even the directions of these effects can be influenced by various factors, especially by the current and historical population size (i.e. inbreeding rate). Using Drosophila littoralis as a model species, we studied the effect of inbreeding rate over a range of inbreeding levels on (i) mean fitness of a population (relative to that of an outbred control population), (ii) within-population inbreeding depression (reduction in fitness of offspring from inbred versus random mating within a population) and (iii) heterosis (increase in fitness of offspring from interpopulation versus within-population random mating). Inbreeding rate was manipulated by using three population sizes (2, 10 and 40), and fitness was measured as offspring survival and fecundity. Fast inbreeding (smaller effective population size) resulted in greater reduction in population mean fitness than slow inbreeding, when populations were compared over similar inbreeding coefficients. Correspondingly, populations with faster inbreeding expressed more heterosis upon interpopulation hybridization. Inbreeding depression within the populations did not have a clear relationship with either the rate or the level of inbreeding.

Polymorphism in developmental mode and its effect on population genetic structure of a spionid polychaete, Pygospio elegans.

Population genetic structure of sedentary marine species is expected to be shaped mainly by the dispersal ability of their larvae. Long-lived planktonic larvae can connect populations through migration and gene flow, whereas species with nondispersive benthic or direct-developing larvae are expected to have genetically differentiated populations. Poecilogonous species producing different larval types are ideal when studying the effect of developmental mode on population genetic structure and connectivity. In the spionid polychaete Pygospio elegans, different larval types have been observed between, and sometimes also within, populations. We used microsatellite markers to study population structure of European P. elegans from the Baltic Sea (BS) and North Sea (NS). We found that populations with planktonic larvae had higher genetic diversity than did populations with benthic larvae. However, this pattern may not be related to developmental mode, since in P. elegans, developmental mode may be associated with geography. Benthic larvae were more commonly seen in the brackish BS and planktonic larvae were predominant in the NS, although both larval types also are found from both areas. Significant isolation-by-distance (IBD) was found overall and within regions. Most of the pair-wise F(ST) comparisons among populations were significant, although some geographically close populations with planktonic larvae were found to be genetically similar. However, these results, together with the pattern of IBD, autocorrelation within populations, as well as high estimated local recruitment, suggest that dispersal is limited in populations with planktonic larvae as well as in those with benthic larvae. The decrease in salinity between the NS and BS causes a barrier to gene flow in many marine species. In P. elegans, low, but significant, differentiation was detected between the NS and BS (3.34% in AMOVA), but no clear transition zone was observed, indicating that larvae are not hampered by the change in salinity.

Introduction to symposium: Poecilogony--a window on larval evolutionary transitions in marine invertebrates.

Poecilogony is the intraspecific variation in developmental mode that has been described in some marine invertebrates. Poecilogonous species produce different larval forms (e.g., free-swimming planktotrophic larvae as well as brooded lecithotrophic or adelphophagic larvae). Poecilogony can be a controversial topic, since it is difficult to identify and characterize the phenomenon with certainty. It has been challenging to determine whether poecilogony represents developmental polymorphism with a genetic basis or developmental polyphenism reflecting plastic responses to environmental cues. Other outstanding questions include whether common mechanisms underlie the developmental variation we observe in poecilogonous species, and whether poecilogony is maintained in different taxa through similar mechanisms or selective pressures. Poecilogonous species provide a unique opportunity to elucidate the cellular, developmental, and genetic mechanisms underlying evolutionary transitions in developmental mode, as well as to help clarify the selective pressures and possible ecological circumstances that might be involved. Here, we describe an integrative approach to the study of poecilogony and its role in larval evolutionary transitions highlighted during a symposium held at the 2012 annual meeting of the Society for Integrative and Comparative Biology.

What do we need to know about speciation?

Speciation has been a major focus of evolutionary biology research in recent years, with many important advances. However, some of the traditional organising principles of the subject area no longer provide a satisfactory framework, such as the classification of speciation mechanisms by geographical context into allopatric, parapatric and sympatry classes. Therefore, we have asked where speciation research should be directed in the coming years. Here, we present a distillation of questions about the mechanisms of speciation, the genetic basis of speciation and the relationship between speciation and diversity. Our list of topics is not exhaustive; rather we aim to promote discussion on research priorities and on the common themes that underlie disparate speciation processes.

Adaptation to a seasonally varying environment: a strong latitudinal cline in reproductive diapause combined with high gene flow in Drosophila montana.

Adaptation to seasonal changes in the northern hemisphere includes an ability to predict the forthcoming cold season from gradual changes in environmental cues early enough to prepare for the harsh winter conditions. The magnitude and speed of changes in these cues vary between the latitudes, which induces strong selection pressures for local adaptation.We studied adaptation to seasonal changes in Drosophila montana, a northern maltfly, by defining the photoperiodic conditions leading to adult reproductive diapause along a latitudinal cline in Finland and by measuring genetic differentiation and the amount of gene flow between the sampling sites with microsatellites. Our data revealed a clear correlation between the latitude and the critical day length (CDL), in which half of the females of different cline populations enter photoperiodic reproductive diapause. There was no sign of limited gene flow between the cline populations, even though these populations showed isolation by distance. Our results show that local adaptation may occur even in the presence of high gene flow, when selection for locally adaptive life-history traits is strong. A wide range of variation in the CDLs of the fly strains within and between the cline populations may be partly due to gene flow and partly due to the opposing selection pressures for fly reproduction and overwinter survival. This variation in the timing of diapause will enhance populations' survival over the years that differ in the severity of the winter and in the length of the warm period and may also help them respond to long-term changes in environmental conditions.

Genetic variability and drift load in populations of an aquatic snail.

Population genetic theory predicts that in small populations, random genetic drift will fix and accumulate slightly deleterious mutations, resulting in reduced reproductive output. This genetic load due to random drift (i.e., drift load) can increase the extinction risk of small populations. We studied the relationship between genetic variability (indicator of past population size) and reproductive output in eight isolated, natural populations of the hermaphroditic snail Lymnaea stagnalis. In a common laboratory environment, snails from populations with the lowest genetic variability mature slower and have lower fecundity than snails from genetically more variable populations. This result suggests that past small population size has resulted in increased drift load, as predicted. The relationship between genetic variability and reproductive output is independent of the amount of nonrandom mating within populations. However, reproductive output and the current density of snails in the populations were not correlated. Instead, data from the natural populations suggest that trematode parasites may determine, at least in part, population densities of the snails.

The effects of mating system and genetic variability on susceptibility to trematode parasites in a freshwater snail, Lymnaea stagnalis.

The amount and distribution of genetic variability in host populations can have significant effects on the outcome of host-parasite interactions. We studied the effect of mating system and genetic variability on susceptibility of Lymnaea stagnalis snails to trematode parasites. Mating system of snails from eight populations differing in the amount of genetic variability was manipulated, and self- and cross-fertilized offspring were exposed to naturally occurring trematode parasites in a controlled lake experiment. Susceptibility of snails varied between populations, but mating-system treatment did not have a significant effect. Heterozygosity of snails was negatively correlated with the probability of trematode infection, however, suggesting that parasitic diseases may pose a serious threat to populations lacking genetic variability.