DNA Sequence Representation methods

DNA sequence representation methods are used to denote a gene structure effectively and help in similarities/dissimilarities analysis of coding sequences. Many different kinds of representations have been proposed in the literature. They can be broadly classified into Numerical, Graphical, Geometrical and Hybrid representation methods. DNA structure and function analysis are made easy with graphical and geometrical representation methods since it gives visual representation of a DNA structure. In numerical method, numerical values are assigned to a sequence and digital signal processing methods are used to analyze the sequence. Hybrid approaches are also reported in the literature to analyze DNA sequences. This paper reviews the latest developments in DNA Sequence representation methods. We also present a taxonomy of various methods. A comparison of these methods where ever possible is also done

Ultrastructure and small-subunit ribosomal DNA sequence of Henneguya lesteri n. sp (Myxosporea), a parasite of sand whiting Sillago analis (Sillaginidae) from the coast of Queensland, Australia

Henneguya lesteri n. sp, (Myxosporea) is described from sand whiting, Sillago analis, from the southern Queensland coast of Australia. H. lesteri displays a preference for the pseudobranchs and is typically positioned along the afferent blood vessels, displacing the adjoining lamellae and disrupting their normal array, The plasmodia appeared as whitish-hyaline, elliptical cysts (mean dimensions 230 x 410 mum) attached to the oral mucosa lining of the hyoid arch on the inner surface of the operculum. Infections of the gills were also found, in which the plasmodia were spherical, averaged 240 x 240 mum in size and were located on the inner hemibranch margin. The parasites lodged in the gill filament crypts and generated a mild hyperplastic response of the branchial epithelium, In histological sections, the plasmodium wall and adjoining ectoplasm appeared as a finely granulated, weakly eosinophilic layer, Ultrastructurally, this section of the host-parasite interface contained an intricate complex of pinocytotic channels. H. lesteri is polysporic, disporoblastic and pansporoblast forming. Sporogenesis is asynchronous, with the earliest developmental stages aligned predominantly along the plasmodium periphery, and maturing sporoblasts and spores toward the center. Ultrastructural details of sporoblast and spore development are in agreement with previously described myxosporeans. The mature spore is drop-shaped, length (mean) 9.1 mum, width 4.7 mum, thickness 2.5 mum, and comprises 2 polar capsules positioned closely together, a binucleated sporoplasm and a caudal process of 12.6 mum. The polar capsules are elongated, 3.2 x 1.6 mum, with 4 turns of the polar filament. Mean length of the everted filament is 23.2 mum, Few studies have analyzed the 18S gene-of marine Myxosporea. In fact, H. lesteri is the first marine species of Henneguya to be characterized at the molecular level: we determined 1966 bp of the small-subunit (18S) rDNA, The results indicated that differences between this and the hitherto studied freshwater Henneguya species are greater than differences among the freshwater Henneguya species.

Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda)

The complete nucleotide sequence of the mitochondrial (mt) DNA molecule of the liverfluke, Fasciola hepatica (phylum Platyhelminthes, class Trematoda, family Fasciolidae), was determined, It comprises 14462 bp, contains 12 protein-encoding, 2 ribosomal and 22 transfer RNA genes, and is the second complete flatworm (and the first trematode) mitochondrial sequence to be described in detail. All of the genes are transcribed from the same strand. Of the genes typically found in mitochondrial genomes of eumetazoans, only atp8 is absent. The nad4L and nad4 genes overlap by 40 nt. Most intergenic sequences are very short. Two larger non-coding regions are present. The longer one (817 nt) is located between trnG and cox3 and consists of 8 identical tandem repeats of 85 nt, rich in G and C, followed by 1 imperfect repeat. The shorter non-coding region (187 nt) exhibits no special features and is separated from the longer region by trnG. The gene arrangement resembles that of some other trematodes including the eastern Asian Schistosoma species (and cyclophyllidean cestode species) but it is strikingly different from that of the African schistosomes, represented by Schistosoma mansoni. The genetic code is as inferred previously for flatworms. Transfer RNA genes range in length from 58 to 70 nt, their products producing characteristic 'clover leaf' structures, except for tRNA(S-VON) and tRNA(S-AGN) lacking the DHU arm.

DNA sequence variation in the ITS-1 rDNA subunit and host relationships in sorghum midge, Stenodiplosis sorghicola (Coquillett) (Diptera : Cecidomyiidae), in Australia

Sequence variation in the internal transcribed spacer (ITS-1) ribosomal DNA subunit was examined for sorghum midge obtained from introduced and native hosts in south-eastern and central Queensland. No variation was observed relative to host plant or geographical distance for midges collected from two introduced hosts, grain sorghum (Sorghum bicolor ) and Johnson grass (S. halepense ); however, sequence differences were observed between midges from introduced and native hosts and among midges from a single native host, slender bluegrass (Dichanthium affine ). No evidence was observed of introduced midges on native hosts, or vice versa. These results agree with previously hypothesised host distributions for native and introduced midges in Australia, and expand the sample of introduced hosts to include Johnson grass. They suggest that Stenodiplosis sorghicola , the principal midge infesting grain sorghum, is also the most common species on Johnson grass. This confirms that Johnson grass plays a role in the population dynamics of S. sorghicola and suggests that midges originating from Johnson grass may influence levels of infestation in grain sorghum.

Are there two species of Pygmy Shrews (Sorex)? Revisiting the question using DNA sequence data

Pygmy Shrews in North America have variously been considered to be one species (Sorex hoyi) or two species (S. hoyi and S. thompsoni). Currently, only S. hoyi is recognized. In this study, we examine mitochondrial DNA sequence data for the cytochrome b gene to evaluate the level of differentiation and phylogeographic relationships among eleven samples of Pygmy Shrews from across Canada. Pygmy Shrews from eastern Canada (i.e., Ontario, Quebec, New Brunswick, Nova Scotia, and Prince Edward Island) are distinct from Pygmy Shrews from western Canada (Alberta, Yukon) and Alaska. The average level of sequence divergence between these clades (3.3%) falls within the range of values for other recognized pairs of sister species of shrews. A molecular clock based on third position transversion substitutions suggests that these two lineages diverged between 0.44 and 1.67 million years ago. These molecular phylogenetic data. combined with a reinterpretation of previously published morphological data, are suggestive of separate species status for S. hoyi and S. thompsoni as has been previously argued by others. Further analysis of specimens from geographically intermediate areas (e.g., Manitoba. northern Ontario) is required to determine if there is secondary contact and/or introgression between these two putative species.

A new repetitive DNA sequence from Trypanosoma cruzi

Tandemly repeated DNA sequences are found in the genome of higher eukaryotes, and have also been demonstrated in Trypanosoma cruzi. Repeated DNA sequences are potentially useful for the diagnostic detection of T. cruzi (A. Gonzales et al., 1984, Proc. Natl. Acad. Sci. USA, 81: 3356-3360). We have isoleted two clones from a genomic library of T. cruzi (Y strain) that contain, in one clone a family of at least seven copies of a repetitive sequence of approximately 600 base pairs, and in the other an independent copy of the same sequence. One copy of the repetition (HSP) and the independent clone (HCR) were sequenced by the Sanger procedure (Fig.). This sequence hybridized to four strains of T. cruzi tested and did not hybridize to eleven species of trypanosotids from five different Genera, being a good candidate for diagnostic assays. GenBank accession numbers: HSP#m31919, HCR#31920.

Mitochondrial DNA sequence variation within the remnant populations of the endangered numbat (Marsupialia: Myrmecobiidae: Myrmecobius fasciatus).

The numbat has been reduced to two populations in Western Australia. To better understand the effects of range reduction on gene flow and genetic variation, and to address questions crucial for the species' management, we analysed mitochondrial DNA (mtDNA) sequences of free-ranging individuals and museum specimens. The results suggest recent connectivity between the remnant populations, although one of those may have lost significant amounts of genetic diversity during the recent population size reduction. We propose that for management purposes the remnant populations should be treated as a single historical lineage and that, subject to certain caveats, consideration should be given to population augmentation by translocation.

CA88, a nuclear repetitive DNA sequence identified in Schistosoma mansoni, aids in the genotyping of nine Schistosoma species of medical and veterinary importance

CA88 is the first long nuclear repetitive DNA sequence identified in the blood fluke, Schistosoma mansoni. The assembled S. mansoni sequence, which contains the CA88 repeat, has 8,887 nucleotides and at least three repeat units of approximately 360 bp. In addition, CA88 also possesses an internal CA microsatellite, identified as SmBr18. Both PCR and BLAST analysis have been used to analyse and confirm the CA88 sequence in other S. mansoni sequences in the public database. PCR-acquired nuclear repetitive DNA sequence profiles from nine Schistosoma species were used to classify this organism into four genotypes. Included among the nine species analysed were five sequences of both African and Asian lineages that are known to infect humans. Within these genotypes, three of them refer to recognised species groups. A panel of four microsatellite loci, including SmBr18 and three previously published loci, has been used to characterise the nine Schistosoma species. Each species has been identified and classified based on its CA88 DNA fingerprint profile. Furthermore, microsatellite sequences and intra-specific variation have also been observed within the nine Schistosoma species sequences. Taken together, these results support the use of these markers in studying the population dynamics of Schistosoma isolates from endemic areas and also provide new methods for investigating the relationships between different populations of parasites. In addition, these data also indicate that Schistosoma magrebowiei is not a sister taxon to Schistosoma mattheei, prompting a new designation to a basal clade.

Estimation of protein coding density in a corpus of DNA sequence data

A number of experimental methods have been reported for estimating the number of genes in a genome, or the closely related coding density of a genome, defined as the fraction of base pairs in codons. Recently, DNA sequence data representative of the genome as a whole have become available for several organisms, making the problem of estimating coding density amenable to sequence analytic methods. Estimates of coding density for a single genome vary widely, so that methods with characterized error bounds have become increasingly desirable. We present a method to estimate the protein coding density in a corpus of DNA sequence data, in which a ‘coding statistic’ is calculated for a large number of windows of the sequence under study, and the distribution of the statistic is decomposed into two normal distributions, assumed to be the distributions of the coding statistic in the coding and noncoding fractions of the sequence windows. The accuracy of the method is evaluated using known data and application is made to the yeast chromosome III sequence and to C.elegans cosmid sequences. It can also be applied to fragmentary data, for example a collection of short sequences determined in the course of STS mapping.

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

The Vi element. A transposon-like repeated DNA sequence interspersed in the vitellogenin locus of Xenopus laevis.

A repeated DNA element in Xenopus laevis is described that is present in about 7500 copies dispersed throughout the genome. It was first identified in the 5' flanking region of one vitellogenin gene and was therefore named the Vi element. Seven copies are present within the vitellogenin gene region, three of them within introns of the genes A1, A2 and B2, and the other four copies in the gene flanking regions. Four of these copies have been sequenced. The Vi element is bounded by a well-conserved 13 base-pair inverted repeat; in addition, it is flanked by a three base-pair direct repeat that appears to be site-specific. The length of these four copies varies from 112 to 469 base-pairs; however, sequence homology between the different copies is very high. Their structural characteristics suggest that length heterogeneity may have arisen by either unequal recombinations, deletions or tandem duplications. Altogether, the characteristics and properties of the Vi element indicate that it might represent a mobile genetic element. One of the four copies sequenced is inserted close (position -535) to the transcription initiation site of the vitellogenin gene B2 in a region otherwise showing considerable homology with the closely related gene B1. Nevertheless, the presence of the Vi element does not seem to influence significantly the estrogen-controlled expression of gene B2. In addition, three alleles of this gene created by length polymorphism in intron 3 and in the Vi element inserted near the transcription initiation site are described.

The estrogen-responsive element as an inducible enhancer: DNA sequence requirements and conversion to a glucocorticoid-responsive element.

The estrogen-responsive element (ERE) present in the 5'-flanking region of the Xenopus laevis vitellogenin (vit) gene B1 has been characterized by transient expression analysis of chimeric vit-tk-CAT (chloramphenicol acetyltransferase) gene constructs transfected into the human estrogen-responsive MCF-7 cell line. The vit B1 ERE behaves like an inducible enhancer, since it is able to confer estrogen inducibility to the heterologous HSV thymidine kinase (tk) promoter in a relative position- and orientation-independent manner. In this assay, the minimal B1 ERE is 33 bp long and consists of two 13 bp imperfect palindromic elements both of which are required for the enhancer activity. A third imperfect palindromic element is present further upstream within the 5'-flanking region of the gene but is unable to confer hormone responsiveness by itself. Similarly, neither element forming the B1 ERE can alone confer estrogen inducibility to the tk promoter. However, in combinations of two, all three imperfect palindromes can act cooperatively to form a functional ERE. In contrast a single 13 bp perfect palindromic element, GGTCACTGTGACC, such as the one found upstream of the vit gene A2, is itself sufficient to act as a fully active ERE. Single point mutations within this element abolish estrogen inducibility, while a defined combination of two mutations converts this ERE into a glucocorticoid-responsive element.

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data.

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.