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.

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.

POLYMORPHISM REGENERATION FOR A NEUTRALIZED SELFISH B CHROMOSOME

Long-run evolution of B chromosomes is mainly made up by an evolutionary arms race between these selfish genetic elements and the standard genome. The suppression of B drive is one of the clearest expressions of genome defense against B chromosomes. After drive neutralization, the B is condemned to extinction unless a new variant showing drive can emerge and replace it. This paper reports the first empirical evidence for the substitution of a neutralized B variant by a new selfish B variant. Such a polymorphism regeneration has recently taken place in a natural population of the grasshopper Eyprepocnemis plorans.

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.

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.

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.

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.

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.

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.

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.

Histone H2AX phosphorylation is associated with most meiotic events in grasshopper

It is widely accepted that the H2AX histone in its phosphorylated form (gamma-H2AX) is related to the repair of DNA double-strand breaks (DSBs). In several organisms, gamma-H2AX presence has been demonstrated in meiotic processes such as recombination and sex chromosome inactivation during prophase I (from leptotene to pachytene). To test whether gamma-H2AX is present beyond pachytene, we have analysed the complete sequence of changes in H2AX phosphorylation during meiosis in grasshopper, a model organism for meiotic studies at the cytological level. We show the presence of phosphorylated H2AX during most of meiosis, with the exception only of diplotene and the end of each meiotic division. During the first meiotic division, gamma-H2AX is associated with i) recombination, as deduced from its presence in leptotene-zygotene over all chromosome length, ii) X chromosome inactivation, since at pachytene gamma-H2AX is present in the X chromosome only, and iii) chromosome segregation, as deduced from gamma-H2AX presence in centromere regions at first metaphase-anaphase. During second meiotic division, gamma-H2AX was very abundant at most chromosome lengths from metaphase to telophase, suggesting its possible association with the maintenance of chromosome condensation and segregation.

Development of a SCAR marker for the analysis of B chromosome presence in Partamona helleri (Hymenoptera, Apidae)

Chromosomes in hymenopteran insects cannot currently be analysed in adult individuals. The only available cytogenetic techniques need to be performed in larvae. Here we develop and implement a SCAR (Sequence Characterized Amplified Region) marker, associated with B chromosomes in the bee Partamona helleri, which has proven to be very useful to reveal B chromosome presence in adults from natural populations. The marker was tested in ten different colonies simultaneously analysed by both molecular (ten adults per colony) and cytogenetic (20 larvae per colony) techniques. The presence of the SCAR marker always showed the same pattern as B chromosome presence: both were present or absent in all individuals from a same colony, or both were present in only part of the individuals from a same colony. This molecular marker is thus a useful tool for analysing new aspects of this B chromosome system such as B frequency and geographical distribution, B transmission, or B effects in adult individuals.

Causes of B chromosome variant substitution in the grasshopper Eyprepocnemis plorans

We have analysed B chromosome frequency for three consecutive years, B transmission rate at population and individual levels, clutch size, egg fertility and embryo-adult viability in a natural population of the grasshopper Eyprepocnemis plorans containing two different B chromosome variants, i.e. B(2) and B(24), the second being derived from the first and having replaced it in nearby populations. From 2002 to 2003 the relative frequency of both variants changed, although the differences did not reach significance. A mother-offspring analysis showed no significant effect of any of the two B variants on clutch size, egg fertility or embryo-adult viability, but B(24) was more efficiently transmitted than B(2) through males from the 2002 season, which explains the observed frequency change. Controlled crosses, at individual level, showed significant drive through some females for B(24) but not for B(2), suggesting that this difference in transmission rate might also be important for the substitution process. The analysis of relative fitness for B(2) and B(24) carriers for all fitness components, as a whole, showed a significantly better performance of B(24)-carrying individuals, suggesting that the cumulative effect of these slight differences might contribute to the replacement of B(2) by B(24).

Detection of B chromosomes in interphase hemolymph nuclei from living specimens of the grasshopper Eyprepocnemis plorans

The two most important evolutionary properties of B chromosomes are their transmission rate (which suggests their selfishness when significantly higher than 0.5) and their net effects on carrier fitness (usually negative for parasitic Bs). The study of transmission rate unavoidably requires the analysis of many controlled crosses in order to accurately measure population average transmission rate. Therefore, getting a marker closely associated to B presence is of crucial importance to alleviate the load of performing many useless crosses between lacking B individuals. After investigating several cytogenetic techniques on several tissues that may be sampled without drastically damaging live specimens of the grasshopper Eyprepocnemis plorans, we report here the excellent results provided by the CMA3 fluorescence and C-banding techniques applied to hemocyte nuclei. These cells may be easily obtained from both males and females and provide information on B presence even during the interphase stage. The two cytogenetic techniques take advantage of the heterochromatic nature of the B chromosomes, so that Bs made predominantly of ribosomal DNA are revealed by CMA3 as bright bodies in the interphase hemocytes, and Bs mostly made of satellite DNA are visualized by C-banding as intensely stained bodies in these cells.

Are the dot-like chromosomes in Trinomys iheringi (Rodentia, Echimyidae) B chromosomes?

In this article we review the existing cytogenetic information on the polymorphic dot-like chromosomes in Trinomys iheringi, the only species in the family Echimyidae harboring them, and provide new data on the frequency, banding properties, meiotic behavior and DNA composition of these minute chromosomes. Since no individuals lacking these chromosomes have hitherto been found, one of the main properties of B chromosomes, i.e. dispensability, has not yet been tested, so that some reasonable doubt might exist on whether they are true B chromosomes. The dot-like chromosomes were also present in the twelve new individuals analyzed, showed intraindividual variation in number, most likely due to mitotic instability during development, failed to show C-bands, showed late-replication, paired among them in meiosis, but not with the large chromosomes, and appeared to be mainly composed of telomeric DNA. These results suggest that these dot-like chromosomes might actually be mitotically unstable micro B chromosomes showing very high frequency in the natural populations thus far analyzed. But, to be confident of this conclusion, individuals lacking the dot-like chromosomes should actively be searched in future research to test their dispensability.

The odd-even effect in mitotically unstable B chromosomes in grasshoppers

The odd-even effect, by which B chromosomes are more detrimental in odd numbers, has been reported in plants and animals. In grasshoppers, there are only a few reports of this effect and all were referred to as traits related to the formation of aberrant meiotic products (AMPs). Here we review the existing information about B chromosome effects on AMPs, chiasma frequency and the number of active nucleolus organizer regions (NORs) per cell. Polysomy for A chromosomes and B chromosomes are two kinds of chromosome polymorphism frequently found in grasshoppers. In some aspects, e.g. meiotic behaviour and mitotic instability leading to individual mosaicism (in the case of mitotically unstable Bs), polysomic As show similar characteristics to B chromosomes. In fact, polysomy is regarded as one of the main mechanisms for B chromosome origin. Here we review some features of meiotic behaviour in known cases of polysomy and mitotically unstable Bs in grasshoppers, in looking for possible causes for the odd-even effect. In all these traits, the odd-even effect was apparent, although its appearance was not universal in any case, with variation among species or populations within the same species. The equational division and lagging of the extra chromosomes, when univalents, could favour the appearance of abnormal meiotic products, and the formation of bivalents, when there are two or more extra chromosomes, inhibits this process. Therefore, the odd-even effect might be a consequence of the concomitant operation of both aspects of extra chromosome meiotic behaviour. The possibility that the odd-even effect might result from an increase in cell stress generated by odd numbers is suggested.

Rapid suppression of drive for a parasitic B chromosome

The persistence of parasitic B chromosomes in natural populations depends on both B ability to drive and host response to counteracting it. In the grasshopper Eyprepocnemis plorans, the B24 chromosome is the most widespread B chromosome variant in the Torrox area (Málaga, Spain). Its evolutionary success, replacing its ancestral neutralized B variant, B2, was based on meiotic drive in females, as we showed in a sample caught in 1992. In females collected six years later, mean B24 transmission ratio (k(B)) was 0.523, implying a very rapid decrease from the 0.696 observed in 1992. This shows that B24 neutralization is running very fast and suggests that it might most likely be based on a single gene of major effect.

Spatio-temporal dynamics of a neutralized B chromosome in the grasshopper Eyprepocnemis plorans

Spatial and temporal patterns of frequency variation for a neutralized B chromosome in the grasshopper Eyprepocnemis plorans were analyzed along six transects in the east of Spain to explore possible factors affecting the population dynamics of this polymorphism. Three parameters were employed to quantify B frequency: prevalence, load and mean frequency. Of them, load seemed to be the less sensitive parameter, probably due to its small range of variation. Prevalence, however, shows ample variation, but the mean frequency of B chromosomes per individual is the best parameter to characterize B frequency. Only river transects revealed significant differences among populations, and the use of two geographic explicit approaches (Mantel test and distograms) revealed significant isolation by distance (IBD), especially at the Segura River mouth, presumably due to low gene flow and drift. No temporal trend was found in the Segura River transects, which is consistent with the slow changes in B frequency expected during the random walk for neutralized B chromosomes. But these transects showed a clear spatial pattern, with B1 showing lower frequency in the upper course of this river. The present results provide the first empirical evidence of IBD in the evolution of a neutralized B chromosome, and support the notion that B dynamics at this evolutionary stage is best explained by a metapopulation approach.

The B chromosome polymorphism of the grasshopper Eyprepocnemis plorans in North Africa. IV. Transmission of rare B chromosome variants

In addition to the principal B chromosome (B(1)) in Moroccan populations of the grasshopper Eyprepocnemis plorans, nine B chromosome variants appeared at low frequency. The transmission of five of these rare B chromosome variants through females was analysed in three natural populations. Sixteen controlled crosses provided useful information on the transmission of B(M2), B(M6) and B(M7) in Smir, B(M3) and B(M6) in SO.DE.A. (Société de Développement Agricole lands near Ksar-el-Kebir city), and B(M2) and B(M10) in Mechra, all located in Morocco. Since six female parents carried two different B variants, a total of 22 progeny analyses could be studied. Intraindividual variation in B transmission rate (k(B)) was observed among the successive egg pods in 26.7 % of the females, but this variation did not show a consistent temporal pattern. Only the B(M2) and B(M6) variants in Smir showed net drive, although variation was high among crosses, especially for B(M2). These two variants are thus good candidates for future regenerations (the replacement of a neutralized B, B(1) in this case, by a new driving variant, B(M2) or B(M6)) in Smir, the northern population where the B polymorphism is presumably older. The analysis of all crosses performed in the three populations, including those reported previously for the analysis of B(1) transmission, showed that the largest variance in k(B) among crosses stands at the individual level, and not at population or type of B levels. The implications of these findings for the occurrence of possible regeneration processes in Moroccan populations are discussed.

Population variation in the A chromosome distribution of satellite DNA and ribosomal DNA in the grasshopper Eyprepocnemis plorans

The double FISH analysis of two repetitive DNAs (a satellite DNA and ribosomal DNA) in 12 natural populations of the grasshopper Eyprepocnemis plorans collected at the south (Granada and Málaga provinces) and south-east (Albacete and Murcia provinces) of the Iberian Peninsula has shown their wide-spread presence throughout the whole genome as well as extensive variation among populations. Both DNAs are found in most A chromosomes. Regularly, both DNAs occurred in the S11 and X chromosomes, rDNA in the S10 and satDNA in the L2 and M3. No correlation was found between the number of satDNA and rDNA clusters in the A genomes of the 12 populations analysed, and both figures were independent of the presence of B chromosomes. The genomic distribution of both DNAs showed no association with the geographical localization of the populations analysed. Finally, we provide evidence that the supernumerary chromosome segment proximally located on the S11 chromosome is, in most cases, the result of satDNA amplification but, in some cases, it might also derive from amplification of both satDNA and rDNA.

The B chromosome polymorphism of the grasshopper Eyprepocnemis plorans in North Africa: III. Mutation rate of B chromosomes

B chromosome variation in nine Moroccan populations of the grasshopper Eyprepocnemis plorans was analysed for 3 consecutive years. In addition to B1, which was the predominant B chromosome in all nine populations, we found 15 other B variants, albeit at very low frequency. Eight variants were found in adults caught in the wild, four appeared in adults reared in the laboratory and seven were found in embryo progeny of controlled crosses between a 0B male and a B-carrying female. Some variants were found in more than one kind of material. At least the seven B variants that appeared in embryo progeny of females carrying a different B type arose de novo through mutation of the maternal B chromosome. The mutation rate of B chromosomes was 0.73%, on average, which explains the high variety of morphs and banding patterns found. The most frequent de novo mutations observed in these chromosomes were centromere misdivision with or without chromatid nondisjunction, which generates iso-B-chromosomes or telocentric Bs, respectively, as well as translocations with A and B chromosomes and deletions. But the whole variation observed, including that found in adult individuals, suggests that other mutations such as duplications, inversions and centric fusions do usually affect B chromosomes. Finally, B chromosome mutation rate was remarkably similar in both Moroccan and Spanish populations, which suggests that it might be dependent on B chromosome intrinsic factors.

Evolutionary dynamics of a B chromosome invasion in island populations of the grasshopper Eyprepocnemis plorans

Four natural populations of the grasshopper Eyprepocnemis plorans in the Mallorca island were analysed for several years revealing the recent invasion of the B1 chromosome from the south-west part of the island (Palma region) towards the north and to the east. In only 10 years, the mean number of Bs in the northern population at Pollença increased from 0.053 to 0.692. Therefore, B chromosome invasion seems to be very rapid and has recently arrived to the north of the island. The south-west (close to Palma) is the most likely point at which B invasion started in the Mallorca Island. Finally, the number of B chromosomes was significantly associated to an increase in chiasma frequency (and thus recombination) in A chromosomes.

B-A interchanges are an unlikely pathway for B chromosome integration into the standard genome

One of the conceivable evolutionary pathways of a parasitic B chromosome is its integration into the host genome through translocation to an A chromosome. To investigate this possibility, we analyze here the nature, meiotic behavior and genetic effects of a spontaneous interchange between a medium-sized autosome and a B chromosome, found in one male caught in a Moroccan population of the grasshopper Eyprepocnemis plorans and the offspring of controlled crosses with six different females. Most metaphase I cells analyzed (115 out of 118) showed a trivalent with chiasmata in both the interstitial and pairing segments, which predicted about half of genetically unbalanced spermatozoa. The analysis of 234 embryos sired by this male on six females showed the lethality of some meiotic products paralleled to a decrease in egg fertility (0.541 +/- 0.051, compared to the 0.879 +/- 0.017 shown by females mated to standard males). These results suggest that the cost of the B-A interchange on host fitness, in terms of gametic inviability, highly diminishes the possibility of frequency increase for the interchange to reach a polymorphic status, which is the first and indispensable step to reach fixation and thus integration of the B chromosome DNA into the host genome.

Multiregional origin of B chromosomes in the grasshopper Eyprepocnemis plorans

Analysis of chromosome localization of three molecular markers, 18S-5.8S-28S rDNA, 5S rDNA and a 180 bp satDNA, showed that B chromosomes in the grasshopper Eyprepocnemis plorans originated independently in Eastern (Caucasus) and Western (Spain and Morocco) populations. Eastern B chromosomes are most likely derived from the smallest autosome, which is the only A chromosome carrying the three markers, in coincidence with Caucasian B chromosomes. Western B chromosomes, however, lack 5S rDNA and are most likely derived from the X chromosome, which is the only A chromosome carrying the two remaining markers, always in the same order with respect to the centromere, as the B chromosome.

Genetic load caused by variation in the amount of rDNA in a wasp

Extensive variation in the size of the short (heterochromatic) arm of chromosome 14 was found in the wasp Trypoxylon (Trypargilum) albitarse. Ten different variants were differentiated by size and C-banding pattern. Fluorescent in-situ hybridization (FISH) revealed that ribosomal DNA in this species is clustered in the darkly C-banded parts of the heterochromatic short arm of chromosome 14. On this basis, we got an indirect estimate of the amount of rDNA from the area of these dark C-bands. The significant absence in males of the three chromosome variants with lower amounts of rDNA indicates that these three variants are lethal in this sex, and suggests the existence of a threshold marking the minimum amount of rDNA which is tolerable in haploidy. This implies about 4% genetic load in the population caused by variation in rDNA amount.

Comparative FISH analysis in five species of Eyprepocnemidine grasshoppers

The chromosomal localization of ribosomal DNA, and a 180 bp satellite DNA isolated from Spanish Eyprepocnemis plorans specimens, has been analysed in five Eyprepocnemidinae species collected in Russia and Central Asia. Caucasian E. plorans individuals carried each of the two DNAs, but the rDNA was limited to only two chromosomes (S(9) and S(11)) in sharp contrast to Spanish specimens that show 4-8 rDNA clusters and to Moroccan specimens which carry rDNA in almost all chromosomes. The four remaining species, however, lacked the 180 bp tandem repeat, and showed rDNA clusters in one (S(9) in Thisoicetrinus pterostichus), two (S(9) and S(10) in Eyprepocnemis unicolor; M(8) and S(11) in Heteracris adspersa), or three (S(9), S(10), and S(11) in Shirakiacris shirakii) chromosome pairs. The implications of these findings for the evolution of these two chromosome markers in this group of species are discussed.

The B chromosomes of the grasshopper Eyprepocnemis plorans and the intragenomic conflict

The grasshopper Eyprepocnemis plorans harbours an extremely widespread polymorphism for supernumerary (B) chromosomes, which is found in almost all circum-Mediterranean and Caucasian populations hitherto analysed. B chromosomes in this species have been shown to evolve through several stages of parasitic and near-neutral nature, presumably because of an arms race between the standard (A) and B chromosomes. This intragenomic conflict can either be solved with the extinction of the neutralised B chromosome or, more interestingly, with the replacement of the neutralised B by a mutant version being parasitic again and thus prolonging B chromosome life. This species thus provides a complete view of the long-term life-cycle of parasitic B chromosomes.

Host recombination is dependent on the degree of parasitism

Parasites and hosts are involved in a continuous coevolutionary process leading to genetic changes in both counterparts. To understand this process, it is necessary to track host responses, one of which could be an increase in sex and recombination, such as is proposed by the Red Queen hypothesis. In this theoretical framework, the inducible recombination hypothesis states that B-chromosomes (genome parasites that prosper in natural populations of many living beings) elicit an increase in host chiasma frequency that is favoured by natural selection because it increases the proportion of recombinant progeny, some of which could be resistant to both B-chromosome effects and B-accumulation in the germline. We have found a clear parallelism between host recombination and the evolutionary status of the B-chromosome polymorphism, which provides explicit evidence for inducible recombination and strong support for the Red Queen hypothesis.

Integration of a B chromosome into the A genome of a wasp, revisited

A previous study showed that in the haplodiploid solitary wasp Trypoxylon albitarse, most individuals carry one B chromosome per haploid genome, the same dosage as the standard (A) chromosomes, indicating a possible regularization of B-chromosome meiotic behaviour and its integration into the A genome. In a new sampling, we have analysed 15 populations (including 9 out of the 10 previously analysed) to test the evolution of this integration process. The new results provide a direct report of the invasion process in the Porto Firme population, where B frequency has dramatically increased in only four generations. In the populations from the Viçosa region, however, B frequency has remained stable, although the principal B type, the metacentric one, has increased in frequency at the expense of the acrocentric one in several populations. The implications of these new results on the hypothesis of the integration of these B chromosomes, as regular members of the A genome, are discussed.

The B-chromosome polymorphism of the grasshopper Eyprepocnemis plorans in North Africa: II. Parasitic and neutralized B1 chromosomes

The transmission of the B1 chromosome through females has been analysed in three Moroccan populations (Smir, SO.DE.A. and Mechra) of the grasshopper Eyprepocnemis plorans. We analysed transmission ratio (kB) variation at two levels: intra-individual (to test female age effects) and inter-individual (to test for A chromosome effects). In 81.8% of females, kB did not differ among successive egg-pods, suggesting no effect of female age. The remaining females (18.2%), showed significant differences in kB values among egg-pods, but without clear temporal patterns. In Smir, kB ranged between B elimination (0.244) and B accumulation (0.689) but there was no net accumulation (mean +/- s.e. = 0.463 +/- 0.045). In SO.DE.A., all females analysed transmitted B1 at a Mendelian rate, with a mean kBequal to 0.512 +/- 0.020. In Mechra, kB ranged from 0.341 to 0.972, with mean kB (0.575 +/- 0.029) showing a net B accumulation. All these observations suggest that the B1 chromosome could be at a drive-suppression stage in Smir and Mechra, but that it has already been neutralised in SO.DE.A.

Transmission analysis of mitotically unstable B chromosomes in Locusta migratoria

Seventeen controlled crosses in which the mitotically unstable B chromosome of Locusta migratoria was carried by one parent only have provided evidence that B chromosomes are significantly eliminated during sexual transmission in males, at a mean rate that almost counteracts the premeiotic accumulation derived from mitotic instability during germ line development. On the other hand, B chromosomes are significantly accumulated in females, presumably by their preferential migration to the secondary oocyte during the first meiotic division. These results substantially change the current knowledge about this B chromosome system, because the main B accumulation occurs in females and not in males, as was hitherto thought. Furthermore, this case shows that the maintenance of a single B system in natural populations may be the result of many different forces and mechanisms acting for and against B chromosomes.Key words: locust, Locusta migratoria, B chromosomes, accumulation mechanisms, meiotic drive, B elimination, parthenogenesis.