- Open Access
Functional constraints and evolutionary dynamics of the repeats in the rDNA internal transcribed spacer 2 of members of the Anopheles barbirostris group
Parasites & Vectorsvolume 7, Article number: 106 (2014)
The Anopheles barbirostris group is widely distributed in Southeast Asia. Although seven species have been formally described, a molecular analysis of the rDNA ITS2 and the mitochondrial cytochrome oxidase I gene suggests that the group includes species that are morphologically very similar or identical.
We have previously shown that species in the Anopheles barbirostris Subgroup have an exceptionally large ITS2 (>1.5 kb), greater than in any other Anopheline group. However, the molecular processes responsible for generating such a large ITS2 have not previously been explored.
To determine the processes by which this large ITS2 is generated, we examined the sequence and secondary structure of the ITS2 of 51 specimens from five species of the Anopheles barbirostris Subgroup. These include the anthropophilic species An. campestris and three morphospecies of the Barbirostris Complex: An. vanderwulpi, An. barbirostris I and III, together with a previously undescribed member of this group (Clade IV).
Results and conclusions
All the specimens were found to have an ITS2 greater than 1.5 kb in length. The possibility that the spacer sequences amplified were pseudogenes was examined and discarded. The large size of ITS2 in the species studied is due to the presence of internal repeats of approximately 110 bp in length, confined to the central region of the spacer. Repeats varied markedly between the species examined, with respect to their organization, number and sequence similarity. The nucleotide diversity increased in direct relation to size variation and the presence of non-repeated elements.
A secondary structure analysis showed that the repeats form hairpin structures with a wide range of free energy values. These hairpin structures are known to facilitate the subsequent processing of mature rRNA. An analysis of the repeats from the different species suggests they originate from a common ancestor, with the repeats appearing before speciation of the Barbirostris Group.
In the ribosomal genome, the first and second internal transcribed spacers (ITS1 and ITS2) are located between the genes coding for 18S, 5.8S and 28S in the ribosomal DNA cluster. Although spacer sequences are not present in the mature ribosomal RNA molecules, they mediate the cleavage of the great RNA precursor  and form stable secondary structures . There is evidence that the deletion of ITS2 affects the maturation of both the small and large subunit rRNA . Whilst many mutations become fixed in a population through natural selection and genetic drift, the non-coding spacer regions in ITS2 evolve through a pattern of concerted evolution . Typically, ITS2 exhibits a low level of sequence variation within species and a high level of divergence between species. It is this feature that renders the region useful to infer phylogenetic relationships in different taxa, including insects . In Anopheles species, ITS2 has proved valuable for distinguishing members of sibling species complexes [6–8] and is the most commonly used marker for the design of species-specific PCR primers for species discrimination .
The Anopheles barbirostris group is distributed in the oriental region . Some of its members have been implicated in malaria and filariasis transmission in Southeast Asia [11–13]. Seven species have been formerly described: An. barbirostris, An. campestris, An. donaldi, An. pollicaris, An. franciscoi, An. hodgkini and An. vanderwulpi, the latter species being found in sympatry with An. barbirostris sensu stricto in Eastern Java. Anopheles campestris is considered the most anthropophilic of these species. Variations in habitat, resting behaviour and feeding preferences have been reported in the Anopheles barbirostris complex [16–18]. A molecular analysis of ITS2 and the cytochrome oxidase I (COI) gene region confirmed that the morphospecies An. barbirostris comprises at least three morphologically identical species, including An. vanderwulpi. In addition, one species, (Clade IV of Paredes-Esquivel, 2009), shares morphological characters with An. barbirostris and An. campestris.
All species in the Barbirostris Subgroup have a large ITS2 (>1.5 kb). The size varies between Clades, namely: Clade I, 1545 bp; Clade II, 1727 bp; Clade III, 1730 bp; Clade IV, 1583 bp and An. campestris, 1519 bp. ITS2 sequences of considerable size have been documented in other Anopheles species. Anopheles crucians species B has an ITS2 of 1012 bp . This length was the result of a complex array of repeats in different combinations. In Anopheles beklemishevi, a member of the Maculipennis Group, ITS2 comprises 638 bp, including two repeats of about 140 bp of similar though not identical sequence . Similarly, An. fluminensis from Bolivia has been reported to have three repeats of 125 bp each within an ITS2 of 596 bp . In Anopheles species of the genus Cellia, a series of internal repeats has led to a large ITS1. ITS2 of considerable length have also been seen in other insect species [8, 22, 23].
In this paper we establish that the outstanding size of ITS2 in the Barbirostris Subgroup is due to the presence of DNA insertions that occur as repeated units located in the central region of the internal transcribed spacer. These vary in copy numbers and organization. To determine whether these mutations affect the stability of the ITS2 and to discard the possibility that they were pseudogenes, i.e. non-functional relatives of genes, we analysed their secondary structure. We have carried out a comprehensive analysis of these repeated elements and discuss their origin in relation to the evolutionary history of the Barbirostris Group members.
The species examined were members of the Barbirostris Group studied by Paredes-Esquivel et al. (2009): Anopheles barbirostris clade I collected from Kalimantan (Indonesia) and Mae Hong Son (Thailand); An. vanderwulpi from the island of Sumatra (Indonesia); An. barbirostris clade III, collected from Mae Hong Son, Sa Kaeo, Tak and Kanchanaburi (Thailand); an unknown species clade IV, with mixed characters between An. barbirostris and An. campestris found in Sumatra and Trat and Sa Kaeo (Thailand) and finally An. campestris collected in Sa Kaeo. These species were identified based on the analysis of the COI, ITS2, combined with morphological examination of adult specimens. GenBank accession numbers, sequence sizes and place of collection are summarized in Table 1. Detailed DNA extraction and amplification procedures are described in Paredes-Esquivel et al. (2009).
Sequence alignment was complicated by the presence of internal repeats. BioEdit v 5.0.6. (Hall 1999) was used to align sequences manually and to determine GC content and length. Boundaries of the ITS2 region were identified in comparison to sequences from Anopheles gambiae (GenBank accession number X67157.1) . When the regions containing the repeats were removed, sequences could be aligned with ClustalW . The substitutional rate in the conserved 5,8S and 28 S regions was analysed in detail to discriminate functional genes from pseudogenes. Repeats were found using the Tandem repeats finder program: http://tandem.bu.edu/trf/trf.html, although an exhaustive visual examination was also required. Nucleotide diversity (π) was calculated with program DNASP5.101  to determine the degree of polymorphism within each type of repeat. This is defined as the average number of pairwise nucleotide substitutions divided by the length of the sequence (π = П/L) .
Sequences were folded using RNA Folding Form (version 2.3 energies), from the Mfold web server , located at: http://mfold.rna.albany.edu/?q=mfold/RNA-Folding-Form2.3. Default parameters (37°C with 5% suboptimal folding) were used to fold sequences. This program provides with several tentative secondary structures for a wide range of free energy values. Only the ITS2 region was included in the analysis.
BLAST searches were carried out to identify sequence similarities. A Bayesian analysis using Mr Bayes 3.2  was carried out to determine phylogenetic relationships among repeats. For this analysis, the Monte Carlo (MCMC) chain length involved 500,000 generations, with trees sampled every 100 generations. The analysis was repeated twice to confirm topologies. Posterior probabilities were employed to test statistical support for clades. Trees were visualized using the FigTree program http://tree.bio.ed.ac.uk/software/figtree/).
Sequence analysis and occurrence of repeat elements
The fifty-one specimens identified as five members of the Barbirostris Group included 47 haplotypes, thus the majority of haplotypes contained a single specimen. The G + C content varied between 52.3% and 55.9%, similar to that found in other Anopheles species. All specimens showed an exceptionally large ITS2 amplicon. The size of ITS2 in the members of the Barbirostris Group examined varied markedly in size: Clade I c. 1543 bp., Clade II 1727 bp, Clade-III had a similar ITS2 size ~ 1730 bp but differed in 30 fixed substitutions. In Anopheles campestris ITS2 was 1519 bp long whilst in the unnamed member of the Barbirostris Group, Clade IV, ITS2 was 1583 bp. (Table 1).
The exceptional length of the ITS2 region was found to be due to the presence of repeated elements, located at the centre of the ITS2, occupying 55% to 61% of the total length of the spacer. The number and organization of the repeats varied between species (Figure 1 and Table 2). In Clade I, eight repeats, organized in two groups comprising four repeats each, were present. There was no length variation in these repeats, being 112 and 108 bp long for types 1 and 2, respectively. Repeats in Anopheles campestris were arranged in a similar manner, except that the last repeat at the 3′ end turned into a non-repeated element. The presence of a single insertion/deletion event resulted in length variation in these repeats (Table 2). Anopheles vanderwulpi and Clade III showed a similar arrangement of nine repeated elements of three types. In both of these, sequences differed by several nucleotides. Non-repeated elements were also present in these species. Finally, the unknown species of the Barbirostris Group (Clade IV) contained five repeats of two types and varied sizes and three non-repeated elements (Figure 1). Most indels were located at the 3′ end of the repeat sequences in species where indels were found.
Table 2 shows the nucleotide diversity of the different type of repeats found in Barbirostris Group members. The nucleotide diversity π increased in direct relation to size variation and the presence of non-repeated elements. The lowest π values were observed in repeats of Clade I and An. campestris, whereas the unknown species Clade IV showed a high degree of polymorphism, particularly in type 2 repeats.
A GGGTGTG motif occurred at the 5′ end of most repeats. Although this motif was common to all species, in some cases variations of this motif were also observed (GGGTGGT in An. campestris, GGGTGGG in clade-I and An. vanderwulpi, GGGTGCG and TGGTGTG in the unknown member of the Barbirostris Subgroup). Using default parameters, most repeats formed long hairpin structures with a wide range of energy values. Results varied only in the number of hairpins that repeats formed, these being more evident in Clade I, Clade III and in An. campestris. At the base of the hairpin stem there was a motif GGGTGTG (or similar) at one side and a palindromic sequence on the other side of the stem. As an example Figure 2 shows the disposition of the hairpins in Clade I, at the lowest free energy value. Interestingly in this species each hairpin has a similar structure with three internal loops located close to the top of the hairpin (Figure 2). Stable RNA secondary structures are characterised by having low energy values. Energy values ranged from 593.6 to 629.4. Nevertheless, repeats formed at least 3 hairpin structures in all 25 resultant secondary structures, the only difference being the number of hairpins formed.
Among members of the Barbirostris Group, An. barbirostris Clade III has the largest ITS2 recorded so far in Anopheles species (~1730 bp). The large size of ITS2 in members of the Barbirostris Group (>1.5 kb) results from the presence of DNA insertions, comprising tandem-repeated elements of c.100 bp located in the centre of the spacer region. In the genus Anopheles the ITS2 is characterized by a very low level of polymorphism between species, even between cryptic species . Despite reports of intra-specific and intra-individual variation , ITS2 is still widely used to infer phylogenetic relations among species [30, 31]. In the Barbirostris Group, the presence of repeats did not affect tree topology, and this was consistent with results from COI analysis (10).
Given the high rate of substitution, we examined the possibility that the ITS2 of members of the Barbirostris Group were pseudogenes. This hypothesis was discarded for the following reasons: 1) There was a low rate of substitution in the flanking 5.8 and 28S regions which were amplified, these comprise 90 nucleotides of the alignment; 2) Hairpin structures were formed for all species in high number, at the lowest free energy values. Pseudogenes in ribosomal spacers tend to accumulate random substitutions at high rates, whereas functional ITS regions show many hairpins with compact stable secondary structures . Even when a point mutation was observed, it can be assumed that the structural integrity was not affected by the insertion of long repetitive elements, as large spacers seem to be less affected by mutations than shorter ones .
Functional constraints on repeats
Long ITS2 regions have been found in distant taxa including oomycetes  and amphipods . They have also been reported in several arthropod species, including rhipicephaline ticks  and insects of various orders [36, 37]. Repetitive DNA regions of considerable size have also been documented in the ITS2 of other Anopheles species [8, 19, 20]. Little has been written to explain the origin and function of repeated elements in ribosomal spacer regions. Secondary structure reconstructions show that repeats found in members of the Barbirostris Group form hairpin structures at a wide range of energy values. Since the Mfold server does not produce a single ‘correct folding’, many configurations obtained with the minimum free energy values were examined. It is known that hairpin structures are related to functional ITS regions , although their functionality has not been fully elucidated. As stated by Paskewitz et al. in their pioneering study of the secondary structure of ITS1 in the Anopheles gambiae complex, conservation of hairpin structures in a wide range of energy values may indicate the conservation of functional constraints, possibly related to the maturation of ribosomal subunits. Long stem-loop secondary structures are also formed by the repeats found in the ITS1 of members of the An. punctulatus group . However, the repeats found in the Barbirostris group showed no similarity to repeats in An. punctulatus or An. gambiae. Furthermore, a Blast search failed to reveal any similar sequence in GenBank. A possible explanation for this is that similarities to other insertional elements, e.g. transposable elements, may be obscured by the high rate of substitutions present in this spacer, as has been observed in other taxa .
Variation in other regions of the ribosomal gene has been reported as a result of adaptations to changes in the environment. Thus in the parthenogenetic Daphnia pulex, longer IGS regions provide these organisms with the plasticity required to adapt to different environments . In Bryophytes, repetitive sequences in the rDNA appear to result from exposure to heavy metals. Whether insertion elements or adaptation to new environments affect the ITS2 region of members of the Barbirostris Group has yet to be determined.
ITS2 in members of the Barbirostris Group is subject to a high rate of evolutionary change. We observed a high rate of substitution, even at an intraspecific level and incomplete homogenization of repeats. Nevertheless, when comparing the rate of substitutions at an interspecies level, ITS2 seems to evolve through a pattern of concerted evolution. Accordingly, repeated elements found were imperfect copies of an original one. This has also been reported in the long spacers of other taxa [22, 39, 40] and seems to be the result of replication slippage events .
Hairpin structures appear to facilitate subsequent processing . Thus rRNA processing in yeast requires sequences and/or higher order structures within ITS2, an example of which would be the extensive folding of ITS2 sequences, bringing into juxtaposition those regions of mature 5.8S and 25S rRNA that must interact during processing. Thus rRNA processing requires sequences and/or higher order structures within ITS2. One example may be the extensive folding of the ITS2 sequences (see Figure 2), bringing into juxtaposition those regions in mature 5.8S and 28S rRNA that must interact but that are at a distance in the primary sequence.
Nucleotide diversity varied considerably among repeat types. It was higher in repeats from An.vanderwulpi, Clade III and Clade IV, where non-repeated elements were observed and the length of the repeats varied considerably. On the other hand, the degree of polymorphism was much lower in An.campestris and particularly in Clade I, where indels were absent. Concerted evolution is a potent force in the formation of ITS2, but there are cases where concerted evolution is incomplete, as seen in Anopheles longirostris from Papua New Guinea . It is conceivable that concerted evolution is more effective in Clade I, leading to the stabilization of the repeats. This species is to date the most widely distributed member of the Anopheles barbirostris complex, having been reported in Thailand, Borneo (Indonesia), Vietnam  and in the islands of Sulawesi and Java, in the Indonesian archipelago . Indels did not affect hairpin structures, as these were located at the 3′ end, outside of hairpin structures.
We postulate that repeat elements were present in the common ancestor of the species of the Barbirostris Group, before members of the Group underwent subsequent speciation. This conclusion is based on two principal facts; firstly repeats in all species were flanked by a common motif (GGGTGTG and variants) and more importantly, the Bayesian analysis shows that type I repeats present in An. campestris form a monophyletic clade with their homologous repeats in Clade I.
We conclude that repeat elements were present in the common ancestor of species of the Barbirostris Group, before members of the Group underwent subsequent speciation. This conclusion is based on two principal observations; firstly repeats in all species are flanked by a common motif (GGGTGTG and variants) and more importantly, the Bayesian analysis shows that type I repeats present in An. campestris form a monophyletic clade with their homologous repeats in Clade I. Repeats form hairpin structures that are linked to the functional constraints on the internal transcribed spacer, since they are known to facilitate the processing of mature rRNA.
Cote CA, Greer CL, Peculis BA: Dynamic conformational model for the role of ITS2 in pre-rRNA processing in yeast. RNA. 2002, 8 (6): 786-797. 10.1017/S1355838202023063.
Coleman AW: ITS2 is a double-edged tool for eukaryote evolutionary comparisons. Trends Genet. 2003, 19 (7): 370-375. 10.1016/S0168-9525(03)00118-5.
Good L, Intine RVA, Nazar RN: Interdependence in the processing of ribosomal RNAs in Schizosaccharomyces pombe. J Mol Biol. 1997, 273 (4): 782-788. 10.1006/jmbi.1997.1351.
Dover G: Molecular drive - a cohesive mode of species evolution. Nature. 1982, 299 (5879): 111-117. 10.1038/299111a0.
Young I, Coleman AW: The advantages of the ITS2 region of the nuclear rDNA cistron for analysis of phylogenetic relationships of insects: a drosophila example. Mol Phyl Evol. 2004, 30 (1): 236-242. 10.1016/S1055-7903(03)00178-7.
Cornel AJ, Porter CH, Collins FH: Polymerase chain reaction species diagnostic assay for Anopheles quadrimaculatus cryptic species (Diptera: Culicidae) based on ribosomal DNA ITS2 sequences. J Med Ent. 1996, 33 (1): 109-116.
Manonmani A, Townson H, Adeniran T, Jambulingam P, Sahu S, Vijayakumar T: rDNA-ITS2 polymerase chain reaction assay for the sibling species of Anopheles fluviatilis. Acta Trop. 2001, 78 (1): 3-9. 10.1016/S0001-706X(00)00154-6.
Wilkerson RC, Reinert JF, Cong L: Ribosomal DNA ITS2 sequences differentiate six species in the Anopheles crucians complex (Diptera: Culicidae). J Med Ent. 2004, 41 (3): 392-401. 10.1603/0022-2585-41.3.392.
Loaiza JR, Bermingham E, Sanjur OI, Scott ME, Bickersmith SA, Conn JE: Review of genetic diversity in malaria vectors (Culicidae: Anophelinae). J Med Epid Inf Gen Evol Infect Dis. 2012, 12 (1): 1-12. 10.1016/j.meegid.2011.08.004.
Paredes-Esquivel C, Donnelly MJ, Harbach RE, Townson H: A molecular phylogeny of mosquitoes in the Anopheles barbirostris subgroup reveals cryptic species: implications for identification of disease vectors. Mol Phyl Evol. 2009, 50 (1): 141-151. 10.1016/j.ympev.2008.10.011.
Reid JA: Anophelines of Malaya and Borneo. Stud Instit Med Res Mal. 31:
Harrison B, Scanlon J: The subgenus Anopheles in Thailand (Diptera: Culicidae). Am Ent Instit. 1975, 12: 1-307. 1, pp iv +
Apiwathnasorn C, Prommongkol S, Samung Y, Limrat D, Rojruthai B: Potential for Anopheles campestris (Diptera : Culicidae) to transmit malaria parasites in Pa Rai subdistrict (Aranyaprathet, Sa Kaeo Province), Thailand. J Med Ent. 2002, 39 (4): 583-586. 10.1603/0022-2585-39.4.583.
Reid JA: The Anopheles barbirostris group (Diptera, Culicidae). Bull Ent Res. 1962, 53: 1-57. 10.1017/S0007485300047945.
Townson H, Dyer N, McAlister E, Satoto TBT, Bangs MJ, Harbach RE: Systematics of Anopheles barbirostris van der Wulp and a sibling species of the barbirostris complex (Diptera: Culicidae) in eastern Java, Indonesia. Syst Ent. 2013, 38 (1): 180-191. 10.1111/j.1365-3113.2012.00653.x.
Lien JC, Kawengian BA, Partono F, Lami B, Cross JH: Brief survey of mosquitos of South Sulawesi, Indonesia, with special reference to identity of Anopheles barbirostris (Diptera-Culicidae) from Margolembo area. J Med Ent. 1977, 13 (6): 719-727.
Reid JA, Harrison BA, Atmosoedjono S: Variation and vector status in Anopheles barbirostris. Mosquito Systematics. 1979, 11: 235-250.
Fu Q, Matsumoto Y, Matsumura M, Hirai Y, Sato Y, Noda H: Presence of a short repeat sequence in internal transcribed spacer (ITS) 1 of the rRNA gene of Sogatella furcifera (Hemiptera: Delphacidae) from geographically different populations in Asia. Appl Ent Zool. 2012, 47 (2): 95-101. 10.1007/s13355-012-0093-y.
Kampen H: The ITS2 ribosomal DNA of Anopheles beklemishevi and further remarks on the phylogenetic relationships within the Anopheles maculipennis group of species (Diptera : Culicidae). Parasit Res. 2005, 97 (2): 118-128. 10.1007/s00436-005-1393-8.
Brelsfoard CL, Fritz GN, Rodriguez R: Sequence analysis of the rDNA internal transcribed spacer 2 and polymerase chain reaction identification of Anopheles fluminensis (Diptera: Culicidae: Anopheles) in Bolivia. J Med Ent. 2006, 43 (3): 460-466. 10.1603/0022-2585(2006)43[460:SAOTRI]2.0.CO;2.
Bower JE, Cooper RD, Beebe NW: Internal repetition and intraindividual variation in the rDNA ITS1 of the Anopheles punctulatus group (Diptera: Culicidae): multiple units and rates of turnover. J Mol Evol. 2009, 68 (1): 66-79. 10.1007/s00239-008-9188-z.
Murrell A, Campbell NJH, Barker SC: Recurrent gains and losses of large (84–109 bp) repeats in the rDNA internal transcribed spacer 2 (ITS2) of rhipicephaline ticks. Insect Mol Biol. 2001, 10 (6): 587-596. 10.1046/j.0962-1075.2001.00298.x.
Hlinka O, Murrell A, Barker SC: Evolution of the secondary structure of the rRNA internal transcribed spacer 2 (ITS2) in hard ticks (Ixodidae, Arthropoda). Heredity. 2002, 88: 275-279. 10.1038/sj.hdy.6800040.
Paskewitz SM, Wesson DM, Collins FH: The internal transcribed spacers of ribosomal DNA in five members of the anopheles gambiae species complex. Insect Mol Biol. 1993, 2 (4): 247-257.
Thompson JD, Higgins DG, Gibson TJ: CLUSTALW: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.
Rozas J, Librado P, Sanchez-Del B, Messeguer X, Rozas R: DnaSP version 5.10.1, 4 March 2010: an integrated program for molecular population genetics and molecular evolution analysis. [http://www.ub.edu/dnasp/]
Nei M, Li WH: Mathematical-model for studying genetic-variation in terms of restriction endonucleases. Proc Natl Acad Sci U S A. 1979, 76 (10): 5269-5273. 10.1073/pnas.76.10.5269.
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.
Paredes-Esquivel C, Harbach RE, Townson H: Molecular taxonomy of members of the Anopheles hyrcanus group from Thailand and Indonesia. Med Vet Ent. 2011, 25 (3): 348-352. 10.1111/j.1365-2915.2010.00937.x.
Ndo C, Simard F, Kengne P, Awono-Ambene P, Morlais I, Sharakhov I, Fontenille D, Antonio-Nkondjio C: Cryptic genetic diversity within the Anopheles nili group of malaria vectors in the equatorial forest area of Cameroon (central Africa). PLoS One. 2013, 8 (3): e58862-10.1371/journal.pone.0058862. doi:10.1371/journal.pone.0058862
Buckler ES, Ippolito A, Holtsford TP: The evolution of ribosomal DNA: divergent paralogues and phylogenetic implications. Genetics. 1997, 145 (3): 821-832.
Alquezar DE, Hemmerter S, Cooper RD, Beebe NW: Incomplete concerted evolution and reproductive isolation at the rDNA locus uncovers nine cryptic species within Anopheles longirostris from Papua New Guinea. BMC Evol Biol. 2010, 10: 392-10.1186/1471-2148-10-392. doi:10.1186/1471-2148-10-392
Thines M: Characterisation and phylogeny of repeated elements giving rise to exceptional length of ITS2 in several downy mildew genera (Peronosporaceae). Fungal Genet Biol. 2007, 44 (3): 199-207. 10.1016/j.fgb.2006.08.002.
Kornobis E, Palsson S: The ITS region of groundwater amphipods: length, secondary structure and phylogenetic information content in Crangonyctoids and Niphargids. J Zool Syst & Evol Res. 2013, 51 (1): 19-28. 10.1111/jzs.12006.
von der Schulenburg JHG, Hancock JM, Pagnamenta A, Sloggett JJ, Majerus MEN, Hurst GDD: Extreme length and length variation in the first ribosomal internal transcribed spacer of ladybird beetles (Coleoptera: Coccinellidae). Mol Biol Evol. 2001, 18 (4): 648-660. 10.1093/oxfordjournals.molbev.a003845.
Oh HK, Yoon HJ, Kim MJ, Jeong HU, Kim SR, Hwang JS, Bae CH, Kim I: ITS2 Ribosomal DNA sequence variation of the bumblebee, Bombus ardens (Hymenoptera: Apidae). Genes & Genomics. 2009, 31 (4): 293-303. 10.1007/BF03191202.
Gorokhova E, Dowling TE, Weider LJ, Crease TJ, Elser JJ: Functional and ecological significance of rDNA intergenic spacer variation in a clonal organism under divergent selection for production rate. Proc of Roy Soc London Ser B-Biol Sci. 2002, 269 (1507): 2373-2379. 10.1098/rspb.2002.2145.
Thines M, Komjati H, Spring O: Exceptional length of ITS in Plasmopara halstedii is due to multiple repetitions in the ITS2 region. Eur J Plant Pathol. 2005, 112 (4): 395-398. 10.1007/s10658-005-6606-6.
Choi YJ, Hong SB, Shin HD: Extreme size and sequence variation in the ITS rDNA of Bremia lactucae. Mycopathologia. 2007, 163 (2): 91-95. 10.1007/s11046-007-0092-7.
Satoto TBT: Cryptic species within Anopheles barbirostris Van der Wulp, 1884 inferred from nuclear and mitochondrial gene sequence variation (Ph.D thesis). 2001, University of Liverpool
The authors declare that they have no competing interests.
CPE carried out the molecular studies and prepared draft Figures and Tables. HT assisted in the analysis of the data and their presentation. CPE provided an initial draft of the text. HT and CPE contributed equally to the final article and both authors read and approved the final manuscript.