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.

Satellite DNA Is an Inseparable Fellow Traveler of B Chromosomes.

Next-Generation Sequencing (NGS) has revealed that B chromosomes in several species are enriched in repetitive DNA, mostly satellite DNA (satDNA). This raises the question of whether satDNA is important to B chromosomes for functional reasons or else its abundance on Bs is simply a consequence of properties of B chromosomes such as their dispensability and late replication. Here we review current knowledge in this respect and contextualize it within the frame of practical difficulties to perform this kind of research, the most important being the absence of good full genome sequencing for B-carrying species, which is an essential requisite to ascertain the intragenomic origin of B chromosomes. Our review analysis on 16 species revealed that 38% of them showed B-specific satDNAs whereas only one of them (6%) carried an inter-specifically originated B chromosome. This shows that B-specific satDNA families can eventually evolve in intraspecifically arisen B chromosomes. Finally, the possibility of satDNA accumulation on B chromosomes for functional reasons is exemplified by B chromosomes in rye, as they contain B-specific satDNAs which are transcribed and occupy chromosome locations where they might facilitate the kind of drive shown by this B chromosome during pollen grain mitosis.

Gene expression changes elicited by a parasitic B chromosome in the grasshopper Eyprepocnemis plorans are consistent with its phenotypic effects.

Parasitism evokes adaptive physiological changes in the host, many of which take place through gene expression changes. This response can be more or less local, depending on the organ or tissue affected by the parasite, or else systemic when the parasite affects the entire host body. The most extreme of the latter cases is intragenomic parasitism, where the parasite is present in all host nuclei as any other genomic element. Here, we show the molecular crosstalk between a parasitic chromosome (also named B chromosome) and the host genome, manifested through gene expression changes. The transcriptome analysis of 0B and 1B females of the grasshopper Eyprepocnemis plorans, validated by a microarray experiment performed on four B-lacking and five B-carrying females, revealed changes in gene expression for 188 unigenes being consistent in both experiments. Once discarded B-derived transcripts, there were 46 differentially expressed genes (30 up- and 16 downregulated) related with the adaptation of the host genome to the presence of the parasitic chromosome. Interestingly, the functions of these genes could explain some of the most important effects of B chromosomes, such as nucleotypic effects derived from the additional DNA they represent, chemical defense and detoxification, protein modification and response to stress, ovary function, and regulation of gene expression. Collectively, these changes uncover an intimate host-parasite interaction between A and B chromosomes during crucial steps of gene expression and protein function.

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.

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.

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.

Satellite DNA content illuminates the ancestry of a supernumerary (B) chromosome.

B chromosomes are supernumerary genomic elements most likely derived from the standard (A) chromosomes, whose dispensability has freed their DNA sequences to evolve fast, thus making it difficult to uncover their ancestry. Here, we show the ancestry of a B chromosome in the grasshopper Eumigus monticola by means of the high-throughput analysis of the satellitome, i.e., the whole collection of satellite DNA (satDNA). The satellitome found in this species consists of 27 satDNA families, with monomer length between 5 and 325 nt and A + T content between 42.9 and 83.3 %. Two out of the 20 clustered satDNA families (EmoSat26-41 and EmoSat27-102) were observed only on the B chromosome. The A chromosome carrying the highest number of satDNA families was the megameric S8 (13 families), six of which were also present in the B chromosome, and three of these were exclusive of the S8 and B chromosomes. The absence in the B chromosome of the H3 histone gene cluster (located interstitially on S8) and three satDNA families (located distally on S8) allowed delimiting the possible origin of the B chromosome to the proximal third of the S8 autosome, through a breakpoint between EmoSat11-122 and the H3 cluster. Interestingly, bioinformatic analysis revealed the presence of seeds for the two B-specific satDNAs in the A chromosomes, suggesting their massive amplification in the B chromosome after its origin. Therefore, intraspecifically arisen B chromosomes can harbor DNA sequences apparently being B-specific.

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.

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

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

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.

Transient Microgeographic Clines during B Chromosome Invasion.

The near-neutral model of B chromosome evolution predicts that the invasion of a new population should last some tens of generations, but the details on how it proceeds in real populations are mostly unknown. Trying to fill this gap, we analyze here a natural population of the grasshopper Eyprepocnemis plorans at three time points during the last 35 years. Our results show that B chromosome frequency increased significantly during this period and that a cline observed in 1992 had disappeared in 2012 once B chromosome frequency reached an upper limit at all sites sampled. This indicates that, during B chromosome invasion, transient clines for B chromosome frequency are formed at the invasion front on a microgeographic scale. Computer simulation experiments showed that the pattern of change observed for genotypic frequencies is consistent with the existence of B chromosome drive through females and selection against individuals with a high number of B chromosomes.

HP1 knockdown is associated with abnormal condensation of almost all chromatin types in a grasshopper (Eyprepocnemis plorans).

Heterochromatin protein 1 (HP1) is a highly conserved family of eukaryotic proteins required for heterochromatic gene silencing and euchromatic gene transcription regulation. In addition, HP1 is involved in chromatin organization and protection of chromosome integrity during cell division. Here, we present a cytological and molecular analysis of the effects of HP1 knockdown in Eyprepocnemis plorans, a grasshopper species polymorphic for supernumerary heterochromatic chromosomes. Our results revealed contrasting effects of HP1 knockdown on gene activity. While the Bub1 gene decreased in expression level in HP1 knockdown animals, NOR activity, rRNA and, contrarily to previous reports in Drosophila, Hsp70 gene expression remained unchanged. Furthermore, HP1 knockdown resulted in abnormal chromatin condensation, chromosomal bridges, higher frequency of macrospermatids, loss of muscle mass and hemolymph amount as well as a low number of dividing cells and survival reduction. All these phenotypes are very likely due to the chromatin condensation disruption observed for almost all kinds of chromatin.

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.

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.

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.

Tandem Repeat DNA Provides Many Cytological Markers for Hybrid Zone Analysis in Two Subspecies of the Grasshopper Chorthippus parallelus.

Recent advances in next generation sequencing (NGS) have greatly increased our understanding of non-coding tandem repeat (TR) DNA. Here we show how TR DNA can be useful for the study of hybrid zones (HZ), as it serves as a marker to identify introgression in areas where two biological entities come in contact. We used Illumina libraries to analyse two subspecies of the grasshopper Chorthippus parallelus, which currently form a HZ in the Pyrenees. We retrieved a total of 152 TR sequences, and used fluorescent in situ hybridization (FISH) to map 77 families in purebred individuals from both subspecies. Our analysis revealed 50 TR families that could serve as markers for analysis of this HZ, using FISH. Differential TR bands were unevenly distributed between chromosomes and subspecies. Some of these TR families yielded FISH bands in only one of the subspecies, suggesting the amplification of these TR families after the geographic separation of the subspecies in the Pleistocene. Our cytological analysis of two TR markers along a transect of the Pyrenean hybrid zone showed asymmetrical introgression of one subspecies into the other, consistent with previous findings using other markers. These results demonstrate the reliability of TR-band markers for hybrid zone studies.

Evolutionary dynamics of 5S rDNA location in acridid grasshoppers and its relationship with H3 histone gene and 45S rDNA location.

We analyze the chromosomal location of 5S rDNA clusters in 29 species of grasshoppers belonging to the family Acrididae. There was extensive variation among species for the number and location of 5S rDNA sites. Out of 148 sites detected, 75% were proximally located, 21.6% were interstitial, and only 3.4% were distal. The number of 5S rDNA sites per species varied from a single chromosome pair (in six species) to all chromosome pairs (in five species), with a range of intermediate situations. Thirteen chromosomes from eight species carried two 5S rDNA clusters. At intraspecific level, differences among populations were detected in Eyprepocnemis plorans, and some heteromorphisms have also been observed in some species. Double FISH for 5S rDNA and H3 histone gene DNA, performed on 17 of these 29 species, revealed that both markers are sometimes placed in a same chromosome but at different location, whereas they appeared to co-localize in five species (Calliptamus barbarus, Heteracris adpersa, Aiolopus strepens, Oedipoda charpentieri and O. coerulescens). Double fiber-FISH in A. strepens and O. coerulescens showed that the two DNAs are closely interspersed with variable relative amounts of both classes of DNA. Finally, no correlation was observed between the number of 5S and 45S rDNA clusters in 23 species where this information was available. These results are discussed in the light of possible mechanisms of spread that led to the extensive variation in the number of clusters observed for both rDNA types in acridid grasshoppers.

Transcription of a B chromosome CAP-G pseudogene does not influence normal Condensin Complex genes in a grasshopper.

Parasitic B chromosomes invade and persist in natural populations through several mechanisms for transmission advantage (drive). They may contain gene-derived sequences which, in some cases, are actively transcribed. A further interesting question is whether B-derived transcripts become functional products. In the grasshopper Eyprepocnemis plorans, one of the gene-derived sequences located on the B chromosome shows homology with the gene coding for the CAP-G subunit of condensin I. We show here, by means of fluorescent in situ hybridization coupled with tyramide signal amplification (FISH-TSA), that this gene is located in the distal region of the B24 chromosome variant. The DNA sequence located in the B chromosome is a pseudogenic version of the CAP-G gene (B-CAP-G). In two Spanish populations, we found active transcription of B-CAP-G, but it did not influence the expression of CAP-D2 and CAP-D3 genes coding for corresponding condensin I and II subunits, respectively. Our results indicate that the transcriptional regulation of the B-CAP-G pseudogene is uncoupled from the standard regulation of the genes that constitute the condensin complex, and suggest that some of the B chromosome known effects may be related with its gene content and transcriptional activity, thus opening new exciting avenues for research.

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.