Comparative FISH analysis in five species of Eyprepocnemidine grasshoppers

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

Origin and Evolution of the Neo-Sex Chromosomes in Pamphagidae Grasshoppers through Chromosome Fusion and Following Heteromorphization

In most phylogenetic lineages, the evolution of sex chromosomes is accompanied by their heteromorphization and degradation of one of them. The neo-sex chromosomes are useful model for studying early stages of these processes. Recently two lineages of the neo-sex chromosomes on different stages of heteromorphization was discovered in Pamphagidae family. The neo-sex chromosome heteromorphization was analyzed by generation of DNA probes derived from the neo-Xs and neo-Ys followed with chromosome painting in nineteen species of Pamphagidae family. The homologous regions of the neo-sex chromosomes were determined in closely related species with the painting procedure and image analysis with application of the Visualization of the Specific Signal in Silico software package. Results of these analyses and distribution of C-positive regions in the neo-sex chromosomes revealed details of the heteromorphization of the neo-sex chromosomes in species from both phylogenetic lineages of Pamphagidae grasshoppers. The hypothetical mechanism of the neo-Y degradation was suggested. It includes expansion of different repeats from the proximal neo-Y chromosome region by inversions, spreading them towards distal region. Amplification of these repeats leads to formation of C-positive regions and elimination of the C-negative regions located between them.

De Novo Assembly and Characterization of the Transcriptome of Grasshopper Shirakiacris shirakii

Background: The grasshopper Shirakiacris shirakii is an important agricultural pest and feeds mainly on gramineous plants, thereby causing economic damage to a wide range of crops. However, genomic information on this species is extremely limited thus far, and transcriptome data relevant to insecticide resistance and pest control are also not available. Methods: The transcriptome of S. shirakii was sequenced using the Illumina HiSeq platform, and we de novo assembled the transcriptome. Results: Its sequencing produced a total of 105,408,878 clean reads, and the de novo assembly revealed 74,657 unigenes with an average length of 680 bp and N50 of 1057 bp. A total of 28,173 unigenes were annotated for the NCBI non-redundant protein sequences (Nr), NCBI non-redundant nucleotide sequences (Nt), a manually-annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Based on the Nr annotation results, we manually identified 79 unigenes encoding cytochrome P450 monooxygenases (P450s), 36 unigenes encoding carboxylesterases (CarEs) and 36 unigenes encoding glutathione S-transferases (GSTs) in S. shirakii. Core RNAi components relevant to miroRNA, siRNA and piRNA pathways, including Pasha, Loquacious, Argonaute-1, Argonaute-2, Argonaute-3, Zucchini, Aubergine, enhanced RNAi-1 and Piwi, were expressed in S. shirakii. We also identified five unigenes that were homologous to the Sid-1 gene. In addition, the analysis of differential gene expressions revealed that a total of 19,764 unigenes were up-regulated and 4185 unigenes were down-regulated in larvae. In total, we predicted 7504 simple sequence repeats (SSRs) from 74,657 unigenes. Conclusions: The comprehensive de novo transcriptomic data of S. shirakii will offer a series of valuable molecular resources for better studying insecticide resistance, RNAi and molecular marker discovery in the transcriptome.

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.

U1 snDNA clusters in grasshoppers: chromosomal dynamics and genomic organization

The spliceosome, constituted by a protein set associated with small nuclear RNA (snRNA), is responsible for mRNA maturation through intron removal. Among snRNA genes, U1 is generally a conserved repetitive sequence. To unveil the chromosomal/genomic dynamics of this multigene family in grasshoppers, we mapped U1 genes by fluorescence in situ hybridization in 70 species belonging to the families Proscopiidae, Pyrgomorphidae, Ommexechidae, Romaleidae and Acrididae. Evident clusters were observed in all species, indicating that, at least, some U1 repeats are tandemly arrayed. High conservation was observed in the first four families, with most species carrying a single U1 cluster, frequently located in the third or fourth longest autosome. By contrast, extensive variation was observed among Acrididae, from a single chromosome pair carrying U1 to all chromosome pairs carrying it, with occasional occurrence of two or more clusters in the same chromosome. DNA sequence analysis in Eyprepocnemis plorans (species carrying U1 clusters on seven different chromosome pairs) and Locusta migratoria (carrying U1 in a single chromosome pair) supported the coexistence of functional and pseudogenic lineages. One of these pseudogenic lineages was truncated in the same nucleotide position in both species, suggesting that it was present in a common ancestor to both species. At least in E. plorans, this U1 snDNA pseudogenic lineage was associated with 5S rDNA and short interspersed elements (SINE)-like mobile elements. Given that we conclude in grasshoppers that the U1 snDNA had evolved under the birth-and-death model and that its intragenomic spread might be related with mobile elements.

Chromosomal divergence and evolutionary inferences in Rhodniini based on the chromosomal location of ribosomal genes

In this study, we used fluorescence in situ hybridisation to determine the chromosomal location of 45S rDNA clusters in 10 species of the tribe Rhodniini (Hemiptera: Reduviidae: Triatominae). The results showed striking inter and intraspecific variability, with the location of the rDNA clusters restricted to sex chromosomes with two patterns: either on one (X chromosome) or both sex chromosomes (X and Y chromosomes). This variation occurs within a genus that has an unchanging diploid chromosome number (2n = 22, including 20 autosomes and 2 sex chromosomes) and a similar chromosome size and genomic DNA content, reflecting a genome dynamic not revealed by these chromosome traits. The rDNA variation in closely related species and the intraspecific polymorphism in Rhodnius ecuadoriensis suggested that the chromosomal position of rDNA clusters might be a useful marker to identify recently diverged species or populations. We discuss the ancestral position of ribosomal genes in the tribe Rhodniini and the possible mechanisms involved in the variation of the rDNA clusters, including the loss of rDNA loci on the Y chromosome, transposition and ectopic pairing. The last two processes involve chromosomal exchanges between both sex chromosomes, in contrast to the widely accepted idea that the achiasmatic sex chromosomes of Heteroptera do not interchange sequences.

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.

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

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

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.

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.

Cytogenetic and molecular markers reveal the complexity of the genus Piabina Reinhardt, 1867 (Characiformes: Characidae)

Cytogenetic and molecular analyses were carried out in fish representative of the genus Piabina. This study specifically involved the species P. argentea and P. anhembi collected from areas of the Paranapanema and Tietê River basins, Brazil. Our findings suggest that fish classified as Piabina argentea in the Paranapanema and Tietê Rivers may represent more than one species. The samples analyzed differed by cytogenetic particularities and molecular analyses using partial sequences of the genes COI and CytB as genetic markers revealed three distinct groups of P. argentea with genetic distances sufficient to support the conclusion that the three samples analyzed are three distinct taxonomic units.

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

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

Location and expression of ribosomal RNA genes in grasshoppers: abundance of silent and cryptic loci

We investigate regularities and restrictions in chromosome location of ribosomal RNA genes, analysed by fluorescent in situ hybridization (FISH), and their phenotypic expression assessed by nucleolus formation at first meiotic prophase cells, analysed by silver impregnation, in 49 grasshopper species. High variation was found for rDNA location between species within most genera analysed. The mean haploid number of rDNA loci detected by FISH was 2.47, but some species had up to 10 loci. Chromosome distribution of rDNA loci differed between the Gomphocerinae and Oedipodinae subfamilies, most loci being proximal to the centromere in the former and distal to it in the latter. Chromosomes 2, 3 and X frequently carried rDNA in Gomphocerinae species with 2n male symbol=17 chromosomes, whereas chromosomes 6 and 9 were the most frequent rDNA locations in the Oedipodinae. About 13% of the 126 rDNA loci detected by FISH were silent, although this figure might be even higher. The comparison of FISH and silver-impregnation results also suggested the existence of cryptic NORs, i.e. those forming small nucleoli with no apparent presence of rDNA revealed by FISH. This was especially clear after the same cells in two species were sequentially treated with both silver impregnation and FISH. The abundance of silent and cryptic loci might thus suggest that rDNA spreads through grasshopper genomes by the Dubcovsky and Dvorak mechanism-that is, the transposition of a few rRNA genes to new chromosome locations, their amplification giving rise to new NORs, and the elimination of the old NORs. The cryptic NORs might correspond to nascent NORs, i.e. a few rRNA gene copies moved to new locations, whereas the inactive rDNA loci might correspond to those being in the process of elimination.

Physical mapping of rDNA and satDNA in A and B chromosomes of the grasshopper Eyprepocnemis plorans from a Greek population

Adult males and females of the grasshopper Eyprepocnemis plorans from a Greek population were analysed by C-banding, silver impregnation and double FISH for two DNA probes, i.e. ribosomal DNA (rDNA) and a 180-bp tandem repeat DNA (satDNA). This population shows characteristics of rDNA location in A chromosomes that are intermediate between those previously reported for eastern (Caucasus) and western (Spain and Morocco) populations. The four rDNA clusters revealed by FISH in chromosomes X, 9, 10 and 11 in Greek specimens imply two more than the two observed in chromosomes 9 and 11 in the Caucasus, but less than the 12 observed in all chromosomes in Morocco. Remarkably, the X chromosome bears one of the new rDNA locations in Greece with respect to the Caucasus, but it appears to be inactive, in contrast to X chromosomes in western populations, which are usually active. B chromosomes were very frequent in the Greek population, and three variants differing in size were observed, all of these being largely composed of rDNA, with the exception of a small pericentromeric satDNA cluster. The high B frequency suggests that B chromosomes in this population might behave parasitically, in resemblance to Bs in western populations.

Comparative analysis of rDNA location in five Neotropical gomphocerine grasshopper species

We report here, for the first time, the chromosome complement, number and location of the nucleolar organizer regions (NORs) revealed by silver staining (AgNO(3)) and fluorescent in situ hybridization (FISH) in five Neotropical gomphocerine species: Rhammatocerus brasiliensis, R. brunneri, R. palustris, R. pictus and Amblytropidia sp. The objective of this study was to summarize available data and propose a model of chromosome evolution in Neotropical gomphocerines. All five species studied showed chromosome numbers consisting of 2n = 23,X0 in males and 2n = 24,XX in females. Amblytropidia sp. was the only species showing a bivalent (M(8)) with megameric behavior during meiosis. The rDNA sites were restricted to autosomal pairs, i.e. the pericentromeric region of the S(9) chromosome, the consensus NOR location in all five species. R. brasiliensis was the only species showing additional NORs on M(4) and M(6) pairs which, likewise the S(9) NOR, were active in all cells analyzed. Comparison of these results with those reported previously in Palearctic gomphocerine species suggests higher resemblance of Neotropical species with the Old World species also possessing 23/24 chromosomes. Evolutionary mechanisms responsible for the observed interspecific variation in NOR location in this group are discussed.

Molecular marker systems in insects: current trends and future avenues

Insects comprise the largest species composition in the entire animal kingdom and possess a vast undiscovered genetic diversity and gene pool that can be better explored using molecular marker techniques. Current trends of application of DNA marker techniques in diverse domains of insect ecological studies show that mitochondrial DNA (mtDNA), microsatellites, random amplified polymorphic DNA (RAPD), expressed sequence tags (EST) and amplified fragment length polymorphism (AFLP) markers have contributed significantly for progresses towards understanding genetic basis of insect diversity and for mapping medically and agriculturally important genes and quantitative trait loci in insect pests. Apart from these popular marker systems, other novel approaches including transposon display, sequence-specific amplification polymorphism (S-SAP), repeat-associated polymerase chain reaction (PCR) markers have been identified as alternate marker systems in insect studies. Besides, whole genome microarray and single nucleotide polymorphism (SNP) assays are becoming more popular to screen genome-wide polymorphisms in fast and cost effective manner. However, use of such methodologies has not gained widespread popularity in entomological studies. The current study highlights the recent trends of applications of molecular markers in insect studies and explores the technological advancements in molecular marker tools and modern high throughput genotyping methodologies that may be applied in entomological researches for better understanding of insect ecology at molecular level.

Population variation in the A chromosome distribution of satellite DNA and ribosomal DNA in the grasshopper Eyprepocnemis plorans

The double FISH analysis of two repetitive DNAs (a satellite DNA and ribosomal DNA) in 12 natural populations of the grasshopper Eyprepocnemis plorans collected at the south (Granada and Málaga provinces) and south-east (Albacete and Murcia provinces) of the Iberian Peninsula has shown their wide-spread presence throughout the whole genome as well as extensive variation among populations. Both DNAs are found in most A chromosomes. Regularly, both DNAs occurred in the S11 and X chromosomes, rDNA in the S10 and satDNA in the L2 and M3. No correlation was found between the number of satDNA and rDNA clusters in the A genomes of the 12 populations analysed, and both figures were independent of the presence of B chromosomes. The genomic distribution of both DNAs showed no association with the geographical localization of the populations analysed. Finally, we provide evidence that the supernumerary chromosome segment proximally located on the S11 chromosome is, in most cases, the result of satDNA amplification but, in some cases, it might also derive from amplification of both satDNA and rDNA.

Multiregional origin of B chromosomes in the grasshopper Eyprepocnemis plorans

Analysis of chromosome localization of three molecular markers, 18S-5.8S-28S rDNA, 5S rDNA and a 180 bp satDNA, showed that B chromosomes in the grasshopper Eyprepocnemis plorans originated independently in Eastern (Caucasus) and Western (Spain and Morocco) populations. Eastern B chromosomes are most likely derived from the smallest autosome, which is the only A chromosome carrying the three markers, in coincidence with Caucasian B chromosomes. Western B chromosomes, however, lack 5S rDNA and are most likely derived from the X chromosome, which is the only A chromosome carrying the two remaining markers, always in the same order with respect to the centromere, as the B chromosome.