Mutation screening of the medium-chain acyl-CoA dehydrogenase (MCAD) and the ornithine transcarbamylase (OTC) genes by multiplex PCR amplification and sequencing

BACKGROUND: Sequencing based mutation screening assays of genes encompassing large numbers of exons could be substantially optimized by multiplex PCR, which enables simultaneous amplification of many targets in one reaction. In the present study, a multiplex PCR protocol originally developed for fragment analysis was evaluated for sequencing based mutation screening of the ornithine transcarbamylase (OTC) and the medium-chain acyl-CoA dehydrogenase (MCAD) genes. METHODS: Single exon and multiplex PCR protocols were applied to generate PCR templates for subsequent DNA sequencing of all exons of the OTC and the MCAD genes. For each PCR protocol and using the same DNA samples, 66 OTC and 98 MCAD sequence reads were generated. The sequences derived from the two different PCR methods were compared at the level of individual signal-to-noise ratios of the four bases and the proportion of high-quality base-signals. RESULTS: The single exon and the multiplex PCR protocol gave qualitatively comparable results for the two genes. CONCLUSIONS: Many existing sequencing based mutation analysis protocols may be easily optimized with the proposed method, since the multiplex PCR protocol was successfully applied without any re-design of the PCR primers and other optimization steps for generating sequencing templates for the OTC and MCAD genes, respectively.

DNA-based diagnosis of rare diseases in veterinary medicine: a 4.4 kb deletion of ITGB4 is associated with epidermolysis bullosa in Charolais cattle.

BACKGROUND Rare diseases in livestock animals are traditionally poorly diagnosed. Other than clinical description and pathological examination, the underlying causes have, for the most part, remained unknown. A single case of congenital skin fragility in cattle was observed, necropsy, histological and ultrastructural examinations were carried out and whole genome sequencing was utilized to identify the causative mutation. RESULTS A single purebred female Charolais calf with severe skin lesions was delivered full-term and died spontaneously after birth. The clinical and pathological findings exactly matched the gross description given by previous reports on epitheliogenesis imperfecta and epidermolysis bullosa (EB) in cattle. Histological and ultrastructural changes were consistent with EB junctionalis (EBJ). Genetic analysis revealed a previously unpublished ITGB4 loss-of-function mutation; the affected calf was homozygous for a 4.4 kb deletion involving exons 17 to 22, and the dam carried a single copy of the deletion indicating recessive inheritance. The homozygous mutant genotype did not occur in healthy controls of various breeds but some heterozygous carriers were found among Charolais cattle belonging to the affected herd. The mutant allele was absent in a representative sample of unrelated sires of the German Charolais population. CONCLUSION This is the first time in which a recessively inherited ITGB4 associated EBJ has been reported in cattle. The identification of heterozygous carriers is of importance in avoiding the transmission of this defect in future. Current DNA sequencing methods offer a powerful tool for understanding the genetic background of rare diseases in domestic animals having a reference genome sequence available.

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Direct sequencing of bacterial and P1 artificial chromosome-nested deletions for identifying position-specific single-nucleotide polymorphisms

A loxP-transposon retrofitting strategy for generating large nested deletions from one end of the insert DNA in bacterial artificial chromosomes and P1 artificial chromosomes was described recently [Chatterjee, P. K. & Coren, J. S. (1997) Nucleic Acids Res. 25, 2205–2212]. In this report, we combine this procedure with direct sequencing of nested-deletion templates by using primers located in the transposon end to illustrate its value for position-specific single-nucleotide polymorphism (SNP) discovery from chosen regions of large insert clones. A simple ampicillin sensitivity screen was developed to facilitate identification and recovery of deletion clones free of transduced transposon plasmid. This directed approach requires minimal DNA sequencing, and no in vitro subclone library generation; positionally oriented SNPs are a consequence of the method. The procedure is used to discover new SNPs as well as physically map those identified from random subcloned libraries or sequence databases. The deletion templates, positioned SNPs, and markers are also used to orient large insert clones into a contig. The deletion clone can serve as a ready resource for future functional genomic studies because each carries a mammalian cell-specific antibiotic resistance gene from the transposon. Furthermore, the technique should be especially applicable to the analysis of genomes for which a full genome sequence or radiation hybrid cell lines are unavailable.

Fluorescence-based directed termination PCR: direct mutation characterization without sequencing

We describe a fluorescence-based directed termination PCR (fluorescent DT–PCR) that allows accurate determination of actual sequence changes without dideoxy DNA sequencing. This is achieved using near infrared dye-labeled primers and performing two PCR reactions under low and unbalanced dNTP concentrations. Visualization of resulting termination fragments is accomplished with a dual dye Li-cor DNA sequencer. As each DT–PCR reaction generates two sets of terminating fragments, a pair of complementary reactions with limiting dATP and dCTP collectively provide information on the entire sequence of a target DNA, allowing an accurate determination of any base change. Blind analysis of 78 mutants of the supF reporter gene using fluorescent DT–PCR not only correctly determined the nature and position of all types of substitution mutations in the supF gene, but also allowed rapid scanning of the signature sequences among identical mutations. The method provides simplicity in the generation of terminating fragments and 100% accuracy in mutation characterization. Fluorescent DT–PCR was successfully used to generate a UV-induced spectrum of mutations in the supF gene following replication on a single plate of human DNA repair-deficient cells. We anticipate that the automated DT–PCR method will serve as a cost-effective alternative to dideoxy sequencing in studies involving large-scale analysis for nucleotide sequence changes.

Cloning and expression in Escherichia coli of the OGG1 gene of Saccharomyces cerevisiae, which codes for a DNA glycosylase that excises 7,8-dihydro-8-oxoguanine and 2,6-diamino-4-hydroxy-5-N-methylformamidopyrimidine.

A spontaneous mutator strain of Escherichia coli (fpg mutY) was used to clone the OGG1 gene of Saccharomyces cerevisiae, which encodes a DNA glycosylase activity that excises 7,8-dihydro-8-oxoguanine (8-OxoG). E. coli (fpg mutY) was transformed by a yeast DNA library, and clones that showed a reduced spontaneous mutagenesis were selected. The antimutator activity was associated with pYSB10, an 11-kbp recombinant plasmid. Cell-free extracts of E. coli (fpg mutY) harboring pYSB10 possess an enzymatic activity that cleaves a 34-mer oligonucleotide containing a single 8-oxoG opposite a cytosine (8-OxoG/C). The yeast DNA fragment of 1.7 kbp that suppresses spontaneous mutagenesis and overproduces the 8-OxoG/C cleavage activity was sequenced and mapped to chromosome XIII. DNA sequencing identified an open reading frame, designated OGG1, which encodes a protein of 376 amino acids with a molecular mass of 43 kDa. The OGG1 gene was inserted in plasmid pUC19, yielding pYSB110. E. coli (fpg) harboring pYSB110 was used to purify the Ogg1 protein of S. cerevisiae to apparent homogeneity. The Ogg1 protein possesses a DNA glycosylase activity that releases 8-OxoG and 2,6-diamino-4-hydroxy-5-N-methylformamidopyrimidine. The Ogg1 protein preferentially incises DNA that contains 8-OxoG opposite cytosine (8-OxoG/C) or thymine (8-OxoG/T). In contrast, Ogg1 protein does not incise the duplex where an adenine is placed opposite 8-OxoG (8-OxoG/A). The mechanism of strand cleavage by Ogg1 protein is probably due to the excision of 8-OxoG followed by a beta-elimination at the resulting apurinic/apyrimidinic site.

Transposon tools for recombinant DNA manipulation: characterization of transcriptional regulators from yeast, Xenopus, and mouse.

Transposon Tn1000 has been adapted to deliver novel DNA sequences for manipulating recombinant DNA. The transposition procedure for these "tagged" Tn1000s is simple and applicable to most plasmids in current use. For yeast molecular biology, tagged Tn1000s introduce a variety of yeast selective markers and replication origins into plasmids and cosmids. In addition, the beta-globin minimal promoter and lacZ gene of Tn(beta)lac serve as a mobile reporter of eukaryotic enhancer activity. In this paper, Tn(beta)lac was used to localize a mouse HoxB-complex enhancer in transgenic mice. Other tagged transposons create Gal4 DNA-binding-domain fusions, in either Escherichia coli or yeast plasmids, for use in one- and two-hybrid tests of transcriptional activation and protein-protein interaction, respectively. With such fusions, the Saccharomyces cerevisiae Swi6 G1/S-phase transcription factor and the Xenopus laevis Pintallavis developmental regulator are shown to activate transcription. Furthermore, the same transposon insertions also facilitated mapping of the Swi6 and Pintallavis domains responsible for transcriptional activation. Thus, as well as introducing novel sequences, tagged transposons share the numerous other applications of transposition such as producing insertional mutations, creating deletion series, or serving as mobile primer sites for DNA sequencing.

Identification of a DNA transformation gene required for com101A+ expression and supertransformer phenotype in Haemophilus influenzae.

DNA sequencing, RNA mapping, and protein expression experiments revealed the presence of a gene, tfoX+, encoding a 24.9-kDa polypeptide, that is transcribed divergently from a common promoter region with the Haemophilus influenzae rec-1+ gene. H. influenzae strains mutant for tfoX failed to bind transforming DNA and were transformation deficient. Primer extension experiments utilizing in vivo total RNA from precompetent and competent H. influenzae cells demonstrated that transcription of tfoX+ increased immediately upon competence induction, suggesting that tfoX+ is an early competence gene. Similar experiments showed that the expression of the late competence-specific gene, com101A+, was tfoX+ dependent. Moreover, expression of plasmid-borne tfoX+ in H. influenzae resulted in constitutive competence. The addition of cyclic adenosine monophosphate (cAMP) to strains carrying a tfoX::lacZ operon fusion resulted in an immediate increase in beta-galactosidase activity that correlated with an increase in genetic transformability. Collectively, our results suggest that TfoX may play a key role in the development of genetic competence by regulating the expression of late competence-specific genes.

Learning Models of Gene Expression from Synthetic DNA Sequences

Thesis (Ph.D.)--University of Washington, 2016-06

Understanding next generation sequencing in oncology: A guide for oncologists

DNA sequencing is now faster and cheaper than ever before, due to the development of next generation sequencing (NGS) technologies. NGS is now widely used in the research setting and is becoming increasingly utilised in clinical practice. However, due to evolving clinical commitments, increased workload and lack of training opportunities, many oncologists may be unfamiliar with the terminology and technology involved. This can lead to oncologists feeling daunted by issues such as how to interpret the vast amounts of data generated by NGS and the differences between sequencing platforms. This review article explains common concepts and terminology, summarises the process of DNA sequencing (including data analysis) and discusses the main factors to consider when deciding on a sequencing method. This article aims to improve oncologists' understanding of the most commonly used sequencing platforms and the ongoing challenges faced in expanding the use of NGS into routine clinical practice.

Tecnología Nanopore para la secuenciación de DNA

Nanopore-based sequencer will open the path to the fourth-generation DNA sequencing technology. The main differences between this technique and the previous ones are: DNA molecule that will be sequenced does not need a previous amplification step, is not necessary any type of specific label both molecular adaptors, and it has been abolished enzymatic process in the nucleotide sequence identification event. These differences have as result a more economic method since don’t spend the necessary reagents for the previous techniques, furthermore it lets to sequence samples with a low DNA concentration. This technique is based in the use of a membrane with a biologic nanopore inserted in it whereby the molecule to analyze (analyte) it made to pass, this membrane is placed between two reservoirs containing ions, when an external volatage is applied in both sides this lead to an ion current through the nanopore. When an analyte cross the nanopore the ion current is modified, that modification in the amplitude and duration of ion current determine the physical and chemical properties of that analyte. By means of subsequent statistical analyzes it can be determined to what sequence own this ion current blockade patterns. More used nanopores are the biologic ones, although they are working to develop synthetic nanopores. The main biologic nanopores are: α-Hemolysin from Staphylococcus aureus (α-HL), Mycobacterium smegmatis porin A (MspA) and bacteriophage phi29 pore (phi29). Α-HL and MspA have in their narrowest point a diameter similar to nucleotide size, they are functional at high temperature both wide range of pH (2-12) but MspA is able to read four nucleotide at the same time while α- HL just can read one by one. Finally, phi29 present a bigger diameter what let to get information about DNA spatial conformation and their interaction with proteins (Feng et al., 2015). Nowaday Oxford Nanopore Technologies (ONT) is the only company which has developed Nanopore technology; they have two devices available to sequencing (PromethION and MinION). The MinION is a single-use DNA sequencing device with the size of a USB memory with a total of 3000 nanopores that can sequence until 200kb. The PrometheION is big size sequencer that own 48 different cells, what let to sequence different samples at the same time, with a total of 144.000 nanopores and reading of several megabases (https://www.nanoporetech.com/). The high processivity and low cost become this technique in a great option to massive- sequencing.

Genetic Structure of the Florida Key Tree Cactus, Pilosocereus robinii, using Restriction Site associated DNA (RAD) markers

Rare plant conservation efforts must utilize current genetic methods to ensure the evolutionary potential of populations is preserved. One such effort involves the Key Tree Cactus, Pilosocereus robinii, which is an endangered columnar cactus native to the Florida Keys. The populations have precipitously declined over the past decade because of habitat loss and increasing soil salinity from rising sea levels and storm surge. Next-generation DNA sequencing was used to assess the genetic structure of the populations. Twenty individuals representative of both wild and extirpated cacti were chosen for Restriction Site Associated DNA (RAD) analysis. Samples processed using the HindIII and NotIII restriction enzymes produced 82,382,440 high quality reads used for genetic mapping, from which 5,265 Single Nucleotide Polymorphisms (SNPs) were discovered. The analysis revealed that the Keys’ populations are closely related with little population differentiation. In addition, the populations display evidence of inbreeding and low genetic diversity.

The development of 10 novel polymorphic microsatellite markers through next generation sequencing and a preliminary population genetic analysis for the endangered Glenelg spiny crayfish, Euastacus bispinosus

The Glenelg spiny crayfish, Euastacus bispinosus, is an iconic freshwater invertebrate of south eastern Australia and listed as 'endangered' under the Environment Protection and Biodiversity Conservation Act 1999, and 'vulnerable' under the International Union for Conservation of Nature's Red List. The species has suffered major population declines as a result of over-fishing, low environmental flows, the introduction of invasive fish species and habitat degradation. In order to develop an effective conservation strategy, patterns of gene flow, genetic structure and genetic diversity across the species distribution need to be clearly understood. In this study we develop a suite of polymorphic microsatellite markers by next generation sequencing. A total of 15 polymorphic loci were identified and 10 characterized using 22 individuals from the lower Glenelg River. We observed low to moderate genetic variation across most loci (mean number of alleles per locus = 2.80; mean expected heterozygosity = 0.36) with no evidence of individual loci deviating significantly from Hardy-Weinberg equilibrium. Marker independence was confirmed with tests for linkage disequilibrium, and analyses indicated no evidence of null alleles across loci. Individuals from two additional sites (Crawford River, Victoria; Ewens Ponds Conservation Park, South Australia) were genotyped at all 10 loci and a preliminary investigation of genetic diversity and population structure was undertaken. Analyses indicate high levels of genetic differentiation among sample locations (F ST = 0.49), while the Ewens Ponds population is genetically homogeneous, indicating a likely small founder group and ongoing inbreeding. Management actions will be needed to restore genetic diversity in this and possibly other at risk populations. These markers will provide a valuable resource for future population genetic assessments so that an effective framework can be developed for implementing conservation strategies for E. bispinosus.

A Kallikrein 15 (KLK15) single nucleotide polymorphism located close to a novel exon shows evidence of association with poor ovarian cancer survival

KLK15 over-expression is reported to be a significant predictor of reduced progression-free survival and overall survival in ovarian cancer. Our aim was to analyse the KLK15 gene for putative functional single nucleotide polymorphisms (SNPs) and assess the association of these and KLK15 HapMap tag SNPs with ovarian cancer survival. Results In silico analysis was performed to identify KLK15 regulatory elements and to classify potentially functional SNPs in these regions. After SNP validation and identification by DNA sequencing of ovarian cancer cell lines and aggressive ovarian cancer patients, 9 SNPs were shortlisted and genotyped using the Sequenom iPLEX Mass Array platform in a cohort of Australian ovarian cancer patients (N = 319). In the Australian dataset we observed significantly worse survival for the KLK15 rs266851 SNP in a dominant model (Hazard Ratio (HR) 1.42, 95% CI 1.02-1.96). This association was observed in the same direction in two independent datasets, with a combined HR for the three studies of 1.16 (1.00-1.34). This SNP lies 15bp downstream of a novel exon and is predicted to be involved in mRNA splicing. The mutant allele is also predicted to abrogate an HSF-2 binding site. Conclusions We provide evidence of association for the SNP rs266851 with ovarian cancer survival. Our results provide the impetus for downstream functional assays and additional independent validation studies to assess the role of KLK15 regulatory SNPs and KLK15 isoforms with alternative intracellular functional roles in ovarian cancer survival.