Computational analysis of epigenetic information in human DNA sequences

Over the last few years, investigations of human epigenetic profiles have identified key elements of change to be Histone Modifications, stable and heritable DNA methylation and Chromatin remodeling. These factors determine gene expression levels and characterise conditions leading to disease. In order to extract information embedded in long DNA sequences, data mining and pattern recognition tools are widely used, but efforts have been limited to date with respect to analyzing epigenetic changes, and their role as catalysts in disease onset. Useful insight, however, can be gained by investigation of associated dinucleotide distributions. The focus of this paper is to explore specific dinucleotides frequencies across defined regions within the human genome, and to identify new patterns between epigenetic mechanisms and DNA content. Signal processing methods, including Fourier and Wavelet Transformations, are employed and principal results are reported.

Mitochondrial DNA sequence analysis from multiple gene fragments reveals genetic heterogeneity of Crassostrea ariakensis in East Asia

The native Asian oyster, Crassostrea ariakensis is one of the most common and important Crassostrea species that occur naturally along the coast of East Asia. Molecular species diagnosis is a prerequisite for population genetic analysis of wild oyster populations because oyster species cannot be discriminated reliably using external morphological characters alone due to character ambiguity. To date there have been few phylogeographic studies of natural edible oyster populations in East Asia, in particular this is true of the common species in Korea C. ariakensis. We therefore assessed the levels and patterns of molecular genetic variation in East Asian wild populations of C. ariakensis from Korea, Japan, and China using DNA sequence analysis of five concatenated mtDNA regions namely; 16S rRNA, cytochrome oxidase I, cytochrome oxidase II, cytochrome oxidase III, and cytochrome b. Two divergent C. ariakensis clades were identified between southern China and remaining sites from the northern region. In addition, hierarchical AMOVA and pairwise UST analyses showed that genetic diversity was discontinuous among wild populations of C. ariakensis in East Asia. Biogeographical and historical sea level changes are discussed as potential factors that may have influenced the genetic heterogeneity of wild C. ariakensis stocks across this region.

Analysis of a Genomic DNA Expression Library of Mycobacterium tuberculosis Using Tuberculosis Patient Sera: Evidence for Modulation of Host Immune Response

DNA obtained from a human sputum isolate of Mycobacterium tuberculosis, NTI-64719, which showed extensive dissemination in the guinea pig model resulting in a high score for virulence was used to construct an expression library in the lambda ZAP vector. The size of DNA inserts in the library ranged from 1 to 3 kb, and recombinants represented 60% of the total plaques obtained. When probed with pooled serum from chronically infected tuberculosis patients, the library yielded 176 recombinants with a range of signal intensities. Among these, 93 recombinants were classified into 12 groups on the basis of DNA hybridization experiments, The polypeptides synthesized by the recombinants were predominantly LacZ fusion proteins, Serum obtained from patients who were clinically diagnosed to be in the early phase of M. tuberculosis infection was used to probe the 176 recombinants obtained. interestingly, some recombinants that gave very strong signals in the original screen did not react with early-phase serum; conversely, others whose signals were extremely weak in the original screen gave very intense signals with serum from recently infected patients, This indicates the differential nature of either the expression of these antigens or the immune response elicited by them as a function of disease progression.

Crystalline A-DNA: The X-Ray Analysis of the Fragment d(G-G-T-A-T-A-C-C)

An A-DNA type double helical conformation was observed in the single crystal X-ray structure of the octamer d(G-G-T-A-T-A-C-C), 1, and its 5-bromouracil-containing analogue, 2. The structure of the isomorphous crystals (space group P61) was solved by a search technique based on packing criteria and R-factor calculations, with use of only low order data. At the present stage of refinement the R factors are 31 % for 1 and 28 % for 2 at a resolution of 2.25 A (0.225 nm). The molecules interact through their minor grooves by hydrogen bonding and base to sugar van der Waals contacts. The stable A conformation observed in the crystal may have some structural relevance to promoter regions where the T-A-T-A sequence is frequently found.

DNA amplification fingerprinting analysis of genetic variation within Fusarium oxysporum f.sp zingiberi

Genetic variation among 29 isolates of Fusarium oxysporum f.sp. zingiberi (Foz) collected from diseased ginger rhizome in production regions throughout Queensland was analysed using DNA amplification fingerprinting (DAF). Eight isolates of other Fusarium species and/or formae speciales were included for comparative analysis. Within the Foz isolates, three haplotypes were identified based on 17 polymorphic bands generated with five primers. Two groups showed very little genetic variation (98.6% similarity), whereas the third single isolate was quite distinct in terms of its molecular profile (77.2% similarity). Genetic similarity among the Fusarium solani, F. oxysporum f.sp. lycopersici and F. oxysporum f.sp. cubense races 1, 3 and 4 isolates compared well with the published literature.

High-resolution melt-curve analysis of random-amplified-polymorphic-DNA markers, for the characterisation of pathogenic Leptospira.

A new test for pathogenic Leptospira isolates, based on RAPD-PCR and high-resolution melt (HRM) analysis (which measures the melting temperature of amplicons in real time, using a fluorescent DNA-binding dye), has recently been developed. A characteristic profile of the amplicons can be used to define serovars or detect genotypes. Ten serovars, of leptospires from the species Leptospira interrogans (serovars Australis, Robinsoni, Hardjo, Pomona, Zanoni, Copenhageni and Szwajizak), L. borgpetersenii (serovar Arborea), L. kirschneri (serovar Cynopteri) and L. weilii (serovar Celledoni), were typed against 13 previously published RAPD primers, using a real-time cycler (the Corbett Life Science RotorGene 6000) and the optimised reagents from a commercial kit (Quantace SensiMix). RAPD-HRM at specific temperatures generated defining amplicon melt profiles for each of the tested serovars. These profiles were evaluated as difference-curve graphs generated using the RotorGene software package, with a cut-off of at least 8 'U' (plus or minus). The results demonstrated that RAPD-HRM can be used to measure serovar diversity and establish identity, with a high degree of stability. The characterisation of Leptospira serotypes using a DNA-based methodology is now possible. As an objective and relatively inexpensive and rapid method of serovar identification, at least for cultured isolates, RAPD-HRM assays show convincing potentia.

High-resolution melt-curve analysis of random amplified polymorphic DNA (RAPD-HRM) for the characterisation of pathogenic leptospires: intra-serovar divergence, interserovar convergence, and evidence of attenuation in Leptospira reference collections.

High-resolution melt-curve analysis of random amplified polymorphic DNA (RAPD-HRM) is a novel technology that has emerged as a possible method to characterise leptospires to serovar level. RAPD-HRM has recently been used to measure intra-serovar convergence between strains of the same serovar as well as inter-serovar divergence between strains of different serovars. The results indicate that intra-serovar heterogeneity and inter-serovar homogeneity may limit the application of RAPD-HRM in routine diagnostics. They also indicate that genetic attenuation of aged, high-passage-number isolates could undermine the use of RAPD-HRM or any other molecular technology. Such genetic attenuation may account for a general decrease seen in titres of rabbit hyperimmune antibodies over time. Before RAPD-HRM can be further advanced as a routine diagnostic tool, strains more representative of the wild-type serovars of a given region need to be identified. Further, RAPD-HRM analysis of reference strains indicates that the routine renewal of reference collections, with new isolates, may be needed to maintain the genetic integrity of the collections.

Comparative analysis of Tritrichomonas foetus (Riedmuller, 1928) cat genotype, T. foetus (Riedmuller, 1928) cattle genotype and Tritrichomonas suis (Davaine, 1875) at 10 DNA loci

The parasitic protists in the genus Tritrichomonas cause significant disease in domestic cattle and cats. To assess the genetic diversity of feline and bovine isolates of Tritrichomonas foetus (Riedmuller, 1928) Wenrich and Emmerson, 1933, we used 10 different genetic regions, namely the protein coding genes of cysteine proteases 1,2 and 4-9 (CP1, 2, 4-9) involved in the pathogenesis of the disease caused by the parasite. The cytosolic malate dehydrogenase 1 (MDH1) and internal transcribed spacer region 2 of the rDNA unit (ITS2) were included as additional markers. The gene sequences were compared with those of Tritrichomonas suis (Davaine. 1875) Morgan and Hawkins, 1948 and Tritrichomonas mobilensis Culberson et al., 1986. The study revealed 100% identity for all 10 genes among all feline isolates (=T. foetus cat genotype), 100% identity among all bovine isolates (=T. foetus cattle genotype) and a genetic distinctness of 1% between the cat and cattle genotypes of T. foetus. The cattle genotype of T. foetus was 100% identical to T. suis at nine loci (CP1, 2,4-8, ITS2, MDH1). At CP9, three out of four T. suis isolates were identical to the T. foetus cattle genotype, while the T. suis isolate SUI-H3B sequence contained a single unique nucleotide substitution. Tritrichomonas mobilensis was 0.4% and 0.7% distinct from the cat and cattle genotypes of T. foetus, respectively. The genetic differences resulted in amino acid changes in the CP genes, most pronouncedly in CP2, potentially providing a platform for elucidation of genotype-specific host-pathogen interactions of T. foetus. On the basis of this data we judge T. suis and T. foetus to be subjective synonyms. For the first time, on objective nomenclatural grounds, the authority of T. suis is given to Davaine, 1875, rather than the commonly cited Gruby and Delafond, 1843. To maintain prevailing usage of T. foetus, we are suppressing the senior synomym T. suis Davaine, 1875 according to Article 23.9, because it has never been used as a valid name after 1899 and T. foetus is widely discussed as the cause of bovine trichomonosis. Thus bovine, feline and porcine isolates should all be given the name T. foetus. This promotes the stability of T. foetus for the veterinary and economically significant venereal parasite causing bovine trichomonosis. (C) 2012 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Sequence Specific Interaction between Promoter DNA and Escherichia coli RNA Polymerase: Comparative Thermodynamic Analysis with One Immobilized Partner

Sequence specific interaction between DNA and protein molecules has been a subject of active investigation for decades now. Here, we have chosen single promoter containing bacteriophage Delta D-III T7 DNA and Escherichia coli RNA polymerase and followed their recognition at the air-water interface by using the surface plasmon resonance (SPR) technique, where the movement of one of the reacting species is restricted by way of arraying them on an immobilized support. For the Langmuir monolayer studies, we used a RNA polymerase with a histidine tag attached to one of its subunits, thus making it an xcellent substrate for Ni(II) ions, while the SPR Studies were done using biotin-labeled DNA immobilized on a streptavidin-coated chip. Detailed analysis of the thermodynamic parameters as a function of concentration and temperature revealed that the interaction of RNA polymerase with T7 DNA is largely entropy driven (83 (+/- 12) kcal mol(-1)) with a positive enthalpy of 13.6 (+/- 3.6) kcal mol(-1), The free energy of reaction determined by SPR and Langmuir-Blodgett technique was -11 (+/- 2) and -15.6 kcal mol(-1), respectively. The ability of these methods to retain the specificity of the recognition process was also established.

Mutational analysis of active-site residues in the Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease: Asp(122) and Asp(193) are crucial to the double-stranded DNA cleavage activity whereas Asp(218) is not

Mycobacterium leprae recA harbors an in-frame insertion sequence that encodes an intein homing endonuclease (PI-MleI). Most inteins (intein endonucleases) possess two conserved LAGLIDADG (DOD) motifs at their ctive center. A common feature of LAGLIDADG-type homing endonucleases is that they recognize and cleave the same or very similar DNA sequences. However, PI-MleI is distinctive from other members of the family of LAGLIDADG-type HEases for its modular structure with functionally separable domains for DNA-binding and cleavage, each with distinct sequence preferences. Sequence alignment analyses of PI-MleI revealed three putative LAGLIDADG motifs; however, there is conflicting bioinformatics data in regard to their identity and specific location within the intein polypeptide. To resolve this conflict and to determine the active-site residues essential for DNA target site recognition and double-stranded DNA cleavage, we performed site-directed mutagenesis of presumptive catalytic residues in the LAGLIDADG motifs. Analysis of target DNA recognition and kinetic parameters of the wild-type PI-MleI and its variants disclosed that the two amino acid residues, Asp(122) (in Block C) and Asp(193) (in functional Block E), are crucial to the double-stranded DNA endonuclease activity, whereas Asp(218) (in pseudo-Block E) is not. However, despite the reduced catalytic activity, the PI-MleI variants, like the wild-type PI-MleI, generated a footprint of the same length around the insertion site. The D122T variant showed significantly reduced catalytic activity, and D122A and D193A mutations although failed to affect their DNA-binding affinities, but abolished the double-stranded DNA cleavage activity. On the other hand, D122C variant showed approximately twofold higher double-stranded DNA cleavage activity, compared with the wild-type PI-MleI. These results provide compelling evidence that Asp(122) and Asp(193) in DOD motif I and II, respectively, are bona fide active-site residues essential for DNA cleavage activity. The implications of these results are discussed in this report.

Analysis of DNA sequence transformations on grids

Study of the evolution of species or organisms is essential for various biological applications. Evolution is typically studied at the molecular level by analyzing the mutations of DNA sequences of organisms. Techniques have been developed for building phylogenetic or evolutionary trees for a set of sequences. Though phylogenetic trees capture the overall evolutionary relationships among the sequences, they do not reveal fine-level details of the evolution. In this work, we attempt to resolve various fine-level sequence transformation details associated with a phylogenetic tree using cellular automata. In particular, our work tries to determine the cellular automata rules for neighbor-dependent mutations of segments of DNA sequences. We also determine the number of time steps needed for evolution of a progeny from an ancestor and the unknown segments of the intermediate sequences in the phylogenetic tree. Due to the existence of vast number of cellular automata rules, we have developed a grid system that performs parallel guided explorations of the rules on grid resources. We demonstrate our techniques by conducting experiments on a grid comprising machines in three countries and obtaining potentially useful statistics regarding evolutions in three HIV sequences. In particular, our work is able to verify the phenomenon of neighbor-dependent mutations and find that certain combinations of neighbor-dependent mutations, defined by a cellular automata rule, occur with greater than 90% probability. We also find the average number of time steps for mutations for some branches of phylogenetic tree over a large number of possible transformations with standard deviations less than 2.

Analysis of the impact of a uracil DNA glycosylase attenuated in AP-DNA binding in maintenance of the genomic integrity in Escherichia coli

Uracil DNA glycosylase (Ung)initiates the uracil excision repair pathway. We have earlier characterized the Y66W and Y66H mutants of Ung and shown that they are compromised by similar to 7- and similar to 170-fold, respectively in their uracil excision activities. In this study, fluorescence anisotropy measurements show that compared with the wild-type, the Y66W protein is moderately compromised and attenuated in binding to AP-DNA. Allelic exchange of ung in Escherichia coli with ung::kan, ungY66H:amp or ungY66W:amp alleles showed similar to 5-, similar to 3.0- and similar to 2.0-fold, respectively increase in mutation frequencies. Analysis of mutations in the rifampicin resistance determining region of rpoB revealed that the Y66W allele resulted in an increase in A to G (or T to C) mutations. However, the increase in A to G mutations was mitigated upon expression of wild-type Ung from a plasmid borne gene. Biochemical and computational analyses showed that the Y66W mutant maintains strict specificity for uracil excision from DNA. Interestingly, a strain deficient in AP-endonucleases also showed an increase in A to G mutations. We discuss these findings in the context of a proposal that the residency of DNA glycosylase(s) onto the AP-sites they generate shields them until recruitment of AP-endonucleases for further repair.

A self-consistent formulation for analysis and generation of non-uniform DNA structure

A fully self-consistent formulation is described here for the analysis and generation of base-pairs in non-uniform DNA structures, in terms of various local parameters. It is shown that the internal "wedge parameters" are mathematically related to the parameters describing the base-pair orientation with respect to an external helix axis. Hence any one set of three translation and three rotation parameters are necessary and sufficient to completely describe the relative orientation of the base-pairs comprising a step (or doublet). A general procedure is outlined for obtaining an average or global helix axis from the local helix axes for each step. A graphical representation of the local helix axes in the form of a polar plot is also shown and its application for estimating the curvature of oligonucleotide structures is illustrated, with examples of both A and B type structures.

Polypurine-polypyrimidine sequences adopt unwound structure in pBR322 form V DNA as probed by single-hit analysis of HpaII sites.

Structure at the polypurine-polypyrimidine sequences flanking the HpaII sites (CCGG) in pBR322 form V DNA was probed employing single-hit analysis using HpaII restriction endonuclease. Reduced cleavage efficiency of HpaII sites flanked by polypurine-polypyrimidine sequences suggested that under high torsional stress these sequences adopt unwound structures rendering these sites insensitive to restriction enzyme cleavage. In addition to polypurine-polypyrimidine sequences. HpaII sites flanked by alternating purine-pyrimidine sequence, a potential motif of left handed Z-DNA, were also found to be resistant to HpaII cleavage. Results obtained from various studies implicating structure sensitivity of restriction endonucleases and methylases were compiled and a direct correlation was observed between the occurrence of altered sites in a domain and its G/C content in pBR322 form V DNA.