Next Article in Journal
The Snakeskin Gourami (Trichopodus pectoralis) Tends to Exhibit XX/XY Sex Determination
Previous Article in Journal
Influence of Fishery Subsidies on Fishing: Empirical Test Based on China’s Provincial Panel Data
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Fishes
Volume 6
Issue 3
10.3390/fishes6030042
fishes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of the nrDNA Repeat Unit of Manila Clam Ruditapes philippinarum and Quahog Mercenaria mercenaria

by
Zhansheng Guo
1,
Zhen Wang
2 and
Xuguang Hou
1,*
1
Marine College, Shandong University, Weihai 264209, China
2
Binzhou Institute of Scientific and Technical Information, Binzhou 256600, China
*
Author to whom correspondence should be addressed.
Fishes 2021, 6(3), 42; https://0-doi-org.brum.beds.ac.uk/10.3390/fishes6030042
Submission received: 16 August 2021 / Revised: 11 September 2021 / Accepted: 16 September 2021 / Published: 17 September 2021
Browse Figures
Review Reports Versions Notes

Abstract

:
Ruditapes philippinarum and Mercenaria mercenaria are economically important bivalve species. The complete ribosomal DNA (rDNA) unit sequences of R. philippinarum and M. mercenaria, with as-sembled rDNA unit lengths of 12,910 and 12,100 bp, respectively, were obtained in this study for the first time. The rDNA unit structural organisation was similar to that in other eukaryotes, in-cluding the following elements in order: 18S rRNA-internal transcribed spacer 1 (ITS1); 5.8S rRNA-ITS2-28S rRNA-intergenic spacer (IGS) (3′ external transcribed spacer (ETS); non-transcribed spacer (NTS)-5′ ETS). The genetic differences between R. philippinarum and M. mercenaria were mainly attributable to non-coding regions (ITS1, ITS2 and IGS), especially the IGS region. The boundaries of putative 3′ ETS, NTS and 5′ ETS were confirmed. Seven and three sub-repeat fragments were found in R. philippinarum and M. mercenaria, respectively. These frag-ments ranged from 4 to 154 bp in length, and were located at the NTS and 5′ ETS regions. Five and six cytosine–guanine (CpG) islands were detected in R. philippinarum and M. mercenaria, respec-tively, and these covered 85.58% and 79.29% of the entire IGS sequence, respectively. The phylo-genetic tree was constructed based on Veneridae ITS and 18S rRNA sequences using the maxi-mum likelihood (ML) method. The ML tree based on ITS revealed that species within the same genus clearly clustered together with relatively high supporting values, and all the genera were recovered as monophyletic. The phylogenetic analyses using 18S rRNA provided a weaker phy-logenetic signal than ITS.

Graphical Abstract

1. Introduction

Veneridae is one of the most diverse families of molluscs, with approximately 800 species [1,2]. Its members are distributed in marine environments worldwide, where they burrow in muddy or sandy habitats. Many Veneridae species are economically important due to their huge abundance in benthic environments, such as the Manila clam Ruditapes philippinarum, the northern quahog Mercenaria mercenaria, the western Pacific hard clam Meretrix petechialis, the equilateral clam Gomphina aequilatera and the Taca clam Protothaca theca [1,2,3]. The classification system and phylogenetic re-lationships of Veneridae have been based mainly on morphological characteristics and molecular data (mitochondrial and nuclear sequences) [1,2,4,5,6], which are currently un-der debate. Problems with existing classification methods are attributed to (1) the con-vergence of the evolutionary and phenotypic plasticity of morphological traits in most venerids [6]; (2) double uniparental inheritance (DUI) of mitochondrial DNA in some Veneridae species [7,8,9]; (3) highly conserved in some nuclear markers, such as the nu-clear protein-coding gene (H3) and ribosomal 18S or 28S rRNA [5].
In eukaryotes, nuclear ribosomal DNA (nrDNA) is usually organised in long tan-dem repeats that form nucleolar organising regions in chromosomes [10,11]. In animals, each nrDNA repeat unit has a similar structure, which consists of conserved coding genes (18S, 5.8S and 28S rRNA) separated by two internal transcribed spacers (ITS1 and ITS2) and an intergenic spacer (IGS) located between 28S rRNA and 18S rRNA [11]. The IGS region comprises two external transcribed spacers (5′ETS and 3′ETS) and a non-transcribed spacer (NTS) [11]. Each nrDNA region has evolved at a different rate and can be used as a genetic marker for phylogenetic analyses under different lev-els. Coding genes (18S and 28S rRNA) are highly conserved between organisms and are suitable for phylogenetic comparisons at the family and higher levels [12]. By con-trast, non-coding regions (ITS1, ITS2 and IGS) evolve rapidly and have high structural variability and can, thus, be used as genetic markers for phylogenetic comparisons at the genus and lower levels, including between individuals [12]. Salvi and Mariottini (2012) validated ITS2 as a suitable marker for venerid phylogenetic and taxonomic studies and showed a higher resolution below the subfamily level [5]. Despite recent studies of rDNA regions, only 11 complete ITS sequences from venerids are available in the NCBI database (Table 1). A complete rDNA unit and IGS region have never been reported.
The Manila clam Ruditapes philippinarum (Adams and Reeve, 1850) and quahog Mercenaria mercenaria (Linnaeus, 1758) are widespread and commercially important bivalve species. R. philippinarum originated in the Indo-Pacific region, but has been in-troduced to the west coast of North America and the Atlantic and Mediterranean coasts of Europe [13]. M. mercenaria is native to the east coast of the US and Canada and has been introduced to China for commercial cultivation [14]. Although R. philip-pinarum and M. mercenaria have been commercially exploited and deeply manipulated by human activities, their molecular genetics remain understudied [13,14]. The com-plete nrDNA repeat units of R. philippinarum and M. mercenaria are sequenced for the first time in this study. The two aims of this study are to analyse the structural features of each rDNA region, especially the IGS region, and to compare the phylogenetic rela-tionships with those of publicly available venerids based on ITS and 18S rRNA se-quences.

2. Materials and Methods

2.1. Sampling and DNA Extraction

Each individual of R. philippinarum and M. mercenaria was collected in May 2019 from Shuangdao Bay, Weihai, China. The adductor muscles were dissected and pre-served in 75% ethanol until further analysis. Total genomic DNA was extracted using an EasyPure® Marine Animal Genomic DNA Kit (TransGen Biotech Co., LTD., Beijing, China) according to the manufacturer’s instructions.

2.2. PCR Amplification, Cloning and Sequencing

The nrDNA fragments were amplified using the primers listed in Table 2. The primers used to amplify the 18S rRNA gene were designed based on the conserved se-quences of the 18S rRNA genes of R. philippinarum (EF426293), Meretrix meretrix (EF426291) and M. mercenaria (JN996711). Because the lengths of the complete 28S rRNA sequences in molluscs generally exceed 3.5 kb, the 28S rRNA genes of R. philip-pinarum and M. mercenaria were divided into three parts for primer design to ensure successful amplification. These primers were then designed based on the 28S rRNA gene sequences of M. mercenaria (KX713401), Mysia undata (KX713408) and R. decus-satus (KX713427). The primers used to amplify the 18S-ITS-28S region were designed based on the 18S rRNA genes of R. philippinarum (EF426293) and M. mercenaria (JN996711) and the 28S rRNA genes of M. mercenaria (KX713401) and M. undata (KX713408). The 28S-IGS-18S region was designed based on the newly sequenced 5′ end of 18S and 3′ end of 28S rRNA gene sequences from R. philippinarum and M. mer-cenaria, respectively. For PCR amplification, each reaction had a total volume of 25 μL and contained 1.0 μL of DNA template (100 ng/μl), 0.3 μL of each primer (10 μmol/L), 12.5 μL of 2× SanTaq PCR Mix (Sangon Biotech Co., Ltd., Shanghai, China) and 10.9 μL of ddH2O. The PCR amplification protocols were as follows: initial denaturation at 95 °C for 3 min; 35 cycles of denaturation at 95 °C for 30 s; annealing at 53–59 °C for 45 s; extension at 72 °C for 1 min–6 min 30 s (Table 2). All reactions included a final exten-sion at 72 °C for 7 min.
All PCR products were electrophoresed on 1.0% agarose gels and then extracted and purified using SanPrep Column DNA Gel Extraction Kit (Sangon Biotech Co., Ltd., Shanghai, China) according to the manufacturer’s instructions. As the lengths of 28S-IGS-18S region of R. philippinarum and M. mercenaria both exceeded 6 kb, which were hard to cloned successfully. The purified PCR products, except for the IGS region, were cloned as described by Guo et al. (2019) [11]. Three positive colonies were picked and sequenced by Sangon Biotech Co., Ltd. (Shanghai, China). The purified IGS region was quantified using the Qubit® 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA), and submitted to Sangon Biotech (Shanghai, China) for library construction and high-throughput sequencing on the Illumina HiSeq 2500 platform.

2.3. Sequence Analysis

The high-throughput sequencing reads were filtered by trimming based on low quality. Cleaned reads were assembled using Geneious v7.0 (Biomatters, Auckland, New Zealand). Annotations were performed, and the IGS region was confirmed using the BLAST tool from Geneious software and the NCBI database [15]. The complete rDNA unit was assembled using DNAMAN v10.0 [16]. The boundaries of each region were confirmed according to Guo et al. (2018) [11]. The regions were submitted to the GenBank database under the following accession numbers: MZ227551-MZ227552 (18S rRNA), MZ241537-MZ241538 (ITS), MZ227553-MZ227554 (28S rRNA) and MZ274272-MZ274273 (IGS). Multiple sequence alignment was performed using MAFFT online tool, and the parameters were set as default [17]. The general molecular features of R. philippinarum and M. mercenaria rDNA were calculated using MEGA 11 [18]. Sub-repeat (SR) fragments of the IGS region were detected using Tandem Repeats Finder software (https://tandem.bu.edu/trf/trf.html, accessed on 8 June 2021) [19]. The putative transcription initiation sites (TISs) and transcription termination sites (TTSs) were predicted based on other mollusc species [11,20,21]. CpG islands were confirmed with EMBOSS CpGPlot software (https://www.ebi.ac.uk/Tools/seqstats/emboss_cpgplot/, accessed on 2 June 2021) [17]. The sequence identity was confirmed using BioEdit v7.2.
For phylogenetic reconstruction, the complete ITS sequences (ITS1-5.8S-ITS2) and 18S rRNA sequences of R. philippinarum and M. mercenaria were determined. The re-maining sequences of Veneridae species were picked up from the GenBank database (Table 1). Perna viridis (MH279802 for ITS and KY081342 for 18S rRNA) was selected as the outgroup. Maximum likelihood (ML) phylogenetic trees were constructed with MEGA 11 using pairwise distances based on the GTR+G+I and GTR+G model of nucle-otide substitution for ITS and 18S rRNA, respectively. Bootstrapping (1000 replicates) was used to gauge support for individual clades in the resulting phylogenetic tree.

3. Results

3.1. Complete nrDNA Sequences of R. philippinarum and M. mercenaria

The complete rDNA units of R. philippinarum and M. mercenaria were determined after amplification, sequencing and assembly. In agreement with data from other eukaryotes, the structural organisation of the rDNA repeat units included the following elements in order: 18S rRNA-ITS1, 5.8S rRNA-ITS2 and 28S rRNA-IGS. The lengths of the complete rDNA repeat unit were 12,910 and 12,100 bp in R. philippinarum and M. mercenaria, respectively, and the respective GC contents were 59.14% and 57.52%, respectively.
The length, GC content, pairwise identity and variable sites of each region are presented in Table 3. The GC content of each region exceeded 51% in both R. philippinarum and M. mercenaria, and this value exceeded 63% in ITS1 and ITS2. The encoding genes (18S, 5.8S and 28S rRNA) in R. philippinarum and M. mercenaria had low levels of genetic difference, in addition to pairwise identities of 99.13%, 100% and 96.65%, respectively. A total of 12 and 90 variable sites were detected within 18S rRNA and 28S rRNA, respectively. The differences between R. philippinarum and M. mercenaria were observed mainly in non-coding regions (ITS1, ITS2 and IGS). The lengths of ITS1 and ITS2 had no evident differences between R. philippinarum and M. mercenaria; however, these regions contained 254 and 133 variable sites, respectively, which could have resulted from the relatively low sequence identities (46.24% in ITS1, 61.66% in ITS2). Each region of the rDNA units was compared, and the lowest pairwise identity (40.47%) and highest size difference (6216 bp in R. philippinarum vs. 5402 bp in M. mercenaria) were found in the IGS region. The contribution of variable sites of the IGS sequences between two species reached up to 43.66%.

3.2. Informative Characteristics of the IGS Region

The IGS region contains several functional elements, including TTS, TIS and SR fragments and CpG islands (Figure 1). A Thymine-rich tract was identified in the IGS region of R. philippinarum from bp 198 to 267, and M. mercenaria had two ploy (T) tracts at positions bp 191–217 and 274–292. These sequences were considered as putative TTSs. Therefore, the positions of 3′ ETS in the R. philippinarum and M. mercenaria IGS regions were at bp 1–267 and 1–292, respectively. The sequences TAAGAAGGCAGGCTACGGGCGGG and GGGAGGG in R. philippinarum and M. mercenaria were considered the putative TISs for RNA polymerase I, and a nucleotide at positions 3760 bp and 3742 bp, respectively, was identified as the + 1 position of the TIS. The 5′ ETS was located between the putative TIS and 5′ 18S rRNA, and the remainder of the IGS contained the NTS. The respective sizes of the 5′ ETS and NTS were 3492 and 2457 bp in R. philippinarum and 3449 and 1661 bp in M. mercenaria.
The DNA repeats in the IGS sequences of R. philippinarum and M. mercenaria were analysed, and the length, copy number, location, GC content and entropy of each sub-repeat fragment are listed in Table 4. Ruditapes philippinarum and M. mercenaria contained seven and three SRs, respectively, which ranged from 4 to 154 bp and were located at the NTS and 5′ ETS, respectively. The numbers and lengths of SR copies varied widely between R. philippinarum and M. mercenaria. In R. philippinarum, SR1 (length: 36 bp) was present in 3.8 copies; SR2 (length: 154 bp) was repeated three times; SR5 was found in only 1.9 repeats. By contrast, in M. mercenaria, SR1 and SR3 were both 48 bp in size and repeated three times. SR2 (96 bp) was present in 1.8 copies.
The nucleotide compositions of the IGS region had a GC bias, and the GC content was slightly higher in R. philippinarum (61.18%) than in M. mercenaria (58.33%). Five and six CpG islands were detected in R. philippinarum and M. mercenaria, respectively, and were both located at the NTS and 5′ ETS. The CpG islands covered 85.58% and 79.29% of the entire IGS sequences, respectively (Figure 1). The CpG islands in R. philippinarum ranged from 206 to 2540 bp in size. In M. mercenaria they varied from 215 to 1585 bp in size.

3.3. Phylogenetic Analysis of Veneridae Based on ITS and 18S rRNA Sequence

Due to the lack of sufficient IGS and 28S rRNA data, ITS and 18S rRNA sequences were used to reveal the phylogenetic relationships among Veneridae species. Specifically, the ML method was used to construct a phylogenetic tree (Figure 2). The ML tree based on ITS revealed that species within the same genus clearly clustered together with relatively high supporting values. At the generic level, all the genera were recovered as monophyletic. While at the subfamilial level, Chioninae and Meretricinae formed a monophyletic and well supported clade. Tapetinae was paraphyletic, P. undulata was more closely related to S. gigantea than genus Ruditapes. The phylogenetic tree based on 18S rRNA produced similar topologies with those of ITS except the position of R. decussatus and P. undulata, but with relatively low supporting values, which was in accordance with previous molecular studies [3,22].

4. Discussion

The complete rDNA units of R. philippinarum and M. mercenaria were sequenced for the first time in this study. The structural composition and order of the rDNA units were similar to those of other reported eukaryotes [10,11,12,20,21,23,24,25,26,27,28]. The rDNA unit lengths of the reported marine invertebrates ranged from 7.6 to 12.91 kb [11,20,21,23,24,25,26], and R. philippinarum had the longest rDNA (12.91 kb) among these taxa. The rDNA units of marine invertebrates were significantly shorter than those of mammals (43 kb in human, 45 kb in mouse) [27,28]. The differences in rDNA lengths between marine invertebrates were mainly attributable to the IGS, which varied in length from 1399 bp in Perna canaliculus [11] to 6216 bp in R. philippinarum.
Different copies of rDNA usually evolve in a concerted manner within species but independently between species [29]. However, if the mutation rate of a particular rDNA sequence occurs faster than that of homogenisation during evolution, multiple variations will cause rDNA polymorphisms, namely, non-concerted or incomplete concerted evolution [29,30,31]. rDNA polymorphism has been reported in Haliotis tuberculate [32], Pleuronichthys cornutus [29], Cynoglossus zanzibarensis [30], Soleidae [31] and Fusarium proliferatum [33]. A comparison of the ITS1 and ITS2 sequences between samples in this study and others that were submitted to the NCBI GenBank database led to the detection of three ITS1 types in R. philippinarum. Type one was the identity of the whole ITS1 sequence, which exceeded 99% with no mismatch fragments, e.g., MZ241537. Type two was a mismatch fragment detected from 485 to 517 bp, e.g., AY498755. Type three was a mismatch fragment detected from 485 to 524 bp, e.g., GU358316. The identities of the other type two and three sequences exceeded 99%. Except for ITS1 in R. philippinarum, no second ITS1 or ITS2-type was detected in R. philippinarum or M. mercenaria.
The 11 ITS sequences used to construct the phylogenetic tree belonged to five sub-families. Species assigned to the same sub-family were placed in the same clade except for P. undulata, which clustered with the genus Meretrix with a high supporting value (99%). This finding agreed with Salvi and Mariottini (2012) [5]. The mitochondrial genome analysis showed that P. undulata was clustered with other Tapetinae species [3,6]. Just four partial ITS sequences of P. undulata were deposited in the GenBank database (JN996788, JN996787, JN996750 and EU183531), and the identity of ITS1 and ITS2 of them were just 36.16% and 58.25%, respectively. More P. undulata samples should be cloned and sequenced to confirm the accuracy and phylogenetic position of the sequence.
A Thymine rich tract at the beginning of the 5′ IGS region could be considered as a putative TTS, a feature commonly found in marine invertebrates [11,20,21,23,24,25,26]. The typical structure of putative TTS was also observed in R. philippinarum and M. mercenaria. The core promoter contained TATA and GGGNGGG boxes, which appeared to be a general feature of TIS, and a AT-rich sequence has been reported upstream [11], but no TATA box and AT-rich sequence were detected in the sequences from R. philippinarum and M. mercenaria. Miscanthus sinensis, Sorghum bicolor, Fagus sylvatica and Quercus suber are also missing a canonical TATA box upstream of the putative TIS sites [34,35,36]. Precise transcription factor interactions are key to promoter recognition; the TATA box may not be a direct binding site for transcription factors, but is probably conserved only for DNA to be easily melted [37].
SR fragments of different sizes are a functional element of the IGS. Both intra- and inter-species variations in the copy number and length were observed for each SR type, and these variations were rarely predictable and did not demonstrate expected similarities [38]. Compared with each IGS region in R. philippinarum and M. mercenaria, the alignments showed that the relative conserved regions were at the beginning of the 3′ ETS and end of the 5′ ETS sequences. SR fragments were interspersed in the NTS and upstream region of the 5′ ETS, and heterogeneity in the IGS length could be attributed to the presence of SRs. IGS SR fragments have been demonstrated to serve as rRNA gene transcription regulatory elements, particularly as enhancers that increase the rRNA transcription rate [20,39]. Due to these differences in the SR sequences (including copy number, sequence and length) among different species, the IGS region could be a suitable genetic marker for revealing phylogenetic inter- and intra-species relationships and could be used to distinguish venerid species, as described previously in plant and other animal studies [35,37,40,41].

5. Conclusions

We performed, for the first time, a comparative analysis of the complete rDNA units of R. philippinarum and M. mercenaria. The structural organisation of these rDNA units was identical to that in other eukaryotes, with assembled rDNA unit lengths of 12,910 and 12,100 bp, respectively. Identified variations were mainly attributable to the non-coding regions, especially IGS. We also analysed the main features of IGS functional elements, namely, the TIS, TTS and SR and CpG islands. Our results suggest that R. philippinarum and M. mercenaria can provide a structural nuclear rDNA model for molecular comparisons, and this model may also contribute to analyses of Veneridae population genetics and evolution. Notably, the complete rDNA unit sequence was only determined from these two species, which likely limits the obtained phylogenetic information. Additional Veneridae species should be sequenced to clarify the phylogenetic relationships in this family.

Author Contributions

Conceptualization, Z.G.; Data curation, Z.W.; Methodology, Z.G.; Project administration, X.H.; Resources, X.H.; Software, Z.G. and Z.W.; Writing—original draft, Z.G. and Z.W.; Writing—review and editing, X.H. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Fundamental Research Funds for the Central Universities (2019ZRJC006).

Data Availability Statement

MDPI Research Data Policies.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Canapa, A.; Schiaparelli, S.; Marota, I.; Barucca, M. Molecular data from the 16S rRNA gene for the phylogeny of Veneridae (Mollusca: Bivalvia). Mar. Biol. 2003, 142, 1125–1130. [Google Scholar] [CrossRef]
  2. Kappner, I.; Bieler, R. Phylogeny of venus clams (Bivalvia: Venerinae) as inferred from nuclear and mitochondrial gene sequences. Mol. Phylogenet. Evol. 2006, 40, 317–331. [Google Scholar] [CrossRef]
  3. Chen, J.; Li, Q.; Kong, L.; Zheng, X. Molecular phylogeny of venus clams (Mollusca, Bivalvia, Veneridae) with emphasis on the systematic position of taxa along the coast of mainland China. Zool. Scr. 2011, 40, 260–271. [Google Scholar] [CrossRef]
  4. Mikkelsen, P.M.; Bieler, R.; Kappner, I.; Rawlings, T.A. Phylogeny of veneroidea (Mollusca: Bivalvia) based on morphology and molecules. Zool. J. Linn. Soc. 2006, 148, 439–521. [Google Scholar] [CrossRef] [Green Version]
  5. Salvi, D.; Mariottini, P. Molecular phylogenetics in 2D: ITS2 rRNA evolution and sequence-structure barcode from Veneridae to Bivalvia. Mol. Phylogenet. Evol. 2012, 65, 792–798. [Google Scholar] [CrossRef]
  6. Wang, Y.; Yang, Y.; Liu, H.Y.; Kong, L.F.; Yu, H.; Liu, S.K.; Li, Q. Phylogeny of Veneridae (Bivalvia) based on mitochondrial genomes. Zool. Scr. 2021, 50, 58–70. [Google Scholar] [CrossRef]
  7. Capt, C.; Bouvet, K.; Guerra, D.; Robicheau, B.M.; Stewart, D.T.; Pante, E.; Breton, S. Unorthodox features in two venerid bivalves with doubly uniparental inheritance of mitochondria. Sci. Rep. 2020, 10, 1087. [Google Scholar] [CrossRef]
  8. Chacón, G.M.; Arias-Pérez, A.; Freire, R.; Martínez, L.; Nóvoa, S.; Naveira, H.; Insua, A. Evidence of doubly uniparental inheritance of the mitochondrial DNA in Polititapes rhomboides (Bivalvia, Veneridae): Evolutionary and population genetic analysis of F and M mitotypes. J. Zool. Syst. Evol. Res. 2020, 58, 541–560. [Google Scholar] [CrossRef]
  9. Passamonti, M.; Plazzi, F. Doubly Uniparental Inheritance and beyond: The contribution of the Manila clam Ruditapes philippinarum. J. Zool. Syst. Evol. Res. 2020, 58, 529–540. [Google Scholar] [CrossRef]
  10. Dyomin, A.; Galkina, S.; Fillon, V.; Cauet, S.; Lopez-Roques, C.; Rodde, N.; Klopp, C.; Vignal, A.; Sokolovskaya, A.; Saiftdinova, A.; et al. Structure of the intergenic spacers in chicken ribosomal DNA. Genet. Sel. Evol. 2019, 51, 59. [Google Scholar] [CrossRef] [Green Version]
  11. Guo, Z.; Han, L.; Liang, Z.; Hou, X. Comparative analysis of the ribosomal DNA repeat unit (rDNA) of Perna viridis (Linnaeus, 1758) and Perna canaliculus (Gmelin, 1791). PeerJ 2019, 7, e7644. [Google Scholar] [CrossRef] [Green Version]
  12. Li, X.; Xu, J.; He, Y.; Shen, S.; Zhu, J.; Shen, Z. The complete nuclear ribosomal DNA (nrDNA) cistron sequence of Pyropia yezoensis (Bangiales, Rhodophyta). J. Appl. Phycol. 2016, 28, 663–669. [Google Scholar] [CrossRef]
  13. Chiesa, S.; Lucentini, L.; Freitas, R.; Marzano, F.N.; Minello, F.; Ferrari, C.; Filonzi, L.; Figueira, E.; Breda, S.; Baccarani, G.; et al. Genetic diversity of introduced Manila clam Ruditapes philippinarum populations inferred by 16S rDNA. Biochem. Syst. Ecol. 2014, 57, 52–59. [Google Scholar] [CrossRef]
  14. Hu, Z.; Song, H.; Zhou, C.; Yu, Z.L.; Yang, M.J.; Zhang, T. Complete mitochondrial genome of the hard clam (Mercenaria mercenaria). Mitochondrial DNA B 2019, 4, 3738–3739. [Google Scholar] [CrossRef]
  15. NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2017, 45, D12–D17. [Google Scholar] [CrossRef] [Green Version]
  16. Hao, X.Y.; Li, J.P.; Gao SQTuerxun, A.; Chang, X.C.; Hu, W.R.; Chen, G.; Huang, Q.S. SsPsaH, a H subunit of the photosystem I reaction center of Suaeda salsa, confers the capacity of osmotic adjustment in tobacco. Genes Genom. 2020, 42, 1455–1465. [Google Scholar] [CrossRef]
  17. Madeira, F.; Park, Y.M.; Lee, J.; Buso, N.; Gur, T.; Madhusoodanan, N.; Basutkar, P.; Tivey, A.R.N.; Potter, S.C.; Finn, R.D.; et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019, 47, W636–W641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  19. Benson, G. Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999, 27, 573–580. [Google Scholar] [CrossRef] [Green Version]
  20. Guo, Z.; Han, L.; Ding, Y.; Hou, X.; Liang, Z. Molecular characterization of the complete nuclear ribosomal DNA sequence of the blacklip abalone Haliotis rubra. N. Z. J. Mar. Fresh. 2018, 52, 430–443. [Google Scholar] [CrossRef]
  21. Guo, Z.; Hou, X.; Han, L. Complete nuclear ribosomal DNA sequence analyses of the black-footed abalone Haliotis iris. N. Z. J. Mar. Fresh. 2018, 52, 232–246. [Google Scholar] [CrossRef]
  22. Cheng, H.L.; Peng, Y.X.; Wang, F.; Meng, X.P.; Yan, B.L.; Dong, Z.G. Sequence analysis of 18S rRNA gene of six Veneridae clams (Mollusca: Bivalvia). J. Fish. Sci. Chin. 2008, 15, 559–567, (In Chinese with English abstract). [Google Scholar]
  23. Chae, J.; Seo, Y.; Yu, W.B.; Yoon, W.D.; Lee, H.E.; Chang, S.J.; Ki, J.S. Comprehensive analysis of the jellyfish Chrysaora pacifca (Goette, 1886) (Semaeostomeae: Pelagiidae) with description of the complete rDNA sequence. Zool. Stud. 2018, 57, 51. [Google Scholar]
  24. Guo, Z.; Ding, Y.; Zhang, X.; Hou, X. Complete nuclear ribosomal DNA sequence analysis of Pacific abalone Haliotis discus hannai. Fish. Sci. 2017, 83, 777–784. [Google Scholar] [CrossRef]
  25. Ki, J.S.; Kim, I.C.; Lee, J.S. Comparative analysis of nuclear ribosomal DNA from the moon jelly Aurelia sp.1 (Cnidaria: Scyphozoa) with characterizations of the 18S, 28S genes, and the intergenic spacer (IGS). Hydrobiologia 2009, 616, 229–239. [Google Scholar] [CrossRef]
  26. Ki, J.S.; Park, H.G.; Lee, J.S. Extensive analysis of nuclear cistron rDNA sequence of Paracyclopina nana (Cyclopoida: Cyclopettidae). Hydrobiologia 2011, 666, 3–9. [Google Scholar] [CrossRef]
  27. Gonzalez, I.L.; Sylvester, J.E. Complete sequence of the 43 kb human ribosomal DNA repeat: Analysis of the intergenic spacer. Genomics 1995, 27, 320–328. [Google Scholar] [CrossRef]
  28. Grozdanov, P.; Georgiev, O.; Karagyozov, L. Complete sequence of the 45-kb mouse ribosomal DNA repeat: Analysis of the intergenic spacer. Genomics 2003, 82, 637–643. [Google Scholar] [CrossRef]
  29. Yang, M.; Kong, X.Y.; Shi, W.; Gong, L.; Luo, H.R.; Wu, B.S. Remarkable sequence polymorphisms in 18S rDNA of Pleuronichthys cornutus (Pleuronectiformes: Pleuronectidae). Gene 2018, 677, 251–258. [Google Scholar] [CrossRef]
  30. Gong, L.; Shi, W.; Yang, M.; Si, L.Z.; Kong, X.Y. Non-concerted evolution in ribosomal ITS2 sequence in Cynoglossus zanzibarensis (Pleuronectiformes: Cynoglossidae). Biochem. Syst. Ecol. 2016, 66, 181–187. [Google Scholar] [CrossRef]
  31. Gong, L.; Shi, W.; Yang, M.; Kong, X.Y. Characterization of 18S-ITS1-5.8S rDNA in eleven species in Soleidae: Implications for phylogenetic analysis. Hydrobiologia 2018, 819, 161–175. [Google Scholar] [CrossRef]
  32. Van Wjormhoudt, A.; Gaume, B.; Le Bras, Y.; Roussel, V.; Huchette, S. Two different and functional nuclear rDNA genes in the abalone Haliotis tuberculata: Tissue differential expression. Genetica 2011, 139, 1217–1227. [Google Scholar] [CrossRef]
  33. Gálvez, L.; Clarkson, J.P.; Palmero, D. IGS region polymorphisms are responsible for failure of commonly used species-specific primers in Fusarium proliferatum isolates from diseased garlic. Plant. Pathol. 2020, 69, 713–722. [Google Scholar] [CrossRef]
  34. Chou, C.H.; Chiang, Y.C.; Chiang, T.Y. Within- and between-individual length heterogeneity of the rDNA IGS in Miscanthus sinensis var. glaber (Poaceae): Phylogenetic analyses. Genome 1999, 42, 1088–1093. [Google Scholar]
  35. Krawczyk, K.; Nobis, M.; Nowak, A.; Szczecińska, M.; Sawicki, J. Phylogenetic implications of nuclear rRNA IGS variation in Stipa L. (Poaceae). Sci. Rep. 2017, 7, 11506. [Google Scholar] [CrossRef] [Green Version]
  36. Inacio, V.; Rocheta, M.; Morais-Cecilio, L. Molecular organization of the 25S-18S rDNA IGS of Fagus sylvatica and Quercus suber: A comparative analysis. PLoS ONE 2014, 9, e98678. [Google Scholar] [CrossRef] [Green Version]
  37. Huang, Y.; Yu, F.; Li, X.; Luo, L.; Wu, J.; Yang, Y.; Deng, Z.; Chen, R.; Zhang, M. Comparative genetic analysis of the 45S rDNA intergenic spacers from three Saccharum species. PLoS ONE 2017, 12, e0183447. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Madani, M.; Ward, L.; Vierstraete, A.; De Boer, S.H.; Moens, M. The ribosomal intergenic spacer (IGS) in the potato and tobacco cyst nematodes, Globodera pallida, G. rostochiensis and G. tabacum. Mol. Cell. Probe. 2019, 48, 101441. [Google Scholar] [CrossRef]
  39. Kakou, B.; Angers, B.; Glémet, H. Extensive length variation in the ribosomal DNA intergenic spacer of yellow perch (Perca flavescens). Genome 2015, 59, 149–158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. He, Y.; Shen, S.D.; Shen, Z. Cloning and Application of the Complete Nuclear Ribosomal DNA (nrDNA) Cistron Sequence of Pyropia haitanensis (Bangiales, Rhodophyta). Bot. Mar. 2017, 60, 1515–1526. [Google Scholar] [CrossRef]
  41. Qiu, Y.Y.; Gao, Y.; Li, Y.; Ma, X.X.; Lv, Q.B.; Hu, Y.; Qiu, H.Y.; Chang, Q.C.; Wang, C.R. Comparative analyses of complete ribosomal DNA sequences of Clonorchis sinensis and Metorchis orientalis: IGS sequences may provide a novel genetic marker for intraspecific variation. Infect. Genet. Evol. 2020, 78, 104125. [Google Scholar] [CrossRef] [PubMed]
Fishes 06 00042 g001 550
Figure 1. Structural organisation scheme of the IGS region in R. philippinarum and M. mercenaria. ETS, external transcribed spacer; NTS, non-transcribed spacer; TTS, transcription termination site; TIS, transcription initiation site; CpG island, cytosine-guanine island.
Figure 1. Structural organisation scheme of the IGS region in R. philippinarum and M. mercenaria. ETS, external transcribed spacer; NTS, non-transcribed spacer; TTS, transcription termination site; TIS, transcription initiation site; CpG island, cytosine-guanine island.
Fishes 06 00042 g001
Fishes 06 00042 g002a 550Fishes 06 00042 g002b 550
Figure 2. Phylogenetic trees constructed using the maximum likelihood (ML) method based on ITS (A) and 18S rRNA sequences (B) of Veneridae species. The asterisk (★) means the sample collected in the present study.
Figure 2. Phylogenetic trees constructed using the maximum likelihood (ML) method based on ITS (A) and 18S rRNA sequences (B) of Veneridae species. The asterisk (★) means the sample collected in the present study.
Fishes 06 00042 g002aFishes 06 00042 g002b
Table
Table 1. Collection information, GenBank accession numbers and the lengths of each ribosomal region for Veneridae species used in the present study.
Table 1. Collection information, GenBank accession numbers and the lengths of each ribosomal region for Veneridae species used in the present study.
SubfamilyGenusSpeciesITS18S28S
GenBank Accession NumberLength (bp)GenBank Accession NumberLength (bp)GenBank Accession NumberLength (bp)
MeretricinaeMeretrixMeretrix lusoriaDQ3891081567JN9967141899MW221234739
MeretricinaeMeretrixMeretrix meretrixAY6958011480EF4262911900//
TapetinaeRuditapesRuditapes philippinarumDQ3994041303JN8073431833AM7797421476
TapetinaeRuditapesRuditapes variegatusJN8073811265JN8073491832DQ3438581388
TapetinaeRuditapesRuditapes decussatusHQ6341391259KX7133401779KX7134272078
TapetinaeParatapesPaphia undulataEU1835311236JN9967221831JQ277802743
ChioninaeLeukomaProtothaca stamineaEF0350841379AM7745701779AM7797441478
ChioninaeLeukomaProtothaca jedoensisDQ2202911437EF4262921831DQ3438561390
ChioninaeMercenariaMercenaria mercenariaDQ1904451332AF1060731804KX7134012085
CallistinaeSaxidomusSaxidomus giganteaEF0351101332///
CyclininaeCyclinaCyclina sinensisDQ1904461263EF4262891838DQ3438491396
outgroupPerna viridisMH279802735KY0813421698MK4191063679
Table
Table 2. PCR amplification primers.
Table 2. PCR amplification primers.
Amplified RegionPrimersSequences (5′–3′)Annealing TemperatureExtension TimeLength (bp)
18S rRNA18S-arTCAAATGTCTGCCCTATC53 °C for R. philippinarum,
55 °C for M. mercenaria		1 min 30 s1501–1573
18S-brTTCACCTACGGATACCTTG
18S rRNA–ITS-28S rRNAITS-arTAACAAGGTATCCGTAGGTG53 °C for R. philippinarum,
55 °C for M. mercenaria		1 min1436–1841
ITS-brCGTGCCAGTATTTAGCC
28S rRNA28S1-arAGTCGGGTTGTTTGGGAATG55 °C1 min1240–1251
28S1-brTTGATTCGGCAGGTGAGTTG
28S2-arCTGTGGGATGAACCAAACGC55 °C1 min 30 s1649–1691
28S2-brACCTTAGGACACCTGCGTTA
28S3-arTCACCCACTAATAGGGAACG55 °C1 min691–702
28S3-brAAGCACCTAAACCAAATGTC
IGS for R. philippinarumIGS-ar1CCAAATGCCTCGTCATCTAA59 °C6 min 30 s7829
IGS-br1CTGCCTTCCTTGGATGTG
IGS for M. mercenariaIGS-ar2GAATACAGACCGTGAAAGCG56 °C5 min 30 s6434
IGS-br1CTGCCTTCCTTGGATGTG
Table
Table 3. Characterisation of rDNA from R. philippinarum and M. mercenaria.
Table 3. Characterisation of rDNA from R. philippinarum and M. mercenaria.
RegionRuditapes philippinarumMercenaria mercenariaAlignment Length (bp)Pairwise Identity (%)Pairwise DistanceVariable Site (V)
Length (bp)GC Content (%)Length (bp)GC Content (%)
18S rRNA183251.85183051.75183399.130.00712
ITS163864.4267963.0469246.240.599254
5.8S rRNA15757.9615757.96157100.000.0000
ITS243465.6742065.7244661.660.432133
28S rRNA363357.68361257.22363896.650.02590
IGS621661.18540258.33623940.470.8782724
Table
Table 4. Characterisation of repeated regions defined for R. philippinarum and M. mercenaria.
Table 4. Characterisation of repeated regions defined for R. philippinarum and M. mercenaria.
SpeciesSub-RepeatWidth (bp)Copy NumberLocationGC Content (%)Entropy (0–2)
Ruditapes philippinarumR-S1363.8592–72755.651.99
R-S21543700–116154.551.98
R-S34292635–275275.001.49
R-S46223164–328463.641.90
R-S5721.93320–345465.751.89
R-S614124338–461965.961.91
R-S7721.94637–477165.751.89
Mercenaria mercenariaM-S14833615–376166.671.89
M-S2961.83599–377664.291.88
M-S34834054–419762.501.91
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, Z.; Wang, Z.; Hou, X. Comparative Analysis of the nrDNA Repeat Unit of Manila Clam Ruditapes philippinarum and Quahog Mercenaria mercenaria. Fishes 2021, 6, 42. https://0-doi-org.brum.beds.ac.uk/10.3390/fishes6030042

AMA Style

Guo Z, Wang Z, Hou X. Comparative Analysis of the nrDNA Repeat Unit of Manila Clam Ruditapes philippinarum and Quahog Mercenaria mercenaria. Fishes. 2021; 6(3):42. https://0-doi-org.brum.beds.ac.uk/10.3390/fishes6030042

Chicago/Turabian Style

Guo, Zhansheng, Zhen Wang, and Xuguang Hou. 2021. "Comparative Analysis of the nrDNA Repeat Unit of Manila Clam Ruditapes philippinarum and Quahog Mercenaria mercenaria" Fishes 6, no. 3: 42. https://0-doi-org.brum.beds.ac.uk/10.3390/fishes6030042

Article Metrics

Back to TopTop