Robust Estimates of Divergence Times and Selection with a Poisson Random Field Model: A Case Study of Comparative Phylogeographic Data

Mutation frequencies can be modeled as a Poisson random field (PRF) to estimate speciation times and the degree of selection on newly arisen mutations. This approach provides a quantitative theory for comparing intraspecific polymorphism with interspecific divergence in the presence of selection and can be used to estimate population genetic parameters. Although the original PRF model has been extended to more general biological settings to make statistical inference about selection and divergence among model organisms, it has not been incorporated into phylogeographic studies that focus on estimating population genetic parameters for nonmodel organisms. Here, we modified a recently developed time-dependent PRF model to independently estimate genetic parameters from a nuclear and mitochondrial DNA data set of 22 sister pairs of birds that have diverged across a biogeographic barrier. We found that species that inhabit humid habitats had more recent divergence times and larger effective population sizes than those that inhabit drier habitats, and divergence time estimated from the PRF model were similar to estimates from a coalescent species-tree approach. Selection coefficients were higher in sister pairs that inhabited drier habitats than in those in humid habitats, but overall the mitochondrial DNA was under weak selection. Our study indicates that PRF models are useful for estimating various population genetic parameters and serve as a framework for incorporating estimates of selection into comparative phylogeographic studies.     

The Great American Interchange in birds: a phylogenetic perspective with the genus Trogon

The 'Great American Interchange' (GAI) is recognized as having had a dramatic effect on biodiversity throughout the Neotropics. However, investigation of patterns in Neotropical avian biodiversity has generally been focused on South American taxa in the Amazon Basin, leaving the contribution of Central American taxa under-studied. More rigorous studies of lineages distributed across the entire Neotropics are needed to uncover phylogeographical patterns throughout the area, offering insights into mechanisms that contribute to overall Neotropical biodiversity. Here we use mitochondrial DNA sequence data and intensive geographical sampling from the widespread Neotropical avian genus Trogon to investigate the role of the GAI in shaping its phylogeographical history. Our results show that genetic diversity in Trogon exceeds the perceived biodiversity, and that the GAI resulted in lineage diversification within the genus. Despite greater diversity in South America, a Central American centre of origin with multiple and independent dispersals into South America is indicated. These dispersals were followed by the evolution of divergent lineages associated with the Andes Mountains and other South American geographical features. According to our phylogenetic reconstructions, several species, which were originally defined by morphological characters, are nonmonophyletic. In sum, our results elucidate the evolutionary history of Trogon , reveal patterns obscured by extant biodiversity, and serve as a biogeographical model to consider in future studies.     

Defining reproductively isolated units in a cryptic and syntopic species complex using mitochondrial and nuclear markers: the brooding brittle star, Amphipholis squamata (Ophiuroidea)

At a time when biodiversity is threatened, we are still discovering new species, and particularly in the marine realm. Delimiting species boundaries is the first step to get a precise idea of diversity. For sympatric species which are morphologically undistinguishable, using a combination of independent molecular markers is a necessary step to define separate species. Amphipholis squamata , a cosmopolitan brittle star, includes several very divergent mitochondrial lineages. These lineages appear totally intermixed in the field and studies on morphology and colour polymorphism failed to find any diagnostic character. Therefore, these mitochondrial lineages may be totally interbreeding presently. To test this hypothesis, we characterized the genetic structure of the complex in the French Mediterranean coast using sequences of mitochondrial DNA (16S) and for the first time, several nuclear DNA markers (introns and microsatellites). The data revealed six phylogenetic lineages corresponding to at least four biological species. These sibling species seem to live in syntopy. However, they seem to display contrasted levels of genetic diversity, suggesting they have distinct demographic histories and/or life-history traits. Genetic differentiation and isolation-by-distance within the French Mediterranean coasts are revealed in three lineages, as expected for a species without a free larval phase. Finally, although recombinant nuclear genotypes are common within mitochondrial lineages, the data set displays a total lack of heterozygotes, suggesting a very high selfing rate, a feature likely to have favoured the formation of the species complex.     

Genetic differentiation in the genus Lithops L. (Ruschioideae, Aizoaceae) reveals a high level of convergent evolution and reflects geographic distribution

Southern Africa is one of the hot spots for plant biodiversity, with ca. 80% of species endemic to this area. Rapid and recent radiations in Southern African plant genera were triggered by fine-scale differences in climate, topography and geology. The genus Lithops (Ruschioideae, Aizoaceae) contains 37 species and is widely distributed in Southern Africa. Species delimitation within the genus is challenging because the limited number of morphological characters in these reduced succulents varies intensely between populations, presumably as adaptations to local geological environments. We analysed phylogenetic relationships within Lithops using non-coding chloroplast DNA (trnS-trnG intergenic spacer), nuclear ribosomal internal transcribed spacer (nrITS) sequences and AFLP data. Genetic variability of the sequence data was very low, but AFLP data detected nine clades within Lithops that do not fit current morphology-based taxonomy. Two of these clades are separated by their distribution on the northern and eastern border of the distribution area, and four clades are found in the Gariep Centre in the estuary of the Orange River. Morphological similarities, especially colour of leaves, evolved repeatedly within the clades, thus we hypothesise that closely related species became adapted to different soil types in a mosaic-like geological environment. One-third of the species are found in the Gariep Centre, characterised by extremely diverse edaphic habitats.     

Mechanisms of radiation in a bat group from the genus Pipistrellus inferred by phylogeography, demography and population genetics

Here, we present a study of the Pipistrellus pipistrellus species complex, a highly diversified bat group with a radiation centre in the Mediterranean biodiversity hotspot. The study sample comprised 583 animals from 118 localities representatively covering the bats' range in the western Palearctic. We used fast-evolving markers (the mitochondrial D-loop sequence and 11 nuclear microsatellites) to describe the phylogeography, demography and population structure of this model taxon and address details of its diversification. The overall pattern within this group includes a mosaic of phylogenetically basal, often morphologically distant, relatively small and mostly allopatric demes in the Mediterranean Basin, as well as two sympatric sibling species in the large continental part of the range. The southern populations exhibit constant size, whereas northern populations show a demographic trend of growth associated with range expansion during the Pleistocene climate oscillations. There is evidence of isolation by distance and female philopatry in P. pipistrellus sensu stricto. Although the northern populations are reproductively isolated, we detected introgression events among several Mediterranean lineages. This pattern implies incomplete establishment of reproductive isolating mechanisms in these populations as well as the existence of a past reinforcement stage in the continental siblings. The occurrence of reticulations in the radiation centre among morphologically and ecologically derived relict demes suggests that adaptive unequal gene exchange within hybridizing populations could play a role in speciation and adaptive radiation within this group.     

Transitions between self-compatibility and self-incompatibility and the evolution of reproductive isolation in the large and diverse tropical genus Dendrobium (Orchidaceae)

Background and Aims The evolution of interspecific reproductive barriers is crucial to understanding species evolution. This study examines the contribution of transitions between self-compatibility (SC) and self-incompatibility (SI) and genetic divergence in the evolution of reproductive barriers in Dendrobium, one of the largest orchid genera. Specifically, it investigates the evolution of pre- and postzygotic isolation and the effects of transitions between compatibility states on interspecific reproductive isolation within the genus.Methods The role of SC and SI changes in reproductive compatibility among species was examined using fruit set and seed viability data available in the literature from 86 species and ∼2500 hand pollinations. The evolution of SC and SI in Dendrobium species was investigated within a phylogenetic framework using internal transcribed spacer sequences available in GenBank.Key Results Based on data from crossing experiments, estimations of genetic distance and the results of a literature survey, it was found that changes in SC and SI significantly influenced the compatibility between species in interspecific crosses. The number of fruits produced was significantly higher in crosses in which self-incompatible species acted as pollen donor for self-compatible species, following the SI × SC rule. Maximum likelihood and Bayesian tests did not reject transitions from SI to SC and from SC to SI across the Dendrobium phylogeny. In addition, postzygotic isolation (embryo mortality) was found to evolve gradually with genetic divergence, in agreement with previous results observed for other plant species, including orchids.Conclusions Transitions between SC and SI and the gradual accumulation of genetic incompatibilities affecting postzygotic isolation are important mechanisms preventing gene flow among Dendrobium species, and may constitute important evolutionary processes contributing to the high levels of species diversity in this tropical orchid group.     

Exploring neutral and adaptive processes in expanding populations of gilthead sea bream, Sparus aurata L., in the North-East Atlantic

Recent studies in empirical population genetics have highlighted the importance of taking into account both neutral and adaptive genetic variation in characterizing microevolutionary dynamics. Here, we explore the genetic population structure and the footprints of selection in four populations of the warm-temperate coastal fish, the gilthead sea bream (Sparus aurata), whose recent northward expansion has been linked to climate change. Samples were collected at four Atlantic locations, including Spain, Portugal, France and the South of Ireland, and genetically assayed using a suite of species-specific markers, including 15 putatively neutral microsatellites and 23 expressed sequence tag-linked markers, as well as a portion of the mitochondrial DNA (mtDNA) control region. Two of the putatively neutral markers, Bld-10 and Ad-10, bore signatures of strong directional selection, particularly in the newly established Irish population, although the potential 'surfing effect' of rare alleles at the edge of the expansion front was also considered. Analyses after the removal of these loci suggest low but significant population structure likely affected by some degree of gene flow counteracting random genetic drift. No signal of historic divergence was detected at mtDNA. BLAST searches conducted with all 38 markers used failed to identify specific genomic regions associated to adaptive functions. However, the availability of genomic resources for this commercially valuable species is rapidly increasing, bringing us closer to the understanding of the interplay between selective and neutral evolutionary forces, shaping population divergence of an expanding species in a heterogeneous milieu.     

Out of sight, out of mind: high cryptic diversity obscures the identities and histories of geminate species in the marine bivalve subgenus Acar

Aim The rise of the Isthmus of Panama and the formation of ‘geminate’ species pairs serves as an important model of allopatric speciation. However, to function as a model system, hypothesized geminates must first be shown to be each other’s closest living relatives. If the presence of cryptic taxa obscures true relationships, the biogeographical histories of transisthmian taxa are likely to be misinterpreted. We have therefore completed a phylogeographic survey of the transisthmian bivalve subgenus Acar in the genus Barbatia to characterize patterns of tropical American diversity and to place transisthmian taxa in a regional phylogeographic context. Location Tropical America. Methods Mitochondrial cytochrome c oxidase I (COI) and nuclear internal transcribed spacer (ITS) sequences were obtained from 233 specimens of Acar. Sequences were analysed using cladistic and Bayesian methods. Divergence times between species were inferred from net nucleotide divergences and a coalescence‐based method. Results The survey revealed 22 COI clades that were also monophyletic at ITS, indicating that the taxonomy of Acar is potentially greater than a fivefold diversity underestimate. The lone previously recognized geminate [Barbatia (Acar) gradata and Barbatia (Acar) domingensis] is composed of 15 clades. Among the four transisthmian lineages identified, two diverged more than 14 Ma; the two other geminates split just prior to the time of final seaway closure. In addition to a fourfold increase in the number of known geminates, our data show that within‐basin diversification has been more impressive, with one geminate splitting into five monophyletic clades in the Western Atlantic alone since seaway closure. Electron microscopy of the larval shells of specimens indicates that the transisthmian lineage with the greatest rate of post‐Isthmian diversification possesses non‐planktonic larvae, a life‐history feature linked to high speciation rates. Main conclusions Our analyses revealed that the identities of geminate pairs split by the Isthmus of Panama were obscured by extremely high tropical American cryptic diversity. Although we have identified four geminates, only two appear to have been split by the Isthmus. Our uncovering of extensive post‐Isthmian diversification is consistent with the palaeontological perspective that the final closure of the Central American Seaway was followed by high rates of subgeneric diversification, particularly in the tropical Western Atlantic.     

Two different satellite DNAs in Beta vulgaris L.: evolution,quantification and distribution in the genus

Summary Two highly repeated EcoRI (0.45 × 10⁶) and BamHI (0.17 × 10⁶) fragments per haploid genome were found in sugar beet genomic DNA. Both fragments were located by 6% acrylamide-gel electrophoresis, purified and cloned in pUC18. Four of the inserts corresponding to each family were chosen for further study. Both fragment families display the main characteristics of the satellite DNA of animals and plants. The EcoRI and BamHI fragment families are arranged in long tandem arrays. Fragments of the EcoRI family (pBVE) were analyzed. They vary both in sequence and in length (158–160 nt) in comparison with the consensus sequence of 159nt. Both families are A-T rich; pBVE is 59% rich while pBVB is 69% rich. The BVESAT family is present in all the members of the section Vulgares. It is conserved in the section Procumbentes with 80% homology and the same length, but is not detectable in the Corollinae. The sequence variation rate and the variation in length (330±5 nt) are of the same order in comparison with those of the BVESAT family. However, the BVBSAT family is present in species of the section Vulgares only. As regards other plant satellite DNAs, the BVESAT family shares homology with Allium cepa satellite DNA, with three of the yeast centromeric sequences, and with three Arabidopsis thaliana sequences. The BVBSAT family is unique to the Vulgares and does not share any homology with other plant or animal satellite DNAs sequences so far.     

Dispersal, recruitment and migratory behaviour in a hawksbill sea turtle aggregation

We investigated the dispersal, recruitment and migratory behaviour of the hawksbill sea turtle ( Eretmochelys imbricata ), among different life-history stages and demographic segments of the large hawksbill turtle aggregation at Mona Island, Puerto Rico. There were significant differences in both mitochondrial DNA (mtDNA) haplotype diversity and haplotype frequencies among the adult males, females and juveniles examined, but little evidence for temporal heterogeneity within these same groups sampled across years. Consistent with previous studies and the hypothesis of strong natal homing, there were striking mtDNA haplotype differences between nesting females on Mona Island and nesting females in other major Caribbean rookeries. Breeding males also showed strong, albeit weaker, genetic evidence of natal homing. Overall, Bayesian mixed-stock analysis suggests that Mona Island was the natal rookery for 79% (65–94%) of males in the aggregation. In contrast, the Mona Island rookery accounted for only a small subset of the new juvenile recruits to the foraging grounds or in the population of older juvenile hawksbills turtles on Mona. Instead, both new recruits and the older juvenile hawksbill turtles on Mona more likely recruited from other Caribbean rookeries, suggesting that a mechanism besides natal homing must be influencing recruitment to feeding habitats. The difference in the apparent degree of natal homing behaviour among the different life-history stages of hawksbill turtles at Mona Island underscores the complexity of the species' life-history dynamics and highlights the need for both local and regional conservation efforts.     

The distribution, organization and evolution of two abundant and widespread repetitive DNA sequences in the genus Hordeum

The genomic organization and chromosomal distributions of two abundant tandemly repeated DNA sequences, dpTa1 and pSc119.2, were examined in six wild Hordeum taxa, representing the four basic genomes of the genus, by Southern and fluorescence in situ hybridization. The dpTa1 probe hybridized to between 30 and 60 sites on the chromosomes of all five diploid species studied, but hybridization patterns differed among the species. Hybridization of the pSc119.2 sequence to the chromosomes and Southern blots of digested DNA detected signals in Hordeum bulbosum, Hordeum chilense, Hordeum marinum and Hordeum murinum 4x, but not in Hordeum murinum 2x and Hordeum vulgare ssp. spontaneum. A maximum of one pSc119.2 signal was observed in the terminal or subterminal region of each chromosome arm in the species carrying this sequence. The species carrying the same I-genome differed in the presence (Hordeum bulbosum) or absence (Hordeum spontaneum) of pSc119.2. The presence of pSc119.2 in the tetraploid cytotype of Hordeum murinum, but its absence in the diploid cytotype, suggests that the tetraploid is not likely to be a simple autotetraploid of the diploid. Data about the inter- and intra-specific variation of the two independent repetitive DNA sequences give information about both the interrelationships of the species and the evolution of the repetitive sequences. Received: 17 March 1999 / Accepted: 16 June 1999     

Range‐wide genetic homogeneity in the California sea mussel (Mytilus californianus): a comparison of allozymes,nuclear DNA markers,and mitochondrial DNA sequences

We tested for genetic differentiation among six populations of California sea mussels (Mytilus californianus) sampled across 4000 km of its geographical range by comparing patterns of variation at four independent types of genetic markers: allozymes, single‐copy nuclear DNA markers, and DNA sequences from the male and female mitochondrial genomes. Despite our extensive sampling and genotyping efforts, we detected no significant differences among localities and no signal of isolation by distance suggesting that M. californianus is genetically homogeneous throughout its range. This concordance differs from similar studies on other mytilids, especially in the role of postsettlement selection generating differences between exposed coastal and estuarine habitats. To assess if this homogeneity was due to M. californianus not inhabiting estuarine environments, we reviewed studies comparing allozymes with other classes of nuclear DNA markers. Although both types of markers gave broadly consistent results, there was a bias favouring studies in which allozymes were more divergent than DNA markers (nine to three) and a disproportionate number of these cases involved marine taxa (seven). Furthermore, allozymes were significantly more heterogeneous than DNA markers in three of the four studies that sampled coastal and estuarine habitats. We conclude that the genetic uniformity exhibited by M. californianus may result from a combination of extensive gene flow and the lack of exposure to strong selective gradients across its range.     

Historical introgression and the role of selective vs. neutral processes in structuring nuclear genetic variation (AFLP) in a circumpolar marine fish, the capelin (Mallotus villosus)

The capelin (Mallotus villosus) is a widespread marine fish species for which previous work has identified geographically distinct mtDNA clades, the frontiers of which are well within adult and larval dispersal capabilities. Here, we use AFLPs to test for the presence of nuclear gene flow among clades. In addition, we evaluate genetic structuring within one clade, the Northwest Atlantic (NWA). We found that each of the mtDNA clades corresponds with a unique nuclear DNA genetic cluster. Within the NWA clade, we detected individuals with small but significant amounts of genetic ancestry from other clades, likely due to historical introgression. Further support for historical introgression comes from analyses of variance in locus-specific differentiation, which support introgression between some clades and divergence without gene flow between others. Within the NWA, we identified two genetic clusters that correspond to sites in geographically adjacent areas. However, these clusters differ primarily at 'outlier' loci, and a genetic subdivision (K=2) was not supported by genetic clustering programs using neutral loci. Significant neutral F(ST) differentiation was found only between sites that otherwise differed at outlier loci. Thus, these populations may be in the initial stages of 'isolation by adaptation'. These results suggest strong between-clade reproductive isolation despite opportunities for gene flow and support the hypothesis that selection can contribute to divergence in otherwise 'open' systems.     

The evolutionary origin of a novel karyotype in Timarcha (Coleoptera, Chrysomelidae) and general trends of chromosome evolution in the genus

In this work, we have analysed the karyotypes of six species of Timarcha for the first time and updated the cytological information for two additional taxa, for one of them confirming previous results ( Timarcha erosa vermiculata ), but not for the other ( T. scabripennis ). We describe the remarkable karyotype of T. aurichalcea , the lowest chromosome number in the genus (2 n  = 18), distinctive as well for the presence of an unusual chiasmatic sexual bivalent hitherto unreported for Timarcha . This study increases the number of species studied cytologically in this genus to forty. Additional cytogenetic analyses are performed on several species, including Ag-NOR staining and fluorescent in situ hybridization (FISH) studies with ribosomal DNA probes. Karyotype evolution is analysed by tracing different karyotype coding strategies on a published independent phylogenetic hypothesis for Timarcha based on the study of three genetic markers. The implementation of a likelihood model of character change optimized onto the phylogeny is tentatively used to detect possible drifts in chromosome changes. These analyses show that karyotype is conservative in the evolution of the genus and that there is an apparent trend to reducing chromosome number. Cytological and phylogenetic data are used to explain the evolutionary origin of the karyotype of T. aurichalcea by two centric fusions involving one pair of acrocentric autosomes and the sexual chromosomes.     

Demographic history,geographical distribution and reproductive isolation of distinct lineages of blue rockfish (Sebastes mystinus), a marine fish with a high dispersal potential

Understanding the barriers to genetic exchange in taxonomic groups that have a high dispersal potential will provide critical information on speciation in general. Blue rockfish (Sebastes mystinus) are good taxa to examine speciation because they are nonmigratory inhabitants of shallow rocky reef habitats along the eastern North Pacific with a pelagic larval stage lasting 3–5 months. The goal of this study was to analyse the evolutionary history and distribution patterns of different lineages within S. mystinus described previously and use this information to understand the speciation process in this group of high dispersal fish. The molecular data derived from specimens sampled over approximately 1650 km of the S. mystinus range revealed a northerly and southerly distribution for the two lineages. Almost equal frequencies of both lineages occurred at centrally located sample locations, with evidence of reproductive isolation between the lineages. A demographic analysis showed that the two lineages diverged and experienced sudden expansion prior to the last glacial maximum, which affected the observed pattern of genetic structure. The spatial distribution, demographic history and degree of genetic distinctiveness found from the genetic analysis, despite the high potential for dispersal in S. mystinus, suggest both lineages diverged in allopatry.     

Reticulate evolution and phenotypic diversity in North American ciscoes, Coregonus ssp. (Teleostei: Salmonidae): implications for the conservation of an evolutionary legacy

Coregonine fishes are notorious taxonomicproblems due to their extreme morphological andecological variation. In North America, diversity is particularly baffling among ciscoes, and both morphological and phylogenetic analyses have resulted in major polytomy among the 8 taxa of the ``Coregonus artedi' species complex. Ciscoes arealso a devastated group, accounting for 10% ofthe fish species listed by the Committee on theStatus of Endangered Wildlife in Canada. Here,we complete the genetic characterization ofNorth American ciscoes with mitochondrial andmicrosatellites markers previously used toanalyse populations of C. artedi in orderto elucidate the evolutionary history andidentify appropriate conservation units. Ourresults revealed a complex evolutionary historymarked by postglacial reticulation eventscoupled with recent and independent evolutionof similar phenotypes (taxa). Genetic variationreflects geography rather than taxonomy, andconsequently, we recommend that a single taxon,C. artedi (sensu lato) be recognized.Local genetic differentiation is often coupledwith ecophenotypic diversification, and gillraker polymorphisms, depth-related habitatpreference and reproductive behaviour areconsidered as phenotypic traits with probableadaptive value contributing to the nicheexpansion of ciscoes. Ecomorphotypes of eachparticular locale thus represent a uniqueexpression of a diverse genetic pool stillundergoing divergence and sorting.Consequently, ciscoes from lakes with distinctecomorphotypes are recognized as ESUs, as wellas each of sympatric forms when they aregenetically differentiated. We recommend thatan ESU strategy focusing at a very local levelbe adopted for continental ciscoes as a validalternative to protect significant evolutionaryprocesses of divergence encountered inpolytypic species of newly colonized habitats.     

Seasonal fluctuations in the population of denitrifying and N2-fixing bacteria in an acid soil of a Norway spruce forest

The seasonal fluctuations in the concentration of cultured denitrifying and N₂-fixing bacteria were followed in an ammonium fertilised and a control soil of a Norway spruce forest near Villingen/Black Forest from December 1994 to August 1998. The horizontal distribution of bacteria in three layers was determined by the MPN-method and by molecular probing (colony hybridisation) using specific 0.4–0.7 kb DNA probes for denitrification steps (narG, nirS, nirK and nosZ) and for N₂-fixation (nifH). The data showed that highest bacterial counts and higher numbers of denitrifying and N₂-fixing bacteria were generally detectable in the upper (= 5 cm) soil layer and that their amount decreased with soil depth. The concentration of these cultured bacteria showed seasonal fluctuations with highest numbers in autumn/winter/early spring and with low counts in summer. Denitrifying and N₂-fixing bacteria amounted to less than 10% of the total number of cultured bacteria determined by the MPN-method. Fertilisation with ammonium did not cause a shift in the population of these bacteria. These findings were corroborated by hybridisation experiments with genomic DNA isolated from the different layers. Strongest DNA–DNA hybridisation band intensities were obtained in the upper soil layer and their intensities decreased with soil depth. Soil samples from Villingen assayed in the laboratory produced N₂O (in dependence of nitrate and C₂H₂ added to the vessels) and utilised this gas with higher activities in the assays with the fertilised soil. It is concluded that molecular techniques can successfully be applied for assessing seasonal fluctuations of bacterial populations in soil. Relative abundance of denitrifying and N₂-fixing bacteria can be determined from experiments with DNA isolated from soils. Attempts to transform these results to the total population of soil bacteria on a single cell basis are faced with many uncertainties.     

Cloning,expression and characterization of a gene which encodes a topoisomerase I with positive supercoiling activity in pea

We have isolated and sequenced the full length cDNA for topoisomerase I. Using degenerate primers, based on the conserved amino acid sequences of five eukaryotic topoisomerase I, a 386 bp fragment was PCR amplified using pea cDNA as template. This fragment was used as a probe to screen a pea cDNA library. Two partial cDNA clones were isolated which were truncated at the 5 end. RACE-PCR was employed to isolate the remaining portion of the gene. The total size of the gene was 3055 bp with an open reading frame of 2676 bp. The deduced structure of pea topoisomerase I contain 892 amino acids with a calculated molecular weight of 100 kDa and an estimated pI of 9.3. A comparison of the deduced amino acid sequences of the pea topo I with the other eukaryotic topoisomerases clearly suggested that they are all related. Pea topoisomerase I has been overexpressed in E. coli system and the recombinant topoisomerase purified to homogeneity. The purified protein relaxes both positive and negative supercoiled DNA in the absence of divalent cation Mg²⁺. In the presence of Mg²⁺ ions the purified enzyme introduces positive supercoils a unique property not reported in any other organism except in archaebacterial topoisomerase I. Polyclonal antibodies were raised against recombinant topoisomerase I and western blotting with sub-cellular fractions indicated the localization of this topoisomerase in pea nuclei.