Short CommunicationMolecular phylogenetics in 2D: ITS2 rRNA evolution and sequence-structure barcode from Veneridae to Bivalvia
Graphical abstract
Highlights
► We analyzed ITS2 rRNA sequence and secondary structure in Veneridae and Bivalvia. ► Phylogenetic trees based on ITS2 sequences and on sequence-structure are congruent. ► ITS2 sequence-structure approach is suitable for phylogeny and taxonomy of venerids. ► We discuss important sequence-structure features of the ITS2 folding unique to venerids and Bivalvia.
Introduction
Molecular phylogenetics and molecular taxonomy have traditionally focused on analyzing DNA (or RNA) sequences for phylogenetic and taxonomic assessments. However, in the last few years the exploitation of secondary structure information of ribosomal RNA (rRNA) molecules such as the nuclear ribosomal internal transcribed spacer 2 (ITS2) has revealed a promising approach not only in phylogenetic reconstruction (Coleman, 2003, Schultz et al., 2006, Schultz and Wolf, 2009, Salvi et al., 2010) but also in species diagnosis (Coleman, 2003, Coleman, 2009, Müller et al., 2007).
The ITS2 is generally organized in four main helix domains (DI–IV) of secondary structure and the combined information of both folding and primary sequence has been exploited for phylogenetic inference, taxonomic classification, and species delimitation, providing an increased phylogenetic resolution and reliability than approaches based on primary sequence data alone (Coleman, 2003, Schultz et al., 2006, Bologna et al., 2008, Song et al., 2008, Schultz and Wolf, 2009, Keller et al., 2010, Salvi et al., 2010). Compensatory base changes (CBCs) – mutations in both nucleotides of a paired position in a double-stranded structure of the transcribed RNA – in conserved regions of the eukaryote ITS2 sequence-structure have been shown to correlate with interbreeding incompatibility between species. Several studies on animal, plants, and fungi demonstrated that the presence of at least one CBC in ITS2 secondary structures is a good classifier (reliability higher than 90%) of two organisms belonging to distinct species (Coleman, 2000, Coleman, 2003, Coleman, 2009, Müller et al., 2007). Moreover, some conserved features of the ITS2 secondary structure such as stem-loops domains have shown to be diagnostic of higher taxonomic grouping such as tribes, subfamilies, families, and orders (Oliverio et al., 2002, Bologna et al., 2008, Keller et al., 2008, Salvi et al., 2010). Thus, the information from the ITS2 folding and sequence-structure variation is a useful tool in molecular taxonomy for species and higher-taxa diagnosis.
In this study we analyzed the primary sequence and the secondary structure information from the nuclear ITS2 rRNA in the phylogenetic context of the family Veneridae. This family is a large group of bivalves with approximately 800 species, with a largely unresolved phylogenetic history and unstable taxonomy despite recent molecular and morphological studies (Canapa et al., 1996, Canapa et al., 2003, Kappner and Bieler, 2006, Mikkelsen et al., 2006, Chen et al., 2011). Indeed, while on the one hand the use of morphological traits is problematic since most of them are either uniform or homoplastic (Mikkelsen et al., 2006), on the other hand the use of molecular markers can be problematic due to (i) double uniparental inheritance of mitochondrial genomes in some venerids (Mikkelsen et al., 2006); (ii) inherent difficulties in amplifying and sequencing “the barcoding gene” COI in venerids (Kappner and Bieler, 2006), and (iii) the very low resolution of nuclear markers such as the H3 gene and the ribosomal 18S/28S rRNAs (Chen et al., 2011). To overcome these problems, the use of an extremely variable and easy to amplify nuclear marker such as the ITS2 is highly advisable, but despite a promising test on four species (Cheng et al., 2006), ITS2 has never been used in phylogenetic and taxonomic studies of Veneridae.
The main aim of this study was to integrate the information of individual RNA secondary structures with primary sequences in venerids phylogeny and taxonomy, by employing a combined model of sequence-structure evolution as well as identifying conserved elements of the ITS2 folding. By including several representatives from Bivalvia subclasses we further test the phylogenetic resolution of the ITS2 marker and its suitability as a barcode tool for molecular taxonomy of venerid species to higher taxonomic levels among bivalves.
Section snippets
Sampling, DNA extraction and sequencing
ITS2 sequences were obtained from live specimens (60%), frozen specimens from trade (4%), or retrieved from Genbank (36%). A total of 45 specimens from 24 venerid species were employed in the molecular analyses. Twenty additional Bivalvia ITS2 sequences from the subclasses Heterodonta, Paleoheterodonta, and Pteriomorphia were also analyzed (see Supplementary Table S1 for taxonomic references, locality data, and Genbank accession numbers). Taxonomic nomenclature in the text, tables and figures
ITS2 sequence variation
The venerid ITS2 sequences ranged in length from 241 (Saxidomus gigantea) to 417 (Meretrix meretrix) base pairs (bp). For the eight species for which multiple specimens were available no differences in ITS2 sequences length were observed except in Ruditapes philippinarum. In this species, specimens from four different localities showed ITS2 sequences length from 310 to 379 bp with an average length of 363 bp. Intra-individual variation of the ITS2 sequences was observed only in two species,
Discussion
The suitability of the ITS2 as a barcode marker for venerids is highly supported by the present study as assessed against several criteria: (1) the presence of a barcode gap between inter- and intra-specific variation; (2) the clustering of conspecific sequences in distinct groups with high bootstrap supports; and (3) the presence of sequence-structure diagnostic character differences between congeneric species. Based on the ITS2 sequence-structure analyses, it was possible to identify all the
Acknowledgments
We thank the Associate Editor Suzanne Williams, Mattias Wolf and Marco Oliverio for useful suggestions and people from “Gruppo Malacologico Mediterraneo” (Rome, Italy) for their help during the sampling. PM wishes to thank the University of Roma Tre for financial support. DS is supported by the FCT post-doctoral grant SFRH/BPD/66592/2009 and by the SYNTHESYS Project http://www.synthesys.info/ which is financed by European Community Research Infrastructure Action under the FP7 Capacities
References (33)
- et al.
Phylogeny and evolutionary history of the blister beetles (Coleoptera, Meloidae)
Mol. Phylogenet. Evol.
(2008) - et al.
A Molecular Phylogeny of Heterodonta (Bivalvia) Based on Small Ribosomal Subunit RNA Sequences
Mol. Phylogenet. Evol.
(2001) - et al.
Study on sequences of ribosomal DNA internal transcribed spacers of clams belonging to the Veneridae family (Mollusca: Bivalvia)
Acta Genet. Sinic.
(2006) The significance of a coincidence between evolutionary landmarks found in mating affinity and a DNA sequence
Protist
(2000)ITS2 is a double-edged tool for eukaryote evolutionary comparisons
Trends Genet.
(2003)Is there a molecular key to the level of “biological species” in eukaryotes? A DNA guide
Mol. Phylogenet. Evol.
(2009)- et al.
Phylogeny of venus clams (Bivalvia: Venerinae) as inferred from nuclear and mitochondrial gene sequences
Mol. Phylogenet. Evol.
(2006) - et al.
Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure
J. Mol. Biol.
(1999) - et al.
ITS2 rRNA evolution and its congruence with the phylogeny of muricid neogastropods (Caenogastropoda, Muricoidea)
Mol. Phylogenet. Evol.
(2002) - et al.
The analysis of rRNA sequence-structure in phylogenetics: an application to the family pectinidae (Mollusca: Bivalvia)
Mol. Phylogenet. Evol.
(2010)
ITS2 sequence-structure analysis in phylogenetics: a how-to manual for molecular systematics
Mol. Phylogenet. Evol.
A molecular phylogeny of the bivalve mollusks
Mol. Biol. Evol.
Phylogenetic analysis of Veneridae (Bivalvia): comparison of molecular and paleontological data
J. Mol. Evol.
Molecular data from the 16S rRNA gene for the phylogeny of Veneridae (Mollusca: Bivalvia)
Mar. Biol.
Molecular phylogeny of venus clams (Mollusca, Bivalvia, Veneridae) with emphasis on the systematic position of taxa along the coast of mainland China
Zool. Scripta
MrBayes: Bayesian inference of phylogenetic trees
Bioinformatics
Cited by (31)
Sperm ultrastructure of Ruditapes variegata and Tapes literatus (Mollusca, Bivalvia, Veneridae, Tapetinae) from Pescadores, Taiwan
2021, Tissue and CellCitation Excerpt :Tapes literatus, a similar shell morphology, is another common Tapetinae species in Taiwan but little research had been reported, despite this species is widely spread in the Indo-West Pacific (Chen et al., 2009, 2010, 2011; Simone and Wilkinson, 2008). Several genetic analyses of the three species did not support the current congeneric assignment of R. variegata and R. philippinarum (Chen et al., 2010, 2011; Salvi and Mariottini, 2012; Aranda-Burgos et al., 2014). The phylogenetic molecular datasets of R. variegata and R. philippinarum, plus Tapes, were well resolved into three major clades and supported robustly by high bootstrapping values (Salvi et al., 2010; Salvi and Mariottini, 2012; Morton and Leung, 2018).
Evaluation of ITS2 molecular morphometrics effectiveness in species delimitation of Ascomycota – A pilot study
2019, Fungal BiologyCitation Excerpt :Hence, documenting the variations (evolutionary events) in the ITS2 sequences and their distribution in the secondary structure was attempted to define species boundaries. Evolutionary events had been documented to be phylogenetically informative (Fritz et al., 1994; Caisová et al., 2011; Nagy et al., 2012; Salvi and Mariottini, 2012). In our dataset, among substitution events, the frequency of transition was observed to be higher than the frequency of transversion.
An Antarctic flock under the Thorson's rule: Diversity and larval development of Antarctic Velutinidae (Mollusca: Gastropoda)
2019, Molecular Phylogenetics and EvolutionCitation Excerpt :To infer the species tree we used the multi-species coalescent model implemented in the *BEAST extension (Heled and Drummond, 2010) of the BEAST package. ITS2 has proven to be a valuable marker for mollusc phylogenetics and taxonomy (e.g. Oliverio et al., 2002; Puillandre et al., 2011), especially when the information from both the sequence and the secondary structure are exploited (Salvi et al., 2010; Salvi et al., 2014; Salvi and Mariottini, 2012, 2017). Including RNA secondary structures improves accuracy and robustness in reconstructing phylogenetic trees (e.g. Keller et al., 2010).
A consensus secondary structure of ITS2 for the diatom Order Cymatosirales (Mediophyceae, Bacillariophyta) and reappraisal of the order based on DNA, morphology, and reproduction
2018, Molecular Phylogenetics and EvolutionCitation Excerpt :In spite of the highly divergent primary sequences of the fast-evolving ITS2, the use of its sequence-structure alignment and resulting phylogenies has been increasing, at least in part due to current advances in bioinformatics (Keller et al., 2010; Wolf et al., 2014). This approach has contributed to disentangling the complicated evolutionary relationships in many eukaryotic groups, either solely or in combination with other molecular markers (Biffin et al., 2007; Buchheim et al., 2012, 2011; Caisová et al., 2013, 2011; Coleman and Vacquier, 2002; Heeg and Wolf, 2015; Merget and Wolf, 2010; Miao et al., 2008; Salvi et al., 2014; Salvi and Mariottini, 2012; Schoch et al., 2012; Wiemers et al., 2009). So far, in diatoms such an approach has only targeted troublesome species or genera (Kaczmarska et al., 2014; Lim et al., 2018; Sorhannus et al., 2010; Wolf, 2004).
Does GenBank provide a reliable DNA barcode reference to identify small alien oysters invading the Mediterranean Sea?
2015, Journal of the Marine Biological Association of the United KingdomPhylogenomic resolution of Imparidentia (Mollusca: Bivalvia) diversification through mitochondrial genomes
2023, Marine Life Science and Technology