Evidence for B chromosome drive suppression in the grasshopper Eyprepocnemis plorans

Evolution of B Chromosomes: From Dispensable Parasitic Chromosomes to Essential Genomic Players

B chromosomes represent additional chromosomes found in many eukaryotic organisms. Their origin is not completely understood but recent genomic studies suggest that they mostly arise through rearrangements and duplications from standard chromosomes. They can occur in single or multiple copies in a cell and are usually present only in a subset of individuals in the population. Because B chromosomes frequently show unstable inheritance, their maintenance in a population is often associated with meiotic drive or other mechanisms that increase the probability of their transmission to the next generation. For all these reasons, B chromosomes have been commonly considered to be nonessential, selfish, parasitic elements. Although it was originally believed that B chromosomes had little or no effect on an organism's biology and fitness, a growing number of studies have shown that B chromosomes can play a significant role in processes such as sex determination, pathogenicity and resistance to pathogens. In some cases, B chromosomes became an essential part of the genome, turning into new sex chromosomes or germline-restricted chromosomes with important roles in the organism's fertility. Here, we review such cases of "cellular domestication" of B chromosomes and show that B chromosomes can be important genomic players with significant evolutionary impact.

Sex-Dependent Inheritance of B Chromosomes in Psalidodon paranae (Teleostei, Characiformes) Revealed by Directed Crossings

B chromosomes are additional dispensable elements to the standard chromosomal set of an organism. In most cases, their transmission differs from Mendelian patterns, leading to their accumulation or extinction. The present study aimed to describe, for the first time, the transmission pattern of B chromosome in a population of Psalidodon paranae through directed crosses, as well as to analyze the populational dynamics of B chromosome. Our results revealed the possible elimination of B chromosome in crossings where only females were B-carriers, with a mean transmission rate (k_B) of 0.149; however, k_B was significantly higher in crossings involving male B-carriers (k_B = 0.328-0.450). Moreover, we observed an increase in the frequency of B chromosomes in the natural population of P. paranae in the last two decades. These apparently contradictory results can make sense if the B chromosome provides adaptive advantages to their carriers. Here, we observed a differential transmission of B chromosomes in each sex of parental individuals, with higher transmission rates in crossing involving males B-carriers, in addition to describe the temporal changes of B chromosome frequency in P. paranae.

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.

Long-term monitoring of B-chromosome invasion and neutralization in a population of Prospero autumnale (Asparagaceae)

B chromosomes have been reported in about 15% of eukaryotes, but long-term dynamics of B chromosomes in a single natural population has rarely been analyzed. Prospero autumnale plants collected in 1981 and 1983 at Cuesta de La Palma population had shown the presence of B chromosomes. We analyze here seven additional samples collected between 1987 and 2015, and show that B frequency increased significantly during the 1980s and showed minor fluctuations between 2005 and 2015. A mother-offspring analysis of B chromosome transmission, at population level, showed significant drive on the male side (k_B  = 0.65) and significant drag on the female side (k_B  = 0.33), with average B transmission rate being very close to the Mendelian rate (0.5). No significant effects of B chromosomes were observed on a number of vigor and fertility-related traits. Within a parasite/host framework, these results suggest that B chromosomes' drive on the male side is the main pathway for B chromosome invasion, whereas B chromosome drag on the female side might be the main manifestation of host genome resistance in this species. Prospero autumnale thus illuminates a novel evolutionary pathway for B chromosome neutralization by means of a decrease in B transmission through the nondriving sex.

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.

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.

Post-meiotic B chromosome expulsion, during spermiogenesis, in two grasshopper species

Most supernumerary (B) chromosomes are parasitic elements carrying out an evolutionary arms race with the standard (A) chromosomes. A variety of weapons for attack and defense have evolved in both contending elements, the most conspicuous being B chromosome drive and A chromosome drive suppression. Here, we show for the first time that most microspermatids formed during spermiogenesis in two grasshopper species contain expulsed B chromosomes. By using DNA probes for B-specific satellite DNAs in Eumigus monticola and Eyprepocnemis plorans, and also 18S rDNA in the latter species, we were able to count the number of B chromosomes in standard spermatids submitted to fluorescence in situ hybridization, as well as visualizing B chromosomes inside most microspermatids. In E. plorans, the presence of B-carrying microspermatids in 1B males was associated with a significant decrease in the proportion of B-carrying standard spermatids. The fact that this decrease was apparent in elongating spermatids but not in round ones demonstrates that meiosis yields 1:1 proportions of 0B and 1B spermatids and hence that B elimination takes place post-meiotically, i.e., during spermiogenesis, implying a 5-25% decrease in B transmission rate. In E. monticola, the B chromosome is mitotically unstable and B number varies between cells within a same individual. A comparison of B frequency between round and elongating spermatids of a same individual revealed a significant 12.3% decrease. We conclude that B chromosome elimination during spermiogenesis is a defense weapon of the host genome to get rid of parasitic chromosomes.

Geographical Barriers Impeded the Spread of a Parasitic Chromosome

Parasitic supernumerary (B) chromosomes show high capability to spread across populations. But the existence of abrupt discontinuities in their distribution demands an explanation. The grasshopper Eyprepocnemis plorans plorans harbour supernumerary chromosomes in all natural populations hitherto analyzed from the Circum-Mediterranean region, with the single exception of the headwaters of the Iberian Segura River and several of its tributaries. To ascertain the causes of this distribution pattern, we analyze here the genetic structure of five natural populations collected in this zone (two +B and three -B), by means of ISSR markers. We found significant population structure, with two kinds of populations coinciding with +B and -B ones, separated by strong barriers to gene flow. This gives strong support to the hypothesis that the non-B populations precede B origin, and that B-carrying individuals from coastal zones have been able to colonize upstream areas, until geographical barriers (usually narrow canyons and arid areas surrounding them) impeded their advance.

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.

Transmission rate variation among three B chromosome variants in the fish Prochilodus lineatus (Characiformes, Prochilodontidae)

Cytogenetic studies were developed in Prochilodus lineatus (Valenciennes 1836), describing an interesting system of small supernumerary chromosomes. The purpose of this work is to study the frequency and morphology of B chromosomes in individuals from the parental line and the inheritance patterns of these elements in individuals obtained from controlled crosses in the species P. lineatus. The transmission rate of B chromosomes revealed a k(B)=0.388 for the acrocentric type, a kB=0.507 for the metacentric type and a k(B)=0.526 for the submetacentric type. The obtained results raise hypothesis that B-acrocentric chromosomes are involved in an extinction process in this species, while the metacentric and submetacentric supernumerary elements comprises a neutral mechanism and follows a Mendelian transmission rate.

Population genetic structure of the grasshopper Eyprepocnemis plorans in the south and east of the Iberian Peninsula

The grasshopper Eyprepocnemis plorans subsp. plorans harbors a very widespread polymorphism for supernumerary (B) chromosomes which appear to have arisen recently. These chromosomes behave as genomic parasites because they are harmful for the individuals carrying them and show meiotic drive in the initial stages of population invasion. The rapid increase in B chromosome frequency at intrapopulation level is thus granted by meiotic drive, but its spread among populations most likely depends on interpopulation gene flow. We analyze here the population genetic structure in 10 natural populations from two regions (in the south and east) of the Iberian Peninsula. The southern populations were coastal whereas the eastern ones were inland populations located at 260–655 m altitude. The analysis of 97 ISSR markers revealed significant genetic differentiation among populations (average G_ST =  0.129), and the Structure software and AMOVA indicated a significant genetic differentiation between southern and eastern populations. There was also significant isolation by distance (IBD) between populations. Remarkably, these results were roughly similar to those found when only the markers showing low or no dropout were included, suggesting that allelic dropout had negligible effects on population genetic analysis. We conclude that high gene flow helped this parasitic B chromosome to spread through most of the geographical range of the subspecies E. plorans plorans.

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.

Population dynamics of a selfish B chromosome neutralized by the standard genome in the grasshopper Eyprepocnemis plorans

Effects of the B chromosome polymorphism of the grasshopper Eyprepocnemis plorans were analyzed in two natural populations. Postmating sexual selection, female fertility, and survival were studied. The B chromosome lacks drive and has no detectable effects on fitness. A neutral B cannot invade a population and establish a polymorphism, but the confidence limits on our estimates cannot exclude the possibility that the polymorphism is maintained by a balance between weak drive and weak selection against individuals with two and three B's. However, other lines of evidence favor the following model of the dynamics of the B in E. plorans. In a newly invaded population, the B has substantial drive, but the evolution of drive suppressor genes in the A chromosomes neutralizes the B drive so that it becomes near-neutral and begins a random walk toward extinction by stochastic loss. Because the B is common by the time drive disappears, the random walk is likely to continue for a long time. If in the course of the random walk a variant B with greater drive appears, then it will displace the original variant, and a new cycle of drive suppression and drift to extinction occurs. A simulation model of this process suggested that the mean time to extinction is proportional to the two-thirds power of the population size; it is much less affected by subpopulation size or the number of populations in a subdivided population.

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.

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.

Parasitic exploitation as an engine of diversity

Parasitic exploitation occurs within and between a wide variety of taxa in a plethora of diverse contexts. Theoretical and empirical analyses indicate that parasitic exploitation can generate substantial genetic and phenotypic polymorphism within species. Under some circumstances, parasitic exploitation may also be an important factor causing reproductive isolation and promoting speciation. Here we review research relevant to the relationship between parasitic exploitation, within species-polymorphism, and speciation in some of the major arenas in which such exploitation has been studied. This includes research on the vertebrate major histocompatibility loci, plant-pathogen interactions, the evolution of sexual reproduction, intragenomic conflict, sexual conflict, kin mimicry and social parasitism, tropical forest diversity and the evolution of language. We conclude by discussing some of the issues raised by comparing the effect of parasitic exploitation on polymorphism and speciation in different contexts.

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.

Polymorphism of attack and defense

Coevolution of attack and defense occurs in host-parasite systems and various forms of genomic conflict. Attackers benefit when their specific molecules allow entry past host defenses. Defenders gain when their matching biochemical specificities aid recognition. Selection continually favors new attack specificities that avoid matching defense and, in turn, new defense specificities that match novel attackers. The introduction of novel specificities strongly influences the spatial and temporal dynamics of conflict. Lack of reciprocally matching diversity in a particular system suggests biochemical constraints that prevent diversification. New work on cytoplasmic male sterility, B chromosomes and meiotic drive suggests that varying biochemical constraints on recognition cause varying patterns of diversity and spatiotemporal dynamics

B-chromosome evolution

B chromosomes are extra chromosomes to the standard complement that occur in many organisms. They can originate in a number of ways including derivation from autosomes and sex chromosomes in intra- and interspecies crosses. Their subsequent molecular evolution resembles that of univalent sex chromosomes, which involves gene silencing, heterochromatinization and the accumulation of repetitive DNA and transposons. B-chromosome frequencies in populations result from a balance between their transmission rates and their effects on host fitness. Their long-term evolution is considered to be the outcome of selection on the host genome to eliminate B chromosomes or suppress their effects and on the B chromosome's ability to escape through the generation of new variants. Because B chromosomes interact with the standard chromosomes, they can play an important role in genome evolution and may be useful for studying molecular evolutionary processes.

Common origin of B chromosome variants in the grasshopper eyprepocnemis plorans

Nine B chromosome variants, from seven different populations of Eyprepocnemis plorans collected at four localities in Spain and three in Morocco, have been shown to be mainly composed of two DNA sequences, i.e. a 180-bp tandem repeat and ribosomal DNA. B types, however, differ in the relative amounts of the two sequences. The most widespread one (B1) bears about the same amount of rDNA and 180 bp repeat, but three other variants that have reached a polymorphism by replacing B1 in smaller areas (B2, B5 and B24) carry a conspicuously larger amount of the 180 bp repeat. In Morocco, the most widespread B variant is also B1, and a rare variant that appeared in a single individual is also built with the same two DNA sequences. All these data point to a common origin for these B chromosomes, with B1 probably being the original one. The origin of the different B types and the possible relationship of the relative amount of 180 bp DNA repeat with B drive are discussed.