Simultaneous quantitative and allele-specific expression analysis with real competitive PCR

Detection of DNA copy number alterations by matrix-assisted laser desorption/ionization time-of-flight mass spectrometric analysis of single nucleotide polymorphisms

OBJECTIVES

Copy number alterations (CNAs) are frequently found in malignant tissues. Different approaches have been used for CNA detection. However, it is not easy to detect a large panel of CNA targets in heterogenous tumors.

METHODS

We have developed a CNAs detection approach through quantitatively analyzed allelic imbalance by allelotyping single nucleotide polymorphisms (SNPs) by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). Furthermore, the copy number changes were quantified by real-competitive PCR (rcPCR) to distinguish loss of heterozygosity (LOH) and genomic amplification. The approach was used to validate the CNA regions detected by next generation sequencing (NGS) in early-stage lung carcinoma.

RESULTS

CNAs were detected in heterogeneous DNA samples where tumor DNA is present at only 10% through the SNP based allelotyping. In addition, two different types of CNAs (loss of heterozygosity and chromosome amplification) were able to be distinguished quantitatively by rcPCR. Validation on a total of 41 SNPs from the selected CNA regions showed that copy number changes did occur, and the tissues from early-stage lung carcinoma were distinguished from normal.

CONCLUSIONS

CNA detection by MALDI-TOF MS can be used for validating potentially interesting genomic regions identified from next generation sequencing, and for detecting CNAs in tumor tissues consisting of a mixture of neoplastic and normal cells.

Whole-Genome Sequencing and iPLEX MassARRAY Genotyping Map an EMS-Induced Mutation Affecting Cell Competition in Drosophila melanogaster

Cell competition, the conditional loss of viable genotypes only when surrounded by other cells, is a phenomenon observed in certain genetic mosaic conditions. We conducted a chemical mutagenesis and screen to recover new mutations that affect cell competition between wild-type and RpS3 heterozygous cells. Mutations were identified by whole-genome sequencing, making use of software tools that greatly facilitate the distinction between newly induced mutations and other sources of apparent sequence polymorphism, thereby reducing false-positive and false-negative identification rates. In addition, we utilized iPLEX MassARRAY for genotyping recombinant chromosomes. These approaches permitted the mapping of a new mutation affecting cell competition when only a single allele existed, with a phenotype assessed only in genetic mosaics, without the benefit of complementation with existing mutations, deletions, or duplications. These techniques expand the utility of chemical mutagenesis and whole-genome sequencing for mutant identification. We discuss mutations in the Atm and Xrp1 genes identified in this screen.

Research progress in allele-specific expression and its regulatory mechanisms

Although the majority of genes are expressed equally from both alleles, some genes are differentially expressed. Organisms possess characteristics to preferentially express a particular allele under regulatory factors, which is termed allele-specific expression (ASE). It is one of the important genetic factors that lead to phenotypic variation and can be used to identify the variance of gene regulation factors. ASE indicates mechanisms such as DNA methylation, histone modifications, and non-coding RNAs function. Here, we review a broad survey of progress in ASE studies, and what this simple yet very effective approach can offer in functional genomics, and possible implications toward our better understanding of the underlying mechanisms of complex traits.

Subunit Interactions and cooperativity in the microtubule-severing AAA ATPase spastin

Spastin is a hexameric ring AAA ATPase that severs microtubules. To see whether the ring complex funnels the energy of multiple ATP hydrolysis events to the site of mechanical action, we investigate here the cooperativity of spastin. Several lines of evidence indicate that interactions among two subunits dominate the cooperative behavior: (i) the ATPase activity shows a sigmoidal dependence on the ATP concentration; (ii) ATPγS displays a mixed-inhibition behavior for normal ATP turnover; and (iii) inactive mutant subunits inhibit the activity of spastin in a hyperbolic dependence, characteristic for two interacting species. A quantitative model based on neighbor interactions fits mutant titration experiments well, suggesting that each subunit is mainly influenced by one of its neighbors. These observations are relevant for patients suffering from SPG4-type hereditary spastic paraplegia and explain why single amino acid exchanges lead to a dominant negative phenotype. In severing assays, wild type spastin is even more sensitive toward the presence of inactive mutants than in enzymatic assays, suggesting a weak coupling of ATPase and severing activity.

Increased de novo copy number variants in the offspring of older males

The offspring of older fathers have an increased risk of neurodevelopmental disorders, such as schizophrenia and autism. In light of the evidence implicating copy number variants (CNVs) with schizophrenia and autism, we used a mouse model to explore the hypothesis that the offspring of older males have an increased risk of de novo CNVs. C57BL/6J sires that were 3- and 12-16-months old were mated with 3-month-old dams to create control offspring and offspring of old sires, respectively. Applying genome-wide microarray screening technology, 7 distinct CNVs were identified in a set of 12 offspring and their parents. Competitive quantitative PCR confirmed these CNVs in the original set and also established their frequency in an independent set of 77 offspring and their parents. On the basis of the combined samples, six de novo CNVs were detected in the offspring of older sires, whereas none were detected in the control group. Two of the CNVs were associated with behavioral and/or neuroanatomical phenotypic features. One of the de novo CNVs involved Auts2 (autism susceptibility candidate 2), and other CNVs included genes linked to schizophrenia, autism and brain development. This is the first experimental demonstration that the offspring of older males have an increased risk of de novo CNVs. Our results support the hypothesis that the offspring of older fathers have an increased risk of neurodevelopmental disorders such as schizophrenia and autism by generation of de novo CNVs in the male germline.

A functional role for ADAM10 in human immunodeficiency virus type-1 replication

BACKGROUND

Gene trap insertional mutagenesis was used as a high-throughput approach to discover cellular genes participating in viral infection by screening libraries of cells selected for survival from lytic infection with a variety of viruses. Cells harboring a disrupted ADAM10 (A Disintegrin and Metalloprotease 10) allele survived reovirus infection, and subsequently ADAM10 was shown by RNA interference to be important for replication of HIV-1.

RESULTS

Silencing ADAM10 expression with small interfering RNA (siRNA) 48 hours before infection significantly inhibited HIV-1 replication in primary human monocyte-derived macrophages and in CD4⁺ cell lines. In agreement, ADAM10 over-expression significantly increased HIV-1 replication. ADAM10 down-regulation did not inhibit viral reverse transcription, indicating that viral entry and uncoating are also independent of ADAM10 expression. Integration of HIV-1 cDNA was reduced in ADAM10 down-regulated cells; however, concomitant 2-LTR circle formation was not detected, suggesting that HIV-1 does not enter the nucleus. Further, ADAM10 silencing inhibited downstream reporter gene expression and viral protein translation. Interestingly, we found that while the metalloprotease domain of ADAM10 is not required for HIV-1 replication, ADAM15 and γ-secretase (which proteolytically release the extracellular and intracellular domains of ADAM10 from the plasma membrane, respectively) do support productive infection.

CONCLUSIONS

We propose that ADAM10 facilitates replication at the level of nuclear trafficking. Collectively, our data support a model whereby ADAM10 is cleaved by ADAM15 and γ-secretase and that the ADAM10 intracellular domain directly facilitates HIV-1 nuclear trafficking. Thus, ADAM10 represents a novel cellular target class for development of antiretroviral drugs.

Allele-specific transcription of the asthma-associated PHD finger protein 11 gene (PHF11) modulated by octamer-binding transcription factor 1 (Oct-1)

BACKGROUND

Asthma is a common, chronic inflammatory airway disease of major public health importance with multiple genetic determinants. Previously, we found by positional cloning that PHD finger protein 11 (PHF11) on chromosome 13q14 modifies serum immunoglobulin E (IgE) concentrations and asthma susceptibility. No coding variants in PHF11 were identified.

OBJECTIVE

Here we investigate the 3 single nucleotide polymorphisms (SNPs) in this gene most significantly associated with total serum IgE levels--rs3765526, rs9526569, and rs1046295--for a role in transcription factor binding.

METHODS

We used electrophoretic mobility shift assays to examine the effect of the 3 SNPs on transcription factor binding in 3 cell lines relevant to asthma pathogenesis. Relative preferential expression of alleles was investigated by using the allelotyping method.

RESULTS

Electrophoretic mobility shift assays show that rs1046295 modulates allele-specific binding by the octamer-binding transcription factor 1 (Oct-1). Analysis of the relative expression levels of the 2 alleles of this SNP in heterozygous individuals showed a modest, but highly significant (P = 6.5 × 10(-16)), preferential expression of the A allele consistent with a functional role for rs1046295.

CONCLUSION

These results suggest a mechanism by which rs1046295 may act as a regulatory variant modulating transcription at this locus and altering asthma susceptibility.

Allele-specific differences in ryanodine receptor 1 mRNA expression levels may contribute to phenotypic variability in malignant hyperthermia

BACKGROUND

Malignant hyperthermia (MH) is a dominantly inherited skeletal muscle disorder that can cause a fatal hypermetabolic reaction to general anaesthetics. The primary locus of MH (MHS1 locus) in humans is linked to chromosome 19q13.1, the position of the gene encoding the ryanodine receptor skeletal muscle calcium release channel (RyR1).

METHODS

In this study, an inexpensive allele-specific PCR (AS-PCR) assay was designed that allowed the relative quantification of the two RyR1 transcripts in heterozygous samples found to be susceptible to MH (MHS). Allele-specific differences in RyR1 expression levels can provide insight into the observed variable penetrance and variations in MH phenotypes between individuals. The presence/absence of the H4833Y mutation in RYR1 transcripts was employed as a marker that allowed discrimination between the two alleles.

RESULTS

In four skeletal muscle samples and two lymphoblastoid cell lines (LCLs) from different MHS patients, the wild type allele was found to be expressed at higher levels than the mutant RyR1 allele. For both LCLs, the ratios between the wild type and mutant RYR1 alleles did not change after different incubation times with actinomycin D. This suggests that there are no allele-specific differences in RyR1 mRNA stability, at least in these cells.

CONCLUSION

The data presented here revealed for the first time allele-specific differences in RYR1 mRNA expression levels in heterozygous MHS samples, and can at least in part contribute to the observed variable penetrance and variations in MH clinical phenotypes.

Non-invasive prenatal detection of fetal trisomy 18 by RNA-SNP allelic ratio analysis using maternal plasma SERPINB2 mRNA: a feasibility study

OBJECTIVE

Non-invasive prenatal diagnosis of chromosome aneuploidies has been achieved by measuring the ratio of two alleles of a single nucleotide polymorphism (SNP) in circulating placental mRNA (the RNA-SNP allelic ratio approach) in maternal plasma. We investigated the feasibility of applying this approach for the non-invasive prenatal detection of fetal trisomy 18.

METHOD

We targeted serpin peptidase inhibitor, clade B (ovalbumin), membrane 2 (SERPINB2) mRNA, which is transcribed from chromosome 18 and is preferentially expressed by the placenta. We developed a mass-spectrometric assay to measure the SERPINB2 RNA-SNP allelic ratios in the placental samples and maternal plasma obtained from pregnancies involving euploid and trisomy 18 fetuses.

RESULTS

We were able to separate all the euploid and trisomy 18 placentas by their SERPINB2 RNA-SNP allelic ratios. The allelic ratios of the trisomy 18 placentas deviated from the reference interval established from the euploid placentas. Due to the relatively low concentrations of SERPINB2 mRNA in maternal plasma, we used pooled maternal plasma samples for analysis. We were able to identify three of the four pooled trisomy 18 plasma samples by their deviated allelic ratios when compared with the reference interval obtained from pooled euploid plasma samples.

CONCLUSION

It is feasible to detect fetal trisomy 18 non-invasively by maternal plasma SERPINB2 RNA-SNP analysis provided that sufficient quantities of plasma samples are used.

PCR with quenching probes enables the rapid detection and identification of ganciclovir-resistance-causing U69 gene mutations in human herpesvirus 6

A single-nucleotide polymorphism detection assay using PCR with quenching probes (QP-PCR) was developed for the rapid detection of antiviral drug-resistance mutations of human herpesvirus 6 (HHV-6). The mutations examined were in the HHV-6 U69 gene, and were single-base mutations in sequences known to be associated with ganciclovir (GCV) resistance in HCMV. We previously confirmed that they conferred GCV resistance to recombinant baculoviruses (Nakano et al., J. Virol. Methods 161:223-230, 2009). Six characterized mutations, including a previously reported one that encodes a GCV-sensitive kinase-activity mutant (Isegawa et al., J. Clin. Virol. 44:15-19, 2009), were used. The six mutations were separated into three groups based on their location in the U69 protein, and detected by the hybridization of three probes. We developed and validated a set of assays for these mutations using PCR followed by differential melting of a fluorescently labeled oligo probe, on a Roche Light Cycler platform. Nucleobase quenching was used to detect the hybridized probe. The optimized assay could distinguish the different mutants, and easily detected mutants representing 30% of the DNA in a mixed sample. This QP-PCR assay permitted the rapid (1.5 h), objective, and reproducible detection of drug-resistant mutations of HHV-6.

Allelic imbalance (AI) identifies novel tissue-specific cis-regulatory variation for human UGT2B15

Allelic imbalance (AI) is a powerful tool to identify cis-regulatory variation for gene expression. UGT2B15 is an important enzyme involved in the metabolism of multiple endobiotics and xenobiotics. In this study, we measured the relative expression of two alleles at this gene by using SNP rs1902023:G>T. An excess of the G over the T allele was consistently observed in liver (P<0.001), but not in breast (P=0.06) samples, suggesting that SNPs in strong linkage disequilibrium with G253T regulate UGT2B15 expression in liver. Seven such SNPs were identified by resequencing the promoter and exon 1, which define two distinct haplotypes. Reporter gene assays confirmed that one haplotype displayed approximately 20% higher promoter activity compared to the other major haplotype in liver HepG2 (P<0.001), but not in breast MCF-7 (P=0.540) cells. Reporter gene assays with additional constructs pointed to rs34010522:G>T and rs35513228:C>T as the cis-regulatory variants; both SNPs were also evaluated in LNCaP and Caco-2 cells. By ChIP, we showed that the transcription factor Nrf2 binds to the region spanning rs34010522:G>T in all four cell lines. Our results provide a good example for how AI can be used to identify cis-regulatory variation and gain insights into the tissue specific regulation of gene expression.

Detection and identification of U69 gene mutations encoded by ganciclovir-resistant human herpesvirus 6 using denaturing high-performance liquid chromatography

A denaturing high-performance liquid chromatography (dHPLC) assay was developed to detect antiviral drug-resistance mutations of human herpesvirus 6 (HHV-6). Recombinant baculoviruses were created that contained wild-type and mutant forms of the HHV-6 U69 gene, which determines sensitivity to the antiviral drug ganciclovir (GCV). The mutations causing GCV resistance in HHV-6 U69 were single-base mutations adapted from known GCV-resistant DNA sequences of HCMV, and their ability to confer GCV resistance on recombinant baculoviruses was confirmed. Six characterized mutant sequences, including one reported previously that encodes a GCV-sensitive kinase-activity mutant, were used. DNA was extracted, and the levels of homoduplex and heteroduplex DNA in the PCR products from mixed wild-type and mutant viral DNAs were determined using dHPLC. The optimized assay could distinguish the different mutants, and could detect mutants representing only 10% of the DNAs. The new assay with dHPLC readout permitted the rapid (4 h), objective, and reproducible detection of HHV-6 drug-resistance mutations.

Qualitative and quantitative genotyping using single base primer extension coupled with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MassARRAY)

Matrix-assisted laser desorption/ionization (MALDI) time-of-flight (TOF) mass spectrometry (MS) has developed over the past decade into a versatile tool for the analysis of nucleic acids and especially as a reliable genotyping platform. This chapter summarizes its use in the context of the most widely used MALDI-TOF MS genomics platform, the Sequenom MassARRAY system. MassARRAY genotyping is based upon region-specific PCR followed by allele-specific single base primer extension reactions which are then desalted, dispensed onto a silica array preloaded with matrix, and the genotyping products are resolved on the basis of mass using MALDI-TOF MS. The platform is versatile in that it can resolve multiplexed reactions (40+ separate loci per reaction), acquires and interprets data quickly, gives a quantitative output which reflects the amount of product generated for each allele within an assay for multiplexed reactions, and is highly sensitive. These characteristics coupled with integrated software for sequence annotation, assay design, data interpretation, and data storage have lead to its wide adoption and use for multiple nucleic acid analysis applications in both the realm of genomics research and molecular diagnostics.

Effect of spastic paraplegia mutations in KIF5A kinesin on transport activity

Hereditary spastic paraplegia (HSP) is a neurodegenerative disease caused by motoneuron degeneration. It is linked to at least 30 loci, among them SPG10, which causes dominant forms and originates in point mutations in the neuronal Kinesin-1 gene (KIF5A). Here, we investigate the motility of KIF5A and four HSP mutants. All mutations are single amino-acid exchanges and located in kinesin's motor or neck domain. The mutation in the neck (A361V) did not change the gliding properties in vitro, the others either reduced microtubule affinity or gliding velocity or both. In laser-trapping assays, none of the mutants moved more than a few steps along microtubules. Motility assays with mixtures of homodimeric wild-type, homodimeric mutant and heterodimeric wild-type/mutant motors revealed that only one mutant (N256S) reduces the gliding velocity at ratios present in heterozygous patients, whereas the others (K253N, R280C) do not. Attached to quantum dots as artificial cargo, mixtures involving N256S mutants produced slower cargo populations lagging behind in transport, whereas mixtures with the other mutants led to populations of quantum dots that rarely bound to microtubules. These differences indicate that the dominant inheritance of SPG10 is caused by two different mechanisms that both reduce the gross cargo flux, leading to deficient supply of the synapse.

‘Other’ applications of single nucleotide polymorphisms

Single nucleotide polymorphisms (SNPs) are the most frequent form of sequence variation in the human genome, occurring on average every 300 base pairs. Owing to their high density, SNPs are considered useful for identifying the genes associated with complex diseases. The focus of this Opinion article is the recent applications of SNPs that fall outside of disease association studies. These applications are diverse, ranging from using SNPs as qualitative markers for distinguishing individuals in mixed samples to using SNPs for quantitative genomic DNA and RNA transcript dosage assessment. These are possible owing to the millions of validated SNPs and the variety of robust SNP analysis platforms.

Allelic imbalance in gene expression as a guide to cis-acting regulatory single nucleotide polymorphisms in cancer cells

Using the relative expression levels of two SNP alleles of a gene in the same sample is an effective approach for identifying cis-acting regulatory SNPs (rSNPs). In the current study, we established a process for systematic screening for cis-acting rSNPs using experimental detection of AI as an initial approach. We selected 160 expressed candidate genes that are involved in cancer and anticancer drug resistance for analysis of AI in a panel of cell lines that represent different types of cancers and have been well characterized for their response patterns against anticancer drugs. Of these genes, 60 contained heterozygous SNPs in their coding regions, and 41 of the genes displayed imbalanced expression of the two cSNP alleles. Genes that displayed AI were subjected to bioinformatics-assisted identification of rSNPs that alter the strength of transcription factor binding. rSNPs in 15 genes were subjected to electrophoretic mobility shift assay, and in eight of these genes (APC, BCL2, CCND2, MLH1, PARP1, SLIT2, YES1, XRCC1) we identified differential protein binding from a nuclear extract between the SNP alleles. The screening process allowed us to zoom in from 160 candidate genes to eight genes that may contain functional rSNPs in their promoter regions.

MALDI-TOF mass spectrometry for quantitative, specific, and sensitive analysis of DNA and RNA

Cell-free fetal DNA and RNA released into the maternal circulation offer new opportunities to study fetal and pregnancy-associated abnormalities. Similarly, tumor cells can release cell-free DNA and RNA into the peripheral circulation, and these cell-free DNA and RNA can be used for cancer diagnosis, monitoring, and prognosis. However, these DNA and RNA often exist at very low concentrations (for fetal DNA, approximately 20 genome-equivalents (G.E.)/mL of plasma in the first trimester). The analysis is further complicated by the predominant amount of blood cell-derived DNA and RNA. MALDI-TOF mass spectrometry can provide quantitative, specific, and sensitive analysis of DNA and RNA, and thus may be a useful technology for the field.

Plasma placental RNA allelic ratio permits noninvasive prenatal chromosomal aneuploidy detection

Current methods for prenatal diagnosis of chromosomal aneuploidies involve the invasive sampling of fetal materials using procedures such as amniocentesis or chorionic villus sampling and constitute a finite risk to the fetus. Here, we outline a strategy for fetal chromosome dosage assessment that can be performed noninvasively through analysis of placental expressed mRNA in maternal plasma. We achieved noninvasive prenatal diagnosis of fetal trisomy 21 by determining the ratio between alleles of a single-nucleotide polymorphism (SNP) in PLAC4 mRNA, which is transcribed from chromosome 21 and expressed by the placenta, in maternal plasma. PLAC4 mRNA in maternal plasma was fetal derived and cleared after delivery. The allelic ratios in maternal plasma correlated with those in the placenta. Fetal trisomy 21 was detected noninvasively in 90% of cases and excluded in 96.5% of controls.

Quantitative Subtyping of Hepatitis B Virus Reveals Complex Dynamics of Ymdd Motif Mutants Development during Long-Term Lamivudine Therapy

Background/aims

Treatment of chronic hepatitis B (CHB) with lamivudine (3TC) is limited by development of drug-resistant mutants at the YMDD motif. We aimed to validate the use of mass spectrometry to detect YMDD mutants and quantify viral subpopulations.

Methods

A total of 21 Chinese patients with severe acute exacerbation of CHB treated with 3TC were studied. Serial serum samples were tested for wild-type and YMDD mutants using matrix assisted laser desorption/ionization time-of-flight mass spectrometry. INNO-LiPA assay was performed for comparison.

Results

At a median follow-up of 192 weeks, 11 patients developed YMDD mutants (six had YIDD, four had YVDD and one had YV/IDD). Mass spectrometry was concordant with INNO-LiPA in all but one patient, in which INNO-LiPA detected coexistence of YIDD and YVDD but mass spectrometry and direct sequencing detected YVDD only. Mass spectrometry was able to reliably detect a minor hepatitis B virus (HBV) subtype at 5% or above. By serial quantitative measurement, several patterns of viral dynamics were observed. In some cases, YMDD mutants dominated the whole viral population. In other cases, the proportion of YMDD mutants fluctuated with time. When more than one mutant was present (that is, YIDD and YVDD), the different mutants might dominate during different time periods.

Conclusions

Mass spectrometry is an accurate and cheap method for the detection of YMDD mutants, even in the presence of overwhelming wild-type HBV. We observed some intriguing mutant viral dynamics during 3TC treatment. Further studies are needed to clarify whether they have any clinical significance.

Matrix-assisted laser desorption/ionisation, time-of-flight mass spectrometry in genomics research

The beginning of this millennium has seen dramatic advances in genomic research. Milestones such as the complete sequencing of the human genome and of many other species were achieved and complemented by the systematic discovery of variation at the single nucleotide (SNP) and whole segment (copy number polymorphism) level. Currently most genomics research efforts are concentrated on the production of whole genome functional annotations, as well as on mapping the epigenome by identifying the methylation status of CpGs, mainly in CpG islands, in different tissues. These recent advances have a major impact on the way genetic research is conducted and have accelerated the discovery of genetic factors contributing to disease. Technology was the critical driving force behind genomics projects: both the combination of Sanger sequencing with high-throughput capillary electrophoresis and the rapid advances in microarray technologies were keys to success. MALDI-TOF MS–based genome analysis represents a relative newcomer in this field. Can it establish itself as a long-term contributor to genetics research, or is it only suitable for niche areas and for laboratories with a passion for mass spectrometry? In this review, we will highlight the potential of MALDI-TOF MS–based tools for resequencing and for epigenetics research applications, as well as for classical complex genetic studies, allele quantification, and quantitative gene expression analysis. We will also identify the current limitations of this approach and attempt to place it in the context of other genome analysis technologies.

Investigating cis-acting regulatory variation using assays of relative allelic expression

Assays of relative allelic expression offer a powerful means of exploring the influence of cis-acting regulatory variation. In this article, we provide a general assay protocol and discuss practical and theoretical considerations for the design and interpretation of relative allelic expression analysis.

DNA analysis by mass spectrometry-past, present and future

The analysis of deoxyribose nucleic acid (DNA) by mass spectrometry (MS) has evolved to where it can be used to analyze most known types of DNA and ribose nucleic acid (RNA) situations. It can efficiently deal with the analysis of DNA polymorphisms, sequences, haplotypes, human leukocyte antigen (HLA) typing, DNA methylation and RNA expression. Implementations of MS for these forms of DNA analyses are reviewed. The use of DNA analysis by MS is compared with competing technologies. Finally, an overview is given of worthwhile applications where the know-how gained so far could be used for future developments.

Computational tradeoffs in multiplex PCR assay design for SNP genotyping

Background

Multiplex PCR is a key technology for detecting infectious microorganisms, whole-genome sequencing, forensic analysis, and for enabling flexible yet low-cost genotyping. However, the design of a multiplex PCR assays requires the consideration of multiple competing objectives and physical constraints, and extensive computational analysis must be performed in order to identify the possible formation of primer-dimers that can negatively impact product yield.

Results

This paper examines the computational design limits of multiplex PCR in the context of SNP genotyping and examines tradeoffs associated with several key design factors including multiplexing level (the number of primer pairs per tube), coverage (the % of SNP whose associated primers are actually assigned to one of several available tube), and tube-size uniformity. We also examine how design performance depends on the total number of available SNPs from which to choose, and primer stringency criterial. We show that finding high-multiplexing/high-coverage designs is subject to a computational phase transition, becoming dramatically more difficult when the probability of primer pair interaction exceeds a critical threshold. The precise location of this critical transition point depends on the number of available SNPs and the level of multiplexing required. We also demonstrate how coverage performance is impacted by the number of available snps, primer selection criteria, and target multiplexing levels.

Conclusion

The presence of a phase transition suggests limits to scaling Multiplex PCR performance for high-throughput genomics applications. Achieving broad SNP coverage rapidly transitions from being very easy to very hard as the target multiplexing level (# of primer pairs per tube) increases. The onset of a phase transition can be "delayed" by having a larger pool of SNPs, or loosening primer selection constraints so as to increase the number of candidate primer pairs per SNP, though the latter may produce other adverse effects. The resulting design performance tradeoffs define a benchmark that can serve as the basis for comparing competing multiplex PCR design optimization algorithms and can also provide general rules-of-thumb to experimentalists seeking to understand the performance limits of standard multiplex PCR.

A single nucleotide polymorphism based approach for the identification and characterization of gene expression modulation using MassARRAY

Single nucleotide polymorphisms (SNPs) are the most common form of genetic variation. Their abundance and the ease with which they can be assayed have lead to their use in applications beyond simple genotyping. One such application is the quantitative determination of transcript levels associated with distinct alleles or haplotypes found in promoters and coding regions of genes. These changes in expression due to allelic variation are often associated with additional genomic or transcript modifications such as DNA methylation or RNA editing. Here, we describe the use of an integrated genetic analysis platform, based on matrix-assisted laser desorption/ionisation-time-of-flight (MALDI-TOF) to first, discover coding SNPs (cSNPs); second, use these cSNPs to identify and analyze allele-specific expression; and third, from this knowledge to further analyze methylation patterns as a putative cause for the allele-specific expression. An established model involving allele-specific expression profiles of the human tumor protein 73 (TP73) gene is presented as an example to outline and validate data obtained from the MassARRAY platform. The availability of a single integrated platform to assay stable and dynamic variation at the genomic and transcript level greatly simplifies complex functional genomic studies.