Cirripede phylogeny using a novel approach: molecular morphometrics

Detection and Molecular Phylogenetic-Morphometric Characterization of Rhizoctonia tuliparum, Causal Agent of Gray Bulb Rot of Tulips and Bulbous Iris

Gray bulb rot of tulips and bulbous iris is caused by the soil-borne fungal pathogen, Rhizoctonia tuliparum (Rtul). Sclerotia present in infected bulbs, as well as overwintering sclerotia in soil and field debris, are the primary sources of infection. A method for accurate and sensitive detection of Rtul from soil and infected bulbs, and estimation of inoculum threshold levels, is needed for the management of disease caused by this pathogen. We designed a unique set of primers targeting the ITS2 region of the Rtul genome and developed a highly sensitive quantitative PCR (qPCR)-based method for Rtul identification using these primers, where the threshold of detection was approximately 1 fg Rtul DNA. The assay was more sensitive with sclerotia collected from the field (natural) than with those grown in the lab, and more sensitive with natural-light than natural-dark sclerotia. Also, the detection method was more sensitive when sclerotia were extracted from soil than from bulb tissue. The qPCR method was highly specific, as no PCR amplification was detected when genomic DNA from 62 non-Rtul Rhizoctonia isolates from a wide range of anastomosis groups were tested. To understand the evolutionary relationships and genomic diversity of Rtul, we performed phylogenetics of the ITS1-5.8S-ITS2 region and ITS2-molecular morphometric characterization (MMC) of Rtul isolates. The three Rtul isolates whose ITS sequences were available in GenBank formed a distinct phylogenetic clade with Ceratobasidium anceps as the nearest relative. Furthermore, MMC analysis revealed genetic divergence among these three Rtul isolates.

Morphopathological and Molecular Morphometric Characterization of Waitea circinata var. prodigus Causing a Novel Sheath Spot Disease of Maize in India

Maize brown sheath spot (MBSS), a new disease of maize, was discovered while surveying for maize leaf and sheath blight diseases in the Indian states of Assam, Jharkhand, Meghalaya, Manipur, and Odisha. Maize is the third most important cereal after rice and wheat in India. Unlike banded leaf and sheath blight disease caused by Rhizoctonia solani, MBSS symptoms on maize were discrete and limited to sheaths only. Symptoms of MBSS in the field were initially water-soaked necrotic lesions of 1 to 2 cm in diameter on the lowermost leaf sheaths, which then progressed to the upper sheaths. Lesions coalesced and covered approximately 2 to 5% of the sheath area. Infected dried lower leaves were shed, whereas infected upper leaves remained on the stem. The pathogen was isolated, characterized morphologically, pathologically, and molecularly, and identified as Waitea circinata var. prodigus, a basidiomycete known to cause basal leaf blight of seashore paspalum. The internal transcribed spacer (ITS) sequence 2 (ITS2) of rDNA from MBSS isolates formed a well supported clade with known W. circinata var. prodigus isolates. Molecular morphometric analysis of the ITS2 regions of the five known varieties of W. circinata detected distinguishing variations in GC content, compensatory base changes (CBCs), hemi- CBCs, indels, and altered base-pairing of helices. Variation in these characteristics may indicate that varieties are distinct biological species within W. circinata sensu lato. The geographical distribution and potential impacts of MBSS on the maize crop in India necessitate further investigations of pathogen identification and disease management.

Molecular identification and phylogenetic analysis of Papaver based on ITS2 barcoding

In forensic cases suspected to involve Papaver somniferum, species identification is key to the investigation. To accurately detect and identify P. somniferum as well as common adulterants of the same genus, 19 internal transcribed spacer 2 (ITS2) sequences of P. somniferum (256 bp), Papaver canescens (254 bp), Papaver nudicaule (254 bp), Papaver pavoninum (250 bp), Papaver radicatum (254 bp), and Papaver rhoeas (256 bp) were obtained. Based on the ITS2 sequence, similarity analysis via BLAST, the nearest Kimura-2-parameter (K2P) genetic distances were calculated, and a phylogenetic tree was constructed using MEGA X software for the identification of six species of Papaver. Finally, differences in the ITS2 secondary structure between species were analyzed. The best matches of the P. somniferum ITS2 sequence were of other P. somniferum from different sources. The nearest K2P genetic distances between P. somniferum and its counterparts from other sources were zero, which was the smallest pairwise genetic distance among distances from the other five Papaver species. Various sources of P. somniferum clustered into an independent branch in the phylogenetic tree. The secondary structures of P. somniferum and P. rhoeas were significantly different from those of the other four species of Papaver. In summary, P. somniferum can be effectively distinguished from five closely related plants of the same genus by using ITS2 as a DNA barcode.

Menzerath-Altmann's Law of Syntax in RNA Accretion History

RNA evolves by adding substructural parts to growing molecules. Molecular accretion history can be dissected with phylogenetic methods that exploit structural and functional evidence. Here, we explore the statistical behaviors of lengths of double-stranded and single-stranded segments of growing tRNA, 5S rRNA, RNase P RNA, and rRNA molecules. The reconstruction of character state changes along branches of phylogenetic trees of molecules and trees of substructures revealed strong pushes towards an economy of scale. In addition, statistically significant negative correlations and strong associations between the average lengths of helical double-stranded stems and their time of origin (age) were identified with the Pearson's correlation and Spearman's rho methods. The ages of substructures were derived directly from published rooted trees of substructures. A similar negative correlation was detected in unpaired segments of rRNA but not for the other molecules studied. These results suggest a principle of diminishing returns in RNA accretion history. We show this principle follows a tendency of substructural parts to decrease their size when molecular systems enlarge that follows the Menzerath-Altmann's law of language in full generality and without interference from the details of molecular growth.

Genetic diversity of Paramecium species on the basis of multiple loci analysis and ITS secondary structure models

Among ciliates, Paramecium has become a privileged model for the study of "species problem" particularly in the case of the "Paramecium aurelia complex" that has been intensely investigated. Despite extensive studies, the taxonomy of Paramecium is still challenging. The major problem is an uneven sampling of Paramecium with relatively few representatives of each species. To investigate species from the less discovered region (Pakistan), 10 isolates of Paramecium species including a standing-alone FT8 strain previously isolated by some of us were subjected to molecular characterization. Fragments of 18S recombinant DNA (rDNA), ITS1-5.8S-ITS2-5'LSU rDNA, cytochrome c oxidase subunit II, and hsp70 genes were used as molecular markers for phylogenetic analysis of particular isolates. The nucleotide sequences of polymerase chain reaction products of all markers were compared with the available sequences of relevant markers of other Paramecium species from GenBank. Phylogenetic trees based on all molecular markers showed that all the nine strains had a very close relationship with Paramecium primaurelia except for the FT8 strain. FT8 consistently showed its unique position in comparison to all other species in the phylogenetic trees. Available sequences of internal transcribed spacer 1 (ITS1) and ITS2 and some other ciliate sequences from GenBank were used for the construction of secondary models. Two highly conserved helices supported by compensatory base changes among all ciliates of ITS2 secondary structures were found similar to other eukaryotes. Therefore, the most conserved 120 to 180 base pairs regions were identified for their comparative studies. We found that out of the three helices in ITS1 structure, helix B was more conserved in Paramecium species. Despite various substitutions in the primary sequence, it was observed that secondary structures of ITS1 and ITS2 could be helpful in interpreting the phylogenetic relationships both at species as well as at generic level.

Evaluation of ITS2 molecular morphometrics effectiveness in species delimitation of Ascomycota - A pilot study

Exploring the secondary structure information of nuclear ribosomal internal transcribed spacer 2 (ITS2) has been a promising approach in species delimitation. However, Compensatory base changes (CBC) concept employed in this approach turns futile when CBC is absent. This prompted us to investigate the utility of insertion/deletion (INDELs) and substitutions in fungal delineation at species level. Upon this rationale, 116 strains representing 97 species, belonging to 6 genera (Colletotrichum, Boeremia, Leptosphaeria, Peyronellaea, Plenodomus and Stagonosporopsis) of Ascomycota were retrieved from Q-bank for molecular morphometric analysis. CBC, INDELs and substitutions between the species of their respective genus were recorded. Most species combinations lacked CBC. Among the substitution events, transitions were predominant. INDELs were less frequent than the substitutions. These evolutionary events were mapped upon the helices to discern species specific variation sites. In 68 species unique variation sites were recognised. The remaining 29 species shared absolute similarity with distinctly named species. The variation sites catalogued in them overlapped with other distinct species and resulted in the blurring of species boundaries. Species specific variation sites recognized in this study are the preliminary results and they could be discerned with absolute confidence when larger datasets encompassing all described species of genera were investigated. They could be of potential use in barcoding fungi at species level. This study also concludes that the ITS2 molecular morphometric analysis is an efficient third dimensional study of the fungal species delimitation. This may help to avoid the false positives in species delimitations and to alleviate the challenges in molecular characterization.

Internal transcribed spacer 2 (ITS2) molecular morphometric analysis based species delimitation of foliar endophytic fungi from Aglaia elaeagnoidea, Flacourtia inermis and Premna serratifolia

Molecular morphometrics is an emerging third dimensional aspect of fungal species delimitation. They have been demonstrated to be more informative than conventional barcoding methods. Hence in this study, foliar endophytic fungal (FEF) assemblages in three Magnoliopsida plants were delimited using nuclear ribosomal internal transcribed spacer 2 (ITS2) sequence-secondary structural features based phylogenetic analysis, also known as molecular morphometrics. A total of 392 FEF isolates were obtained from the Aglaia elaeagnoidea, Flacourtia inermis, and Premna serratifolia leaves and grouped into 98 morphotypes. Among these host plants, P. serratifolia showed the maximum percentage of colonization frequency. Representatives of each morphotype was sequenced and subjected to further molecular characterization. The results revealed that morphotypes were belonged to the phylum of Ascomycota, distributed over two classes (Sordariomycetes (68.59%) and Dothideomycetes (31.41%)), 6 orders and 19 genera. Based on compensatory base changes (CBC) analysis and absolute identity of ITS2 structure, 21, 20 and 23 species were recognized from A. elaeagnoidea, F. inermis, and P. serratifolia respectively. Diversity indices were higher in A. elaeagnoidea, despite it accounted for a modest 16.8% of total isolates recorded in this study. The genus Colletotrichum was predominant in A. elaeagnoidea (39%) and P. serratifolia (48%). Similarly, Diaporthe (43%) was dominant in F. inermis. Several host-specific species were also observed. This study concludes that these plants host diverse species of Ascomycota. To the best of our knowledge, this is the first detailed report on FEF diversity from these plants. Also, the inclusion of ITS2 secondary structure information along with the sequence provides a further dimension to resolve the inherent problems in identification of fungal species.

Computing the origin and evolution of the ribosome from its structure - Uncovering processes of macromolecular accretion benefiting synthetic biology

Accretion occurs pervasively in nature at widely different timeframes. The process also manifests in the evolution of macromolecules. Here we review recent computational and structural biology studies of evolutionary accretion that make use of the ideographic (historical, retrodictive) and nomothetic (universal, predictive) scientific frameworks. Computational studies uncover explicit timelines of accretion of structural parts in molecular repertoires and molecules. Phylogenetic trees of protein structural domains and proteomes and their molecular functions were built from a genomic census of millions of encoded proteins and associated terminal Gene Ontology terms. Trees reveal a ‘metabolic-first’ origin of proteins, the late development of translation, and a patchwork distribution of proteins in biological networks mediated by molecular recruitment. Similarly, the natural history of ancient RNA molecules inferred from trees of molecular substructures built from a census of molecular features shows patchwork-like accretion patterns. Ideographic analyses of ribosomal history uncover the early appearance of structures supporting mRNA decoding and tRNA translocation, the coevolution of ribosomal proteins and RNA, and a first evolutionary transition that brings ribosomal subunits together into a processive protein biosynthetic complex. Nomothetic structural biology studies of tertiary interactions and ancient insertions in rRNA complement these findings, once concentric layering assumptions are removed. Patterns of coaxial helical stacking reveal a frustrated dynamics of outward and inward ribosomal growth possibly mediated by structural grafting. The early rise of the ribosomal ‘turnstile’ suggests an evolutionary transition in natural biological computation. Results make explicit the need to understand processes of molecular growth and information transfer of macromolecules.

ITS2 Secondary Structure Improves Discrimination between Medicinal "Mu Tong" Species when Using DNA Barcoding

DNA barcoding is a promising species identification method, but it has proved difficult to find a standardized DNA marker in plant. Although the ITS/ITS2 RNA transcript has been proposed as the core barcode for seed plants, it has been criticized for being too conserved in some species to provide enough information or too variable in some species to align it within the different taxa ranks. We selected 30 individuals, representing 16 species and four families, to explore whether ITS2 can successfully resolve species in terms of secondary structure. Secondary structure was predicted using Mfold software and sequence-structure was aligned by MARNA. RNAstat software transformed the secondary structures into 28 symbol code data for maximum parsimony (MP) analysis. The results showed that the ITS2 structures in our samples had a common four-helix folding type with some shared motifs. This conserved structure facilitated the alignment of ambiguous sequences from divergent families. The structure alignment yielded a MP tree, in which most topological relationships were congruent with the tree constructed using nucleotide sequence data. When the data was combined, we obtained a well-resolved and highly supported phylogeny, in which individuals of a same species were clustered together into a monophyletic group. As a result, the different species that are often referred to as the herb “Mu tong” were successfully identified using short fragments of 250 bp ITS2 sequences, together with their secondary structure. Thus our analysis strengthens the potential of ITS2 as a promising DNA barcode because it incorporates valuable secondary structure information that will help improve discrimination between species.

The complete mitogenome of Apocheima cinerarius (Lepidoptera: Geometridae: Ennominae) and comparison with that of other lepidopteran insects

The complete mitochondrial genome (mitogenome) of a female flightless geometrid moth Apocheima cinerarius was found to be 15,722 bp in length, containing 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a control region. The A+T content of the complete mitogenome is 80.83%. The AT skew value ([A-T]/[A+T]) is 0.027. The 13 PCGs of the mitogenome start with typical ATN codons, except for cox1 with the start codon CGA. All the tRNA genes have typical cloverleaf secondary structures, except for trnSer(AGN). The secondary structures of rrnL and rrnS were predicted. Six structural domains including conserved regions (IV, V) and variable regions (I, II, III, VI) were identified in the secondary structure of rrnL. The secondary structure of rrnS consists of 3 structural domains. The control region of A. cinerarius begins with conserved motifs of "ATAGA"+19-bp poly T. It also contains a microsatellite-like (TA)26, a stem-and-loop structure, and a poly-A stretch. Phylogenetic analysis showed that Geometroidea is more closely related to Bombycoidea than to Noctuoidea. A. cinerarius is more closely related to Biston panterinaria than to Phthonandria atrilineata, which is in accordance with the conventional morphology-based classification.

Archaea: the first domain of diversified life

The study of the origin of diversified life has been plagued by technical and conceptual difficulties, controversy, and apriorism. It is now popularly accepted that the universal tree of life is rooted in the akaryotes and that Archaea and Eukarya are sister groups to each other. However, evolutionary studies have overwhelmingly focused on nucleic acid and protein sequences, which partially fulfill only two of the three main steps of phylogenetic analysis, formulation of realistic evolutionary models, and optimization of tree reconstruction. In the absence of character polarization, that is, the ability to identify ancestral and derived character states, any statement about the rooting of the tree of life should be considered suspect. Here we show that macromolecular structure and a new phylogenetic framework of analysis that focuses on the parts of biological systems instead of the whole provide both deep and reliable phylogenetic signal and enable us to put forth hypotheses of origin. We review over a decade of phylogenomic studies, which mine information in a genomic census of millions of encoded proteins and RNAs. We show how the use of process models of molecular accumulation that comply with Weston's generality criterion supports a consistent phylogenomic scenario in which the origin of diversified life can be traced back to the early history of Archaea.

Secondary structure models of 18S and 28S rRNAs of the true bugs based on complete rDNA sequences of Eurydema maracandica Oshanin, 1871 (Heteroptera, Pentatomidae)

The sequences of 18S and 28S rDNAs have been used as molecular markers to resolve phylogenetic relationships of Heteroptera for two decades. The complete sequences of 18S rDNAs have been used in many studies, while in most studies only partial sequences of 28S rDNAs have been used due to technical difficulties of amplifying the complete lengths. In this study, we amplified the complete 18S and 28S rDNA sequences of Eurydema maracandica Oshanin, 1871, and reconstructed the secondary structure models of the corresponding rRNAs. In addition, and more importantly, all of the length variable regions of 18S rRNA were compared among 37 families of Heteroptera based on 140 sequences, and the D3 region of 28S rRNA was compared among 51 families based on 84 sequences. It was found that 8 length variable regions could potentially serve as molecular synapomorphies for some monophyletic groups. Therefore discoveries of more molecular synapomorphies for specific clades can be anticipated from amplification of complete 18S and 28S rDNAs of more representatives of Heteroptera.

Secondary structure and phylogenetic utility of the ribosomal large subunit (28S) in monogeneans of the genus Thaparocleidus and Bifurcohaptor (Monogenea: Dactylogyridae)

Present communication deals with secondary structure of 28S rDNA of two already known species of monogeneans viz., Bifurcohaptor indicus and Thaparocleidus parvulus parasitizing gill filaments of a freshwater fish, Mystus vittatus for phylogenetic inference. Secondary structure data are best used as accessory taxonomic characters as their phylogenetic resolving power and confidence in validity. Secondary structure of the 28S rDNA transcript could provide information for identifying homologous nucleotide characters, useful for cladistic inference of relationships. Such structure data could be used as taxonomic character. The study supports that species-level sequence variability renders 28S sequence as a unique window for examining the behavior of fast evolving, non-coding DNA sequences. Apart from this it also confirms that molecular similarity present in various species could be host-induced.

Diversity of Streptomyces spp. in Eastern Himalayan region - computational RNomics approach to phylogeny

Isolation and characterization of actinomycetes from soil samples from altitudinal gradient of North-East India were investigated for computational RNomics based phylogeny. A total of 52 diverse isolates of Streptomyces from the soil samples were isolated on four different media and from these 6 isolates were selected on the basis of cultural characteristics, microscopic and biochemical studies. Sequencing of 16S rDNA of the selected isolates identified them to belong to six different species of Streptomyces. The molecular morphometric and physico-kinetic analysis of 16S rRNA sequences were performed to predict the diversity of the genus. The computational RNomics study revealed the significance of the structural RNA based phylogenetic analysis in a relatively diverse group of Streptomyces.

The ancient history of the structure of ribonuclease P and the early origins of Archaea

BACKGROUND

Ribonuclease P is an ancient endonuclease that cleaves precursor tRNA and generally consists of a catalytic RNA subunit (RPR) and one or more proteins (RPPs). It represents an important macromolecular complex and model system that is universally distributed in life. Its putative origins have inspired fundamental hypotheses, including the proposal of an ancient RNA world.

RESULTS

To study the evolution of this complex, we constructed rooted phylogenetic trees of RPR molecules and substructures and estimated RPP age using a cladistic method that embeds structure directly into phylogenetic analysis. The general approach was used previously to study the evolution of tRNA, SINE RNA and 5S rRNA, the origins of metabolism, and the evolution and complexity of the protein world, and revealed here remarkable evolutionary patterns. Trees of molecules uncovered the tripartite nature of life and the early origin of archaeal RPRs. Trees of substructures showed molecules originated in stem P12 and were accessorized with a catalytic P1-P4 core structure before the first substructure was lost in Archaea. This core currently interacts with RPPs and ancient segments of the tRNA molecule. Finally, a census of protein domain structure in hundreds of genomes established RPPs appeared after the rise of metabolic enzymes at the onset of the protein world.

CONCLUSIONS

The study provides a detailed account of the history and early diversification of a fundamental ribonucleoprotein and offers further evidence in support of the existence of a tripartite organismal world that originated by the segregation of archaeal lineages from an ancient community of primordial organisms.

A coleopteran triosephosphate isomerase: X-ray structure and phylogenetic impact of insect sequences

A coleopteran triosephosphate isomerase (TIM) from Tenebrio molitor (yellow mealworm beetle) was recombinantly expressed in Escherichia coli and characterized with respect to thermal stability, kinetic parameters and oligomeric state. The enzyme was successfully crystallized and the structure determined by X-ray analysis to 2.0 A resolution. This is the first example of an invertebrate TIM. We compare structural features with known structures of TIMs from microorganisms, plants and vertebrates, and discuss the utility of the Tenebrio TIM sequence, together with several newly sequenced insect TIMs, for molecular phylogenetic analysis.

Phylogenetic reconstruction using secondary structures and sequence motifs of ITS2 rDNA of Paragonimus westermani (Kerbert, 1878) Braun, 1899 (Digenea: Paragonimidae) and related species

BACKGROUND

Most phylogenetic studies using current methods have focused on primary DNA sequence information. However, RNA secondary structures are particularly useful in systematics because they include characteristics that give "morphological" information, not found in the primary sequence. In several mountainous regions of Northeastern India, foci of Paragonimus (lung fluke) infection reportedly involve species that are known to prevail in neighbouring countries. The present study was undertaken to demonstrate the sequence analysis of the ribosomal DNA (ITS2) of the infective (metacercarial) stage of the lung fluke collected from the edible crab hosts that are abundant in a mountain stream of the area (Miao, Changlang District in Arunachal Pradesh) and to construct its phylogeny. Using the approach of molecular morphometrics that is based on ITS2 secondary structure homologies, phylogenetic relationships of the various isolates of Paragonimus species that are prevalent in the neighbouring Near-eastern countries have been discussed.

RESULTS

Initially, ten predicted RNA secondary structures were reconstructed and the topology based only on the predicted RNA secondary structure of the ITS2 region resolved most relationships among the species studied. We obtained three similar topologies for seven species of the genus Paragonimus on the basis of traditional primary sequence analysis using MEGA and a Bayesian analysis of the combined data. The latter approach allowed us to include both primary sequence and RNA molecular morphometrics; each data partition was allowed to have a different evolution rate. Paragonimus westermani was found to group with P. siamensis of Thailand; this was best supported by both the molecular morphometrics and combined analyses. P. heterotremus, P. proliferus, P. skrjabini, P. bangkokensis and P. harinasutai formed a separate clade in the molecular phylogenies, and were reciprocally monophyletic with respect to other species. ITS2 sequence motifs allowed an accurate in-silico distinction of lung flukes.

CONCLUSION

Data indicate that ITS2 motifs (<or= 50 bp in size) can be considered a promising tool for trematode species identification. RNA secondary structure analysis could be a valuable tool for distinguishing new species and completing Paragonimus systematics, more so because ITS2 secondary structure contains more information than the usual primary sequence alignment.

Use of sequence motifs as barcodes and secondary structures of internal transcribed spacer 2 (ITS2, rDNA) for identification of the Indian liver fluke, Fasciola (Trematoda: Fasciolidae)

Most phylogenetic studies using current methods have focused on primary DNA sequence information. However, RNA secondary structures are particularly useful in systematics because they include characteristics that give “morphological” information which is not found in the primary sequence. Also DNA sequence motifs from the internal transcribed spacer (ITS) of the nuclear rRNA repeat are useful for identification of trematodes. The species of liver flukes of the genus Fasciola (Platyhelminthes: Digenea: Fasciolidae) are obligate parasitic trematodes residing in the large biliary ducts of herbivorous mammals. While Fasciola hepatica has a cosmopolitan distribution, the other major species, i.e., F. gigantica is reportedly prevalent in the tropical and subtropical regions of Africa and Asia. To determine the Fasciola sp. of Assam (India) origin based on rDNA molecular data, ribosomal ITS2 region was sequenced (EF027103) and analysed. NCBI databases were used for sequence homology analysis and the phylogenetic trees were constructed based upon the ITS2 using MEGA and a Bayesian analysis of the combined data. The latter approach allowed us to include both primary sequence and RNA molecular morphometrics and revealed a close relationship with isolates of F. gigantica from China, Indonesia and Japan, the isolate from China with significant bootstrap values being the closest. ITS2 sequence motifs allowed an accurate in silico distinction of liver flukes. The data indicate that ITS2 motifs (≤ 50 bp in size) can be considered promising tool for trematode species identification. Using the novel approach of molecular morphometrics that is based on ITS2 secondary structure homologies, phylogenetic relationships of the various isolates of fasciolid species have been discussed.

Comparative analysis of sequences and secondary structures of the rRNA internal transcribed spacer 2 (ITS2) in pollen beetles of the subfamily Meligethinae (Coleoptera, Nitidulidae): potential use of slippage-derived sequences in molecular systematics

A comparative analysis of ITS2 sequences and secondary structures in 89 species of pollen beetles of the subfamily Meligethinae (Coleoptera, Nitidulidae) was performed. The ITS2 folding pattern was highly conserved and comparable with the general model proposed for eukaryotes. Simple sequence repeats (SSRs) were responsible for most of the observed nucleotide variability (approximately 1-3%) and length variation (359-459bp). When plotted on secondary structures, SSRs mapped in expansion segments positioned at the apices of three ITS2 helices ('A', 'B' and 'D1') and appeared to have evolved under mechanisms of compensatory slippage. Homologies among SSRs nucleotides could not be unambiguously assigned, and thus were not useful to resolve phylogeny. However, slippage-derived motifs provided some preliminary genetic support for newly proposed taxonomic arrangements of several genera and subgenera of Meligethinae, corroborating existing morphological and ecological datasets.

Phylogenetic comparison of local length plasticity of the small subunit of nuclear rDNAs among all Hexapoda orders and the impact of hyper-length-variation on alignment

The SSU nrDNA (18S), is one of the most frequently sequenced molecular markers in phylogenetic studies. However, the length-hyper-variation at multiple positions of this gene can affect the accuracy of alignment greatly and this length variation makes alignment across arthropod orders a serious problem. The analyses of Hexapoda phylogeny is such a case. A more clear recognition of the distribution of the length-variable-regions is needed. In this study, the secondary structure of some length-variable-regions in the SSU nrRNA of Arthropoda was adjusted by the principle of co-variation. It is found that the extent of plasticity of some length-variable-region can extraordinarily be higher than 600 bases in hexapods. And the numbers of hyper length-variable-regions are largest in Strepsiptera and Sternorrhyncha (Hemiptera). Our study shows that some length-variable-regions can serve as synapomorphies for some groups. The phylogenetic comparison also suggested that the expansion of a lateral bulge could be the origin of a helix.

Phylogenetic analysis of the internal transcribed spacer (ITS) region in Menyanthaceae using predicted secondary structure

Sequences of the nuclear internal transcribed spacer (ITS) regions ITS1 and ITS2 have been used widely in molecular phylogenetic studies because of their relatively high variability and facility of amplification. For phylogenetic applications, most researchers use sequence alignments that are based on nucleotide similarity. However, confidence in the alignment often deteriorates at taxonomic levels above genus, due to increasing variability among sequences. Like ribosomal RNA (rRNA) and other RNA molecules, the ITS transcripts consist in part of conserved secondary structures ('stems' and 'loops') that can be predicted by mathematical algorithm. Researchers have long considered the evolutionary conservation of rRNA secondary structure, but until recently few phylogenetic analyses of the ITS regions specifically incorporated structural data. We outline a novel method by which to derive additional phylogenetic data from ITS secondary structure in order to evaluate support for relationships at higher taxonomic levels. To illustrate the method, we describe an example from the plant family Menyanthaceae. Using predicted ITS secondary structure data, we obtained a well-resolved and moderately supported phylogeny, in which most topological relationships were congruent with the tree constructed using ITS nucleotide sequence data. Furthermore, the explicit encoding of ITS structural data in a phylogenetic framework allowed for the reconstruction of putative ancestral states and structural evolution in the functional but highly variable ITS region.

Evolutionary patterns in the sequence and structure of transfer RNA: a window into early translation and the genetic code

Transfer RNA (tRNA) molecules play vital roles during protein synthesis. Their acceptor arms are aminoacylated with specific amino acid residues while their anticodons delimit codon specificity. The history of these two functions has been generally linked in evolutionary studies of the genetic code. However, these functions could have been differentially recruited as evolutionary signatures were left embedded in tRNA molecules. Here we built phylogenies derived from the sequence and structure of tRNA, we forced taxa into monophyletic groups using constraint analyses, tested competing evolutionary hypotheses, and generated timelines of amino acid charging and codon discovery. Charging of Sec, Tyr, Ser and Leu appeared ancient, while specificities related to Asn, Met, and Arg were derived. The timelines also uncovered an early role of the second and then first codon bases, identified codons for Ala and Pro as the most ancient, and revealed important evolutionary take-overs related to the loss of the long variable arm in tRNA. The lack of correlation between ancestries of amino acid charging and encoding indicated that the separate discoveries of these functions reflected independent histories of recruitment. These histories were probably curbed by co-options and important take-overs during early diversification of the living world.

NOBAI: a web server for character coding of geometrical and statistical features in RNA structure

The Numeration of Objects in Biology: Alignment Inferences (NOBAI) web server provides a web interface to the applications in the NOBAI software package. This software codes topological and thermodynamic information related to the secondary structure of RNA molecules as multi-state phylogenetic characters, builds character matrices directly in NEXUS format and provides sequence randomization options. The web server is an effective tool that facilitates the search for evolutionary history embedded in the structure of functional RNA molecules. The NOBAI web server is accessible at 'http://www.manet.uiuc.edu/nobai/nobai.php'. This web site is free and open to all users and there is no login requirement.

Evolutionary patterns in the sequence and structure of transfer RNA: early origins of archaea and viruses

Transfer RNAs (tRNAs) are ancient molecules that are central to translation. Since they probably carry evolutionary signatures that were left behind when the living world diversified, we reconstructed phylogenies directly from the sequence and structure of tRNA using well-established phylogenetic methods. The trees placed tRNAs with long variable arms charging Sec, Tyr, Ser, and Leu consistently at the base of the rooted phylogenies, but failed to reveal groupings that would indicate clear evolutionary links to organismal origin or molecular functions. In order to uncover evolutionary patterns in the trees, we forced tRNAs into monophyletic groups using constraint analyses to generate timelines of organismal diversification and test competing evolutionary hypotheses. Remarkably, organismal timelines showed Archaea was the most ancestral superkingdom, followed by viruses, then superkingdoms Eukarya and Bacteria, in that order, supporting conclusions from recent phylogenomic studies of protein architecture. Strikingly, constraint analyses showed that the origin of viruses was not only ancient, but was linked to Archaea. Our findings have important implications. They support the notion that the archaeal lineage was very ancient, resulted in the first organismal divide, and predated diversification of tRNA function and specificity. Results are also consistent with the concept that viruses contributed to the development of the DNA replication machinery during the early diversification of the living world.

The origin and evolution of tRNA inferred from phylogenetic analysis of structure

The evolutionary history of the two structural and functional domains of tRNA is controversial but harbors the secrets of early translation and the genetic code. To explore the origin and evolution of tRNA, we reconstructed phylogenetic trees directly from molecular structure. Forty-two structural characters describing the geometry of 571 tRNAs and three statistical parameters describing thermodynamic and mechanical features of molecules quantitatively were used to derive phylogenetic trees of molecules and molecular substructures. Trees of molecules failed to group tRNA according to amino acid specificity and did not reveal the tripartite nature of life, probably due to loss of phylogenetic signal or because tRNA diversification predated organismal diversification. Trees of substructures derived from both structural and statistical characters support the origin of tRNA in the acceptor arm and the hypothesis that the top half domain composed of acceptor and pseudouridine (TPsiC) arms is more ancient than the bottom half domain composed of dihydrouridine (DHU) and anticodon arms. This constitutes the cornerstone of the genomic tag hypothesis that postulates tRNAs were ancient telomeres in the RNA world. The trees of substructures suggest a model for the evolution of the major functional and structural components of tRNA. In this model, short RNA hairpins with stems homologous to the acceptor arm of present day tRNAs were extended with regions homologous to TPsiC and anticodon arms. The DHU arm was then incorporated into the resulting three-stemmed structure to form a proto-cloverleaf structure. The variable region was the last structural addition to the molecular repertoire of evolving tRNA substructures.

Phylogenetic reconstruction using secondary structures of Internal Transcribed Spacer 2 (ITS2, rDNA): finding the molecular and morphological gap in Caribbean gorgonian corals

BACKGROUND

Most phylogenetic studies using current methods have focused on primary DNA sequence information. However, RNA secondary structures are particularly useful in systematics because they include characteristics, not found in the primary sequence, that give "morphological" information. Despite the number of recent molecular studies on octocorals, there is no consensus opinion about a region that carries enough phylogenetic resolution to solve intrageneric or close species relationships. Moreover, intrageneric morphological information by itself does not always produce accurate phylogenies; intra-species comparisons can reveal greater differences than intra-generic ones. The search for new phylogenetic approaches, such as by RNA secondary structure analysis, is therefore a priority in octocoral research.

RESULTS

Initially, twelve predicted RNA secondary structures were reconstructed to provide the basic information for phylogenetic analyses; they accorded with the 6 helicoidal ring model, also present in other groups of corals and eukaryotes. We obtained three similar topologies for nine species of the Caribbean gorgonian genus Eunicea (candelabrum corals) with two sister taxa as outgroups (genera Plexaura and Pseudoplexaura) on the basis of molecular morphometrics of ITS2 RNA secondary structures only, traditional primary sequence analyses and maximum likelihood, and a Bayesian analysis of the combined data. The latter approach allowed us to include both primary sequence and RNA molecular morphometrics; each data partition was allowed to have a different evolution rate. In addition, each helix was partitioned as if it had evolved at a distinct rate. Plexaura flexuosa was found to group within Eunicea; this was best supported by both the molecular morphometrics and combined analyses. We suggest Eunicea flexuosa (Lamouroux, 1821) comb. nov., and we present a new species description including Scanning Electron Microscopy (SEM) images of morphological characteristics (sclerites). Eunicea flexuosa, E. pallida, E. laxispica and E. mammosa formed a separate clade in the molecular phylogenies, and were reciprocally monophyletic with respect to other Eunicea (Euniceopsis subgenus, e.g. E. tourneforti and E. laciniata) in the molecular morphometrics tree, with the exception of E. fusca. Moreover, we suggest a new diagnostic character for Eunicea, also present in E. flexuosa: middle layer sclerites > 1 mm in length.

CONCLUSION

ITS2 was a reliable sequence for intrageneric studies in gorgonian octocorals because of the amount of phylogenetic signal, and was corroborated against morphological characters separating Eunicea from Plexaura. The ITS2 RNA secondary structure approach to phylogeny presented here did not rely on alignment methods such as INDELS, but provided clearly homologous characters for partition analysis and RNA molecular morphometrics. These approaches support the divergence of Eunicea flexuosa comb. nov. from the outgroup Plexaura, although it has been considered part of this outgroup for nearly two centuries because of morphological resemblance.

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.

Analysis of the secondary structure of ITS1 in Pectinidae: implications for phylogenetic reconstruction and structural evolution

It is at present difficult to accurately position gaps in sequence alignment and to determine substructural homology in structure alignment when reconstructing phylogenies based on highly divergent sequences. Therefore, we have developed a new strategy for inferring phylogenies based on highly divergent sequences. In this new strategy, the whole secondary structure presented as a string in bracket notation is used as phylogenetic characters to infer phylogenetic relationships. It is no longer necessary to decompose the secondary structure into homologous substructural components. In this study, reliable phylogenetic relationships of eight species in Pectinidae were inferred from the structure alignment, but not from sequence alignment, even with the aid of structural information. The results suggest that this new strategy should be useful for inferring phylogenetic relationships based on highly divergent sequences. Moreover, the structural evolution of ITS1 in Pectinidae was also investigated. The whole ITS1 structure could be divided into four structural domains. Compensatory changes were found in all four structural domains. Structural motifs in these domains were identified further. These motifs, especially those in D2 and D3, may have important functions in the maturation of rRNAs.

Application of the secondary structure model of rRNA for phylogeny: D2-D3 expansion segments of the LSU gene of plant-parasitic nematodes from the family Hoplolaimidae Filipjev, 1934

Knowledge of rRNA structure is increasingly important to assist phylogenetic analysis through reconstructing optimal alignment, utilizing molecule features as an additional source of data and refining appropriate models of evolution of the molecule. We describe a procedure of optimization for alignment and a new coding method for nucleotide sequence data using secondary structure models of the D2 and D3 expansion fragments of the LSU-rRNA gene reconstructed for fifteen nematode species of the agriculturally important and diverse family Hoplolaimidae, order Tylenchida. Using secondary structure information we converted the original sequence data into twenty-eight symbol codes and submitted the transformed data to maximum parsimony analysis. We also applied the original sequence data set for Bayesian inference. This used the doublet model with sixteen states of nucleotide doublets for the stem region and the standard model of DNA substitution with four nucleotide states for loops and bulges. By this approach, we demonstrate that using structural information for phylogenetic analyses led to trees with lower resolved relationships between clades and likely eliminated some artefactual support for misinterpreted relationships, such as paraphyly of Helicotylenchus or Rotylenchus. This study as well as future phylogenetic analyses is herein supported by the development of an on-line database, NEMrRNA, for rRNA molecules in a structural format for nematodes. We also have developed a new computer program, RNAstat, for calculation of nucleotide statistics designed and proposed for phylogenetic studies.

Higher phylogeny of zygaenid moths (Insecta: Lepidoptera) inferred from nuclear and mitochondrial sequence data and the evolution of larval cuticular cavities for chemical defence

Zygaenid moths are capable of releasing hydrogen cyanide in their defense by enzymatic break-down of cyanoglucosides, but only larvae of chalcosiine and zygaenine moths store cyanogenic compounds in cuticular cavities and thus are able to discharge defense droplets, which effectively deter potential predators. A previously proposed phylogeny of Zygaenidae hypothesized a sister group relationship of chalcosiine and zygaenine moths because of their similar larval defense system. Not all chalcosiine taxa possess cuticular cavities, however, and a comparable defense mechanism has been reported in larvae of the zygaenoid family Heterogynidae. Considering sequence data of seven molecular loci, the present study estimates the posterior probability of phylogenetic hypotheses explaining the occurrence of larval cuticular cavities. The molecular data confirm the previous exclusion of Himantopteridae from Zygaenidae and suggest their close affinity to Somabrachyidae. The sequence data also corroborate the recently proposed exclusion of the Phaudinae from the Zygaenidae, because this subfamily is recovered in a reasonably well supported species cluster consisting of members of the families Lacturidae, Limacodidae, Himantopteridae, and Somabrachyidae. We consequently agree to raise Phaudinae to family rank. Within Zygaenidae, the subfamilies Callizygaeninae, Chalcosiinae, and Procridinae most likely constitute a monophyletic group, which is sister to the Zygaeninae. Our results imply that cuticular cavities were probably present in the larvae of the most recent common ancestor of Zygaenidae. Heterogynidae cannot be confirmed as sister taxon to this family, but appear at the very first split of the Zygaenoidea, although with poor support. The specific pattern of taxa in the molecular phylogeny showing larval cuticular cavities opens the possibility that these structures could have been already present in the most recent common ancestor of the Zygaenoidea.

Parasite histories and novel phylogenetic tools: Alternative approaches to inferring parasite evolution from molecular markers

Parasitological research is often contingent on the knowledge of the phylogeny/genealogy of the studied group. Although molecular phylogenetics has proved to be a powerful tool in such investigations, its application in the traditional fashion, based on a tree inference from the primary nucleotide sequences may, in many cases, be insufficient or even improper. These limitations are due to a number of factors, such as a scarcity/ambiguity of phylogenetic information in the sequences, an intricacy of gene relationships at low phylogenetic levels, or a lack of criteria when deciding among several competing coevolutionary scenarios. With respect to the importance of a precise and reliable phylogenetic background in many biological studies, attempts are being made to extend molecular phylogenetics with a variety of new data sources and methodologies. In this review, selected approaches potentially applicable to parasitological research are presented and their advantages as well as drawbacks are discussed. These issues include the usage of idiosyncratic markers (unique features with presumably low probability of homoplasy), such as insertion of mobile elements, gene rearrangements and secondary structure features; the problem of ancestral polymorphism and reticulate relationships at low phylogenetic levels; and the utility of a molecular clock to facilitate discrimination among alternative scenarios in host-parasite coevolution.

A Hexapod nuclear SSU rRNA secondary-structure model and catalog of taxon-specific structural variation

RNA molecules and in particular the nuclear SSU RNA play an important role in molecular systematics. With the advent of increasingly parameterized substitution models in systematic research, the incorporation of secondary-structure information became a realistic option compensating interdependence of character variation. As a prerequisite, consensus structures of eukaryotic SSU RNA molecules have become available through extensive comparative analyses and crystallographic studies. Despite extensive research in hexapod phylogenetics, consensus SSU RNA secondary structures focusing on hexapods have not yet been explored. In this study, we compiled a representative hexapod SSU data set of 261 sequences and inferred a specific consensus SSU secondary-structure model. Our search for conserved structural motives relied on a combined approach of thermodynamic and covariation analyses. The hexapod consensus-structure model deviates from the canonical eukaryotic model in a number of helices. Additionally, in several helices the hexapod sequences did not support a single consensus structure. We provide consensus structures of these sections of single less-inclusive taxa, thus facilitating the adaptation of the consensus hexapod model to less-inclusive phylogenetic questions. The secondary-structure catalog will foster the application of RNA structure models in phylogenetic analyses using the SSU rRNA molecule, and it will improve the realism of substitution models and the reliability of reconstructions based on rRNA sequences.

Exploration of phylogenetic data using a global sequence analysis method

Background

Molecular phylogenetic methods are based on alignments of nucleic or peptidic sequences. The tremendous increase in molecular data permits phylogenetic analyses of very long sequences and of many species, but also requires methods to help manage large datasets.

Results

Here we explore the phylogenetic signal present in molecular data by genomic signatures, defined as the set of frequencies of short oligonucleotides present in DNA sequences. Although violating many of the standard assumptions of traditional phylogenetic analyses – in particular explicit statements of homology inherent in character matrices – the use of the signature does permit the analysis of very long sequences, even those that are unalignable, and is therefore most useful in cases where alignment is questionable. We compare the results obtained by traditional phylogenetic methods to those inferred by the signature method for two genes: RAG1, which is easily alignable, and 18S RNA, where alignments are often ambiguous for some regions. We also apply this method to a multigene data set of 33 genes for 9 bacteria and one archea species as well as to the whole genome of a set of 16 γ-proteobacteria. In addition to delivering phylogenetic results comparable to traditional methods, the comparison of signatures for the sequences involved in the bacterial example identified putative candidates for horizontal gene transfers.

Conclusion

The signature method is therefore a fast tool for exploring phylogenetic data, providing not only a pretreatment for discovering new sequence relationships, but also for identifying cases of sequence evolution that could confound traditional phylogenetic analysis.

A Secondary Structural Model of the 28S rRNA Expansion Segments D2 and D3 for Chalcidoid Wasps (Hymenoptera: Chalcidoidea)

We analyze the secondary structure of two expansion segments (D2, D3) of the 28S ribosomal (rRNA)-encoding gene region from 527 chalcidoid wasp taxa (Hymenoptera: Chalcidoidea) representing 18 of the 19 extant families. The sequences are compared in a multiple sequence alignment, with secondary structure inferred primarily from the evidence of compensatory base changes in conserved helices of the rRNA molecules. This covariation analysis yielded 36 helices that are composed of base pairs exhibiting positional covariation. Several additional regions are also involved in hydrogen bonding, and they form highly variable base-pairing patterns across the alignment. These are identified as regions of expansion and contraction or regions of slipped-strand compensation. Additionally, 31 single-stranded locales are characterized as regions of ambiguous alignment based on the difficulty in assigning positional homology in the presence of multiple adjacent indels. Based on comparative analysis of these sequences, the largest genetic study on any hymenopteran group to date, we report an annotated secondary structural model for the D2, D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related apocritan taxa.

A secondary structural model of the 28S rRNA expansion segments D2 and D3 from rootworms and related leaf beetles (Coleoptera: Chrysomelidae; Galerucinae)

We analysed the secondary structure of two expansion segments (D2, D3) of the 28S rRNA gene from 229 leaf beetles (Coleoptera: Chrysomelidae), the majority of which are in the subfamily Galerucinae. The sequences were compared in a multiple sequence alignment, with secondary structure inferred primarily from the compensatory base changes in the conserved helices of the rRNA molecules. This comparative approach yielded thirty helices comprised of base pairs with positional covariation. Based on these leaf beetle sequences, we report an annotated secondary structural model for the D2 and D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related beetle taxa. This predicted structure, consisting of seven major compound helices, is mostly consistent with previously proposed models for the D2 and D3 expansion segments in insects. Despite a lack of conservation in the primary structure of these regions of insect 28S rRNA, the evolution of the secondary structure of these seven major motifs may be informative above the nucleotide level for higher-order phylogeny reconstruction of major insect lineages.

Unraveling the evolutionary radiation of the thoracican barnacles using molecular and morphological evidence: a comparison of several divergence time estimation approaches

The Thoracica includes the ordinary barnacles found along the sea shore and is the most diverse and well-studied superorder of Cirripedia. However, although the literature abounds with scenarios explaining the evolution of these barnacles, very few studies have attempted to test these hypotheses in a phylogenetic context. The few attempts at phylogenetic analyses have suffered from a lack of phylogenetic signal and small numbers of taxa. We collected DNA sequences from the nuclear 18S, 28S, and histone H3 genes and the mitochondrial 12S and 16S genes (4,871 bp total) and data for 37 adult and 53 larval morphological characters from 43 taxa representing all the extant thoracican suborders (except the monospecific Brachylepadomorpha). Four Rhizocephala (highly modified parasitic barnacles) taxa and a Rhizocephala + Acrothoracica (burrowing barnacles) hypothetical ancestor were used as the outgroup for the molecular and morphological analyses, respectively. We analyzed these data separately and combined using maximum likelihood (ML) under "hill-climbing" and genetic algorithm heuristic searches, maximum parsimony procedures, and Bayesian inference coupled with Markov chain Monte Carlo techniques under mixed and homogeneous models of nucleotide substitution. The resulting phylogenetic trees answered key questions in barnacle evolution. The four-plated Iblomorpha were shown as the most primitive thoracican, and the plateless Heteralepadomorpha were placed as the sister group of the Lepadomorpha. These relationships suggest for the first time in an invertebrate that exoskeleton biomineralization may have evolved from phosphatic to calcitic. Sessilia (nonpedunculate) barnacles were depicted as monophyletic and appear to have evolved from a stalked (pedunculate) multiplated (5+) scalpelloidlike ancestor rather than a five-plated lepadomorphan ancestor. The Balanomorpha (symmetric sessile barnacles) appear to have the following relationship: (Chthamaloidea(Coronuloidea(Tetraclitoidea, Balanoidea))). Thoracican divergence times were estimated under ML-based local clock, Bayesian, and penalized likelihood approaches using an 18S data set and three calibration points: Heteralepadomorpha = 530 million years ago (MYA), Scalpellomorpha = 340 MYA, and Verrucomorpha = 120 MYA. Estimated dates varied considerably within and between approaches depending on the calibration point. Highly parameterized local clock models that assume independent rates (r > or = 15) for confamilial or congeneric species generated the most congruent estimates among calibrations and agreed more closely with the barnacle fossil record. Reasonable estimates were also obtained under the Bayesian procedure of Kishino et al. (2001, Mol. Biol. Evol. 18:352-361) but using multiple calibrations. Most of the dates estimated under the Bayesian procedure of Aris-Brosou and Yang (2002, Syst. Biol. 51:703-714) and the penalized likelihood method using single and/or multiple calibrations were inconsistent among calibrations and did not fit the fossil record.

Comparative analysis of mt LSU rRNA secondary structures of Odonates: structural variability and phylogenetic signal

Secondary structures of the most conserved part of the mt 16S rRNA gene, domains IV and V, have been recently analysed in a comparative study. However, full secondary structures of the mt LSU rRNA molecule are published for only a few insect species. The present study presents full secondary structures of domains I, II, IV and V of Odonates and one representative of mayflies, Ephemera sp. The reconstructions are based on a comparative approach and minimal consensus structures derived from sequence alignments. The inferred structures exhibit remarkable similarities to the published Drosophila melanogaster model, which increases confidence in these structures. Structural variance within Odonates is homoplastic, and neighbour-joining trees based on tree edit distances do not correspond to any of the phylogenetically expected patterns. However, despite homoplastic quantitative structural variation, many similarities between Odonates and Ephemera sp. suggest promising character sets for higher order insect systematics that merit further investigations.

Phylogenetic analyses among octocorals (Cnidaria): mitochondrial and nuclear DNA sequences (lsu-rRNA, 16S and ssu-rRNA, 18S) support two convergent clades of branching gorgonians

Gorgonian octocorals lack corroborated hypotheses of phylogeny. This study reconstructs genealogical relationships among some octocoral species based on published DNA sequences from the large ribosomal subunit of the mitochondrial RNA (lsu-rRNA, 16S: 524bp and 21 species) and the small subunit of the nuclear RNA (ssu-rRNA, 18S: 1815bp and 13 spp) using information from insertions-deletions (INDELS) and the predicted secondary structure of the lsu-rRNA (16S). There were seven short (3-10bp) INDELS in the 18S with consistent phylogenetic information. The INDELS in the 16S corresponded to informative signature sequences homologous to the G13 helix found in Escherichia coli. We found two main groups of gorgonian octocorals using a maximum parsimony analysis of the two genes. One group corresponds to deep-water taxa including species from the suborders Calcaxonia and Scleraxonia characterized by an enlargement of the G13 helix. The second group has species from Alcyoniina, Holaxonia and again Scleraxonia characterized by insertions in the 18S. Gorgonian corals, branching colonies with a gorgonin-containing flexible multilayered axis (Holaxonia and Calcaxonia), do not form a monophyletic group. These corroborated results from maternally inherited (16S) and biparentally inherited (18S) genes support a hypothesis of independent evolution of branching in the two octocoral clades.

Structural rRNA characters support monophyly of raptorial limbs and paraphyly of limb specialization in water fleas

The evolutionary success of arthropods has been attributed partly to the diversity of their limb morphologies. Large morphological diversity and increased specialization are observed in water flea (Cladocera) limbs, but it is unclear whether the increased limb specialization in different cladoceran orders is the result of shared ancestry or parallel evolution. We inferred a robust among-order cladoceran phylogeny using small-subunit and large-subunit rRNA nuclear gene sequences, signature sequence regions, novel stem-loops and secondary structure morphometrics to assess the phylogenetic distribution of limb specialization. The sequence-based and structural rRNA morphometric phylogenies were congruent and suggested monophyly of orders with raptorial limbs, but paraphyly of orders with reduced numbers of specialized limbs. These results highlight the utility of complex molecular structural characters in resolving ancient rapid radiations.

Louse (Insecta: Phthiraptera) mitochondrial 12S rRNA secondary structure is highly variable

Lice are ectoparasitic insects hosted by birds and mammals. Mitochondrial 12S rRNA sequences obtained from lice show considerable length variation and are very difficult to align. We show that the louse 12S rRNA domain III secondary structure displays considerable variation compared to other insects, in both the shape and number of stems and loops. Phylogenetic trees constructed from tree edit distances between louse 12S rRNA structures do not closely resemble trees constructed from sequence data, suggesting that at least some of this structural variation has arisen independently in different louse lineages. Taken together with previous work on mitochondrial gene order and elevated rates of substitution in louse mitochondrial sequences, the structural variation in louse 12S rRNA confirms the highly distinctive nature of molecular evolution in these insects.