Alkaline-mediated differential interaction (AMDI): A simple automatable single-nucleotide polymorphism assay

The key requirements for high-throughput single-nucleotide polymorphism (SNP) typing of DNA samples in large-scale disease case-control studies are automatability, simplicity, and robustness, coupled with minimal cost. In this paper we describe a fluorescence technique for the detection of SNPs that have been amplified by using the amplification refractory mutation system (ARMS)-PCR procedure. Its performance was evaluated using 32 sequence-specific primer mixes to assign the HLA-DRB alleles to 80 lymphoblastoid cell line DNAs chosen from our database for their diversity. All had been typed previously by alternative methods, either direct sequencing or gel electrophoresis. We believe the detection system that we call AMDI (alkaline-mediated differential interaction) satisfies the above criteria and is suitable for general high-throughput SNP typing.

A single nucleotide polymorphism in the RAD51 gene modifies cancer risk in BRCA2 but not BRCA1 carriers

BRCA1 and BRCA2 carriers are at increased risk for both breast and ovarian cancer, but estimates of lifetime risk vary widely, suggesting their penetrance is modified by other genetic and/or environmental factors. The BRCA1 and BRCA2 proteins function in DNA repair in conjunction with RAD51. A preliminary report suggested that a single nucleotide polymorphism in the 5′ untranslated region of RAD51 (135C/G) increases breast cancer risk in BRCA1 and BRCA2 carriers. To investigate this effect we studied 257 female Ashkenazi Jewish carriers of one of the common BRCA1 (185delAG, 5382insC) or BRCA2 (6174delT) mutations. Of this group, 164 were affected with breast and/or ovarian cancer and 93 were unaffected. RAD51 genotyping was performed on all subjects. Among BRCA1 carriers, RAD51-135C frequency was similar in healthy and affected women [6.1% (3 of 49) and 9.9% (12 of 121), respectively], and RAD-135C did not influence age of cancer diagnosis [Hazard ratio (HR) = 1.18 for disease in RAD51-135C heterozygotes, not significant]. However, in BRCA2 carriers, RAD51-135C heterozygote frequency in affected women was 17.4% (8 of 46) compared with 4.9% (2 of 41) in unaffected women (P = 0.07). Survival analysis in BRCA2 carriers showed RAD51-135C increased risk of breast and/or ovarian cancer with an HR of 4.0 [95% confidence interval 1.6–9.8, P = 0.003]. This effect was largely due to increased breast cancer risk with an HR of 3.46 (95% confidence interval 1.3–9.2, P = 0.01) for breast cancer in BRCA2 carriers who were RAD51-135C heterozygotes. RAD51 status did not affect ovarian cancer risk. These results show RAD51-135C is a clinically significant modifier of BRCA2 penetrance, specifically in raising breast cancer risk at younger ages.

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.

The human XPG gene: gene architecture, alternative splicing and single nucleotide polymorphisms

Defects in the XPG DNA repair endonuclease gene can result in the cancer-prone disorders xeroderma pigmentosum (XP) or the XP–Cockayne syndrome complex. While the XPG cDNA sequence was known, determination of the genomic sequence was required to understand its different functions. In cells from normal donors, we found that the genomic sequence of the human XPG gene spans 30 kb, contains 15 exons that range from 61 to 1074 bp and 14 introns that range from 250 to 5763 bp. Analysis of the splice donor and acceptor sites using an information theory-based approach revealed three splice sites with low information content, which are components of the minor (U12) spliceosome. We identified six alternatively spliced XPG mRNA isoforms in cells from normal donors and from XPG patients: partial deletion of exon 8, partial retention of intron 8, two with alternative exons (in introns 1 and 6) and two that retained complete introns (introns 3 and 9). The amount of alternatively spliced XPG mRNA isoforms varied in different tissues. Most alternative splice donor and acceptor sites had a relatively high information content, but one has the U12 spliceosome sequence. A single nucleotide polymorphism has allele frequencies of 0.74 for 3507G and 0.26 for 3507C in 91 donors. The human XPG gene contains multiple splice sites with low information content in association with multiple alternatively spliced isoforms of XPG mRNA.

Genetic epidemiology of single-nucleotide polymorphisms

On the causal hypothesis, most genetic determinants of disease are single-nucleotide polymorphisms (SNPs) that are likely to be selected as markers for positional cloning. On the proximity hypothesis, most disease determinants will not be included among markers but may be detected through linkage disequilibrium with other SNPs. In that event, allelic association among SNPs is an essential factor in positional cloning. Recent simulation based on monotonic population expansion suggests that useful association does not usually extend beyond 3 kb. This is contradicted by significant disequilibrium at much greater distances, with corresponding reduction in the number of SNPs required for a cost-effective genome scan. A plausible explanation is that cyclical expansions follow population bottlenecks that establish new disequilibria. Data on more than 1,000 locus pairs indicate that most disequilibria trace to the Neolithic, with no apparent difference between haplotypes that are random or selected through a major disease gene. Short duration may be characteristic of alleles contributing to disease susceptibility and haplotypes characteristic of particular ethnic groups. Alleles that are highly polymorphic in all ethnic groups may be older, neutral, or advantageous, in weak disequilibrium with nearby markers, and therefore less useful for positional cloning of disease genes. Significant disequilibrium at large distance makes the number of suitably chosen SNPs required for genome screening as small as 30,000, or 1 per 100 kb, with greater density (including less common SNPs) reserved for candidate regions.

A single nucleotide in the SMN gene regulates splicing and is responsible for spinal muscular atrophy

SMN1 and SMN2 (survival motor neuron) encode identical proteins. A critical question is why only the homozygous loss of SMN1, and not SMN2, results in spinal muscular atrophy (SMA). Analysis of transcripts from SMN1/SMN2 hybrid genes and a new SMN1 mutation showed a direct relationship between presence of disease and exon 7 skipping. We have reported previously that the exon-skipped product SMNΔ7 is partially defective for self-association and SMN self-oligomerization correlated with clinical severity. To evaluate systematically which of the five nucleotides that differ between SMN1 and SMN2 effect alternative splicing of exon 7, a series of SMN minigenes was engineered and transfected into cultured cells, and their transcripts were characterized. Of these nucleotide differences, the exon 7 C-to-T transition at codon 280, a translationally silent variance, was necessary and sufficient to dictate exon 7 alternative splicing. Thus, the failure of SMN2 to fully compensate for SMN1 and protect from SMA is due to a nucleotide exchange (C/T) that attenuates activity of an exonic enhancer. These findings demonstrate the molecular genetic basis for the nature and pathogenesis of SMA and illustrate a novel disease mechanism. Because individuals with SMA retain the SMN2 allele, therapy targeted at preventing exon 7 skipping could modify clinical outcome.

A common polymorphism associated with antibiotic-induced cardiac arrhythmia

Drug-induced long QT syndrome (LQTS) is a prevalent disorder of uncertain etiology that predisposes to sudden death. KCNE2 encodes MinK-related peptide 1 (MiRP1), a subunit of the cardiac potassium channel IKr that has been associated previously with inherited LQTS. Here, we examine KCNE2 in 98 patients with drug-induced LQTS, identifying three individuals with sporadic mutations and a patient with sulfamethoxazole-associated LQTS who carried a single-nucleotide polymorphism (SNP) found in ≈1.6% of the general population. While mutant channels showed diminished potassium flux at baseline and wild-type drug sensitivity, channels with the SNP were normal at baseline but inhibited by sulfamethoxazole at therapeutic levels that did not affect wild-type channels. We conclude that allelic variants of MiRP1 contribute to a significant fraction of cases of drug-induced LQTS through multiple mechanisms and that common sequence variations that increase the risk of life-threatening drug reactions can be clinically silent before drug exposure.

Detection of mitochondrial single nucleotide polymorphisms using a primer elongation reaction on oligonucleotide microarrays

We have developed a novel allele-specific primer elongation protocol using a DNA polymerase on oligonucleotide chips. Oligonucleotide primers carrying polymorphic sites at their free 3́end were covalently bound to glass slides. The generation of single-stranded targets of genomic DNA containing single nuclotide polymorphisms (SNPs) to be typed was achieved by an asymmetric PCR reaction or exonuclease treatment of phosphothioate (PTO)-modified PCR products. In the presence of DNA polymerase and all four dNTPs, with Cy3-dUTP replacing dTTP, allele-specific extension of the immobilized primers took place along a stretch of target DNA sequence. The yield of elongated products was increased by repeated reaction cycles. We performed multiplexed assays with many small DNA targets, or used single targets of up to 4.4 kb mitochondrial DNA (mtDNA) sequence to detect multiple SNPs in one reaction. The latter approach greatly simplifies preamplification of SNP-containing regions, thereby providing a framework for typing hundreds of mtDNA polymorphisms.

Operational RNA code for amino acids: species-specific aminoacylation of minihelices switched by a single nucleotide.

The genetic code is based on aminoacylation reactions where specific amino acids are attached to tRNAs bearing anticodon trinucleotides. However, the anticodon-independent specific aminoacylation of RNA minihelix substrates by bacterial and yeast tRNA synthetases suggested an operational RNA code for amino acids whereby specific RNA sequences/structures in tRNA acceptor stems correspond to specific amino acids. Because of the possible significance of the operational RNA code for the development of the genetic code, we investigated aminoacylation of synthetic RNA minihelices with a human enzyme to understand the sequences needed for that aminoacylation compared with those needed for a microbial system. We show here that the species-specific aminoacylation of glycine tRNAs is recapitulated by a species-specific aminoacylation of minihelices. Although the mammalian and Escherichia coli minihelices differ at 6 of 12 base pairs, two of the three nucleotides essential for aminoacylation by the E. coli enzyme are conserved in the mammalian minihelix. The two conserved nucleotides were shown to be also important for aminoacylation of the mammalian minihelix by the human enzyme. A simple interchange of the differing nucleotide enabled the human enzyme to now charge the bacterial substrate and not the mammalian minihelix. Conversely, this interchange made the bacterial enzyme specific for the mammalian substrate. Thus, the positional locations (if not the actual nucleotides) for the operational RNA code for glycine appear conserved from bacteria to mammals.

Differential subcellular mRNA targeting: deletion of a single nucleotide prevents the transport to axons but not to dendrites of rat hypothalamic magnocellular neurons.

It has previously been shown that mRNA encoding the arginine vasopressin (AVP) precursor is targeted to axons of rat magnocellular neurons of the hypothalamo-neurohypophyseal tract. In the homozygous Brattle-boro rat, which has a G nucleotide deletion in the coding region of the AVP gene, no such targeting is observed although the gene is transcribed. RNase protection and heteroduplex analyses demonstrate that, in heterozygous animals, which express both alleles of the AVP gene, the wild-type but not the mutant transcript is subject to axonal compartmentation. In contrast, wild-type and mutant AVP mRNAs are present in dendrites. These data suggest the existence of different mechanisms for mRNA targeting to the two subcellular compartments. Axonal mRNA localization appears to take place after protein synthesis; the mutant transcript is not available for axonal targeting because it lacks a stop codon preventing its release from ribosomes. Dendritic compartmentation, on the other hand, is likely to precede translation and, thus, would be unable to discriminate between the two mRNAs.

Adaptive evolution of a candidate gene for aging in Drosophila

Examination of the phenotypic effects of specific mutations has been extensively used to identify candidate genes affecting traits of interest. However, such analyses do not reveal anything about the evolutionary forces acting at these loci, or whether standing allelic variation contributes to phenotypic variance in natural populations. The Drosophila gene methuselah (mth) has been proposed as having major effects on organismal stress response and longevity phenotype. Here, we examine patterns of polymorphism and divergence at mth in population level samples of Drosophila melanogaster, D. simulans, and D. yakuba. Mth has experienced an unusually high level of adaptive amino acid divergence concentrated in the intra- and extracellular loop domains of the receptor protein, suggesting the historical action of positive selection on those regions of the molecule that modulate signal transduction. Further analysis of single nucleotide polymorphisms (SNPs) in D. melanogaster provided evidence for contemporary and spatially variable selection at the mth locus. In ten surveyed populations, the most common mth haplotype exhibited a 40% cline in frequency that coincided with population level differences in multiple life-history traits including lifespan. This clinal pattern was not associated with any particular SNP in the coding region, indicating that selection is operating at a closely linked site that may be involved in gene expression. Together, these consistently nonneutral patterns of inter- and intraspecific variation suggest adaptive evolution of a signal transduction pathway that may modulate lifespan in nature.

Population structure and recent evolution of Plasmodium falciparum

Plasmodium falciparum is the agent of malignant malaria, one of mankind's most severe maladies. The parasite exhibits antigenic polymorphisms that have been postulated to be ancient. We have proposed that the extant world populations of P. falciparum have derived from one single parasite, a cenancestor, within the last 5,000–50,000 years. This inference derives from the virtual or complete absence of synonymous nucleotide polymorphisms at genes not involved in immune or drug responses. Seeking to conciliate this claim with extensive antigenic polymorphism, we first note that allele substitutions or polymorphisms can arise very rapidly, even in a single generation, in large populations subject to strong natural selection. Second, new alleles can arise not only by single-nucleotide mutations, but also by duplication/deletion of short simple-repeat DNA sequences, a process several orders of magnitude faster than single-nucleotide mutation. We analyze three antigenic genes known to be extremely polymorphic: Csp, Msp-1, and Msp-2. We identify regions consisting of tandem or proximally repetitive short DNA sequences, including some previously unnoticed. We conclude that the antigenic polymorphisms are consistent with the recent origin of the world populations of P. falciparum inferred from the analysis of nonantigenic genes.