Satellitome comparison of two oedipodine grasshoppers highlights the contingent nature of satellite DNA evolution

BACKGROUND

The full catalog of satellite DNA (satDNA) within a same genome constitutes the satellitome. The Library Hypothesis predicts that satDNA in relative species reflects that in their common ancestor, but the evolutionary mechanisms and pathways of satDNA evolution have never been analyzed for full satellitomes. We compare here the satellitomes of two Oedipodine grasshoppers (Locusta migratoria and Oedaleus decorus) which shared their most recent common ancestor about 22.8 Ma ago.

RESULTS

We found that about one third of their satDNA families (near 60 in every species) showed sequence homology and were grouped into 12 orthologous superfamilies. The turnover rate of consensus sequences was extremely variable among the 20 orthologous family pairs analyzed in both species. The satDNAs shared by both species showed poor association with sequence signatures and motives frequently argued as functional, except for short inverted repeats allowing short dyad symmetries and non-B DNA conformations. Orthologous satDNAs frequently showed different FISH patterns at both intra- and interspecific levels. We defined indices of homogenization and degeneration and quantified the level of incomplete library sorting between species.

CONCLUSIONS

Our analyses revealed that satDNA degenerates through point mutation and homogenizes through partial turnovers caused by massive tandem duplications (the so-called satDNA amplification). Remarkably, satDNA amplification increases homogenization, at intragenomic level, and diversification between species, thus constituting the basis for concerted evolution. We suggest a model of satDNA evolution by means of recursive cycles of amplification and degeneration, leading to mostly contingent evolutionary pathways where concerted evolution emerges promptly after lineages split.

Occasional paternal inheritance of the germline-restricted chromosome in songbirds

Songbirds have one special accessory chromosome, the so-called germline-restricted chromosome (GRC), which is only present in germline cells and absent from all somatic tissues. Earlier work on the zebra finch (Taeniopygia guttata castanotis) showed that the GRC is inherited only through the female line-like the mitochondria-and is eliminated from the sperm during spermatogenesis. Here, we show that the GRC has the potential to be paternally inherited. Confocal microscopy using GRC-specific fluorescent in situ hybridization probes indicated that a considerable fraction of sperm heads (1 to 19%) in zebra finch ejaculates still contained the GRC. In line with these cytogenetic data, sequencing of ejaculates revealed that individual males from two families differed strongly and consistently in the number of GRCs in their ejaculates. Examining a captive-bred male hybrid of the two zebra finch subspecies (T. g. guttata and T. g. castanotis) revealed that the mitochondria originated from a castanotis mother, whereas the GRC came from a guttata father. Moreover, analyzing GRC haplotypes across nine castanotis matrilines, estimated to have diverged for up to 250,000 y, showed surprisingly little variability among GRCs. This suggests that a single GRC haplotype has spread relatively recently across all examined matrilines. A few diagnostic GRC mutations that arose since this inferred spreading suggest that the GRC has continued to jump across matriline boundaries. Our findings raise the possibility that certain GRC haplotypes could selfishly spread through the population via occasional paternal transmission, thereby outcompeting other GRC haplotypes that were limited to strict maternal inheritance, even if this was partly detrimental to organismal fitness.

Out of patterns, the euchromatic B chromosome of the grasshopper Abracris flavolineata is not enriched in high-copy repeats

In addition to the normal set of standard (A) chromosomes, some eukaryote species harbor supernumerary (B) chromosomes. In most cases, B chromosomes show differential condensation with respect to A chromosomes and display dark C-bands of heterochromatin, and some of them are highly enriched in repetitive DNA. Here we perform a comprehensive NGS (next-generation sequencing) analysis of the repeatome in the grasshopper Abracris flavolineata aimed at uncovering the molecular composition and origin of its B chromosome. Our results have revealed that this B chromosome shows a DNA repeat content highly similar to the DNA repeat content observed for euchromatic (non-C-banded) regions of A chromosomes. Moreover, this B chromosome shows little enrichment for high-copy repeats, with only a few elements showing overabundance in B-carrying individuals compared to the 0B individuals. Consequently, the few satellite DNAs (satDNAs) mapping on the B chromosome were mostly restricted to its centromeric and telomeric regions, and they displayed much smaller bands than those observed on the A chromosomes. Our data support the intraspecific origin of the B chromosome from the longest autosome by misdivision, isochromosome formation, and additional restructuring, with accumulation of specific repeats in one or both B chromosome arms, yielding a submetacentric B. Finally, the absence of B-specific satDNAs, which are frequent in other species, along with its euchromatic nature, suggest that this B chromosome arose recently and might still be starting a heterochromatinization process. On this basis, it could be a good model to investigate the initial steps of B chromosome evolution.

Transposable element landscapes illuminate past evolutionary events in the endangered fern Vandenboschia speciosa

Vandenboschia speciosa is an endangered tetraploid fern species with a large genome (10.5 Gb). Its geographical distribution is characterized by disjoined tertiary flora refuges, with relict populations that survived past climate crises. Here, we analyzed the transposable elements (TEs) and found that they comprise approximately 76% of the V. speciosa genome, thus being the most abundant type of DNA sequence in this gigantic genome. The V. speciosa genome is composed of 51% and 5.6% of Class I and Class II elements, respectively. LTR retrotransposons were the most abundant TEs in this species (at least 42% of the genome), followed by non-LTR retrotransposons, which constituted at least 8.7% of the genome of this species. We introduce an additional analysis to identify the nature of non-annotated elements (19% of the genome). A BLAST search of the non-annotated contigs against the V. speciosa TE database allowed for the identification of almost half of them, which were most likely diverged sequence variants of the annotated TEs. In general, the TE composition in V. speciosa resembles the TE composition in seed plants. In addition, repeat landscapes revealed three episodes of amplification for all TEs, most likely due to demographic changes associated with past climate crises.

Long-term persistence of supernumerary B chromosomes in multiple species of Astyanax fish

BACKGROUND

Eukaryote genomes frequently harbor supernumerary B chromosomes in addition to the "standard" A chromosome set. B chromosomes are thought to arise as byproducts of genome rearrangements and have mostly been considered intraspecific oddities. However, their evolutionary transcendence beyond species level has remained untested.

RESULTS

Here we reveal that the large metacentric B chromosomes reported in several fish species of the genus Astyanax arose in a common ancestor at least 4 million years ago. We generated transcriptomes of A. scabripinnis and A. paranae 0B and 1B individuals and used these assemblies as a reference for mapping all gDNA and RNA libraries to quantify coverage differences between B-lacking and B-carrying genomes. We show that the B chromosomes of A. scabripinnis and A. paranae share 19 protein-coding genes, of which 14 and 11 were also present in the B chromosomes of A. bockmanni and A. fasciatus, respectively. Our search for B-specific single-nucleotide polymorphisms (SNPs) identified the presence of B-derived transcripts in B-carrying ovaries, 80% of which belonged to nobox, a gene involved in oogenesis regulation. Importantly, the B chromosome nobox paralog is expressed > 30× more than the A chromosome paralog. This indicates that the normal regulation of this gene is altered in B-carrying females, which could potentially facilitate B inheritance at higher rates than Mendelian law prediction.

CONCLUSIONS

Taken together, our results demonstrate the long-term survival of B chromosomes despite their lack of regular pairing and segregation during meiosis and that they can endure episodes of population divergence leading to species formation.

Eight Million Years of Satellite DNA Evolution in Grasshoppers of the Genus Schistocerca Illuminate the Ins and Outs of the Library Hypothesis

Satellite DNA (satDNA) is an abundant class of tandemly repeated noncoding sequences, showing high rate of change in sequence, abundance, and physical location. However, the mechanisms promoting these changes are still controversial. The library model was put forward to explain the conservation of some satDNAs for long periods, predicting that related species share a common collection of satDNAs, which mostly experience quantitative changes. Here, we tested the library model by analyzing three satDNAs in ten species of Schistocerca grasshoppers. This group represents a valuable material because it diversified during the last 7.9 Myr across the American continent from the African desert locust (Schistocerca gregaria), and this thus illuminates the direction of evolutionary changes. By combining bioinformatic and cytogenetic, we tested whether these three satDNA families found in S. gregaria are also present in nine American species, and whether differential gains and/or losses have occurred in the lineages. We found that the three satDNAs are present in all species but display remarkable interspecies differences in their abundance and sequences while being highly consistent with genus phylogeny. The number of chromosomal loci where satDNA is present was also consistent with phylogeny for two satDNA families but not for the other. Our results suggest eminently chance events for satDNA evolution. Several evolutionary trends clearly imply either massive amplifications or contractions, thus closely fitting the library model prediction that changes are mostly quantitative. Finally, we found that satDNA amplifications or contractions may influence the evolution of monomer consensus sequences and by chance playing a major role in driftlike dynamics.

Programmed DNA elimination of germline development genes in songbirds

In some eukaryotes, germline and somatic genomes differ dramatically in their composition. Here we characterise a major germline–soma dissimilarity caused by a germline-restricted chromosome (GRC) in songbirds. We show that the zebra finch GRC contains >115 genes paralogous to single-copy genes on 18 autosomes and the Z chromosome, and is enriched in genes involved in female gonad development. Many genes are likely functional, evidenced by expression in testes and ovaries at the RNA and protein level. Using comparative genomics, we show that genes have been added to the GRC over millions of years of evolution, with embryonic development genes bicc1 and trim71 dating to the ancestor of songbirds and dozens of other genes added very recently. The somatic elimination of this evolutionarily dynamic chromosome in songbirds implies a unique mechanism to minimise genetic conflict between germline and soma, relevant to antagonistic pleiotropy, an evolutionary process underlying ageing and sexual traits.

Songbirds have extensive germline–soma genome differences due to developmental elimination of a germline-specific chromosome (GRC). Here, the authors show that the GRC contains dozens of expressed developmental genes, some of which have been on the GRC since the ancestor of all songbirds.

Satellite DNA content of B chromosomes in the characid fish Characidium gomesi supports their origin from sex chromosomes

The origin of supernumerary (B) chromosomes is clearly conditioned by their ancestry from the standard (A) chromosomes. Sequence similarity between A and B chromosomes is thus crucial to determine B chromosome origin. For this purpose, we compare here the DNA sequences from A and B chromosomes in the characid fish Characidium gomesi using two main approaches. First, we found 59 satellite DNA (satDNA) families constituting the satellitome of this species and performed FISH analysis for 18 of them. This showed the presence of six satDNAs on the B chromosome: one shared with sex chromosomes and autosomes, two shared with sex chromosomes, one shared with autosomes and two being B-specific. This indicated that B chromosomes most likely arose from the sex chromosomes. Our second approach consisted of the analysis of five repetitive DNA families: 18S and 5S ribosomal DNA (rDNA), the H3 histone gene, U2 snDNA and the most abundant satDNA (CgoSat01-184) on DNA obtained from microdissected B chromosomes and from B-lacking genomes. PCR and sequence analysis of these repetitive sequences was successful for three of them (5S rDNA, H3 histone gene and CgoSat01-184), and sequence comparison revealed that DNA sequences obtained from the B chromosomes displayed higher identity with C. gomesi genomic DNA than with those obtained from other Characidium species. Taken together, our results support the intraspecific origin of B chromosomes in C. gomesi and point to sex chromosomes as B chromosome ancestors, which raises interesting prospects for future joint research on the genetic content of sex and B chromosomes in this species.

DNA metabarcoding of Neotropical ichthyoplankton: Enabling high accuracy with lower cost

Knowledge of ichthyoplankton dynamics is extremely important for conservation management as it can provide information about preferential spawning sites, reproductive period, migratory routes and recruitment success, which can be used to guide management and conservation efforts. However, identification of the eggs and larvae of Neotropical freshwater fish is a difficult task. DNA barcodes have emerged as an alternative and highly accurate approach for species identification, but DNA barcoding can be time-consuming and costly. To solve this problem, we aimed to develop a simple protocol based on DNA metabarcoding, to investigate whether it is possible to detect and quantify all species present in a pool of organisms. To do this, 230 larvae were cut in half, one half was sequenced by the Sanger technique and the other half was used to compose six arrays with a pool of larvae that were sequenced using a next-generation technique (NGS). The results of the Sanger sequencing allowed the identification of almost all larvae at species level, and the results from NGS showed high accuracy in species detection, ranging from 83% to 100%, with an average of 95% in all samples. No false positives were detected. The frequency of organisms in the two methods was positively correlated (Pearson), with low variation among species. In conclusion, this protocol represents a considerable advance in ichthyoplankton studies, allowing a rapid, cost-effective, quali-quantitative approach that improves the accuracy of identification.

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.

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.

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.

High-throughput analysis of satellite DNA in the grasshopper Pyrgomorpha conica reveals abundance of homologous and heterologous higher-order repeats

Satellite DNA (satDNA) constitutes an important fraction of repetitive DNA in eukaryotic genomes, but it is barely known in most species. The high-throughput analysis of satDNA in the grasshopper Pyrgomorpha conica revealed 87 satDNA variants grouped into 76 different families, representing 9.4% of the genome. Fluorescent in situ hybridization (FISH) analysis of the 38 most abundant satDNA families revealed four different patterns of chromosome distribution. Homology search between the 76 satDNA families showed the existence of 15 superfamilies, each including two or more families, with the most abundant superfamily representing more than 80% of all satDNA found in this species. This also revealed the presence of two types of higher-order repeats (HORs), one showing internal homologous subrepeats, as conventional HORs, and an additional type showing non-homologous internal subrepeats, the latter arising by the combination of a given satDNA family with a non-annotated sequence, or with telomeric DNA. Interestingly, the heterologous subrepeats included in these HORs showed higher divergence within the HOR than outside it, suggesting that heterologous HORs show poor homogenization, in high contrast with conventional (homologous) HORs. Finally, heterologous HORs can show high differences in divergence between their constituent subrepeats, suggesting the possibility of regional homogenization.

First complete female mitochondrial genome in four bivalve species genus Donax and their phylogenetic relationships within the Veneroida order

BACKGROUND

Four species of the genus Donax (D. semistriatus, D. trunculus, D. variegatus and D. vittatus) are common on Iberian Peninsula coasts. Nevertheless, despite their economic importance and overexploitation, scarce genetic resources are available. In this work, we newly determined the complete mitochondrial genomes of these four representatives of the family Donacidae, with the aim of contributing to unveil phylogenetic relationships within the Veneroida order, and of developing genetic markers being useful in wedge clam identification and authentication, and aquaculture stock management.

PRINCIPAL FINDINGS

The complete female mitochondrial genomes of the four species vary in size from 17,044 to 17,365 bp, and encode 13 protein-coding genes (including the atp8 gene), 2 rRNAs and 22 tRNAs, all located on the same strand. A long non-coding region was identified in each of the four Donax species between cob and cox2 genes, presumably corresponding to the Control Region. The Bayesian and Maximum Likelihood phylogenetic analysis of the Veneroida order indicate that all four species of Donax form a single clade as a sister group of other bivalves within the Tellinoidea superfamily. However, although Tellinoidea is actually monophyletic, none of its families are monophyletic.

CONCLUSIONS

Sequencing of complete mitochondrial genomes provides highly valuable information to establish the phylogenetic relationships within the Veneroida order. Furthermore, we provide here significant genetic resources for further research and conservation of this commercially important fishing resource.

Quantitative sequence characterization for repetitive DNA content in the supernumerary chromosome of the migratory locust

Repetitive DNA is a major component in most eukaryotic genomes but is ignored in most genome sequencing projects. Here, we report the quantitative composition in repetitive DNA for a supernumerary (B) chromosome, in the migratory locust (Locusta migratoria), by Illumina sequencing of genomic DNA from B-carrying and B-lacking individuals and DNA obtained from a microdissected B chromosome, as well as the physical mapping of some elements. B chromosome DNA of 94.9% was repetitive, in high contrast with the 64.1% of standard (A) chromosomes. B chromosomes are enriched in satellite DNA (satDNA) (65.2% of B-DNA), with a single satellite (LmiSat02-176) comprising 55% of the B. Six satDNAs were visualized by FISH on the B chromosome, and the only A chromosome carrying all these satellites was autosome 9, pointing to this chromosome, along with autosome 8 (which shares histone genes with the B) as putative ancestors of the B chromosome. We found several transposable elements (TEs) showing nucleotidic variation specific to B-carrying individuals, which was also present in B-carrying transcriptomes. Remarkably, an interstitial region of the B chromosome included a 17 kb chimera composed of 29 different TEs, suggesting reiterative TE insertion in this B chromosome region.

A Glimpse into the Satellite DNA Library in Characidae Fish (Teleostei, Characiformes)

Satellite DNA (satDNA) is an abundant fraction of repetitive DNA in eukaryotic genomes and plays an important role in genome organization and evolution. In general, satDNA sequences follow a concerted evolutionary pattern through the intragenomic homogenization of different repeat units. In addition, the satDNA library hypothesis predicts that related species share a series of satDNA variants descended from a common ancestor species, with differential amplification of different satDNA variants. The finding of a same satDNA family in species belonging to different genera within Characidae fish provided the opportunity to test both concerted evolution and library hypotheses. For this purpose, we analyzed here sequence variation and abundance of this satDNA family in ten species, by a combination of next generation sequencing (NGS), PCR and Sanger sequencing, and fluorescence in situ hybridization (FISH). We found extensive between-species variation for the number and size of pericentromeric FISH signals. At genomic level, the analysis of 1000s of DNA sequences obtained by Illumina sequencing and PCR amplification allowed defining 150 haplotypes which were linked in a common minimum spanning tree, where different patterns of concerted evolution were apparent. This also provided a glimpse into the satDNA library of this group of species. In consistency with the library hypothesis, different variants for this satDNA showed high differences in abundance between species, from highly abundant to simply relictual variants.

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.

High-throughput analysis of the satellitome illuminates satellite DNA evolution

Satellite DNA (satDNA) is a major component yet the great unknown of eukaryote genomes and clearly underrepresented in genome sequencing projects. Here we show the high-throughput analysis of satellite DNA content in the migratory locust by means of the bioinformatic analysis of Illumina reads with the RepeatExplorer and RepeatMasker programs. This unveiled 62 satDNA families and we propose the term "satellitome" for the whole collection of different satDNA families in a genome. The finding that satDNAs were present in many contigs of the migratory locust draft genome indicates that they show many genomic locations invisible by fluorescent in situ hybridization (FISH). The cytological pattern of five satellites showing common descent (belonging to the SF3 superfamily) suggests that non-clustered satDNAs can become into clustered through local amplification at any of the many genomic loci resulting from previous dissemination of short satDNA arrays. The fact that all kinds of satDNA (micro- mini- and satellites) can show the non-clustered and clustered states suggests that all these elements are mostly similar, except for repeat length. Finally, the presence of VNTRs in bacteria, showing similar properties to non-clustered satDNAs in eukaryotes, suggests that this kind of tandem repeats show common properties in all living beings.

Comment on Schielzeth et al. (2014): "Genome size variation affects song attractiveness in grasshoppers: Evidence for sexual selection against large genomes"

Schielzeth et al. (2014) concluded that attractive grasshopper singers have significantly smaller genomes thus suggesting a possible role for sexual selection on genome size. Whereas this conclusion could still be conceivably valid, it is not supported by the data presented due to some technical flaws. In addition, the interpretation of the results, speculating on the possible presence of B chromosomes, is not justified.

Origin of B chromosomes in the genus Astyanax (Characiformes, Characidae) and the limits of chromosome painting

Eukaryote genomes are frequently burdened with the presence of supernumerary (B) chromosomes. Their origin is frequently investigated by chromosome painting, under the hypothesis that sharing the repetitive DNA sequences contained in the painting probes is a sign of common descent. However, the intragenomic mobility of many anonymous DNA sequences contained in these probes (e.g., transposable elements) adds high uncertainty to this conclusion. Here we test the validity of chromosome painting to investigate B chromosome origin by comparing its results for seven B chromosome types in two fish species genus Astyanax, with those obtained (1) by means of the physical mapping of 18S ribosomal DNA (rDNA), H1 histone genes, the As51 satellite DNA and the (AC)15 microsatellite, and (2) by comparing the nucleotide sequence of one of these families (ITS regions from ribosomal DNA) between genomic DNA from B-lacking individuals in both species and the microdissected DNA from two metacentric B chromosomes found in these same species. Intra- and inter-specific painting suggested that all B chromosomes that were assayed shared homologous DNA sequences among them, as well as with a variable number of A chromosomes in each species. This finding would be consistent with a common origin for all seven B chromosomes analyzed. By contrast, the physical mapping of repetitive DNA sequences failed to give support to this hypothesis, as no more than two B-types shared a given repetitive DNA. Finally, sequence analysis of the ITS regions suggested that at least some of the B chromosomes could have had a common origin.

Uncovering the Ancestry of B Chromosomes in Moenkhausia sanctaefilomenae (Teleostei, Characidae)

B chromosomes constitute a heterogeneous mixture of genomic parasites that are sometimes derived intraspecifically from the standard genome of the host species, but result from interspecific hybridization in other cases. The mode of origin determines the DNA content, with the B chromosomes showing high similarity with the A genome in the first case, but presenting higher similarity with a different species in the second. The characid fish Moenkhausia sanctaefilomenae harbours highly invasive B chromosomes, which are present in all populations analyzed to date in the Parana and Tietê rivers. To investigate the origin of these B chromosomes, we analyzed two natural populations: one carrying B chromosomes and the other lacking them, using a combination of molecular cytogenetic techniques, nucleotide sequence analysis and high-throughput sequencing (Illumina HiSeq2000). Our results showed that i) B chromosomes have not yet reached the Paranapanema River basin; ii) B chromosomes are mitotically unstable; iii) there are two types of B chromosomes, the most frequent of which is lightly C-banded (similar to euchromatin in A chromosomes) (B₁), while the other is darkly C-banded (heterochromatin-like) (B₂); iv) the two B types contain the same tandem repeat DNA sequences (18S ribosomal DNA, H3 histone genes, MS3 and MS7 satellite DNA), with a higher content of 18S rDNA in the heterochromatic variant; v) all of these repetitive DNAs are present together only in the paracentromeric region of autosome pair no. 6, suggesting that the B chromosomes are derived from this A chromosome; vi) the two B chromosome variants show MS3 sequences that are highly divergent from each other and from the 0B genome, although the B₂-derived sequences exhibit higher similarity with the 0B genome (this suggests an independent origin of the two B variants, with the less frequent, B₂ type presumably being younger); and vii) the dN/dS ratio for the H3.2 histone gene is almost 4–6 times higher for B chromosomes than for A chromosome sequences, suggesting that purifying selection is relaxed for the DNA sequences located on the B chromosomes, presumably because they are mostly inactive.

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.

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.

Intragenomic distribution of RTE retroelements suggests intrachromosomal movement

Much is known about the abundance of transposable elements (TEs) in eukaryotic genomes, but much is still unknown on their behaviour within cells. We employ here a combination of cytological, molecular and genomic approaches providing information on the intragenomic distribution and behaviour of non-long terminal repeat (LTR) retrotransposon-like elements (RTE). We microdissected every chromosome in a single first meiotic metaphase cell of the grasshopper Eyprepocnemis plorans and polymerase chain reaction (PCR) amplified a fragment of the RTE reverse transcriptase gene with specific primers. PCR products were cloned and 139 clones were sequenced. Analysis of molecular variance (AMOVA) showed significant intragenomic structure for these elements, with 4.6 % of molecular variance being found between chromosomes. A maximum likelihood tree built with the RTE sequences revealed the frequent presence of two or more elements showing very high similarity and being located on the same chromosome, thus suggesting intrachromosome movement. The 454 pyrosequencing of genomic DNA gave strong support to the microdissection results and provided evidence for the existence of 5' truncated elements. Our results thus indicate a tendency of RTE elements to reinsert into the same chromosome from where they were transcribed, which could be achieved if retrotranscription and insertion takes place immediately after transcription.

Possible introgression of B chromosomes between bee species (Genus Partamona )

The origin of supernumerary (B) chromosomes is still a debated topic, with intra- and interspecific origins being the most plausible options. In the bee Partamona helleri, a sequence-characterized amplified region (SCAR) marker being specific to B chromosomes suggested the possibility of interspecific origin. Here, we search for this marker in 3 close relative species and perform DNA sequence comparison between species. The SCAR sequence does not show homology with other sequences in the databases, but does contain an open reading frame with sequence homology with a reverse transcriptase. Dot-blot hybridization using the SCAR marker as a probe confirmed that it is present in B-carrying, but not B-lacking larvae of P. helleri, and indicated its presence in adult individuals of P. cupira and P. criptica. Additionally, PCR amplification of the SCAR marker was successful on genomic DNA obtained from P. helleri and P. rustica larvae carrying B chromosomes, and on genomic DNA obtained from adult individuals of P. cupira, P. criptica and P. rustica. Finally, a comparison of the DNA sequence of the SCAR markers amplified from these 4 species showed very few nucleotide differences between the species. The complete association between B chromosome and SCAR presence and the scarce divergence observed for this DNA sequence between the 4 species analyzed suggest the possibility that this B chromosome has recently been transferred between species through several episodes of interspecific hybridization.

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.

Next generation sequencing and FISH reveal uneven and nonrandom microsatellite distribution in two grasshopper genomes

Simple sequence repeats (SSRs), also known as microsatellites, are one of the prominent DNA sequences shaping the repeated fraction of eukaryotic genomes. In spite of their profuse use as molecular markers for a variety of genetic and evolutionary studies, their genomic location, distribution, and function are not yet well understood. Here we report the first thorough joint analysis of microsatellite motifs at both genomic and chromosomal levels in animal species, by a combination of 454 sequencing and fluorescent in situ hybridization (FISH) techniques performed on two grasshopper species. The in silico analysis of the 454 reads suggested that microsatellite expansion is not driving size increase of these genomes, as SSR abundance was higher in the species showing the smallest genome. However, the two species showed the same uneven and nonrandom location of SSRs, with clear predominance of dinucleotide motifs and association with several types of repetitive elements, mostly histone gene spacers, ribosomal DNA intergenic spacers (IGS), and transposable elements (TEs). The FISH analysis showed a dispersed chromosome distribution of microsatellite motifs in euchromatic regions, in coincidence with chromosome location patterns previously observed for many mobile elements in these species. However, some SSR motifs were clustered, especially those located in the histone gene cluster.

Single origin of sex chromosomes and multiple origins of B chromosomes in fish genus Characidium

Chromosome painting with DNA probes obtained from supernumerary (B) and sex chromosomes in three species of fish genus Characidium (C. gomesi, C. pterostictum and C. oiticicai) showed a close resemblance in repetitive DNA content between B and sex chromosomes in C. gomesi and C. pterostictum. This suggests an intraspecific origin for B chromosomes in these two species, probably deriving from sex chromosomes. In C. oiticicai, however, a DNA probe obtained from its B chromosome hybridized with the B but not with the A chromosomes, suggesting that the B chromosome in this species could have arisen interspecifically, although this hypothesis needs further investigation. A molecular phylogenetic analysis performed on nine Characidium species, with two mtDNA genes, showed that the presence of heteromorphic sex chromosomes in these species is a derived condition, and that their origin could have been unique, a conclusion also supported by interspecific chromosome painting with a CgW probe derived from the W chromosome in C. gomesi. Summing up, our results indicate that whereas heteromorphic sex chromosomes in the genus Characidium appear to have had a common and unique origin, B chromosomes may have had independent origins in different species. Our results also show that molecular phylogenetic analysis is an excellent complement for cytogenetic studies by unveiling the direction of evolutionary chromosome changes.

B chromosomes in the grasshopper Eyprepocnemis plorans are present in all body parts analyzed and show extensive variation for rDNA copy number

B chromosomes in the grasshopper Eyprepocnemis plorans are considered to be mitotically stable, because all meiotic (primary spermatocytes and oocytes) or mitotic (embryos, ovarioles, and gastric caecum) cells analyzed within the same individual show the same B chromosome number. Nothing is known, however, about body parts with somatic tissues with no mitotic activity in adult individuals, constituting the immense majority of their body. Therefore, we investigated whether B chromosomes are present in 8 non-mitotically active somatic body parts from both sexes in addition to ovarioles and testes by PCR analysis of 2 B-specific molecular markers. We also elucidated the number of B chromosomes that an individual carried through quantifying the B-located rDNA copy number by qPCR. Our results indicated the amplification of both B-specific markers in all analyzed body parts. However, we found high variation between males for the estimated number of rDNA units in the B chromosomes. These results demonstrate the presence of B chromosomes in all body parts from the same individual and suggest a high variation in the rDNA content of the B chromosomes carried by different individuals from the same population, presumably due to unequal crossovers during meiosis.

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.

Delimiting the origin of a B chromosome by FISH mapping, chromosome painting and DNA sequence analysis in Astyanax paranae (Teleostei, Characiformes)

Supernumerary (B) chromosomes have been shown to contain a wide variety of repetitive sequences. For this reason, fluorescent in situ hybridisation (FISH) is a useful tool for ascertaining the origin of these genomic elements, especially when combined with painting from microdissected B chromosomes. In order to investigate the origin of B chromosomes in the fish species Astyanax paranae, these two approaches were used along with PCR amplification of specific DNA sequences obtained from the B chromosomes and its comparison with those residing in the A chromosomes. Remarkably, chromosome painting with the one-arm metacentric B chromosome probe showed hybridization signals on entire B chromosome, while FISH mapping revealed the presence of H1 histone and 18S rDNA genes symmetrically placed in both arms of the B chromosome. These results support the hypothesis that the B chromosome of A. paranae is an isochromosome. Additionally, the chromosome pairs Nos. 2 or 23 are considered the possible B chromosome ancestors since both contain syntenic H1 and 18S rRNA sequences. The analysis of DNA sequence fragments of the histone and rRNA genes obtained from the microdissected B chromosomes showed high similarity with those obtained from 0B individuals, which supports the intraspecific origin of B chromosomes in A. paranae. Finally, the population hereby analysed showed a female-biased B chromosome presence suggesting that B chromosomes in this species could influence sex determinism.

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.

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.

Ribosomal DNA is active in different B chromosome variants of the grasshopper Eyprepocnemis plorans

B chromosomes are considered to be genetically inert elements. However, some of them are able to show nucleolus organizer region (NOR) activity, as detected by both cytological and molecular means. The grasshopper Eyprepocnemis plorans shows a B chromosome polymorphism characterized by the existence of many B variants. One of them, B24, shows NOR activity in about half of B-carrying males in the Torrox population. Molecular data have suggested the recent origin for B chromosomes in this species, and on this basis it would be expected that NOR activity was widespread among the different B variants. Here we test this hypothesis in four different B chromosome variants (B1, B2, B5, and B24) from 11 natural populations of the grasshopper E. plorans covering the south and east of the Iberian Peninsula plus the Balearic Islands. We used two different approaches: (1) the cytological observation of nucleoli attached to the distal region of the B chromosome (where the rDNA is located), and (2) the molecular detection of the rDNA transcripts carrying an adenine insertion characteristic of B chromosome ITS2 sequences. The results showed NOR expression not only for B24 but also for the B1 and B2 variants. However, the level of B-NOR expression in these latter variants, measured by the proportion of cells showing nucleoli attached to the B chromosomes, was much lower than that previously reported for B24. This suggests the possibility that structural or genetic background conditions are enhancing the expressivity of the rDNA in the B24 variant.

The Ku70 DNA-repair protein is involved in centromere function in a grasshopper species

The Ku70 protein is involved in numerous cell functions, the nonhomologous end joining (NHEJ) DNA repair pathway being the best known. Here, we report a novel function for this protein in the grasshopper Eyprepocnemis plorans. We observed the presence of large Ku70 foci on the centromeres of meiotic and mitotic chromosomes during the cell cycle stages showing the highest centromeric activity (i.e., metaphase and anaphase). The fact that colchicine treatment prevented centromeric location of Ku70, suggests a microtubule-dependent centromeric function for Ku70. Likewise, the absence of Ku70 at metaphase-anaphase centromeres from three males whose Ku70 gene had been knocked down using interference RNA, and the dramatic increase in the frequency of polyploid spermatids observed in these males, suggest that the centromeric presence of Ku70 is required for normal cytokinesis in this species. The centromeric function of Ku70 was not observed in 14 other grasshopper and locust species, or in the mouse, thus suggesting that it is an autapomorphy in E. plorans.

Gypsy, RTE and Mariner transposable elements populate Eyprepocnemis plorans genome

We analyze here the presence and abundance of three types of transposable elements (TEs), i.e. Gypsy, RTE and Mariner, in the genome of the grasshopper Eyprepocnemis plorans. PCR experiments allowed amplification, cloning and sequencing of these elements (EploGypI, EploRTE5, EploMar20) from the E. plorans genome. Fluorescent in situ hybridization (FISH) showed that all three elements are restricted to euchromatic regions, thus being absent from the pericentromeric region of all A chromosomes, which contain a satellite DNA (satDNA) and ribosomal DNA (rDNA), and being very scarce in B chromosomes mostly made up of these two types of repetitive DNA. FISH suggested that EploGypI is the most abundant and EploMar20 is the least abundant, with EploRTE5 showing intermediate abundance. An estimation of copy number, by means of quantitative PCR, showed that EploGypI is, by far, the most abundant element, followed by EploRTE5 and EploMar20, in consistency with FISH results. RNA isolation and PCR experiments on complementary DNA (cDNA) showed the presence of transcripts for the three TE elements. The implications of the preferential location of these TE elements into euchromatin, the significance of TE abundance in the giant genome of this species, and a possible relationship between TEs and B chromosome mutability, are discussed.

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.

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.

New insights on the origin of B chromosomes in Astyanax scabripinnis obtained by chromosome painting and FISH

Chromosome painting (CP) with a probe of B chromosome obtained by microdissection and fluorescence in situ hybridization (FISH) with probes of As51 satellite DNA, C( o )t-1 DNA, and 18S and 5S rDNA confirmed sharing of some repetitive DNA but not rDNA between A and B chromosomes in the fish Astyanax scabripinnis. Meiotic analysis revealed a pachytene B chromosome bivalent nearly half the size of its mitotic configuration, suggesting a self-pairing of B chromosome arms. Such an isochromosome nature of somatic B chromosome was further evidenced by CP and FISH. All the findings obtained suggest (i) intraspecific origin of B chromosome, and (ii) evolutionary enrichment of repetitive DNA classes, especially those contained in the C( o )t-1 and the As51 probes, in B chromosome. However, the precise origin of B chromosome in the present species remains to be elucidated by further molecular cytogenetic analysis because of painting of some A chromosome regions with the B chromosome-derived probe.

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.

Prevalence of B chromosomes in Orthoptera is associated with shape and number of A chromosomes

We analyze the prevalence of B chromosomes in 1,601 species of orthopteran insects where chromosome number and shape are known. B chromosomes have been reported in 191 of these species. Bs are not uniformly distributed among orthopteran superfamilies, with evident hotspots in the Pyrgomorphoidea (32.3% of species carrying Bs), Grylloidea (14.9%), Acridoidea (14.6%) and Tetrigoidea (14.3%). As expected under the theory of centromeric drive, we found a correlation between B chromosome presence and A chromosome shape-Bs are more frequent in karyotypes with more acrocentric A chromosomes. We also found that Bs are less common in species with high chromosome numbers and appear to be most common at the modal chromosome number (2n = 24). Study effort, measured for each genus, was not associated with B prevalence, A chromosome shape or A chromosome number. Our results thus provide support for centromeric drive as an important and prevalent force in the karyotypic evolution of Orthoptera, just as it appears to be in mammals. We suggest that centromeric drive may provide a mechanistic explanation for White's principle of karyotypic orthoselection.