Discovering Frequent Poly-Regions in DNA Sequences

The problem of discovering frequent arrangements of regions of high occurrence of one or more items of a given alphabet in a sequence is studied, and two efficient approaches are proposed to solve it. The first approach is entropy-based and uses an existing recursive segmentation technique to split the input sequence into a set of homogeneous segments. The key idea of the second approach is to use a set of sliding windows over the sequence. Each sliding window keeps a set of statistics of a sequence segment that mainly includes the number of occurrences of each item in that segment. Combining these statistics efficiently yields the complete set of regions of high occurrence of the items of the given alphabet. After identifying these regions, the sequence is converted to a sequence of labeled intervals (each one corresponding to a region). An efficient algorithm for mining frequent arrangements of temporal intervals on a single sequence is applied on the converted sequence to discover frequently occurring arrangements of these regions. The proposed algorithms are tested on various DNA sequences producing results with significant biological meaning.

Separating DNA with different topologies by atomic force microscopy in comparison with gel electrophoresis.

Atomic force microscopy, which is normally used for DNA imaging to gain qualitative results, can also be used for quantitative DNA research, at a single-molecular level. Here, we evaluate the performance of AFM imaging specifically for quantifying supercoiled and relaxed plasmid DNA fractions within a mixture, and compare the results with the bulk material analysis method, gel electrophoresis. The advantages and shortcomings of both methods are discussed in detail. Gel electrophoresis is a quick and well-established quantification method. However, it requires a large amount of DNA, and needs to be carefully calibrated for even slightly different experimental conditions for accurate quantification. AFM imaging is accurate, in that single DNA molecules in different conformations can be seen and counted. When used carefully with necessary correction, both methods provide consistent results. Thus, AFM imaging can be used for DNA quantification, as an alternative to gel electrophoresis.

Mitochondrial DNA polymorphism A4917G is independently associated with age-related macular degeneration.

The objective of this study was to determine if MTND2*LHON4917G (4917G), a specific non-synonymous polymorphism in the mitochondrial genome previously associated with neurodegenerative phenotypes, is associated with increased risk for age-related macular degeneration (AMD). A preliminary study of 393 individuals (293 cases and 100 controls) ascertained at Vanderbilt revealed an increased occurrence of 4917G in cases compared to controls (15.4% vs.9.0%, p = 0.11). Since there was a significant age difference between cases and controls in this initial analysis, we extended the study by selecting Caucasian pairs matched at the exact age at examination. From the 1547 individuals in the Vanderbilt/Duke AMD population association study (including 157 in the preliminary study), we were able to match 560 (280 cases and 280 unaffected) on exact age at examination. This study population was genotyped for 4917G plus specific AMD-associated nuclear genome polymorphisms in CFH, LOC387715 and ApoE. Following adjustment for the listed nuclear genome polymorphisms, 4917G independently predicts the presence of AMD (OR = 2.16, 95%CI 1.20-3.91, p = 0.01). In conclusion, a specific mitochondrial polymorphism previously implicated in other neurodegenerative phenotypes (4917G) appears to convey risk for AMD independent of recently discovered nuclear DNA polymorphisms.

Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing.

We used ultra-deep sequencing to obtain tens of thousands of HIV-1 sequences from regions targeted by CD8+ T lymphocytes from longitudinal samples from three acutely infected subjects, and modeled viral evolution during the critical first weeks of infection. Previous studies suggested that a single virus established productive infection, but these conclusions were tempered because of limited sampling; now, we have greatly increased our confidence in this observation through modeling the observed earliest sample diversity based on vastly more extensive sampling. Conventional sequencing of HIV-1 from acute/early infection has shown different patterns of escape at different epitopes; we investigated the earliest escapes in exquisite detail. Over 3-6 weeks, ultradeep sequencing revealed that the virus explored an extraordinary array of potential escape routes in the process of evading the earliest CD8 T-lymphocyte responses--using 454 sequencing, we identified over 50 variant forms of each targeted epitope during early immune escape, while only 2-7 variants were detected in the same samples via conventional sequencing. In contrast to the diversity seen within epitopes, non-epitope regions, including the Envelope V3 region, which was sequenced as a control in each subject, displayed very low levels of variation. In early infection, in the regions sequenced, the consensus forms did not have a fitness advantage large enough to trigger reversion to consensus amino acids in the absence of immune pressure. In one subject, a genetic bottleneck was observed, with extensive diversity at the second time point narrowing to two dominant escape forms by the third time point, all within two months of infection. Traces of immune escape were observed in the earliest samples, suggesting that immune pressure is present and effective earlier than previously reported; quantifying the loss rate of the founder virus suggests a direct role for CD8 T-lymphocyte responses in viral containment after peak viremia. Dramatic shifts in the frequencies of epitope variants during the first weeks of infection revealed a complex interplay between viral fitness and immune escape.

Association between DNA damage response and repair genes and risk of invasive serous ovarian cancer.

BACKGROUND: We analyzed the association between 53 genes related to DNA repair and p53-mediated damage response and serous ovarian cancer risk using case-control data from the North Carolina Ovarian Cancer Study (NCOCS), a population-based, case-control study. METHODS/PRINCIPAL FINDINGS: The analysis was restricted to 364 invasive serous ovarian cancer cases and 761 controls of white, non-Hispanic race. Statistical analysis was two staged: a screen using marginal Bayes factors (BFs) for 484 SNPs and a modeling stage in which we calculated multivariate adjusted posterior probabilities of association for 77 SNPs that passed the screen. These probabilities were conditional on subject age at diagnosis/interview, batch, a DNA quality metric and genotypes of other SNPs and allowed for uncertainty in the genetic parameterizations of the SNPs and number of associated SNPs. Six SNPs had Bayes factors greater than 10 in favor of an association with invasive serous ovarian cancer. These included rs5762746 (median OR(odds ratio)(per allele) = 0.66; 95% credible interval (CI) = 0.44-1.00) and rs6005835 (median OR(per allele) = 0.69; 95% CI = 0.53-0.91) in CHEK2, rs2078486 (median OR(per allele) = 1.65; 95% CI = 1.21-2.25) and rs12951053 (median OR(per allele) = 1.65; 95% CI = 1.20-2.26) in TP53, rs411697 (median OR (rare homozygote) = 0.53; 95% CI = 0.35 - 0.79) in BACH1 and rs10131 (median OR( rare homozygote) = not estimable) in LIG4. The six most highly associated SNPs are either predicted to be functionally significant or are in LD with such a variant. The variants in TP53 were confirmed to be associated in a large follow-up study. CONCLUSIONS/SIGNIFICANCE: Based on our findings, further follow-up of the DNA repair and response pathways in a larger dataset is warranted to confirm these results.

Comparative analyses of seven algorithms for copy number variant identification from single nucleotide polymorphism arrays.

Determination of copy number variants (CNVs) inferred in genome wide single nucleotide polymorphism arrays has shown increasing utility in genetic variant disease associations. Several CNV detection methods are available, but differences in CNV call thresholds and characteristics exist. We evaluated the relative performance of seven methods: circular binary segmentation, CNVFinder, cnvPartition, gain and loss of DNA, Nexus algorithms, PennCNV and QuantiSNP. Tested data included real and simulated Illumina HumHap 550 data from the Singapore cohort study of the risk factors for Myopia (SCORM) and simulated data from Affymetrix 6.0 and platform-independent distributions. The normalized singleton ratio (NSR) is proposed as a metric for parameter optimization before enacting full analysis. We used 10 SCORM samples for optimizing parameter settings for each method and then evaluated method performance at optimal parameters using 100 SCORM samples. The statistical power, false positive rates, and receiver operating characteristic (ROC) curve residuals were evaluated by simulation studies. Optimal parameters, as determined by NSR and ROC curve residuals, were consistent across datasets. QuantiSNP outperformed other methods based on ROC curve residuals over most datasets. Nexus Rank and SNPRank have low specificity and high power. Nexus Rank calls oversized CNVs. PennCNV detects one of the fewest numbers of CNVs.

Computational Inference of Genome-Wide Protein-DNA Interactions Using High-Throughput Genomic Data

Transcriptional regulation has been studied intensively in recent decades. One important aspect of this regulation is the interaction between regulatory proteins, such as transcription factors (TF) and nucleosomes, and the genome. Different high-throughput techniques have been invented to map these interactions genome-wide, including ChIP-based methods (ChIP-chip, ChIP-seq, etc.), nuclease digestion methods (DNase-seq, MNase-seq, etc.), and others. However, a single experimental technique often only provides partial and noisy information about the whole picture of protein-DNA interactions. Therefore, the overarching goal of this dissertation is to provide computational developments for jointly modeling different experimental datasets to achieve a holistic inference on the protein-DNA interaction landscape.

We first present a computational framework that can incorporate the protein binding information in MNase-seq data into a thermodynamic model of protein-DNA interaction. We use a correlation-based objective function to model the MNase-seq data and a Markov chain Monte Carlo method to maximize the function. Our results show that the inferred protein-DNA interaction landscape is concordant with the MNase-seq data and provides a mechanistic explanation for the experimentally collected MNase-seq fragments. Our framework is flexible and can easily incorporate other data sources. To demonstrate this flexibility, we use prior distributions to integrate experimentally measured protein concentrations.

We also study the ability of DNase-seq data to position nucleosomes. Traditionally, DNase-seq has only been widely used to identify DNase hypersensitive sites, which tend to be open chromatin regulatory regions devoid of nucleosomes. We reveal for the first time that DNase-seq datasets also contain substantial information about nucleosome translational positioning, and that existing DNase-seq data can be used to infer nucleosome positions with high accuracy. We develop a Bayes-factor-based nucleosome scoring method to position nucleosomes using DNase-seq data. Our approach utilizes several effective strategies to extract nucleosome positioning signals from the noisy DNase-seq data, including jointly modeling data points across the nucleosome body and explicitly modeling the quadratic and oscillatory DNase I digestion pattern on nucleosomes. We show that our DNase-seq-based nucleosome map is highly consistent with previous high-resolution maps. We also show that the oscillatory DNase I digestion pattern is useful in revealing the nucleosome rotational context around TF binding sites.

Finally, we present a state-space model (SSM) for jointly modeling different kinds of genomic data to provide an accurate view of the protein-DNA interaction landscape. We also provide an efficient expectation-maximization algorithm to learn model parameters from data. We first show in simulation studies that the SSM can effectively recover underlying true protein binding configurations. We then apply the SSM to model real genomic data (both DNase-seq and MNase-seq data). Through incrementally increasing the types of genomic data in the SSM, we show that different data types can contribute complementary information for the inference of protein binding landscape and that the most accurate inference comes from modeling all available datasets.

This dissertation provides a foundation for future research by taking a step toward the genome-wide inference of protein-DNA interaction landscape through data integration.

Comparison of the thermodynamics and base-pair dynamics of a full LNA:DNA duplex and of the isosequential DNA:DNA duplex

Locked nucleic acids (LNA), conformationally restricted nucleotide analogues, are known to enhance pairing stability and selectivity toward complementary strands. With the aim to contribute to a better understanding of the origin of these effects, the structure, thermal stability, hybridization thermodynamics, and base-pair dynamics of a full-LNA:DNA heteroduplex and of its isosequential DNA:DNA homoduplex were monitored and compared. CD measurements highlight differences in the duplex structures: the homoduplex and heteroduplex present B-type and A-type helical conformations, respectively. The pairing of the hybrid duplex is characterized, at all temperatures monitored (between 15 and 37 degrees C), by a larger stability constant but a less favorable enthalpic term. A major contribution to this thermodynamic profile emanates from the presence of a hairpin structure in the LNA single strand which contributes favorably to the entropy of interaction but leads to an enthalpy penalty upon duplex formation. The base-pair opening dynamics of both systems was monitored by NMR spectroscopy via imino protons exchange measurements. The measurements highlight that hybrid G-C base-pairs present a longer base-pair lifetime and higher stability than natural G-C base-pairs, but that an LNA substitution in an A-T base-pair does not have a favorable effect on the stability. The thermodynamic and dynamic data confirm a more favorable stacking of the bases in the hybrid duplex. This study emphasizes the complementarities between dynamic and thermodynamical studies for the elucidation of the relevant factors in binding events.

Evidence for distinct dehydrogenase and isomerase sites within a single 3β-hydroxysteroid dehydrogenase/5-ene-4-ene isomerase protein

Complementary DNA encoding human 3β-hydroxysteroid dehydrogenase/5-ene-4-ene isomerase (30-HSD) has been expressed in transfected GH4C1 with use of the cytomegalovirus promoter. The activity of the expressed protein clearly shows that both dehydrogenase and isomerase enzymatic activities are present within a single protein. However, such findings do not indicate whether the two activities reside within one or two closely related catalytic sites. With use of [3H]-5-androstenedione, the intermediate compound in dehydroepiandrosterone (DHEA) transformation into 4-androstenedione by 3β-HSD, the present study shows that 4MA (N,N-diethyl-4-rnethyl-3-oxo-4-aza-5α-androstane-17β-carboxamide) and its analogues inhibit DHEA oxidation competitively while they exert a noncompetitive inhibition of the isomerization of 5-androstenedione to 4-androstenedione with an approximately 1000-fold higher Ki value. The present results thus strongly suggest that dehydrogenase and isomerase activities are present at separate sites on the 3β-HSD protein. In addition, using 5α-dihydrotestosterone (DHT) and 5α-androstane-3β,17β-diol as substrates for dehydrogenase activity only, we have found that dehydrogenase activity is reversibly and competitively inhibited by 4MA. Such data suggest that the irreversible step in the transformation of DHEA to 4-androstenedione is due to a separate site possessing isomerase activity that converts the 5-ene-3-keto to a much more stable 4-ene-3-keto configuration. © 1991 American Chemical Society.

The application of randomly amplified polymorphic DNA (RAPD's) to the study of Lasaea rubra (Montagu).

This paper describes the random amplification of polymorphic DNA markers (RAPDs) in Lasaea rubra (Erycinidae: Bivalvia). Present evidence suggests that L. rubra is an asexual species; however, the exact mode of clonal reproduction in this species is still a matter of debate. In this preliminary study, four of the primers used generated polymorphic RAPDs. One primer was able to distinguish between individuals from the same or different crevice population. This same primer also resolved a single band difference between otherwise identical RAPD patterns of a parent and its offspring. No familial differences have been detected in several previous studies using allozyme electrophoresis. This paper suggests that many polymorphic markers could be obtained with this species using the RAPD technique. Population genetic analysis of L. rubra has long been hampered by a dearth of polymorphic markers due to its small size. These findings suggest that this technique has the potential to further the study of population genetics in this asexual species.

Differences in mitochondrial and nuclear DNA status of high-density and low-density sperm fractions after density centrifugation preparation

Background: In order to isolate the â??bestâ?? sperm for assisted conception a discontinuous two-step density gradient centrifugation is usually employed. This technique is known to isolate a subpopulation with good motility, morphology and nuclear DNA (nDNA) integrity. As yet its ability to isolate sperm with unfragmented mitochondrial DNA (mtDNA) is unknown. Methods: Semen was obtained from men (n=28) attending our Regional Fertility Centre for infertility investigations. We employed a modified long polymerase chain reaction to study mtDNA and a modified alkaline Comet assay to determine nDNA fragmentation. Results: The high- density fraction displayed significantly more wild type mtDNA (75% of samples) than that of the low- density fraction (25% of samples). In the high-density fraction, there was a higher incidence of single, rather than double or multiple deletions and the deletions were predominantly small scale (0.1-4.0kb). There was a strong correlation between nDNA fragmentation, the number of mtDNA deletions (r=0.7, p

Retinal Pigment Epithelial Cell DNA is Damaged by Exposure to Benzo[a]pyrene, a Constituent of Cigarette Smoke.

This study examined the effect of exogenous benzo[ a ]pyrene (BaP), an important constituent of cigarette smoke, on cultured bovine retinal pigment epithelial (RPE) cells. Evidence is presented for its metabolic conversion into benzo[ a ]pyrene diol epoxide (BPDE) and the consequent formation of potentially cytotoxic nucleobase adducts in DNA. Cultured RPE cells were treated with BaP at concentrations in the range of 0–100 µm. The presence of BaP was found to cause inhibition of cell growth and replication. BaP induced the expression of a phase I drug metabolizing enzyme which was identified as cytochrome P450 1A1 (CYP 1A1) by RT–PCR and by Western blotting. Coincident with the increased expression of CYP 1A1, covalent adducts between the mutagenic metabolite BPDE and DNA could be detected within RPE cells by immunocytochemical staining. Additional support for their formation was afforded by nuclease P1 enhanced 32P-postlabelling assays on cellular DNA. Single-cell gel electrophoresis (comet) assays showed that exposure of RPE cells to BaP rendered them markedly more susceptible to DNA damage induced by broad band UVB or blue light laser irradiation. In the case of UVB, this is consistent with the photosensitization of DNA cleavage by nucleobase adducts of BPDE. Collectively, these findings imply that BaP has a significant impact on RPE cell pathophysiology and suggest mechanisms whereby exposure to cigarette smoke might cause RPE dysfunction and cell death, thus possibly contributing to degenerative disorders of the retina.

Cryopreservation of human semen and prepared sperm: effects on motility parameters and DNA integrity

Objective: To investigate effects of cryopreservation on sperm motility and DNA integrity. Design: Pre-cryopreservation and post-cryopreservation analysis of motility and DNA integrity of semen and prepared sperm samples. Setting: A hospital andrology laboratory. Patient(s): Forty men attending the Regional Fertility Centre, Belfast, Northern Ireland. Intervention(s): Each sample was divided, and an aliquot was frozen unprepared. Remaining aliquots were prepared by Percoll density centrifugation (95.0:47.5) or direct swim-up procedure and divided into aliquots to allow direct comparison of fresh and frozen semen and prepared sperm (frozen with or without the addition of seminal plasma) from the same ejaculate. Samples were frozen by static-phase vapor cooling and being plunged into liquid nitrogen. Thawing was carried out at room temperature. Main Outcome Measure(s): Sperm DNA integrity was determined using a modified alkaline single cell gel electrophoresis (comet) assay, and motility was determined using computer-assisted semen analysis. Result(s): Sperm frozen unprepared in seminal fluid appeared more resistant to freezing damage than frozen prepared sperm. Further improvements can be achieved by selecting out the subpopulation of sperm with best motility and DNA integrity and freezing these sperm in seminal plasma, making this the optimal procedure. Conclusion(s): Freezing sperm in seminal plasma improves postthaw motility and DNA integrity.

Assessment of DNA integrity and morphology of ejaculated sperm from fertile and infertile men before and after cryopreservation

Cryopreservation of human spermatozoa is extensively used in artifical insemination and IVF programmes. Despite various advances in cryopreservation methodology, the recovery rate of functional post thaw spermatozoa remains mediocre, with sperm motility being significantly decreased after freezing. The aim of this study was to investigate the effects of cryopreservation on both DNA integrity and morphology of spermatozoa from fertile and infertile men. Semen samples were obtained from 17 fertile men and 40 infertile men. All samples were prepared by discontinuous Percoll density centrifugation ( 95.0:47.5). Samples were divided into aliquots to allow direct comparison of fresh and frozen spermatozoa from the same ejaculate. Aliquots for cryopreservation were mixed with a commercial cryoprotectant and frozen by static phase vapour cooling before plunging into liquid nitrogen. Thawing was carried out slowly at room temperature. Sperm DNA integrity was determined using a modified alkaline single cell gel electrophoresis ( comet ) assay and sperm morphology analysed using the Tygerberg criteria. DNA of semen and prepared sperm from fertile men was found to be unaffected by cryopreservation. In marked contrast, spermatozoa from infertile men were significantly damaged by freeze- thawing. Cryopreservation had a detrimental effect on morphology of semen and prepared sperm from fertile and infertile men.

A comparison of mitochondrial and nuclear DNA status in testicular sperm from fertile men and those with obstructive azoospermia

Background: Mitochondria are vital to sperm as their motility powerhouses. They are also the only animal organelles with their own unique genome; encoding subunits for the complexes required for the electron transfer chain. Methods: A modified long PCR technique was used to study mitochondrial DNA (mtDNA) in ejaculated and testicular sperm samples from fertile men (n=11) and testicular sperm from men with obstructive azoospermia (n=25). Nuclear DNA fragmentation was measured by an alkaline single cell gel electrophoresis (COMET) assay. Results: Wild-type mtDNA was detected in only 60% of fertile mens�??�?�¢?? testicular sperm, 50% of their ejaculated sperm and 46% of testicular sperm from men with obstructive azoospermia. The incidence of mitochondrial deletions in testicular sperm of fertile and infertile men was not significantly different but the mean size of the deletions was significantly less in testicular sperm from fertile men compared with men with obstructive azoospermia (p<0.02). Nuclear DNA fragmentation in testicular sperm from fertile men and men with obstructive azoospermia was not significantly different. Conclusion: Multiple mtDNA deletions are common in testicular and ejaculated sperm from both fertile and infertile men. However, in males with obstructive azoospermia the mtDNA deletions in testicular sperm are of a larger scale.