Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

New Contig Creation Algorithm for the de novo DNA Assembly Problem

DNA assembly is among the most fundamental and difficult problems in bioinformatics. Near optimal assembly solutions are available for bacterial and small genomes, however assembling large and complex genomes especially the human genome using Next-Generation-Sequencing (NGS) technologies is shown to be very difficult because of the highly repetitive and complex nature of the human genome, short read lengths, uneven data coverage and tools that are not specifically built for human genomes. Moreover, many algorithms are not even scalable to human genome datasets containing hundreds of millions of short reads. The DNA assembly problem is usually divided into several subproblems including DNA data error detection and correction, contig creation, scaffolding and contigs orientation; each can be seen as a distinct research area. This thesis specifically focuses on creating contigs from the short reads and combining them with outputs from other tools in order to obtain better results. Three different assemblers including SOAPdenovo [Li09], Velvet [ZB08] and Meraculous [CHS+11] are selected for comparative purposes in this thesis. Obtained results show that this thesis’ work produces comparable results to other assemblers and combining our contigs to outputs from other tools, produces the best results outperforming all other investigated assemblers.

Estimate of corrected values and the effect of correction for measurement error in dietary data obtained by the Food Frequency Questionnaire for Adolescents (AFFQ)

The scope of this study was to estimate calibrated values for dietary data obtained by the Food Frequency Questionnaire for Adolescents (FFQA) and illustrate the effect of this approach on food consumption data. The adolescents were assessed on two occasions, with an average interval of twelve months. In 2004, 393 adolescents participated, and 289 were then reassessed in 2005. Dietary data obtained by the FFQA were calibrated using the regression coefficients estimated from the average of two 24-hour recalls (24HR) of the subsample. The calibrated values were similar to the the 24HR reference measurement in the subsample. In 2004 and 2005 a significant difference was observed between the average consumption levels of the FFQA before and after calibration for all nutrients. With the use of calibrated data the proportion of schoolchildren who had fiber intake below the recommended level increased. Therefore, it is seen that calibrated data can be used to obtain adjusted associations due to reclassification of subjects within the predetermined categories.

Functional correction in mouse models of muscular dystrophy using exon-skipping tricyclo-DNA oligomers

Antisense oligonucleotides (AONs) hold promise for therapeutic correction of many genetic diseases via exon skipping, and the first AON-based drugs have entered clinical trials for neuromuscular disorders1, 2. However, despite advances in AON chemistry and design, systemic use of AONs is limited because of poor tissue uptake, and recent clinical reports confirm that sufficient therapeutic efficacy has not yet been achieved. Here we present a new class of AONs made of tricyclo-DNA (tcDNA), which displays unique pharmacological properties and unprecedented uptake by many tissues after systemic administration. We demonstrate these properties in two mouse models of Duchenne muscular dystrophy (DMD), a neurogenetic disease typically caused by frame-shifting deletions or nonsense mutations in the gene encoding dystrophin3, 4 and characterized by progressive muscle weakness, cardiomyopathy, respiratory failure5 and neurocognitive impairment6. Although current naked AONs do not enter the heart or cross the blood-brain barrier to any substantial extent, we show that systemic delivery of tcDNA-AONs promotes a high degree of rescue of dystrophin expression in skeletal muscles, the heart and, to a lesser extent, the brain. Our results demonstrate for the first time a physiological improvement of cardio-respiratory functions and a correction of behavioral features in DMD model mice. This makes tcDNA-AON chemistry particularly attractive as a potential future therapy for patients with DMD and other neuromuscular disorders or with other diseases that are eligible for exon-skipping approaches requiring whole-body treatment.

The DNA Helicase Activity of BLM Is Necessary for the Correction of the Genomic Instability of Bloom Syndrome Cells

Bloom syndrome (BS) is a rare autosomal recessive disorder characterized by growth deficiency, immunodeficiency, genomic instability, and the early development of cancers of many types. BLM, the protein encoded by BLM, the gene mutated in BS, is localized in nuclear foci and absent from BS cells. BLM encodes a DNA helicase, and proteins from three missense alleles lack displacement activity. BLM transfected into BS cells reduces the frequency of sister chromatid exchanges and restores BLM in the nucleus. Missense alleles fail to reduce the sister chromatid exchanges in transfected BS cells or restore the normal nuclear pattern. BLM complements a phenotype of a Saccharomyces cerevisiae sgs1 top3 strain, and the missense alleles do not. This work demonstrates the importance of the enzymatic activity of BLM for its function and nuclear localization pattern.

Roles of DNA polymerases V and II in SOS-induced error-prone and error-free repair in Escherichia coli

DNA polymerase V, composed of a heterotrimer of the DNA damage-inducible UmuC and UmuD\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{2}^{^{\prime}}}}\end{equation*}\end{document} proteins, working in conjunction with RecA, single-stranded DNA (ssDNA)-binding protein (SSB), β sliding clamp, and γ clamp loading complex, are responsible for most SOS lesion-targeted mutations in Escherichia coli, by catalyzing translesion synthesis (TLS). DNA polymerase II, the product of the damage-inducible polB (dinA ) gene plays a pivotal role in replication-restart, a process that bypasses DNA damage in an error-free manner. Replication-restart takes place almost immediately after the DNA is damaged (≈2 min post-UV irradiation), whereas TLS occurs after pol V is induced ≈50 min later. We discuss recent data for pol V-catalyzed TLS and pol II-catalyzed replication-restart. Specific roles during TLS for pol V and each of its accessory factors have been recently determined. Although the precise molecular mechanism of pol II-dependent replication-restart remains to be elucidated, it has recently been shown to operate in conjunction with RecFOR and PriA proteins.

Targeted gene correction of episomal DNA in mammalian cells mediated by a chimeric RNA.DNA oligonucleotide.

An experimental strategy to facilitate correction of single-base mutations of episomal targets in mammalian cells has been developed. The method utilizes a chimeric oligonucleotide composed of a contiguous stretch of RNA and DNA residues in a duplex conformation with double hairpin caps on the ends. The RNA/DNA sequence is designed to align with the sequence of the mutant locus and to contain the desired nucleotide change. Activity of the chimeric molecule in targeted correction was tested in a model system in which the aim was to correct a point mutation in the gene encoding the human liver/bone/kidney alkaline phosphatase. When the chimeric molecule was introduced into cells containing the mutant gene on an extrachromosomal plasmid, correction of the point mutation was accomplished with a frequency approaching 30%. These results extend the usefulness of the oligonucleotide-based gene targeting approaches by increasing specific targeting frequency. This strategy should enable the design of antiviral agents.

Correction of chromosomal instability and sensitivity to diverse mutagens by a cloned cDNA of the XRCC3 DNA repair gene.

The mutagen-sensitive CHO line irs1SF was previously isolated on the basis of hypersensitivity to ionizing radiation and was found to be chromosomally unstable as well as cross-sensitive to diverse kinds of DNA-damaging agents. The analysis of somatic cell hybrids formed between irs1SF and human lymphocytes implicated a human gene (defined as XRCC3; x-ray repair cross-complementing), which partially restored mitomycin C resistance to the mutant. A functional cDNA that confers mitomycin C resistance was transferred to irs1SF cells by transforming them with an expression cDNA library and obtaining primary and secondary transformants. Functional cDNA clones were recovered from a cosmid library prepared from a secondary transformant. Transformants also showed partial correction of sensitivity to cisplatin and gamma-rays, efficient correction of chromosomal instability, and substantially improved plating efficiency and growth rate. The XRCC3 cDNA insert is approximately 2.5 kb and detects an approximately 3.0-kb mRNA on Northern blots. The cDNA was mapped by fluorescence in situ hybridization to human chromosome 14q32.3, which was consistent with the chromosome concordance data of two independent hybrid clone panels.

New Approach for Correction of Error Associated With Keratometric Estimation of Corneal Power in Keratoconus

The aim of this study was to obtain the exact value of the keratometric index (nkexact) and to clinically validate a variable keratometric index (nkadj) that minimizes this error. Methods: The nkexact value was determined by obtaining differences (DPc) between keratometric corneal power (Pk) and Gaussian corneal power (PGauss c ) equal to 0. The nkexact was defined as the value associated with an equivalent difference in the magnitude of DPc for extreme values of posterior corneal radius (r2c) for each anterior corneal radius value (r1c). This nkadj was considered for the calculation of the adjusted corneal power (Pkadj). Values of r1c ∈ (4.2, 8.5) mm and r2c ∈ (3.1, 8.2) mm were considered. Differences of True Net Power with PGauss c , Pkadj, and Pk(1.3375) were calculated in a clinical sample of 44 eyes with keratoconus. Results: nkexact ranged from 1.3153 to 1.3396 and nkadj from 1.3190 to 1.3339 depending on the eye model analyzed. All the nkadj values adjusted perfectly to 8 linear algorithms. Differences between Pkadj and PGauss c did not exceed 60.7 D (Diopter). Clinically, nk = 1.3375 was not valid in any case. Pkadj and True Net Power and Pk(1.3375) and Pkadj were statistically different (P , 0.01), whereas no differences were found between PGauss c and Pkadj (P . 0.01). Conclusions: The use of a single value of nk for the calculation of the total corneal power in keratoconus has been shown to be imprecise, leading to inaccuracies in the detection and classification of this corneal condition. Furthermore, our study shows the relevance of corneal thickness in corneal power calculations in keratoconus.

Error detection and correction in binary parallel digital computers /

"Retyped October, 1964"

Correction of the misclassification of species in the Portuguese collection of Cucurbita pepo L. using DNA markers

In this study, the genetic variability among 130 accessions of the Portuguese germplasm collection of Cucurbita pepo L. maintained at the Banco Portugues de Germoplasma Vegetal was assessed using AFLP (amplified fragment length polymorphism) and RAPD (random amplified polymorphic DNA) techniques for the identification of a genetically diverse core group of accessions for field phenotypic analysis. The surprisingly completely different molecular patterns exhibited by multiple accessions was later confirmed in the distribution of the putative C. pepo plants into two clusters drastically separated at a very low level of genetic similarity (DICE coefficient = 0.37). Additional analyses with RAPD and ISSR (inter single sequence repeat) markers and the introduction of standard genotypes of C. maxima L. and C. moschata L. into the analyses allowed the identification of multiple accessions of the last species wrongly included in the C. pepo collection. This study is a good example of the usefulness of DNA markers in the establishment and management of plant germplasm collections.

Correction of the misclassification of species in the Portuguese collection of Cucurbita pepo L. using DNA markers

In this study, the genetic variability among 130 accessions of the Portuguese germplasm collection of Cucurbita pepo L. maintained at the Banco Portugues de Germoplasma Vegetal was assessed using AFLP (amplified fragment length polymorphism) and RAPD (random amplified polymorphic DNA) techniques for the identification of a genetically diverse core group of accessions for field phenotypic analysis. The surprisingly completely different molecular patterns exhibited by multiple accessions was later confirmed in the distribution of the putative C. pepo plants into two clusters drastically separated at a very low level of genetic similarity (DICE coefficient = 0.37). Additional analyses with RAPD and ISSR (inter single sequence repeat) markers and the introduction of standard genotypes of C. maxima L. and C. moschata L. into the analyses allowed the identification of multiple accessions of the last species wrongly included in the C. pepo collection. This study is a good example of the usefulness of DNA markers in the establishment and management of plant germplasm collections.

Patient-specific computational biomechanics for simulating adolescent scoliosis surgery : predicted vs clinical correction for a preliminary series of six patients

Scoliosis is a spinal deformity that requires surgical correction in progressive cases. In order to optimize surgical outcomes, patient-specific finite element models are being developed by our group. In this paper, a single rod anterior correction procedure is simulated for a group of six scoliosis patients. For each patient, personalised model geometry was derived from low-dose CT scans, and clinically measured intra-operative corrective forces were applied. However, tissue material properties were not patient-specific, being derived from existing literature. Clinically, the patient group had a mean initial Cobb angle of 47.3 degrees, which was corrected to 17.5 degrees after surgery. The mean simulated post-operative Cobb angle for the group was 18.1 degrees. Although this represents good agreement between clinical and simulated corrections, the discrepancy between clinical and simulated Cobb angle for individual patients varied between -10.3 and +8.6 degrees, with only three of the six patients matching the clinical result to within accepted Cobb measurement error of +-5 degrees. The results of this study suggest that spinal tissue material properties play an important role in governing the correction obtained during surgery, and that patient-specific modelling approaches must address the question of how to prescribe patient-specific soft tissue properties for spine surgery simulation.

The measurement of applied forces during anterior single rod surgical correction of adolescent idiopathic scoliosis

INTRODUCTION. Following anterior thoracoscopic instrumentation and fusion for the treatment of thoracic AIS, implant related complications have been reported as high as 20.8%. Currently the magnitudes of the forces applied to the spine during anterior scoliosis surgery are unknown. The aim of this study was to measure the segmental compressive forces applied during anterior single rod instrumentation in a series of adolescent idiopathic scoliosis patients. METHODS. A force transducer was designed, constructed and retrofitted to a surgical cable compression tool, routinely used to apply segmental compression during anterior scoliosis correction. Transducer output was continuously logged during the compression of each spinal joint, the output at completion converted to an applied compression force using calibration data. The angle between adjacent vertebral body screws was also measured on intra-operative frontal plane fluoroscope images taken both before and after each joint compression. The difference in angle between the two images was calculated as an estimate for the achieved correction at each spinal joint. RESULTS. Force measurements were obtained for 15 scoliosis patients (Aged 11-19 years) with single thoracic curves (Cobb angles 47˚- 67˚). In total, 95 spinal joints were instrumented. The average force applied for a single joint was 540 N (± 229 N)ranging between 88 N and 1018 N. Experimental error in the force measurement, determined from transducer calibration was ± 43 N. A trend for higher forces applied at joints close to the apex of the scoliosis was observed. The average joint correction angle measured by fluoroscope imaging was 4.8˚ (±2.6˚, range 0˚-12.6˚). CONCLUSION. This study has quantified in-vivo, the intra-operative correction forces applied by the surgeon during anterior single rod instrumentation. This data provides a useful contribution towards an improved understanding of the biomechanics of scoliosis correction. In particular, this data will be used as input for developing patient-specific finite element simulations of scoliosis correction surgery.