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.

Non-invasive assessment of ventriculo-arterial coupling using aortic wave intensity analysis combining central blood pressure and phase-contrast cardiovascular magnetic resonance.

BACKGROUND: Wave intensity analysis (WIA) in the aorta offers important clinical and mechanistic insight into ventriculo-arterial coupling, but is difficult to measure non-invasively. We performed WIA by combining standard cardiovascular magnetic resonance (CMR) flow-velocity and non-invasive central blood pressure (cBP) waveforms. METHODS AND RESULTS: Two hundred and six healthy volunteers (age range 21-73 years, 47% male) underwent sequential phase contrast CMR (Siemens Aera 1.5 T, 1.97 × 1.77 mm2, 9.2 ms temporal resolution) and supra-systolic oscillometric cBP measurement (200 Hz). Velocity (U) and central pressure (P) waveforms were aligned using the waveform foot, and local wave speed was calculated both from the PU-loop (c) and the sum of squares method (cSS). These were compared with CMR transit time derived aortic arch pulse wave velocity (PWVtt). Associations were examined using multivariable regression. The peak intensity of the initial compression wave, backward compression wave, and forward decompression wave were 69.5 ± 28, -6.6 ± 4.2, and 6.2 ± 2.5 × 104 W/m2/cycle2, respectively; reflection index was 0.10 ± 0.06. PWVtt correlated with c or cSS (r = 0.60 and 0.68, respectively, P < 0.01 for both). Increasing age decade and female sex were independently associated with decreased forward compression wave (-8.6 and -20.7 W/m2/cycle2, respectively, P < 0.01) and greater wave reflection index (0.02 and 0.03, respectively, P < 0.001). CONCLUSION: This novel non-invasive technique permits straightforward measurement of wave intensity at scale. Local wave speed showed good agreement with PWVtt, and correlation was stronger using the cSS than the PU-loop. Ageing and female sex were associated with poorer ventriculo-arterial coupling in healthy individuals.

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.

Evolutionary conservation and function of the human embryonic stem cell specific miR-302/367 cluster.

miRNA clusters define a group of related miRNAs closely localized in the genome with an evolution that remains poorly understood. The miR-302/367 cluster represents a single polycistronic transcript that produces five precursor miRNAs. The cluster is highly expressed and essential for maintenance of human embryonic stem cells. We found the cluster to be highly conserved and present in most mammals. In primates, seed sequence and miRNA structure are conserved, but inter-precursor sequences are evolving. Insertions of new miRNAs, deletions of individual miRNAs, and a cluster duplication observed in different species suggest an actively evolving cluster. Core transcriptional machinery consisting of NANOG and OCT-4 transcription factors that define stem cells are present upstream of the miR-302/367 cluster. Interestingly, we found the miR-302/367 cluster flanking region to be enriched as a target site of other miRNAs suggesting a mechanism for feedback control. Analysis of miR-302 and miR-367 targets demonstrated concordance of gene set enrichment groups at high gene ontology levels. This cluster also expresses isomiRs providing another means of establishing sequence diversity. Finally, using three different kidney tumor datasets, we observed consistent expression of miR-302 family members in normal tissue while adjacent tumor tissue showed a significant lack of expression. Clustering expression levels of miR-302 validated target genes showed a significant correlation between miR-302/367 cluster miRNAs and a subset of validated gene targets in healthy and adjacent tumor tissues. Taken together, our data show a highly conserved and still evolving miRNA cluster that may have additional unrecognized functions.

Nutrient and mineral composition during shoot growth in seven species of Phyllostachys and Pseudosasa bamboo consumed by giant panda.

During the annual period of bamboo shoot growth in spring, free-ranging giant pandas feed almost exclusively on the shoots while ignoring the leaves and full- height culm. Little is known about the nutritional changes that occur during bamboo shoot growth, if nutritional changes differ among species, or how these changes might influence forage selection. Our objective was to examine the nutrient and mineral composition during three phases of shoot growth (<60, 90-150 and >180 cm) for seven species of bamboo (Phyllostachys (P.) aurea, P. aureosulcata, P. bissetii, P. glauca, P. nuda, P. rubromarginata, Pseudosasa japonica) fed to captive giant pandas at the Memphis Zoo. Total dietary fiber content of bamboo shoots increased (p < 0.0001) from an overall species average of 61% dry matter (DM) at < 60 cm to 75% DM at shoot heights > 180 cm, while crude protein, fat and ash exhibited significant declines (p < 0.05). Phyllostachys nuda had the overall greatest (p = 0.007) crude protein (21% DM) and fat (4% DM) content, and lowest overall total fibre (61% DM) content compared to the other species examined. In contrast, Pseudosasa japonica had the overall lowest crude protein and fat, and relatively higher fibre content (9%, 3% and 74% respectively). Concentrations of Zn and Fe were highest in shoots <60 cm (10-50 µg/g DM) and decreased (p < 0.05) during growth in all species examined. Concentrations of Ca, Cu, Mn, Na and K varied among species and were largely unaffected by growth stage. Due to their higher concentrations of nutrients and lower fibre content in comparison to culm and leaf, bamboo shoots should be a major component of captive giant panda diets when available.

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.

Metallothionein gene expression differs in earthworm populations with different exposure history.

Metals are persistent pollutants in soils that can harm soil organisms and decrease species diversity. Animals can cope with metal contamination with the help of metallothioneins, small metal-binding proteins involved in homeostasis and detoxification of metals. We studied the expression of metallothionein with qPCR in a small, epigeic earthworm, Dendrobaena octaedra. We compared expression patterns and metal body content in earthworms collected from two sites with different metal contamination histories: Harjavalta, contaminated by a Cu-Ni smelter operational for over 50 years, and Jyväskylä, an uncontaminated site. Earthworms from both sites were also experimentally exposed to different concentrations of Cu (control, 50, 100 or 200 mg/kg) or Zn (control, 75, 150 or 300 mg/kg) for 7, 14 or 28 days to determine if there is a time related dose-response in gene expression. Population comparison showed that metallothionein expression was higher in earthworms from the contaminated site. In the exposure experiment, exposure time affected expression, but only in the earthworms from the uncontaminated site, suggesting that there is a delay in the metallothionein response of earthworms in this population. In contrast, earthworms from the contaminated site showed higher and constant levels of metallothionein expression at all exposure concentrations and durations. The constant metallothionein expression in earthworms from the contaminated site suggests that inducibility of metallothionein response could be lost in earthworms with metal exposure history. Adaptation of D. octaedra to metal exposure could explain the differences between the populations and explain the persistence of this species in contaminated forest soils.

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.

The benefits of interpopulation hybridization diminish with increasing divergence of small populations.

Interpopulation hybridization can increase the viability of small populations suffering from inbreeding and genetic drift, but it can also result in outbreeding depression. The outcome of hybridization can depend on various factors, including the level of genetic divergence between the populations, and the number of source populations. Furthermore, the effects of hybridization can change between generations following the hybridization. We studied the effects of population divergence (low vs. high level of divergence) and the number of source populations (two vs. four source populations) on the viability of hybrid populations using experimental Drosophila littoralis populations. Population viability was measured for seven generations after hybridization as proportion of populations facing extinction and as per capita offspring production. Hybrid populations established at the low level of population divergence were more viable than the inbred source populations and had higher offspring production than the large control population. The positive effects of hybridization lasted for the seven generations. In contrast, at the high level of divergence, the viability of the hybrid populations was not significantly different from the inbred source populations, and offspring production in the hybrid populations was lower than in the large control population. The number of source populations did not have a significant effect at either low or high level of population divergence. The study shows that the benefits of interpopulation hybridization may decrease with increasing divergence of the populations, even when the populations share identical environmental conditions. We discuss the possible genetic mechanisms explaining the results and address the implications for conservation of populations.

Inbreeding rate modifies the dynamics of genetic load in small populations.

The negative fitness consequences of close inbreeding are widely recognized, but predicting the long-term effects of inbreeding and genetic drift due to limited population size is not straightforward. As the frequency and homozygosity of recessive deleterious alleles increase, selection can remove (purge) them from a population, reducing the genetic load. At the same time, small population size relaxes selection against mildly harmful mutations, which may lead to accumulation of genetic load. The efficiency of purging and the accumulation of mutations both depend on the rate of inbreeding (i.e., population size) and on the nature of mutations. We studied how increasing levels of inbreeding affect offspring production and extinction in experimental Drosophila littoralis populations replicated in two sizes, N = 10 and N = 40. Offspring production and extinction were measured over 25 generations concurrently with a large control population. In the N = 10 populations, offspring production decreased strongly at low levels of inbreeding, then recovered only to show a consistent subsequent decline, suggesting early expression and purging of recessive highly deleterious alleles and subsequent accumulation of mildly harmful mutations. In the N = 40 populations, offspring production declined only after inbreeding reached higher levels, suggesting that inbreeding and genetic drift pose a smaller threat to population fitness when inbreeding is slow. Our results suggest that highly deleterious alleles can be purged in small populations already at low levels of inbreeding, but that purging does not protect the small populations from eventual genetic deterioration and extinction.