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

A New Variant B Chromosome in Auchenipteridae: The Role of (GATA)n and (TTAGGG)n Sequences in Understanding the Evolution of Supernumeraries in Trachelyopterus

Basic and molecular cytogenetic techniques were carried out in 3 Neotropical region populations of catfishes, two of Trachelyopterus galeatus (one from the marshlands of Paraguay River basin and another from Lago Catalão, Amazon River basin) and one of Trachelyopterus porosus, a sympatric population to T. galeatus from the Amazon River basin. This study aimed to describe and understand the structure and evolution of Trachelyopterus B chromosomes, mainly through physical mapping of repetitive elements. A diploid number of 58 chromosomes was found for all individuals, as well as the presence of B chromosomes. For T. porosus this is the first report of a supernumerary. The sympatric species of T. galeatus and T. porosus from Amazon River had 1-3 B chromosomes and T. galeatus from Paraguay River had 1-2 B chromosomes, all of them showed intra- and interindividual numerical variation. Two females of T. porosus exhibited a new variant B chromosome (B2), previously not seen in Auchenipteridae, which might have originated from B1 chromosomes. All B chromosomes were entirely heterochromatic. In contrast to all complement A and B2 chromosomes, in which the telomeric sequences were found in the telomeric regions, B1 chromosomes of all populations were totally marked by (TTAGGG)n probes. (GATA)n sequence sites were found through all complement A chromosomes, but B1 and B2 chromosomes exhibited only a clustered block in one of the chromosome arms. The most frequent B chromosomes (B1) in all populations/species, including those previously studied in Auchenipteridae catfishes, share the following characteristics: totally heterochromatic, small, metacentric, with accumulation of repetitive (TTAGGG)n sequences, and a low number of (GATA)n copies, which might suggest a common ancient origin in Trachelyopterus species/populations.

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.

Genome size variation affects song attractiveness in grasshoppers: evidence for sexual selection against large genomes

Genome size is largely uncorrelated to organismal complexity and adaptive scenarios. Genetic drift as well as intragenomic conflict have been put forward to explain this observation. We here study the impact of genome size on sexual attractiveness in the bow-winged grasshopper Chorthippus biguttulus. Grasshoppers show particularly large variation in genome size due to the high prevalence of supernumerary chromosomes that are considered (mildly) selfish, as evidenced by non-Mendelian inheritance and fitness costs if present in high numbers. We ranked male grasshoppers by song characteristics that are known to affect female preferences in this species and scored genome sizes of attractive and unattractive individuals from the extremes of this distribution. We find that attractive singers have significantly smaller genomes, demonstrating that genome size is reflected in male courtship songs and that females prefer songs of males with small genomes. Such a genome size dependent mate preference effectively selects against selfish genetic elements that tend to increase genome size. The data therefore provide a novel example of how sexual selection can reinforce natural selection and can act as an agent in an intragenomic arms race. Furthermore, our findings indicate an underappreciated route of how choosy females could gain indirect benefits.

B chromosomes showing active ribosomal RNA genes contribute insignificant amounts of rRNA in the grasshopper Eyprepocnemis plorans

The genetic inertness of supernumerary (B) chromosomes has recently been called into question after finding several cases of gene activity on them. The grasshopper Eyprepocnemis plorans harbors B chromosomes containing large amounts of ribosomal DNA (rDNA) units, some of which are eventually active, but the amount of rRNA transcripts contributed by B chromosomes, compared to those of the standard (A) chromosomes, is unknown. Here, we address this question by means of quantitative PCR (qPCR) for two different ITS2 amplicons, one coming from rDNA units located in both A and B chromosomes (ITS2(A+B)) and the other being specific to B chromosomes (ITS2(B)). We analyzed six body parts in nine males showing rDNA expression in their B chromosomes in the testis. Amplification of the ITS2(B) amplicon was successful in RNA extracted from all six body parts analyzed, but showed relative quantification (RQ) values four orders of magnitude lower than those obtained for the ITS(A+B) amplicon. RQ values differed significantly between body parts for the two amplicons, with testis, accessory gland and wing muscle showing threefold higher values than head, gastric cecum and hind leg. We conclude that the level of B-specific rDNA expression is extremely low even in individuals where B chromosome rDNA is not completely silenced. Bearing in mind that B chromosomes carry the largest rDNA cluster in the E. plorans genome, we also infer that the relative contribution of B chromosome rRNA genes to ribosome biogenesis is insignificant, at least in the body parts analyzed.

Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans

Wide arrays of repetitive DNA sequences form an important part of eukaryotic genomes. These repeats appear to evolve as coherent families, where repeats within a family are more similar to each other than to other orthologous representatives in related species. The continuous homogenization of repeats, through selective and non-selective processes, is termed concerted evolution. Ascertaining the level of variation between repeats is crucial to determining which evolutionary model best explains the homogenization observed for these sequences. Here, for the grasshopper Eyprepocnemis plorans, we present the analysis of intragenomic diversity for two repetitive DNA sequences (a satellite DNA (satDNA) and the 45S rDNA) resulting from the independent microdissection of several chromosomes. Our results show different homogenization patterns for these two kinds of paralogous DNA sequences, with a high between-chromosome structure for rDNA but no structure at all for the satDNA. This difference is puzzling, considering the adjacent localization of the two repetitive DNAs on paracentromeric regions in most chromosomes. The disparate homogenization patterns detected for these two repetitive DNA sequences suggest that several processes participate in the concerted evolution in E. plorans, and that these mechanisms might not work as genome-wide processes but rather as sequence-specific ones.

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.

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.

Human supernumeraries: are they B chromosomes?

In humans, the presence of supernumerary chromosomes is an unusual phenomenon, which is often associated with developmental abnormalities and malformations. In contrast to most animal and plant species, the extensive knowledge of the human genome and the ample set of molecular and cytogenetic tools available have permitted to ascertain not only that most human supernumerary chromosomes (HSCs) derive from the A chromosome set, but also the specific A chromosome from which most of them arose. These extra chromosomes are classified into six types on the basis of morphology and size. There are both heterochromatic and euchromatic HSCs, the latter being more detrimental. Most are mitotically stable, except some producing individual mosaicism. No information is available on the HSC transmission rate since extensive familial studies are not usually performed generally because of death of the relatives or lack of cooperation. The main B chromosome property failing in HSCs seems to be their population spread as polymorphisms, since most HSCs seem to correspond to extra A chromosomes or centric fragments spontaneously arisen in the analysed individual or one of his/her parents. However, we cannot rule out at this moment, that more intensive studies on population distribution and frequency of those HSCs most closely resembling B chromosomes (i.e. those heterochromatic and thus less detrimental) would reveal possible HSCs polymorphisms. Although HSCs cannot be considered B chromosomes, some of them might be a source for future B chromosomes. The best candidates would be heterochromatic HSCs, which might manage to drive in either sex. To ascertain this possibility, research on inheritance and population studies would be very helpful in combination with the powerful cytogenetic and molecular tools available for our species.

FISH detection of ribosomal cistrons and assortment-distortion for X and B chromosomes in Dichroplus pratensis (Acrididae)

Assortment-distortion with respect to the X and NOR activity of a rare mitotically stable B chromosome (B(N)), was examined in 16 males of Dichroplus pratensis (Acrididae: Melanoplinae) from Argentine populations. In 1B individuals, the X and B associate preferentially during prophase I reaching a maximum level of association at zygotene. Frequency of X/B association remains relatively high up to diplotene-diakinesis and decreases steeply towards metaphase I. The percent X/B association at each stage is positively influenced by association at the previous stage, and interindividual variability in X/B association decreases as the frequency of association increases. Both chromosomes tended to preferentially orientate toward the same pole at MI (mean ratio of 16 individuals, 1.50:1) which determined an excess of XB and 00 second spermatocytes over X0 and 0B ones (1.39:1). No significant differences occurred between the MI, AI and MII assortment ratios. Fluorescent in situ hybridisation (FISH) confirmed that the B chromosome carries ribosomal genes and helped to establish that, during spermiogenesis, both the B and the normal NOR-bearing chromosome (S8) are clustered near the centriole adjunct region of spermatids. However, FISH failed to reveal the existence of inactive ribosomal cistrons in the X chromosome, as previously suggested, thus providing no support to a simple origin of the B from the X.

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.