- Open Access
Functional constraints and evolutionary dynamics of the repeats in the rDNA internal transcribed spacer 2 of members of the Anopheles barbirostris group
© Paredes-Esquivel and Townson; licensee BioMed Central Ltd. 2014
Received: 29 January 2014
Accepted: 7 March 2014
Published: 19 March 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.
List of specimens and size of ITS2 (bp)
GenBank Acc N°
Mae Hong Son (Thailand)
Mae Hong Son (Thailand)
Sa Kaeo (Thailand)
Sa Kaeo (Thailand)
Sa Kaeo (Thailand)
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).
Nucleotide diversity and sequence variation in repeats
N° of SS
N° sites compared*
Clade I type 1
A, C, E, G
Clade I type 2
D, B, F, H
An. campestris type 1
A, C, E, G
110 (G), 111 (A,C,E)
An. campestris type 2
D, B, F
112 (F), 113 (D,B)
An. vanderwulpi type 1
A, E, H
103 (E), 108 (H), 117 (A)
An. vanderwulpi type 2
C, D, G
102 (D), 105 (G), 107 (C)
An. vanderwulpi type 3
B, F, I
99 (B,F), 110 (I)
Clade III type 1
A, E, H
107 (E), 113 (H), 116 (A)
Clade III type 2
C, D, G
103 (D,G), 105 (C)
Clade III type 3
B, F, I
99 (B), 101 (F), 113 (I)
Clade IV type 1
A, D, F
97 (F), 98 (D), 116 (A)
Clade IV type 2
114 (E), 119 (G)
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.
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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Dover G: Molecular drive - a cohesive mode of species evolution. Nature. 1982, 299 (5879): 111-117. 10.1038/299111a0.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Reid JA: Anophelines of Malaya and Borneo. Stud Instit Med Res Mal. 31:Google Scholar
- Harrison B, Scanlon J: The subgenus Anopheles in Thailand (Diptera: Culicidae). Am Ent Instit. 1975, 12: 1-307. 1, pp iv +Google Scholar
- 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.View ArticleGoogle Scholar
- Reid JA: The Anopheles barbirostris group (Diptera, Culicidae). Bull Ent Res. 1962, 53: 1-57. 10.1017/S0007485300047945.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Reid JA, Harrison BA, Atmosoedjono S: Variation and vector status in Anopheles barbirostris. Mosquito Systematics. 1979, 11: 235-250.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.PubMed CentralView ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.0058862PubMed CentralView ArticlePubMedGoogle Scholar
- Buckler ES, Ippolito A, Holtsford TP: The evolution of ribosomal DNA: divergent paralogues and phylogenetic implications. Genetics. 1997, 145 (3): 821-832.PubMedGoogle Scholar
- 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-392PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 LiverpoolGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.