Rapid suppression of drive for a parasitic B chromosome

Interpopulation spread of a parasitic B chromosome is unlikely through males in the grasshopper Eyprepocnemis plorans

The near-neutral model of B chromosome evolution predicts that population invasion is quite fast. To test this prediction, in 1994, we introduced males of the grasshopper Eyprepocnemis plorans from a B-carrying population into a B-lacking population and monitored the evolution of B-chromosome frequency up to 2013. We observed fluctuating very low B frequency across years but, remarkably, the B chromosome introduced (the B2 variant) was found up to 1996 only, whereas the B1 variant was present from 1996 onwards, presumably introduced by fishermen using E. plorans males as bait. Effective introgression of genetic material from the donor population was evidenced by the presence of a satellite DNA on autosome 9 (up to 1999) and the presence of one individual in 2006 showing an ISSR marker profile being highly similar to that found in the donor population. This indicated that the males introduced by us effectively mated with resident females, but donor genes rapidly decreased in frequency after this non-recurrent migration event. Taken together, our results indicated: (i) that the non-recurrent migration event had a slight, transient genetic effect on the recipient population, which was diluted in only a few generations; and (ii) that even with recurrent migration (forced by fishermen) the B chromosome failed to increase in frequency. Bearing in mind that B chromosomes in this species drive through females only, we hypothesize that B chromosomes most likely failed invasion in both migration events because the migrating sex shows no B-drive.

Protein-coding genes in B chromosomes of the grasshopper Eyprepocnemis plorans

For many years, parasitic B chromosomes have been considered genetically inert elements. Here we show the presence of ten protein-coding genes in the B chromosome of the grasshopper Eyprepocnemis plorans. Four of these genes (CIP2A, GTPB6, KIF20A, and MTG1) were complete in the B chromosome whereas the six remaining (CKAP2, CAP-G, HYI, MYCB2, SLIT and TOP2A) were truncated. Five of these genes (CIP2A, CKAP2, CAP-G, KIF20A, and MYCB2) were significantly up-regulated in B-carrying individuals, as expected if they were actively transcribed from the B chromosome. This conclusion is supported by three truncated genes (CKAP2, CAP-G and MYCB2) which showed up-regulation only in the regions being present in the B chromosome. Our results indicate that B chromosomes are not so silenced as was hitherto believed. Interestingly, the five active genes in the B chromosome code for functions related with cell division, which is the main arena where B chromosome destiny is played. This suggests that B chromosome evolutionary success can lie on its gene content.

B-chromosome effects on Hsp70 gene expression does not occur at transcriptional level in the grasshopper Eyprepocnemis plorans

As intragenomic parasites, B chromosomes can elicit stress in the host genome, thus inducing a response for host adaptation to this kind of continuous parasitism. In the grasshopper Eyprepocnemis plorans, B-chromosome presence has been previously associated with a decrease in the amount of the heat-shock protein 70 (HSP70). To investigate whether this effect is already apparent at transcriptional level, we analyze the expression levels of the Hsp70 gene in gonads and somatic tissues of males and females with and without B chromosomes from two populations, where the predominant B chromosome variants (B2 and B24) exhibit different levels of parasitism, by means of quantitative real-time PCR (qPCR) on complementary DNA (cDNA). The results revealed the absence of significant differences for Hsp70 transcripts associated with B-chromosome presence in virtually all samples. This indicates that the decrease in HSP70 protein levels, formerly reported in this species, may not be a consequence of transcriptional down-regulation of Hsp70 genes, but the result of post-transcriptional regulation. These results will help to design future studies oriented to identifying factors modulating Hsp70 expression, and will also contribute to uncover the biological role of B chromosomes in eukaryotic genomes.

Preferential occupancy of R2 retroelements on the B chromosomes of the grasshopper Eyprepocnemis plorans

R2 non-LTR retrotransposons exclusively insert into the 28S rRNA genes of their host, and are expressed by co-transcription with the rDNA unit. The grasshopper Eyprepocnemis plorans contains transcribed rDNA clusters on most of its A chromosomes, as well as non-transcribed rDNA clusters on the parasitic B chromosomes found in many populations. Here the structure of the E. plorans R2 element, its abundance relative to the number of rDNA units and its retrotransposition activity were determined. Animals screened from five populations contained on average over 12,000 rDNA units on their A chromosomes, but surprisingly only about 100 R2 elements. Monitoring the patterns of R2 insertions in individuals from these populations revealed only low levels of retrotransposition. The low rates of R2 insertion observed in E. plorans differ from the high levels of R2 insertion previously observed in insect species that have many fewer rDNA units. It is proposed that high levels of R2 are strongly selected against in E. plorans, because the rDNA transcription machinery in this species is unable to differentiate between R2-inserted and uninserted units. The B chromosomes of E. plorans contain an additional 7,000 to 15,000 rDNA units, but in contrast to the A chromosomes, from 150 to over 1,500 R2 elements. The higher concentration of R2 in the inactive B chromosomes rDNA clusters suggests these chromosomes can act as a sink for R2 insertions thus further reducing the level of insertions on the A chromosomes. These studies suggest an interesting evolutionary relationship between the parasitic B chromosomes and R2 elements.

Spread of a new parasitic B chromosome variant is facilitated by high gene flow

The B₂₄ chromosome variant emerged several decades ago in a Spanish population of the grasshopper Eyprepocnemis plorans and is currently reaching adjacent populations. Here we report, for the first time, how a parasitic B chromosome (a strictly vertically transmitted parasite) expands its geographical range aided by high gene flow in the host species. For six years we analyzed B frequency in several populations to the east and west of the original population and found extensive spatial variation, but only a slight temporal trend. The highest B₂₄ frequency was found in its original population (Torrox) and it decreased closer to both the eastern and the western populations. The analysis of Inter Simple Sequence Repeat (ISSR) markers showed the existence of a low but significant degree of population subdivision, as well as significant isolation by distance (IBD). Pairwise N_em estimates suggested the existence of high gene flow between the four populations located in the Torrox area, with higher values towards the east. No significant barriers to gene flow were found among these four populations, and we conclude that high gene flow is facilitating B₂₄ diffusion both eastward and westward, with minor role for B₂₄ drive due to the arrival of drive suppressor genes which are also frequent in the donor population.

B-chromosome ribosomal DNA is functional in the grasshopper Eyprepocnemis plorans

B-chromosomes are frequently argued to be genetically inert elements, but activity for some particular genes has been reported, especially for ribosomal RNA (rRNA) genes whose expression can easily be detected at the cytological level by the visualization of their phenotypic expression, i.e., the nucleolus. The B(24) chromosome in the grasshopper Eyprepocnemis plorans frequently shows a nucleolus attached to it during meiotic prophase I. Here we show the presence of rRNA transcripts that unequivocally came from the B(24) chromosome. To detect these transcripts, we designed primers specifically anchoring at the ITS-2 region, so that the reverse primer was complementary to the B chromosome DNA sequence including a differential adenine insertion being absent in the ITS2 of A chromosomes. PCR analysis carried out on genomic DNA showed amplification in B-carrying males but not in B-lacking ones. PCR analyses performed on complementary DNA showed amplification in about half of B-carrying males. Joint cytological and molecular analysis performed on 34 B-carrying males showed a close correspondence between the presence of B-specific transcripts and of nucleoli attached to the B chromosome. In addition, the molecular analysis revealed activity of the B chromosome rDNA in 10 out of the 13 B-carrying females analysed. Our results suggest that the nucleoli attached to B chromosomes are actively formed by expression of the rDNA carried by them, and not by recruitment of nucleolar materials formed in A chromosome nucleolar organizing regions. Therefore, B-chromosome rDNA in E. plorans is functional since it is actively transcribed to form the nucleolus attached to the B chromosome. This demonstrates that some heterochromatic B chromosomes can harbour functional genes.

A single, recent origin of the accessory B chromosome of the grasshopper Eyprepocnemis plorans

B chromosomes are dispensable chromosomes found in >2000 eukaryotic species, usually behaving as genomic parasites. Most B chromosomes seem to be made up of the same kind of DNA sequences present in the A chromosomes. This sequence similarity makes it difficult to obtain specific molecular probes that may permit B-presence diagnosis without cytogenetic analysis. We have developed a sequence-characterized amplified region (SCAR) marker for B chromosomes in the grasshopper Eyprepocnemis plorans, which specifically amplifies a 1510-bp DNA fragment exclusively in B-carrying individuals. Fluorescent in situ hybridization and fiber FISH analyses showed that this marker is a tandemly repeated DNA sequence closely intermingled with 45S rDNA. PCR reactions showed the presence of SCAR-like sequences in the A chromosomes, but in two separate fragments, supporting the intraspecific origin of B chromosomes in this species. SCAR marker DNA sequence showed to be identical in B chromosome variants from several localities from Spain and Morocco, and it was very similar to those found in B chromosome variants from Greece and Armenia. This strongly suggests that this sequence was already present in the ancestral B chromosome of this species. In addition, the scarce sequence variation observed among several B variants from very distant populations suggests either a functional constraint or, more likely, a recent and unique origin for B chromosomes in this species.

A- and B-chromosome pairing and recombination in male meiosis of the silver fox (Vulpes vulpes L., 1758, Carnivora, Canidae)

We examined A- and B-chromosome pairing and recombination in 12 males from the farm-bred population of the silver fox (2n = 34 + 0-10 Bs) by means of electron and immunofluorescent microscopy. To detect recombination at A and B chromosomes, we used immunolocalisation of MLH1, a mismatch repair protein of mature recombination nodules, at synaptonemal complexes. The mean total number of MLH1 foci at A-autosomes was 29.6 foci per cell. The XY bivalent had one MLH1 focus at the pairing region. Total recombination length of the male fox genome map was estimated as 1,530 centimorgans. We detected single MLH1 foci at 61% of linear synaptic configurations involving B chromosomes. The distribution of the foci along B- and A-bivalents was the same. This may be considered as a first molecular evidence that meiotic recombination does occur in mammalian B chromosomes. There was no correlation between the number of synaptic configurations involving B chromosomes per cell and the recombination rate of the A-genome.

Quantitative analysis of NOR expression in a B chromosome of the grasshopper Eyprepocnemis plorans

The B24 chromosome in the Torrox population of the grasshopper Eyprepocnemis plorans is recurrently attached to a nucleolus in diplotene cells, indicating the activity of its distally located ribosomal DNA (rDNA). The frequency of males expressing the B chromosome nucleolus organizer region (B-NOR) almost doubled in 4 years. The likelihood of expressing the B-NOR increased with the B number and, in males expressing it, about 20% of their cells showed a nucleolus attached to the B. When active, the B-NOR contributed more than 25% of total cell nucleolar area (NA). Within males expressing the B-NOR, total cell NA did not differ between cells showing the active or inactive B-NOR, suggesting that total cell NA is tightly regulated in this species. However, this parameter tended to increase in this population from 1999 to 2004, in parallel to the neutralization process which is taking place in this population. Finally, an analysis of A chromosome NOR interdependence for activity revealed a positive correlation among autosomes but a negative correlation between autosomes and the X chromosome, the manifestation of which depends on B-NOR activity. These results are discussed in the context of the nucleolus as a sensor of the stress caused by parasitic B chromosomes.

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.