U1 snDNA clusters in grasshoppers: chromosomal dynamics and genomic organization

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Karyotypes of Manatees: New Insights into Hybrid Formation (Trichechus inunguis × Trichechus m. manatus) in the Amazon Estuary

Great efforts have been made to preserve manatees. Recently, a hybrid zone was described between Trichechus inunguis (TIN) and the Trichechus manatus manatus (TMM) in the Amazon estuary. Cytogenetic data on these sirenians are limited, despite being fundamental to understanding the hybridization/introgression dynamics and genomic organization in Trichechus. We analyzed the karyotype of TMM, TIN, and two hybrid specimens ("Poque" and "Vitor") by classical and molecular cytogenetics. G-band analysis revealed that TMM (2n = 48) and TIN (2n = 56) diverge by at least six Robertsonian translocations and a pericentric inversion. Hybrids had 2n = 50, however, with Autosomal Fundamental Number (FNA) = 88 in "Poque" and FNA = 74 in "Vitor", and chromosomal distinct pairs in heterozygous; additionally, "Vitor" exhibited heteromorphisms and chromosomes whose pairs could not be determined. The U2 snDNA and Histone H3 multi genes are distributed in small clusters along TIN and TMM chromosomes and have transposable Keno and Helitron elements (TEs) in their sequences. The different karyotypes observed among manatee hybrids may indicate that they represent different generations formed by crossing between fertile hybrids and TIN. On the other hand, it is also possible that all hybrids recorded represent F1 and the observed karyotype differences must result from mechanisms of elimination.

Major and minor U small nuclear RNAs genes characterization in a neotropical fish genome: Chromosomal remodeling and repeat units dispersion in Parodontidae

In association with many proteins, small nuclear RNAs (snRNAs) organize the spliceosomes that play a significant role in processing precursor mRNAs during gene expression. According to snRNAs genic arrangements, two kinds of spliceosomes (major and minor) can be organized into eukaryotic cells. Although in situ localization of U1 and U2 snDNAs have been performed in fish karyotypes, studies with genomic characterization and functionality of U snRNAs integrated into chromosomal changes on Teleostei are still scarce. This study aimed to achieve a genomic characterization of the U snRNAs genes in Apareiodon sp. (2n = 54, ZZ/ZW), apply these data to recognize functional/defective copies, and map chromosomal changes involving snDNAs in Parodontidae species karyotype diversification. Nine snRNA multigene families (U1, U2, U4, U5, U6, U11, U12, U4atac and U6atac) arranged in putatively functional copies in the genome were analyzed. Proximal Sequence Elements (PSE) and TATA-box promoters occurrence, besides an entire transcribed region and conserved secondary structures, qualify them for spliceosome activity. In addition, several defective copies or pseudogenes were identified for the snRNAs that make up the major spliceosome. In situ localization of snDNAs in five species of Parodontidae demonstrated that U1, U2, and U4 snDNAs were involved in chromosomal location changes or units dispersion. The U snRNAs defective/pseudogenes units dispersion could be favored by the probable occurrence of active retrotransposition enzymes in the Apareiodon genome. The U2 and U4 snDNAs sites were involved in independent events in the differentiation of sex chromosomes among Parodontidae lineages. The study characterized U snRNA genes that compose major and minor spliceosomes in the Apareiodon sp. genome and proposes that their defective copies trigger chromosome differentiation and diversification events in Parodontidae.

New insights into the six decades of Mesa’s hypothesis of chromosomal evolution in Ommexechinae grasshoppers (Orthoptera: Acridoidea)

In Acridoidea grasshoppers, chromosomal rearrangements are frequently found as deviations from the standard acrocentric karyotype (2n = 23♂/24♀, FN = 23♂/24♀) in either phylogenetically unrelated species or shared by closely related ones, i.e. genus. In the South American subfamily Ommexechinae, most of the species show a unique karyotype (2n = 23♂/24♀, FN = 25♂/26♀) owing to the occurrence of a large autosomal pair (L1) with submetacentric morphology. In the early 1960s, Alejo Mesa proposed the hypothesis of an ancestral pericentric inversion to explain this karyotype variation. Furthermore, in Ommexechinae, extra chromosomal rearrangements (e.g. centric fusions) are recorded between the ancestral X chromosome and autosomes that originated the so-called neo-sex chromosomes. However, the evolutionary significance of the pericentric inversions and centric fusions in Ommexechinae remains poorly explored. Aiming for a better understanding of chromosomal evolution in Ommexechinae, we performed a detailed cytogenetic analysis in five species. Our findings support the hypothesis about the occurrence of an early pericentric inversion in the ancestor of Ommexechinae. Moreover, our results show a complex karyotype diversification pattern due to several chromosome rearrangements, variations in heterochromatin and repetitive DNA dynamics. Finally, the chromosomal mapping of U2 snDNA in L1 provided new insights about the morphological evolution of this autosomal pair and revealed unnoticed chromosome reorganizations.

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.

Large-scale comparative analysis of cytogenetic markers across Lepidoptera

Fluorescence in situ hybridization (FISH) allows identification of particular chromosomes and their rearrangements. Using FISH with signal enhancement via antibody amplification and enzymatically catalysed reporter deposition, we evaluated applicability of universal cytogenetic markers, namely 18S and 5S rDNA genes, U1 and U2 snRNA genes, and histone H3 genes, in the study of the karyotype evolution in moths and butterflies. Major rDNA underwent rather erratic evolution, which does not always reflect chromosomal changes. In contrast, the hybridization pattern of histone H3 genes was well conserved, reflecting the stable organisation of lepidopteran genomes. Unlike 5S rDNA and U1 and U2 snRNA genes which we failed to detect, except for 5S rDNA in a few representatives of early diverging lepidopteran lineages. To explain the negative FISH results, we used quantitative PCR and Southern hybridization to estimate the copy number and organization of the studied genes in selected species. The results suggested that their detection was hampered by long spacers between the genes and/or their scattered distribution. Our results question homology of 5S rDNA and U1 and U2 snRNA loci in comparative studies. We recommend the use of histone H3 in studies of karyotype evolution.

Universal fluorescence in situ hybridization (FISH) protocol for mapping repetitive DNAs in insects and other arthropods

Repetitive DNAs comprise large portion of eukaryote genomes. In genome projects, the assembly of repetitive DNAs is challenging due to the similarity between repeats, which generate ambiguities for alignment. Fluorescence in situ hybridization (FISH) is a powerful technique for the physical mapping of various sequences on chromosomes. This technique is thus very helpful in chromosome-based genome assemblies, providing information on the fine architecture of genomes and their evolution. However, various protocols are currently used for FISH mapping, most of which are relatively laborious and expensive, or work properly only with a specific type of probes or sequences, and there is a need for a universal and affordable FISH protocol. Here we tested a FISH protocol for mapping of different DNA repeats, such as multigene families (rDNAs, U snDNAs, histone genes), satellite DNAs, microsatellites, transposable elements, DOP-PCR products, and telomeric motif (TTAGG)_n, on the chromosomes of various insects and other arthropods. Different cell types and stages obtained from diverse tissues were used. The FISH procedure proved high quality and reliable results in all experiments performed. We obtained data on the chromosomal distribution of DNA repeats in representatives of insects and other arthropods. Thus, our results allow us to conclude that the protocol is universal and requires only time adjustment for chromosome/DNA denaturation. The use of this FISH protocol will facilitate studies focused on understanding the evolution and role of repetitive DNA in arthropod genomes.

Evolutionary Dynamics of the Repetitive DNA in the Karyotypes of Pipa carvalhoi and Xenopus tropicalis (Anura, Pipidae)

The large amphibian genomes contain numerous repetitive DNA components that have played an important role in the karyotypic diversification of this vertebrate group. Hypotheses based on the presumable primitive karyotype (2n = 20) of the anurans of the family Pipidae suggest that they have evolved principally through intrachromosomal rearrangements. Pipa is the only South American pipid, while all the other genera are found in Africa. The divergence of the South American lineages from the African ones occurred at least 136 million years ago and is thought to have had a strong biogeographic component. Here, we tested the potential of the repetitive DNA to enable a better understanding of the differentiation of the karyotype among the family Pipidae and to expand our capacity to interpret the chromosomal evolution in this frog family. Our results indicate a long history of conservation in the chromosome bearing the H3 histone locus, corroborating inferences on the chromosomal homologies between the species in pairs 6, 8, and 9. The chromosomal distribution of the microsatellite motifs also provides useful markers for comparative genomics at the chromosome level between Pipa carvalhoi and Xenopus tropicalis, contributing new insights into the evolution of the karyotypes of these species. We detected similar patterns in the distribution and abundance of the microsatellite arrangements, which reflect the shared organization in the terminal/subterminal region of the chromosomes between these two species. By contrast, the microsatellite probes detected a differential arrangement of the repetitive DNA among the chromosomes of the two species, allowing longitudinal differentiation of pairs that are identical in size and morphology, such as pairs 1, 2, 4, and 5. We also found evidence of the distinctive composition of the repetitive motifs of the centromeric region between the species analyzed in the present study, with a clear enrichment of the (CA) and (GA) microsatellite motifs in P. carvalhoi. Finally, microsatellite enrichment in the pericentromeric region of chromosome pairs 6, 8, and 9 in the P. carvalhoi karyotype, together with interstitial telomeric sequences (ITS), validate the hypothesis that pericentromeric inversions occurred during the chromosomal evolution of P. carvalhoi and reinforce the role of the repetitive DNA in the remodeling of the karyotype architecture of the Pipidae.

Conserved number of U2 snDNA sites in Piabina argentea, Piabarchus stramineus and two Bryconamericus species (Characidae, Stevardiinae)

ABSTRACT The chromosomal location of 5S rRNA and U2 snRNA genes of Piabina argentea, Piabarchus stramineus and two Bryconamericus species from two different Brazilian river basins were investigated, in order to contribute to the understanding of evolutionary characteristics of these repetitive DNAs in the subfamily Stevardiinae. The diploid chromosome number was 2n = 52 for Bryconamericus cf. iheringii, Bryconamericus turiuba, Piabarchus stramineus and Piabina argentea. The 5S rDNA clusters were located on one chromosome pair in P. stramineus and B. cf. iheringii, and on two pairs in B. turiuba and P. argentea. The U2 snDNA clusters were located on the one pair in all species. Two-color FISH experiments showed that the co-localization between 5S rDNA and U2 snDNA in P. stramineus can represent a marker for this species. Thus, the present study demonstrated that the number of U2 snDNA clusters observed for the four species was conserved, but particular characteristics can be found in the genome of each species.

U1 snDNA chromosomal mapping in ten spittlebug species (Cercopidade, Auchenorrhyncha, Hemiptera)

Spittlebugs, which belong to the family Cercopidae (Auchenorrhyncha, Hemiptera), form a large group of xylem-feeding insects that are best known for causing damage to plantations and pasture grasses. The holocentric chromosomes of these insects remain poorly studied in regards to the organization of different classes of repetitive DNA. To improve chromosomal maps based on repetitive DNAs and to better understand the chromosomal organization and evolutionary dynamics of multigene families in spittlebugs, we physically mapped the U1 snRNA gene with fluorescence in situ hybridization (FISH) in 10 species of Cercopidae belonging to three different genera. All the U1 snDNA clusters were autosomal and located in interstitial position. In seven species, they were restricted to one autosome per haploid genome, while three species of the genus Mahanarva showed two clusters in two different autosomes. Although it was not possible to precisely define the ancestral location of this gene, it was possible to observe the presence of at least one cluster located in a small bivalent in all karyotypes. The karyotype stability observed in Cercopidae is also observed in respect to the distribution of U1 snDNA. Our data are discussed in light of possible mechanisms for U1 snDNA conservation and compared with the available data from other species.

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.

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.

Karyotypes and Repetitive DNA Evolution in Six Species of the Genus Mahanarva (Auchenorrhyncha: Cercopidae)

Insects of the Cercopidae family are widely distributed and comprise 59 genera and 431 species in the New World. They are xylemophagous, causing losses in agricultural and pasture grasses, and are considered as emerging pests. Chromosomally, these insects have been studied by standard techniques, revealing variable diploid numbers and primarily X0 sex chromosome systems (males). We performed chromosome studies in 6 Mahanarva (Cercopidae) species using standard and differential chromosome staining as well as mapping of repetitive DNAs. Moreover, the relationship between the repetitive DNAs was analyzed at the interspecific level. A diploid chromosome number of 2n = 19,X0 was documented, with chromosomes gradually decreasing in size. Neutral or GC-rich regions were detected which varied depending on the species. Fluorescence in situ hybridization with a (TTAGG)n telomeric motif probe revealed terminal signals, matching those of the Cot DNAs obtained from each species, that were also restricted to the terminal regions of all chromosomes. Dot blot analysis with the Cot fraction from M. quadripunctata showed that at least part of the repetitive genome is shared among the 6 species. Our data highlight the conservation of chromosomal features and organization of repetitive DNAs in the genus Mahanarva, suggesting a low differentiation for chromosomes and repetitive DNAs in most of the 6 species studied.

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.

The 5S rDNA in two Abracris grasshoppers (Ommatolampidinae: Acrididae): molecular and chromosomal organization

The 5S ribosomal DNA (rDNA) sequences are subject of dynamic evolution at chromosomal and molecular levels, evolving through concerted and/or birth-and-death fashion. Among grasshoppers, the chromosomal location for this sequence was established for some species, but little molecular information was obtained to infer evolutionary patterns. Here, we integrated data from chromosomal and nucleotide sequence analysis for 5S rDNA in two Abracris species aiming to identify evolutionary dynamics. For both species, two arrays were identified, a larger sequence (named type-I) that consisted of the entire 5S rDNA gene plus NTS (non-transcribed spacer) and a smaller (named type-II) with truncated 5S rDNA gene plus short NTS that was considered a pseudogene. For type-I sequences, the gene corresponding region contained the internal control region and poly-T motif and the NTS presented partial transposable elements. Between the species, nucleotide differences for type-I were noticed, while type-II was identical, suggesting pseudogenization in a common ancestor. At chromosomal point to view, the type-II was placed in one bivalent, while type-I occurred in multiple copies in distinct chromosomes. In Abracris, the evolution of 5S rDNA was apparently influenced by the chromosomal distribution of clusters (single or multiple location), resulting in a mixed mechanism integrating concerted and birth-and-death evolution depending on the unit.

Organization of some repetitive DNAs and B chromosomes in the grasshopper Eumastusia koebelei koebelei (Rehn, 1909) (Orthoptera, Acrididae, Leptysminae)

B chromosomes occur in approximately 15% of eukaryotes and are usually heterochromatic and rich in repetitive DNAs. Here we describe characteristics of a B chromosome in the grasshopper Eumastusia koebelei koebelei (Rehn, 1909) through classical cytogenetic methods and mapping of some repetitive DNAs, including multigene families, telomeric repeats and a DNA fraction enriched with repetitive DNAs obtained from DOP-PCR. Eumastusia koebelei koebelei presented 2n=23, X0 and, in one individual, two copies of the same variant of a B chromosome were noticed, which are associated during meiosis. The C-positive blocks were located in the pericentromeric regions of the standard complement and along the entire length of the B chromosomes. Some G+C-rich heterochromatic blocks were noticed, including conspicuous blocks in the B chromosomes. The mapping of 18S rDNA and U2 snDNA revealed only autosomal clusters, and the telomeric probe hybridized in terminal regions. Finally, the DOP-PCR probe obtained from an individual without a B chromosome revealed signals in the heterochromatic regions, including the entire length of the B chromosome. The possible intraspecific origin of the B chromosomes, due to the shared pool of repetitive DNAs between the A and B chromosomes and the possible consequences of their association are discussed.

Molecular Cytogenetic Analysis of the European Hake Merluccius merluccius (Merlucciidae, Gadiformes): U1 and U2 snRNA Gene Clusters Map to the Same Location

The European hake (Merluccius merluccius) is a highly valuable and intensely fished species in which a long-term alive stock has been established in captivity for aquaculture purposes. Due to their huge economic importance, genetic studies on hakes were mostly focused on phylogenetic and phylogeographic aspects; however chromosome numbers are still not described for any of the fifteen species in the genus Merluccius. In this work we report a chromosome number of 2n = 42 and a karyotype composed of three meta/submetacentric and 18 subtelo/telocentric chromosome pairs. Telomeric sequences appear exclusively at both ends of every single chromosome. Concerning rRNA genes, this species show a single 45S rDNA cluster at an intercalary location on the long arm of subtelocentric chromosome pair 12; the single 5S rDNA cluster is also intercalary to the long arm of chromosome pair 4. While U2 snRNA gene clusters map to a single subcentromeric position on chromosome pair 13, U1 snRNA gene clusters seem to appear on almost all chromosome pairs, but showing bigger clusters on pairs 5, 13, 16, 17 and 19. The brightest signals on pair 13 are coincident with the single U2 snRNA gene cluster signals. Therefore, the use of these probes allows the unequivocal identification of at least 7 of the chromosome pairs that compose the karyotype of Merluccius merluccius thus opening the way to integrate molecular genetics and cytological data on the study of the genome of this important species.

Contrasting the Chromosomal Organization of Repetitive DNAs in Two Gryllidae Crickets with Highly Divergent Karyotypes

A large percentage of eukaryotic genomes consist of repetitive DNA that plays an important role in the organization, size and evolution. In the case of crickets, chromosomal variability has been found using classical cytogenetics, but almost no information concerning the organization of their repetitive DNAs is available. To better understand the chromosomal organization and diversification of repetitive DNAs in crickets, we studied the chromosomes of two Gryllidae species with highly divergent karyotypes, i.e., 2n(♂) = 29,X0 (Gryllus assimilis) and 2n = 9, neo-X1X2Y (Eneoptera surinamensis). The analyses were performed using classical cytogenetic techniques, repetitive DNA mapping and genome-size estimation. Conserved characteristics were observed, such as the occurrence of a small number of clusters of rDNAs and U snDNAs, in contrast to the multiple clusters/dispersal of the H3 histone genes. The positions of U2 snDNA and 18S rDNA are also conserved, being intermingled within the largest autosome. The distribution and base-pair composition of the heterochromatin and repetitive DNA pools of these organisms differed, suggesting reorganization. Although the microsatellite arrays had a similar distribution pattern, being dispersed along entire chromosomes, as has been observed in some grasshopper species, a band-like pattern was also observed in the E. surinamensis chromosomes, putatively due to their amplification and clustering. In addition to these differences, the genome of E. surinamensis is approximately 2.5 times larger than that of G. assimilis, which we hypothesize is due to the amplification of repetitive DNAs. Finally, we discuss the possible involvement of repetitive DNAs in the differentiation of the neo-sex chromosomes of E. surinamensis, as has been reported in other eukaryotic groups. This study provided an opportunity to explore the evolutionary dynamics of repetitive DNAs in two non-model species and will contribute to the understanding of chromosomal evolution in a group about which little chromosomal and genomic information is known.