Direct molecular haplotyping by melting curve analysis of hybridization probes: beta 2-adrenergic receptor haplotypes as an example

Direct determination of the association of multiple genetic polymorphisms, or haplotyping, in individual samples is challenging because of chromosome diploidy. Here, we describe the ability of hybridization probes, commonly used as genotyping tools, to establish single nucleotide polymorphism (SNP) haplotypes in a single step. Three haplotypes found in the beta 2-adrenergic receptor (β2AR) gene and characterized by three different SNPs combinations are presented as examples. Each combination of SNPs has a unique stability, recorded by its melting temperature, even when intervening sequences from the template must loop out during probe hybridization. In the course of this study, two haplotypes in β2AR not described previously were discovered. This approach provides a tool for molecular haplotyping that should prove useful in clinical molecular genetics diagnostics and pharmacogenetic research where methods for direct haplotyping are needed.     

Algorithm for haplotype inference via galled-tree networks with simple galls

The problem of determining haplotypes from genotypes has gained considerable prominence in the research community. Here the focus is on determining sets of SNP values on individual chromosomes since such information captures the genetic causes of diseases. The most efficient algorithmic tool for haplotyping is based on perfect phylogenetic trees. A drawback of this method is that it cannot be applied in situations when the data contains homoplasies (multiple mutations of the same character) or recombinations. Recently, Song et al. ( 2005 ) studied the two cases: haplotyping via imperfect phylogenies with a single homoplasy and via galled-tree networks with one gall. In Gupta et al. ( 2010 ), we have shown that the haplotyping via galled-tree networks is NP-hard, even if we restrict to the case when every gall contains at most 3 mutations. We present a polynomial algorithm for haplotyping via galled-tree networks with simple galls (each having two mutations) for genotype matrices which satisfy a natural condition which is implied by presence of at least one 1 in each column that contains a 2. In the end, we give the experimental results comparing our algorithm with PHASE on simulated data.     

Direct micro-haplotyping by multiple double PCR amplifications of specific alleles (MD-PASA)

Analysis of haplotypes is an important tool in population genetics, familial heredity and gene mapping. Determination of haplotypes of multiple single nucleotide polymorphisms (SNPs) or other simple mutations is time consuming and expensive when analyzing large populations, and often requires the help of computational and statistical procedures. Based on double PCR amplification of specific alleles, described previously, we have developed a simple, rapid and low-cost method for direct haplotyping of multiple SNPs and simple mutations found within relatively short specific regions or genes (micro-haplotypes). Using this method, it is possible to directly determine the physical linkage of multiple heterozygous alleles, by conducting a series of double allele-specific PCR amplification sets with simple analysis by gel electrophoresis. Application of the method requires prior information as to the sequence of the segment to be haplotyped, including the polymorphic sites. We applied the method to haplotyping of nine sites in the chicken HSP108 gene. One of the haplotypes in the population apparently arose by recombination between two existing haplotypes, and we were able to locate the point of recombination within a segment of 19 bp. We anticipate rapidly growing needs for SNP haplotyping in human (medical and pharmacogenetics), animal and plant genetics; in this context, the multiple double PCR amplifications of specific alleles (MD-PASA) method offers a useful haplotyping tool.     

SNPs detection by a single-strand specific nuclease on a PNA zip-code microarray

In this report, a reliable peptide nucleic acid (PNA) microarray-based method for accurately detecting single nucleotide polymorphism (SNP) in human genes is described. The technique relies on the mismatched cleavage activity of a single-strand specific (SSS) nuclease. PCR amplification was performed to prepare gene fragments containing the mutation sites. The amplified fragments were then employed as templates for the SSS nuclease reaction using chimeric probes, modified with biotin at the 5' end and extended with a unique anchoring zip-code complement sequence at the 3' end. The SSS nuclease promotes cleavage of heteroduplex DNAs at base mismatched positions to produce crumbled chimeric probes in the presence of imperfectly matching template strands. In contrast, the probes remain intact when they interact with perfectly matched template strands. Only the non-fragmented probes generate fluorescence signals after treatment with streptavidin-Cy3 on the PNA zip-code array. This methodology was used to successfully genotype selected Korean-specific BRCA mutation sites with wild type and mutant samples. The investigation has led to the development of a reliable SSS nuclease-based system for the diagnosis of human genetic mutations or SNPs.     

DNA suspension arrays: silencing discrete artifacts for high-sensitivity applications

Detection of low frequency single nucleotide polymorphisms (SNPs) has important implications in early screening for tumorgenesis, genetic disorders and pathogen drug resistance. Nucleic acid arrays are a powerful tool for genome-scale SNP analysis, but detection of low-frequency SNPs in a mixed population on an array is problematic. We demonstrate a model assay for HIV-1 drug resistance mutations, wherein ligase discrimination products are collected on a suspension array. In developing this system, we discovered that signal from multiple polymorphisms was obscured by two discrete hybridization artifacts. Specifically: 1) tethering of unligated probes on the template DNA elicited false signal and 2) unpredictable probe secondary structures impaired probe capture and suppressed legitimate signal from the array. Two sets of oligonucleotides were used to disrupt these structures; one to displace unligated reporter labels from the bead-bound species and another to occupy sequences which interfered with array hybridization. This artifact silencing system resulted in a mean 21-fold increased sensitivity for 29 minority variants of 17 codons in our model assay for mutations most commonly associated with HIV-1 drug resistance. Furthermore, since the artifacts we characterized are not unique to our system, their specific inhibition might improve the quality of data from solid-state microarrays as well as from the growing number of multiple analyte suspension arrays relying on sequence-specific nucleic acid target capture.     

Multiplex locked nucleic acid probes for analysis of hepatitis B virus mutants using real-time PCR

Current methods of detecting hepatitis B virus (HBV) mutations are time consuming, labor intensive, and not suitable for screening large numbers of samples. A multiplex real-time PCR approach presented in this article is a hepatitis B virus quantification method that employs the SYBR Green I dye in conjunction with wild-type HBV sequence-specific locked nucleic acid (LNA) probes. The three short LNA probes distinguished the wild-type strain or three groups of mutants (rt173, rt180/rt181, and rt202/rt204) depending on perfect-match hybrids or mismatch within one template simultaneously. Primers labeled with quencher minimized the background signals. This sensitive approach could quantify 10(2) copies of HBV virus, and as low as 1% mutants among 10(4) copies of wild-type HBV could be identified. The technique is handy and convenient, requiring only 3.5 h to analyze 30 hepatitis B surface antigen-positive serum samples. The HBV isolates were confirmed by direct sequencing. Our data indicate that real-time PCR with SYBR Green I dye is a reliable, rapid, and convenient technique for HBV quantification. Furthermore, by incorporating fluorescent LNA probes, this technique becomes handy in identifying and classifying mutations in the HBV polymerase gene. Being sensitive, specific, accurate, rapid, and convenient in nature, this technique could be a suitable diagnostic tool with wide application particularly in cases in which large volumes of clinical samples are handled.     

Does the location of a mutation determine the ability to form amyloid fibrils?

We have previously reported studies of fibril formation by a set of protein G B1 domain (beta1) variants, with mutations located around the central parallel beta-strands. In this study, we designed multiple mutations in the edge strands of beta1 to create proteins with a stability range comparable to that of the set of central mutants. All the edge variants are able to form amyloid fibrils when they are incubated at their melting temperatures. This result suggests that overall protein stability is the key determinant for amyloid formation and not the specific location of destabilizing mutations. The edge strand and variants cross-seed with each other and with members of the central variant family. Interesting fibrillar morphology was observed in some cross-seeding cases and its implications for a better understanding of nucleation and elongation events are discussed.     

Simultaneous MLPA-based multiplex point mutation and deletion analysis of the <Emphasis Type="Italic">Dystrophin</Emphasis> gene

The Multiplex Ligation-dependent Probe Amplification assay (MLPA) is the method of choice for the initial mutation screen in the analysis of a large number of genes where partial or total gene deletion is part of the mutation spectrum. Although MLPA dosage probes are usually designed to bind to normal DNA sequence to identify dosage imbalance, point mutation-specific MLPA probes can also be made. Using the dystrophin gene as a model, we have designed two MLPA probe multiplexes that are specific to a number of commonly listed point mutations in the Leiden dystrophin point mutation database (http://www.dmd.nl). The point mutation probes are designed to work simultaneously with two widely used dystrophin MLPA multiplexes, allowing both full dosage analysis and partial point mutation analysis in a single test. This approach may be adapted for other syndromes with well defined common point mutations or polymorphisms.     

High tolerance to simultaneous active‐site mutations in TEM‐1 β‐lactamase: Distinct mutational paths provide more generalized β‐lactam recognition

The diversity in substrate recognition spectra exhibited by various β‐lactamases can result from one or a few mutations in the active‐site area. Using Escherichia coli TEM‐1 β‐lactamase as a template that efficiently hydrolyses penicillins, we performed site‐saturation mutagenesis simultaneously on two opposite faces of the active‐site cavity. Residues 104 and 105 as well as 238, 240, and 244 were targeted to verify their combinatorial effects on substrate specificity and enzyme activity and to probe for cooperativity between these residues. Selection for hydrolysis of an extended‐spectrum cephalosporin, cefotaxime (CTX), led to the identification of a variety of novel mutational combinations. In vivo survival assays and in vitro characterization demonstrated a general tendency toward increased CTX and decreased penicillin resistance. Although selection was undertaken with CTX, productive binding (K_M) was improved for all substrates tested, including benzylpenicillin for which catalytic turnover (k_cat) was reduced. This indicates broadened substrate specificity, resulting in more generalized (or less specialized) variants. In most variants, the G238S mutation largely accounted for the observed properties, with additional mutations acting in an additive fashion to enhance these properties. However, the most efficient variant did not harbor the mutation G238S but combined two neighboring mutations that acted synergistically, also providing a catalytic generalization. Our exploration of concurrent mutations illustrates the high tolerance of the TEM‐1 active site to multiple simultaneous mutations and reveals two distinct mutational paths to substrate spectrum diversification.     

Transient hypermutability, chromothripsis and replication-based mechanisms in the generation of concurrent clustered mutations

Clustered mutations may be broadly defined as the presence of two or more mutations within a spatially localized genomic region on a single chromosome. Known instances vary in terms of both the number and type of the component mutations, ranging from two closely spaced point mutations to tens or even hundreds of genomic rearrangements. Although clustered mutations can represent the observable net result of independent lesions sequentially acquired over multiple cell cycles, they can also be generated in a simultaneous or quasi-simultaneous manner within a single cell cycle. This review focuses on those mechanisms known to underlie the latter type. Both gene conversion and transient hypermutability are capable of generating closely spaced multiple mutations. However, a recently described phenomenon in human cancer cells, known as 'chromothripsis', has provided convincing evidence that tens to hundreds of genomic rearrangements can sometimes be generated simultaneously via a single catastrophic event. The distinctive genomic features observed in the derivative chromosomes, together with the highly characteristic junction sequences, point to non-homologous end joining (NHEJ) as being the likely underlying mutational mechanism. By contrast, replication-based mechanisms such as microhomology-mediated break-induced replication (MMBIR) which involves serial replication slippage or serial template switching probably account for those complex genomic rearrangements that comprise multiple duplications and/or triplications.     

Restriction fragment length polymorphism analysis of the kappa-casein locus in cattle

The two common genetic variants (A and B) of bovine kappa-casein originate from two point mutations in the codons for the aminoacids in position 136 and 148. These mutations give rise to polymorphic sites for the restriction endonucleases Hin dIII, AluI, HinfI, Mbo II and TaqI. We have examined DNAs of several Italian Friesian cows and bulls of known and unknown genotype by Southern analyses using kappa-casein cDNA probes. Restriction fragment length polymorphisms (RFLPs) specific for the A and B alleles were identified for each of the above enzymes, except for AluI, which has a non-polymorphic site 12bp away from the polymorphic one. We have also found two new polymorphic sites for MboII and TaqI in the non-coding regions. These sites differentiate the A allele into two new variants, named A1 and A2. The RFLP analysis permits the characterization of kappa-casein alleles even in the absence of their expression. This should facilitate selective breeding programmes aimed at increasing the frequency of the kappa-casein B allele whose product improves the cheesemaking properties of milk.     

Transient hypermutability, chromothripsis and replication-based mechanisms in the generation of concurrent clustered mutations

Clustered mutations may be broadly defined as the presence of two or more mutations within a spatially localized genomic region on a single chromosome. Known instances vary in terms of both the number and type of the component mutations, ranging from two closely spaced point mutations to tens or even hundreds of genomic rearrangements. Although clustered mutations can represent the observable net result of independent lesions sequentially acquired over multiple cell cycles, they can also be generated in a simultaneous or quasi-simultaneous manner within a single cell cycle. This review focuses on those mechanisms known to underlie the latter type. Both gene conversion and transient hypermutability are capable of generating closely spaced multiple mutations. However, a recently described phenomenon in human cancer cells, known as ‘chromothripsis’, has provided convincing evidence that tens to hundreds of genomic rearrangements can sometimes be generated simultaneously via a single catastrophic event. The distinctive genomic features observed in the derivative chromosomes, together with the highly characteristic junction sequences, point to non-homologous end joining (NHEJ) as being the likely underlying mutational mechanism. By contrast, replication-based mechanisms such as microhomology-mediated break-induced replication (MMBIR) which involves serial replication slippage or serial template switching probably account for those complex genomic rearrangements that comprise multiple duplications and/or triplications.     

Deciphering the genetic architecture and ethnographic distribution of IRD in three ethnic populations by whole genome sequence analysis

Patients with inherited retinal dystrophies (IRDs) were recruited from two understudied populations: Mexico and Pakistan as well as a third well-studied population of European Americans to define the genetic architecture of IRD by performing whole-genome sequencing (WGS). Whole-genome analysis was performed on 409 individuals from 108 unrelated pedigrees with IRDs. All patients underwent an ophthalmic evaluation to establish the retinal phenotype. Although the 108 pedigrees in this study had previously been examined for mutations in known IRD genes using a wide range of methodologies including targeted gene(s) or mutation(s) screening, linkage analysis and exome sequencing, the gene mutations responsible for IRD in these 108 pedigrees were not determined. WGS was performed on these pedigrees using Illumina X10 at a minimum of 30X depth. The sequence reads were mapped against hg19 followed by variant calling using GATK. The genome variants were annotated using SnpEff, PolyPhen2, and CADD score; the structural variants (SVs) were called using GenomeSTRiP and LUMPY. We identified potential causative sequence alterations in 61 pedigrees (57%), including 39 novel and 54 reported variants in IRD genes. For 57 of these pedigrees the observed genotype was consistent with the initial clinical diagnosis, the remaining 4 had the clinical diagnosis reclassified based on our findings. In seven pedigrees (12%) we observed atypical causal variants, i.e. unexpected genotype(s), including 4 pedigrees with causal variants in more than one IRD gene within all affected family members, one pedigree with intrafamilial genetic heterogeneity (different affected family members carrying causal variants in different IRD genes), one pedigree carrying a dominant causative variant present in pseudo-recessive form due to consanguinity and one pedigree with a de-novo variant in the affected family member. Combined atypical and large structural variants contributed to about 20% of cases. Among the novel mutations, 75% were detected in Mexican and 50% found in European American pedigrees and have not been reported in any other population while only 20% were detected in Pakistani pedigrees and were not previously reported. The remaining novel IRD causative variants were listed in gnomAD but were found to be very rare and population specific. Mutations in known IRD associated genes contributed to pathology in 63% Mexican, 60% Pakistani and 45% European American pedigrees analyzed. Overall, contribution of known IRD gene variants to disease pathology in these three populations was similar to that observed in other populations worldwide. This study revealed a spectrum of mutations contributing to IRD in three populations, identified a large proportion of novel potentially causative variants that are specific to the corresponding population or not reported in gnomAD and shed light on the genetic architecture of IRD in these diverse global populations.     

Multiplex automated primer extension analysis: simultaneous genotyping of several polymorphisms.

Accurate and fast genotyping of single nucleotide polymorphisms (SNPs) is of significant scientific importance for linkage and association studies. We report here an automated fluorescent method we call multiplex automated primer extension analysis (MAPA) that can accurately genotype multiple known SNPs simultaneously. This is achieved by substantially improving a commercially available protocol (SNaPshot). This protocol relies on the extension of a primer that ends one nucleotide 5'of a given SNP with fluorescent dideoxy-NTPs (minisequencing), followed by analysis on an ABI PRisMS 377 Semi-Automated DNA Sequencer Our modification works by multiplexing the initial reaction that produces the DNA template for primer extension and/or multiplexing several primers (corresponding to several SNPs) in the same primer extension reaction. Then, we run each multiplexed reaction on a single gel lane. We demonstrate that MAPA can be used to genotype up to four SNPs simultaneously, even in compound heterozygote samples, with complete accuracy (based on concordance with sequencing results). We also show that primer design, unlike the DNA template purification method, can significantly affect genotyping accuracy, and we suggest useful guidelines for quick optimization.     

AGBE: a dual deaminase-mediated base editor by fusing CGBE with ABE for creating a saturated mutant population with multiple editing patterns

Establishing saturated mutagenesis in a specific gene through gene editing is an efficient approach for identifying the relationships between mutations and the corresponding phenotypes. CRISPR/Cas9-based sgRNA library screening often creates indel mutations with multiple nucleotides. Single base editors and dual deaminase-mediated base editors can achieve only one and two types of base substitutions, respectively. A new glycosylase base editor (CGBE) system, in which the uracil glycosylase inhibitor (UGI) is replaced with uracil-DNA glycosylase (UNG), was recently reported to efficiently induce multiple base conversions, including C-to-G, C-to-T and C-to-A. In this study, we fused a CGBE with ABE to develop a new type of dual deaminase-mediated base editing system, the AGBE system, that can simultaneously introduce 4 types of base conversions (C-to-G, C-to-T, C-to-A and A-to-G) as well as indels with a single sgRNA in mammalian cells. AGBEs can be used to establish saturated mutant populations for verification of the functions and consequences of multiple gene mutation patterns, including single-nucleotide variants (SNVs) and indels, through high-throughput screening.     

Evolution in vitro: analysis of a lineage of ribozymes

Background: Catalytic RNAs, or ribozymes, possessing both a genotype and a phenotype, are ideal molecules for evolution experiments in vitro. A large, heterogeneous pool of RNAs can be subjected to multiple rounds of selection, amplification and mutation, leading to the development of variants that have some desired phenotype. Such experiments allow the investigator to correlate specific genetic changes with quantifiable alterations of the catalytic properties of the RNA. In addition, patterns of evolutionary change can be discerned through a detailed examination of the genotypic composition of the evolving RNA population. Results: Beginning with a pool of 10(13) variants of the Tetrahymena ribozyme, we carried out in vitro evolution experiments that led to the generation of ribozymes with the ability to cleave an RNA substrate in the presence of Ca2+ ions, an activity that does not exist for the wild-type molecule. Over the course of 12 generations, a seven-error variant emerged that has substantial Ca(2+)-dependent RNA-cleavage activity. Advantageous mutations increased in frequency in the population according to three distinct dynamics--logarithmic, linear and transient. Through a comparative analysis of 31 individual variants, we infer how certain mutations influence the catalytic properties of the ribozyme. Conclusions: In vitro evolution experiments make it possible to elucidate important aspects of both evolutionary biology and structural biochemistry on a reasonable short time scale.     

Rapid single-base mismatch detection in genotyping for phenylketonuria

Phenylketonuria (PKU) is a metabolic disorder that results from a deficiency of hepatic phenylalanine hydroxylase (PAH). Identification of the PKU genotype is useful for predicting clinical PKU phenotype. More than 400 mutations resulting in PAH deficiency have been reported worldwide. We used a genedetecting instrument to identify the nine prevalent Japanese mutations in the PAH gene among 31 PKU patients as a preliminary study. This instrument can automatically detect mutations through the use of allele-specific oligonucleotide (ASO) capture probes, and gave results comparable to those of sequencing studies. Each country has uniquely prevalent and specific mutations causing PKU, and less than 50 types of such mutations are generally present in each country. Early genotyping of PKU makes it possible to identify the phenotype and select the optimal therapy for the disease. For early genotyping, the instrumental method described here shortens the time required for genotyping based on mRNA and/or genomic DNA of PKU parents.     

Massively parallel haplotyping on microscopic beads for the high-throughput phase analysis of single molecules

In spite of the many advances in haplotyping methods, it is still very difficult to characterize rare haplotypes in tissues and different environmental samples or to accurately assess the haplotype diversity in large mixtures. This would require a haplotyping method capable of analyzing the phase of single molecules with an unprecedented throughput. Here we describe such a haplotyping method capable of analyzing in parallel hundreds of thousands single molecules in one experiment. In this method, multiple PCR reactions amplify different polymorphic regions of a single DNA molecule on a magnetic bead compartmentalized in an emulsion drop. The allelic states of the amplified polymorphisms are identified with fluorescently labeled probes that are then decoded from images taken of the arrayed beads by a microscope. This method can evaluate the phase of up to 3 polymorphisms separated by up to 5 kilobases in hundreds of thousands single molecules. We tested the sensitivity of the method by measuring the number of mutant haplotypes synthesized by four different commercially available enzymes: Phusion, Platinum Taq, Titanium Taq, and Phire. The digital nature of the method makes it highly sensitive to detecting haplotype ratios of less than 1:10,000. We also accurately quantified chimera formation during the exponential phase of PCR by different DNA polymerases.