Microgeographic variation in recruitment under adult trees: arrival of new genotypes or perpetuation of the existing ones?

Investigating spatial variation in the relative importance of sexual reproduction and clonal propagation is critical to obtain more accurate estimates of future effective population sizes and genetic diversity, as well as to identify ecological correlates of clonality. We combined a stratified sampling scheme with microsatellite genetic analyses to estimate variation in the proportion of sexual versus clonal recruits among saplings in five populations of the tree Pyrus bourgaeana. Using a likelihood framework, we identified clones among the genotypes analysed and examined variation among populations regarding the proportion of saplings coming from clonal propagation. We also examined the relationship between the relative abundance of clonal shoots across the studied populations and their herbivory levels. Our results revealed that one third of the saplings examined (N = 225 saplings) had a probability above 0.9 of being clones of nearby (<10 m) trees, with the ratio between clonal propagation and sexual recruitment varying up to eight-fold among populations. A small portion of these putative clonal shoots reached sexual maturity. Relative abundance of clonal shoots did not significantly relate to the herbivory by ungulates. Our results call into question optimistic expectations of previous studies reporting sufficient levels of recruitment under parental trees without animal seed dispersal services. Nevertheless, given that some of these clonal shoots reach sexual maturity, clonal propagation can ultimately facilitate the long-term persistence of populations during adverse periods (e.g. environmental stress, impoverished pollinator communities, seed dispersal limitation).

Genetics at the verge of extinction: insights from the Iberian lynx.

Population viability might become compromised by the loss of genetic diversity and the accumulation of inbreeding resulting from population decline and fragmentation. The Iberian lynx (Lynx pardinus) provides a paradigmatic example of a species at the verge of extinction, and because of the well-documented and different demographic histories of the two remaining populations (Doñana and Andújar), it provides the opportunity to evaluate the performance of analytical methods commonly applied to recently declined populations. We used mitochondrial sequences and 36 microsatellite markers to evaluate the current genetic status of the species and to assess the genetic signatures of its past history. Mitochondrial diversity was extremely low with only two haplotypes, alternatively fixed in each population. Both remnant populations have low levels of genetic diversity at microsatellite markers, particularly the population from Doñana, and genetic differentiation between the two populations is high. Bayesian coalescent-based methods suggest an earlier decline starting hundreds of years ago, while heterozygosity excess and M-ratio tests did not provide conclusive and consistent evidence for recent bottlenecks. Also, a model of gene flow received overwhelming support over a model of pure drift. Results that are in conflict with the known recent demography of the species call for caution in the use of these methods, especially when no information on previous demographic history is available. Overall, our results suggest that current genetic patterns in the Iberian lynx are mainly the result of its recent decline and fragmentation and alerts on possible genetic risks for its persistence. Conservation strategies should explicitly consider this threat and incorporate an integrated genetic management of wild, captive and re-introduced populations, including genetic restoration through translocations.

Identification of neotropical felid faeces using RCP-PCR.

Faeces similarity among sympatric felid species has generally hampered their use in distributional, demographic and dietary studies. Here, we present a new and simple approach based on a set of species-specific primers, for the unambiguous identification of faeces from sympatric neotropical felids (i.e. puma, jaguar, jaguarundi and ocelot/ margay). This method, referred to as rapid classificatory protocol-PCR (RCP-PCR), consists of a single-tube multiplex PCR yielding species-specific banding patterns on agarose gel. The method was optimized with samples of known origin (14 blood and 15 fresh faeces) and validated in faecal samples of unknown origin (n = 138), for some of which (n = 40) we also obtained species identification based on mtDNA sequencing. This approach proved reliable and provides high identification success rates from faeces. Its simplicity and cost effectiveness should facilitate its application for routine surveys of presence and abundance of these species.

Searching for DNA in museum specimens: a comparison of sources in a mammal species.

The number of genetic studies that use preserved specimens as sources of DNA has been steadily increasing during the last few years. Therefore, selecting the sources that are more likely to provide a suitable amount of DNA of enough quality to be amplified and at the minimum cost to the original specimen is an important step for future research. We have compared different types of tissue (hides vs. bones) from museum specimens of Iberian lynx and multiple alternative sources within each type (skin, footpad, footpad powder, claw, diaphysis, maxilloturbinal bone, mastoid process and canine) for DNA yield and probability of amplification of both mitochondrial and nuclear targets. Our results show that bone samples yield more and better DNA than hides, particularly from sources from skull, such as mastoid process and canines. However, claws offer an amplification success as high as bone sources, which makes them the preferred DNA source when no skeletal pieces have been preserved. Most importantly, these recommended sources can be sampled incurring minimal damage to the specimens while amplifying at a high success rate for both mitochondrial and microsatellite markers.

Phylogeography of Southern Water Vole (Arvicola sapidus): evidence for refugia within the Iberian glacial refugium?

The role of Southern European peninsulas as glacial refugia for temperate species has been widely established, but phylogeographic patterns within refugia are being only recently addressed. Here we describe the phylogeographic patterns for Southern water vole (Arvicola sapidus) in its whole distribution across Iberia and France. Control region and cytochrome b sequences were obtained for 228 samples from 130 localities across Iberia and France. Eighty-five haplotypes were found in total yielding a high overall mitochondrial diversity (pi = 0.027; H = 0.974). Phylogeographic structure was relatively shallow (3.1% average intraspecific divergence) with few supported clades and 95% and 90% maximum parsimony unconnected networks, but significant, as reflected in increased pairwise nucleotide divergences with distance (r = 0.197, P = 0.03) and significant autocorrelation up to approximately 500 km. Spatial analysis of molecular variance analysis detected seven geographical groups explaining 43.73% of the total mitochondrial variation. We detected demographic expansions in three of these groups. A recent colonization of France from Iberia was suggested and estimated around 62 000 years bp by an isolation-with-migration model. Our results suggest the contribution of episodes of isolation in glacial subrefugia in Iberia, but seem to exclude a long-term isolation over successive glacial cycles. Phylogeographic divergence was probably tempered by relatively large population sizes and rapid and extensive mixing among subrefugia during interglacials, that might have eroded the phylogeographic structure accumulated at glacial peaks. Phenotypic differences in A. sapidus do not delineate historically isolated intraspecific divisions and do not warrant subspecific delimitations. Our results do support the existence of subrefugia within Iberia and their role in promoting intraspecific divergences.

Development and characterization of 11 microsatellite loci in a historically introduced carnivoran, the common genet (Genetta genetta).

Microsatellite markers were developed to assess population structure and patterns of translocation in the introduced European common genet (Genetta genetta). Primer pairs were designed for 60 microsatellite sequences enriched for CA, GA, CATC and TAGA repeat motifs. Eleven loci that proved to be polymorphic were genotyped in 33 individuals from southwestern France. The number of alleles per locus and observed heterozygosities varied from three to seven and from 0.2121 to 0.7576, respectively. One locus (B103) showed significant departure from Hardy-Weinberg equilibrium, probably due to the presence of null alleles. Tests of linkage disequilibrium did not detect significant associations among loci.

Contemporary pollen and seed dispersal in a Prunus mahaleb population: patterns in distance and direction.

Pollination and seed dispersal determine the spatial pattern of gene flow in plant populations and, for those species relying on pollinators and frugivores as dispersal vectors, animal activity plays a key role in determining this spatial pattern. For these plant species, reported dispersal patterns are dominated by short-distance movements with a significant amount of immigration. However, the contribution of seed and pollen to the overall contemporary gene immigration is still poorly documented for most plant populations. In this study we investigated pollination and seed dispersal at two spatial scales in a local population of Prunus mahaleb (L.), a species pollinated by insects and dispersed by frugivorous vertebrates. First, we dissected the relative contribution of pollen and seed dispersal to gene immigration from other parts of the metapopulation. We found high levels of gene immigration (18.50%), due to frequent long distance seed dispersal events. Second, we assessed the distance and directionality for pollen and seed dispersal events within the local population. Pollen and seed movement patterns were non-random, with skewed distance distributions: pollen tended moved up to 548 m along an axis approaching the N-S direction, and seeds were dispersed up to 990 m, frequently along the SW and SE axes. Animal-mediated dispersal contributed significantly towards gene immigration into the local population and had a markedly nonrandom pattern within the local population. Our data suggest that animals can impose distinct spatial signatures in contemporary gene flow, with the potential to induce significant genetic structure at a local level.

Differential contribution of frugivores to complex seed dispersal patterns.

Frugivores are highly variable in their contribution to fruit removal in plant populations. However, data are lacking on species-specific variation in two central aspects of seed dispersal, distance of dispersal and probability of dispersal among populations through long-distance transport. We used DNA-based genotyping techniques on Prunus mahaleb seeds dispersed by birds (small- and medium-sized passerines) and carnivorous mammals to infer each seed's source tree, dispersal distance, and the probability of having originated from outside the study population. Small passerines dispersed most seeds short distances (50% dispersed <51 m from source trees) and into covered microhabitats. Mammals and medium-sized birds dispersed seeds long distances (50% of mammals dispersed seeds >495 m, and 50% of medium-sized birds dispersed seeds to >110 m) and mostly into open microhabitats. Thus, dispersal distance and microhabitat of seed deposition were linked through the contrasting behaviors of different frugivores. When the quantitative contribution to fruit removal was accounted for, mammals were responsible for introducing two-thirds of the immigrant seeds into the population, whereas birds accounted for one-third. Our results demonstrate that frugivores differ widely in their effects on seed-mediated gene flow. Despite highly diverse coteries of mutualistic frugivores dispersing seeds, critical long-distance dispersal events might rely on a small subset of large species. Population declines of these key frugivore species may seriously impair seed-mediated gene flow in fragmented landscapes by truncating the long-distance events and collapsing seed arrival to a restricted subset of available microsites.

Population fragmentation leads to spatial and temporal genetic structure in the endangered Spanish imperial eagle.

The fragmentation of a population may have important consequences for population genetic diversity and structure due to the effects of genetic drift and reduced gene flow. We studied the genetic consequences of the fragmentation of the Spanish imperial eagle (Aquila adalberti) population into small patches through a temporal analysis. Thirty-four museum individuals representing the population predating the fragmentation were analysed for a 345-bp segment of the mitochondrial control region and a set of 10 nuclear microsatellite loci. Data from a previous study on the current population (N = 79) were re-analysed for this subset of 10 microsatellite markers and results compared to those obtained from the historical sample. Three shared mitochondrial haplotypes were found in both populations, although fluctuations in haplotype frequencies and the occurrence of a fourth haplotype in the historical population resulted in lower current levels of haplotype and nucleotide diversity. However, microsatellite markers revealed undiminished levels of nuclear diversity. No evidence for genetic structure was observed for the historical Spanish imperial eagle population, suggesting that the current pattern of structure is the direct consequence of population fragmentation. Temporal fluctuations in mitochondrial and microsatellite allelic frequencies were found between the historical and the current population as well as for each pairwise comparison between historical and current Centro and historical and current Parque Nacional de Doñana nuclei. Our results indicate an ancestral panmictic situation for the species that management policies should aim to restore. A historical analysis like the one taken here provides the baseline upon which the relative role of recent drift in shaping current genetic patterns in endangered species can be evaluated and this knowledge is used to guide conservation actions.

Mating patterns, pollen dispersal, and the ecological maternal neighbourhood in a Prunus mahaleb L. population.

Gender polymorphism, plant-animal interactions, and environmental heterogeneity are the three important sources of variation in mating system and pollen dispersal patterns. We used progeny arrays and paternity analysis to assess the effects of gender type and density level on variation in mating patterns within a highly isolated population of Prunus mahaleb, a gynodioecious species. All the adult trees in the population were sampled and located. The direct estimate of long-distance insect-mediated pollination events was low (< 10%). Gender expression deeply influenced the mating system, decreasing the outcrossing rates (t(m)) and the pollen pool diversity in hermaphrodite trees. Long intermate distances (> 250 m) were significantly more frequent among female mother trees. Variation in local tree density also affected pollen pool diversity and intermate distance, with a higher effective number of fathers (k(e)) and longer intermate distances for female trees in low-density patches. A canonical correlation analysis showed significant correlations between mating variables and the maternal ecological neighbourhood. Only the first canonical variable was significant and explained 78% of variation. Outcrossing rates tended to decrease, and the relatedness among the fathers tended to increase, when mother trees grew in dense patches with high cover of other woody species and taller vegetation away from the pine forest edge. We highlight the relevance of considering maternal ecological neighbourhood effects on mating system and gene flow studies as maternal trees act simultaneously as receptors of pollen and as sources of the seeds to be dispersed.

Organization and variation of the mitochondrial control region in two vulture species, Gypaetus barbatus and Neophron percnopterus.

We report the first entire mitochondrial DNA (mtDNA) control region sequences in two endangered vulture species, the bearded vulture (Gypaetus barbatus) and the Egyptian vulture (Neophron percnopterus). Results showed that the general organization of vulture control regions was very similar to other birds, with three distinct domains: a left variable domain (DI), a central conserved one (DII) including the F, E, D, and C boxes, and a right domain (DIII) containing the CSB1 sequence. However, due to the presence of long tandem repeats, vulture control regions differed from other avian control regions both in size and nucleotide composition. The Egyptian vulture control region was found to be the largest sequenced so far (2031 bp), due to the simultaneous presence of repeats in both DI (80 bp) and DIII (77 bp). Low variation was found in vulture control regions, particularly in G. barbatus, as the probable result of populations declines in the last few centuries.

Population genetics after fragmentation: the case of the endangered Spanish imperial eagle (Aquila adalberti).

The highly endangered Spanish imperial eagle, Aquila adalberti, has suffered from both population decline and fragmentation during the last century. Here we describe the current genetic status of the population using an extensive sampling of its current distribution range and both mitochondrial control region sequences and nuclear microsatellite markers. Results were evaluated in comparison to those obtained for the Eastern imperial eagle, Aquila heliaca, its nearest extant relative. Mitochondrial haplotype diversity was lower in the Spanish than in the Eastern species whereas microsatellite allelic richness and expected heterozygosity did not differ. Both allelic richness and expected heterozygosity were lower in the small Parque Nacional de Doñana breeding nucleus compared to the remaining nuclei. A signal for a recent genetic bottleneck was not detected in the current Spanish imperial eagle population. We obtained low but significant pairwise FST values that were congruent with a model of isolation by distance. FST and exact tests showed differentiation among the peripheral and small Parque Nacional de Doñana population and the remaining breeding subgroups. The centrally located Montes de Toledo population did not differ from the surrounding Centro, Extremadura and Sierra Morena populations whereas the latter were significantly differentiated. On the other hand, a Bayesian approach identified two groups, Parque Nacional de Doñana and the rest of breeding nuclei. Recent migration rates into and from Parque Nacional de Doñana and the rest of breeding nuclei were detected by assignment methods and estimated as 2.4 and 5.7 individuals per generation, respectively, by a Bayesian approach. We discuss how management strategies should aim at the maintenance of current genetic variability levels and the avoidance of inbreeding depression through the connection of the different nuclei.

Phylogenetic and phylogeographic analysis of Iberian lynx populations.

The Iberian lynx (Lynx pardinus), one of the world's most endangered cat species, is vulnerable due to habitat loss, increased fragmentation of populations, and precipitous demographic reductions. An understanding of Iberian lynx evolutionary history is necessary to develop rational management plans for the species. Our objectives were to assess Iberian lynx genetic diversity at three evolutionary timescales. First we analyzed mitochondrial DNA (mtDNA) sequence variation to position the Iberian lynx relative to other species of the genus LYNX: We then assessed the pattern of mtDNA variation of isolated populations across the Iberian Peninsula. Finally we estimated levels of gene flow between two of the most important remaining lynx populations (Doñana National Park and the Sierra Morena Mountains) and characterized the extent of microsatellite locus variation in these populations. Phylogenetic analyses of 1613 bp of mtDNA sequence variation supports the hypothesis that the Iberian lynx, Eurasian lynx, and Canadian lynx diverged within a short time period around 1.53-1.68 million years ago, and that the Iberian lynx and Eurasian lynx are sister taxa. Relative to most other felid species, genetic variation in mtDNA genes and nuclear microsatellites were reduced in Iberian lynx, suggesting that they experienced a fairly severe demographic bottleneck. In addition, the effects of more recent reductions in gene flow and population size are being manifested in local patterns of molecular genetic variation. These data, combined with recent studies modeling the viability of Iberian lynx populations, should provide greater urgency for the development and implementation of rational in situ and ex situ conservation plans.

Activation of Wnt signaling rescues neurodegeneration and behavioral impairments induced by beta-amyloid fibrils.

Alzheimer's disease (AD) is a progressive neurodegenerative disorder, which is probably caused by the cytotoxic effect of the amyloid beta-peptide (Abeta). We report here molecular changes induced by Abeta, both in neuronal cells in culture and in rats injected in the dorsal hippocampus with preformed Abeta fibrils, as an in vivo model of the disease. Results indicate that in both systems, Abeta neurotoxicity resulted in the destabilization of endogenous levels of beta-catenin, a key transducer of the Wnt signaling pathway. Lithium chloride, which mimics Wnt signaling by inhibiting glycogen synthase kinase-3beta promoted the survival of post-mitotic neurons against Abeta neurotoxicity and recovered cytosolic beta-catenin to control levels. Moreover, the neurotoxic effect of Abeta fibrils was also modulated with protein kinase C agonists/inhibitors and reversed with conditioned medium containing the Wnt-3a ligand. We also examined the spatial memory performance of rats injected with preformed Abeta fibrils in the Morris water maze paradigm, and found that chronic lithium treatment protected neurodegeneration by rescuing beta-catenin levels and improved the deficit in spatial learning induced by Abeta. Our results are consistent with the idea that Abeta-dependent neurotoxicity induces a loss of function of Wnt signaling components and indicate that lithium or compounds that mimic this signaling cascade may be putative candidates for therapeutic intervention in Alzheimer's patients.

Faecal genetic analysis to determine the presence and distribution of elusive carnivores: design and feasibility for the Iberian lynx.

Noninvasive methods using genetic markers have been suggested as ways to overcome difficulties associated with documenting the presence of elusive species. We present and assess a novel, reliable and effective molecular genetic technique for the unequivocal genetic identification of faeces from the endangered Iberian lynx (Lynx pardinus). From mitochondrial DNA (mtDNA) cytochrome b and D-loop region sequences, we designed four species-specific primers (for products 130-161 bp long) that were considered to be likely to amplify degraded DNA. We compared two DNA extraction methods, various DNA amplification conditions and the robustness and specificity of the primer pairs with 87 lynx samples from 5 potentially different lynx populations and with 328 samples of other carnivore species. The utility of the identification technique was tested with faeces of different ages, with faeces from controlled field experiments, and with faeces collected from locales with possible lynx populations from throughout the state of Andalusia, Spain (8052 km2). Faecal mtDNA extraction was more efficient using PBS wash of the faeces instead of a faeces homogenate. Our assay increased from 92.6 to 99% efficiency with a second amplification and a reduction in template concentration to overcome polymerase chain reaction (PCR) inhibition. Our assay never produced false positives, and correctly identified all lynx faeces. Of 252 faeces samples of unknown species collected throughout Andalusia, 26.6% (from three different areas) were classified as Iberian lynx, 1.4% showed evidence of PCR inhibition and 1.2% were of uncertain origin. This method has proven to be a reliable technique that can be incorporated into large-scale surveys of Iberian lynx populations and exemplifies an approach that can easily be extended to other species.

Seed dispersal by animals: exact identification of source trees with endocarp DNA microsatellites.

A long-standing challenge in studies of seed dispersal by animal frugivores has been the characterization of the spatial relationships between dispersed seeds and the maternal plants, i.e. the seed shadow. The difficulties to track unambiguously the origin of frugivore-dispersed seeds in natural communities has been considered an unavoidable limitation of the research field and precluded a robust analysis of the direct consequences of zoochory. Here we report that the multilocus genotype at simple sequence repeat (SSR; microsatellite) loci of the woody endocarp, a tissue of maternal origin, provides an unequivocal genetic fingerprint of the source tree. By comparing the endocarp genotype against the complete set of genotypes of reproductive trees in the population, we could unambiguously identify the source tree for 82.1% of the seeds collected in seed traps and hypothesize that the remaining 17.9% of sampled seeds come from other populations. Identification of the source tree for Prunus mahaleb seeds dispersed by frugivores revealed a marked heterogeneity in the genetic composition of the seed rain in different microhabitats, with a range of 1-5 distinct maternal trees contributing seeds to a particular landscape patch. Within-population dispersal distances ranged between 0 and 316 m, with up to 62% of the seeds delivered within 15 m of the source trees. Long distance dispersal events, detected by the exclusion of all reproductive trees in the population, accounted for up to 17.9% of the seeds sampled. Our results indicate strong distance limitation of seed delivery combined with infrequent long-distance dispersal events, extreme heterogeneity in the landscape pattern of genetic makeup, and a marked mosaic of multiple parentage for the seeds delivered to a particular patch.

RAPD variation and population genetic structure in Prunus mahaleb (Rosaceae), an animal-dispersed tree.

We examined the patterns of random amplified polymorphic DNA (RAPD) variation among seven Prunus mahaleb (Rosaceae) populations extending over approximately 100 km2 to examine local differentiation in relation to spatial isolation due to both geographical distance and differences in elevation. No less than 51. 4% of the RAPD loci were polymorphic, but very few were fixed and among-population variation accounted for 16.46% of variation in RAPD patterns. Mean gene diversity was 0.1441, with mean Nei's genetic diversity for individual populations ranging between 0.089 and 0.149. Mean GST value across loci was 0.1935 (range, 0.0162-0.4685), giving an average estimate for Nm of 1.191. These results suggest extensive gene flow among populations, but higher GST and lower Nm values relative to other outcrossing, woody species with endozoochorous dispersal, also suggest a process of isolation by distance. The combined effect of both geographical and elevation distances and nonoverlapping flowering and fruiting phenophases on the GST matrix was partially significant, revealing only marginal isolation of the P. mahaleb populations. The matrix correlation between estimated Nm values among populations and the geographical + elevation distance matrices (r = -0.4623, P = 0.07), suggests a marginal trend for more isolated populations to exchange less immigrants. Long-distance seed dispersal by efficient medium-sized frugivorous birds and mammals is most likely associated to the high levels of within-population genetic diversity. However, vicariance factors and demographic bottlenecks (high postdispersal seed and seedling mortality) explain comparatively high levels of local differentiation.

Structure of the dehydrin tas14 gene of tomato and its developmental and environmental regulation in transgenic tobacco.

We have isolated a genomic clone encoding tomato TAS14, a dehydrin that accumulates in response to mannitol, NaCl or abscisic acid (ABA) treatment. A fragment of tas14 gene containing the region from -2591 to +162 fused to beta-glucuronidase gene drives ABA- and osmotic stress-induced GUS expression in transgenic tobacco. Histochemical analysis of salt-, mannitol- and ABA-treated plants showed GUS activity mainly localized to vascular tissues, outer cortex and adventitious root meristems, coinciding with the previously observed distribution of TAS14 protein in salt-stressed tomato plants. In addition, GUS activity was also observed in guard cells, trichomes and leaf axils. Developmentally regulated gus expression was studied in unstressed plants and found to occur not only in embryos, but also in flowers and pollen. Tas14 expression in floral organs was confirmed by northern blots of tomato flowers.

Regulation of type I plasminogen activator inhibitor by fibrin degradation products in rat lung fibroblasts.

Persistent fibrin deposition in tissues characterizes the early pathology of many types of injury. In an animal model of bleomycin-induced lung fibrosis, increased expression of type 1 plasminogen activator inhibitor (PAI-1) is associated with accumulation of fibrin in fibroproliferative lesions. Plasmin proteolysis of cross-linked fibrin generates fibrin degradation products (FDPs) with multiple biological activities in several cell types. We reasoned that fibrin fragments may also regulate fibroblast-mediated fibrinolysis. In this study, we describe induction of PAI-1 mRNA, protein, and activity by soluble FDPs and fibrinogen in rat lung fibroblast monolayers. FDPs are more potent than fibrinogen, inducing a concentration-dependent, maximal 3.7 (+/- 0.9)-fold increase in PAI-1 mRNA as measured by northern blotting and a 9.0 (+/- 1.3)-fold induction of PAI-1 antigen levels. Active PAI-1 is demonstrated in fibrinogen- and FDP-stimulated conditioned media. Further characterization of this response shows that PAI-1 expression is induced by the DD/D fragments, but not by immunopurified fragment E. Experiments using Actinomycin D and puromycin indicate that the induction appears to be transcriptionally regulated and is not dependent on new protein synthesis. FDP induction of PAI-1 suggests a matrix-cell feedback process in which a fibrin fragment modulates expression of an important regulator of fibrinolysis.

Cellular and molecular events in a newly organizing lateral root meristem.

Spontaneous or auxin-induced lateral root formation in radish and Arabidopis provides an efficient system in which to examine molecular and cellular events associated with the initiation of a new meristem. Subtracted cDNA libraries made at different times in lateral root initiation were used as a source of genes that are expressed differentially during this developmental process, and expression studies on a small gene family of ribosomal protein genes were conducted. From analysis of cell division patterns in pericycle cells the number of founder cells for lateral roots was established. By the use of in vitro growth assays lateral root formation was determined to be a two-stage process. First a primordium is formed, and subsequently a subset of primordial cells begins to function as the lateral root apical meristem. This mode of root development has implications for pattern formation in newly organizing organs.