Secondary structure models of the nuclear internal transcribed spacer regions and 5.8S rRNA in Calciodinelloideae (Peridiniaceae) and other dinoflagellates

Secondary structure models of the 5.8S rRNA and both internal transcribed spacers (ITS1 and ITS2) are proposed for Calciodinelloideae (Peridiniaceae) and are also plausible for other dinoflagellates. The secondary structure of the 5.8S rRNA corresponds to previously developed models, with two internal paired regions and at least one 5.8S rRNA–28S rRNA interaction. A general secondary structure model of ITS1 for Calciodinelloideae (and other dinoflagellates), consisting of an open multibranch loop with three major helices, is proposed. The homology of these paired regions with those found in other taxa, published in previous studies (e.g. yeast, green algae and Platyhelmithes) remains to be determined. Finally, a general secondary structure model of ITS2 for Calciodinelloideae (and other dinoflagellates) is reconstructed. Based on the 5.8S rRNA–28S rRNA interaction, it consists of a closed multibranch loop, with four major helices. At least helix III and IV have homology with paired regions found in other eukaryotic taxa (e.g. yeast, green algae and vertebrates). Since the secondary structures of both ITS regions are more conserved than the nucleotide sequences, their analysis helps in understanding molecular evolution and increases the number of structural characters. Thus, the structure models developed in this study may be generally useful for future phylogenetic analyses.     

Secondary structure and phylogenetic analysis of the internal transcribed spacers 1 and 2 of bush crickets (Orthoptera: Tettigoniidae: Barbitistini)

We inferred secondary structure models of the internal transcribed spacers (ITS) 1 and 2 of bush crickets using a combined comparative and thermodynamic approach. The inferred secondary structure models were used to account for interdependency of interacting nucleotides in a phylogenetic analysis of the bush cricket genus Poecilimon. Our analysis indicates that the two previously reported conformational structures (i.e., hairpin and ring) of ITS2 are likely to fold in bush crickets as well and that both predicted structures are similar to those proposed for other eukaryotes. Comparing predicted ITS1 secondary structure models proved to be difficult because of substantial variation in their nucleotide sequence length. Our study revealed that the phylogenetic signal of ITS1 and ITS2 is largely congruent with that preserved in the mitochondrial genes 16S rRNA, tRNA‐Val and 12S rRNA. The phylogenetic signal in both the nuclear and the mitochondrial genome question the monophyly of the genus Poecilimon: species of the genera Poecilimonella, Parapoecilimon, Polysarcus and Phonochorion consistently cluster within Poecilimon.     

Molecular Evolution and Phylogenetic Utility of the Internal Transcribed Spacer 2 (ITS2) in Calyptratae (Diptera: Brachycera)

The resolution potential of internal transcribed spacer 2 (ITS2) at deeper levels remains controversial. In this study, 105 ITS2 sequences of 55 species in Calyptratae were analyzed to examine the phylogenetic utility of the spacer above the subfamily level and to further understand its evolutionary characteristics. We predicted the secondary structure of each sequence using the minimum-energy algorithm and constructed two data matrixes for phylogenetic analysis. The ITS2 regions of Calyptratae display strong A-T bias and slight variation in length. The tandem and dispersed repeats embedded in the spacers possibly resulted from replication slippage or transposition. Most foldings conformed to the four-domain model. Sequence comparison in combination with the secondary structures revealed six conserved motifs. Covariation analysis from the conserved motifs indicated that the secondary structure restrains the sequence evolution of the spacer. The deep-level phylogeny derived from the ITS2 data largely agreed with the phylogenetic hypotheses from morphologic and other molecular evidence. Our analyses suggest that the accordant resolutions generated from different analyses can be used to infer deep-level phylogenetic relations.     

Phylogenetic hypotheses of gorgoniid octocorals according to ITS2 and their predicted RNA secondary structures

Gorgoniid octocorals taxonomy (Cnidaria; Octocorallia; Gorgoniidae) includes diagnostic characters not well defined at the generic level, and based on the family diagnosis some species could be classified in either Gorgoniidae or Plexauridae. In this study, we used sequences from the Internal Transcribed Spacer 2 (ITS2) and their predicted RNA secondary structure to both correct the alignment and reconstruct phylogenies using molecular morphometrics for 24 octocorals mostly from the Atlantic. ITS2 exhibited the six-helicoidal ring-model structure found in eukaryotes, and provided 38 parsimony-informative characters. The proposed phylogenies, though differing between sequence- and structure-base results, provided consistent support for several clades. Genera considered part of the polyphyletic genus Leptogorgia, such as Filigorgia, were distantly related to the former in all phylogenetic hypotheses. Main differences among the hypotheses consisted in the placement of Muriceopsis (previously considered from the Plexauridae family) and Filigorgia. Excluding Muriceopsis and an undescribed octocoral from Tobago, Plexaurella and Pterogorgia grouped together as a sister branch of Pinnigorgia spp. but long-branch attraction was evident for the grouping of Plexaurella nutans (another plexaurid) and Pterogorgia citrina. Unexpected results were the divergence between Caribbean genera, Gorgonia and Pseudopterogorgia, which were placed basal respect to Pacifigorgia and Leptogorgia (=Lophogorgia). ITS2 provided support to corroborate observations based on sclerite morphology: species with "capstan sclerites" (e.g., Pacifigorgia and Leptogorgia) were characterized by a long helix IV with one internal loop and a helix V with four internal loops; "scaphoid sclerites" had a predominantly long helix V if compared to helix IV; "asymmetric spiny sclerites" (Muriceopsis, Pinnigorgia and the undescribed octocoral) exhibited one or two lateral bulges in the V helix. Remarkably, Muriceopsis and Pinnigorgia were supported by a complete Compensatory Base Change (CBC) (A-U to G-C) in helix V. Filigorgia with simple "spindles" had a short helix IV and a large central ring. DNA sequences from the nuclear ITS2 region, including information from predicted RNA secondary structure, despite their reduced length, provided numerous characters and phylogenetic information among Gorgoniidae genera and species.     

Phylogenetic relationships among Myrceugenia, Blepharocalyx, and Luma (Myrtaceae) based on paired-sites models and the secondary structures of ITS and ETS sequences

To establish relationships among the genera Blepharocalyx, Luma, and Myrceugenia, the phylogeny of tribe Myrteae was reconstructed based on secondary structures of the sequences of ITS (ITS1-5.8S-ITS2) and ETS regions from 93 taxa belonging to 29 genera. The DNA sequences were aligned according to their secondary structure and analysed with Bayesian inference (BI) using base substitution models for RNA, which can differentiate between the regions consisting of base pairs (helices) and unpaired bases (loops). The best-fit models were RNA7C for ITS and RNA7B for ETS, which differ in how they predict double substitutions and symmetry of the frequency of base pairs. The phylogenetic hypothesis was compared with results obtained using classical 4-state models, the results of maximum parsimony, and using sequence alignment without considering the secondary structure. Analyses show low support along the backbone of the consensus trees, those using substitution models of paired sites showing more highly supported derived clades than substitution models of independent sites. All tests were consistent and differed in the positioning of certain groups, but they also showed that Blepharocalyx, Luma, and Myrceugenia arise from three independent lineages that are not closely related. The substitution model for RNA shows that Luma is a monophyletic group and sister to the other genera of Myrteae, except for Myrtus communis. Myrceugenia appears as a monophyletic group, except for M. fernandeziana, which is related to Blepharocalyx, a genus that appears to be biphyletic.     

Phylogenetic reconstruction of endophytic fungal isolates using internal transcribed spacer 2 (ITS2) region

Endophytic fungi are inhabitants of plants, living most part of their lifecycle asymptomatically which mainly confer protection and ecological advantages to the host plant. In this present study, 48 endophytic fungi were isolated from the leaves of three medicinal plants and characterized based on ITS2 sequence – secondary structure analysis. ITS2 secondary structures were elucidated with minimum free energy method (MFOLD version 3.1) and consensus structure of each genus was generated by 4SALE. ProfDistS was used to generate ITS2 sequence structure based phylogenetic tree respectively. Our elucidated isolates were belonging to Ascomycetes family, representing 5 orders and 6 genera. Colletotrichum/Glomerella spp., Diaporthae/Phomopsis spp., and Alternaria spp., were predominantly observed while Cochliobolus sp., Cladosporium sp., and Emericella sp., were represented by singletons. The constructed phylogenetic tree has well resolved monophyletic groups with >50% bootstrap value support. Secondary structures based fungal systematics improves not only the stability; it also increases the precision of phylogenetic inference. Above ITS2 based phylogenetic analysis was performed for our 48 isolates along with sequences of known ex-types taken from GenBank which confirms the efficiency of the proposed method. Further, we propose it as superlative marker for reconstructing phylogenetic relationships at different taxonomic levels due to their lesser length.     

Phylogenetics of Anthyllis (Leguminosae: Papilionoideae: Loteae): Partial incongruence between nuclear and plastid markers, a long branch problem and implications for morphological evolution

Phylogenetic relationships in the genus Anthyllis (Leguminosae: Papilionoideae: Loteae) were investigated using data from the nuclear ribosomal internal transcribed spacer regions (ITS) and three plastid regions (psbA-trnH intergenic spacer, petB-petD region and rps16 intron). Bayesian and maximum parsimony (MP) analysis of a concatenated plastid dataset recovered well-resolved trees that are topologically similar, with many clades supported by unique indels. MP and Bayesian analyses of the ITS sequence data recovered trees that have several well-supported topological differences, both among analyses, and to trees inferred from the plastid data. The most substantial of these concerns A. vulneraria and A. lemanniana, whose placement in the parsimony analysis of the ITS data appears to be due to a strong long-branch effect. Analysis of the secondary structure of the ITS1 spacer showed a strong bias towards transitions in A. vulneraria and A. lemanniana, many of which were also characteristic of certain outgroup taxa. This may contribute to the conflicting placement of this clade in the MP tree for the ITS data. Additional conflicts between the plastid and ITS trees were more taxonomically focused. These differences may reflect the occurrence of reticulate evolution between closely related species, including a possible hybrid origin for A. hystrix. The patterns of incongruence between the plastid and the ITS data seem to correlate with taxon ranks. All of our phylogenetic analyses supported the monophyly of Anthyllis (incl. Hymenocarpos). Although they are often taxonomically associated with Anthyllis, the genera Dorycnopsis and Tripodion are shown here to be more closely related to other genera of Loteae. We infer up to six major clades in Anthyllis that are morphologically well-characterized, and which could be recognized as sections. Four of these agree with various morphology-based classifications, while the other two are novel. We reconstruct the evolution of several morphological characteristics found only in Anthyllis or tribe Loteae. Some of these characters support major clades, while others show evidence of homoplasy within Anthyllis.     

ITS secondary structure derived from comparative analysis: implications for sequence alignment and phylogeny of the Asteraceae

An RNA secondary structure model is presented for the nuclear ribosomal internal transcribed spacers (ITS) based on comparative analysis of 340 sequences from the angiosperm family Asteraceae. The model based on covariation analysis agrees with structural features proposed in previous studies using mainly thermodynamic criteria and provides evidence for additional structural motifs within ITS1 and ITS2. The minimum structure model suggests that at least 20% of ITS1 and 38% of ITS2 nucleotide positions are involved in base pairing to form helices. The sequence alignment enabled by conserved structural features provides a framework for broadscale molecular evolutionary studies and the first family-level phylogeny of the Asteraceae based on nuclear DNA data. The phylogeny based on ITS sequence data is very well resolved and shows considerable congruence with relationships among major lineages of the family suggested by chloroplast DNA studies, including a monophyletic subfamily Asteroideae and a paraphyletic subfamily Cichorioideae. Combined analyses of ndhF and ITS sequences provide additional resolution and support for relationships in the family.     

Genetic diversity of nuclear ITS1-5.8S-ITS2 rDNA sequence in Clonorchis sinensis Cobbold, 1875 (Trematoda: Opisthorchidae) from the Russian Far East

The present study examined the molecular organisation and sequence variation in the nuclear ribosomal DNA (rDNA) region, including the two internal transcribed spacers (ITS1 and ITS2) and the 5.8S gene of the Clonorchis sinensis from the Russian Far East. The relevant sequences from other parts of this species' area were downloaded from GenBank. The results showed 100% identity for all investigated 5.8S-ITS2 rDNA sequences. In contrast, two levels of intraspecific variations were revealed in the complete ITS1 sequences. The intra-genomic variation resulted from a C/T polymorphism in a single position. The inter-individual differences between the ITS1 sequences were both due to nucleotide and size polymorphisms resulting from a varying number of five-nucleotide repeats and followed by two ITS1 length variants. These variant frequencies correlate with the clonorchiasis level in some geographical localities. ITS1 differences, both in the mutation profile and mutation localisation, were revealed between northern and southern geographical samples. The presence of GC boxes that are identical to known regulatory motifs in eukaryotes was detected within the ITS1 sub-repeats. The predicted secondary structures for ITS1 consist of two large branches, one of which was invariable, while another depended on ITS1 length. The predicted secondary structure for ITS2 includes four helices around the core. The main differences between C. sinensis and other opisthorchids were localised on the tops of helices 2, 3, and 4. A phylogenetic MST reconstruction subdivided all ITS1 sequences into two well differentiated clusters, each with the major widespread ribotype, and showed that ribotype diversity in both Russia and Korea is much lower than in China. The results obtained demonstrate the feasibility of complete ITS1 sequences in C. sinensis population genetics and can be considered as a basis for further studies of the parasite infection because they may help to elucidate the molecular mechanisms of pathogen evolution and adaptation.     

Evolution of Helix Formation in the Ribosomal Internal Transcribed Spacer 2 (ITS2) and Its Significance for RNA Secondary Structures

Helices are the most common elements of RNA secondary structure. Despite intensive investigations of various types of RNAs, the evolutionary history of the formation of new helices (novel helical structures) remains largely elusive. Here, by studying the nuclear ribosomal Internal Transcribed Spacer 2 (ITS2), a fast-evolving part of the eukaryotic nuclear ribosomal operon, we identify two possible types of helix formation: one type is “dichotomous helix formation”—transition from one large helix to two smaller helices by invagination of the apical part of a helix, which significantly changes the shape of the original secondary structure but does not increase its complexity (i.e., the total length of the RNA). An alternative type is “lateral helix formation”—origin of an extra helical region by the extension of a bulge loop or a spacer in a multi-helix loop of the original helix, which does not disrupt the pre-existing structure but increases RNA size. Moreover, we present examples from the RNA sequence literature indicating that both types of helix formation may have implications for RNA evolution beyond ITS2.     

Genetic heterogeneity in internal transcribed spacer genes of Balantidium coli (Litostomatea, Ciliophora)

The species Balantidium coli is the only ciliate that parasitizes humans. It has been described in other primates, and it has been proposed that the species B. suis from pigs and B. struthionis from ostriches are synonyms of B. coli. Previous genetic analysis of pig and ostrich Balantidium isolates found a genetic polymorphism in the ITS region but its taxonomic relevance was not established. We have extended the genetic analysis to Balantidium isolates of pig, gorilla, human and ostrich origin. We have PCR-amplified and sequenced the ITS region of individual Balantidium cells. The predicted ITS secondary structures of the sequences obtained were transferred by homology modelling to the sequences of other Trichostomatia ciliates (Isotricha, Troglodytella, Lacrymaria and Spathidium) and compared to determine the importance of the differences in the primary sequences. The results show that the ITS2 secondary structure of the species considered follows the general pattern of other ciliates, although with some deviations. There are at least two main types of ITS sequence variants in B. coli which could be present in the same cell and they are common to the mammal and avian hosts studied. These data do not support B. suis and B. struthionis as distinct species.     

Conservation of sequence in the internal transcribed spacers and 5.8S ribosomal RNA among geographically separated isolates of parasitic scuticociliates (Ciliophora, Orchitophryidae)

Nucleotide sequence from the internal transcribed spacers (ITS1 and ITS2) and the 5.8S gene from the ribosomal RNA gene cluster of isolates of the scuticociliate Orchitophrya stellarum from 4 asteroid hosts were compared. Surprisingly, these data (495 bp) were identical for O. stellarum isolated from the testes of Asterias amurensis from Japan; Pisaster ochraceus from British Columbia, Canada; Asterias rubens from The Netherlands; and Asterias vulgaris from Prince Edward Island, Canada. These sequence data were compared to those from 3 scuticociliates which parasitise crustaceans: Mesanophrys pugettensis, M. chesapeakensis and Anophryoides haemophila. No difference was found in this region between the nucleotide sequence of M. pugettensis and M. chesapeakensis. The sequence of Mesanophrys spp. differed by 9.2% in the ITS1 and 4.7% in the ITS2 from that of O. stellarum. The sequence from the ITS1 (135 bp) and ITS2 (233 bp) of A. haemophila differed by 42.6 and 20.5% respectively from those of O. stellarum. Therefore, nucleotide sequence of the ITS regions in these scuticociliates is highly conserved.     

Ovalopodium desertum n. sp. and the phylogenetic relationships of Cochliopodiidae (Amoebozoa)

An amoeba isolated from a weakly saline semi-desert pond in Kazakhstan (Central Asia) resembles a small Cochliopodium in the light microscope, but has a dorsal fibrous cell coat without scales. Thus it can be identified morphologically as a new species of Ovalopodium Sawyer, 1980, and it is herein named O. desertum. Phylogenetic analysis of the SSU rRNA gene sequences of the new species and four Cochliopodium spp. sequenced additionally shows that Ovalopodium desertum is a sister clade to a robustly monophyletic Cochliopodium. The close relationship between Ovalopodium and Cochliopodium is also confirmed by the analysis of SSU rRNA secondary structure showing the specific helices in the region V5 in all species of both genera. Analysis of actin gene sequences fails to resolve the position of Ovalopodium but demonstrates that Parvamoeba Rogerson, 1993 is probably related to Cochliopodium. The position of Cochliopodiidae within Amoebozoa remains unresolved, despite our efforts to resolve it using broader taxonomic sampling of Amoebozoa, testing alternative tree topologies and removing the fast-evolving sites. Among sequenced genera, Parvamoeba and Endostelium Olive et al., 1984 are probable relatives to Cochliopodiidae. Molecular trees weakly support an inclusion of the family in Flabellinia (Discosea), but more phylogenomic data are necessary to test this hypothesis.     

Secondary Structure of the ITS1 Transcript and its Application in a Reconstruction of the Phylogeny of Boraginales

Abstract: In this study, we present the possibilities for calculating systematic relationships on higher taxonomical levels based exclusively on ITS1. This is demonstrated for Boraginales (Boraginaceae s.str., Cordiaceae, Ehretiaceae, Heliotropiaceae, Hydrophyllaceae and Lennoaceae). Secondary structure of the ITS1 region is more conserved than the primary structure (i.e., sequence itself) and is therefore a tool for optimizing alignments. It increases the number of structural characters. Information inferred from the secondary structure enables us to construct well-resolved phylogenetic trees at higher taxonomic levels. A general secondary structure of ITS1 for Boraginales, with four major helices, is proposed here. In each subtaxon, derivations from this structure are found. The paraphyly of Boraginaceae s.l. is evident from both a comparison of secondary structures and from bootstrap analysis. Boraginaceae s.str. are the sister group of a clade formed by Hydrophyllaceae, Cordiaceae, Ehretiaceae, Heliotropiaceae and Lennoaceae. The last four taxa constitute a monophyletic group which is the sister group of Hydrophyllaceae. Lennoaceae are closely related to Ehretiaceae, and these two taxa in turn are the sister group of Cordiaceae.     

Phylogenetic analysis and predicted secondary structures of the rDNA internal transcribed spacers of the powdery mildew fungi (Erysiphaceae)

The nucleotide sequences of the internal transcribed spacer (ITS) regions of the ribosomal DNA including the 5.8S rRNA gene and the 5′ end of the 28S rRNA gene have been determined for 19 species in 10 genera of the powdery mildew fungi in order to analyze their phylogenetic relationship. These fungi were divided into two large groups based on the nucleotide length of the ITS regions, and this grouping was in line with that based on the morphological characters of the anamorphic stage rather than the teleomorphic stage. Although the variable ITS sequences were often ambiguously aligned, conserved sites were also found. Thus, a neighbor-joining tree was constructed using the nucleotide sequence data of the conserved sites of the ITS regions, the 5.8S rRNA gene, and the 5′ end of the 28S rRNA gene. The phylogenetic tree displayed the presence of four groups in the powdery mildews, which were distinguished by their morphology and/or host ranges. In the ITS2 region, the presence of a common secondary structure having four hairpin domains was suggested, in spite of the highly variable nucleotide sequences of this region. The predicted secondary structure was supported by the compensatory mutations as well as compensatory conserved sequences and high G+C content in the predicted stem regions. Contribution No. 142 from the Laboratory of Plant Pathology, Mie University.     

Evolution of secondary heads in Nassauviinae (Asteraceae, Mutisieae)

The evolution of the inflorescence head in Asteraceae is important in the diversification of this largest angiosperm family. The aggregation of heads into higher-order capitulescences (secondary heads or syncephalia) is considered evolutionarily advanced. The genera Moscharia, Nassauvia, Polyachyrus, and Triptilion of the subtribe Nassauviinae (Mutisieae) have syncephalia with differing degrees of capitula condensation. ITS and plastid trnL-trnF regions were analyzed separately and together using maximum parsimony and maximum likelihood to examine the evolution of syncephalia in the Nassauviinae. The four genera displaying syncephalia do not form a clade minus taxa without syncephalia, indicating that secondary heads in Nassauviinae have either convergently evolved twice in the subtribe (or, very unlikely) once with multiple reversions. Strong support was obtained for a sister relationship between Leucheria (without syncephalium) and Polyachyrus, and both sister to Moscharia. Nassauvia and Triptilion form a distinct clade but are sister to other genera, Perezia and Panphalea, without syncephalium. Previous hypotheses postulated the evolution from simple to more complex secondary heads. We show that the ancestor of Moscharia, Polyachyrus, and Leucheria, in a more arid habitat, had a complex type of secondary head, and loss of complexity occurred in response to a shift from arid to mesic conditions.     

Internal transcribed spacer 2 (nu ITS2 rRNA) sequence-structure phylogenetics: towards an automated reconstruction of the green algal tree of life

Background

Chloroplast-encoded genes (matK and rbcL) have been formally proposed for use in DNA barcoding efforts targeting embryophytes. Extending such a protocol to chlorophytan green algae, though, is fraught with problems including non homology (matK) and heterogeneity that prevents the creation of a universal PCR toolkit (rbcL). Some have advocated the use of the nuclear-encoded, internal transcribed spacer two (ITS2) as an alternative to the traditional chloroplast markers. However, the ITS2 is broadly perceived to be insufficiently conserved or to be confounded by introgression or biparental inheritance patterns, precluding its broad use in phylogenetic reconstruction or as a DNA barcode. A growing body of evidence has shown that simultaneous analysis of nucleotide data with secondary structure information can overcome at least some of the limitations of ITS2. The goal of this investigation was to assess the feasibility of an automated, sequence-structure approach for analysis of IT2 data from a large sampling of phylum Chlorophyta.

Methodology/Principal Findings

Sequences and secondary structures from 591 chlorophycean, 741 trebouxiophycean and 938 ulvophycean algae, all obtained from the ITS2 Database, were aligned using a sequence structure-specific scoring matrix. Phylogenetic relationships were reconstructed by Profile Neighbor-Joining coupled with a sequence structure-specific, general time reversible substitution model. Results from analyses of the ITS2 data were robust at multiple nodes and showed considerable congruence with results from published phylogenetic analyses.

Conclusions/Significance

Our observations on the power of automated, sequence-structure analyses of ITS2 to reconstruct phylum-level phylogenies of the green algae validate this approach to assessing diversity for large sets of chlorophytan taxa. Moreover, our results indicate that objections to the use of ITS2 for DNA barcoding should be weighed against the utility of an automated, data analysis approach with demonstrated power to reconstruct evolutionary patterns for highly divergent lineages.     

ITS-2 secondary structures and phylogeny of Anopheles culicifacies species

Background: Second internal transcribed spacer (ITS2) has proven to contain useful biological information at higher taxonomic levels. Objectives: This study was carried out to unravel the biological information in the ITS2 region of An. culicifacies and the internal relationships between the five species of Anopheles culicifacies. Methodology: In achieving these objectives, twenty two ITS2 sequences (~370bp) of An. culicifacies species were retrieved from GenBank and secondary structures were generated. For the refinement of the primary structures, i.e. nucleotide sequence of ITS2 sequences, generated secondary structures were used. The improved ITS2 primary structures sequences were then aligned and used for the construction of phylogenetic trees. Results and discussions: ITS2 secondary structures of culicifacies closely resembled near universal eukaryotes secondary structure and had three helices, and the structures of helix II and distal region of helix III of ITS2 of An. culicifacies were strikingly similar to those regions of other organisms strengthening possible involvement of these regions in rRNA biogenesis. Phylogenetic analysis of improved ITS2 sequences revealed two main clades one representing sibling B, C and E and A and D in the other. Conclusions: Near sequence identity of ITS2 regions of the members in a particular clade indicate that this region is undergoing parallel evolution to perform clade specific RNA biogenesis. The divergence of certain isolates of An. culicifacies from main clades in phylogenetic analyses suggests the possible existence of camouflaged sub-species within the complex of culicifacies. Using the fixed nucleotide differences, we estimate that these two clades have diverged nearly 3.3 million years ago, while the sibling species in clade 2 are under less evolutionary pressure, which may have evolved much later than the members in clade 1.