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.

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.

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.

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.

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.

Sensitive detection of DNA polymorphisms by the serial invasive signal amplification reaction

The invasive signal amplification reaction has been previously developed for quantitative detection of nucleic acids and discrimination of single-nucleotide polymorphisms. Here we describe a method that couples two invasive reactions into a serial isothermal homogeneous assay using fluorescence resonance energy transfer detection. The serial version of the assay generates more than 107 reporter molecules for each molecule of target DNA in a 4-h reaction; this sensitivity, coupled with the exquisite specificity of the reaction, is sufficient for direct detection of less than 1,000 target molecules with no prior target amplification. Here we present a kinetic analysis of the parameters affecting signal and background generation in the serial invasive signal amplification reaction and describe a simple kinetic model of the assay. We demonstrate the ability of the assay to detect as few as 600 copies of the methylene tetrahydrofolate reductase gene in samples of human genomic DNA. We also demonstrate the ability of the assay to discriminate single base differences in this gene by using 20 ng of human genomic DNA.

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.

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.

DNA transport by a micromachined Brownian ratchet device

We have micromachined a silicon-chip device that transports DNA with a Brownian ratchet that rectifies the Brownian motion of microscopic particles. Transport properties for a DNA 50-mer agree with theoretical predictions, and the DNA diffusion constant agrees with previous experiments. This type of micromachine could provide a generic pump or separation component for DNA or other charged species as part of a microscale lab-on-a-chip. A device with reduced feature size could produce a size-based separation of DNA molecules, with applications including the detection of single-nucleotide polymorphisms.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed “reference” duplex (usually a Watson-Crick duplex) and a related “test” duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.