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.

Structural characterization of an engineered tandem repeat contrasts the importance of context and sequence in protein folding

To test a different approach to understanding the relationship between the sequence of part of a protein and its conformation in the overall folded structure, the amino acid sequence corresponding to an α-helix of T4 lysozyme was duplicated in tandem. The presence of such a sequence repeat provides the protein with “choices” during folding. The mutant protein folds with almost wild-type stability, is active, and crystallizes in two different space groups, one isomorphous with wild type and the other with two molecules in the asymmetric unit. The fold of the mutant is essentially the same in all cases, showing that the inserted segment has a well-defined structure. More than half of the inserted residues are themselves helical and extend the helix present in the wild-type protein. Participation of additional duplicated residues in this helix would have required major disruption of the parent structure. The results clearly show that the residues within the duplicated sequence tend to maintain a helical conformation even though the packing interactions with the remainder of the protein are different from those of the original helix. It supports the hypothesis that the structures of individual α-helices are determined predominantly by the nature of the amino acids within the helix, rather than the structural environment provided by the rest of the protein.

Sequence-tagged microsatellite profiling (STMP): a rapid technique for developing SSR markers

We describe a technique, sequence-tagged microsatellite profiling (STMP), to rapidly generate large numbers of simple sequence repeat (SSR) markers from genomic or cDNA. This technique eliminates the need for library screening to identify SSR-containing clones and provides an ∼25-fold increase in sequencing throughput compared to traditional methods. STMP generates short but characteristic nucleotide sequence tags for fragments that are present within a pool of SSR amplicons. These tags are then ligated together to form concatemers for cloning and sequencing. The analysis of thousands of tags gives rise to a representational profile of the abundance and frequency of SSRs within the DNA pool, from which low copy sequences can be identified. As each tag contains sufficient nucleotide sequence for primer design, their conversion into PCR primers allows the amplification of corresponding full-length fragments from the pool of SSR amplicons. These fragments permit the full characterisation of a SSR locus and provide flanking sequence for the development of a microsatellite marker. Alternatively, sequence tag primers can be used to directly amplify corresponding SSR loci from genomic DNA, thereby reducing the cost of developing a microsatellite marker to the synthesis of just one sequence-specific primer. We demonstrate the utility of STMP by the development of SSR markers in bread wheat.

Targeted development of informative microsatellite (SSR) markers

We describe a novel approach, selectively amplified microsatellite (SAM) analysis, for the targeted development of informative simple sequence repeat (SSR) markers. A modified selectively amplified microsatellite polymorphic loci assay is used to generate multi-locus SSR fingerprints that provide a source of polymorphic DNA markers (SAMs) for use in genetic studies. These polymorphisms capture the repeat length variation associated with SSRs and allow their chromosomal location to be determined prior to the expense of isolating and characterising individual loci. SAMs can then be converted to locus-specific SSR markers with the design and synthesis of a single primer specific to the conserved region flanking the repeat. This approach offers a cost-efficient and rapid method for developing SSR markers for predetermined chromosomal locations and of potential informativeness. The high recovery rate of useful SSR markers makes this strategy a valuable tool for population and genetic mapping studies. The utility of SAM analysis was demonstrated by the development of SSR markers in bread wheat.

Sequence tag identification of intact proteins by matching tanden mass spectral data against sequence data bases.

Molecular and fragment ion data of intact 8- to 43-kDa proteins from electrospray Fourier-transform tandem mass spectrometry are matched against the corresponding data in sequence data bases. Extending the sequence tag concept of Mann and Wilm for matching peptides, a partial amino acid sequence in the unknown is first identified from the mass differences of a series of fragment ions, and the mass position of this sequence is defined from molecular weight and the fragment ion masses. For three studied proteins, a single sequence tag retrieved only the correct protein from the data base; a fourth protein required the input of two sequence tags. However, three of the data base proteins differed by having an extra methionine or by missing an acetyl or heme substitution. The positions of these modifications in the protein examined were greatly restricted by the mass differences of its molecular and fragment ions versus those of the data base. To characterize the primary structure of an unknown represented in the data base, this method is fast and specific and does not require prior enzymatic or chemical degradation.

Strained DNA is kinked by low concentrations of Zn2+

A novel atomic force microscope with a magnetically oscillated tip has provided unprecedented resolution of small DNA fragments spontaneously adsorbed to mica and imaged in situ in the presence of divalent ions. Kinks (localized bends of average angle 78°) were observed in axially strained minicircles consisting of tandemly repeated d(A)5 and d(GGGCC[C]) sequences. The frequency of kinks in identical minicircles increased 4-fold in the presence of 1 mM Zn2+ compared with 1 mM Mg2+. Kinking persisted in mixed Mg2+/Zn2+ electrolytes until the Zn2+ concentration dropped below 100 μM, indicating that this type of kinking may occur under physiological conditions. Kinking appears to replace intrinsic bending, and statistical analysis shows that kinks are not localized within any single sequence element. A surprisingly small free energy is associated with kink formation.

The onset and extent of genomic instability in sporadic colorectal tumor progression

Cancer cell genomes contain alterations beyond known etiologic events, but their total number has been unknown at even the order of magnitude level. By sampling colorectal premalignant polyp and carcinoma cell genomes through use of the technique inter-(simple sequence repeat) PCR, we have found genomic alterations to be considerably more abundant than expected, with the mean number of genomic events per carcinoma cell totaling approximately 11,000. Colonic polyps early in the tumor progression pathway showed similar numbers of events. These results indicate that, as with certain hereditary cancer syndromes, genomic destabilization is an early step in sporadic tumor development. Together these results support the model of genomic instability being a cause rather than an effect of malignancy, facilitating vastly accelerated somatic cell evolution, with the observed orderly steps of the colon cancer progression pathway reflecting the consequences of natural selection.

Quantitative assessment of enzyme specificity in vivo: P2 recognition by Kex2 protease defined in a genetic system

The specificity of the yeast proprotein-processing Kex2 protease was examined in vivo by using a sensitive, quantitative assay. A truncated prepro-α-factor gene encoding an α-factor precursor with a single α-factor repeat was constructed with restriction sites for cassette mutagenesis flanking the single Kex2 cleavage site (-SLDKR↓EAEA-). All of the 19 substitutions for the Lys (P2) residue in the cleavage site were made. The wild-type and mutant precursors were expressed in a yeast strain lacking the chromosomal genes encoding Kex2 and prepro-α-factor. Cleavage of the 20 sites by Kex2, expressed at the wild-type level, was assessed by using a quantitative-mating assay with an effective range greater than six orders of magnitude. All substitutions for Lys at P2 decreased mating, from 2-fold for Arg to >106-fold for Trp. Eviction of the Kex2-encoding plasmid indicated that cleavage of mutant sites by other cellular proteases was not a complicating factor. Mating efficiencies of strains expressing the mutant precursors correlated well with the specificity (kcat/KM) of purified Kex2 for comparable model peptide substrates, validating the in vivo approach as a quantitative method. The results support the conclusion that KM, which is heavily influenced by the nature of the P2 residue, is a major determinant of cleavage efficiency in vivo. P2 preference followed the rank order: Lys > Arg > Thr > Pro > Glu > Ile > Ser > Ala > Asn > Val > Cys > AsP > Gln > Gly > His > Met > Leu > Tyr > Phe > Trp.

Comparative mutational analysis of cis-acting RNA signals for translational frameshifting in HIV-1 and HTLV-2

Human immunodeficiency virus type 1 (HIV-1) and human T cell leukemia virus type II (HTLV-2) use a similar mechanism for –1 translational frameshifting to overcome the termination codon in viral RNA at the end of the gag gene. Previous studies have identified two important RNA signals for frameshifting, the slippery sequence and a downstream stem–loop structure. However, there have been somewhat conflicting reports concerning the individual contributions of these sequences. In this study we have performed a comprehensive mutational analysis of the cis-acting RNA sequences involved in HIV-1 gag–pol and HTLV-2 gag–pro frameshifting. Using an in vitro translation system we determined frameshifting efficiencies for shuffled HIV-1/HTLV-2 RNA elements in a background of HIV-1 or HTLV-2 sequences. We show that the ability of the slippery sequence and stem–loop to promote ribosomal frameshifting is influenced by the flanking upstream sequence and the nucleotides in the spacer element. A wide range of frameshift efficiency rates was observed for both viruses when shuffling single sequence elements. The results for HIV-1/HTLV-2 chimeric constructs represent strong evidence supporting the notion that the viral wild-type sequences are not designed for maximal frameshifting activity but are optimized to a level suited to efficient viral replication.

A statistical mechanical model for β-hairpin kinetics

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a β-hairpin, the minimal structural element of the antiparallel β-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, β-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 215 possible conformations and 15 × 215 possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the β-hairpin peptide.

Incorporation of glutamine repeats makes protein oligomerize: implications for neurodegenerative diseases.

Many transcription factors and some other proteins contain glutamine repeats; their abnormal expansion has been linked to several dominantly inherited neuro-degenerative diseases. Having found that poly(L-glutamine) alone forms beta-strands held together by hydrogen bonds between their amide groups, we surmised that glutamine repeats may form polar zippers, an unusual motif for protein-protein interactions. To test this hypothesis, we have engineered a Gly-Gln10-Gly peptide into the inhibitory loop of truncated chymotrypsin inhibitor 2 (CI2), a small protein from barley seeds, by both insertion and replacement. Gel filtration resolved both mutant inhibitors into at least three fractions, which analytical ultracentrifugation identified as monomers, dimers, and trimers of the recombinant protein; the truncated wild-type CI2 formed only monomers. CD difference spectra of the dimers and trimers versus wild type indicated that their glutamine repeats formed beta-pleated sheets, while those of the monomers versus wild type were more suggestive of type I beta-turns. The CD spectra of all three fractions remained unchanged even after incubation at 70 degrees C; neither the dimers nor the trimers dissociated at this temperature. We argue that the stability of all three fractions is due to the multiplicity of hydrogen bonds between extended strands of glutamine repeats in the oligomers or within a beta-hairpin formed by the single glutamine repeat of each monomer. Pathological effects may arise when expanded glutamine repeats cause proteins to acquire excessively high affinities for each other or for other proteins with glutamine repeats.

Functional copies of a human gene can be directly isolated by transformation-associated recombination cloning with a small 3′ end target sequence

Unique, small sequences (sequence tag sites) have been identified at the 3′ ends of most human genes that serve as landmarks in genome mapping. We investigated whether a single copy gene could be isolated directly from total human DNA by transformation-associated recombination (TAR) cloning in yeast using a short, 3′ unique target. A TAR cloning vector was constructed that, when linearized, contained a small amount (381 bp) of 3′ hypoxanthine phosphoribosyltransferase (HPRT) sequence at one end and an 189-bp Alu repeat at the other end. Transformation with this vector along with human DNA led to selective isolations of the entire HPRT gene as yeast artificial chromosomes (YACs) that extended from the 3′ end sequence to various Alu positions as much as 600 kb upstream. These YACs were retrofitted with a NeoR and a bacterial artificial chromosome (BAC) sequence to transfer the YACs to bacteria and subsequently the BACs to mouse cells by using a Neo selection. Most of the HPRT isolates were functional, demonstrating that TAR cloning retains the functional integrity of the isolated material. Thus, this modified version of TAR cloning, which we refer to as radial TAR cloning, can be used to isolate large segments of the human genome accurately and directly with only a small amount of sequence information.

Three functional luciferase domains in a single polypeptide chain

We report a unique case of a gene containing three homologous and contiguous repeat sequences, each of which, after excision, cloning, and expression in Escherichia coli, is shown to code for a peptide catalyzing the same reaction as the native protein, Gonyaulax polyedra luciferase (Mr = 137). This enzyme, which catalyzes the light-emitting oxidation of a linear tetrapyrrole (dinoflagellate luciferin), exhibits no sequence similarities to other luciferases in databases. Sequence analysis also reveals an unusual evolutionary feature of this gene: synonymous substitutions are strongly constrained in the central regions of each of the repeated coding sequences.

Dual roles for DNA sequence identity and the mismatch repair system in the regulation of mitotic crossing-over in yeast

Sequence divergence acts as a potent barrier to homologous recombination; much of this barrier derives from an antirecombination activity exerted by mismatch repair proteins. An inverted repeat assay system with recombination substrates ranging in identity from 74% to 100% has been used to define the relationship between sequence divergence and the rate of mitotic crossing-over in yeast. To elucidate the role of the mismatch repair machinery in regulating recombination between mismatched substrates, we performed experiments in both wild-type and mismatch repair defective strains. We find that a single mismatch is sufficient to inhibit recombination between otherwise identical sequences, and that this inhibition is dependent on the mismatch repair system. Additional mismatches have a cumulative negative effect on the recombination rate. With sequence divergence of up to approximately 10%, the inhibitory effect of mismatches results mainly from antirecombination activity of the mismatch repair system. With greater levels of divergence, recombination is inefficient even in the absence of mismatch repair activity. In both wild-type and mismatch repair defective strains, an approximate log-linear relationship is observed between the recombination rate and the level of sequence divergence.

Evolution of the Friedreich’s ataxia trinucleotide repeat expansion: Founder effect and premutations

Friedreich’s ataxia, the most frequent inherited ataxia, is caused, in the vast majority of cases, by large GAA repeat expansions in the first intron of the frataxin gene. The normal sequence corresponds to a moderately polymorphic trinucleotide repeat with bimodal size distribution. Small normal alleles have approximately eight to nine repeats whereas a more heterogeneous mode of large normal alleles ranges from 16 to 34 GAA. The latter class accounts for ≈17% of normal alleles. To identify the origin of the expansion mutation, we analyzed linkage disequilibrium between expansion mutations or normal alleles and a haplotype of five polymorphic markers within or close to the frataxin gene; 51% of the expansions were associated with a single haplotype, and the other expansions were associated with haplotypes that could be related to the major one by mutation at a polymorphic marker or by ancient recombination. Of interest, the major haplotype associated with expansion is also the major haplotype associated with the larger alleles in the normal size range and was almost never found associated with the smaller normal alleles. The results indicate that most if not all large normal alleles derive from a single founder chromosome and that they represent a reservoir for larger expansion events, possibly through “premutation” intermediates. Indeed, we found two such alleles (42 and 60 GAA) that underwent cataclysmic expansion to pathological range in a single generation. This stepwise evolution to large trinucleotide expansions already was suggested for myotonic dystrophy and fragile X syndrome and may relate to a common mutational mechanism, despite sequence motif differences.