Román Vilas

Identification of the Major Sex-Determining Region of Turbot ( Scophthalmus maximus )

Sex determination in fish is a labile character in evolutionary terms. The sex-determining (SD) master gene can differ even between closely related fish species. This group is an interesting model for studying the evolution of the SD region and the gonadal differentiation pathway. The turbot (Scophthalmus maximus) is a flatfish of great commercial value, where a strong sexual dimorphism exists for growth rate. Following a QTL and marker association approach in five families and a natural population, we identified the main SD region of turbot at the proximal end of linkage group (LG) 5, close to the SmaUSC-E30 marker. The refined map of this region suggested that this marker would be 2.6 cM and 1.4 Mb from the putative SD gene. This region appeared mostly undifferentiated between males and females, and no relevant recombination frequency differences were detected between sexes. Comparative genomics of LG5 marker sequences against five model species showed no similarity of this chromosome to the sex chromosomes of medaka, stickleback, and fugu, but suggested a similarity to a sex-associated QTL from Oreochromis spp. The segregation analysis of the closest markers to the SD region demonstrated a ZW/ZZ model of sex determination in turbot. A small proportion of families did not fit perfectly with this model, which suggests that other minor genetic and/or environmental factors are involved in sex determination in this species.

Generalized lattice graphs for 2D-visualization of biological information

Abstract Several graph representations have been introduced for different data in theoretical biology. For instance, complex networks based on Graph theory are used to represent the structure and/or dynamics of different large biological systems such as protein–protein interaction networks. In addition, Randic, Liao, Nandy, Basak, and many others developed some special types of graph-based representations. This special type of graph includes geometrical constrains to node positioning in space and adopts final geometrical shapes that resemble lattice-like patterns. Lattice networks have been used to visually depict DNA and protein sequences but they are very flexible. However, despite the proved efficacy of new lattice-like graph/networks to represent diverse systems, most works focus on only one specific type of biological data. This work proposes a generalized type of lattice and illustrates how to use it in order to represent and compare biological data from different sources. We exemplify the following cases: protein sequence; mass spectra (MS) of protein peptide mass fingerprints (PMF); molecular dynamic trajectory (MDTs) from structural studies; mRNA microarray data; single nucleotide polymorphisms (SNPs); 1D or 2D-Electrophoresis study of protein polymorphisms and protein-research patent and/or copyright information. We used data available from public sources for some examples but for other, we used experimental results reported herein for the first time. This work may break new ground for the application of Graph theory in theoretical biology and other areas of biomedical sciences.

Characterization of EST‐derived microsatellites for gene mapping and evolutionary genomics in turbot

The detection of microsatellite sequences within expressed sequence tags (ESTs) connects potential markers with specific genes, generating type I markers. We have developed and mapped by linkage analysis a set of EST-derived microsatellites in the turbot, Scophthalmus maximus. One hundred and ninety-one microsatellites were identified from 9256 turbot ESTs. Primer design was possible with 98 microsatellites. After genotyping 25 wild turbot and the parents of two reference families for linkage analysis, 43 EST-derived microsatellites were selected because they met technical and polymorphism criteria. A final set of 31 EST-derived microsatellites could be mapped to 17 linkage groups of the turbot consensus map based on 242 anonymous microsatellites. Twenty-four microsatellite-containing ESTs were functionally annotated, confirming them as type I markers. Nineteen were mapped in the turbot consensus map. These EST-derived microsatellites constitute useful tools for genome scanning of turbot populations, marker-assisted selection programmes and comparative mapping.

Genetic variability of natural populations of trematodes of the genus Lecithochirium parasites of eels.

Allozyme variation within and among populations of 3 species of the genus Lecithochirium (Trematoda: Hemiuridae) was studied by starch gel electrophoresis. In total, 19 loci were analysed in 7 populations. The level of genetic variability was relatively high in all populations. The percentage of polymorphic loci (0.95 criterion) observed per population varied from 21.0% to 55.5%, and expected heterozygosity levels varied from 0.082 to 0.197. All populations showed significant heterozygote deficiencies. In Lecithochirium fusiforme most of the deviations from Hardy-Weinberg proportions were within the populations and this species showed moderate population structuring (F(IS)=0.486, F(ST)=0.142, Nm= 1.51) and accordingly low intraspecific genetic distances (D=0.003 to 0.027). A significant lack of heterozygotes for several polymorphic loci was revealed in Lecithochirium rufoviride and Lecithochirium musculus. The most probable cause of the population genetic subdivision in L. rufoviride is the presence of at least 1 cryptic species in the populations studied. Although the lowest percentage of fixed genetic differences was that between L. fusiforme and L. musculus, two different algorithms for the construction of evolutionary trees on a matrix of genetic distances confirmed that L. fusiforme and L. rufoviride are phenetically the most closely related species.

More than meets the eye: detecting cryptic microgeographic population structure in a parasite with a complex life cycle

Nonrandom recruitment of parasites among hosts can lead to genetic differentiation among hosts and mating dynamics that promote inbreeding. It has been hypothesized that strictly aquatic parasites with intermediate hosts will behave as panmictic populations among hosts because ample opportunity exists for random mixing of unrelated individuals during transmission to the definitive host. A previous allozyme study on the marine trematode Lecithochirium fusiforme did not support this hypothesis; in that, there was genetic differentiation among, and significant heterozygote deficiencies within, definitive hosts. We revisit this system and use microsatellites to obtain multilocus genotypes. Our goal was to determine whether cryptic subgroups and/or the presence of clones could account for the apparent deviation from ‘panmixia’. We find strong evidence for cryptic subdivision (three genetic clusters) that causes the Wahlund effect and differentiation among definitive hosts. After accounting for these cryptic groups, we see panmictic genetic structure among definitive hosts that is consistent with the ‘high mixing in aquatic habitats’ hypothesis. We see evidence for cotransmission of clones in all three clusters, but this level of clonal structure did not have a major impact in causing deviations from Hardy–Weinberg equilibrium, and only affected genetic differentiation among hosts in one cluster. A cursory examination of the data may have led to incorrect conclusions about nonrandom transmission. However, it is obvious in this system that there is more than meets the eye in relation to the actual make‐up of parasite populations. In general, the methods we employ will be useful for elucidating hidden patterns in other organisms where cryptic structure may be common (e.g. those with limited morphology or complex life histories).

Development and Validation of Single Nucleotide Polymorphisms (SNPs) Markers from Two Transcriptome 454-Runs of Turbot (Scophthalmus maximus) Using High-Throughput Genotyping

The turbot (Scophthalmus maximus) is a commercially valuable flatfish and one of the most promising aquaculture species in Europe. Two transcriptome 454-pyrosequencing runs were used in order to detect Single Nucleotide Polymorphisms (SNPs) in genes related to immune response and gonad differentiation. A total of 866 true SNPs were detected in 140 different contigs representing 262,093 bp as a whole. Only one true SNP was analyzed in each contig. One hundred and thirteen SNPs out of the 140 analyzed were feasible (genotyped), while III were polymorphic in a wild population. Transition/transversion ratio (1.354) was similar to that observed in other fish studies. Unbiased gene diversity (He) estimates ranged from 0.060 to 0.510 (mean = 0.351), minimum allele frequency (MAF) from 0.030 to 0.500 (mean = 0.259) and all loci were in Hardy-Weinberg equilibrium after Bonferroni correction. A large number of SNPs (49) were located in the coding region, 33 representing synonymous and 16 non-synonymous changes. Most SNP-containing genes were related to immune response and gonad differentiation processes, and could be candidates for functional changes leading to phenotypic changes. These markers will be useful for population screening to look for adaptive variation in wild and domestic turbot.

Genetic variation within and among infrapopulations of the marine digenetic trematode Lecithochirium fusiforme.

Allozyme markers were used to study genetic variation in Lecithochirium fusiforme within a natural population of Conger conger. Six of 16 enzyme-coding loci studied were found to be polymorphic. These loci were surveyed in 12 infrapopulations of adult flukes. High levels of genetic variation were detected (P = 0.375); Ho = 0.048; He = 0.085). However, the population did not conform to Hardy-Weinberg expectations, as it showed a significant deficit of heterozygotes. L. fusiforme also exhibited low differentiation between infrapopulations (FST = 0.064). Despite significant linkage disequilibrium at Pgm-1 and Pgm-2 (P < 0.05), mating system does not appear to be the principal reason for the deficit of heterozygotes detected, because some polymorphic loci were in Hardy-Weinberg equilibrium. Association between FIS and FST statistics suggests the existence of the Wahlund effect. However, all infrapopulations showed a strong deficit of heterozygotes for most polymorphic loci (FIS = 0.409). Detection of significant genetic differentiation among temporal samples and the existence of paratenic hosts in the life-cycle suggests the Wahlund effect, caused by the mixture of genetically distinct temporal samples in the infrapopulations. Occasional temporal gene flow also might explain the high estimated genetic polymorphism.

Regional environmental pressure influences population differentiation in turbot (Scophthalmus maximus)

Unravelling the factors shaping the genetic structure of mobile marine species is challenging due to the high potential for gene flow. However, genetic inference can be greatly enhanced by increasing the genomic, geographical or environmental resolution of population genetic studies. Here, we investigated the population structure of turbot (Scophthalmus maximus) by screening 17 random and gene‐linked markers in 999 individuals at 290 geographical locations throughout the northeast Atlantic Ocean. A seascape genetics approach with the inclusion of high‐resolution oceanographical data was used to quantify the association of genetic variation with spatial, temporal and environmental parameters. Neutral loci identified three subgroups: an Atlantic group, a Baltic Sea group and one on the Irish Shelf. The inclusion of loci putatively under selection suggested an additional break in the North Sea, subdividing southern from northern Atlantic individuals. Environmental and spatial seascape variables correlated marginally with neutral genetic variation, but explained significant proportions (respectively, 8.7% and 10.3%) of adaptive genetic variation. Environmental variables associated with outlier allele frequencies included salinity, temperature, bottom shear stress, dissolved oxygen concentration and depth of the pycnocline. Furthermore, levels of explained adaptive genetic variation differed markedly between basins (3% vs. 12% in the North and Baltic Sea, respectively). We suggest that stable environmental selection pressure contributes to relatively strong local adaptation in the Baltic Sea. Our seascape genetic approach using a large number of sampling locations and associated oceanographical data proved useful for the identification of population units as the basis of management decisions.

Contrasting patterns of population genetic structure of Fasciola hepatica from cattle and sheep: Implications for the evolution of anthelmintic resistance

Twelve polymorphic genetic markers, eight allozymic loci and four microsatellites, were used to characterize 20 infrapopulations of Fasciola hepatica (all flukes from 10 individual cattle and 10 sheep) from 11 farms in Northwest Spain. Results suggest different patterns of population genetic structure depending on the host species. Individuals identified as clones were much more frequent in sheep. The common presence of clones and its nonrandom occurrence among individual hosts suggests clumped transmission of liver flukes in sheep. After reducing significant repeated multilocus genotypes to one unique copy within infrapopulations, results show relatively high levels of gene diversity within infrapopulations from cattle and sheep (0.411 and 0.360 on average, respectively). However, parasites of sheep appear to show significantly more structured variation at the infrapopulation level (Standardized F(ST)=0.087 and 0.170 for parasites of cattle and sheep, respectively). Compared to the parasites from cattle, results suggest that populations from sheep show lower levels of gene flow, higher degree of aggregate transmission, higher probability of mating within clones, and lower parasitic load. These differences have implications for the evolution of anthelmintic resistance because they affect the effective population size and the degree of inbreeding. The development and rapid spread of resistance seems likely in the parasites of cattle because populations from the study area are characterized by high gene flow. However, results also suggest that the efficient selection of a new recessive advantageous mutation would be favored in parasites of sheep due to a greater potential for inbreeding.

Very low microsatellite polymorphism and large heterozygote deficits suggest founder effects and cryptic structure in the parasite Perkinsus olseni

Twelve microsatellite markers were used to characterize 130 clonal cultures of Perkinsus olseni derived from 30 clams from six different geographic locations. Only two loci were polymorphic in the four populations studied from Spanish coast (mean sample size = 31.2), and a third locus was variable in only two populations. In contrast, five parasites isolated from five clams from Japan and New Zealand showed variation at nine loci. Low genetic variation (2.08 ± 0.64 alleles per locus; mean genetic diversity: 0.101 ± 0.022), and very high F(IS) values (0.857 on average) were observed in Spanish populations. A total of 39 multilocus genotypes (MLGs) were identified in the whole sample (121 clonal isolates after excluding incomplete MLGs due to missing data). A three-level hierarchical analysis of molecular variance found significant levels of genetic variation within infrapopulations (all the parasites in a single host; Φ(IS) = 0.679) and among infrapopulations within the component population (all the parasites among a host population; Φ(SC) = 0.579). Differences among the component population from different geographic locations were not significant (Φ(CT) = 0.057). These results suggest that an important fraction of F(IS) is explained by the Wahlund effect, but also strong inbreeding within infrapopulations. Another explanation for the high F(IS) within infrapopulations is the presence of haploid and diploid stages in the clam. Although fully aquatic system provides many opportunities for mixing of parasites from different clams, results are consistent with the consideration of all P. olseni in a clam as a cohesive genetic unit (i.e., deme). If the parasite was introduced into the Spanish coast with the importation of infected clams from Asia and Oceania, the low microsatellite polymorphism could be reflecting founder effects in the recent evolutionary history of P. olseni. The loss of alleles would be intensified in a scenario structured in numerous demes because of recurrent founder effects at microgeographic level.