Quantitative analysis of NOR expression in a B chromosome of the grasshopper Eyprepocnemis plorans

Transcribing the enigma: the B chromosome as a territory of uncharted RNAs

B chromosomes are supernumerary elements found in several groups of eukaryotes, including fungi, plants, and animals. Typically, these chromosomes either originate from their hosts through errors in meiosis or interspecifically through horizontal transfer. While many B chromosomes are primarily heterochromatic and possess a low number of coding genes, these additional elements are still capable of transcribing sequences and exerting influence on the expression of host genes. How B chromosomes escape elimination and which impacts can be promoted in the cell always intrigued the cytogeneticists. In pursuit of understanding the behavior and functional impacts of these extra elements, cytogenetic studies meet the advances of molecular biology, incorporating various techniques into investigating B chromosomes from a functional perspective. In this review, we present a timeline of studies investigating B chromosomes and RNAs, highlighting the advances and key findings throughout their history. Additionally, we identified which RNA classes are reported in the B chromosomes and emphasized the necessity for further investigation into new perspectives on the B chromosome functions. In this context, we present a phylogenetic tree that illustrates which branches either report B chromosome presence or have functional RNA studies related to B chromosomes. We propose investigating other unexplored RNA classes and conducting functional analysis in conjunction with cytogenetic studies to enhance our understanding of the B chromosome from an RNA perspective.

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.

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.

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.

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.

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.

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.

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.

Chromosomal Mapping of Repetitive DNAs in the Grasshopper Abracris flavolineata Reveal Possible Ancestry of the B Chromosome and H3 Histone Spreading

Supernumerary chromosomes (B chromosomes) occur in approximately 15% of eukaryote species. Although these chromosomes have been extensively studied, knowledge concerning their specific molecular composition is lacking in most cases. The accumulation of repetitive DNAs is one remarkable characteristic of B chromosomes, and the occurrence of distinct types of multigene families, satellite DNAs and some transposable elements have been reported. Here, we describe the organization of repetitive DNAs in the A complement and B chromosome system in the grasshopper species Abracris flavolineata using classical cytogenetic techniques and FISH analysis using probes for five multigene families, telomeric repeats and repetitive C0t-1 DNA fractions. The 18S rRNA and H3 histone multigene families are highly variable and well distributed in A. flavolineata chromosomes, which contrasts with the conservation of U snRNA genes and less variable distribution of 5S rDNA sequences. The H3 histone gene was an extensively distributed with clusters occurring in all chromosomes. Repetitive DNAs were concentrated in C-positive regions, including the pericentromeric region and small chromosomal arms, with some occurrence in C-negative regions, but abundance was low in the B chromosome. Finally, the first demonstration of the U2 snRNA gene in B chromosomes in A. flavolineata may shed light on its possible origin. These results provide new information regarding chromosomal variability for repetitive DNAs in grasshoppers and the specific molecular composition of B chromosomes.

Comparative cytogenetics in four species of Palinuridae: B chromosomes, ribosomal genes and telomeric sequences

The evolutionary pathway of Palinuridae (Crustacea, Decapoda) is still controversial, uncertain and unexplored, expecially from a karyological point of view. Here we describe the South African spiny lobster Jasus lalandii karyotype: n and 2n values, heterochromatin distribution, nucleolar organizer region (NOR) location and telomeric repeat structure and location. To compare the genomic and chromosomal organization in Palinuridae we located NORs in Panulirus regius, Palinurus gilchristi and Palinurus mauritanicus: all species showed multiple NORs. In J. lalandii NORs were located on three chromosome pairs, with interindividual polymorphism. In P. regius and in the two Palinurus species NORs were located on two chromosome pairs. In the two last species 45S ribosomal gene loci were also found on B chromosomes. In addition, the nature and location of telomeric repeats were investigated by FISH in J. lalandii, P. gilchristi, P. mauritanicus Palinurus elephas, and P. regius (Palinuridae, Achelata), and in Scyllarus arctus (Scyllaridae, Achelata): all these Achelata species showed the (TTAGG)n pentameric repeats. Furthermore, in J. lalandii these repeats occurred in all the telomeres and in some interstitial chromosomal sites, associated with NORs.

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.

Cytogenetic analysis of B chromosomes in one population of the fish Moenkhausia sanctaefilomenae (Steindachner, 1907) (Teleostei, Characiformes)

The aim of this study was to characterize cytogenetically one population of the fish Moenkhausia sanctaefilomenae (Steindachner, 1907), with emphasis on the analysis of B chromosomes. The nucleolar activity in the B microchromosomes was characterized, and an analysis of mitotic instability of these microchromosomes was accomplished. The results showed a diploid chromosome number of 50 chromosomes. In all individuals, we observed the presence of B microchromosomes with intra- and inter-individual variability. The analysis of the nucleolus organizing regions (NORs) by silver nitrate staining demonstrated multiple NORs. We observed active sites of ribosomal DNA in the B microchromosomes, with a frequency of 20% in the analyzed cells, which shows gene activity in these chromosomal elements. The analysis of constitutive heterochromatin patterns showed that the B microchromosomes are heterochromatic or euchromatic, which demonstrates differentiation of DNA composition between these genomic elements. The calculation of the mitotic instability index implied that B chromosomes in this species might be in a final stage of instability.

Chromosomal mapping of rDNAs and H3 histone sequences in the grasshopper rhammatocerus brasiliensis (acrididae, gomphocerinae): extensive chromosomal dispersion and co-localization of 5S rDNA/H3 histone clusters in the A complement and B chromosome

BACKGROUND

Supernumerary B chromosomes occur in addition to standard karyotype and have been described in about 15% of eukaryotes, being the repetitive DNAs the major component of these chromosomes, including in some cases the presence of multigene families. To advance in the understanding of chromosomal organization of multigene families and B chromosome structure and evolution, the distribution of rRNA and H3 histone genes were analyzed in the standard karyotype and B chromosome of three populations of the grasshopper Rhammatocerus brasiliensis.

RESULTS

The location of major rDNA was coincident with the previous analysis for this species. On the other hand, the 5S rDNA mapped in almost all chromosomes of the standard complement (except in the pair 11) and in the B chromosome, showing a distinct result from other populations previously analyzed. Besides the spreading of 5S rDNA in the genome of R. brasiliensis it was also observed multiple sites for H3 histone genes, being located in the same chromosomal regions of 5S rDNAs, including the presence of the H3 gene in the B chromosome.

CONCLUSIONS

Due to the intense spreading of 5S rRNA and H3 histone genes in the genome of R. brasiliensis, their chromosomal distribution was not informative in the clarification of the origin of B elements. Our results indicate a linked organization for the 5S rRNA and H3 histone multigene families investigated in R. brasiliensis, reinforcing previous data concerning the association of both genes in some insect groups. The present findings contribute to understanding the organization/evolution of multigene families in the insect genomes.