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.

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.

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.

42Sp48 in previtellogenic Xenopus oocytes is structurally homologous to EF-1 alpha and may be a stage-specific elongation factor.

We have isolated the cDNA for 42Sp48 and EF-1 alpha from mixed stage oocytes and tailbud (stage 22) Xenopus laevis cDNA libraries by use of the cDNA for human elongation factor-1 alpha (EF-1 alpha) as probe. The nucleotide and deduced amino acid sequences of the entire coding region of 42Sp48 and EF-1 alpha cDNA were established. The proposed functional homology of the proteins is reflected in highly conserved amino acid sequences (91% identity), while the large number of silent mutations at the gene level may serve to prevent recombination at their loci. 42Sp48 is apparently encoded by two genes in Xenopus, while no sequences corresponding to 42Sp48 could be found in murine or human genomic DNA. 42Sp48 has been proposed to act as a stage-specific elongation factor in Xenopus. Comparison of the deduced amino acid sequences of 42Sp48 and EF-1 alpha with that of elongation factor Tu from E. coli, for which the three-dimensional structure including that of the GTP binding sites have been determined, supports this hypothesis.

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.

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.

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.

Isolation and characterization of the gene encoding EF-1 alpha O, an elongation factor 1-alpha expressed during early development of Xenopus laevis.

In Xenopus laevis, the gene encoding the elongation factor 1-alpha variant EF-1 alpha O, where O stands for oocyte, is expressed in oocytes and early embryos. A genomic library from X. laevis was screened with a cDNA probe coding for EF-1 alpha O. Two recombinant phages were isolated, one of which carries an entire EF-1 alpha O gene. This clone was characterized by restriction enzyme mapping and sequencing. Comparison of cDNA and genomic sequences revealed that EF-1 alpha O consists of seven exons spanning about 6.5 kb. The structure of the gene is very homologous to the human EF-1 alpha gene, as all locations of the splice junctions are conserved between the two genes. The sequence immediately upstream from the transcription start point (tsp) contains a CCAAT box, but does not contain either a TATA box or a Sp1-binding site. Interestingly, this sequence has a sequence homologous to the negative regulatory element from the TFIIIA promoter. A region located about 400 bp upstream from the tsp contains an additional number of possible regulatory sequence elements. The first intron contains G + C-rich elements which exist both isolated and as part of longer inverted repeats. Furthermore, one octamer and four Sp1-binding sites are found in this intron.

Identification of nuclear factors which interact with the 5' flanking region of the EF-1 alpha O gene in Xenopus laevis.

The EF-1 alpha O gene of Xenopus laevis is a stage-specific gene, being transcribed in oogonia and oocytes, but not in postmeiotic germ cells and terminally differentiated cells. We found that two trans-acting factors from oocyte nuclear extract are able to interact with a DNA sequence in the 5'-upstream region of the EF-1 alpha O gene. Methylation interference experiments suggested that the two factors recognised the same DNA element. Gel retardation assays indicated that part of the protein binding site could be confined to a 21 bp sequence, located between -51 and -72, relative to the cap site. Interestingly, this region shares great homology to a negative regulatory segment in the promoter of the TFIIIA gene, another developmentally regulated gene.

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.

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.

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.

The sequence of 16S rRNA from Mycoplasma strain PG50.

The DNA sequence of one of the 16S rRNA genes (cistron rrnA) of Mycoplasma strain PG50 was determined. It is 1523 bp long and has about 70% homology to the sequence of Escherichia coli 16S rRNA (rrnB). The G + C content of the sequence is 48% compared with 56% G + C in the E. coli sequence. The secondary structure is formed and it is determined that most of the differences between the two sequences are seen in stems while loops in the secondary structure are conserved. A detailed description of differences and similarities to known sequences and rRNA oligonucleotide catalogues is given, and this information is used to discuss functional properties and phylogenetic relations of mycoplasma 16S rRNA.

Cloning of Mycoplasma pneumoniae DNA and expression of P1-epitopes in Escherichia coli.

A genomic library of Mycoplasma pneumoniae was constructed by cloning partial Sau3A-digested genomic DNA into the expression plasmids pEX1 to pEX3. The recombinant clones were screened for production of M. pneumoniae P1-antigen by an in situ colony enzyme-linked immunosorbent assay (ELISA) blot method with a monospecific rabbit antiserum raised against the surface protein P1. The length of the translated P1-sequence and the size of the inserted DNA were determined. By comparison it was shown that six clones contained DNA fragments coding for an internal part of the P1-protein and eight clones code for the C-terminal part of terminated P1-protein. In reactions against sera from patients suffering from M. pneumoniae infection and sera from healthy persons, one of the internal clones and five of the C-terminal clones reacted with one or two of the patient sera, but only one clone reacted with all patient sera.

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.

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.