Evidence for 5S rDNA horizontal transfer in the toadfish Halobatrachus didactylus (Schneider, 1801) based on the analysis of three multigene families

Long-read de novo genome assembly of Gulf toadfish (Opsanus beta)

BACKGROUND

The family Batrachoididae are a group of ecologically important teleost fishes with unique life histories, behavior, and physiology that has made them popular model organisms. Batrachoididae remain understudied in the realm of genomics, with only four reference genome assemblies available for the family, with three being highly fragmented and not up to current assembly standards. Among these is the Gulf toadfish, Opsanus beta, a model organism for serotonin physiology which has recently been bred in captivity.

RESULTS

Here we present a new, de novo genome and transcriptome assemblies for the Gulf toadfish using PacBio long read technology. The genome size of the final assembly is 2.1 gigabases, which is among the largest teleost genomes. This new assembly improves significantly upon the currently available reference for Opsanus beta with a final scaffold count of 62, of which 23 are chromosome scale, an N50 of 98,402,768, and a BUSCO completeness score of 97.3%. Annotation with ab initio and transcriptome-based methods generated 41,076 gene models. The genome is highly repetitive, with ~ 70% of the genome composed of simple repeats and transposable elements. Satellite DNA analysis identified potential telomeric and centromeric regions.

CONCLUSIONS

This improved assembly represents a valuable resource for future research using this important model organism and to teleost genomics more broadly.

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.

Chromosomal mapping of repetitive sequences in Hyphessobrycon eques (Characiformes, Characidae): a special case of the spreading of 5S rDNA clusters in a genome

Cytogenetic data showed a variation in diploid chromosome number in the genus Hyphessobrycon ranging from 2n = 46 to 52, and studies involving repetitive DNA sequences are scarce in representatives of this genus. The purpose of this paper was the chromosomal mapping of repetitive sequences (rDNA, histone genes, U snDNA and microsatellites) and investigation of the amplification of 5S rDNA clusters in the Hyphessobrycon eques genome. Two H. eques populations displayed 2n = 52 chromosomes, with the acrocentric pair No. 24 bearing Ag-NORs corresponding with CMA₃+/DAPI-. FISH with a 18S rDNA probe identified the NORs on the short (p) arms of the acrocentric pairs Nos. 22 and 24. The 5S rDNA probe visualized signals on almost all chromosomes in genomes of individuals from both populations (40 signals); FISH with H3 histone probe identified two chromosome pairs, with the pericentromeric location of signals; FISH with a U2 snDNA probe identified one chromosome pair bearing signals, on the interstitial chromosomal region. The mononucleotide (A), dinucleotide (CA) and tetranucleotide (GATA) repeats were observed on the centromeric/pericentromeric and/or terminal positions of all chromosomes, while the trinucleotide (CAG) repeat showed signals on few chromosomes. Molecular analysis of 5S rDNA and non-transcribed spacers (NTS) showed microsatellites (GATA and A repeats) and a fragment of retrotransposon (SINE3/5S-Sauria) inside the sequences. This study expanded the available cytogenetic data for H. eques and demonstrated to the dispersion of the 5S rDNA sequences on almost all chromosomes.

Integrative genetic map of repetitive DNA in the sole Solea senegalensis genome shows a Rex transposon located in a proto-sex chromosome

Repetitive sequences play an essential role in the structural and functional evolution of the genome, particularly in the sexual chromosomes. The Senegalese sole (Solea senegalensis) is a valuable flatfish in aquaculture albeit few studies have addressed the mapping and characterization of repetitive DNA families. Here we analyzed the Simple Sequence Repeats (SSRs) and Transposable elements (TEs) content from fifty-seven BAC clones (spanning 7.9 Mb) of this species, located in chromosomes by multiple fluorescence in situ hybridization (m-BAC-FISH) technique. The SSR analysis revealed an average density of 675.1 loci per Mb and a high abundance (59.69%) of dinucleotide coverage was observed, being ‘AC’ the most abundant. An SSR-FISH analysis using eleven probes was also carried out and seven of the 11 probes yielded positive signals. ‘AC’ probes were present as large clusters in almost all chromosomes, supporting the bioinformatic analysis. Regarding TEs, DNA transposons (Class II) were the most abundant. In Class I, LINE elements were the most abundant and the hAT family was the most represented in Class II. Rex/Babar subfamily, observed in two BAC clones mapping to chromosome pair 1, showed the longest match. This chromosome pair has been recently reported as a putative sexual proto-chromosome in this species, highlighting the possible role of the Rex element in the evolution of this chromosome. In the Rex1 phylogenetic tree, the Senegalese sole Rex1 retrotransposon could be associated with one of the four major ancient lineages in fish genomes, in which it is included O. latipes.

Chromosomal and molecular characterization of 5S rRNA genes in the North American abalones Haliotis rufescens Swainson (red abalone) and H. fulgens Philippi (blue abalone)

Abalone is an extremely valuable food source derived from cultured and wild animals, the later from populations under intense fishing exploitation and of high conservation value. As part of a long-term study to characterize genes from abalone that can be used as markers for hybrids certification, we characterised 5S ribosomal DNA (5S rDNA) in red abalone (Haliotis rufescens) and blue abalone (H. fulgens). The 5S rDNA arrays occur to a single pair of metacentric chromosomes at interstitial positions in both species. Two types of 5S genes were found, named types I and II, each associated with different non-transcribed spacer (NTS) sequences. The structure of the 5S rRNA genes and the NTS indicate incomplete homogenisation of the 5S rDNA arrays. The divergence of the 5S genes between species provide polymorphisms which can be used to distinguish red from blue abalone in forensic analysis of commercial production.

Chromosomal Mapping of Repetitive DNA Sequences in the Genus Bryconamericus (Characidae) and DNA Barcoding to Differentiate Populations

The mapping of repetitive DNA sites by fluorescence in situ hybridization has been widely used for karyotype studies in different species of fish, especially when dealing with related species or even genera presenting high chromosome variability. This study analyzed three populations of Bryconamericus, with diploid number preserved, but with different karyotype formulae. Bryconamericus ecai, from the Forquetinha river/RS, presented three new cytotypes, increasing the number of karyotype forms to seven in this population. Other two populations of Bryconamericus sp. from the Vermelho stream/PR and Cambuta river/PR exhibited interpopulation variation. The chromosome mapping of rDNA sites revealed unique markings among the three populations, showing inter- and intrapopulation variability located in the terminal region. The molecular analysis using DNA barcoding complementing the cytogenetic analysis also showed differentiation among the three populations. The U2 small nuclear DNA repetitive sequence exhibited conserved features, being located in the interstitial region of a single chromosome pair. This is the first report on its occurrence in the genus Bryconamericus. Data obtained revealed a karyotype variability already assigned to the genus, along with polymorphism of ribosomal sites, demonstrating that this group of fish can be undergoing a divergent evolutionary process, constituting a substantive model for studies of chromosomal evolution.

Evolutionary Dynamics of rDNAs and U2 Small Nuclear DNAs in Triportheus (Characiformes, Triportheidae): High Variability and Particular Syntenic Organization

Multigene families correspond to a group of genes tandemly repeated, showing enormous diversity in both number of units and genomic organization. In fishes, unlike rDNAs that have been well explored in cytogenetic studies, U2 small nuclear RNA (snRNA) genes are poorly investigated concerning their chromosomal localization. All Triportheus species (Characiformes, Triportheidae) studied so far carry a ZZ/ZW sex chromosomes system, where the W chromosome contains a huge 18S rDNA cistron. In some species the syntenic organization of rDNAs on autosomes was also verified. To explore this particular organization, we performed three-color-fluorescence in situ hybridization using 5S, 18S rDNA, and U2 snRNA genes as probes in eight Triportheus species. This work represents the first one analyzing the chromosomal distribution of U2 snRNA genes in genomes of Triportheidae. The variability in number of rDNA clusters, and the divergent syntenies for these three multigene families, put in evidence their evolutionary dynamism, revealing a much more complex organization of these genes than previously supposed for closely related species. Our study also provides additional data on the accumulation of repetitive sequences in the sex-specific chromosome. Besides, the chromosomal organization of U2 snDNAs among fish species is also reviewed.

Evidence of birth-and-death evolution of 5S rRNA gene in Channa species (Teleostei, Perciformes)

In higher eukaryotes, minor rDNA family codes for 5S rRNA that is arranged in tandem arrays and comprises of a highly conserved 120 bp long coding sequence with a variable non-transcribed spacer (NTS). Initially the 5S rDNA repeats are considered to be evolved by the process of concerted evolution. But some recent reports, including teleost fishes suggested that evolution of 5S rDNA repeat does not fit into the concerted evolution model and evolution of 5S rDNA family may be explained by a birth-and-death evolution model. In order to study the mode of evolution of 5S rDNA repeats in Perciformes fish species, nucleotide sequence and molecular organization of five species of genus Channa were analyzed in the present study. Molecular analyses revealed several variants of 5S rDNA repeats (four types of NTS) and networks created by a neighbor net algorithm for each type of sequences (I, II, III and IV) did not show a clear clustering in species specific manner. The stable secondary structure is predicted and upstream and downstream conserved regulatory elements were characterized. Sequence analyses also shown the presence of two putative pseudogenes in Channa marulius. Present study supported that 5S rDNA repeats in genus Channa were evolved under the process of birth-and-death.

Evolutionary Dynamics of 5S rDNA and Recurrent Association of Transposable Elements in Electric Fish of the Family Gymnotidae (Gymnotiformes): The Case of Gymnotus mamiraua

Gymnotidae is a family of electric fish endemic to the Neotropics consisting of 2 genera: Electrophorus and Gymnotus. The genus Gymnotus is widely distributed and is found in all of the major Brazilian river systems. Physical and molecular mapping data for the ribosomal DNA (rDNA) in this genus are still scarce, with its chromosomal location known in only 11 species. As other species of Gymnotus with 2n = 54 chromosomes from the Paraná-Paraguay basin, G. mamiraua was found to have a large number of 5S rDNA sites. Isolation and cloning of the 5S rDNA sequences from G. mamiraua identified a fragment of a transposable element similar to the Tc1/mariner transposon associated with a non-transcribed spacer. Double fluorescence in situ hybridization analysis of this element and the 5S rDNA showed that they were colocalized on several chromosomes, in addition to acting as nonsyntenic markers on others. Our data show the association between these sequences and suggest that the Tc1 retrotransposon may be the agent that drives the spread of these 5S rDNA-like sequences in the G. mamiraua genome.

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.

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.

Structure and Organization of the Engraulidae Family U2 snRNA: An Evolutionary Model Gene?

The U2 snRNA multigene family has been analyzed in four species of the Engraulidae family--Engraulis encrasicolus, Engraulis mordax, Engraulis ringens, and Engraulis japonicas--with the object of understanding more about the structure of this multigene family in these pelagic species and studying their phylogenetic relationships. The results showed that the cluster of this gene family in the Engraulis genus is formed by the U2-U5 snRNA with highly conserved sequences of mini- and micro-satellites, such as (CTGT)n, embedded downstream of the transcription unit; findings indicate that this gene family evolved following the concerted model. The phylogenetic analysis of the non-transcribed spacer of cluster U2-U5 snDNA in the 4 species showed that the sequences of the species E. encrasicolus and E. japonicus are closely related; these two are genetically close to E. mordax and slightly more distant from E. ringens. The data obtained by molecular analysis of U2-U5 snDNA and their secondary structure, with the presence of the micro-satellite (CTGT)n and mini-satellites, show clearly that the species E. encrasicolus and E. japonicus are closely related and would be older than E. mordax and E. ringens.

Expression profiling of the sex-related gene Dmrt1 in adults of the Lusitanian toadfish Halobatrachus didactylus (Bloch and Schneider, 1801)

Doublesex and mab-3 related transcription factor 1 (Dmrt1) gene is a widely conserved gene involved in sex determination and differentiation across phyla. To gain insights on Dmrt1 implication for fish gonad cell differentiation and gametogenesis development, its mRNA was isolated from testis and ovary from the Lusitanian toadfish (Halobatrachus didactylus). The cDNA from Dmrt1 was synthesized and cloned, whereas its quantitative and qualitative gene expression, as well as its protein immunolocalization, were analyzed. A main product of 1.38 kb, which encodes a protein of 295 aa, was reported, but other minority Dmrt1 products were also identified by RACE-PCR. This gene is predominantly expressed in testis (about 20 times more than in other organs or tissues), specially in spermatogonia, spermatocytes and spermatids, as well as in somatic Sertoli cells, indicating that Dmrt1 plays an important role in spermatogenesis. Although Dmrt1 transcripts also seem to be involved in oogenesis development, and it cannot be excluded that toadfish Dmrt1 could be functionally involved in other processes not related to sex.

Karyotype and cytogenetic mapping of 9 classes of repetitive DNAs in the genome of the naked catfish Mystus bocourti (Siluriformes, Bagridae)

BACKGROUND

In the present study, conventional and molecular cytogenetic studies were performed in the naked catfish Mystus bocourti (Siluriformes, Bagridae). Besides the conventional Giemsa staining, fluorescence in situ hybridization (FISH) using nine classes of repetitive DNAs namely 5S and 18S rDNAs, U2 snRNA, the microsatellites (CA)15 and (GA)15, telomeric repeats, and the retrotransposable elements Rex1, 3 and 6. was also performed.

RESULTS

M. bocourti had 2n = 56 chromosomes with a karyotype composed by 11 m + 11 sm + 6 st/a and a fundamental number (NF) equal to 100 in both sexes. Heteromorphic sex chromosome cannot be identified. The U2 snRNA, 5S and 18S rDNA were present in only one pair of chromosomes but none of them in a syntenic position. Microsatellites (CA)15 and (GA)15 showed hybridization signals at subtelomeric regions of all chromosomes with a stronger accumulation into one specific chromosomal pair. FISH with the telomeric probe revealed hybridization signals on each telomere of all chromosomes and interstitial telomeric sites (ITS) were not detected. The retrotransposable elements Rex1, 3 and 6 were generally spread throughout the genome.

CONCLUSIONS

In general, the repetitive sequences were not randomly distributed in the genome, suggesting a pattern of compartmentalization on the heterochromatic region of the chromosomes. Little is known about the structure and organization of bagrid genomes and the knowledge of the chromosomal distribution of repetitive DNA sequences in M. bocourti represents the first step for achieving an integrated view of their genomes.

High evolutionary dynamism in 5S rDNA of fish: state of the art

The 5S ribosomal DNA (rDNA) consists of one transcriptional unit of about 120 base pairs, which is separated from the next unit by a non-transcribed spacer (NTS). The coding sequence and the NTS together form a repeat unit which can be found in hundreds to thousands of copies tandemly repeated in the genomes. The NTS regions seem to be subject to rapid evolution. The first general model of evolution of these multigene families was referred to as divergent evolution, based on studies using hemoglobin and myoglobin as model systems. Later studies showed that nucleotide sequences of different multigene family members are more closely related within species than between species. This observation led to a new model of multigene family evolution, termed concerted evolution. Another model of evolution, named the birth-and-death model, has been found to be more suitable to explain the long-term evolution of these multigene families. According to this model, new genes originate by successive duplications, and these new genes are either maintained for a long time or are lost, or else degenerate into pseudogenes. In this review we describe different sources of variability in the 5S rDNA genes observed in several distinct fish species. This variability is mainly referred to NTSs and includes the presence of other multigene families (mainly LINEs, SINEs, non-LTR retrotransposons, and U snRNA families). Different types of microsatellites have also been found to contribute to the increase of variability in this region. Our recent results suggest that horizontal transfer contributes to the increase of diversity in the NTSs of some species. Variability in the 5S rDNA coding region affecting the stability of the structure, but without effects on the function of the 5S rRNA, is also described. Retrotransposons seem to be responsible for the high dynamism of 5S rDNA, while microsatellites acting as recombination hot spots could stabilize a wide variety of unusual DNA structures, affecting DNA replication and enhancing or decreasing promoter activity in gene expression. The relationship between the high variability found at molecular level and the low variability found at chromosomal level is also discussed.

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.