Full plastome sequence of the fern Vandenboschia speciosa (Hymenophyllales): structural singularities and evolutionary insights.

We provide here the first full chloroplast genome sequence, i.e., the plastome, for a species belonging to the fern order Hymenophyllales. The phylogenetic position of this order within leptosporangiate ferns, together with the general scarcity of information about fern plastomes, places this research as a valuable study on the analysis of the diversity of plastomes throughout fern evolution. Gene content of V. speciosa plastome was similar to that in most ferns, although there were some characteristic gene losses and lineage-specific differences. In addition, an important number of genes required U to C RNA editing for proper protein translation and two genes showed start codons alternative to the canonical AUG (AUA). Concerning gene order, V. speciosa shared the specific 30-kb inversion of euphyllophytes plastomes and the 3.3-kb inversion of fern plastomes, keeping the ancestral gene order shared by eusporangiate and early leptosporangiate ferns. Conversely, V. speciosa has expanded IR regions comprising the rps7, rps12, ndhB and trnL genes in addition to rRNA and other tRNA genes, a condition shared with several eusporangiate ferns, lycophytes and hornworts, as well as most seed plants.

Characterization of the satellitome in lower vascular plants: the case of the endangered fern Vandenboschia speciosa.

BACKGROUND AND AIMS: Vandenboschia speciosa is a highly vulnerable fern species, with a large genome (10.5 Gb). Haploid gametophytes and diploid sporophytes are perennial, can reproduce vegetatively, and certain populations are composed only of independent gametophytes. These features make this fern a good model: (1) for high-throughput analysis of satellite DNA (satDNA) to investigate possible evolutionary trends in satDNA sequence features; (2) to determine the relative contribution of satDNA and other repetitive DNAs to its large genome; and (3) to analyse whether the reproduction mode or phase alternation between long-lasting haploid and diploid stages influences satDNA abundance or divergence. METHODS: We analysed the repetitive fraction of the genome of this species in three different populations (one comprised only of independent gametophytes) using Illumina sequencing and bioinformatic analysis with RepeatExplorer and satMiner. KEY RESULTS: The satellitome of V. speciosa is composed of 11 satDNA families, most of them showing a short repeat length and being A + T rich. Some satDNAs had complex repeats composed of sub-repeats, showing high similarity to shorter satDNAs. Three families had particular structural features and highly conserved motifs. SatDNA only amounts to approx. 0.4 % of its genome. Likewise, microsatellites do not represent more than 2 %, but transposable elements (TEs) represent approx. 50 % of the sporophytic genomes. We found high resemblance in satDNA abundance and divergence between both gametophyte and sporophyte samples from the same population and between populations. CONCLUSIONS: (1) Longer (and older) satellites in V. speciosa have a higher A + T content and evolve from shorter ones and, in some cases, microsatellites were a source of new satDNAs; (2) the satellitome does not explain the huge genome size in this species while TEs are the major repetitive component of the V. speciosa genome and mostly contribute to its large genome; and (3) reproduction mode or phase alternation between gametophytes and sporophytes does not entail accumulation or divergence of satellites.

Non-random expression of ribosomal DNA units in a grasshopper showing high intragenomic variation for the ITS2 region.

We analyse intragenomic variation of the ITS2 internal transcribed spacer of ribosomal DNA (rDNA) in the grasshopper Eyprepocnemis plorans, by means of tagged PCR 454 amplicon sequencing performed on both genomic DNA (gDNA) and RNA-derived complementary DNA (cDNA), using part of the ITS2 flanking coding regions (5.8S and 28S rDNA) as an internal control for sequencing errors. Six different ITS2 haplotypes (i.e. variants for at least one nucleotide in the complete ITS2 sequence) were found in a single population, one of them (Hap4) being specific to a supernumerary (B) chromosome. The analysis of both gDNA and cDNA from the same individuals provided an estimate of the expression efficiency of the different haplotypes. We found random expression (i.e. about similar recovery in gDNA and cDNA) for three haplotypes (Hap1, Hap2 and Hap5), but significant underexpression for three others (Hap3, Hap4 and Hap6). Hap4 was the most extremely underexpressed and, remarkably, it showed the lowest sequence conservation for the flanking 5.8-28S coding regions in the gDNA reads but the highest conservation (100%) in the cDNA ones, suggesting the preferential expression of mutation-free rDNA units carrying this ITS2 haplotype. These results indicate that the ITS2 region of rDNA is far from complete homogenization in this species, and that the different rDNA units are not expressed at random, with some of them being severely downregulated.

A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs.

The desert locust (Schistocerca gregaria) has been used as material for numerous cytogenetic studies. Its genome size is estimated to be 8.55 Gb of DNA comprised in 11 autosomes and the X chromosome. Its X0/XX sex chromosome determinism therefore results in females having 24 chromosomes whereas males have 23. Surprisingly, little is known about the DNA content of this locust's huge chromosomes. Here, we use the Feulgen Image Analysis Densitometry and C-banding techniques to respectively estimate the DNA quantity and heterochromatin content of each chromosome. We also identify three satellite DNAs using both restriction endonucleases and next-generation sequencing. We then use fluorescent in situ hybridization to determine the chromosomal location of these satellite DNAs as well as that of six tandem repeat DNA gene families. The combination of the results obtained in this work allows distinguishing between the different chromosomes not only by size, but also by the kind of repetitive DNAs that they contain. The recent publication of the draft genome of the migratory locust (Locusta migratoria), the largest animal genome hitherto sequenced, invites for sequencing even larger genomes. S. gregaria is a pest that causes high economic losses. It is thus among the primary candidates for genome sequencing. But this species genome is about 50 % larger than that of L. migratoria, and although next-generation sequencing currently allows sequencing large genomes, sequencing it would mean a greater challenge. The chromosome sizes and markers provided here should not only help planning the sequencing project and guide the assembly but would also facilitate assigning assembled linkage groups to actual chromosomes.

U1 snDNA clusters in grasshoppers: chromosomal dynamics and genomic organization.

The spliceosome, constituted by a protein set associated with small nuclear RNA (snRNA), is responsible for mRNA maturation through intron removal. Among snRNA genes, U1 is generally a conserved repetitive sequence. To unveil the chromosomal/genomic dynamics of this multigene family in grasshoppers, we mapped U1 genes by fluorescence in situ hybridization in 70 species belonging to the families Proscopiidae, Pyrgomorphidae, Ommexechidae, Romaleidae and Acrididae. Evident clusters were observed in all species, indicating that, at least, some U1 repeats are tandemly arrayed. High conservation was observed in the first four families, with most species carrying a single U1 cluster, frequently located in the third or fourth longest autosome. By contrast, extensive variation was observed among Acrididae, from a single chromosome pair carrying U1 to all chromosome pairs carrying it, with occasional occurrence of two or more clusters in the same chromosome. DNA sequence analysis in Eyprepocnemis plorans (species carrying U1 clusters on seven different chromosome pairs) and Locusta migratoria (carrying U1 in a single chromosome pair) supported the coexistence of functional and pseudogenic lineages. One of these pseudogenic lineages was truncated in the same nucleotide position in both species, suggesting that it was present in a common ancestor to both species. At least in E. plorans, this U1 snDNA pseudogenic lineage was associated with 5S rDNA and short interspersed elements (SINE)-like mobile elements. Given that we conclude in grasshoppers that the U1 snDNA had evolved under the birth-and-death model and that its intragenomic spread might be related with mobile elements.

B1 was the ancestor B chromosome variant in the western Mediterranean area in the grasshopper Eyprepocnemis plorans.

We analyzed the distribution of 2 repetitive DNAs, i.e. ribosomal DNA (rDNA) and a satellite DNA (satDNA), on the B chromosomes found in 17 natural populations of the grasshopper Eyprepocnemis ploransplorans sampled around the western Mediterranean region, including the Iberian Peninsula, Balearic Islands, Sicily, and Tunisia. Based on the amount of these repetitive DNAs, 4 types of B variants were found: B1, showing an equal or higher amount of rDNA than satDNA, and 3 other variants, B2, B24 and B5, bearing a higher amount of satDNA than rDNA. The variants B1 and B2 varied in size among populations: B1 was about half the size of the X chromosome in Balearic Islands, but two-thirds of the X in Iberian populations at Alicante, Murcia and Albacete provinces. Likewise, B2 was about one-third the size of the X chromosome in populations from the Granada province but half the size of the X in the populations collected at Málaga province. The widespread geographical distribution of the B1 variant makes it the best candidate for being the ancestor B chromosome in the whole western Mediterranean region.

Common descent of B chromosomes in two species of the fish genus Prochilodus (Characiformes, Prochilodontidae).

To ascertain the origin of B chromosomes in 2 fish species of the genus Prochilodus, i.e. P. lineatus and P. nigricans, we microdissected them and generated B-specific DNA probes. These probes were used to perform chromosome painting in both species and in 3 further ones belonging to the same genus (P. argenteus, P. brevis and P. costatus). Both probes hybridized with the B chromosomes in P. lineatus and P. nigricans, but with none of the chromosomes in the 5 species. This indicates that the B chromosomes have low similarity with DNAs located in the A chromosomes and suggests the possibility that the B chromosomes in the 2 species have a common origin. The most parsimonious explanation would imply intergeneric hybridization in an ancestor of P. lineatus and P. nigricans yielding the B chromosome as a byproduct, which remained in these 2 species after their phylogenetic origin, but was perhaps lost in other Prochilodus species. This hypothesis predicts that B chromosomes are old genomic elements in this genus, and this could be tested once a species from a relative genus would be found showing homology of its A chromosomes with the B-probes employed here, through a comparison of B chromosome DNA sequences with those in the A chromosomes of this other species.

Chromosomal localization of ribosomal and telomeric DNA provides new insights on the evolution of gomphocerinae grasshoppers.

Chromosome location of ribosomal DNA (rDNA) and telomeric repeats was analysed in mitotic chromosomes of 15 species of Gomphocerinae grasshoppers belonging to the tribes Arcypterini, Gomphocerini, Stenobothrini, and Chrysochraontini. Two types of rDNA distribution were found in the Gomphocerini tribe. Type 1, found in 9 species, was characterized by the presence of rDNA in the short arm of the long biarmed chromosomes 2 and 3 and, in some species, also in the X chromosome. Type 2 was found only in Aeropus sibiricus and Stauroderus scalaris and consisted in the presence of pericentromeric rDNA blocks in all chromosomes. A comparison of rDNA distribution in Gomphocerini species with 2n ♂ = 23, 2n ♂ = 21, and 2n ♂ = 17 suggested the possible involvement of chromosome 6 in the ancestral karyotype (2n ♂ = 23) in 1 of the 3 centric fusions that decreased the chromosome number in these species. In the tribe Stenobothrini, Stenobothrus eurasius carried a single rDNA cluster in the X chromosome, likewise 2 Spanish species previously analysed, but Omocestus viridulus unusually showed a single rDNA cluster in the longest autosome. Telomeric repeats were located primarily on the ends of chromosome arms. In 2 species, however, we observed the presence of interstitial clusters outside telomeric regions. The first one, Aeropus sibiricus, exhibited a polymorphic interstitial site of telomeric repeats in chromosome 6 as a consequence of a paracentric inversion. Most remarkably, Chorthippus jacobsoni showed the presence of telomeric repeats in the pericentric regions of the 3 biarmed chromosome pairs originated by centric fusion, thus suggesting that these rearrangements were not of the Robertsonian type but true centric fusion with a probable generation of dicentric chromosomes.

DNA amount of X and B chromosomes in the grasshoppers Eyprepocnemis plorans and Locusta migratoria.

We analyzed the DNA amount in X and B chromosomes of 2 XX/X0 grasshopper species (Eyprepocnemis plorans and Locusta migratoria), by means of Feulgen image analysis densitometry (FIAD), using previous estimates in L. migratoria as standard (5.89 pg). We first analyzed spermatids of 0B males and found a bimodal distribution of integrated optical densities (IODs), suggesting that one peak corresponded to +X and the other to -X spermatids. The difference between the 2 peaks corresponded to the X chromosome DNA amount, which was 1.28 pg in E. plorans and 0.80 pg in L. migratoria. In addition, the +X peak in E. plorans gave an estimate of the C-value in this species (10.39 pg). We next analyzed diplotene cells from 1B males in E. plorans and +B males in L. migratoria (a species where Bs are mitotically unstable and no integer B number can be defined for an individual) and measured B chromosome IOD relative to X chromosome IOD, within the same cell, taking advantage of the similar degree of condensation for both positively heteropycnotic chromosomes at this meiotic stage. From this proportion, we estimated the DNA amount for 3 different B chromosome variants found in individuals from 3 E. plorans Spanish populations (0.54 pg for B1 from Saladares, 0.51 pg for B2 from Salobreña and 0.64 for B24 from Torrox). Likewise, we estimated the DNA amount of the B chromosome in L. migratoria to be 0.15 pg. To automate measurements, we wrote a GPL3 licensed Python program (pyFIA). We discuss the utility of the present approach for estimating X and B chromosome DNA amount in a variety of situations, and the meaning of the DNA amount estimates for X and B chromosomes in these 2 species.

B chromosomes and sex in animals.

Supernumerary (B) chromosomes are dispensable elements found in many eukaryote genomes in addition to standard (A) chromosomes. In many respects, B chromosomes resemble sex chromosomes, so that a common ancestry for them has frequently been suggested. For instance, B chromosomes in grasshoppers, and other insects, show a pycnotic cycle of condensation-decondensation during meiosis remarkably similar to that of the X chromosome. In some cases, B chromosome size is even very similar to that of the X chromosome. These resemblances have led to suggest the X as the B ancestor in many cases. In addition, sex chromosome origin from B chromosomes has also been suggested. In this article, we review the existing evidence for both evolutionary pathways, as well as sex differences for B frequency at adult and embryo progeny levels, B chromosome effects or B chromosome transmission. In addition, we review cases found in the literature showing sex-ratio distortion associated with B chromosome presence, the most extreme case being the paternal sex ratio (PSR) chromosomes in some Hymenoptera. We finally analyse the possibility of B chromosome regularisation within the host genome and, as a consequence of it, whether B chromosomes can become regular members of the host genome.

Level of heat shock proteins decreases in individuals carrying B-chromosomes in the grasshopper Eyprepocnemis plorans.

We analyzed the effect of B-chromosome presence on expression level of heat shock protein 70 (Hsp70) in cerebral ganglion and gonad in both males and females of the grasshopper Eyprepocnemis plorans. Two natural Spanish populations, Salobreña (Granada) and Torrox (Málaga) were assayed, the former harbouring a neutralized (non-driving) B-chromosome (B(2)) and the latter a parasitic (driving) B-chromosome (B(24)). The analysis was performed by Western blotting, immunostaining and densitometric measuring expression level of the Hsp70 family in adult individuals. The results showed that Hsp70 levels of testis were significantly higher in Salobreña than Torrox, and were significantly lower in testes of B-carrying males from both populations. A similar effect was observed in the ovary of females from Torrox. No effect was, however, observed in cerebral ganglia in any sex or population. B-chromosome effects in Torrox showed a dose-dependent pattern. The results point to an interesting interaction between B-chromosome and stress protein expression in reproductive tissue.

Repetitive DNAs and differentiation of sex chromosomes in neotropical fishes.

The processes working on sex chromosome differentiation are still not completely understood. However, the accumulation of repetitive DNA sequences has been shown to be one of the first steps in the early stages of such differentiation. In addition, regions with suppressed or no recombination have a potential to accumulate these DNA sequences and, for this reason, the absence of recombination between the sex chromosomes favors, by itself, the accumulation of repetitive sequences on these chromosomes during evolution. The diversity of sex-determining mechanisms in fish, alongside with the absence of heteromorphic sex chromosomes in many species, makes this group a useful model to better understand evolutionary processes of sex chromosomes in vertebrates, considering that fish occupy the basal position in the phylogeny of this group. In this review we draw attention to a preferential accumulation and enrichment in repetitive DNAs in sex chromosomes of many neotropical fish species in comparison with autosomes. This phenomenon has been observed between both morphologically differentiated and nascent sex chromosome systems, which highlight the potential role of these sequences in the differentiation of fish sex chromosomes generating differences in morphology and size between them.

Abnormal spermatid formation in the presence of the parasitic B(24) chromosome in the grasshopper Eyprepocnemis plorans.

Morphology and size of spermatids were analysed in the grasshopper Eyprepocnemis plorans by means of light and electron microscopy. At light microscopy, normal and abnormal (macro- and micro-) spermatids differed in size and number of centriolar adjuncts (CAs): 1 CA in normal spermatids and 2 or more CAs, depending on ploidy level, in macrospermatids. Males carrying the additional B(24) chromosome showed significantly more macro- and microspermatids than 0B males. The frequency of macro- and microspermatids showed an odd-even pattern in respect to the number of B chromosomes, with a higher frequency of abnormal spermatids associated with odd B numbers. Transmission electron microscopy showed that macrospermatids carried more than one axoneme, depending on ploidy level: 2 for diploid, 3 for triploid, and 4 for tetraploid spermatids. In 0B males, the most frequent abnormal spermatids were diploid, whereas in 1B males they were the tetraploid spermatids and, to a lesser extent, triploid ones. This suggests that most macrospermatids derived from cytokinesis failure and nucleus restitution. The implications of aberrant spermatids on B chromosome transmission and male fertility are discussed.

Microdissection and chromosome painting of X and B chromosomes in the grasshopper Eyprepocnemis plorans.

The relative location of 2 repetitive DNAs, i.e. ribosomal (rDNA) and a tandemly repeated satellite DNA (satDNA), with respect to the centromere, suggested that B chromosomes in the grasshopper Eyprepocnemis plorans derived intraspecifically from the X chromosome. To test this hypothesis, we microdissected X and B chromosomes and amplified the obtained DNA by 2 different procedures, the conventional DOP-PCR method and the single-cell whole-genome amplification GenomePlex method. We then generated DNA probes to perform chromosome painting. Our results have confirmed that X and B chromosomes share many DNA sequences between them and with most of the autosomes, especially at locations where the satDNA and rDNA reside, in consistency with previous information. This supports the hypothesis of an intraspecific origin of B chromosomes in E. plorans. Nevertheless, the present results did not help to clarify whether Bs were derived from the X chromosome or else from 1 or more autosomes.

Local adaptation and maladaptation to pollinators in a generalist geographic mosaic.

The Geographic Mosaic Theory of Coevolution predicts the occurrence of mosaics of interaction-mediated local adaptations and maladaptations. Empirical support to this prediction has come mostly from specialist interactions. In contrast, local adaptation is considered highly unlikely in generalist interactions. In this study, we experimentally test local adaptation in a generalist plant-pollinator geographic mosaic, by means of a transplant experiment in which plants coming from two evolutionary hotspots and two coldspots were offered to pollinators at the same four localities. Plants produced in the hotspots attracted more pollinators in all populations, whereas coldspot plants attracted fewer pollinators in all populations. Differences in adaptation were not related to genetic similarity between populations, suggesting that it was mainly due to spatial variation in previous selective regimes. Our experiment provides the first strong support for a spatially structured pattern of adaptation and maladaptation generated by a generalist free-living mutualism.

Differences in ribosomal DNA distribution on A and B chromosomes between eastern and western populations of the grasshopper Eyprepocnemis plorans plorans.

Distribution of ribosomal DNA (rDNA) on standard (A) and supernumerary (B) chromosomes of the grasshopper Eyprepocnemis plorans was analysed in specimens collected in Turkey and Armenia, belonging to the E. p. plorans subspecies, and in South Africa, belonging to the E. p. meridionalis subspecies. The latter individuals showed rDNA loci in chromosomes 9 and 11 only, whereas those from Armenia carried it in chromosomes 9 and 11 or else in chromosomes 9-11, depending on the population. The specimens from Turkey carried it in chromosomes 1, 9-11 and X. A comparison of this pattern with those previously observed in populations from Spain, Morocco, and Greece (belonging to E. p. plorans) suggests the existence of two evolutionary patterns in rDNA chromosome location in A chromosomes of this subspecies: eastern populations showing rDNA restricted to the small (9-11) chromosomes (as in E. p. meridionalis and other closely related taxa within the Eyprepocneminae subfamily) and western populations carrying rDNA in most A chromosomes (Spain) or all of them (Morocco). The intermediate pattern discerned in geographically intermediate populations (in Greece and Turkey), with rDNA also being located on the X chromosome, suggests a possible east-west cline. Additional support for east-west differentiation in the rDNA location pattern comes from the analysis of B chromosomes. In eastern populations, including Daghestan, Armenia, Turkey, and Greece, B chromosomes are composed mostly of rDNA, whereas in western populations (Spain and Morocco) they contain roughly similar amounts of rDNA and a 180-bp tandem repeat (satDNA), the latter being scarce in eastern Bs.

Spatial variation in selection on corolla shape in a generalist plant is promoted by the preference patterns of its local pollinators.

An adaptive role of corolla shape has been often asserted without an empirical demonstration of how natural selection acts on this trait. In generalist plants, in which flowers are visited by diverse pollinator fauna that commonly vary spatially, detecting pollinator-mediated selection on corolla shape is even more difficult. In this study, we explore the mechanisms promoting selection on corolla shape in the generalist crucifer Erysimum mediohispanicum Polatschek (Brassicaceae). We found that the main pollinators of E. mediohispanicum (large bees, small bees and bee flies) discriminate between different corolla shapes when offered artificial flowers without reward. Importantly, different pollinators prefer different shapes: bees prefer flowers with narrow petals, whereas bee flies prefer flowers with rounded overlapping petals. We also found that flowers with narrow petals (those preferred by bees) produce both more pollen and nectar than those with rounded petals. Finally, different plant populations were visited by different faunas. As a result, we found spatial variation in the selection acting on corolla shape. Selection favoured flowers with narrow petals in the populations where large or small bees are the most abundant pollinator groups. Our study suggests that pollinators, by preferring flowers with high reward, exert strong selection on the E. mediohispanicum corolla shape. The geographical variation in the pollinator-mediated selection on E. mediohispanicum corolla shape suggests that phenotypic evolution and diversification can occur in this complex floral trait even without specialization.

Association between floral traits and rewards in Erysimum mediohispanicum (Brassicaceae).

BACKGROUND AND AIMS: Floral rewards may be associated with certain morphological floral traits and thus act as underlying factors promoting selection on these traits. This study investigates whether some traits that are under pollinator-mediated selection (flower number, stalk height, corolla diameter, corolla tube length and corolla tube width) in the Mediterranean herb E. mediohispanicum (Brassicaceae) are associated with rewards (pollen and nectar). METHODS: During 2005 the phenotypic traits and the visitation rate of the main pollinator functional groups were quantified in 720 plants belonging to eight populations in south-east Spain, and during 2006 the same phenotypic traits and the reward production were quantified in 400 additional plants from the same populations. KEY RESULTS: A significant correlation was found between nectar production rate and corolla tube length, and between pollen production and corolla diameter. Visitation rates of large bees and butterflies were significantly higher in plants exhibiting larger flowers with longer corolla tubes. CONCLUSIONS: The association between reward production and floral traits may be a factor underlying the pattern of visitation rate displayed by some pollinators.

Histone H3 lysine 9 acetylation pattern suggests that X and B chromosomes are silenced during entire male meiosis in a grasshopper.

The facultative heterochromatic X chromosome in leptotene spermatocytes of the grasshopper Eyprepocnemis plorans showed marked hypoacetylation for lysine 9 in the H3 histone (H3-K9) with no sign of histone H2AX phosphorylation. Since H3-K9 hypoacetylation precedes the meiotic appearance of phosphorylated H2AX (gamma-H2AX), which marks the beginning of recombinational DNA double-strand breaks (DSBs), it seems that meiotic sex-chromosome inactivation (MSCI) in this grasshopper occurs prior to the beginning of recombination and hence synapsis (which in this species begins later than recombination). In addition, all constitutively heterochromatic chromosome regions harbouring a 180-bp tandem-repeat DNA and rDNA (B chromosomes and pericentromeric regions of A chromosomes) were H3-K9 hypoacetylated at early leptotene even though they will synapse at subsequent stages. This also suggests that meiotic silencing in this grasshopper might be independent of synapsis. The H3-K9 hypoacetylated state of facultative and constitutive heterochromatin persisted during subsequent meiotic stages and was even apparent in round spermatids. Finally, the fact that B chromosomes are differentially hypoacetylated in testis and embryo interphase cells suggests that they might be silenced early in development and remain this way for most (or all) life-cycle stages.

Physical mapping of rDNA and satDNA in A and B chromosomes of the grasshopper Eyprepocnemis plorans from a Greek population.

Adult males and females of the grasshopper Eyprepocnemis plorans from a Greek population were analysed by C-banding, silver impregnation and double FISH for two DNA probes, i.e. ribosomal DNA (rDNA) and a 180-bp tandem repeat DNA (satDNA). This population shows characteristics of rDNA location in A chromosomes that are intermediate between those previously reported for eastern (Caucasus) and western (Spain and Morocco) populations. The four rDNA clusters revealed by FISH in chromosomes X, 9, 10 and 11 in Greek specimens imply two more than the two observed in chromosomes 9 and 11 in the Caucasus, but less than the 12 observed in all chromosomes in Morocco. Remarkably, the X chromosome bears one of the new rDNA locations in Greece with respect to the Caucasus, but it appears to be inactive, in contrast to X chromosomes in western populations, which are usually active. B chromosomes were very frequent in the Greek population, and three variants differing in size were observed, all of these being largely composed of rDNA, with the exception of a small pericentromeric satDNA cluster. The high B frequency suggests that B chromosomes in this population might behave parasitically, in resemblance to Bs in western populations.