Genome-wide single-generation signatures of local selection in the panmictic European eel.

Next-generation sequencing and the collection of genome-wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD-sequenced European eel individuals (glass eels) from eight locations between 34 and 64(o) N, we examined the patterns of genome-wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F(ST) -based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F(ST) = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single-generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand-receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.

Assessing patterns of hybridization between North Atlantic eels using diagnostic single-nucleotide polymorphisms.

The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

A resource of genome-wide single-nucleotide polymorphisms generated by RAD tag sequencing in the critically endangered European eel.

Reduced representation genome sequencing such as restriction-site-associated DNA (RAD) sequencing is finding increased use to identify and genotype large numbers of single-nucleotide polymorphisms (SNPs) in model and nonmodel species. We generated a unique resource of novel SNP markers for the European eel using the RAD sequencing approach that was simultaneously identified and scored in a genome-wide scan of 30 individuals. Whereas genomic resources are increasingly becoming available for this species, including the recent release of a draft genome, no genome-wide set of SNP markers was available until now. The generated SNPs were widely distributed across the eel genome, aligning to 4779 different contigs and 19,703 different scaffolds. Significant variation was identified, with an average nucleotide diversity of 0.00529 across individuals. Results varied widely across the genome, ranging from 0.00048 to 0.00737 per locus. Based on the average nucleotide diversity across all loci, long-term effective population size was estimated to range between 132,000 and 1,320,000, which is much higher than previous estimates based on microsatellite loci. The generated SNP resource consisting of 82,425 loci and 376,918 associated SNPs provides a valuable tool for future population genetics and genomics studies and allows for targeting specific genes and particularly interesting regions of the eel genome.

Altitudinal and seasonal variation in microsatellite allele frequencies of Drosophila buzzatii.

Variation in climate, particularly temperature, is known to affect the genetic composition of populations. Although there have been many studies of latitudinal variation, comparisons of populations across altitudes or seasons, particularly for animal species, are less common. Here, we study genetic variation (microsatellite markers) in populations of Drosophila buzzatii collected along altitudinal gradients and in different seasons. We found no differences in genetic variation between 2 years or between seasons within years. However, there were numerous cases of significant associations between allele frequencies or expected heterozygosities and altitude, with more than half showing nonlinear relationships. While these associations indicate possible selection and local altitudinal adaptation, direct tests gave strong evidence for selection affecting two loci and weaker evidence for five other loci. Two loci that are located within an inversion (including the one with strongest evidence for selection) show a linear increase in genetic diversity with altitude, likely due to thermal selection. Parallel associations with altitude here and with latitude in Australian populations indicate that selection is operating on chromosomal regions marked by some of the loci.

Bottlenecks, population differentiation and apparent selection at microsatellite loci in Australian Drosophila buzzatii.

Species colonizing new areas disjunct from their original habitat may be subject to novel selection pressures, and exhibit adaptive genetic changes. However, if colonization occurs through a small number of founders, the genetic composition of the colonized population may differ from that of the original population simply due to genetic drift. Disentangling the effects of founder drift and selection after colonization is crucial to understanding the adaptive process. Drosophila buzzatii colonized Australia some 600-700 generations ago, and spread rapidly over a wide geographical range. Genetic variation for 15 microsatellite loci in each of nine populations in eastern Australia was used to estimate the size of the bottleneck, and to determine if any of these microsatellites marked genomic regions subject to recent selection. We estimate that on its introduction to Australia, D. buzzatii went through a moderate bottleneck (approximately 30-40 founders). Linkage disequilibrium was common, both intrachromosomal and between loci on different chromosomes. Of the 15 loci, 2 showed evidence of selection, one exhibiting local adaptation in different populations and the other balancing selection. We conclude that linkage disequilibria may be far more common in natural populations than is generally assumed, and the loci apparently affected by selection may well be marking selection in large genome regions including many loci that are not necessarily closely linked.

Distance education and its potential for international co-operative education.

This paper summarises some of the major trends in distance education in the first decade of the 21st Century, and explores the implications of these trends for international collaboration among institutions of veterinary medicine.

Local adaptation of stress related traits in Drosophila buzzatii and Drosophila simulans in spite of high gene flow.

We addressed the question if local adaptation to a thermal gradient is possible in spite of a high gene flow among closely spaced populations of two species of Drosophila from the island of La Gomera (Canary Islands). Variation in multiple traits related to stress resistance in different life stages was measured in both species in flies collected from five localities at different altitudes and thereby with different climatic conditions. Based on microsatellite loci, the populations were not genetically differentiated. However, 18 of the 24 independent traits measured showed significant differentiation among populations of Drosophila buzzatii, but only nine of 25 for Drosophila simulans. This difference in the number of traits might reflect higher habitat specificity and thus higher potential for local adaptation of D. buzzatii than D. simulans. We found clinal variation, as some traits showed significant linear regressions on altitude, but more on altitude cubed.

Nucleotide diversity in the Hsp90 gene in natural populations of Drosophila melanogaster from Australia.

Hsp90 is regarded as one of the best candidates for an evolved mechanism that regulates the expression of genetic and phenotypic variability. We examined nucleotide diversity in both the promoter and coding regions of Hsp90, the gene which encodes Hsp90 in Drosophila, in natural populations of Drosophila melanogaster from eastern Australia. We found that Hsp90 is polymorphic for only two nonsynonymous changes in the coding region, both of which are deletions of a lysine residue. One of these lysine deletions was in complete linkage disequilibrium with the inversion In(3L)P, and showed a significant association with latitude. The other lysine deletion reported here for the first time varied from 0 to 15% in natural populations, but did not show a clinal pattern. The regulatory and coding regions of Hsp90 showed very low nucleotide diversity compared to other nuclear genes, and chromosomes containing In(3L)P had lower levels of nucleotide diversity than the standard arrangements. Non-neutral evolution of Hsp90 was not supported by analyses of either the regulatory or coding regions of the gene. These results are discussed within the context of Hsp90 variation being involved in thermotolerance as well as the expression of genetic and phenotypic variability.

Adaptive divergence in a high gene flow environment: Hsc70 variation in the European flounder (Platichthys flesus L.).

Little is known about local adaptations in marine fishes since population genetic surveys in these species have typically not applied genetic markers subject to selection. In this study, we used a candidate gene approach to investigate adaptive population divergence in the European flounder (Platichthys flesus L.) throughout the northeastern Atlantic. We contrasted patterns of genetic variation in a presumably neutral microsatellite baseline to patterns from a heat-shock cognate protein gene, Hsc70. Using two different neutrality tests we found that the microsatellite data set most likely represented a neutral baseline. In contrast, Hsc70 strongly deviated from neutral expectations. Importantly, when estimating standardized levels of population divergence (F(ST)'), we also found a large discrepancy in the patterns of structuring in the two data sets. Thus, samples grouped according to geographical or historical proximity with regards to microsatellites, but according to environmental similarities with regards to Hsc70. The differences between the data sets were particularly pronounced in pairwise comparisons involving populations in the western and central Baltic Sea. For instance, the genetic differentiation between geographically close Baltic Sea and North Sea populations was found to be 0.02 and 0.45 for microsatellites and Hsc70 respectively. Our results strongly suggest adaptive population divergence and indicate local adaptations at the DNA level in a background of high levels of gene flow, typically found in many marine fish species. Furthermore, this study highlights the usefulness of the candidate gene approach for demonstrating local selection in non-model organisms such as most marine fishes.

Microsatellites provide insight into contrasting mating patterns in arribada vs. non-arribada olive ridley sea turtle rookeries.

Molecular studies of sea turtles have shown that the frequency of multiple paternity (MP) varies between species, and between rookeries of the same species. This study uses nuclear microsatellite markers to compare the incidence of MP in two neighbouring olive ridley rookeries on the Pacific coast of Costa Rica, with contrasting nesting behaviours -- the 'arribada' population nesting at Ostional and the solitary nesters of Playa Hermosa. Using two highly polymorphic microsatellite markers, we tested 13 nests from each location and found a significant difference (P < 0.001) between the level of MP of the arribada rookery (92%- the highest found for marine turtles) and that of the solitary nesting rookery (30%). Additional analyses based on six microsatellite loci revealed no genetic differentiation between nesting females from the two locations, or between nesting females and attendant males from the Ostional breeding area. Sixty-nine per cent of the nests with MP were fathered by a minimum of three different males, and three nests showed evidence of at least four fathers. The results suggest that the differences observed in levels of MP between arribada and solitary rookeries are due to an effect of abundance of individuals on the mating system. This is supported by a regression analysis combining other paternity studies on sea turtles which shows that levels of MP increase with increasing abundance of nesting females.

Use of ITS rDNA for discriminating of larval stages of two microphallid (Digenea) species using Hydrobia ulvae (Pennant, 1777) and Corophium volutator (Pallas, 1766) as intermediate hosts.

Digenean trematodes encompass several species with little morphological differentiation in the larval stages and, as a result, uncertainty prevails regarding species identification. The microphallid trematode Maritrema subdolum occurs widespread geographically in mud snail and crustacean hosts in European marine shallow-water ecosystems. Larval stages of this and other congeneric species are, however, difficult to separate morphologically. To verify the species status and to secure identification of two co-occurring microphallids, M. subdolum and microphallid sp. no. 15 (Deblock, 1980), we examined the nucleotide sequences of the internal transcribed spacer regions (ITS1, ITS2). From fragments consisting of both ITS regions and the 5.8S gene (nearly 1,200 bp), a sequence divergence of 2.9% between the two types was recorded. In accordance with the morphological traits of the cercariae (stylet shape, length), the results support the view that the two types actually represent different species. Species-specific primers were prepared for each species. They proved to be efficient diagnostic tools for identifying single larval-stage individuals. Using these primers, infections in host organisms were also verified without performing a dissection of the host individuals.

Remating and sperm displacement in a natural population of Drosophila buzzatii inferred from mother-offspring analysis of microsatellite loci.

Prospects for estimation of parameters of models of sperm competition from field data have improved recently with the development of methods that employ multilocus genotype data from brood-structured samples. Sperm competition in Drosophila buzzatii is of special interest because it is possible to directly observe the breeding behaviour of this species in its natural habitat of rotting cactus. Previous laboratory experiments showed that this species exhibits an unusual pattern of frequent remating and sperm partitioning. This paper reports the first attempt to estimate the frequency of female remating and sperm competition in natural populations of D. buzzatii. For the Australian population studied, the mean remating frequency was lower (alpha = 2.12-2.20) than previously estimated in laboratory experiments with the same population, whereas mean sperm displacement (beta = 0.69-0.71) fell within the limits of previous laboratory results. The evolution of the D. buzzatii mating system is discussed.

DNA sequence variation and latitudinal associations in hsp23, hsp26 and hsp27 from natural populations of Drosophila melanogaster.

Heat shock genes are considered to be likely candidate genes for environmental stress resistance. Nucleotide variation in the coding sequence of the small heat shock genes (hsps) hsp26 and hsp27 from Drosophila melanogaster was studied in flies originating from the Netherlands and eastern Australia. The hsp26 gene was polymorphic for an insertion/deletion of three extra amino acids and two nonsynonymous changes in all populations. The hsp27 gene exhibited two nonsynonymous changes and three synonymous mutations. The hsp26 polymorphism showed a latitudinal cline along the east coast of Australia. This pattern was not confounded by the fact that the shsps are located in the inversion In(3 L)P which also shows a latitudinal cline in eastern Australia. A similar latitudinal cline was found for the previously described variation in hsp23, while frequencies of hsp27 alleles did not change with latitude. These findings suggest that variation at two of the shsps or closely linked loci are under selection in natural populations of D. melanogaster.

Genetic variation in original and colonizing Drosophila buzzatii populations analysed by microsatellite loci isolated with a new PCR screening method.

A new polymerase chain reaction-based screening method for microsatellites is presented. Using this method, we isolated 12 microsatellite loci from Drosophila buzzatii, two of which were X-linked. We applied the other 10 microsatellite loci to the analysis of genetic variation in five natural populations of D. buzzatii. Two populations were from the species' original distribution in Argentina, whereas the other three were from Europe (two) and Australia that were colonized 200 and 65 years ago, respectively. Allelic variation was much larger in the original populations than in the colonizing ones and there was a tendency to decreased heterozygosity in the colonizing populations. We used three different statistical procedures for detecting population bottlenecks. All procedures suggested that the low variability in the populations in the Old World was not the result of the recent population decline, but was due to a founder effect followed by a population expansion. In fact, one procedure which detects population expansions and declines based on the genealogical history of microsatellite data suggested that an expansion had taken place in all the colonized populations.

Drosophila melanogaster is polymorphic for a specific repeated (CATA) sequence in the regulatory region of hsp23.

To identify sequence variation associated with a selection response for heat tolerance in Drosophila melanogaster, we sequenced 1400bp of the heat shock protein 23 gene (hsp23) promoter region in four heat-selected and two control lines. The region was found to be variable for a specific (CATA) repeated sequence, and the sequence CTT seems to be a hot spot for mutation. The repeated tetranucleotide sequence was located in several short repeats scattered throughout the entire region. Similar variable repeats are also located downstream the of hsp23 gene in the intergenic region between hsp23 and hsp27. We detected nine different hsp23 alleles. Their frequencies in the selection and control lines seemed to be mainly determined by genetic drift. The function of the CATA repeats is not yet known, though these regions have homology to SAR elements located in the intergenic region between two hsp70 genes, suggesting a similar function.

Evolutionary transition from single to multiple mating in fungus-growing ants.

Queens of leafcutter ants exhibit the highest known levels of multiple mating (up to 10 mates per queen) among ants. Multiple mating may have been selected to increase genetic diversity among nestmate workers, which is hypothesized to be critical in social systems with large, long-lived colonies under severe pressure of pathogens. Advanced fungus-growing (leafcutter) ants have large numbers (104-106 workers) and long-lived colonies, whereas basal genera in the attine tribe have small (< 200 workers) colonies with probably substantially shorter lifespans. Basal attines are therefore expected to have lower queen mating frequencies, similar to those found in most other ants. We tested this prediction by analysing queen mating frequency and colony kin structure in three basal attine species: Myrmicocrypta ednaella, Apterostigma collare and Cyphomyrmex longiscapus. Microsatellite marker analyses revealed that queens in all three species were single mated, and that worker-to-worker relatedness in these basal attine species is very close to 0.75, the value expected under exclusively single mating. Fungus growing per se has therefore not selected for multiple queen mating. Instead, the advanced and highly productive social structure of the higher attine ants, which is fully dependent on the rearing of an ancient clonal fungus, may have necessitated high genetic diversity among nestmate workers. This is not the case in the lower attines, which rear fungi that were more recently derived from free-living fungal populations.

CD19-selected B lymphocytes synthesize, secrete and migrate in the presence of IL-8. TNF-alpha and gammaIP-10 are also B lymphocyte migratory factors.

B lymphocytes are responsible for antigen uptake and presentation, as well as antibody production. These reactions require close cell-to-cell contact between B lymphocytes and monocytes. In this study we demonstrate that interleukin 8 (IL-8), gamma-immune protein 10 (gammaIP-10) and tumour necrosis factor alpha (TNF-alpha) all induce a significant chemokinetic response of human B lymphocytes. Among the cytokines tested, rIL-8 was the strongest B lymphocyte migratory factor with a migratory index (MI) of 2.03+/-0.32, (P<0.002), followed by rTNF-alpha (MI=1.89+/-0.17, P<0.001) and rgammaIP-10 (MI=1.63+/-0.17, P<0.001). We did not observe B lymphocyte migration towards rIL-1alpha, rIL-2, rIL-4, rIL-10, interferon gamma (rINF-gamma) or transforming growth factor beta (rTGF-beta). Furthermore, we report that human B lymphocytes have a constitutive IL-8 mRNA expression and protein secretion in vitro. Resting as well as stimulated B lymphocytes secrete on average 1.5 ng IL-8/ml medium/24 h (2x10(6) B lymphocytes). Our data indicate a possible mechanism by which B lymphocytes make contact with other cells, during immuno-inflammatory processes.

Recombinant human growth-regulated oncogene-alpha induces T lymphocyte chemotaxis. A process regulated via IL-8 receptors by IFN-gamma, TNF-alpha, IL-4, IL-10, and IL-13.

The human cytokine growth-regulated oncogene (GRO)-alpha is a small glycoprotein secreted by monocytes, endothelial cells, glycoprotein secreted by monocytes, endothelial cells, fibroblasts, synovial cells, and some tumor cells such as melanoma cells. It is structurally related to IL-8 and can activate neutrophils, whereas it induces chemotaxis, exocytosis, and a respiratory burst in neutrophils. To date, its functions on T lymphocytes have not been well established. We report here that recombinant human (rh)GRO-alpha is a potent chemoattractant for freshly isolated T lymphocytes, but not for anti-CD3 mAb-activated T lymphocytes. It attracts CD4+ and CD8+ T lymphocyte subsets to an equal extent. The migrating T lymphocytes toward rhGRO-alpha are predominantly CD45RO+ memory CD4+ and CD8+ subsets. The chemotactic migration of T lymphocytes toward rhGRO-alpha is stimulated via the IL-8 receptors on the cells. This process can be augmented by IFN-gamma and TNF-alpha, and inhibited by IL-4, IL-10, and IL-13. In addition, we also document that on T lymphocytes there exist IL-8 receptors that can be up-regulated by IFN-gamma, TNF-alpha, and IL-2. Our results demonstrate that rhGRO-alpha gene encodes for an inflammatory mediator that stimulates the directional migration of T lymphocytes. It may thus be another important mediator in the diseases in which T lymphocytes form the major constituent of the cellular infiltration.

In vitro genetically aberrant T-cell clones with continuous growth are associated with atopic dermatitis.

Atopic dermatitis is a disease with a genetic predisposition affecting the immune system, with T lymphocytes participating in the immune dysregulation. Most in vitro T lymphocyte studies of atopic dermatitis have focused on antigen-specific T-cell clones. However, antigen-non-specific regulatory T lymphocytes may also take part in the pathway leading to antigen-specific clonal T-lymphocyte proliferation. T lymphocytes from skin biopsy specimens from three patients with severe atopic dermatitis were cultured in the presence of IL-2 and IL-4, but without antigen added. Initially, proliferation was oligo- or polyclonal, but in all cases overgrowth by T cells with clonal chromosomal aberrations was subsequently observed. These abnormal T-cell clones demonstrated continuous growth and complete or partial phenotypic loss of the T-cell antigen receptor complex. In summary, these findings suggest that a subset of aberrant skin-homing T lymphocytes is associated with atopic dermatitis.

Tissue-dependent variation in the expression of elongation factor-1 alpha isoforms: isolation and characterisation of a cDNA encoding a novel variant of human elongation-factor 1 alpha.

A novel isoform of human elongation factor-1 alpha (EF-1 alpha 2) has been characterised. It shows a high similarity to other EF-1 alpha proteins, especially to a rat EF-1 alpha variant and it has all the characteristics of a functional EF-1 alpha protein. The pattern of expression of both EF-1 alpha 2 and EF-1 alpha was analysed in different human tissues. This showed that the two proteins were differentially expressed, EF-1 alpha 2 was expressed in brain, heart, skeletal muscle and in the transformed cell lines AMA and K14, but was undetectable in other tissues and in both primary and transformed human fibroblasts. EF-1 alpha was expressed in brain, placenta, lung, liver, kidney, pancreas and in all the cell lines that we have analysed but barely detectable in heart and skeletal muscle.