Statistical analysis of shape through triangulation of landmarks: A study of sexual dimorphism in hominids

Two objects with homologous landmarks are said to be of the same shape if the configuration of landmarks of one object can be exactly matched with that of the other by translation, rotation/reflection, and scaling. In an earlier paper, the authors proposed statistical analysis of shape by considering logarithmic differences of all possible Euclidean distances between landmarks. Tests of significance for differences in the shape of objects and methods of discrimination between populations were developed with such data. In the present paper, the corresponding statistical methodology is developed by triangulation of the landmarks and by considering the angles as natural measurements of shape. This method is applied to the study of sexual dimorphism in hominids.

Residual delay maps unveil global patterns of atmospheric nonlinearity and produce improved local forecasts

We use residual-delay maps of observational field data for barometric pressure to demonstrate the structure of latitudinal gradients in nonlinearity in the atmosphere. Nonlinearity is weak and largely lacking in tropical and subtropical sites and increases rapidly into the temperate regions where the time series also appear to be much noisier. The degree of nonlinearity closely follows the meridional variation of midlatitude storm track frequency. We extract the specific functional form of this nonlinearity, a V shape in the lagged residuals that appears to be a basic feature of midlatitude synoptic weather systems associated with frontal passages. We present evidence that this form arises from the relative time scales of high-pressure versus low-pressure events. Finally, we show that this nonlinear feature is weaker in a well regarded numerical forecast model (European Centre for Medium-Range Forecasts) because small-scale temporal and spatial variation is smoothed out in the grided inputs. This is significant, in that it allows us to demonstrate how application of statistical corrections based on the residual-delay map may provide marked increases in local forecast accuracy, especially for severe weather systems.

Multiple-complete-digest restriction fragment mapping: Generating sequence-ready maps for large-scale DNA sequencing

Multiple-complete-digest mapping is a DNA mapping technique based on complete-restriction-digest fingerprints of a set of clones that provides highly redundant coverage of the mapping target. The maps assembled from these fingerprints order both the clones and the restriction fragments. Maps are coordinated across three enzymes in the examples presented. Starting with yeast artificial chromosome contigs from the 7q31.3 and 7p14 regions of the human genome, we have produced cosmid-based maps spanning more than one million base pairs. Each yeast artificial chromosome is first subcloned into cosmids at a redundancy of ×15–30. Complete-digest fragments are electrophoresed on agarose gels, poststained, and imaged on a fluorescent scanner. Aberrant clones that are not representative of the underlying genome are rejected in the map construction process. Almost every restriction fragment is ordered, allowing selection of minimal tiling paths with clone-to-clone overlaps of only a few thousand base pairs. These maps demonstrate the practicality of applying the experimental and software-based steps in multiple-complete-digest mapping to a target of significant size and complexity. We present evidence that the maps are sufficiently accurate to validate both the clones selected for sequencing and the sequence assemblies obtained once these clones have been sequenced by a “shotgun” method.

Congruent Docking of Dimeric Kinesin and ncd into Three-dimensional Electron Cryomicroscopy Maps of Microtubule–Motor ADP Complexes

We present a new map showing dimeric kinesin bound to microtubules in the presence of ADP that was obtained by electron cryomicroscopy and image reconstruction. The directly bound monomer (first head) shows a different conformation from one in the more tightly bound empty state. This change in the first head is amplified as a movement of the second (tethered) head, which tilts upward. The atomic coordinates of kinesin·ADP dock into our map so that the tethered head associates with the bound head as in the kinesin dimer structure seen by x-ray crystallography. The new docking orientation avoids problems associated with previous predictions; it puts residues implicated by proteolysis-protection and mutagenesis studies near the microtubule but does not lead to steric interference between the coiled-coil tail and the microtubule surface. The observed conformational changes in the tightly bound states would probably bring some important residues closer to tubulin. As expected from the homology with kinesin, the atomic coordinates of nonclaret disjunctional protein (ncd)·ADP dock in the same orientation into the attached head in a map of microtubules decorated with dimeric ncd·ADP. Our results support the idea that the observed direct interaction between the two heads is important at some stages of the mechanism by which kinesin moves processively along microtubules.

Coexistence of linear zones and pinwheels within orientation maps in cat visual cortex

Revealing the layout of cortical maps is important both for understanding the processes involved in their development and for uncovering the mechanisms underlying neural computation. The typical organization of orientation maps in the cat visual cortex is radial; complete orientation cycles are mapped around orientation singularities. In contrast, long linear zones of orientation representation have been detected in the primary visual cortex of the tree shrew. In this study, we searched for the existence of long linear sequences and wide linear zones within orientation preference maps of the cat visual cortex. Optical imaging based on intrinsic signals was used. Long linear sequences and wide linear zones of preferred orientation were occasionally detected along the border between areas 17 and 18, as well as within area 18. Adjacent zones of distinct radial and linear organizations were observed across area 18 of a single hemisphere. However, radial and linear organizations were not necessarily segregated; long (7.5 mm) linear sequences of preferred orientation were found embedded within a typical pinwheel-like organization of orientation. We conclude that, although the radial organization is dominant, perfectly linear organization may develop and perform the processing related to orientation in the cat visual cortex.

A novel modifier gene for plasma von Willebrand factor level maps to distal mouse chromosome 11

Type 1 von Willebrand disease (VWD), characterized by reduced levels of plasma von Willebrand factor (VWF), is the most common inherited bleeding disorder in humans. Penetrance of VWD is incomplete, and expression of the bleeding phenotype is highly variable. In addition, plasma VWF levels vary widely among normal individuals. To identify genes that influence VWF level, we analyzed a genetic cross between RIIIS/J and CASA/Rk, two strains of mice that exhibit a 20-fold difference in plasma VWF level. DNA samples from F2 progeny demonstrating either extremely high or extremely low plasma VWF levels were pooled and genotyped for 41 markers spanning the autosomal genome. A novel locus accounting for 63% of the total variance in VWF level was mapped to distal mouse chromosome 11, which is distinct from the murine Vwf locus on chromosome 6. We designated this locus Mvwf for “modifier of VWF.” Additional genotyping of as many as 2407 meioses established a high resolution genetic map with gene order Cola1-Itg3a-Ngfr-Mvwf/Gip-Hoxb9-Hoxb1-Cbx·rs2-Cox5a-Gfap. The Mvwf candidate interval between Ngfr and Hoxb9 is ≈0.5 centimorgan (cM). These results demonstrate that a single dominant gene accounts for the low VWF phenotype of RIIIS/J mice in crosses with several other strains. The pattern of inheritance suggests a gain-of-function mutation in a unique component of VWF biosynthesis or processing. Characterization of the human homologue for Mvwf may have relevance for a subset of type 1 VWD cases and may define an important genetic factor modifying penetrance and expression of mutations at the VWF locus.

Formation of temporal-feature maps by axonal propagation of synaptic learning

Computational maps are of central importance to a neuronal representation of the outside world. In a map, neighboring neurons respond to similar sensory features. A well studied example is the computational map of interaural time differences (ITDs), which is essential to sound localization in a variety of species and allows resolution of ITDs of the order of 10 μs. Nevertheless, it is unclear how such an orderly representation of temporal features arises. We address this problem by modeling the ontogenetic development of an ITD map in the laminar nucleus of the barn owl. We show how the owl's ITD map can emerge from a combined action of homosynaptic spike-based Hebbian learning and its propagation along the presynaptic axon. In spike-based Hebbian learning, synaptic strengths are modified according to the timing of pre- and postsynaptic action potentials. In unspecific axonal learning, a synapse's modification gives rise to a factor that propagates along the presynaptic axon and affects the properties of synapses at neighboring neurons. Our results indicate that both Hebbian learning and its presynaptic propagation are necessary for map formation in the laminar nucleus, but the latter can be orders of magnitude weaker than the former. We argue that the algorithm is important for the formation of computational maps, when, in particular, time plays a key role.

Functional and structural mapping of human cerebral cortex: Solutions are in the surfaces

The human cerebral cortex is notorious for the depth and irregularity of its convolutions and for its variability from one individual to the next. These complexities of cortical geography have been a chronic impediment to studies of functional specialization in the cortex. In this report, we discuss ways to compensate for the convolutions by using a combination of strategies whose common denominator involves explicit reconstructions of the cortical surface. Surface-based visualization involves reconstructing cortical surfaces and displaying them, along with associated experimental data, in various complementary formats (including three-dimensional native configurations, two-dimensional slices, extensively smoothed surfaces, ellipsoidal representations, and cortical flat maps). Generating these representations for the cortex of the Visible Man leads to a surface-based atlas that has important advantages over conventional stereotaxic atlases as a substrate for displaying and analyzing large amounts of experimental data. We illustrate this by showing the relationship between functionally specialized regions and topographically organized areas in human visual cortex. Surface-based warping allows data to be mapped from individual hemispheres to a surface-based atlas while respecting surface topology, improving registration of identifiable landmarks, and minimizing unwanted distortions. Surface-based warping also can aid in comparisons between species, which we illustrate by warping a macaque flat map to match the shape of a human flat map. Collectively, these approaches will allow more refined analyses of commonalities as well as individual differences in the functional organization of primate cerebral cortex.

Imaging neuroscience: Principles or maps?

This article reviews some recent trends in imaging neuroscience. A distinction is made between making maps of functional responses in the brain and discerning the rules or principles that underlie their organization. After considering developments in the characterization of brain imaging data, several examples are presented that highlight the context-sensitive nature of neuronal responses that we measure. These contexts can be endogenous and physiological, reflecting the fact that each cortical area, or neuronal population, expresses its dynamics in the context of interactions with other areas. Conversely, these contexts can be experimental or psychological and can have a profound effect on the regional effects elicited. In this review we consider experimental designs and analytic strategies that go beyond cognitive subtraction and speculate on how functional imaging can be used to address both the details and principles underlying functional integration and specialization in the brain.

A bacterial artificial chromosome-based framework contig map of human chromosome 22q.

We have constructed a physical map of human chromosome 22q using bacterial artificial chromosome (BAC) clones. The map consists of 613 chromosome 22-specific BAC clones that have been localized and assembled into contigs using 452 landmarks, 346 of which were previously ordered and mapped to specific regions of the q arm of the chromosome by means of chromosome 22-specific yeast artificial chromosome clones. The BAC-based map provides immediate access to clones that are stable and convenient for direct genome analysis. The approach to rapidly developing marker-specific BAC contigs is relatively straightforward and can be extended to generate scaffold BAC contig maps of the rest of the chromosomes. These contigs will provide substrates for sequencing the entire human genome. We discuss how to efficiently close contig gaps using the end sequences of BAC clone inserts.

Gene for the catalytic subunit of mouse DNA-dependent protein kinase maps to the scid locus.

The gene encoding the catalytic subunit of DNA-dependent protein kinase (DNA-PKcs) has been proposed recently as a candidate gene for the mouse severe combined immune deficiency (scid) locus. We have used a partial cDNA clone for human DNA-PKcs to map the mouse homologue using a large interspecific backcross panel. We found that the mouse gene for DNA-PKcs does not recombine with scid, consistent with the hypothesis that scid is a mutation in the mouse gene for DNA-PKcs.

A maximum likelihood algorithm for genome mapping of cytogenetic loci from meiotic configuration data.

Frequencies of meiotic configurations in cytogenetic stocks are dependent on chiasma frequencies in segments defined by centromeres, breakpoints, and telomeres. The expectation maximization algorithm is proposed as a general method to perform maximum likelihood estimations of the chiasma frequencies in the intervals between such locations. The estimates can be translated via mapping functions into genetic maps of cytogenetic landmarks. One set of observational data was analyzed to exemplify application of these methods, results of which were largely concordant with other comparable data. The method was also tested by Monte Carlo simulation of frequencies of meiotic configurations from a monotelodisomic translocation heterozygote, assuming six different sample sizes. The estimate averages were always close to the values given initially to the parameters. The maximum likelihood estimation procedures can be extended readily to other kinds of cytogenetic stocks and allow the pooling of diverse cytogenetic data to collectively estimate lengths of segments, arms, and chromosomes.

Gene for the catalytic subunit of the human DNA-activated protein kinase maps to the site of the XRCC7 gene on chromosome 8.

The DNA-activated serine/threonine protein kinase (DNA-PK) is composed of a large (approximately 460 kDa) catalytic polypeptide (DNA-PKcs) and Ku, a heterodimeric DNA-binding component (p70/p80) that targets DNA-PKcs to DNA. A 41-kbp segment of the DNA-PKcs gene was isolated, and a 7902-bp segment was sequenced. The sequence contains a polymorphic Pvu II restriction enzyme site, and comparing the sequence with that of the cDNA revealed the positions of nine exons. The DNA-PKcs gene was mapped to band q11 of chromosome 8 by in situ hybridization. This location is coincident with that of XRCC7, the gene that complements the DNA double-strand break repair and V(D)J recombination defects (where V is variable, D is diversity, and J is joining) of hamster V3 and murine severe combined immunodeficient (scid) cells.

A region of consistent deletion in neuroblastoma maps within human chromosome 1p36.2-36.3.

Deletion of the short arm of human chromosome 1 is the most common cytogenetic abnormality observed in neuroblastoma. To characterize the region of consistent deletion, we performed loss of heterozygosity (LOH) studies on 122 neuroblastoma tumor samples with 30 distal chromosome 1p polymorphisms. LOH was detected in 32 of the 122 tumors (26%). A single region of LOH, marked distally by D1Z2 and proximally by D1S228, was detected in all tumors demonstrating loss. Also, cells from a patient with a constitutional deletion of 1p36, and from a neuroblastoma cell line with a small 1p36 deletion, were analyzed by fluorescence in situ hybridization. Cells from both sources had interstitial deletions of 1p36.2-36.3 which overlapped the consensus region of LOH defined by the tumors. Interstitial deletion in the constitutional case was confirmed by allelic loss studies using the panel of polymorphic markers. Four proposed candidate genes--DAN, ID3 (heir-1), CDC2L1 (p58), and TNFR2--were shown to lie outside of the consensus region of allelic loss, as defined by the above deletions. These results more precisely define the location of a neuroblastoma suppressor gene within 1p36.2-36.3, eliminating 33 centimorgans of proximal 1p36 from consideration. Furthermore, a consensus region of loss, which excludes the four leading candidate genes, was found in all tumors with 1p36 LOH.