Molecular evolutionary characterization of the mussel Mytilus histone multigene family: first record of a tandemly repeated unit of five histone genes containing an H1 subtype with "orphon" features

Weak purifying selection in allelic diversity of the ADSL gene in indigenous and local chicken breeds and red junglefowl in Thailand

High levels of purine and uric acid, which are associated with health issues such as gout and cardiovascular disease, are found in the meat of fast-growing broiler chickens, which raises concerns about the quality of chicken meat and the health of the consumers who consume it. High genetic homogeneity and uniformity, particularly in genes involved in the synthesis of inosine monophosphate (IMP) and subsequent process of purine synthesis, which are associated with the meat quality, are exhibited in commercial broiler chickens owing to intensive inbreeding programs. Adenosine succinate lyase (ADSL) is a key enzyme involved in de novo purine biosynthetic pathway and its genetic polymorphisms affect IMP metabolism and purine content. In this study, we investigated the polymorphism of the ADSL gene in indigenous and local chicken breeds and red junglefowl in Thailand, using metabarcoding and genetic diversity analyses. Five alleles with 73 single nucleotide polymorphisms in exon 2, including missense and silent mutations, which may act on the synthesis efficiency of IMP and purine. Their protein structures revealed changes in amino acid composition that may affect ADSL enzyme activity. Weak purifying selection in these ADSL alleles was observed in the chicken population studied, implying that the variants have minor fitness impacts and a greater probability of fixation of beneficial mutations than strong purifying selection. A potential selective sweep was observed in Mae Hong Son chickens, whose purine content was lower than that in other breeds. This suggests a potential correlation between variations of the ADSL gene and reduced purine content and an impact of ADSL expression on the quality of chicken meat. However, further studies are required to validate its potential availability as a genetic marker for selecting useful traits that are beneficial to human health and well-being.

Histone derived antimicrobial peptides identified from Mytilus coruscus serum by peptidomics

Histones and their N-terminal or C-terminal derived peptides have been studied in vertebrates and presented as potential antimicrobial agents playing important roles in the innate immune defenses. Although histones and their derived peptides had been reported as components of innate immunity in invertebrates, the knowledge about the histone derived antimicrobial peptides (HDAPs) in invertebrates are still limited. Using a peptidomic technique, a set of peptide fragments derived from the histones was identified in this study from the serum of microbes challenged Mytilus coruscus. Among the 85 identified histone-derived-peptides with high confidence, 5 HDAPs were chemically synthesized and the antimicrobial activities were verified, showing strong growth inhibition against Gram-positive bacteria, Gram-negative bacteria, and fungus. The gene expression level of the precursor histones matched by representative HDAPs were further tested using q-PCR, and the results showed a significant upregulation of the histone gene expression levels in hemocytes, gill, and mantle of the mussel after immune stress. In addition, three identified HDAPs were selected for preparation of specific antibodies, and the corresponding histones and their derived C-terminal fragments were detected by Western blotting in the blood cell and serum of immune challenged mussel, respectively, indicating the existence of HDAPs in M. coruscus. Our findings revealed the immune function of histones in Mytilus, and confirmed the existence of HDAPs in the mussel. The identified Mytilus HDAPs represent a new source of immune effector with antimicrobial function in the innate immune system, and thus provide promising candidates for the treatment of microbial infections in aquaculture and medicine.

Genome-Wide Characterization and Analysis of Expression of the Histone Gene Family in Razor Clam, Sinonovacula constricta

The Chinese razor clam (Sinonovacula constricta), a bivalve species widely distributed in estuaries and mudflats, is often exposed to extreme environmental and microbial stresses. Histones are fundamental components of chromatin and play an important role in innate immunity, as demonstrated by its antimicrobial activities in clams. However, little attention has been paid to histones in bivalves. To fill this gap, we investigated the genomic distribution, structural characteristics, conserved motifs, and phylogenetic relationships of histones in S. constricta. A total of 114 histone genes were detected in the S. constricta genome, which were divided into 25 types in phylogenetic analysis. Among them, partial histones exhibited a tissue-dependent expression pattern, indicating that they may be involved in sustaining the homeostasis of organs/tissues in adult S. constricta. Furthermore, mRNA expression of certain histones changed significantly in S. constricta when infected with Vibrio parahaemolyticus, suggesting that histones play a role in the immune defense of S. constricta. All together, this study on histone genes in S. constricta not only greatly expands our knowledge of histone function in the clam, but also histone evolution in molluscs.

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.

Cytogenetic markers using single-sequence probes reveal chromosomal locations of tandemly repetitive genes in scleractinian coral Acropora pruinosa

The short and similar sized chromosomes of Acropora pose a challenge for karyotyping. Conventional methods, such as staining of heterochromatic regions, provide unclear banding patterns that hamper identification of such chromosomes. In this study, we used short single-sequence probes from tandemly repetitive 5S ribosomal RNA (rRNA) and core histone coding sequences to identify specific chromosomes of Acropora pruinosa. Both the probes produced intense signals in fluorescence in situ hybridization, which distinguished chromosome pairs. The locus of the 5S rDNA probe was on chromosome 5, whereas that of core histone probe was on chromosome 8. The sequence of the 5S rDNA probe was composed largely of U1 and U2 spliceosomal small nuclear RNA (snRNA) genes and their interspacers, flanked by short sequences of the 5S rDNA. This is the first report of a tandemly repetitive linkage of snRNA and 5S rDNA sequences in Cnidaria. Based on the constructed tentative karyogram and whole genome hybridization, the longest chromosome pair (chromosome 1) was heteromorphic. The probes also hybridized effectively with chromosomes of other Acropora species and population, revealing an additional core histone gene locus. We demonstrated the applicability of short-sequence probes as chromosomal markers with potential for use across populations and species of Acropora.

Environmental (e)RNA advances the reliability of eDNA by predicting its age

Environmental DNA (eDNA) analysis has advanced conservation biology and biodiversity management. However, accurate estimation of age and origin of eDNA is complicated by particle transport and the presence of legacy genetic material, which can obscure accurate interpretation of eDNA detection and quantification. To understand the state of genomic material within the environment, we investigated the degradation relationships between (a) size of fragments (long vs short), (b) genomic origins (mitochondrial vs nuclear), (c) nucleic acids (eDNA vs eRNA), and (d) RNA types (messenger (m)RNA vs ribosomal (r)RNA) from non-indigenous Dreissena mussels. Initial concentrations of eRNA followed expected transcriptional trends, with rRNAs found at > 1000 × that of eDNA, and a mitosis-associated mRNA falling below detection limits within 24 h. Furthermore, the ratio of eRNA:eDNA significantly decreased throughout degradation, potentially providing an estimate for the age of genomic material. Thus, eRNA quantification can increase detection due to the high concentrations of rRNAs. Furthermore, it may improve interpretation of positive detections through the eRNA:eDNA ratio and/or by detecting low abundant mitosis-associated mRNAs that degrade within ~ 24 h.

A single clonal lineage of transmissible cancer identified in two marine mussel species in South America and Europe

Transmissible cancers, in which cancer cells themselves act as an infectious agent, have been identified in Tasmanian devils, dogs, and four bivalves. We investigated a disseminated neoplasia affecting geographically distant populations of two species of mussels (Mytilus chilensis in South America and M. edulis in Europe). Sequencing alleles from four loci (two nuclear and two mitochondrial) provided evidence of transmissible cancer in both species. Phylogenetic analysis of cancer-associated alleles and analysis of diagnostic SNPs showed that cancers in both species likely arose in a third species of mussel (M. trossulus), but these cancer cells are independent from the previously identified transmissible cancer in M. trossulus from Canada. Unexpectedly, cancers from M. chilensis and M. edulis are nearly identical, showing that the same cancer lineage affects both. Thus, a single transmissible cancer lineage has crossed into two new host species and has been transferred across the Atlantic and Pacific Oceans and between the Northern and Southern hemispheres.

Sequence characterization of the 5S ribosomal DNA and the internal transcribed spacer (ITS) region in four European Donax species (Bivalvia: Donacidae)

Background

The whole repeat unit of 5S rDNA and the internal transcribed spacer (ITS) of four European Donax species were analysed. After amplifying, cloning and sequencing several 5S and ITS units, their basic features and their variation were described. The phylogenetic usefulness of 5S and ITS sequences in the inference of evolutionary relationships among these wedge clams was also investigated.

Results

The length of the 5S repeat presented little variation among species, except D. trunculus that differed from the rest of the Donax species in 170–210 bp. The deduced coding region covered 120 bp, and showed recognizable internal control regions (ICRs) involved in the transcription. The length of non-transcribed spacer region (NTS) ranged from 157 bp to 165 bp in Donax trunculus and from 335 bp to 367 bp in the other three species. The conservation degree of transcriptional regulatory regions was analysed revealing a conserved TATA-like box in the upstream region. Regarding ITS sequences, the four Donax species showed slight size differences among clones due to the variation occurring in the ITS1 and ITS2, except Donax variegatus did not display size differences in the ITS2. The total length of the ITS sequence ranged between 814 and 1014 bp. Resulting phylogenetic trees display that the two ribosomal DNA regions provide well-resolved phylogenies where the four European Donax species form a single clade receiving high support in nodes. The topology obtained with 5S sequences was in agreement with Donax evolutionary relationships inferred from several sequences of different nature in previous studies.

Conclusions

This is not only a basic research work, where new data and new knowledge is provided about Donax species, but also have allowed the authentication of these wedge clams and offers future applications to provide other genetic resources.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0684-x) contains supplementary material, which is available to authorized users.

Oxidative damage of 18S and 5S ribosomal RNA in digestive gland of mussels exposed to trace metals

Numerous studies have shown the ability of trace metals to accumulate in marine organisms and cause oxidative stress that leads to perturbations in many important intracellular processes, including protein synthesis. This study is mainly focused on the exploration of structural changes, like base modifications, scissions, and conformational changes, caused in 18S and 5S ribosomal RNA (rRNA) isolated from the mussel Mytilus galloprovincialis exposed to 40μg/L Cu, 30μg/L Hg, or 100μg/L Cd, for 5 or 15days. 18S rRNA and 5S rRNA are components of the small and large ribosomal subunit, respectively, found in complex with ribosomal proteins, translation factors and other auxiliary components (metal ions, toxins etc). 18S rRNA plays crucial roles in all stages of protein synthesis, while 5S rRNA serves as a master signal transducer between several functional regions of 28S rRNA. Therefore, structural changes in these ribosomal constituents could affect the basic functions of ribosomes and hence the normal metabolism of cells. Especially, 18S rRNA along with ribosomal proteins forms the decoding centre that ensures the correct codon-anticodon pairing. As exemplified by ELISA, primer extension analysis and DMS footprinting analysis, each metal caused oxidative damage to rRNA, depending on the nature of metal ion and the duration of exposure. Interestingly, exposure of mussels to Cu or Hg caused structural alterations in 5S rRNA, localized in paired regions and within loops A, B, C, and E, leading to a continuous progressive loss of the 5S RNA structural integrity. In contrast, structural impairments of 5S rRNA in mussels exposed to Cd were accumulating for the initial 5days, and then progressively decreased to almost the normal level by day 15, probably due to the parallel elevation of metallothionein content that depletes the pools of free Cd. Regions of interest in 18S rRNA, such as the decoding centre, sites implicated in the binding of tRNAs (A- and P-sites) or translation factors, and areas related to translation fidelity, were found to undergo significant metal-induced conformational alterations, leading either to loosening of their structure or to more compact folding. These modifications were associated with parallel alterations in the translation process at multiple levels, a fact suggesting that structural perturbations in ribosomes, caused by metals, pose significant hurdles in translational efficiency and fidelity.

Molecular and Biochemical Methods Useful for the Epigenetic Characterization of Chromatin-Associated Proteins in Bivalve Molluscs

Bivalve molluscs constitute a ubiquitous taxonomic group playing key functions in virtually all ecosystems, and encompassing critical commercial relevance. Along with a sessile and filter-feeding lifestyle in most cases, these characteristics make bivalves model sentinel organisms routinely used for environmental monitoring studies in aquatic habitats. The study of epigenetic mechanisms linking environmental exposure and specific physiological responses (i.e., environmental epigenetics) stands out as a very innovative monitoring strategy, given the role of epigenetic modifications in acclimatization and adaptation. Furthermore, the heritable nature of many of those modifications constitutes a very promising avenue to explore the applicability of epigenetic conditioning and selection in management and restoration strategies. Chromatin provides a framework for the study of environmental epigenetic responses. Unfortunately, chromatin and epigenetic information are very limited in most non-traditional model organisms and even completely lacking in most environmentally and ecologically relevant organisms. The present work aims to provide a comprehensive and reproducible experimental workflow for the study of bivalve chromatin. First, a series of guidelines for the molecular isolation of genes encoding chromatin-associated proteins is provided, including information on primers suitable for conventional PCR, Rapid Amplification of cDNA Ends (RACE), genome walking and quantitative PCR (qPCR) experiments. This section is followed by the description of methods specifically developed for the analysis of histone and SNBP proteins in different bivalve tissues, including protein extraction, purification, separation and immunodetection. Lastly, information about available antibodies, their specificity and performance is also provided. The tools and protocols described here complement current epigenetic analyses (usually limited to DNA methylation) by incorporating the study of structural elements modulating chromatin dynamics.

The cnidarian Hydractinia echinata employs canonical and highly adapted histones to pack its DNA

Background

Cnidarians are a group of early branching animals including corals, jellyfish and hydroids that are renowned for their high regenerative ability, growth plasticity and longevity. Because cnidarian genomes are conventional in terms of protein-coding genes, their remarkable features are likely a consequence of epigenetic regulation. To facilitate epigenetics research in cnidarians, we analysed the histone complement of the cnidarian model organism Hydractinia echinata using phylogenomics, proteomics, transcriptomics and mRNA in situ hybridisations.

Results

We find that the Hydractinia genome encodes 19 histones and analyse their spatial expression patterns, genomic loci and replication-dependency. Alongside core and other replication-independent histone variants, we find several histone replication-dependent variants, including a rare replication-dependent H3.3, a female germ cell-specific H2A.X and an unusual set of five H2B variants, four of which are male germ cell-specific. We further confirm the absence of protamines in Hydractinia.

Conclusions

Since no protamines are found in hydroids, we suggest that the novel H2B variants are pivotal for sperm DNA packaging in this class of Cnidaria. This study adds to the limited number of full histone gene complements available in animals and sets a comprehensive framework for future studies on the role of histones and their post-translational modifications in cnidarian epigenetics. Finally, it provides insight into the evolution of spermatogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0085-1) contains supplementary material, which is available to authorized users.

Tralopyril bioconcentration and effects on the gill proteome of the Mediterranean mussel Mytilus galloprovincialis

Antifouling (AF) systems are used worldwide as one of the most cost-effective ways of protecting submerged structures against heavy biofouling. The emergence of environmentally friendly AF biocides requires knowledge on their environmental fate and toxicity. In this study we measured the bioconcentration of the emerging AF biocide tralopyril (TP) in the Mediterranean mussel Mytilus galloprovincialis and investigated the effects of TP on the mussel gill proteome following acute (2days) and chronic (30days) exposure, as well as after a 10-day depuration period. The experiments were carried out with 1μg/L TP; blank and solvent (5×10(-5)% DMSO) controls were also included. Proteomics analysis was performed by mass spectrometry-based multidimensional protein identification technology (MudPIT). Differentially expressed proteins were identified using a label-free approach based on spectral counts and G-test. Our results show that TP is rapidly accumulated by mussels at concentrations up to 362ng/g dw (whole tissues), reaching steady-state condition within 13days. Ten days of depuration resulted in 80% elimination of accumulated TP from the organism, suggesting that a complete elimination could be reached with longer depuration times. In total, 46 proteins were found to be regulated in the different exposure scenarios. Interestingly, not only TP but also DMSO alone significantly modulated the protein expression in mussel gills following acute and chronic exposure. Both compounds regulated proteins involved in bioenergetics, immune system, active efflux and oxidative stress, often in the opposite way. Alterations of several proteins, notably several cytoskeletal ones, were still observed after the depuration period. These may reflect either the continuing chemical effect due to incomplete elimination or an onset of recovery processes in the mussel gills. Our study shows that exposure of adult mussels to sublethal TP concentration results in the bioconcentration of this biocide in the tissues and modulates the expression of several proteins that may intervene in important metabolic pathways.

Environmental epigenetics: A promising venue for developing next-generation pollution biomonitoring tools in marine invertebrates

Environmental epigenetics investigates the cause-effect relationships between specific environmental factors and the subsequent epigenetic modifications triggering adaptive responses in the cell. Given the dynamic and potentially reversible nature of the different types of epigenetic marks, environmental epigenetics constitutes a promising venue for developing fast and sensible biomonitoring programs. Indeed, several epigenetic biomarkers have been successfully developed and applied in traditional model organisms (e.g., human and mouse). Nevertheless, the lack of epigenetic knowledge in other ecologically and environmentally relevant organisms has hampered the application of these tools in a broader range of ecosystems, most notably in the marine environment. Fortunately, that scenario is now changing thanks to the growing availability of complete reference genome sequences along with the development of high-throughput DNA sequencing and bioinformatic methods. Altogether, these resources make the epigenetic study of marine organisms (and more specifically marine invertebrates) a reality. By building on this knowledge, the present work provides a timely perspective highlighting the extraordinary potential of environmental epigenetic analyses as a promising source of rapid and sensible tools for pollution biomonitoring, using marine invertebrates as sentinel organisms. This strategy represents an innovative, groundbreaking approach, improving the conservation and management of natural resources in the oceans.

Divergent evolutionary behavior of H3 histone gene and rDNA clusters in venerid clams

BACKGROUND

Histone H3 gene clusters have been described as highly conserved chromosomal markers in invertebrates. Surprisingly, in bivalves remarkable interspecific differences were found among the eight mussels and between the two clams in which histone H3 gene clusters have already been located. Although the family Veneridae comprises 10 % of the species of marine bivalves, their chromosomes are poorly studied. The clams belonging to this family present 2n = 38 chromosomes and similar karyotypes showing chromosome pairs gradually decreasing in length. In order to assess the evolutionary behavior of histone and rRNA multigene families in bivalves, we mapped histone H3 and ribosomal RNA probes to chromosomes of ten species of venerid clams.

RESULTS

In contrast with the reported conservation of histone H3 gene clusters and their intercalary location in invertebrates, these loci varied in number and were mostly subterminal in venerid clams. On the other hand, while a single 45S rDNA cluster, highly variable in location, was found in these organisms, 5S rDNA clusters showed interspecific differences in both number and location. The distribution patterns of these sequences were species-specific and mapped to different chromosomal positions in all clams but Ruditapes decussatus, in which one of the minor rDNA clusters and the major rDNA cluster co-located.

CONCLUSION

The diversity in the distribution patterns of histone H3 gene, 5S rDNA and 28S rDNA clusters found in venerid clams, together with their different evolutionary behaviors in other invertebrate taxa, strongly suggest that the control of the spreading of these multigene families in a group of organisms relies upon a combination of evolutionary forces that operate differently depending not only on the specific multigene family but also on the particular taxa. Our data also showed that H3 histone gene and rDNA clusters are useful landmarks to integrate nex-generation sequencing (NGS) and evolutionary genomic data in non-model species.

High-Resolution Physical Chromosome Mapping of Multigene Families in Lagria villosa (Tenebrionidae): Occurrence of Interspersed Ribosomal Genes in Coleoptera

The organization and mapping of multigene families can produce useful genetic markers, and its use may elucidate the mechanisms of karyotype variation and genomic organization in different groups of eukaryotes. To date, few species of Coleoptera have been analyzed using FISH for the location of multigene families. The purpose of this study was to use high-resolution chromosome mapping to establish the genomic organization of the 18S rDNA, 5S rDNA and histone H3 gene families in Lagria villosa. FISH was performed using 18S rDNA, 5S rDNA and histone H3 probes prepared via PCR labeling. Fiber-FISH for 18S and 5S rDNA indicated that both ribosomal elements are colocalized in the short arm of chromosome 4. Additionally, FISH, using the histone H3 probe, revealed that this sequence is found in only one autosomal pair and did not colocalize with rDNA. Fiber-FISH with 5S and 18S probes, used to improve the mapping resolution of these regions, showed that both genes are closely interspersed with varying amounts of both DNA classes.

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.

Karyotypic diversification in Mytilus mussels (Bivalvia: Mytilidae) inferred from chromosomal mapping of rRNA and histone gene clusters

Background

Mussels of the genus Mytilus present morphologically similar karyotypes that are presumably conserved. The absence of chromosome painting probes in bivalves makes difficult verifying this hypothesis. In this context, we comparatively mapped ribosomal RNA and histone gene families on the chromosomes of Mytilus edulis, M. galloprovincialis, M. trossulus and M. californianus by fluorescent in situ hybridization (FISH).

Results

Major rRNA, core and linker histone gene clusters mapped to different chromosome pairs in the four taxa. In contrast, minor rRNA gene clusters showed a different behavior. In all Mytilus two of the 5S rDNA clusters mapped to the same chromosome pair and one of them showed overlapping signals with those corresponding to one of the histone H1 gene clusters. The overlapping signals on mitotic chromosomes became a pattern of alternate 5S rRNA and linker histone gene signals on extended chromatin fibers. Additionally, M. trossulus showed minor and major rDNA clusters on the same chromosome pair.

Conclusion

The results obtained suggest that at least some of the chromosomes bearing these sequences are orthologous and that chromosomal mapping of rRNA and histone gene clusters could be a good tool to help deciphering some of the many unsolved questions in the systematic classification of Mytilidae.

Bivalve omics: state of the art and potential applications for the biomonitoring of harmful marine compounds

The extraordinary progress experienced by sequencing technologies and bioinformatics has made the development of omic studies virtually ubiquitous in all fields of life sciences nowadays. However, scientific attention has been quite unevenly distributed throughout the different branches of the tree of life, leaving molluscs, one of the most diverse animal groups, relatively unexplored and without representation within the narrow collection of well established model organisms. Within this Phylum, bivalve molluscs play a fundamental role in the functioning of the marine ecosystem, constitute very valuable commercial resources in aquaculture, and have been widely used as sentinel organisms in the biomonitoring of marine pollution. Yet, it has only been very recently that this complex group of organisms became a preferential subject for omic studies, posing new challenges for their integrative characterization. The present contribution aims to give a detailed insight into the state of the art of the omic studies and functional information analysis of bivalve molluscs, providing a timely perspective on the available data resources and on the current and prospective applications for the biomonitoring of harmful marine compounds.

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.

Systematic analysis and evolution of 5S ribosomal DNA in metazoans

Several studies on 5S ribosomal DNA (5S rDNA) have been focused on a subset of the following features in mostly one organism: number of copies, pseudogenes, secondary structure, promoter and terminator characteristics, genomic arrangements, types of non-transcribed spacers and evolution. In this work, we systematically analyzed 5S rDNA sequence diversity in available metazoan genomes, and showed organism-specific and evolutionary-conserved features. Putatively functional sequences (12,766) from 97 organisms allowed us to identify general features of this multigene family in animals. Interestingly, we show that each mammal species has a highly conserved (housekeeping) 5S rRNA type and many variable ones. The genomic organization of 5S rDNA is still under debate. Here, we report the occurrence of several paralog 5S rRNA sequences in 58 of the examined species, and a flexible genome organization of 5S rDNA in animals. We found heterogeneous 5S rDNA clusters in several species, supporting the hypothesis of an exchange of 5S rDNA from one locus to another. A rather high degree of variation of upstream, internal and downstream putative regulatory regions appears to characterize metazoan 5S rDNA. We systematically studied the internal promoters and described three different types of termination signals, as well as variable distances between the coding region and the typical termination signal. Finally, we present a statistical method for detection of linkage among noncoding RNA (ncRNA) gene families. This method showed no evolutionary-conserved linkage among 5S rDNAs and any other ncRNA genes within Metazoa, even though we found 5S rDNA to be linked to various ncRNAs in several clades.

The CHROMEVALOA database: a resource for the evaluation of Okadaic Acid contamination in the marine environment based on the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis

Okadaic Acid (OA) constitutes the main active principle in Diarrhetic Shellfish Poisoning (DSP) toxins produced during Harmful Algal Blooms (HABs), representing a serious threat for human consumers of edible shellfish. Furthermore, OA conveys critical deleterious effects for marine organisms due to its genotoxic potential. Many efforts have been dedicated to OA biomonitoring during the last three decades. However, it is only now with the current availability of detailed molecular information on DNA organization and the mechanisms involved in the maintenance of genome integrity, that a new arena starts opening up for the study of OA contamination. In the present work we address the links between OA genotoxicity and chromatin by combining Next Generation Sequencing (NGS) technologies and bioinformatics. To this end, we introduce CHROMEVALOAdb, a public database containing the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis (a sentinel model organism) in response to OA exposure. This resource constitutes a leap forward for the development of chromatin-based biomarkers, paving the road towards the generation of powerful and sensitive tests for the detection and evaluation of the genotoxic effects of OA in coastal areas.

Chromatin specialization in bivalve molluscs: a leap forward for the evaluation of Okadaic Acid genotoxicity in the marine environment

Marine biotoxins synthesized by Harmful Algal Blooms (HABs) represent one of the most important sources of contamination in marine environments as well as a serious threat to fisheries and aquaculture-based industries in coastal areas. Among these biotoxins Okadaic Acid (OA) is of critical interest as it represents the most predominant Diarrhetic Shellfish Poisoning biotoxin in the European coasts. Furthermore, OA is a potent tumor promoter with aneugenic and clastogenic effects on the hereditary material, most notably DNA breaks and alterations in DNA repair mechanisms. Therefore, a great effort has been devoted to the biomonitoring of OA in the marine environment during the last two decades, mainly based on physicochemical and physiological parameters using mussels as sentinel organisms. However, the molecular genotoxic effects of this biotoxin make chromatin structure a good candidate for an alternative strategy for toxicity assessment with faster and more sensitive evaluation. To date, the development of chromatin-based studies to this purpose has been hampered by the complete lack of information on chromatin of invertebrate marine organisms, especially in bivalve molluscs. Our preliminary results have revealed the presence of histone variants involved in DNA repair and chromatin specialization in mussels and clams. In this work we use this information to put forward a proposal focused on the development of chromatin-based tests for OA genotoxicity in the marine environment. The implementation of such tests in natural populations has the potential to provide an important leap in the biomonitoring of this biotoxin. The outcome of such monitoring may have critical implications for the evaluation of DNA damage in these marine organisms. They will provide as well important tools for the optimization of their harvesting and for the elaboration of additional tests designed to evaluate the safety of their consumption and potential implications for consumer's health.

Histone H2A (H2A.X and H2A.Z) variants in molluscs: molecular characterization and potential implications for chromatin dynamics

Histone variants are used by the cell to build specialized nucleosomes, replacing canonical histones and generating functionally specialized chromatin domains. Among many other processes, the specialization imparted by histone H2A (H2A.X and H2A.Z) variants to the nucleosome core particle constitutes the earliest response to DNA damage in the cell. Consequently, chromatin-based genotoxicity tests have been developed in those cases where enough information pertaining chromatin structure and dynamics is available (i.e., human and mouse). However, detailed chromatin knowledge is almost absent in most organisms, specially protostome animals. Molluscs (which represent sentinel organisms for the study of pollution) are not an exception to this lack of knowledge. In the present work we first identified the existence of functionally differentiated histone H2A.X and H2A.Z variants in the mussel Mytilus galloprovincialis (MgH2A.X and MgH2A.Z), a marine organism widely used in biomonitoring programs. Our results support the functional specialization of these variants based on: a) their active expression in different tissues, as revealed by the isolation of native MgH2A.X and MgH2A.Z proteins in gonad and hepatopancreas; b) the evolutionary conservation of different residues encompassing functional relevance; and c) their ability to confer specialization to nucleosomes, as revealed by nucleosome reconstitution experiments using recombinant MgH2A.X and MgH2A.Z histones. Given the seminal role of these variants in maintaining genomic integrity and regulating gene expression, their preliminary characterization opens up new potential applications for the future development of chromatin-based genotoxicity tests in pollution biomonitoring programs.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Chromosomal organization of the 18S and 5S rRNAs and histone H3 genes in Scarabaeinae coleopterans: insights into the evolutionary dynamics of multigene families and heterochromatin

BACKGROUND

Scarabaeinae beetles show a high level of macro-chromosomal variability, although the karyotypic organization of heterochromatin and multigene families (rDNAs and histone genes) is poorly understood in this group. To better understand the chromosomal organization and evolution in this group, we analyzed the karyotypes, heterochromatin distribution and chromosomal locations of the rRNAs and histone H3 genes in beetles belonging to eight tribes from the Scarabaeinae subfamily (Coleoptera, Scarabaeidae).

RESULTS

The number of 18S rRNA gene (a member of the 45S rDNA unit) sites varied from one to 16 and were located on the autosomes, sex chromosomes or both, although two clusters were most common. Comparison of the 45S rDNA cluster number and the diploid numbers revealed a low correlation value. However, a comparison between the number of 45S rDNA sites per genome and the quantity of heterochromatin revealed (i) species presenting heterochromatin restricted to the centromeric/pericentromeric region that contained few rDNA sites and (ii) species with a high quantity of heterochromatin and a higher number of rDNA sites. In contrast to the high variability for heterochromatin and 45S rDNA cluster, the presence of two clusters (one bivalent cluster) co-located on autosomal chromosomes with the 5S rRNA and histone H3 genes was highly conserved.

CONCLUSIONS

Our results indicate that the variability of the 45S rDNA chromosomal clusters is not associated with macro-chromosomal rearrangements but are instead related to the spread of heterochromatin. The data obtained also indicate that both heterochromatin and the 45S rDNA loci could be constrained by similar evolutionary forces regulating spreading in the distinct Scarabaeinae subfamily lineages. For the 5S rRNA and the histone H3 genes, a similar chromosomal organization could be attributed to their association/co-localization in the Scarabaeinae karyotypes. These data provide evidence that different evolutionary forces act at the heterochromatin and the 45S rDNA loci compared to the 5S rRNA and histone H3 genes during the evolution of the Scarabainae karyotypes.

Chromosomal organization of the 18S and 5S rRNAs and histone H3 genes in Scarabaeinae coleopterans: insights into the evolutionary dynamics of multigene families and heterochromatin

Background

Scarabaeinae beetles show a high level of macro-chromosomal variability, although the karyotypic organization of heterochromatin and multigene families (rDNAs and histone genes) is poorly understood in this group. To better understand the chromosomal organization and evolution in this group, we analyzed the karyotypes, heterochromatin distribution and chromosomal locations of the rRNAs and histone H3 genes in beetles belonging to eight tribes from the Scarabaeinae subfamily (Coleoptera, Scarabaeidae).

Results

The number of 18S rRNA gene (a member of the 45S rDNA unit) sites varied from one to 16 and were located on the autosomes, sex chromosomes or both, although two clusters were most common. Comparison of the 45S rDNA cluster number and the diploid numbers revealed a low correlation value. However, a comparison between the number of 45S rDNA sites per genome and the quantity of heterochromatin revealed (i) species presenting heterochromatin restricted to the centromeric/pericentromeric region that contained few rDNA sites and (ii) species with a high quantity of heterochromatin and a higher number of rDNA sites. In contrast to the high variability for heterochromatin and 45S rDNA cluster, the presence of two clusters (one bivalent cluster) co-located on autosomal chromosomes with the 5S rRNA and histone H3 genes was highly conserved.

Conclusions

Our results indicate that the variability of the 45S rDNA chromosomal clusters is not associated with macro-chromosomal rearrangements but are instead related to the spread of heterochromatin. The data obtained also indicate that both heterochromatin and the 45S rDNA loci could be constrained by similar evolutionary forces regulating spreading in the distinct Scarabaeinae subfamily lineages. For the 5S rRNA and the histone H3 genes, a similar chromosomal organization could be attributed to their association/co-localization in the Scarabaeinae karyotypes. These data provide evidence that different evolutionary forces act at the heterochromatin and the 45S rDNA loci compared to the 5S rRNA and histone H3 genes during the evolution of the Scarabainae karyotypes.

Cytogenetic characterization of the invasive mussel species Xenostrobus securis Lmk. (Bivalvia: Mytilidae)

The chromosomes of the invasive black-pigmy mussel (Xenostrobus securis (Lmk. 1819)) were analyzed by means of 4',6-diamidino-2-phenylindole (DAPI) / propidium iodide (PI) and chromomycin A3 (CMA) / DAPI fluorescence staining and fluorescent in situ hybridization using major rDNA, 5S rDNA, core histone genes, linker histone genes, and telomeric sequences as probes. The diploid chromosome number in this species is 2n = 30. The karyotype is composed of seven metacentric, one meta/submetacentric, and seven submetacentric chromosome pairs. Telomeric sequences appear at both ends of every single chromosome. Major rDNA clusters appear near the centromeres on chromosome pairs 1 and 3 and are associated with bright CMA fluorescence and dull DAPI fluorescence. This species shows five 5S rDNA clusters close to the centromeres on four chromosome pairs (2, 5, 6, and 8). Three of the four core histone gene clusters map to centromeric positions on chromosome pairs 7, 10, and 13. The fourth core histone gene cluster occupies a terminal position on chromosome pair 8, also bearing a 5S rDNA cluster. The two linker histone gene clusters are close to the centromeres on chromosome pairs 12 and 14. Therefore, the use of these probes allows the unequivocal identification of 11 of the 15 chromosome pairs that compose the karyotype of X. securis.

Cytogenetic characterization and mapping of rDNAs, core histone genes and telomeric sequences in Venerupis aurea and Tapes rhomboides (Bivalvia: Veneridae)

We describe the chromosomal location of GC-rich regions, 28S and 5S rDNA, core histone genes, and telomeric sequences in the veneroid bivalve species Venerupis aurea and Tapes (Venerupis) rhomboides, using fluorochrome staining with propidium iodide, DAPI and chromomycin A3 (CMA) and fluorescent in situ hybridization (FISH). DAPI dull/CMA bright bands were coincident with the chromosomal location of 28S rDNA in both species. The major rDNA was interstitially clustered at a single locus on the short arms of the metacentric chromosome pair 5 in V. aurea, whereas in T. rhomboides it was subtelomerically clustered on the long arms of the subtelocentric chromosome pair 17. 5S rDNA also was a single subtelomeric cluster on the long arms of subtelocentric pair 17 in V. aurea and on the short arms of the metacentric pair 9 in T. rhomboides. Furthermore, V. aurea showed four telomeric histone gene clusters on three metacentric pairs, at both ends of chromosome 2 and on the long arms of chromosomes 3 and 8, whereas histone genes in T. rhomboides clustered interstitially on the long arms of the metacentric pair 5 and proximally on the long arms of the subtelocentric pair 12. Double and triple FISH experiments demonstrated that rDNA and H3 histone genes localized on different chromosome pairs in the two clam species. Telomeric signals were found at both ends of every single chromosome in both species. Chromosomal location of these three gene families in two species of Veneridae provides a clue to karyotype evolution in this commercially important bivalve family.

Chromosomal mapping of rRNA genes, core histone genes and telomeric sequences in Brachidontes puniceus and Brachidontes rodriguezi (Bivalvia, Mytilidae)

Background

Chromosome rearrangements are an important part of the speciation process in many taxa. The study of chromosome evolution in bivalves is hampered by the absence of clear chromosomal banding patterns and the similarity in both chromosome size and morphology. For this reason, obtaining good chromosome markers is essential for reliable karyotypic comparisons. To begin this task, the chromosomes of the mussels Brachidontes puniceus and B. rodriguezi were studied by means of fluorochrome staining and fluorescent in situ hybridization (FISH).

Results

Brachidontes puniceus and B. rodriguezi both have 2n = 32 chromosomes but differing karyotype composition. Vertebrate-type telomeric sequences appear at both ends of every single chromosome. B. puniceus presents a single terminal major rRNA gene cluster on a chromosome pair while B. rodriguezi shows two. Both mussels present two 5S rDNA and two core histone gene clusters intercalary located on the long arms of two chromosome pairs. Double and triple-FISH experiments demonstrated that one of the 5S rDNA and one of the major rDNA clusters appear on the same chromosome pair in B. rodriguezi but not in B. puniceus. On the other hand, the second 5S rDNA cluster is located in one of the chromosome pairs also bearing one of the core histone gene clusters in the two mussel species.

Conclusion

Knowledge of the chromosomal distribution of these sequences in the two species of Brachidontes is a first step in the understanding of the role of chromosome changes on bivalve evolution.

Histone and ribosomal RNA repetitive gene clusters of the boll weevil are linked in a tandem array

Histones are the major protein component of chromatin structure. The histone family is made up of a quintet of proteins, four core histones (H2A, H2B, H3 & H4) and the linker histones (H1). Spacers are found between the coding regions. Among insects this quintet of genes is usually clustered and the clusters are tandemly repeated. Ribosomal DNA contains a cluster of the rRNA sequences 18S, 5.8S and 28S. The rRNA genes are separated by the spacers ITS1, ITS2 and IGS. This cluster is also tandemly repeated. We found that the ribosomal RNA repeat unit of at least two species of Anthonomine weevils, Anthonomus grandis and Anthonomus texanus (Coleoptera: Curculionidae), is interspersed with a block containing the histone gene quintet. The histone genes are situated between the rRNA 18S and 28S genes in what is known as the intergenic spacer region (IGS). The complete reiterated Anthonomus grandis histone-ribosomal sequence is 16,248 bp.

Evolutionary dynamics of the 5S rDNA gene family in the mussel Mytilus: mixed effects of birth-and-death and concerted evolution

In higher eukaryotes, the gene family encoding the 5S ribosomal RNA (5S rRNA) has been used (together with histones) to showcase the archetypal example of a gene family subject to concerted evolution. However, recent studies have revealed conspicuous features challenging the predictions of this model, including heterogeneity of repeat units, the presence of functional 5S gene variants as well as the existence of 5S rDNA divergent pseudogenes lacking traces of homogenization. In the present work, we have broadened the scope in the evolutionary study of ribosomal gene families by studying the 5S rRNA family in mussels, a model organism which stands out among other animals due to the heterogeneity it displays regarding sequence and organization. To this end, 48 previously unknown 5S rDNA units (coding and spacer regions) were sequenced in five mussel species, leading to the characterization of two new types of units (referred to here as small-beta 5S rDNA and gamma-5S rDNA) coexisting in the genome with alpha and beta rDNA units. The intense genetic dynamics of this family is further supported by the first description of an association between gamma-5S rDNA units and tRNA genes. Molecular evolutionary and phylogenetic analyses revealed an extensive lack of homology among spacer sequences belonging to different rDNA types, suggesting the presence of independent evolutionary pathways leading to their differentiation. Overall, our results suggest that the long-term evolution of the 5S rRNA gene family in mussels is most likely mediated by a mixed mechanism involving the generation of genetic diversity through birth-and-death, followed by a process of local homogenization resulting from concerted evolution in order to maintain the genetic identities of the different 5S units, probably after their transposition to independent chromosomal locations.

Birth-and-death long-term evolution promotes histone H2B variant diversification in the male germinal cell line

The rich diversity within each of the five histone families (H1, H2A, H2B, H3, and H4) can hardly be reconciled with the notion of homogenizing evolution. The prevalence of birth-and-death long-term evolution over concerted evolution has already been demonstrated in the linker histone H1 family as well as for the H2A, H3, and H4 core histone families. However, information about histone H2B is lacking. In the present work, we have analyzed the diversity of the members of this histone family across different eukaryotic genomes and have characterized the mechanisms involved in their long-term evolution. Our results reveal that, quite in contrast with other histones, H2B variants are subject to a very rapid process of diversification that primarily affects the male germinal cell lineage and involves their functional specialization probably as a consequence of neofunctionalization and subfunctionalization events after gene duplication. The overall parallelism observed between the molecular phylogenies and the relationships among the electrostatic potentials of the different variants suggests that the latter may have played a major structural selective constraint during H2B evolution. It thus seems that the reorganization of chromatin structure during spermiogenesis might have affected the evolutionary constraints driving histone H2B evolution, leading to an increase in diversity.

Histone genes of the razor clam Solen marginatus unveil new aspects of linker histone evolution in protostomes

The association of DNA with histones results in a nucleoprotein complex called chromatin that consists of repetitive nucleosomal subunits. Nucleosomes are joined together in the chromatin fiber by short stretches of linker DNA that interact with a wide diversity of linker H1 histones involved in chromatin compaction and dynamics. Although the long-term evolution of the H1 family has been the subject of different studies during the last 5 years, the lack of molecular data on replication-independent (RI) H1 variants from protostomes has been hampering attempts to complete the evolutionary picture of this histone family in eukaryotes, especially as it pertains to the functional specialization they impart to the chromatin structure in members of this bilaterian lineage. In an attempt to fill this gap, the present work characterizes the histone gene complement from the razor clam Solen marginatus. Molecular evolutionary analyses reveal that the H1 gene from this organism represents one of the few protostome RI H1 genes known to date, a notion which is further supported by its location within the monophyletic group encompassing the RI H1 variants in the overall phylogeny of eukaryotic H1 proteins. Although the detailed characterization of the nucleotide substitution patterns in RI H1 variants agrees with the model of birth-and-death evolution under strong purifying selection, maximum-likelihood approaches unveil the presence of adaptive selection during at least part of the evolutionary differentiation between protostomes and deuterostomes. The presence of increased levels of specialization in RI H1 proteins from deuterostomes as well as the significant differences observed in electrostatic properties between protostome and deuterostome RI H1s represent novel and important preliminary results for future studies of the functional differentiation of this histone H1 lineage across bilaterians.

The evolutionary differentiation of two histone H2A.Z variants in chordates (H2A.Z-1 and H2A.Z-2) is mediated by a stepwise mutation process that affects three amino acid residues

Background

The histone H2A family encompasses the greatest number of core histone variants of which the replacement variant H2A.Z is currently one of the most heavily studied. No clear mechanism for the functional variability that H2A.Z imparts to chromatin has yet been proposed. While most of the past studies have referred to H2A.Z generically as a single protein, in vertebrates it is a mixture of two protein forms H2A.Z-1 (previously H2A.Z) and H2A.Z-2 (previously H2A.F/Z or H2A.V) that differ by three amino acids.

Results

We have performed an extensive study on the long-term evolution of H2A.Z across metazoans with special emphasis on the possible selective mechanisms responsible for the differentiation between H2A.Z-1 and H2A.Z-2. Our results reveal a common origin of both forms early in chordate evolution. The evolutionary process responsible for the differentiation involves refined stepwise mutation change within the codons of the three differential residues. This eventually led to differences in the intensity of the selective constraints acting upon the different H2A.Z forms in vertebrates.

Conclusion

The results presented in this work definitively reveal that the existence of H2A.Z-1 and H2A.Z-2 is not a whim of random genetic drift. Our analyses demonstrate that H2A.Z-2 is not only subject to a strong purifying selection but it is significantly more evolutionarily constrained than H2A.Z-1. Whether or not the evolutionary drift between H2A.Z-1 and H2A.Z-2 has resulted in a functional diversification of these proteins awaits further research. Nevertheless, the present work suggests that in the process of their differently constrained evolutionary pathways, these two forms may have acquired new or complementary functions.

FISH mapping and identification of Zhikong scallop (Chlamys farreri) chromosomes

Chromosome identification is the first step in genomic research of a species, but it remains a challenge in scallops. In the present study, fluorescence in situ hybridization (FISH) mapping of 19 fosmid clones was attempted and used for chromosome identification in Zhikong scallop (Chlamys farreri Jones et Preston, 1904). Data showed that 10 clones were successfully mapped, including 7 without and 3 with C ( 0 ) t-1 DNA. Among them, 2 represented multiple signals and made no contribution to chromosome identification. Karyotypic analysis and cohybridization indicated that the remaining 8 clones realized the identification of 8 chromosomes. All 10 clones were sequenced at both ends, which could be developed as sequence-tagged sites and used for the unification of the cytological and genetic linkage maps. This study shows that fosmid clones can benefit chromosome identification and will undoubtedly be useful for cytogenetic research in Zhikong scallop.

The karyotype and 5S rRNA genes from Spanish individuals of the bat species Rhinolophus hipposideros (Rhinolophidae; Chiroptera)

The karyotype of individuals of the species Rhinolophus hipposideros from Spain present a chromosome number of 2n = 54 (NFa = 62). The described karyotype for these specimens is very similar to another previously described in individual from Bulgaria. However, the presence of one additional pair of autosomal acrocentric chromosomes in the Bulgarian karyotype and the differences in X chromosome morphology indicated that we have described a new karyotype variant in this species. In addition, we have analyzed several clones of 1.4 and 1 kb of a PstI repeated DNA sequence from the genome of R. hipposideros. The repeated sequence included a region with high identity with the 5S rDNA genes and flanking regions, with no homology with GenBank sequences. Search for polymerase III regulatory elements demonstrated the presence of type I promoter elements (A-box, Intermediate Element and C-box) in the 5S rDNA region. In addition, upstream regulatory elements, as a D-box and Sp1 binding sequences, were present in flanking regions. All data indicated that the cloned repeated sequences are the functional rDNA genes from this species. Finally, FISH demonstrated the presence of rDNA in nine chromosome pairs, which is surprising as most mammals have only one carrier chromosome pair.