Genomic organization and comparative chromosome mapping of the U1 snRNA gene in cichlid fish, with an emphasis in Oreochromis niloticus

Cytogenetic Analysis of the Fish Genus Carassius Indicates Divergence, Fission, and Segmental Duplication as Drivers of Tandem Repeat and Microchromosome Evolution

Fishes of the genus Carassius are useful experimental vertebrate models for the study of evolutionary biology and cytogenetics. Carassius demonstrates diverse biological characteristics, such as variation in ploidy levels and chromosome numbers, and presence of microchromosomes. Those Carassius polyploids with ≥150 chromosomes have microchromosomes, but the origin of microchromosomes, especially in European populations, is unknown. We used cytogenetics to study evolution of tandem repeats (U1 and U2 small nuclear DNAs and H3 histone) and microchromosomes in Carassius from the Czech Republic. We tested the hypotheses whether the number of tandem repeats was affected by polyploidization or divergence between species and what mechanism drives evolution of microchromosomes. Tandem repeats were found in tetraploid and hexaploid Carassius gibelio, and tetraploid Carassius auratus and Carassius carassius in conserved numbers, with the exception of U1 small nuclear DNA in C. auratus. This conservation indicates reduction and/or loss in the number of copies per locus in hexaploids and may have occurred by divergence rather than polyploidization. To study the evolution of microchromosomes, we used the whole microchromosome painting probe from hexaploid C. gibelio and hybridized it to tetraploid and hexaploid C. gibelio, and tetraploid C. auratus and C. carassius. Our results revealed variation in the number of microchromosomes in hexaploids and indicated that the evolution of the Carassius karyotype is governed by macrochromosome fissions followed by segmental duplication in pericentromeric areas. These are potential mechanisms responsible for the presence of microchromosomes in Carassius hexaploids. Differential efficacy of one or both of these mechanisms in different tetraploids could ensure variability in chromosome number in polyploids in general.

U2 and U4 snDNA Comparative Chromosomal Mapping in the Neotropical Fish Genera Apareiodon and Parodon (Characiformes: Parodontidae)

Small nuclear DNA (snDNA) are valuable cytogenetic markers for comparative studies in chromosome evolution because different distribution patterns were found among species. Parodontidae, a Neotropical fish family, is known to have female heterogametic sex chromosome systems in some species. The U2 and U4 snDNA sites have been found to be involved in Z and W chromosome differentiation in Apareiodon sp., Apareiodon affinis, and Parodon hilarii. However, few studies have evaluated snDNA sites as propulsors of chromosome diversification among closely related fish species. In this study, we investigated the distribution of U2 and U4 snDNA clusters in the chromosomes of 10 populations/species belonging to Apareiodon and Parodon, aiming to identify chromosomal homeologies or diversification. In situ localization data revealed a submetacentric pair carrying the U2 snDNA site among the populations/species analyzed. Furthermore, all studied species demonstrated homeology in the location of U4 snDNA cluster in the proximal region of metacentric pair 1, besides an additional signal showing up with a divergence in Apareiodon. Comparative chromosomal mapping of U4 snDNA also helped to reinforce the proposal of the ZZ/ZW1W2 sex chromosome system origin in an A. affinis population. According to cytogenetic data, the study corroborates the diversification in Parodontidae paired species with uncertain taxonomy.

Mechanisms of Karyotypic Diversification in Ancistrus (Siluriformes, Loricariidae): Inferences from Repetitive Sequence Analysis

Ancistrus is a highly diverse neotropical fish genus that exhibits extensive chromosomal variability, encompassing karyotypic morphology, diploid chromosome number (2n = 34-54), and the evolution of various types of sex chromosome systems. Robertsonian rearrangements related to unstable chromosomal sites are here described. Here, the karyotypes of two Ancistrus species were comparatively analyzed using classical cytogenetic techniques, in addition to isolation, cloning, sequencing, molecular characterization, and fluorescence in situ hybridization of repetitive sequences (i.e., 18S and 5S rDNA; U1, U2, and U5 snDNA; and telomere sequences). The species analyzed here have different karyotypes: Ancistrus sp. 1 (2n = 38, XX/XY) and Ancistrus cirrhosus (2n = 34, no heteromorphic sex chromosomes). Comparative mapping showed different organizations for the analyzed repetitive sequences: 18S and U1 sequences occurred in a single site in all populations of the analyzed species, while 5S and U2 sequences could occur in single or multiple sites. A sequencing analysis confirmed the identities of the U1, U2, and U5 snDNA sequences. Additionally, a syntenic condition for U2-U5 snDNA was found in Ancistrus. In a comparative analysis, the sequences of rDNA and U snDNA showed inter- and intraspecific chromosomal diversification. The occurrence of Robertsonian rearrangements and other dispersal mechanisms of repetitive sequences are discussed.

Chromosomal Mapping of the Histone Multigene Family and U2 snRNA in Hypancistrus Species (Siluriformes, Loricariidae) from the Brazilian Amazon

Loricariidae (Siluriformes) comprises ∼1026 species of neotropical fish, being considered the most diverse among the Siluriformes. Studies on repetitive DNA sequences have provided important data on the evolution of the genomes of members of this family, especially of the Hypostominae subfamily. In this study, the chromosomal mapping of the histone multigene family and U2 snRNA was performed in two species belonging to the Hypancistrus genus, Hypancistrus sp. "pão" (2n = 52, 22m + 18sm +12st) and Hypancistrus zebra (2n = 52, 16m + 20sm +16st). The presence of dispersed signals of histones H2A, H2B, H3, and H4 in the karyotype of both species, with each sequence displaying a varied level of accumulation and dispersion of these sequences between them was observed; in addition, U2 snDNA probe only showed positive results in H. zebra, which present this multigene in the terminal region of three chromosomal pairs. The obtained results resemble data already analyzed in the literature, in which the action of transposable elements interfere in the organization of these multigene families, in addition to other evolutionary processes that shape the evolution of the genome, such as circular or ectopic recombination. This study also shows that the dispersion of the multigene histone family is quite complex, and from this, these data serve as a point of discussion for the evolutionary processes that occur in the Hypancistrus karyotype.

Chromosomal Diversification in Pseudacanthicus Species (Loricariidae, Hypostominae) Revealed by Comparative Mapping of Repetitive Sequences

Simple Summary

The fishes of the Loricariidae family have a huge genetic diversity, mainly involving variations in the number and shape of chromosomes. The recognition of the species genus Pseudacanthicus is complex due to the large diversity of forms and limited knowledge of their genetic diversity. In this study, the karyotypes of three Pseudacanthicus species were comparatively analyzed using classical and molecular methods. They presented the same diploid number, but with different compositions of repetitive DNA sequences. Such information can be useful for the recognition of distinct species, in addition to providing important insights into the real biodiversity of this important group of Neotropical fish.

Abstract

Pseudacanthicus is a genus of Neotropical fish with eight valid species, in addition to numerous lineages not formally identified. It occurs along the Amazon and Tocantins River basins, in Suriname and in the Guiana shield. There are no karyotypic data in the literature for species of this genus. Here, the karyotypes of three Pseudacanthicus species (P. spinosus, P. leopardus and Pseudacanthicus sp.) were comparatively analyzed by classical cytogenetics and fluorescence in situ hybridization using 18S and 5S rDNA probes, U2 snDNA and telomeric sequences. The analyzed species presented 52 chromosomes and KF = 18 m + 34 sm. Constitutive heterochromatin occurred in blocks on a few chromosomes. The 18S rDNA occurred in a single pair; interestingly, P. leopardus presented only one locus of this sequence in its diploid genome. The 5S rDNA sequence occurred in only one pair in P. leopardus, and in multiple sites in Pseudacanthicus sp. and P. spinosus. The snDNA U2 occurred in only one pair in all analyzed species. Telomeric sequences did not show interstitial sites. Although Pseudacanthicus species share the same 2n and KF, repetitive sequence analysis revealed karyotypic diversity among these species. The occurrence of DNA double-strand breaks related to fragile sites, unequal crossing over and transpositions is proposed as the mechanism of karyotypic diversification, suggesting that the conservation of the karyotypic macrostructure is only apparent in this group of fish.

Comparative Cytogenetic and Sequence Analysis of U Small Nuclear RNA Genes in Three Ancistrus Species (Siluriformes: Loricariidae)

Ancistrus presents a wide karyotypic diversity, resulting from numeric and structural chromosomal rearrangements. It has been proposed that some genome-specific regions containing repetitive units could organize prone-to-break DNA sites in Loricariidae, triggering chromosomal rearrangements such as Robertsonian fusions (Rb fusions), centric fissions, translocations, and inversions. The tandemly repeats of the small nuclear RNAs (snRNAs) gene families are considered good cytogenetic markers for understanding chromosomal remodeling events among closely related species, but these snRNAs have been scarcely analyzed in Ancistrus. This study presented the nucleotide sequencing and comparative in situ location of U snRNA sequences from Ancistrus aguaboensis, Ancistrus cf. multispinis, and Ancistrus sp. (2n = 50, 52, and 50, respectively), aiming to provide information about snRNA clusters in the genome and chromosome evolution in Ancistrus. U snRNA nucleotide sequences of Ancistrus presented identity to orthologous copies and folded their secondary structures correctly. In situ localization and karyotyping of the three Ancistrus species revealed clustered copies of U2 and U5 snRNA gene families to a single chromosome site, one chromosome pair bearing U1 snRNA sequence, and one main locus of U4 snRNA sequence, besides scattered signals along the chromosomes. Previous studies related the participation of the rRNA gene families in centric fusion events, contributing to chromosome rearrangements and karyotype plasticity present in Loricariidae. In this study, homeologies in U snRNA loci chromosomal locations were detected, indicating the occurrence of conserved sites of these gene families in these three Ancistrus species with 2n = 50 or 52 chromosomes.

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.

Comparative Cytogenetics Analysis Among Peckoltia Species (Siluriformes, Loricariidae): Insights on Karyotype Evolution and Biogeography in the Amazon Region

Peckoltia is widely distributed genus in the Amazon and Orinoco basins and the Guiana Shield, containing 18 valid species, and distinct morphotypes still needing description in the scientific literature due to its great taxonomic complexity. This study performed a comparative chromosomal analysis of two undescribed Peckoltia species (Peckoltia sp. 3 Jarumã and Peckoltia sp. 4 Caripetuba) from the Brazilian Amazon using conventional chromosome bands methods and in situ localization of the repetitive DNA (5S and 18S rRNA and U1 snRNA genes and telomeric sequences). Both species presented 2n = 52 but differed in their karyotype formula, probably due to inversions or translocations. The nucleolus organizer regions (NORs) showed distal location on a probably homeologous submetacentric pair in both species, besides an extra signal in a subtelocentric chromosome in Peckoltia sp. 4 Caripetuba. Heterochromatin occurred in large blocks, with different distributions in the species. The mapping of the 18S and 5S rDNA, and U1 snDNA showed differences in locations and number of sites. No interstitial telomeric sites were detected using the (TTAGGG)n probes. Despite 2n conservationism in Peckoltia species, the results showed variation in karyotype formulas, chromosomal bands, and locations of repetitive sites, demonstrating great chromosomal diversity. A proposal for Peckoltia karyotype evolution was inferred in this study based on the diversity of location and number of chromosomal markers analyzed. A comparative analysis with other Peckoltia karyotypes described in the literature, their biogeography patterns, and molecular phylogeny led to the hypothesis that the derived karyotype was raised in the left bank of the Amazon River.

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.

Molecular cytogenetic analyses reveal extensive chromosomal rearrangements and novel B chromosomes in Moenkhausia (Teleostei, Characidae)

The cytogenetic characteristics of five fish species of the Moenkhausia are described, based on the analysis of specimens collected in different headwater. All the species analyzed presented 2n=50 chromosomes. The C-banding revealed a similar distribution pattern of heterochromatic blocks in all the species, except Moenkhausia nigromarginata. The 5S rDNA sites were distributed on multiple chromosome pairs in all five species. Single and multiple histone H1 sites were observed in all the species, and histone H1 was shown to be co-located with the 18S rRNA gene in a single chromosome pair. The U2 snDNA gene was distributed at multiple sites in all the Moenkhausia species. The presence of B microchromosomes was confirmed in Moenkhausia forestii, while individuals of the three study populations of Moenkhausia oligolepis presented three morphologically distinct types of B chromosome. The chromosomal mapping of the 18S rDNA sites using the FISH technique revealed signals in the B chromosomes of M. forestii, while clusters of the H1 histone and U2 snDNA genes were found in the B chromosomes of M. forestii and M. oligolepis. The classical and molecular cytogenetic markers used in this study revealed ample variation in the Moenkhausia karyotypes, reflecting the dynamic nature of the chromosomal evolution.

Karyotypic Diversification in Two Megaleporinus Species (Characiformes, Anostomidae) Inferred from In Situ Localization of Repetitive DNA Sequences

Anostomidae species have conserved diploid numbers (2n = 54), although comparative cytogenetic studies have demonstrated chromosomal rearrangements occurrence among them, especially in repetitive DNA rich regions. The location and distribution of ribosomal DNA (rDNA) and small nuclear RNAs (snRNAs) multigene families are highly dynamic in the genomes of several organisms. In this study, we in situ located the rDNA and snRNA sites in two populations of Megaleporinus obtusidens and a sample of Megaleporinus reinhardti to infer their chromosomal changes in the evolutionary lineages. Both species of Megaleporinus shared 2n = 54 chromosomes with the presence of ZZ/ZW sex chromosome system, but they diverged in relationship to the location of 5S and 45S rDNAs as well as the distribution of snRNAs sites. The characterization of the analyzed sequences revealed the presence of complete rDNA and snRNAs sequences as well as snRNAs containing transposable elements (TEs) and microsatellite repeats. After chromosomal mapping, the sequences encompassing TEs proved to be dispersed through autosomes and accumulated on sex chromosomes. The data demonstrate that intra- and interspecific chromosomal changes occurred involving the multigene family's sites in Megaleporinus karyotypes. Furthermore, we detected TE-like sequences in the differentiation of sex chromosome systems in M. obtusidens and M. reinhardti.

Molecular Analysis and Chromosome Mapping of Repetitive DNAs in the Green Terror Andinoacara rivulatus (Cichlidae: Cichlasomatini)

Neotropical cichlids include hundreds of species whose taxonomy has benefited of molecular phylogeny and whose karyotype evolution has been related to the amount and distribution of different classes of repetitive sequences. This study provides the first integrative molecular (cytochrome c oxidase subunit 1 and 16S sequences) and cytogenetic analyses of wild samples of the green terror Andinoacara rivulatus, a cichlid naturally distributed in Ecuador and spread throughout the world as an aquarium pet. Molecular data revealed that sequences of green terror constitute a single monophyletic clade within the genus and allowed species attribution of uncertain samples previously cytogenetically analyzed. Chromosome number (2n = 48) conforms to the general trend observed within neotropical cichlids. However, mapping of different classes of repeated sequences (18S rDNA, 5S rDNA, U1 snDNA and telomeric) revealed the presence of features uncommon among representatives of these fishes, like multiple major rDNA sites, and suggested a recent occurrence of rearrangements (fusion/inversion) in two chromosome pairs.

Taxonomic Diversity Not Associated with Gross Karyotype Differentiation: The Case of Bighead Carps, Genus Hypophthalmichthys (Teleostei, Cypriniformes, Xenocyprididae)

The bighead carps of the genus Hypophthalmichthys (H. molitrix and H. nobilis) are important aquaculture species. They were subjected to extensive multidisciplinary research, but with cytogenetics confined to conventional protocols only. Here, we employed Giemsa-/C-/CMA3- stainings and chromosomal mapping of multigene families and telomeric repeats. Both species shared (i) a diploid chromosome number 2n = 48 and the karyotype structure, (ii) low amount of constitutive heterochromatin, (iii) the absence of interstitial telomeric sites (ITSs), (iv) a single pair of 5S rDNA loci adjacent to one major rDNA cluster, and (v) a single pair of co-localized U1/U2 snDNA tandem repeats. Both species, on the other hand, differed in (i) the presence/absence of remarkable interstitial block of constitutive heterochromatin on the largest acrocentric pair 11 and (ii) the number of major (CMA3-positive) rDNA sites. Additionally, we applied here, for the first time, the conventional cytogenetics in H. harmandi, a species considered extinct in the wild and/or extensively cross-hybridized with H. molitrix. Its 2n and karyotype description match those found in the previous two species, while silver staining showed differences in distribution of major rDNA. The bighead carps thus represent another case of taxonomic diversity not associated with gross karyotype differentiation, where 2n and karyotype structure cannot help in distinguishing between genomes of closely related species. On the other hand, we demonstrated that two cytogenetic characters (distribution of constitutive heterochromatin and major rDNA) may be useful for diagnosis of pure species. The universality of these markers must be further verified by analyzing other pure populations of bighead carps.

De novo genome assembly of the cichlid fish Astatotilapia latifasciata reveals a higher level of genomic polymorphism and genes related to B chromosomes

Supernumerary B chromosomes (Bs) are accessory elements to the regular chromosome set (As) and have been observed in a huge diversity of eukaryotic species. Although extensively investigated, the biological significance of Bs remains enigmatic. Here, we present de novo genome assemblies for the cichlid fish Astatotilapia latifasciata, a well-known model to study Bs. High coverage data with Illumina sequencing was obtained for males and females with 0B (B-), 1B, and 2B (B+) chromosomes to provide information regarding the diversity among these genomes. The draft assemblies comprised 771 Mb for the B- genome and 781 Mb for the B+ genome. Comparative analysis of the B+ and B- assemblies reveals syntenic discontinuity, duplicated blocks and several insertions, deletions, and inversions indicative of rearrangements in the B+ genome. Hundreds of transposable elements and 1546 protein coding sequences were annotated in the duplicated B+ regions. Our work contributes a list of thousands of genes harbored on the B chromosome, with functions in several biological processes, including the cell cycle.

Identification of a New Mullet Species Complex Based on an Integrative Molecular and Cytogenetic Investigation of Mugil hospes (Mugilidae: Mugiliformes)

Mullets are very common fishes included in the family Mugilidae, (Mugiliformes), which are characterized by both a remarkably uniform external morphology and internal anatomy. Recently, within this family, different species complexes were molecularly identified within Mugil, a genus which is characterized by lineages that sometimes show very different karyotypes. Here we report the results of cytogenetic and molecular analyses conducted on Mugil hospes, commonly known as the hospe mullet, from Ecuador. The study aims to verify whether the original described species from the Pacific Ocean corresponds to that identified in the Atlantic Ocean, and to identify species-specific chromosome markers that can add new comparative data about Mugilidae karyotype evolution. The karyotype of M. hospes from Ecuador is composed of 48 acrocentric chromosomes and shows two active nucleolar organizer regions (NORs). In situ hybridization, using different types of repetitive sequences (rDNAs, U1 snDNA, telomeric repeats) as probes, identified species-specific chromosome markers that have been compared with those of other species of the genus Mugil. Cytochrome c oxidase subunit I (COI) sequence analysis shows only 92-93% similarity with sequences previously deposited under this species name in GenBank, all of which were from the Atlantic Ocean. Phylogenetic reconstructions indicate the presence of three well-supported hospe mullet lineages whose molecular divergence is compatible with the presence of distinct species. Indeed, the first lineage includes samples from Ecuador, whereas the other two lineages include the Atlantic samples and correspond to M. brevirostris from Brazil and Mugil sp. R from Belize/Venezuela. Results here provided reiterate the pivotal importance of an integrative molecular and cytogenetic approach in the reconstruction of the relationships within Mugilidae.

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.

Development of chromosomal markers based on next-generation sequencing: the B chromosome of the cichlid fish Astatotilapia latifasciata as a model

BACKGROUND

B chromosomes (Bs) are additional chromosomal elements found in a wide range of eukaryotes including fungi, plants and animals. B chromosomes are still enigmatic despite being the subject of hundreds, even thousands of reports. As yet there is no comprehensive theory for the biological role of B chromsomes thus, new studies are needed. Next-generation sequencing (NGS) holds promise for investigating classical issues in chromosome biology. NGS uses a large-scale approach that is required for advancing classical cytogenetic studies. Based on 454 sequencing data of a microdissected B chromosome and Illumina whole-genome sequencing data generated for 0B, 1B and 2B animals, we developed PCR- and qPCR-based markers for the B chromosomes of the cichlid fish Astatotilapia latifasciata (that possess 0, 1 or 2 B chromosomes).

RESULTS

Specific PCR primers were designed to produce two amplified fragments for B-positive samples and the control fragment for B-negative samples. Thus, PCR markers detected the presence/absence of Bs but did not provide information about the number of Bs. However, quantitative PCR (qPCR) markers clearly discriminated between 1B and 2B samples. The high copy number of the marker identified in the B chromosomes was confirmed by chromosome mapping.

CONCLUSIONS

The analysis of chromosome polymorphisms based on a NGS approach is a powerful strategy to obtain markers that detect the presence/absence of extra chromosomes or the gain or loss of genomic blocks. Further, qPCR can also provide information regarding the relative copy number of specific DNA fragments. These methods are useful to investigate various chromosome polymorphisms, including B and sex chromosomes, as well as chromosomal duplications and deletions. NGS data provide a detailed analysis of the composition of genomic regions that are thought to be present in B chromosomes.

The repetitive DNA element BncDNA, enriched in the B chromosome of the cichlid fish Astatotilapia latifasciata, transcribes a potentially noncoding RNA

Supernumerary chromosomes have been studied in many species of eukaryotes, including the cichlid fish, Astatotilapia latifasciata. However, there are many unanswered questions about the maintenance, inheritance, and functional aspects of supernumerary chromosomes. The cichlid family has been highlighted as a model for evolutionary studies, including those that focus on mechanisms of chromosome evolution. Individuals of A. latifasciata are known to carry up to two B heterochromatic isochromosomes that are enriched in repetitive DNA and contain few intact gene sequences. We isolated and characterized a transcriptionally active repeated DNA, called B chromosome noncoding DNA (BncDNA), highly represented across all B chromosomes of A. latifasciata. BncDNA transcripts are differentially processed among six different tissues, including the production of smaller transcripts, indicating transcriptional variation may be linked to B chromosome presence and sexual phenotype. The transcript lengths and lack of similarity with known protein/gene sequences indicate BncRNA might represent a novel long noncoding RNA family (lncRNA). The potential for interaction between BncRNA and known miRNAs were computationally predicted, resulting in the identification of possible binding of this sequence in upregulated miRNAs related to the presence of B chromosomes. In conclusion, Bnc is a transcriptionally active repetitive DNA enriched in B chromosomes with potential action over B chromosome maintenance in somatic cells and meiotic drive in gametic cells.

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.

Non-Random Distribution of 5S rDNA Sites and Its Association with 45S rDNA in Plant Chromosomes

5S and 45S rDNA sites are the best mapped chromosome regions in eukaryotic chromosomes. In this work, a database was built gathering information about the position and number of 5S rDNA sites in 784 plant species, aiming to identify patterns of distribution along the chromosomes and its correlation with the position of 45S rDNA sites. Data revealed that in most karyotypes (54.5%, including polyploids) two 5S rDNA sites (a single pair) are present, with 58.7% of all sites occurring in the short arm, mainly in the proximal region. In karyotypes of angiosperms with only 1 pair of sites (single sites) they are mostly found in the proximal region (52.0%), whereas in karyotypes with multiple sites the location varies according to the average chromosome size. Karyotypes with multiple sites and small chromosomes (<3 µm) often display proximal sites, while medium-sized (between 3 and 6 µm) and large chromosomes (>6 µm) more commonly show terminal or interstitial sites. In species with holokinetic chromosomes, the modal value of sites per karyotype was also 2, but they were found mainly in a terminal position. Adjacent 5S and 45S rDNA sites were often found in the short arm, reflecting the preferential distribution of both sites in this arm. The high frequency of genera with at least 1 species with adjacent 5S and 45S sites reveals that this association appeared several times during angiosperm evolution, but it has been maintained only rarely as the dominant array in plant genera.

Chromosomal Mapping of Repetitive DNA Sequences in Five Species of Astyanax (Characiformes, Characidae) Reveals Independent Location of U1 and U2 snRNA Sites and Association of U1 snRNA and 5S rDNA

Astyanax is a genus of Characidae fishes currently composed of 155 valid species. Previous cytogenetic studies revealed high chromosomal diversification among them, and several studies have been performed using traditional cytogenetic techniques to investigate karyotypes and chromosomal locations of 18S and 5S rDNA genes. However, only a few studies are currently available about other repetitive sequences. Here, the chromosomal location of small nuclear RNA genes, identified as U1 and U2 snRNA clusters, was established and compared to the distribution of 5S rDNA and histone clusters in 5 Astyanax species (A. paranae, A. fasciatus, A. bockmanni, A. altiparanae, and A. jordani) using FISH. The cytogenetic mapping of U1 and U2 snRNA demonstrated a conserved pattern in the number of sites per genome independent of the location in Astyanax species. The location of the U1 snRNA gene was frequently associated with 5S rDNA sequences, indicating a possible interaction between the distinct repetitive DNA families. Finally, comparisons involving the location of U1 and U2 snRNA clusters in the chromosomes of Astyanax species revealed a very diverse pattern, suggesting that many rearrangements have occurred during the diversification process of this group.

High similarity of U2 snDNA sequence between A and B chromosomes in the grasshopper Abracris flavolineata

B chromosomes are frequently enriched for a wide variety of repetitive DNAs. Among grasshoppers in the species Abracris flavolineata (Ommatolampidinae) the B chromosomes are submetacentric, C-negative and harbor repetitive DNAs such as, U2 snDNA, C 0 t-1 DNA, two Mariner-like elements and some microsatellites. Here, we provide evidence showing the intragenome similarity between the B chromosome and the A complement in A. flavolineata, combining analysis of microdissection and chromosome painting and B chromosome-specific amplification through polymerase chain reaction (PCR) of U2 snDNA. Chromosome painting revealed signals spread through the C-negative regions, including the A and B chromosomes. Moreover, significant clustered signals forming bands were observed in some A chromosomes, and for the B chromosome, significant signals were located on both arms, which could be caused by accumulation of repetitive DNA sequences. The C-positive regions did not reveal any signals. Sequence comparison of U2 snDNA between that obtained from a genome without the B chromosome and that from µB-DNA revealed high similarity with the occurrence of four shared haplotypes, one of them (i.e., Hap1) being highly prevalent and putatively ancestral. The highest divergence from Hap1 was observed for Hap3, which was caused by only six mutational steps. These data support an intraspecific origin of the B chromosome in A. flavolineata that is highly similar with the A complement, and the low U2 snDNA sequence diversity observed in the B chromosome could be related to its recent origin, besides intrachromosomal concerted evolution for U2 snDNA repeats in the B chromosome.

Chromosomal distribution of microsatellite repeats in Amazon cichlids genome (Pisces, Cichlidae)

Fish of the family Cichlidae are recognized as an excellent model for evolutionary studies because of their morphological and behavioral adaptations to a wide diversity of explored ecological niches. In addition, the family has a dynamic genome with variable structure, composition and karyotype organization. Microsatellites represent the most dynamic genomic component and a better understanding of their organization may help clarify the role of repetitive DNA elements in the mechanisms of chromosomal evolution. Thus, in this study, microsatellite sequences were mapped in the chromosomes of Cichla monoculus Agassiz, 1831, Pterophyllum scalare Schultze, 1823, and Symphysodon discus Heckel, 1840. Four microsatellites demonstrated positive results in the genome of Cichla monoculus and Symphysodon discus, and five demonstrated positive results in the genome of Pterophyllum scalare. In most cases, the microsatellite was dispersed in the chromosome with conspicuous markings in the centromeric or telomeric regions, which suggests that sequences contribute to chromosome structure and may have played a role in the evolution of this fish family. The comparative genome mapping data presented here provide novel information on the structure and organization of the repetitive DNA region of the cichlid genome and contribute to a better understanding of this fish family's genome.

A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs

The desert locust (Schistocerca gregaria) has been used as material for numerous cytogenetic studies. Its genome size is estimated to be 8.55 Gb of DNA comprised in 11 autosomes and the X chromosome. Its X0/XX sex chromosome determinism therefore results in females having 24 chromosomes whereas males have 23. Surprisingly, little is known about the DNA content of this locust's huge chromosomes. Here, we use the Feulgen Image Analysis Densitometry and C-banding techniques to respectively estimate the DNA quantity and heterochromatin content of each chromosome. We also identify three satellite DNAs using both restriction endonucleases and next-generation sequencing. We then use fluorescent in situ hybridization to determine the chromosomal location of these satellite DNAs as well as that of six tandem repeat DNA gene families. The combination of the results obtained in this work allows distinguishing between the different chromosomes not only by size, but also by the kind of repetitive DNAs that they contain. The recent publication of the draft genome of the migratory locust (Locusta migratoria), the largest animal genome hitherto sequenced, invites for sequencing even larger genomes. S. gregaria is a pest that causes high economic losses. It is thus among the primary candidates for genome sequencing. But this species genome is about 50 % larger than that of L. migratoria, and although next-generation sequencing currently allows sequencing large genomes, sequencing it would mean a greater challenge. The chromosome sizes and markers provided here should not only help planning the sequencing project and guide the assembly but would also facilitate assigning assembled linkage groups to actual chromosomes.

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.

Sequencing and characterization of the transcriptome of half-smooth tongue sole (Cynoglossus semilaevis)

BACKGROUND

Half-smooth tongue sole (Cynoglossus semilaevis) is a valuable fish for aquaculture in China. This fish exhibits sexual dimorphism, particularly different growth rates and body sizes between two genders. Thus, C. semilaevis is a good model that can be used to investigate mechanisms responsible for such dimorphism, this model can also be utilized to answer fundamental questions in evolution and applied fields of aquaculture. Hence, advances in second-generation sequencing technology, such as 454 pyrosequencing, could provide a robust tool to study the genome characteristics of non-model species.

RESULTS

In this study, C. semilaevis was subjected to de novo transcriptome sequencing and characterization. A total of 749,954 reads were generated using a single 454 sequencing run in a full PicoTiter plate. These reads were then assembled into 62,632 contigs with a 10-fold average sequencing coverage. A total of 26,589 sequences were successfully annotated based on sequence similarities; among these sequences, 3,451 transcripts exhibited gene ontology terms and 2,362 showed enzyme commissions associated with 186 pathways from Kyoto Encyclopedia of Gene and Genomes pathways. A search of repetitive elements was performed, and 1,898 transposable elements were identified. Approximately 7,800 simple-sequence repeats and 21,234 single-nucleotide polymorphisms were also detected.

CONCLUSIONS

Our data provided an integrated and comprehensive transcriptome resource for C. semilaevis. These data could be used for further research in population genetics, gene function, and tissue-specific gene expressions.

Tracking the evolution of sex chromosome systems in Melanoplinae grasshoppers through chromosomal mapping of repetitive DNA sequences

BACKGROUND

The accumulation of repetitive DNA during sex chromosome differentiation is a common feature of many eukaryotes and becomes more evident after recombination has been restricted or abolished. The accumulated repetitive sequences include multigene families, microsatellites, satellite DNAs and mobile elements, all of which are important for the structural remodeling of heterochromatin. In grasshoppers, derived sex chromosome systems, such as neo-XY♂/XX♀ and neo-X1X2Y♂/X1X1X2X2♀, are frequently observed in the Melanoplinae subfamily. However, no studies concerning the evolution of sex chromosomes in Melanoplinae have addressed the role of the repetitive DNA sequences. To further investigate the evolution of sex chromosomes in grasshoppers, we used classical cytogenetic and FISH analyses to examine the repetitive DNA sequences in six phylogenetically related Melanoplinae species with X0♂/XX♀, neo-XY♂/XX♀ and neo-X1X2Y♂/X1X1X2X2♀ sex chromosome systems.

RESULTS

Our data indicate a non-spreading of heterochromatic blocks and pool of repetitive DNAs (C0t-1 DNA) in the sex chromosomes; however, the spreading of multigene families among the neo-sex chromosomes of Eurotettix and Dichromatos was remarkable, particularly for 5S rDNA. In autosomes, FISH mapping of multigene families revealed distinct patterns of chromosomal organization at the intra- and intergenomic levels.

CONCLUSIONS

These results suggest a common origin and subsequent differential accumulation of repetitive DNAs in the sex chromosomes of Dichromatos and an independent origin of the sex chromosomes of the neo-XY and neo-X1X2Y systems. Our data indicate a possible role for repetitive DNAs in the diversification of sex chromosome systems in grasshoppers.

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.

Evolutionary dynamics of rRNA gene clusters in cichlid fish

Background

Among multigene families, ribosomal RNA (rRNA) genes are the most frequently studied and have been explored as cytogenetic markers to study the evolutionary history of karyotypes among animals and plants. In this report, we applied cytogenetic and genomic methods to investigate the organization of rRNA genes among cichlid fishes. Cichlids are a group of fishes that are of increasing scientific interest due to their rapid and convergent adaptive radiation, which has led to extensive ecological diversity.

Results

The present paper reports the cytogenetic mapping of the 5S rRNA genes from 18 South American, 22 African and one Asian species and the 18S rRNA genes from 3 African species. The data obtained were comparatively analyzed with previously published information related to the mapping of rRNA genes in cichlids. The number of 5S rRNA clusters per diploid genome ranged from 2 to 15, with the most common pattern being the presence of 2 chromosomes bearing a 5S rDNA cluster. Regarding 18S rDNA mapping, the number of sites ranged from 2 to 6, with the most common pattern being the presence of 2 sites per diploid genome. Furthermore, searching the Oreochromis niloticus genome database led to the identification of a total of 59 copies of 5S rRNA and 38 copies of 18S rRNA genes that were distributed in several genomic scaffolds. The rRNA genes were frequently flanked by transposable elements (TEs) and spread throughout the genome, complementing the FISH analysis that detect only clustered copies of rRNA genes.

Conclusions

The organization of rRNA gene clusters seems to reflect their intense and particular evolutionary pathway and not the evolutionary history of the associated taxa. The possible role of TEs as one source of rRNA gene movement, that could generates the spreading of ribosomal clusters/copies, is discussed. The present paper reinforces the notion that the integration of cytogenetic data and genomic analysis provides a more complete picture for understanding the organization of repeated sequences in the genome.

Genetic characterization of Plectorhinchus mediterraneus yields important clues about genome organization and evolution of multigene families

Background

Molecular and cytogenetic markers are of great use for to fish characterization, identification, phylogenetics and evolution. Multigene families have proven to be good markers for a better understanding of the variability, organization and evolution of fish species. Three different tandemly-repeated gene families (45S rDNA, 5S rDNA and U2 snDNA) have been studied in Plectorhinchus mediterraneus (Teleostei: Haemulidae), at both molecular and cytogenetic level, to elucidate the taxonomy and evolution of these multigene families, as well as for comparative purposes with other species of the family.

Results

Four different types of 5S rDNA were obtained; two of them showed a high homology with that of Raja asterias, and the putative implication of a horizontal transfer event and its consequences for the organization and evolution of the 5S rDNA have been discussed. The other two types do not resemble any other species, but in one of them a putative tRNA-derived SINE was observed for the first time, which could have implications in the evolution of the 5S rDNA. The ITS-1 sequence was more related to a species of another different genus than to that of the same genus, therefore a revision of the Hamulidae family systematic has been proposed. In the analysis of the U2 snDNA, we were able to corroborate that U2 snDNA and U5 snDNA were linked in the same tandem array, and this has interest for tracing evolutionary lines. The karyotype of the species was composed of 2n = 48 acrocentric chromosomes, and each of the three multigene families were located in different chromosome pairs, thus providing three different chromosomal markers.

Conclusions

Novel data can be extracted from the results: a putative event of horizontal transfer, a possible tRNA-derived SINE linked to one of the four 5S rDNA types characterized, and a linkage between U2 and U5 snDNA. In addition, a revision of the taxonomy of the Haemulidae family has been suggested, and three cytogenetic markers have been obtained. Some of these results have not been described before in any other fish species. New clues about the genome organization and evolution of the multigene families are offered in this study.