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.

Discrimination of the closely related A and D genomes of the hexaploid oat Avena sativa L.

A satellite DNA sequence, As120a, specific to the A-genome chromosomes in the hexaploid oat, Avena sativa L., was isolated by subcloning a fragment with internal tandem repeats from a plasmid, pAs120, that had been obtained from an Avena strigosa (As genome) genomic library. Southern and in situ hybridization showed that sequences with homology to sequences within pAs120 were dispersed throughout the genome of diploid (A and C genomes), tetraploid (AC genomes), and hexaploid (ACD genomes) Avena species. In contrast, sequences homologous to As120a were found in two A-genome species (A. strigosa and Avena longiglumis) and in the hexaploid A. sativa whereas this sequence was little amplified in the tetraploid Avena murphyi and was absent in the remaining A- and C-genome diploid species. In situ hybridization of pAs120a to hexaploid oat species revealed the distribution of elements of the As120a repeated family over both arms of 14 of 42 chromosomes of this species. By using double in situ hybridization with pAs120a and a C genome-specific probe, three sets of 14 chromosomes were revealed corresponding to the A, C, and D genomes of the hexaploid species. Simultaneous in situ hybridizations with pAs120a and ribosomal probes were used to assign the SAT chromosomes of hexaploid species to their correct genomes. This work reports a sequence able to distinguish between the closely related A and D genomes of hexaploid oats. This sequence offers new opportunities to analyze the relationships of Avena species and to explore the possible evolution of various polyploid oat species.

Chromosome-specific molecular organization of maize (Zea mays L.) centromeric regions

A set of oat–maize chromosome addition lines with individual maize (Zea mays L.) chromosomes present in plants with a complete oat (Avena sativa L.) chromosome complement provides a unique opportunity to analyze the organization of centromeric regions of each maize chromosome. A DNA sequence, MCS1a, described previously as a maize centromere-associated sequence, was used as a probe to isolate cosmid clones from a genomic library made of DNA purified from a maize chromosome 9 addition line. Analysis of six cosmid clones containing centromeric DNA segments revealed a complex organization. The MCS1a sequence was found to comprise a portion of the long terminal repeats of a retrotransposon-like repeated element, termed CentA. Two of the six cosmid clones contained regions composed of a newly identified family of tandem repeats, termed CentC. Copies of CentA and tandem arrays of CentC are interspersed with other repetitive elements, including the previously identified maize retroelements Huck and Prem2. Fluorescence in situ hybridization revealed that CentC and CentA elements are limited to the centromeric region of each maize chromosome. The retroelements Huck and Prem2 are dispersed along all maize chromosomes, although Huck elements are present in an increased concentration around centromeric regions. Significant variation in the size of the blocks of CentC and in the copy number of CentA elements, as well as restriction fragment length variations were detected within the centromeric region of each maize chromosome studied. The different proportions and arrangements of these elements and likely others provide each centromeric region with a unique overall structure.

Sequence variation in the Fanconi anemia gene FAA

Fanconi anemia (FA) is a genetically heterogeneous autosomal recessive syndrome associated with chromosomal instability, hypersensitivity to DNA crosslinking agents, and predisposition to malignancy. The gene for FA complementation group A (FAA) recently has been cloned. The cDNA is predicted to encode a polypeptide of 1,455 amino acids, with no homologies to any known protein that might suggest a function for FAA. We have used single-strand conformational polymorphism analysis to screen genomic DNA from a panel of 97 racially and ethnically diverse FA patients from the International Fanconi Anemia Registry for mutations in the FAA gene. A total of 85 variant bands were detected. Forty-five of the variants are probably benign polymorphisms, of which nine are common and can be used for various applications, including mapping studies for other genes in this region of chromosome 16q. Amplification refractory mutation system assays were developed to simplify their detection. Forty variants are likely to be pathogenic mutations. Seventeen of these are microdeletions/microinsertions associated with short direct repeats or homonucleotide tracts, a type of mutation thought to be generated by a mechanism of slipped-strand mispairing during DNA replication. A screening of 350 FA probands from the International Fanconi Anemia Registry for two of these deletions (1115–1118del and 3788–3790del) revealed that they are carried on about 2% and 5% of the FA alleles, respectively. 3788–3790del appears in a variety of ethnic groups and is found on at least two different haplotypes. We suggest that FAA is hypermutable, and that slipped-strand mispairing, a mutational mechanism recognized as important for the generation of germ-line and somatic mutations in a variety of cancer-related genes, including p53, APC, RB1, WT1, and BRCA1, may be a major mechanism for FAA mutagenesis.

Identification and characterization of a single-stranded DNA-binding protein from the archaeon Methanococcus jannaschii

Single-stranded DNA-binding proteins (SSBs) play essential roles in DNA replication, recombination, and repair in bacteria and eukarya. We report here the identification and characterization of the SSB of an archaeon, Methanococcus jannaschii. The M. jannaschii SSB (mjaSSB) has significant amino acid sequence similarity to the eukaryotic SSB, replication protein A (RPA), and contains four tandem repeats of the core single-stranded DNA (ssDNA) binding domain originally defined by structural studies of RPA. Homologous SSBs are encoded by the genomes of other archaeal species, including Methanobacterium thermoautotrophicum and Archaeoglobus fulgidus. The purified mjaSSB binds to ssDNA with high affinity and selectivity. The apparent association constant for binding to ssDNA is similar to that of RPA under comparable experimental conditions, and the affinity for ssDNA exceeds that for double-stranded DNA by at least two orders of magnitude. The binding site size for mjaSSB is ≈20 nucleotides. Given that RPA is related to mjaSSB at the sequence level and to Escherichia coli SSB at the structural level, we conclude that the SSBs of archaea, eukarya, and bacteria share a common core ssDNA-binding domain. This ssDNA-binding domain was presumably present in the common ancestor to all three major branches of life.

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.

Unexpected homology between inducible cell wall protein QID74 of filamentous fungi and BR3 salivary protein of the insect Chironomus

A gene, qid74, of mycoparasitic filamentous fungus Trichoderma harzianum and its allies encodes a cell wall protein that is induced by replacing glucose in the culture medium with chitin (simulated mycoparasitism conditions). Because no trace of this gene can be detected in related species such as Gibberella fujikuroi and Saccharomyces cerevisiae, the qid74 gene appears to have arisen de novo within the genus Trichoderma. Qid74 protein, 687 residues long, is now seen as highly conserved tandem repeats of the 59-residue-long unit. This unit itself, however, may have arisen as tandem repeats of the shorter 13-residue-long basic unit. Within the genus Trichoderma, the amino acid sequence of Qid74 proteins has been conserved in toto. The most striking is the fact that Qid74 shares 25.3% sequence identity with the carboxyl-terminal half of the 1,572-residue-long BR3 protein of the dipteran insect Chironomus tentans. BR3 protein is secreted by the salivary gland of each aquatic larva of Chironomus to form a tube to house itself. Furthermore, the consensus sequence derived from these 59-residue-long repeating units resembles those of epidermal growth factor-like domains found in divergent invertebrate and vertebrate proteins as to the positions of critical cysteine residues and homology of residues surrounding these cysteines.

High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology

The worldwide threat of tuberculosis to human health emphasizes the need to develop novel approaches to a global epidemiological surveillance. The current standard for Mycobacterium tuberculosis typing based on IS6110 restriction fragment length polymorphism (RFLP) suffers from the difficulty of comparing data between independent laboratories. Here, we propose a high-resolution typing method based on variable number tandem repeats (VNTRs) of genetic elements named mycobacterial interspersed repetitive units (MIRUs) in 12 human minisatellite-like regions of the M. tuberculosis genome. MIRU-VNTR profiles of 72 different M. tuberculosis isolates were established by PCR analysis of all 12 loci. From 2 to 8 MIRU-VNTR alleles were identified in the 12 regions in these strains, which corresponds to a potential of over 16 million different combinations, yielding a resolution power close to that of IS6110-RFLP. All epidemiologically related isolates tested were perfectly clustered by MIRU-VNTR typing, indicating that the stability of these MIRU-VNTRs is adequate to track outbreak episodes. The correlation between genetic relationships inferred from MIRU-VNTR and IS6110-RFLP typing was highly significant. Compared with IS6110-RFLP, high-resolution MIRU-VNTR typing has the considerable advantages of being fast, appropriate for all M. tuberculosis isolates, including strains that have a few IS6110 copies, and permitting easy and rapid comparison of results from independent laboratories. This typing method opens the way to the construction of digital global databases for molecular epidemiology studies of M. tuberculosis.

A link between DNA methylation and epigenetic silencing in transgenic Volvox carteri

Epigenetic silencing of foreign genes introduced into plants poses an unsolved problem for transgenic technology. Here we have used the simple multicellular green alga Volvox carteri as a model to analyse the relation of DNA methylation to transgenic silencing. Volvox DNA contains on average 1.1% 5-methylcytosine and 0.3% N6-methyladenine, as revealed by electrospray mass spectrometry and phosphoimaging of chromatographically separated 32P-labelled nucleotides. In two nuclear transformants of V.carteri, produced in 1993 by biolistic bombardment with a foreign arylsulphatase gene (C-ars), the transgene is still expressed in one (Hill 181), but not in the other (Hill 183), after an estimated 500–1000 generations. Each transformant clone contains multiple intact copies of C-ars, most of them integrated into the genome as tandem repeats. When the bisulphite genomic sequencing protocol was applied to examine two select regions of transgenic C-ars, we found that the inactivated copies (Hill 183) exhibited a high-level methylation (40%) of CpG dinucleotides, whereas the active copies (Hill 181) displayed low-level (7%) CpG methylation. These are average values from 40 PCR clones sequenced from each DNA strand in the two portions of C-ars. The observed correlation of CpG methylation and transgene inactivation in a green alga will be discussed in the light of transcriptional silencing.

The complete genome of the crenarchaeon Sulfolobus solfataricus P2

The genome of the crenarchaeon Sulfolobus solfataricus P2 contains 2,992,245 bp on a single chromosome and encodes 2,977 proteins and many RNAs. One-third of the encoded proteins have no detectable homologs in other sequenced genomes. Moreover, 40% appear to be archaeal-specific, and only 12% and 2.3% are shared exclusively with bacteria and eukarya, respectively. The genome shows a high level of plasticity with 200 diverse insertion sequence elements, many putative nonautonomous mobile elements, and evidence of integrase-mediated insertion events. There are also long clusters of regularly spaced tandem repeats. Different transfer systems are used for the uptake of inorganic and organic solutes, and a wealth of intracellular and extracellular proteases, sugar, and sulfur metabolizing enzymes are encoded, as well as enzymes of the central metabolic pathways and motility proteins. The major metabolic electron carrier is not NADH as in bacteria and eukarya but probably ferredoxin. The essential components required for DNA replication, DNA repair and recombination, the cell cycle, transcriptional initiation and translation, but not DNA folding, show a strong eukaryal character with many archaeal-specific features. The results illustrate major differences between crenarchaea and euryarchaea, especially for their DNA replication mechanism and cell cycle processes and their translational apparatus.

Copy-choice recombination mediated by DNA polymerase III holoenzyme from Escherichia coli.

Formation of deletions by recombination between short direct repeats is thought to involve either a break-join or a copy-choice process. The key step of the latter is slippage of the replication machinery between the repeats. We report that the main replicase of Escherichia coli, DNA polymerase III holoenzyme, slips between two direct repeats of 27 bp that flank an inverted repeat of approximately equal 300bp. Slippage was detected in vitro, on a single-stranded DNA template, in a primer extension assay. It requires the presence of a short (8 bp) G+C-rich sequence at the base of a hairpin that can form by annealing of the inverted repeats. It is stimulated by (i) high salt concentration, which might stabilize the hairpin, and (ii) two proteins that ensure the processivity of the DNA polymerase III holoenzyme: the single-stranded DNA binding protein and the beta subunit of the polymerase. Slippage is rather efficient under optimal reaction conditions because it can take place on >50% of template molecules. This observation supports the copy-choice model for recombination between short direct repeats.

Cloning, expression, and properties of the microtubule-stabilizing protein STOP.

Nerve cells contain abundant subpopulations of cold-stable microtubules. We have previously isolated a calmodulin-regulated brain protein, STOP (stable tubule-only polypeptide), which reconstitutes microtubule cold stability when added to cold-labile microtubules in vitro. We have now cloned cDNA encoding STOP. We find that STOP is a 100.5-kDa protein with no homology to known proteins. The primary structure of STOP includes two distinct domains of repeated motifs. The central region of STOP contains 5 tandem repeats of 46 amino acids, 4 with 98% homology to the consensus sequence. The STOP C terminus contains 28 imperfect repeats of an 11-amino acid motif. STOP also contains a putative SH3-binding motif close to its N terminus. In vitro translated STOP binds to both microtubules and Ca2+-calmodulin. When STOP cDNA is expressed in cells that lack cold-stable microtubules, STOP associates with microtubules at 37 degrees C, and stabilizes microtubule networks, inducing cold stability, nocodazole resistance, and tubulin detyrosination on microtubules in transfected cells. We conclude that STOP must play an important role in the generation of microtubule cold stability and in the control of microtubule dynamics in brain.

DNA profile match probabilities in a subdivided population: when can subdivision be ignored?

Li and Chakravarti [Li, C.C. & Chakravarti, A. (1994) Hum. Hered. 44, 100-109] compared the probability (MO) of a random match between the two DNA profiles of a pair of individuals drawn from a random-mating population to the probability (MF) of the match between a pair of random individuals drawn from a subdivided population. The level of heterogeneity in this subdivided population is measured by the parameter F, where there is no subdivision when F = 0 and increasing values of F indicate increasing subdivisions. Li and Chakravarti concluded that it is conservative to use the match probability MO, which is derived under the assumption that the two individuals are drawn from a homogeneous random-mating population without subdivision. However, MO may not be always greater than MF, even for biologically reasonable values of F. We explore here those mathematical conditions under which MO is less than MF, and we find that MO is not conservative mainly when there is an allele with a much higher frequency than all the other alleles. When empirical data for both variable number of tandem repeat (VNTR) and short tandem repeat (STR) systems are evaluated, we find that in the majority of cases MO represents a conservative probability of a match, and so the subdivision of human populations may usually be ignored for a random match, although not, of course, for relatives. Loci for which MO is not conservative should be avoided for forensic inference.

Down-regulation of NF-kappa B protein levels in activated human lymphocytes by 1,25-dihydroxyvitamin D3.

The effect of 1,25-dihydroxyvitamin D3 [1,25(OH)2)D3], a steroid hormone with immunomodulating properties, on nuclear factor kappa B (NF-kappa B) proteins was examined in in vitro activated normal human lymphocytes by Western blot analysis. Over a 72-hr period of activation, the expression of the 50-kDa NF-kappa B, p50, and its precursor, p105, was increased progressively. When cells were activated in the presence of 1,25(OH)2D3, the levels of the mature protein as well as its precursor were decreased. The effect of the hormone on the levels of p50 was demonstrable in the cytosolic and nuclear compartments; it required between 4 and 8 hr and was specific, as 25-hydroxyvitamin D3 and 24,25-dihydroxyvitamin D3 were ineffective. Besides p50, 1,25(OH)2D3 decreased the levels of another NF-kappa B protein, namely c-rel. In addition, 1,25(OH)2D3 decreased the abundance of a specific DNA-protein complex formed upon incubation of nuclear extracts from activated lymphocytes with a labeled NF-kappa B DNA binding motif. Further, 1,25(OH)2D3 inhibited the transcriptional activity of NF-kappa B in Jurkat cells transiently transfected with a construct containing four tandem repeats of the NF-kappa B binding sequence of the immunoglobulin kappa light chain gene linked to the chloramphenicol acetyltransferase reporter gene. These observations demonstrate directly that there is de novo synthesis of NF-kappa B during human lymphocyte activation and suggest that this process is hormonally regulated.

Identification of a retroviroid-like element from plants.

The biological nature of carnation small viroid-like RNA (CarSV RNA), a 275-nt circular molecule with self-cleaving hammerhead structures in its strands of both polarities, was investigated. The lack of infectivity observed in a series of transmission assays in carnation indicates that CarSV RNA, in spite of sharing structural similarities with viroid and viroid-like satellite RNAs from plants, does not belong to either of these two groups. Additional evidence in this direction comes from the observation that CarSV RNA also exists in carnation plants as DNA tandem repeats. In this respect, CarSV RNA is similar to a small transcript of a tandemly repeated DNA sequence of the newt genome. Moreover, CarSV and newt RNAs have similarities in their sequences as well as in some characteristics of their corresponding hammerhead structures. Further analyses have revealed that CarSV DNA is found directly fused to DNA sequences of carnation etched ring caulimovirus, a pararetrovirus, most likely in the form of an extrachromosomal element. The properties of the CarSV RNA/DNA system are those of a retroviroid-like element having some features in common with viroid and viroid-like satellite RNAs from plants and others with the newt transcript.