Long range physical map of the Wilms' tumor - aniridia region on human chromosome 11

Nephroblastoma or Wilms' tumor is a pediatric renal malignancy that is the most frequently occurring childhood solid tumor. Approximately 1-2% of children with Wilms' tumor also present with aniridia, a congenital absence of all or part of the iris of the eye. These children also have high rates of genitourinary anomalies and mental retardation resulting in what is called the WAGR (Wilms' tumor, aniridia, genitourinary anomaly, mental retardation) syndrome. Cytogenetic analysis of metaphase chromosomes from these patients revealed a consistent deletion of band P13 on chromosome 11. These observations suggest close physical linkage between the disease-related loci, and further imply that development of each phenotype results from the loss of normal gene function.^ The objective of this work is to understand the molecular events at chromosome band 11p13 that are essential to the development of sporadic Wilms' tumor and sporadic aniridia. Two human/hamster somatic cell hybrids have been used to identify sixteen independent DNA probes that map to this segment of the human genome. These newly identified DNA probes and four previously reported probes (CAT, FSHB, D11S16, and HBVIS) have been used to subdivide 11p13 into five intervals defined by overlapping constitutional deletions from several WAGR patients. A long-range physical map of 11p13 has been constructed using each of these probes in Southern blot analysis of genomic DNA after digestion with infrequently cutting restriction enzymes and pulse-field gel electrophoresis. This map, established primarily with MluI and NotI, spans approximately 13 $\times$ 10$\sp{6}$ bp and encompasses deletion and translocation breakpoints associated with genitourinary anomalies, aniridia, and sporadic Wilms' tumor. This complete physical map of human chromosome band 11p13 enables us to localize the genes for sporadic Wilms' tumor and sporadic aniridia to a small number of specific NotI fragments. ^

The physical map of {\it Escherichia coli\/} K-12 strain MG1655: Construction, comparison with other strains and implications for mechanisms controlling oligonucleotide frequency

Complete NotI, SfiI, XbaI and BlnI cleavage maps of Escherichia coli K-12 strain MG1655 were constructed. Techniques used included: CHEF pulsed field gel electrophoresis; transposon mutagenesis; fragment hybridization to the ordered $\lambda$ library of Kohara et al.; fragment and cosmid hybridization to Southern blots; correlation of fragments and cleavage sites with EcoMap, a sequence-modified version of the genomic restriction map of Kohara et al.; and correlation of cleavage sites with DNA sequence databases. In all, 105 restriction sites were mapped and correlated with the EcoMap coordinate system.^ NotI, SfiI, XbaI and BlnI restriction patterns of five commonly used E. coli K-12 strains were compared to those of MG1655. The variability between strains, some of which are separated by numerous steps of mutagenic treatment, is readily detectable by pulsed-field gel electrophoresis. A model is presented to account for the difference between the strains on the basis of simple insertions, deletions, and in one case an inversion. Insertions and deletions ranged in size from 1 kb to 86 kb. Several of the larger features have previously been characterized and some of the smaller rearrangements can potentially account for previously reported genetic features of these strains.^ Some aspects of the frequency and distribution of NotI, SfiI, XbaI and BlnI cleavage sites were analyzed using a method based on Markov chain theory. Overlaps of Dam and Dcm methylase sites with XbaI and SfiI cleavage sites were examined. The one XbaI-Dam overlap in the database is in accord with the expected frequency of this overlap. The occurrence of certain types of SfiI-Dcm overlaps are overrepresented. Of the four subtypes of SfiI-Dcm overlap, only one has a partial inhibitory effect on the activity of SfiI. Recognition sites for all four enzymes are rarer than expected based on oligonucleotide frequency data, with this effect being much stronger for XbaI and BlnI than for NotI and SfiI. The latter two enzyme sites are rare mainly due to apparent negative selection against GGCC (both) and CGGCCG (NotI). The former two enzyme sites are rare mainly due to effects of the VSP repair system on certain di-tri- and tetranucleotides, most notably CTAG. Models are proposed to explain several of the anomalies of oligonucleotide distribution in E. coli, and the biological significance of the systems that produce these anomalies is discussed. ^

Mapping and cloning of genes from portions of the human genome using somatic cell hybrids

Human x rodent somatic cell hybrids have played an important role in human genetics research. They have been especially useful for assigning genes to chromosomes and isolating DNA markers from specific regions of the human genome.^ By employing a combination of somatic cell genetic, recombinant DNA, and cytogenetic techniques, human DNA excision repair gene ERCC4 was mapped regionally to human 16p13.13-13.2, even though the gene has not been cloned. Human x Chinese hamster ovary (CHO) cell hybrids selected for human ERCC4 activity and containing 16p13.1-p13.3 as the only human genetic material were identified. These hybrids were used to order DNA markers located in 16p13.1-p13.3. New DNA markers physically close to ERCC4 were isolated from such hybrids. Using amplified human DNA from the hybrids as probe in fluorescent in situ hybridization, the short arm breakpoint in the chromosome 16 inversion associated with acute myelomonocytic leukemia (AMML) was found to be physically close to the ERCC4 gene. The physical mapping and eventually, the cloning of the ERCC4 gene, will benefit the understanding of the DNA repair system and the study of other important biomedical problems such as tumorigenesis.^ To facilitate the cloning of ERCC4 gene and, in general, the cloning of genes from any defined regions of the human genome, a method was developed for the direct isolation of human transcribed genes ffom somatic cell hybrids. cDNA was prepared from human x rodent hybrid by using consensus 5$\sp\prime$ splice site sequences as primers. These primers were designed to select immature, unspliced messenger RNA (still retaining species specific repeat sequences) as templates. Screening of a derived cDNA library for human repeat sequences resulted in the isolation of human clones at the anticipated frequency with characteristics expected of exons of transcribed human genes. The usefulness of the splice site specific primers was analyzed and the cDNA synthesis conditions with these primers were optimized. The procedure was shown to be sensitive enough to clone weakly expressed genes. Studying the expression of the represented genes with the isolated clones was shown to be feasible. Such regional specific human gene fragments will be very valuable for many human genetic studies such as the search of inherited disease genes and the construction of a cDNA map of the human genome. ^

Characterization of the beta-lactamase gene from {\it Enterococcus faecalis\/} strain HH22

The plasmid-encoded, constitutively produced $\beta$-lactamase gene from Enterococcus faecalis strain HH22 was genetically characterized. A restriction endonuclease map of the 5.1 kb EcoRI fragment encoding the enterococcal $\beta$-lactamase was prepared and compared with the restriction map of a cloned staphylococcal $\beta$-lactamase gene (from the naturally-occurring staphylococcal $\beta$-lactamase plasmid pI258). Comparison and hybridization studies showed that there were identical restriction sites in the region of the $\beta$-lactamase structural gene but not in the region surrounding this gene. Also the enterococcal $\beta$-lactamase plasmid did not encode resistance to mercury or cadmium which is encoded by the small, transducible staphylococcal $\beta$-lactamase plasmids. The nucleotide sequence of the enterococcal gene was shown to be identical to the published sequences of three of four staphylococcal type A $\beta$-lactamase genes; more differences were seen with the genes for staphylococcal type C and D enzymes. One hundred-forty nucleotides upstream of the $\beta$-lactamase start codon were also determined for the inducible staphylococcal $\beta$-lactamase gene on pI258; this sequence was identical to that of the constitutively expressed enterococcal gene indicating that the changes resulting in constitutive expression are not due to changes in the promoter or operator region. Moreover, complementation studies indicated that production of the enterococcal enzyme could be repressed. The gene for the enterococcal $\beta$-lactamase and an inducible staphylococcal $\beta$-lactamase were each cloned into a shuttle vector and then transformed into enterococcal and staphylococcal recipients. The major difference between the two host backgrounds was that more enzyme was produced by the staphylococcal host, regardless of the source of the gene but no qualitative difference was seen between the two genera. Also a difference in the level of resistance to ampicillin was seen between the two backgrounds with the cloned enzymes by MIC and time-kill studies. The location of the enzyme was found to be host dependent since each cloned gene generated extracellular (free) enzyme in the staphylococcus and cell bound enzyme in the enterococcus. Based on the identity of the enterococcal $\beta$-lactamase and several staphylococcal $\beta$-lactamases, these data suggest recent spread of $\beta$-lactamase to enterococci and also suggest loss of a functional repressor. ^

Control of chromosomal rearrangements in {\it Escherichia coli\/}

A complete physical map of Escherichia coli K-12 strain MG1655 was constructed by digesting chromosomal DNA with the infrequently cutting restriction enzymes NotI, SfiI and XbaI and separating the fragments by pulsed field gel electrophoresis. The map was used to compare six K-12 strains of E. coli. Although several differences were noted and localized, the map of MG1655 was representative of all the K-12 strains tested. The maps were also used to analyze chromosomal rearrangements in the E. coli strain MG1655. The spontaneous and UV induced frequencies of tandem duplication formation were measured at several loci distributed around the chromosome. The spontaneous duplication frequency varied from 10$\sp{-5}$ to 10$\sp{-3}$ and increased at least ten-fold following mild UV irradiation treatment. Duplications of several regions of the chromosome, including the serA region and the metE region, were mapped using pulsed field gel electrophoresis. Duplications of serA were found to be large, ranging in size from 600 kb to 2100 kb. Several of the duplications isolated at serA were caused by ectopic recombination between IS5 elements and between IS186 elements. Duplications of the metE region, however, were almost exclusively the result of ectopic recombination between ribosomal RNA cistrons. Duplication frequencies were determined at both serA and metE in wild type and mismatch repair mutant strains (mutL, mutS, uvrD and recF). Even though all of the mismatch repair mutations increased duplication frequency of metE, the largest increases were observed in the mutL and mutS strains. Duplication frequency of serA was increased less dramatically by mutations in mismatch repair. Several duplications of metE isolated in a wild type and a mismatch repair mutant were mapped. The results showed that the same repeated sequences were used for duplication formation in the mismatch repair mutant as were used in the wild type strain. Several isolates showed evidence of multiple rearrangements indicating that mismatch repair may play a role in stabilizing the genome by controlling chromosomal rearrangement. ^

The story of RP10: A positional cloning approach to the identification of an autosomal dominant retinitis pigmentosa gene

Retinitis pigmentosa (RP) is an inherited retinal degenerative disease that is the leading cause of inherited blindness worldwide. Characteristic features of the disease include night blindness, progressive loss of visual fields, and deposition of pigment on the retina in a bone spicule-like pattern. RP is marked by extreme genetic heterogeneity with at least 19 autosomal dominant, autosomal recessive and X-linked loci identified. RP10, which maps to chromosome 7q, was the fifth autosomal dominant RP locus identified, and accounts for the early-onset disease in two independent families. Extensive linkage and haplotype analyses have been performed in these two families which have allowed the assignment of the disease locus to a 5-cM region on chromosome 7q31.3. In collaboration with Dr. Eric Green (National Center for Human Genome Research, National Institutes of Health), a well-characterized physical map of the region was constructed which includes YAC, BAC and cosmid coverage. The entire RP10 critical region resides within a 9-Mb well-characterized YAC contig. These physical maps not only provided the resources to undertake the CAIGES (cDNA amplification for identification of genomic expressed sequences) procedure for identification of retinal candidate genes within the critical region, but also identified a number of candidate genes, including transducin-$\gamma$ and blue cone pigment genes. All candidate genes examined were excluded. In addition, a number of ESTs were mapped within the critical region. EST20241, which was isolated from an eye library, corresponded to the 3$\sp\prime$ region of the ADP-ribosylation factor (ARF) 5 gene. ARF5, with its role in vesicle transport and possible participation in the regulation of the visual transduction pathway, became an extremely interesting candidate gene. Using a primer walking approach, the entire 3.2 kb genomic sequence of the ARF5 gene was generated and developed intronic primers to screen for coding region mutations in affected family members. No mutations were found in the ARF5 gene, however, a number of additional ESTs have been mapped to the critical region, and, as the large-scale sequencing projects get underway, megabases of raw sequence data from the RP10 region are becoming available. These resources will hasten the isolation and characterization of the RP10 gene. ^

CHARACTERIZATION OF THE COLONIZATION FACTOR ANTIGEN II PLASMID (CFA/II) FROM ENTEROTOXIGENIC ESCHERICHIA COLI

The etiological role of enterotoxigenic E. coli (ETEC) in diarrheal diseases of man and domestic animals is firmly established. Besides the production of enterotoxins (ST and LT), ETEC produces other important virulence factors; the colonization factor antigens (CFAs). CFAs mediate the attachment of ETEC to the epithelial cells of the small intestine, and this favors colonization by the bacteria and facilitates delivery of the enterotoxins to the intestinal cells.^ The production of enterotoxin and CFA is determined by plasmids and has been found to be restricted to a select number of E. coli serotypes.^ In this work, plasmid DNA analysis was performed in twenty-three CFA/II-producing enterotoxigenic Escherichia coli strains and their spontaneous CFA/II-negative derivatives. In some cases, strains lost the high molecular weight plasmid and also the ability to produce CFA/II, ST and LT. In other cases there was a deletion of the plasmid, which produced strains that were CFA/II('-), ST('-), LT('-) or CFA/II('-), ST('+), LT('+).^ The CFA/II plasmid from strain PB-176 (06:H16:CFA/II('+), ST('+), LT('+)) was transferred by transformation into E. coli K12 with concomitant transfer of the three characteristics: CFA/II, ST and LT.^ A physical map of the prototype CFA/II:ST:LT (pMEP60) plasmid was constructed by restriction endonuclease analysis and compared to plasmids from three other CFA/II-producing strains. A CFA/II-negative (but ST and LT positive) deletion derivative of pMEP60 (pMEP30) was also included in the map. The four CFA/II plasmids analyzed had a common region of approximately 30 kilobase pairs. The toxin genes were approximately 5 kbp apart and about 20 kbp from the common region. The information given by this physical map could be of great value when constructing a clone that will express the CFA/II genes but not the toxin genes. ^

GENETIC ORGANIZATION OF THE STRUCTURAL GENES FOR NITRATE REDUCTASE (TN5)

Nitrate reductase in Escherichia coli is a membrane-bound anaerobic enzyme that is repressed by oxygen and induced by nitrate. The genetic organization of the structural genes for the two larger subunits of nitrate reductase ((alpha) and (beta)) was determined by immunoprecipitation analysis of the formation of these proteins in nitrate reductase-deficient mutants resulting from transposon Tn5 mutagenesis. The results suggested that the genes encoding the (alpha) and (beta) subunits (narG and H) were arranged in an operon with transcription in the direction promoter(--->)(alpha)(--->)(beta). Segments of the chromosome containing the Tn5 inserts from several of the mutants were cloned into plasmid pBR322 and the positions of the transposons determined by restriction mapping. The Tn5 insertion sites were localized on two contiguous EcoRI fragments spanning about 6.6 kilobases of DNA. The narI gene (proposed to encode the (gamma) subunit) was positioned immediately downstream from the (beta)-gene (narH) by Southern analysis of Tn10 insertions into the narI locus. A Tn10 insertion into the narK locus, proposed to encode a nitrate-sensitive repressor of other anaerobic enzymes, was located about 1.5 kilobases upstream from the narGHI operon promoter. The narL locus, proposed to encode a nitrate-sensitive positive regulator of the narGHI operon and known to be genetically linked to the other nar genes, was demonstrated to lie outside a 19.3-kilobase region of the chromosome which encompasses the other nar genes. The physical limit of the narGHI promoter was defined by studying the effect of Tn5 insertions into a hybrid plasmid containing the functional operon. The points of origin of the coding regions for the (alpha) and (beta) genes were deduced by alignment of the chromosomal map of Tn5 insertion sites with the sizes of (alpha) and (beta) subunit fragments produced by plasmids carrying these Tn5 inserts in the nar operon. The coding region for the (alpha) subunit (143,000 daltons) begins about 250 nucleotides downstream from the deduced limit of the promoter region and includes about 4.0 kilobases of DNA; the region encoding (beta) (60,000 daltons) lies immediately downstream from the (alpha)-gene and is approximately 1.6 kilobases in length. The adjacent region encoding the (gamma) subunit (19,000 daltons) is approximately 0.5 kilobase in length. ^

Theoretical and experimental studies of linkage disequilibrium in human populations

Linkage disequilibrium (LD) is defined as the nonrandom association of alleles at two or more loci in a population and may be a useful tool in a diverse array of applications including disease gene mapping, elucidating the demographic history of populations, and testing hypotheses of human evolution. However, the successful application of LD-based approaches to pertinent genetic questions is hampered by a lack of understanding about the forces that mediate the genome-wide distribution of LD within and between human populations. Delineating the genomic patterns of LD is a complex task that will require interdisciplinary research that transcends traditional scientific boundaries. The research presented in this dissertation is predicated upon the need for interdisciplinary studies and both theoretical and experimental projects were pursued. In the theoretical studies, I have investigated the effect of genotyping errors and SNP identification strategies on estimates of LD. The primary importance of these two chapters is that they provide important insights and guidance for the design of future empirical LD studies. Furthermore, I analyzed the allele frequency distribution of 26,530 single nucleotide polymorphisms (SNPs) in three populations and generated the first-generation natural selection map of the human genome, which will be an important resource for explaining and understanding genomic patterns of LD. Finally, in the experimental study, I describe a novel and simple, low-cost, and high-throughput SNP genotyping method. The theoretical analyses and experimental tools developed in this dissertation will facilitate a more complete understanding of patterns of LD in human populations. ^

Modeling Techniques for the High-Resolution Interpretation of Cryo-Electron Microscopy Reconstructions

Essential biological processes are governed by organized, dynamic interactions between multiple biomolecular systems. Complexes are thus formed to enable the biological function and get dissembled as the process is completed. Examples of such processes include the translation of the messenger RNA into protein by the ribosome, the folding of proteins by chaperonins or the entry of viruses in host cells. Understanding these fundamental processes by characterizing the molecular mechanisms that enable then, would allow the (better) design of therapies and drugs. Such molecular mechanisms may be revealed trough the structural elucidation of the biomolecular assemblies at the core of these processes. Various experimental techniques may be applied to investigate the molecular architecture of biomolecular assemblies. High-resolution techniques, such as X-ray crystallography, may solve the atomic structure of the system, but are typically constrained to biomolecules of reduced flexibility and dimensions. In particular, X-ray crystallography requires the sample to form a three dimensional (3D) crystal lattice which is technically di‑cult, if not impossible, to obtain, especially for large, dynamic systems. Often these techniques solve the structure of the different constituent components within the assembly, but encounter difficulties when investigating the entire system. On the other hand, imaging techniques, such as cryo-electron microscopy (cryo-EM), are able to depict large systems in near-native environment, without requiring the formation of crystals. The structures solved by cryo-EM cover a wide range of resolutions, from very low level of detail where only the overall shape of the system is visible, to high-resolution that approach, but not yet reach, atomic level of detail. In this dissertation, several modeling methods are introduced to either integrate cryo-EM datasets with structural data from X-ray crystallography, or to directly interpret the cryo-EM reconstruction. Such computational techniques were developed with the goal of creating an atomic model for the cryo-EM data. The low-resolution reconstructions lack the level of detail to permit a direct atomic interpretation, i.e. one cannot reliably locate the atoms or amino-acid residues within the structure obtained by cryo-EM. Thereby one needs to consider additional information, for example, structural data from other sources such as X-ray crystallography, in order to enable such a high-resolution interpretation. Modeling techniques are thus developed to integrate the structural data from the different biophysical sources, examples including the work described in the manuscript I and II of this dissertation. At intermediate and high-resolution, cryo-EM reconstructions depict consistent 3D folds such as tubular features which in general correspond to alpha-helices. Such features can be annotated and later on used to build the atomic model of the system, see manuscript III as alternative. Three manuscripts are presented as part of the PhD dissertation, each introducing a computational technique that facilitates the interpretation of cryo-EM reconstructions. The first manuscript is an application paper that describes a heuristics to generate the atomic model for the protein envelope of the Rift Valley fever virus. The second manuscript introduces the evolutionary tabu search strategies to enable the integration of multiple component atomic structures with the cryo-EM map of their assembly. Finally, the third manuscript develops further the latter technique and apply it to annotate consistent 3D patterns in intermediate-resolution cryo-EM reconstructions. The first manuscript, titled An assembly model for Rift Valley fever virus, was submitted for publication in the Journal of Molecular Biology. The cryo-EM structure of the Rift Valley fever virus was previously solved at 27Å-resolution by Dr. Freiberg and collaborators. Such reconstruction shows the overall shape of the virus envelope, yet the reduced level of detail prevents the direct atomic interpretation. High-resolution structures are not yet available for the entire virus nor for the two different component glycoproteins that form its envelope. However, homology models may be generated for these glycoproteins based on similar structures that are available at atomic resolutions. The manuscript presents the steps required to identify an atomic model of the entire virus envelope, based on the low-resolution cryo-EM map of the envelope and the homology models of the two glycoproteins. Starting with the results of the exhaustive search to place the two glycoproteins, the model is built iterative by running multiple multi-body refinements to hierarchically generate models for the different regions of the envelope. The generated atomic model is supported by prior knowledge regarding virus biology and contains valuable information about the molecular architecture of the system. It provides the basis for further investigations seeking to reveal different processes in which the virus is involved such as assembly or fusion. The second manuscript was recently published in the of Journal of Structural Biology (doi:10.1016/j.jsb.2009.12.028) under the title Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions. This manuscript introduces the evolutionary tabu search strategies applied to enable a multi-body registration. This technique is a hybrid approach that combines a genetic algorithm with a tabu search strategy to promote the proper exploration of the high-dimensional search space. Similar to the Rift Valley fever virus, it is common that the structure of a large multi-component assembly is available at low-resolution from cryo-EM, while high-resolution structures are solved for the different components but lack for the entire system. Evolutionary tabu search strategies enable the building of an atomic model for the entire system by considering simultaneously the different components. Such registration indirectly introduces spatial constrains as all components need to be placed within the assembly, enabling the proper docked in the low-resolution map of the entire assembly. Along with the method description, the manuscript covers the validation, presenting the benefit of the technique in both synthetic and experimental test cases. Such approach successfully docked multiple components up to resolutions of 40Å. The third manuscript is entitled Evolutionary Bidirectional Expansion for the Annotation of Alpha Helices in Electron Cryo-Microscopy Reconstructions and was submitted for publication in the Journal of Structural Biology. The modeling approach described in this manuscript applies the evolutionary tabu search strategies in combination with the bidirectional expansion to annotate secondary structure elements in intermediate resolution cryo-EM reconstructions. In particular, secondary structure elements such as alpha helices show consistent patterns in cryo-EM data, and are visible as rod-like patterns of high density. The evolutionary tabu search strategy is applied to identify the placement of the different alpha helices, while the bidirectional expansion characterizes their length and curvature. The manuscript presents the validation of the approach at resolutions ranging between 6 and 14Å, a level of detail where alpha helices are visible. Up to resolution of 12 Å, the method measures sensitivities between 70-100% as estimated in experimental test cases, i.e. 70-100% of the alpha-helices were correctly predicted in an automatic manner in the experimental data. The three manuscripts presented in this PhD dissertation cover different computation methods for the integration and interpretation of cryo-EM reconstructions. The methods were developed in the molecular modeling software Sculptor (http://sculptor.biomachina.org) and are available for the scientific community interested in the multi-resolution modeling of cryo-EM data. The work spans a wide range of resolution covering multi-body refinement and registration at low-resolution along with annotation of consistent patterns at high-resolution. Such methods are essential for the modeling of cryo-EM data, and may be applied in other fields where similar spatial problems are encountered, such as medical imaging.

Learning -dependent plasticity of movement representations in the primate primary motor cortex

Primary motor cortex (M1) is involved in the production of voluntary movement and contains a complete functional representation, or map, of the skeletal musculature. This functional map can be altered by pathological experiences, such as peripheral nerve injury or stroke, by pharmacological manipulation, and by behavioral experience. The process by which experience-dependent alterations of cortical function occur is termed plasticity. In this thesis, plasticity of M1 functional organization as a consequence of behavioral experience was examined in adult primates (squirrel monkeys). Maps of movement representations were derived under anesthesia using intracortical microstimulation, whereby a microelectrode was inserted into the cortex to electrically stimulate corticospinal neurons at low current levels and evoke movements of the forelimb, principally of the hand. Movement representations were examined before and at several times after training on behavioral tasks that emphasized use of the fingers. Two behavioral tasks were utilized that dissociated the repetition of motor activity from the acquisition of motor skills. One task was easy to perform, and as such promoted repetitive motor activity without learning. The other task was more difficult, requiring the acquisition of motor skills for successful performance. Kinematic analysis indicated that monkeys used a consistent set of forelimb movements during pellet extractions. Functional mapping revealed that repetitive motor activity during the easier task did not produce plastic changes in movement representations. Instead, map plasticity, in the form of selective expansions of task-related movement representations, was only produced following skill acquisition on the difficult task. Additional studies revealed that, in general, map plasticity persisted without further training for up to three months, in parallel with the retention of task-related motor skills. Also, extensive additional training on the small well task produced further improvements in performance, and further changes in movement maps. In sum, these experiments support the following three conclusions regarding the role of M1 in motor learning. First, behaviorally-driven plasticity is learning-dependent, not activity-dependent. Second, plastic changes in M1 functional representations represent a neural correlate of acquired motor skills. Third, the persistence of map plasticity suggests that M1 is part of the neural substrate for the memory of motor skills. ^