Quantification of single nucleotide polymorphisms by automated DNA sequencing

ABSP: an automated R tool to efficiently analyze region-specific CpG methylation from bisulfite sequencing PCR

MOTIVATION

Nowadays, epigenetic gene regulations are studied in each part of the biology, from embryonic development to diseases such as cancers and neurodegenerative disorders. Currently, to quantify and compare CpG methylation levels of a specific region of interest, the most accessible technique is the bisulfite sequencing PCR (BSP). However, no existing user-friendly tool is able to analyze data from all approaches of BSP. Therefore, the most convenient way to process results from the direct sequencing of PCR products (direct-BSP) is to manually analyze the chromatogram traces, which is a repetitive and prone to error task.

RESULTS

Here, we implement a new R-based tool, called ABSP for analysis of bisulfite sequencing PCR, providing a complete analytic process of both direct-BSP and cloning-BSP data. It uses the raw sequencing trace files (.ab1) as input to compute and compare CpG methylation percentages. It is fully automated and includes a user-friendly interface as a built-in R shiny app, quality control steps and generates publication-ready graphics.

AVAILABILITY AND IMPLEMENTATION

The ABSP tool and associated data are available on GitHub at https://github.com/ABSP-methylation-tool/ABSP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DNA Methylation in Solid Tumors: Functions and Methods of Detection

DNA methylation, i.e., addition of methyl group to 5′-carbon of cytosine residues in CpG dinucleotides, is an important epigenetic modification regulating gene expression, and thus implied in many cellular processes. Deregulation of DNA methylation is strongly associated with onset of various diseases, including cancer. Here, we review how DNA methylation affects carcinogenesis process and give examples of solid tumors where aberrant DNA methylation is often present. We explain principles of methods developed for DNA methylation analysis at both single gene and whole genome level, based on (i) sodium bisulfite conversion, (ii) methylation-sensitive restriction enzymes, and (iii) interactions of 5-methylcytosine (5mC) with methyl-binding proteins or antibodies against 5mC. In addition to standard methods, we describe recent advances in next generation sequencing technologies applied to DNA methylation analysis, as well as in development of biosensors that represent their cheaper and faster alternatives. Most importantly, we highlight not only advantages, but also disadvantages and challenges of each method.

DNA methylation: a timeline of methods and applications

DNA methylation is a biochemical process where a DNA base, usually cytosine, is enzymatically methylated at the 5-carbon position. An epigenetic modification associated with gene regulation, DNA methylation is of paramount importance to biological health and disease. Recently, the quest to unravel the Human Epigenome commenced, calling for a modernization of previous DNA methylation profiling techniques. Here, we describe the major developments in the methodologies used over the past three decades to examine the elusive epigenome (or methylome). The earliest techniques were based on the separation of methylated and unmethylated cytosines via chromatography. The following years would see molecular techniques being employed to indirectly examine DNA methylation levels at both a genome-wide and locus-specific context, notably immunoprecipitation via anti-5'methylcytosine and selective digestion with methylation-sensitive restriction endonucleases. With the advent of sodium bisulfite treatment of DNA, a deamination reaction that converts cytosine to uracil only when unmethylated, the epigenetic modification can now be identified in the same manner as a DNA base-pair change. More recently, these three techniques have been applied to more technically advanced systems such as DNA microarrays and next-generation sequencing platforms, bringing us closer to unveiling a complete human epigenetic profile.

Functional analysis of upstream common polymorphisms of the dopamine transporter gene

The human dopamine transporter (DAT, SLC6A3) has been extensively investigated because of its potential involvement in neuropsychiatric disorders. The core elements responsible for its transcription have been identified. A regulatory role for certain genomic variants upstream to the core promoter is known. Recently, other single-nucleotide polymorphisms (SNPs) have been identified in this region and are thought to be associated with schizophrenia and bipolar I disorder. Hence, we have investigated the impact of common SNPs in a 2.8-kilobase region flanking the core promoter region (-2.7 to +63 base pair) in the neuroblastoma cell line SH-SY5Y. Haplotypes generated by site-directed mutagenesis revealed varying impact of individual SNPs on promoter activity using dual luciferase assays. In silico analyses also predicted allele-specific binding of transcription factors for some of these SNPs. Though electrophoretic mobility shift assays indicated several factors that appeared to bind to specific sites within this region, allele-specific binding was not detected for any SNP apart from rs3756450. We have thus identified novel putative regulatory domains flanking the core promoter of DAT that merit further investigation.

Successful computational prediction of novel imprinted genes from epigenomic features

Approximately 100 mouse genes undergo genomic imprinting, whereby one of the two parental alleles is epigenetically silenced. Imprinted genes influence processes including development, X chromosome inactivation, obesity, schizophrenia, and diabetes, motivating the identification of all imprinted loci. Local sequence features have been used to predict candidate imprinted genes, but rigorous testing using reciprocal crosses validated only three, one of which resided in previously identified imprinting clusters. Here we show that specific epigenetic features in mouse cells correlate with imprinting status in mice, and we identify hundreds of additional genes predicted to be imprinted in the mouse. We used a multitiered approach to validate imprinted expression, including use of a custom single nucleotide polymorphism array and traditional molecular methods. Of 65 candidates subjected to molecular assays for allele-specific expression, we found 10 novel imprinted genes that were maternally expressed in the placenta.

Rapid quantification of DNA methylation by measuring relative peak heights in direct bisulfite-PCR sequencing traces

Various technologies are currently available to quantify DNA methylation. However, rapid and simple methods for determining the DNA methylation status of CpG sites in genes still remain elusive. In this report, we describe a novel method for the rapid quantification of CpG methylation on the basis of direct bisulfite-PCR sequencing method. According to the principles of bisulfite-PCR, converting unmethylated cytosines to thymine while leaving methylated cytosines unchanged, we regard the CpG site as a SNP and estimate the methylation status of cytosines in the given CG dinucleotides by measuring the ratio of the cytosine peak height to the sum of cytosine and thymine peak heights in automated DNA sequencing traces. Furthermore, we take several effective measures to break through the 'bottleneck' problems that render the routine bisulfite sequencing method unsuitable for quantitative methylation. In comparison with pyrosequencing and bisulfite-cloning sequencing, our method is confirmed to be a simple, high-throughput and cost-effective technology for determining the methylation status of specific genes. Accordingly, this novel method is anticipated to be an efficient and economical alternative tool for rapid quantification of methylation patterns in screening large numbers of clinical samples across multiple genes.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Identification of candidate regulatory SNPs by combination of transcription-factor-binding site prediction, SNP genotyping and haploChIP

Disease-associated SNPs detected in large-scale association studies are frequently located in non-coding genomic regions, suggesting that they may be involved in transcriptional regulation. Here we describe a new strategy for detecting regulatory SNPs (rSNPs), by combining computational and experimental approaches. Whole genome ChIP-chip data for USF1 was analyzed using a novel motif finding algorithm called BCRANK. 1754 binding sites were identified and 140 candidate rSNPs were found in the predicted sites. For validating their regulatory function, seven SNPs found to be heterozygous in at least one of four human cell samples were investigated by ChIP and sequence analysis (haploChIP). In four of five cases where the SNP was predicted to affect binding, USF1 was preferentially bound to the allele containing the consensus motif. Allelic differences in binding for other proteins and histone marks further reinforced the SNPs regulatory potential. Moreover, for one of these SNPs, H3K36me3 and POLR2A levels at neighboring heterozygous SNPs indicated effects on transcription. Our strategy, which is entirely based on in vivo data for both the prediction and validation steps, can identify individual binding sites at base pair resolution and predict rSNPs. Overall, this approach can help to pinpoint the causative SNPs in complex disorders where the associated haplotypes are located in regulatory regions. Availability: BCRANK is available from Bioconductor (http://www.bioconductor.org/).

The maternally expressed WRKY transcription factor TTG2 controls lethality in interploidy crosses of Arabidopsis

The molecular mechanisms underlying lethality of F1 hybrids between diverged parents are one target of speciation research. Crosses between diploid and tetraploid individuals of the same genotype can result in F1 lethality, and this dosage-sensitive incompatibility plays a role in polyploid speciation. We have identified variation in F1 lethality in interploidy crosses of Arabidopsis thaliana and determined the genetic architecture of the maternally expressed variation via QTL mapping. A single large-effect QTL, DR. STRANGELOVE 1 (DSL1), was identified as well as two QTL with epistatic relationships to DSL1. DSL1 affects the rate of postzygotic lethality via expression in the maternal sporophyte. Fine mapping placed DSL1 in an interval encoding the maternal effect transcription factor TTG2. Maternal parents carrying loss-of-function mutations in TTG2 suppressed the F1 lethality caused by paternal excess interploidy crosses. The frequency of cellularization in the endosperm was similarly affected by both natural variation and ttg2 loss-of-function mutants. The simple genetic basis of the natural variation and effects of single-gene mutations suggests that F1 lethality in polyploids could evolve rapidly. Furthermore, the role of the sporophytically active TTG2 gene in interploidy crosses indicates that the developmental programming of the mother regulates the viability of interploidy hybrid offspring.

A simple algorithm for quantifying DNA methylation levels on multiple independent CpG sites in bisulfite genomic sequencing electropherograms

DNA methylation at cytosines is a widely studied epigenetic modification. Methylation is commonly detected using bisulfite modification of DNA followed by PCR and additional techniques such as restriction digestion or sequencing. These additional techniques are either laborious, require specialized equipment, or are not quantitative. Here we describe a simple algorithm that yields quantitative results from analysis of conventional four-dye-trace sequencing. We call this method Mquant and we compare it with the established laboratory method of combined bisulfite restriction assay (COBRA). This analysis of sequencing electropherograms provides a simple, easily applied method to quantify DNA methylation at specific CpG sites.

A non-ancestralRPGR missense mutation in families with either recessive or semi-dominant X-linked retinitis pigmentosa

Most X-linked diseases show a recessive pattern of inheritance in which female carriers are unaffected. In X-linked retinitis pigmentosa (XLRP), however, both recessive and semi-dominant inheritance patterns have been reported. We identified an Israeli family with semi-dominant XLRP due to a missense mutation (p.G275S) in the RPGR gene. The mutation was previously reported in two Danish families with recessive XLRP. Obligate carriers from the two Danish families had no visual complaints and normal to slightly reduced retinal function, while those from the Israeli family suffered from high myopia, low visual acuity, constricted visual fields, and severely reduced electroretinogram (ERG) amplitudes. The disease-related RPGR haplotype of the Israeli family was found to be different from the one found in the two Danish families, indicating that the mutation arose twice independently on different X-chromosome backgrounds. A series of genetic analyses excluded skewed X-inactivation pattern, chromosomal abnormalities, distorted RPGR expression level, and mutations in candidate genes as the cause for the differences in disease severity of female carriers. To the best of our knowledge, this is the first detailed analysis of an identical mutation causing either a recessive or a semi-dominant X-linked pattern of disease in different families. Our results indicate that an additional gene (or genes), linked to RPGR, modulate disease expression in severely affected carriers. These may be related to the high myopia concomitantly found in affected carriers from the Israeli family.

Large genomic fibrillin-1 (FBN1) gene deletions provide evidence for true haploinsufficiency in Marfan syndrome

Mutations in the FBN1 gene are the major cause of Marfan syndrome (MFS), an autosomal dominant connective tissue disorder, which displays variable manifestations in the cardiovascular, ocular, and skeletal systems. Current molecular genetic testing of FBN1 may miss mutations in the promoter region or in other noncoding sequences as well as partial or complete gene deletions and duplications. In this study, we tested for copy number variations by successively applying multiplex ligation-dependent probe amplification (MLPA) and the Affymetrix Human Mapping 500 K Array Set, which contains probes for approximately 500,000 single-nucleotide polymorphisms (SNPs) across the genome. By analyzing genomic DNA of 101 unrelated individuals with MFS or related phenotypes in whom standard genetic testing detected no mutation, we identified FBN1 deletions in two patients with MFS. Our high-resolution approach narrowed down the deletion breakpoints. Subsequent sequencing of the junctional fragments revealed the deletion sizes of 26,887 and 302,580 bp, respectively. Surprisingly, both deletions affect the putative regulatory and promoter region of the FBN1 gene, strongly indicating that they abolish transcription of the deleted allele. This expectation of complete loss of function of one allele, i.e. true haploinsufficiency, was confirmed by transcript analyses. Our findings not only emphasize the importance of screening for large genomic rearrangements in comprehensive genetic testing of FBN1 but, importantly, also extend the molecular etiology of MFS by providing hitherto unreported evidence that true haploinsufficiency is sufficient to cause MFS.

Zscan4: a novel gene expressed exclusively in late 2-cell embryos and embryonic stem cells

The first wave of transcription, called zygotic genome activation (ZGA), begins during the 2-cell stage in mouse preimplantation development and marks a vital transition from the maternal genetic to the embryonic genetic program. Utilizing DNA microarray data, we looked for genes that are expressed only during ZGA and found Zscan4, whose expression is restricted to late 2-cell stage embryos. Sequence analysis of genomic DNA and cDNA clones revealed nine paralogous genes tightly clustered in 0.85 Mb on mouse chromosome 7. Three genes are not transcribed and are thus considered pseudogenes. Among the six expressed genes named Zscan4a-Zscan4f, three - Zscan4c, Zscan4d, and Zscan4f - encode full-length ORFs with 506 amino acids. Zscan4d is a predominant transcript at the late 2-cell stage, whereas Zscan4c is a predominant transcript in embryonic stem (ES) cells. No transcripts of any Zscan4 genes are detected in any other cell types. Reduction of Zscan4 transcript levels by siRNAs delays the progression from the 2-cell to the 4-cell stage and produces blastocysts that fail to implant or proliferate in blastocyst outgrowth culture. Zscan4 thus seems to be essential for preimplantation development.

Natural selection for rash-forming genotypes of the varicella-zoster vaccine virus detected within immunized human hosts

The Oka vaccine strain is a live attenuated virus that is routinely administered to children in the United States and Europe to prevent chickenpox. It is effective and safe but occasionally produces a rash. The vaccine virus has accumulated mutations during its attenuation, but the rashes are not explained by their reversion, unlike complications reported for other viral vaccines. Indeed, most of the novel mutations distinguishing the Oka vaccine from the more virulent parental virus have not actually become fixed. Because the parental alleles are still present, the vaccine is polymorphic at >30 loci and therefore contains a mixture of related viruses. The inoculation of >40 million patients has consequently created a highly replicated evolutionary experiment that we have used to assess the competitive ability of these different viral genotypes in a human host. Using virus recovered from rash vesicles, we show that two vaccine mutations, causing amino acid substitutions in the major transactivating protein IE62, are outcompeted by the ancestral alleles. Standard interpretations of varicella disease severity concentrate on the undeniably important effects of host genotype and immune status, yet our results allow us to demonstrate that the viral genotype is associated with virulence and to identify the key sites. We propose that these loci have pleiotropic effects on the immunogenic properties of the virus, rash formation, and its epidemiological spread, which mould the evolution of its virulence. These findings are of practical importance for reducing the incidence of vaccine-associated rash and promoting public acceptance of the vaccine.

DNA pooling: a comprehensive, multi-stage association analysis of ACSL6 and SIRT5 polymorphisms in schizophrenia

Many candidate gene association studies have evaluated incomplete, unrepresentative sets of single nucleotide polymorphisms (SNPs), producing non-significant results that are difficult to interpret. Using a rapid, efficient strategy designed to investigate all common SNPs, we tested associations between schizophrenia and two positional candidate genes: ACSL6 (Acyl-Coenzyme A synthetase long-chain family member 6) and SIRT5 (silent mating type information regulation 2 homologue 5). We initially evaluated the utility of DNA sequencing traces to estimate SNP allele frequencies in pooled DNA samples. The mean variances for the DNA sequencing estimates were acceptable and were comparable to other published methods (mean variance: 0.0008, range 0-0.0119). Using pooled DNA samples from cases with schizophrenia/schizoaffective disorder (Diagnostic and Statistical Manual of Mental Disorders edition IV criteria) and controls (n=200, each group), we next sequenced all exons, introns and flanking upstream/downstream sequences for ACSL6 and SIRT5. Among 69 identified SNPs, case-control allele frequency comparisons revealed nine suggestive associations (P<0.2). Each of these SNPs was next genotyped in the individual samples composing the pools. A suggestive association with rs 11743803 at ACSL6 remained (allele-wise P=0.02), with diminished evidence in an extended sample (448 cases, 554 controls, P=0.062). In conclusion, we propose a multi-stage method for comprehensive, rapid, efficient and economical genetic association analysis that enables simultaneous SNP detection and allele frequency estimation in large samples. This strategy may be particularly useful for research groups lacking access to high throughput genotyping facilities. Our analyses did not yield convincing evidence for associations of schizophrenia with ACSL6 or SIRT5.

Detection of aneuploidies by paralogous sequence quantification

BACKGROUND

Chromosomal aneuploidies are a common cause of congenital disorders associated with cognitive impairment and multiple dysmorphic features. Pre-natal diagnosis of aneuploidies is most commonly performed by the karyotyping of fetal cells obtained by amniocentesis or chorionic villus sampling, but this method is labour intensive and requires about 14 days to complete.

METHODS

We have developed a PCR based method for the detection of targeted chromosome number abnormalities termed paralogous sequence quantification (PSQ), based on the use of paralogous genes. Paralogous sequences have a high degree of sequence identity, but accumulate nucleotide substitutions in a locus specific manner. These sequence differences, which we term paralogous sequence mismatches (PSMs), can be quantified using pyrosequencing technology, to estimate the relative dosage between different chromosomes. We designed 10 assays for the detection of trisomies of chromosomes 13, 18, and 21 and sex chromosome aneuploidies.

RESULTS

We evaluated the performance of this method on 175 DNAs, highly enriched for abnormal samples. A correct and unambiguous diagnosis was given for 119 out of 120 aneuploid samples as well as for all the controls. One sample which gave an intermediate value for the chromosome 13 assays could not be diagnosed.

CONCLUSIONS

Our data suggests that PSQ is a robust, easy to interpret, and easy to set up method for the diagnosis of common aneuploidies, and can be performed in less than 48 h, representing a competitive alternative for widespread use in diagnostic laboratories.