Semidwarf (sd-1), "green revolution" rice, contains a defective gibberellin 20-oxidase gene

The introduction of semidwarf rice (Oryza sativa L.) led to record yield increases throughout Asia in the 1960s. The major semidwarfing allele, sd-1, is still extensively used in modern rice cultivars. The phenotype of sd-1 is consistent with dwarfism that results from a deficiency in gibberellin (GA) plant growth hormones. We propose that the semidwarf (sd-1) phenotype is the result of a deficiency of active GAs in the elongating stem arising from a defective 20-oxidase GA biosynthetic enzyme. Sequence data from the rice genome was combined with previous mapping studies to locate a putative GA 20-oxidase gene (Os20ox2) at the predicted map location of sd-1 on chromosome 1. Two independent sd-1 alleles contained alterations within Os20ox2: a deletion of 280 bp within the coding region of Os20ox2 was predicted to encode a nonfunctional protein in an indica type semidwarf (Doongara), whereas a substitution in an amino acid residue (Leu-266) that is highly conserved among dioxygenases could explain loss of function of Os20ox2 in a japonica semidwarf (Calrose76). The quantification of GAs in elongating stems by GC-MS showed that the initial substrate of GA 20-oxidase activity (GA53) accumulated, whereas the content of the major product (GA20) and of bioactive GA1 was lower in semidwarf compared with tall lines. We propose that the Os20ox2 gene corresponds to the sd-1 locus.

Sequence rearrangement in the AT-rich minisatellite of the novel rice transposable element Basho

In the process of characterizing a rice wx deletion mutant, an AT-rich minisatellite sequence that consisted of units of approximately 80 bp was detected about 2.3 kb downstream of the wx gene. This AT-rich minisatellite was a multiple-copy element (1 x 10(3) to 2 x 10(3) copies per haploid genome) and interspersed in the rice genome. By BLAST homology search it was indicated that not only the tandem repeat but also both flanking sequences were conserved among copies. According to the characteristics of the termini (5'-CHH ... CTAG-3') and a target site preference for T, this AT-rich minisatellite accompanying the flanking sequences was classified into a novel transposon, Basho. The results of direct amplification of Basho showed that relatively large variation in size existed in the Basho family. We estimate the variation to be generated by not only alteration of the number of units in the minisatellite but also by duplications of larger blocks including the conserved flanking sequences caused by single-strand mispairing (SSM) at noncontiguous repeats. Because the AT-rich minisatellite contained in Basho possessed several motifs of the matrix attachment region (MAR) in its repeat unit, the functional role as MAR in the rice genome was discussed.

Sequences upstream of the YRTA core region are essential for transcription of the tobacco atpB NEP promoter in chloroplasts in vivo

Transcription of the plastid atpB gene is accomplished by plastid-encoded (PEP) and nuclear-encoded (NEP) RNA polymerases. In contrast to NEP promoters of many other plastid genes, the tobacco atpB NEP promoter exhibits robust activity in chloroplasts in vivo. Previously, in vitro transcription assays using extracts from non-photosynthetic cells identified two elements required for full atpB NEP promoter activity, a core sequence and an upstream GAA box. Of these, only the core sequence containing the motif YRTA is conserved in the majority of NEP promoters. We used plastid transformation to examine the requirements for atpB NEP promoter activity in chloroplasts. Our results demonstrate that sequences upstream of the core element are essential for promoter activity in vivo, and that transcription of the NEP promoter in chloroplasts is not dependent on activity from an overlapping PEP promoter.

A variable minisatellite sequence in the chloroplast genome of Sorbus L. (Rosaceae: Maloideae)

The chloroplast genome is now known to be more variable than was once thought. Reports of RFLP (restriction fragment length polymorphism) and sequence variation, as well as variation in chloroplast microsatellites, are common. Here, data are presented on the variability of a minisatellite sequence in the chloroplast genome of Sorbus species. RFLP analysis of a PCR product comprising the region between the trnM and rbcL genes of nine Sorbus species identified seven size variants. Sequencing revealed the observed size polymorphism to be due to differences in the number of copies of an imperfect 9-bp motif. A more intensive survey of the variability of the minisatellite was undertaken in populations of Sorbus aucuparia. The potential uses of such regions in chloroplast DNA are discussed and a possible mechanism for the evolution of the minisatellite is presented.

Evolution of 5S rDNA unit arrays in the plant genus Nicotiana (Solanaceae)

Nicotiana tabacum (tobacco, Solanaceae) has two 5S ribosomal DNA (rDNA) families, one of unit length approximately 646 bp and the other -430 bp, that differ in the length of the 5S rDNA non-transcribed spacer (NTS). The long 5S rDNA family, found on the T genome of tobacco and in Nicotiana tomentosiformis, contains a GC-rich subregion that is absent in the short family. We designed primers for this subregion and generated a probe that we used against a range of Nicotiana and related Solanaceous species. We demonstrated the presence of the GC-rich subregion in a range of Nicotiana species, but it was absent in Nicotiana sylvestris, Nicotiana longiflora, and two closely related genera, Petunia and Solanum. These data suggest that this subregion of the NTS is likely to have evolved with the genus Nicotiana. The absence of the subregion in N. sylvestris and N. longiflora is likely to have arisen by a deletion event in the evolution of section alatae. We demonstrate patterns of evolution in the 5S rDNA unit cluster in relation to a phylogenetic reconstruction of species relationships in section tomentosae. Nicotiana glutinosa diverged early from the section and contains a 5S rDNA family based on a 550-bp unit. After this divergence, 430- and 650-bp rDNA unit families evolved. The 650-bp family is found in all species of tomentosae (except N. glutinosa) and in tobacco. The 430-bp family within tomentosae includes the GC-rich subregion and is thus unrelated to the 430-bp family in N. sylvestris. Nicotiana setchellii is unusual in that it has three 5S rDNA loci, including one locus that is exceptionally large. This species, unique to tomentosae, has a very abundant 900-bp unit family. It is possible that this 900-bp family occurs on this one large locus. In N. tomentosa and N. kawakamii, the 650-bp family is predominant, whereas N. tomentosiformis and N. otophora have only the 650-bp family. There is no clear relationship between the number of 5S families and the number of 5S rDNA loci. Certainly the replacement of 5S rDNA units, perhaps by gene conversion, has occurred repeatedly in the evolution of genus Nicotiana.

Molecular evolution at the self-incompatibility locus of Physalis longifolia (Solanaceae)

The simple codominant expression of the alleles at the self-incompatibility locus ( S) of Solanaceae and their extraordinary spectrum of sequence diversity present an ideal case for understanding the molecular forces that shape the sequence polymorphism of genes involved in recognition reactions. Here, with unprecedented conspecific sequences of 33 Physalis longifolia S alleles, including 27 of their single introns, the tempo-spatial nucleotide substitution patterns were first detailed for the S locus and analyzed in genealogical time scales. Three major genealogical clades of the 33 P. longifolia S alleles were used to divide the genealogical time into the within-clade and the between-clade periods. During the within-clade period, the average nonsynonymous substitution rate was 50% higher than the intron substitution rate but the opposite trend emerged for the between-clade period. A new and simple method developed here was utilized to estimate the selection intensity in the coding regions. The magnitudes and the distribution of these estimates, in conjunction with the spatial substitution pattern among closely related S sequences, revealed an initial short-term action of strong positive selection and a continuous but weaker action of negative selection brought by functional/structural constraints on the S alleles. The two modes of selection can significantly modify the branch lengths of an S genealogy and may be ubiquitous to recognition systems.

Identification of repetitive, genome-specific probes in crucifer oilseed species

Direct amplification of minisatellite DNA by PCR (DAMD PCR) was used to amplify and subsequently clone several fragments of DNA from crucifer species. The PCR-derived fragments of DNA were generated using known minisatellite core sequences as PCR primers. Southern hybridization of these putative minisatellite DNA fragments revealed that many were genome-specific; they hybridized with high affinity only to the genomic DNA of the species from which they were cloned. The DNA fragments were believed to be dispersed in the genome, based on smear-like hybridization signals on EcoRI-, BamHI-, and HindIII-digested genomic DNA. Genome-specific probes were specifically isolated from Brassica rapa (A genome), Brassica nigra (B genome), and Sinapis alba in addition to several other crucifer species. The sequence of a B. rapa specific probe (pBr17.1.3A) contained a minisatellite region that could be divided into three tandem repeats; each repeat contained between two and five subrepeats and each subrepeat shared a highly conserved core region of 29 bp. This minisatellite sequence also hybridized with high affinity to the A genome species B. napus and B. juncea. This research showed that dispersed, genome-specific probes can be isolated using DAMD PCR and that these probes could be used to detect and quantify alien DNA present in progeny from intergeneric or interspecific crosses.

The cloning of Ty1-copia-like retrotransposons from 10 varieties of banana (Musa Sp.)

Ty1-copia-like retrotransposons have been identified and investigated in several plant species. Here, the internal region of the reverse transcriptase (RT) gene of Ty1-copia-like retrotransposons was amplified by PCR from total genomic DNA of 10 varieties of banana. Two to four clones from each variety were sequenced. Extreme heterogeneity in the sequences of Ty1-copia-like retrotransposons from all the varieties was revealed following sequence analysis of the reverse transcriptase (RT) fragments. The size of the individual RT gene fragments varied between 213 and 309 bp. Southern blots of genomic DNA digested from Musa acuminata and other banana varieties probed with W8 clone from M. acuminata and A4 clone from Pisang Abu Nipah showed similar strong, multiple restriction fragments together with other faint hybridization band patterns with variable intensities indicating the presence of many copies of the Ty1-copia-like retrotransposons in the genomes. There was no correlation between retroelement sequence and the banana species (with A or B genomes) from which it arose, suggesting that the probes are not useful for tracking genomes through breeding populations.

Coding sequence divergence between two closely related plant species: Arabidopsis thaliana and Brassica rapa ssp. pekinensis

To characterize the coding-sequence divergence of closely related genomes, we compared DNA sequence divergence between sequences from a Brassica rapa ssp. pekinensis EST library isolated from flower buds and genomic sequences from Arabidopsis thaliana. The specific objectives were (i) to determine the distribution of and relationship between K(a) and K(s), (ii) to identify genes with the lowest and highest K(a): K(s) values, and (iii) to evaluate how codon usage has diverged between two closely related species. We found that the distribution of K(a): K(s) was unimodal, and that substitution rates were more variable at nonsynonymous than synonymous sites, and detected no evidence that K(a) and K(s) were positively correlated. Several genes had K(a): K(s) values equal to or near zero, as expected for genes that have evolved under strong selective constraint. In contrast, there were no genes with K(a): K(s) >1 and thus we found no strong evidence that any of the 218 sequences we analyzed have evolved in response to positive selection. We detected a stronger codon bias but a lower frequency of GC at synonymous sites in A. thaliana than B. rapa. Moreover, there has been a shift in the profile of most commonly used synonymous codons since these two species diverged from one another. This shift in codon usage may have been caused by stronger selection acting on codon usage or by a shift in the direction of mutational bias in the B. rapa phylogenetic lineage.

Laurasian migration explains Gondwanan disjunctions: evidence from Malpighiaceae

Explanations for biogeographic disjunctions involving South America and Africa typically invoke vicariance of western Gondwanan biotas or long distance dispersal. These hypotheses are problematical because many groups originated and diversified well after the last known connection between Africa and South America (approximately 105 million years ago), and it is unlikely that "sweepstakes" dispersal accounts for many of these disjunctions. Phylogenetic analyses of the angiosperm clade Malpighiaceae, combined with fossil evidence and molecular divergence-time estimates, suggest an alternative hypothesis to account for such distributions. We propose that Malpighiaceae originated in northern South America, and that members of several clades repeatedly migrated into North America and subsequently moved via North Atlantic land connections into the Old World during episodes starting in the Eocene, when climates supported tropical forests. This Laurasian migration route may explain many other extant lineages that exhibit western Gondwanan distributions.

Recombinant pronapin precursor produced in Pichia pastoris displays structural and immunologic equivalent properties to its mature product isolated from rapeseed

2S albumin storage proteins from rapeseed (Brassica napus), called napins, consist of two different polypeptide chains linked by disulphide bridges, which are derived by proteolytic cleavage from a single precursor. The precursor form of the napin BnIb (proBnIb) has been cloned using a PCR strategy and sequenced. The amino-acid sequence deduced from the clone includes 31 residues of the small chain and 75 of the large chain, which are connected by the peptide Ser-Glu-Asn. Expression of the cDNA encoding proBnIb has been carried out in the methylotrophic yeast Pichia pastoris. The induced protein was secreted to the extracellular medium at a yield of 80 mg.L(-1) of culture and was purified by means of size-exclusion chromatography and reverse phase-HPLC. Recombinant proBnIb appeared properly folded as its molecular and spectroscopic properties were equivalent to those of the mature heterodimeric protein. As 2S albumin storage proteins from Brassicaceae have been shown to be type I allergy inducers, the immunological activity of the recombinant proBnIb was analysed as a measure of its structural integrity. The immunological properties of the recombinant precursor and the natural napin were indistinguishable by immunoblotting and ELISA inhibition using polyclonal antisera and sera of patients allergic to mustard and rapeseed. In conclusion, the recombinant expression of napin precursors in P. pastoris has been shown to be a successful method for high yield production of homogeneous and properly folded proteins whose polymorphism and complex maturation process limited hitherto their availability.

Identification of a tolerant locus on Arabidopsis thaliana to hypervirulent beet curly top virus CFH strain

The infection of hosts by the geminivirus depends on the interactions between host and viral factors for viral DNA replication, viral gene expression, and the movement of virus throughout the hosts. This work reports that a hypervirulent strain of Beet curly top virus (BCTV) is different in its ability to infect several ecotypes of Arabidopsis thaliana. Symptoms appeared on Arabidopsis ecotypes around 7 to 10 d after inoculation with BCTV-CFH. Symptoms were more severe in ecotype SKKU including severe leaf curling and development of severely deformed and stunted boting compared to Col-O as a lab standard ecotype. One ecotype Cen-O was asymptomatic to BCTV-CFH infection. Studies of viral DNA replication and virus movement in three excised organs of asymptomatic Cen-O demonstrated that BCTV-CFH could replicate viral DNA and move systemically in this ecotype, suggesting that tolerance was due to the blocks of interactions between host and viral factors on symptom development. This asymptomatic phenotype is similar to the mutation of leftward ORFs, especially ORF R2. Genetic analysis of this ecotype Cen-O indicated that tolerance is due to a single recessive locus.

Assessing the level of collinearity between Arabidopsis thaliana and Brassica napus for A. thaliana chromosome 5

This study describes a comprehensive comparison of chromosome 5 of the model crucifer Arabidopsis with the genome of its amphidiploid crop relative Brassica napus and introduces the use of in silico sequence homology to identify conserved loci between the two species. A region of chromosome 5, spanning 8 Mb, was found in six highly conserved copies in the B. napus genome. A single inversion appeared to be the predominant rearrangement that had separated the two lineages leading to the formation of Arabidopsis chromosome 5 and its homologues in B. napus. The observed results could be explained by the fusion of three ancestral genomes with strong similarities to modern-day Arabidopsis to generate the constituent diploid genomes of B. napus. This supports the hypothesis that the diploid Brassica genomes evolved from a common hexaploid ancestor. Alignment of the genetic linkage map of B. napus with the genomic sequence of Arabidopsis indicated that for specific regions a genetic distance of 1 cM in B. napus was equivalent to 285 Kb of Arabidopsis DNA sequence. This analysis strongly supports the application of Arabidopsis as a tool in marker development, map-based gene cloning, and candidate gene identification for the larger genomes of Brassica crop species.

The 42-kDa coat protein of Andean potato mottle virus acts as a transcriptional activator in yeast

Interactions of viral proteins play an important role in the virus life cycle, especially in capsid assembly. Andean potato mottle comovirus (APMoV) is a plant RNA virus with a virion formed by two coat proteins (CP42 and CP22). Both APMoV coat protein open reading frames were cloned into pGBT9 and pGAD10, two-hybrid system vectors. HF7c yeast cells transformed with the p9CP42 construct grew on yeast dropout selection media lacking tryptophan and histidine. Clones also exhibited beta-galactosidase activity in both qualitative and quantitative assays. These results suggest that CP42 protein contains an amino acid motif able to activate transcription of His3 and lacZ reporter genes in Saccharomyces cerevisiae. Several deletions of the CP42 gene were cloned into the pGBT9 vector to locate the region involved in this activation. CP42 constructions lacking 12 residues from the C-terminal region and another one with 267 residues deleted from the N-terminus are still able to activate transcription of reporter genes. However, transcription activation was not observed with construction p9CP42deltaC57, which does not contain the last 57 amino acid residues. These results demonstrate that a transcription activation domain is present at the C-terminus of CP42 between residues 267 and 374.

Chromosome painting as a tool for rice genetics and breeding

Chromosome painting and genomic in situ hybridization (GISH) are both effective methods for basic genetic research and practical breeding. These methods were applied even in the typically small chromosomes of rice. This manuscript describes in detail, highly reproducible, complete protocols for chromosome painting and GISH in rice chromosomes. Examples of useful applications of these methods are also presented.

Fluorescent in situ hybridization in plant polytene chromosomes

Polytene chromosomes are found in specialized tissues, with high metabolic activity, of a few angiosperm genera. They differ from Diptera polytenics in several aspects, mainly because their chromatids on each chromosome are not tightly paired, nor are they so highly endoreplicated as those of Diptera. In situ hybridization with isotopic and non-isotopic probes has been successfully used in plant polytene chromosomes, mainly in Phaseolus coccineus and Vigna unguiculata, where they have been best investigated. The results reported for mitotic and polytene chromosomes of these species, and a few others, are compared aiming to ascertain the efficiency and limitations of FISH in plant polytenics. In general, polytene chromosomes either from embryo suspensor cells of P. coccineus or from anther tapetal cells of V. unguiculata proved to be quite a suitable system for localizing DNA sequences by FISH. The partially unsynapsed chromatids, typically found in plant polytenics, seem to be the most important hindrance for a precise chromosome mapping. On the other hand, the interphase polytene nucleus is a valuable system for localizing FISH signals since they conserve a spatial organization similar to that of mitotic interphase and produce much amplified signals.

TGTCACA motif is a novel cis-regulatory enhancer element involved in fruit-specific expression of the cucumisin gene

Cucumisin, a subtilisin-like serine protease, is expressed at high levels in the fruit of melon (Cucumis melo L.) and accumulates in the juice. We investigated roles of the promoter regions and DNA-protein interactions in fruit-specific expression of the cucumisin gene. In transient expression analysis, a chimeric gene construct containing a 1.2-kb cucumisin promoter fused to a beta-glucuronidase (GUS) reporter gene was expressed in fruit tissues at high levels, but the promoter activities in leaves and stems were very low. Deletion analysis indicated that a positive regulatory region is located between nucleotides -234 and -214 relative to the transcriptional initiation site. Gain-of-function experiments revealed that this 20-bp sequence conferred fruit specificity and contained a regulatory enhancer. Gel mobility shift experiments demonstrated the presence of fruit nuclear factors that interact with the cucumisin promoter. A typical G-box (GACACGTGTC) present in the 20-bp sequence did not bind fruit protein, but two possible cis-elements, an I-box-like sequence (AGATATGATAAAA) and an odd base palindromic TGTCACA motif, were identified in the promoter region between positions -254 and -215. The I-box-like sequence bound more tightly to fruit nuclear protein than the TGTCACA motif. The I-box-like sequence functions as a negative regulatory element, and the TGTCACA motif is a novel enhancer element necessary for fruit-specific expression of the cucumisin gene. Specific nucleotides responsible for the binding of fruit nuclear protein in these two elements were also determined.

The organization of Physcomitrella patensRAD51 genes is unique among eukaryotic organisms

Genetic recombination pathways and genes are well studied, but relatively little is known in plants, especially in lower plants. To study the recombination apparatus of a lower land plant, a recombination gene well characterized particularly in yeast, mouse, and man, the RAD51 gene, was isolated from the moss Physcomitrella patens and characterized. Two highly homologous RAD51 genes were found to be present. Duplicated RAD51 genes have been found thus far exclusively in eukaryotes with duplicated genomes. Therefore the presence of two highly homologous genes suggests a recent genome duplication event in the ancestry of Physcomitrella. Comparison of the protein sequences to Rad51 proteins from other organisms showed that both RAD51 genes originated within the group of plant Rad51 proteins. However, the two proteins form a separate clade in a phylogenetic tree of plant Rad51 proteins. In contrast to RAD51 genes from other multicellular eukaryotes, the Physcomitrella genes are not interrupted by introns. Because introns are a common feature of Physcomitrella genes, the lack of introns in the RAD51 genes is unusual and may indicate the presence of an unusual recombination apparatus in this organism. The presence of duplicated intronless RAD51 genes is unique among eukaryotes. Studies of further members of this lineage are needed to determine whether this feature may be typical of lower plants.

A 70-kDa chloroplast DNA polymerase from pea ( Pisum sativum) that shows high processivity and displays moderate fidelity

A 70-kDa chloroplast (ct) DNA polymerase from pea has been purified to apparent homogeneity. The ct DNA polymerase was insensitive to dideoxynucleotides (d(2) NTP) but showed high sensitivity to phosphonoacetic acid. The enzyme lacked any detectable 5'-->3' exonuclease activity but showed 3'-->5' exonuclease activity. The polymerase displayed high processivity (3 kb) and moderate fidelity, which may be sufficient for the faithful replication of the 140-kb pea ct genome. A 43-kDa accessory protein increased the polymerization rate but did not affect the rate of mis-incorporation in vitro, thus indicating that the domains for polymerisation and proof reading may be spatially separate.

Evolutionary analysis of S-RNase genes from Rosaceae species

Eight new cDNA sequences for S-RNases were cloned and analysed from almond (Prunus dulcis) cultivars of European origin, and compared to published sequences from other Rosaceae species. Insertions/deletions of 10-20 amino acid residues were detected in the RC4 and C5 domains of S-RNases from almond and sweet cherry. The S-RNases of the Prunus species and those of the genera Malus and Pyrus formed two distinct groups on phylogenetic analysis. Nucleotide substitutions were analysed in the S-RNase genes of these species. The S-genes of almond and sweet cherry have a lower Ka/Ks value than those of apple, pear and wild apple do. The fact that there is no fixed difference between the S-RNase genes of almond and sweet cherry, or between apple and pear, suggests that nucleotide substitutions only introduce transient polymorphism into the two groups, and rarely became fixed and contribute to divergence. Through the comparative study of 17 S-RNase genes from the genus Prunus and 18 from the genera Malus and Pyrus, some fixed nucleotide differences between the two groups were identified. These differences do not appear to be the result of selection for adaptive mutations, since the number of replacement substitutions is not significantly greater than the number of synonymous substitutions. S-RNase genes of almond and sweet cherry, and of apple and pear, showed little heterogeneity in nucleotide substitution rates. However, heterogeneity was observed between the two groups of S-alleles, with the Prunus alleles exhibiting a lower rate of non-synonymous substitutions than alleles from Malus and Pyrus. The evolutionary relationships between these species are discussed.

The domains required to direct core proteins of hepatitis C virus and GB virus-B to lipid droplets share common features with plant oleosin proteins

In mammalian tissue culture cells, the core protein of hepatitis C virus (HCV) is located at the surface of lipid droplets, which are cytoplasmic structures that store lipid. The critical amino acid sequences necessary for this localization are in a region of core protein that is absent in flavi- and pestiviruses, which are related to HCV. From our sequence comparisons, this region in HCV core was present in the corresponding protein of GBV-B, another virus whose genomic sequence has significant similarity to HCV. Expression of the putative GBV-B core protein revealed that it also was directed to lipid droplets. By extending the comparisons to cellular proteins, there were amino acid sequence similarities between the domains for lipid droplet association in HCV core and plant oleosin proteins. To determine whether these similarities were related functionally, an oleosin encoded by the Brassica napus bniii gene was expressed in different mammalian cell lines, where it retained the capacity to bind to lipid droplets. Analysis of deletion mutants indicated that the critical region within the protein required for this localization was the same for both plant and mammalian cells. A common feature in the viral and plant sequences was a motif containing proline residues. Mutagenesis of these residues in HCV core and plant oleosin abolished lipid droplet association. Finally, the domain within HCV core required for binding to lipid droplets could substitute for the equivalent domain in oleosin, further indicating the functional relatedness between the viral and plant sequences. These studies identify common features in disparate proteins that are required for lipid droplet localization.

Identification and characterization of RAPD markers inferring genetic relationships among Pine species

Total genomic DNAs were extracted from several populations of pine species and amplified using oligonucleotides of random sequences. Polymorphism in random amplified polymorphic DNA (RAPD) markers was high and sufficient in distinguishing each of the species. Genetic relationships among eight pine species (Pinus sylvestris, Pinus strobus, Pinus rigida, Pinus resinosa, Pinus nigra, Pinus contorta, Pinus monticola, and Pinus banksiana) from different provenances were analyzed. The degree of band sharing was used to evaluate genetic distance between species and to construct a phylogenetic tree. In general, the dendrogram corroborated the description of relationships based on morphological characteristics and crossability, but also provided new insights into pine taxonomy. RAPD markers specific to some pine species were cloned and sequenced. PCR amplifications using pairs of designed specific primers revealed that all the cloned sequences were likely genus specific because they were not found in spruce or larch. True species-specific sequences were identified using designed primers flanking cloned RAPD fragments. The analysis of RAPD fragment sequences confirmed the genetic relationships among species. A 2281-bp RAPD band called PI-Mt-Stb-23 from P. strobus was used as a probe in restriction fragment length polymorphism (RFLP) analysis and produced distinct banding patterns for each species examined, consistent with the highly polymorphic character of DNA-fingerprinting probes.