Computerized polymorphic marker identification: Experimental validation and a predicted human polymorphism catalog

A computational system for the prediction of polymorphic loci directly and efficiently from human genomic sequence was developed and verified. A suite of programs, collectively called pompous (polymorphic marker prediction of ubiquitous simple sequences) detects tandem repeats ranging from dinucleotides up to 250 mers, scores them according to predicted level of polymorphism, and designs appropriate flanking primers for PCR amplification. This approach was validated on an approximately 750-kilobase region of human chromosome 3p21.3, involved in lung and breast carcinoma homozygous deletions. Target DNA from 36 paired B lymphoblastoid and lung cancer lines was amplified and allelotyped for 33 loci predicted by pompous to be variable in repeat size. We found that among those 36 predominately Caucasian individuals 22 of the 33 (67%) predicted loci were polymorphic with an average heterozygosity of 0.42. Allele loss in this region was found in 27/36 (75%) of the tumor lines using these markers. pompous provides the genetic researcher with an additional tool for the rapid and efficient identification of polymorphic markers, and through a World Wide Web site, investigators can use pompous to identify polymorphic markers for their research. A catalog of 13,261 potential polymorphic markers and associated primer sets has been created from the analysis of 141,779,504 base pairs of human genomic sequence in GenBank. This data is available on our Web site (pompous.swmed.edu) and will be updated periodically as GenBank is expanded and algorithm accuracy is improved.

PCR-based subtractive hybridization and differences in gene content among strains of Helicobacter pylori

Genes that are characteristic of only certain strains of a bacterial species can be of great biologic interest. Here we describe a PCR-based subtractive hybridization method for efficiently detecting such DNAs and apply it to the gastric pathogen Helicobacter pylori. Eighteen DNAs specific to a monkey-colonizing strain (J166) were obtained by subtractive hybridization against an unrelated strain whose genome has been fully sequenced (26695). Seven J166-specific clones had no DNA sequence match to the 26695 genome, and 11 other clones were mixed, with adjacent patches that did and did not match any sequences in 26695. At the protein level, seven clones had homology to putative DNA restriction-modification enzymes, and two had homology to putative metabolic enzymes. Nine others had no database match with proteins of assigned function. PCR tests of 13 unrelated H. pylori strains by using primers specific for 12 subtracted clones and complementary Southern blot hybridizations indicated that these DNAs are highly polymorphic in the H. pylori population, with each strain yielding a different pattern of gene-specific PCR amplification. The search for polymorphic DNAs, as described here, should help identify previously unknown virulence genes in pathogens and provide new insights into microbial genetic diversity and evolution.

Posttranscriptional modification of retroviral primers is required for late stages of DNA replication

During reverse transcription of retroviral RNA, synthesis of (−) strand DNA is primed by a cellular tRNA that anneals to an 18-nt primer binding site within the 5′ long terminal repeat. For (+) strand synthesis using a (−) strand DNA template linked to the tRNA primer, only the first 18 nt of tRNA are replicated to regenerate the primer binding site, creating the (+) strand strong stop DNA intermediate and providing a 3′ terminus capable of strand transfer and further elongation. On model HIV templates that approximate the (−) strand linked to natural modified or synthetic unmodified tRNA3Lys, we find that a (+) strand strong stop intermediate of the proper length is generated only on templates containing the natural, modified tRNA3Lys, suggesting that a posttranscriptional modification provides the termination signal. In the presence of a recipient template, synthesis after strand transfer occurs only from intermediates generated from templates containing modified tRNA3Lys. Reverse transcriptase from Moloney murine leukemia virus and avian myoblastosis virus shows the same requirement for a modified tRNA3Lys template. Because all retroviral tRNA primers contain the same 1-methyl-A58 modification, our results suggest that 1-methyl-A58 is generally required for termination of replication 18 nt into the tRNA sequence, generating the (+) strand intermediate, strand transfer, and subsequent synthesis of the entire (+) strand. The possibility that the host methyl transferase responsible for methylating A58 may provide a target for HIV chemotherapy is discussed.

Soldier caste-specific gene expression in the mandibular glands of Hodotermopsis japonica (Isoptera: Termopsidae)

Although “polymorphic castes” in social insects are well known as one of the most important phenomena of polyphenism, few studies of caste-specific gene expressions have been performed in social insects. To identify genes specifically expressed in the soldier caste of the Japanese damp-wood termite Hodotermopsis japonica, we employed the differential-display method using oligo(dT) and arbitrary primers, compared mRNA from the heads of mature soldiers and pseudergates (worker caste), and identified a clone (PCR product) 329 bp in length termed SOL1. Northern blot analysis showed that the SOL1 mRNA is about 1.0 kb in length and is expressed specifically in mature soldiers, but not in pseudergates, even in the presoldier induction by juvenile hormone analogue, suggesting that the product is specific for terminally differentiated soldiers. By using the method of 5′- and 3′-rapid amplification of cDNA ends, we isolated the full length of SOL1 cDNA, which contained an ORF with a putative signal peptide at the N terminus. The sequence showed no significant homology with any other known protein sequences. In situ hybridization analysis showed that SOL1 is expressed specifically in the mandibular glands. These results strongly suggest that the SOL1 gene encodes a secretory protein specifically synthesized in the mandibular glands of the soldiers. Histological observations revealed that the gland actually develops during the differentiation into the soldier caste.

A serotonin transporter gene intron 2 polymorphic region, correlated with affective disorders, has allele-dependent differential enhancer-like properties in the mouse embryo

Polymorphic regions consisting of a variable number of tandem repeats within intron 2 of the gene coding for the serotonin transporter protein 5-HTT have been associated with susceptibility to affective disorders. We have cloned two of these intronic polymorphisms, Stin2.10 and Stin2.12, into an expression vector containing a heterologous minimal promoter and the bacterial LacZ reporter gene. These constructs were then used to produce transgenic mice. In embryonic day 10.5 embryos, both Stin2.10 and Stin2.12 produced consistent β-galactosidase expression in the embryonic midbrain, hindbrain, and spinal cord floor plate. However, we observed that the levels of β-galactosidase expression produced by both the Stin2.10 and Stin2.12 within the rostral hindbrain differed significantly at embryonic day 10.5. Our data suggest that these polymorphic variable number of tandem repeats regions act as transcriptional regulators and have allele-dependent differential enhancer-like properties within an area of the hindbrain where the 5-HTT gene is known to be transcribed at this stage of development.

Direct interaction of flagellin termini essential for polymorphic ability of flagellar filament

We report the structures of flagellar filaments reconstituted from various flagellins with small terminal truncations. Flagellins from Salmonella typhimurium strains SJW1103 (wild type), SJW1660, and SJW1655 were used, which form a left-handed supercoil, the L- and R-type straight forms, respectively. Structure analyses were done by electron cryomicroscopy and helical image reconstruction with a help of x-ray fiber diffraction for determining precise helical symmetries. Truncation of either terminal region, irrespective of the original flagellin species, results in a straight filament having a helical symmetry distinct either from the L- or R-type. This filament structure is named Lt-type. Although the local subunit packing is similar in all three types, a close comparison shows that the Lt-type packing is almost identical to the R-type but distinct from the L-type, which demonstrates the strong two-state preference of the subunit interactions. The structure clearly suggests that both termini are located in the inner tube of the concentric double-tubular structure of the filament core, and their proper interaction is responsible for the correct folding of fairly large terminal regions that form the inner tube. The double tubular structure appears to be essential for the polymorphic ability of flagellar filaments, which is required for the swimming–tumbling of bacterial taxis.

Somatic and germ cell parasitism in a colonial ascidian: Possible role for a highly polymorphic allorecognition system

A colonial protochordate, Botryllus schlosseri, undergoes a natural transplantation reaction in the wild that results alternatively in colony fusion (chimera formation) or inflammatory rejection. A single, highly polymorphic histocompatibility locus (called Fu/HC) is responsible for rejection versus fusion. Gonads are seeded and gametogenesis can occur in colonies well after fusion, and involves circulating germ-line progenitors. Buss proposed that colonial organisms might develop self/non-self histocompatibility systems to limit the possibility of interindividual germ cell “parasitism” (GCP) to histocompatible kin [Buss, L. W. (1982) Proc. Natl. Acad. Sci. USA 79, 5337–5341 and Buss, L. W. (1987) The Evolution of Individuality (Princeton Univ. Press, Princeton]. Here we demonstrate in laboratory and field experiments that both somatic cell and (more importantly) germ-line parasitism are a common occurrence in fused chimeras. These experiments support the tenet in Buss’s hypothesis that germ cell and somatic cell parasitism can occur in fused chimeras and that a somatic appearance may mask the winner of a gametic war. They also provide an interesting challenge to develop formulas that describe the inheritance of competing germ lines rather than competing individuals. The fact that fused B. schlosseri have higher rates of GCP than unfused colonies additionally provides a rational explanation for the generation and maintenance of a high degree of Fu/HC polymorphism, largely limiting GCP to sibling offspring.

RNA: a method to specifically inhibit PCR amplification of known members of a multigene family by degenerate primers

The polymerase chain reaction (PCR) is a versatile method to amplify specific DNA with oligonucleotide primers. By designing degenerate PCR primers based on amino acid sequences that are highly conserved among all known gene family members, new members of a multigene family can be identified. The inherent weakness of this approach is that the degenerate primers will amplify previously identified, in addition to new, family members. To specifically address this problem, we synthesized a specific RNA for each known family member so that it hybridized to one strand of the template, adjacent to the 3′-end of the primer, allowing the degenerate primer to bind yet preventing extension by DNA polymerase. To test our strategy, we used known members of the soluble, nitric oxide-sensitive guanylyl cyclase family as our templates and degenerate primers that discriminate this family from other guanylyl cyclases. We demonstrate that amplification of known members of this family is effectively and specifically inhibited by the corresponding RNAs, alone or in combination. This robust method can be adapted to any application where multiple PCR products are amplified, as long as the sequence of the desired and the undesired PCR product(s) is sufficiently distinct between the primers.

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.

Evidence for clonal propagation in natural isolates of Plasmodium falciparum from Venezuela

We have analyzed 75 isolates of Plasmodium falciparum, collected in Venezuela during both the dry (November) and rainy (May–July) seasons, with a range of genetic markers including antigen genes and 14 random amplified polymorphic DNA (RAPD) primers. Thirteen P. falciparum stocks from Kenya and four other Plasmodium species are included in the analysis for comparison. Cross-hybridization shows that the 14 RAPD primers reveal 14 separate regions of the parasite's genome. The P. falciparum isolates are a monophyletic clade, significantly different from the other Plasmodium species. We identify three RAPD characters that could be useful as “tags” for rapid species identification. The Venezuelan genotypes fall into two discrete genetic subdivisions associated with either the dry or the rainy season; the isolates collected in the rainy season exhibit greater genetic diversity. There is significant linkage disequilibrium in each seasonal subsample and in the full sample. In contrast, no linkage disequilibrium is detected in the African sample. These results support the hypothesis that the population structure of P. falciparum in Venezuela, but not in Africa, is predominantly clonal. However, the impact of genetic recombination on Venezuelan P. falciparum seems higher than in parasitic species with long-term clonal evolution like Trypanosoma cruzi, the agent of Chagas' disease. The genetic structure of the Venezuelan samples is similar to that of Escherichia coli, a bacterium that propagates clonally, with occasional genetic recombination.

Multiplex genotype determination at a large number of gene loci.

To facilitate large-scale genotype analysis, an efficient PCR-based multiplex approach has been developed. For simultaneously amplifying the target sequences at a large number of genetic loci, locus-specific primers containing 5' universal tails are used. Attaching the universal tails to the target sequences in the initial PCR steps allows replacement of all specific primers with a pair of primers identical to the universal tails and converts the multiplex amplification into "uniplex." Simultaneous amplification of 26 genetic loci with this approach is described. The multiplex amplification can be coupled with genotype determination. By incorporating a single-base mismatch between a primer and the template into the target sequences, a polymorphic site can be converted into a desirable restriction fragment length polymorphism when it is necessary. In this way, the allelic PCR products for the polymorphic loci can be discriminated by gel electrophoresis after restriction enzyme digestion. In this study, 32 loci were typed in such a multiplex way.

Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines.

Simple sequence repeats (SSRs), consisting of tandemly repeated multiple copies of mono-, di-, tri-, or tetranucleotide motifs, are ubiquitous in eukaryotic genomes and are frequently used as genetic markers, taking advantage of their length polymorphism. We have examined the polymorphism of such sequences in the chloroplast genomes of plants, by using a PCR-based assay. GenBank searches identified the presence of several (dA)n.(dT)n mononucleotide stretches in chloroplast genomes. A chloroplast (cp) SSR was identified in three pine species (Pinus contorta, Pinus sylvestris, and Pinus thunbergii) 312 bp upstream of the psbA gene. DNA amplification of this repeated region from 11 pine species identified nine length variants. The polymorphic amplified fragments were isolated and the DNA sequence was determined, confirming that the length polymorphism was caused by variation in the length of the repeated region. In the pines, the chloroplast genome is transmitted through pollen and this PCR assay may be used to monitor gene flow in this genus. Analysis of 305 individuals from seven populations of Pinus leucodermis Ant. revealed the presence of four variants with intrapopulational diversities ranging from 0.000 to 0.629 and an average of 0.320. Restriction fragment length polymorphism analysis of cpDNA on the same populations previously failed to detect any variation. Population subdivision based on cpSSR was higher (Gst = 0.22, where Gst is coefficient of gene differentiation) than that revealed in a previous isozyme study (Gst = 0.05). We anticipate that SSR loci within the chloroplast genome should provide a highly informative assay for the analysis of the genetic structure of plant populations.

Complete congruence between gene diversity estimates derived from genotypic data at enzyme and random amplified polymorphic DNA loci in black spruce.

Controversy still exists over the adaptive nature of variation of enzyme loci. In conifers, random amplified polymorphic DNAs (RAPDs) represent a class of marker loci that is unlikely to fall within or be strongly linked to coding DNA. We have compared the genetic diversity in natural populations of black spruce [Picea mariana (Mill.) B.S.P.] using genotypic data at allozyme loci and RAPD loci as well as phenotypic data from inferred RAPD fingerprints. The genotypic data for both allozymes and RAPDs were obtained from at least six haploid megagametophytes for each of 75 sexually mature individuals distributed in five populations. Heterozygosities and population fixation indices were in complete agreement between allozyme loci and RAPD loci. In black spruce, it is more likely that the similar levels of variation detected at both enzyme and RAPD loci are due to such evolutionary forces as migration and the mating system, rather than to balancing selection and overdominance. Furthermore, we show that biased estimates of expected heterozygosity and among-population differentiation are obtained when using allele frequencies derived from dominant RAPD phenotypes.

Fluorescence energy transfer dye-labeled primers for DNA sequencing and analysis.

Fluorescent dye-labeled DNA primers have been developed that exploit fluorescence energy transfer (ET) to optimize the absorption and emission properties of the label. These primers carry a fluorescein derivative at the 5' end as a common donor and other fluorescein and rhodamine derivatives attached to a modified thymidine residue within the primer sequence as acceptors. Adjustment of the donor-acceptor spacing through the placement of the modified thymidine in the primer sequence allowed generation of four primers, all having strong absorption at a common excitation wavelength (488 nm) and fluorescence emission maxima of 525, 555, 580, and 605 nm. The ET efficiency of these primers ranges from 65% to 97%, and they exhibit similar electrophoretic mobilities by gel electrophoresis. With argon-ion laser excitation, the fluorescence of the ET primers and of the DNA sequencing fragments generated with ET primers is 2- to 6-fold greater than that of the corresponding primers or fragments labeled with single dyes. The higher fluorescence intensity of the ET primers allows DNA sequencing with one-fourth of the DNA template typically required when using T7 DNA polymerase. With single-stranded M13mp18 DNA as the template, a typical sequencing reaction with ET primers on a commercial sequencer provided DNA sequences with 99.8% accuracy in the first 500 bases. ET primers should be generally useful in the development of other multiplex DNA sequencing and analysis methods.

Identification of mutations in the Wiskott-Aldrich syndrome gene and characterization of a polymorphic dinucleotide repeat at DXS6940, adjacent to the disease gene.

The Wiskott-Aldrich syndrome (WAS) is an X-chromosome-linked recessive disease characterized by eczema, thrombocytopenia, and immunodeficiency. The disease gene has been localized to the proximal short arm of the X chromosome and recently isolated through positional cloning. The function of the encoded protein remains undetermined. In this study we have characterized mutations in 12 unrelated patients to confirm the identity of the disease gene. We have also revised the coding sequence and genomic structure for the WAS gene. To analyze further the transmittance of the disease gene, we have characterized a polymorphic microsatellite at the DXS6940 locus within 30 kb of the gene and demonstrate the inheritance of the affected alleles in families with a history of WAS.