Quantitative Analysis of Protein Expression to Study Lineage Specification in Mouse Preimplantation Embryos

This protocol presents a method to perform quantitative, single-cell in situ analyses of protein expression to study lineage specificationin mouse preimplantation embryos. The procedures necessary for embryo collection, immunofluorescence, imaging on a confocal microscope, and image segmentation and analysis are described. This method allows quantitation of the expression of multiple nuclear markers and the spatial (XYZ) coordinates of all cells in the embryo. It takes advantage of MINS, an image segmentation software tool specifically developed for the analysis of confocal images of preimplantation embryos and embryonic stem cell (ESC) colonies. MINS carries out unsupervised nuclear segmentation across the X, Y and Z dimensions, and produces information on cell position in three-dimensional space, as well as nuclear fluorescence levels for all channels with minimal user input. While this protocol has been optimized for the analysis of images of preimplantation stage mouse embryos, it can easily be adapted to the analysis of any other samples exhibiting a good signal-to-noise ratio and where high nuclear density poses a hurdle to image segmentation (e.g., expression analysis of embryonic stem cell (ESC) colonies, differentiating cells in culture, embryos of other species or stages, etc.).     

Morphogenetic changes during newt tail regeneration under changed gravity conditions

Gravity-dependent shape alterations in newt tail regenerates are described, which were previously noticed in experiments onboard satellites Foton M2, M3 and in corresponding laboratory controls. Laboratory conditions were developed that allow reproducing this phenomenon persistently in the adult newts Pleurodeles waltl (Michahelles, 1830). The newts kept in an aquarium (in partial weightlessness) after 1/3 tail amputation developed normal lanceolate regenerates, while those that stayed on a moist mat (exposed to greater gravity than in aquarium) developed curved tail regenerates. Dynamics of the shape alterations were described using computer morphometric analysis. The curve was shown to develop at stage III of regeneration and to be caused by bending of the developing axial structures: the ependymal tube and the cartilage rode. Cellular processes were described that accompany the tail shape changes, such as cell migration and formation of dense aggregates. Unequal proliferation throughout the wound epidermis and blastema was revealed using BrdU assay. Proliferation increased within dorsal and apical regions of the regenerates in the newts kept on the mat cell compared with the aquarian animals.     

Segmental organisation of the tail region in the embryo of Drosophila melanogaster

Summary The segmental organisation of the tail region in the embryo of Drosophila melanogaster, which is defined here as the epidermal region posterior to the boundary between abdominal segments A7 and A8, has been investigated by means of ultraviolet (UV) laser fate-mapping and phenotypic analysis of embryonic mutants that alter the segmental pattern of the larval cuticle. Wild-type embryos were irradiated in the presumptive tail region with a UV- laser microbeam of 20 m diameter at the blastoderm stage. The ensuing defects were scored in the cuticle pattern of the tail region of the first-instar larva, which is described in detail in this paper. The spatial distribution of defect frequencies was used to construct a blastoderm fate-map of the cuticle structures of the larval tail region. The segmental origin of the larval tail structures was inferred from the phenotypic analysis of segmentation and homoeotic mutants, which revealed pattern repetition throughout the embryonic tail region corresponding to four segment anlagen, A8 to A11, and a non-segmental telson. These data enabled the transformation of the blastoderm fate-map of cuticle structures into a map of tail segment anlagen. The tail anlage occupies about 10% of the egg length (EL), bounded by segment A7 anteriorly at 20% EL and by the proctodaeum posteriorly at 10% EL, as measured from the posterior pole. The anlagen of segments A8 and A9 appear to be narrow dorso-ventral strips of blastoderm cells similar to the anlagen of the trunk segments, whereas the anlagen of A10 and A11 are smaller and produce fewer pattern elements. The telson is represented in the cuticle by the tuft which derives from a very dorsal posterior position. The antero-posterior axis of the entire tail anlage appears curved upward posteriorly. Differences in the mode of development between tail and trunk segments are discussed, as are similarities of larval and imaginal tail development in Drosophila. Comparison with tail development in other insects suggests that, during evolution, the transition from semi-long-germ to long-germ development modified the organisation of the tail region without affecting its primary subdivision into metameric units.     

A method for quantifying rotational symmetry

Here, a new approach for quantifying rotational symmetry based on vector analysis was described and compared with information obtained from a geometric morphometric analysis and a technique based on distance alone. A new method was developed that generates a polygon from the length and angle data of a structure and then quantifies the minimum change necessary to convert that polygon into a regular polygon. This technique yielded an asymmetry score (s) that can range from 0 (perfect symmetry) to 1 (complete asymmetry). Using digital images of Geranium robertianum flowers, this new method was compared with a technique based on lengths alone and with established geometric morphometric methods used to quantify shape variation. Asymmetry scores (s) more clearly described variation in symmetry and were more consistent with a visual assessment of the images than either comparative technique. This procedure is the first to quantify the asymmetry of radial structures accurately, uses easily obtainable measures to calculate the asymmetry score and allows comparisons among individuals and species, even when the comparisons involve structures with different patterns of symmetry. This technique enables the rigorous analysis of polysymmetric structures and provides a foundation for a better understanding of symmetry in nature.     

Cell shape and cell division

The correlation between cell shape elongation and the orientation of the division axis described by early cell biologists is still used as a paradigm in developmental studies. However, analysis of early embryo development and tissue morphogenesis has highlighted the role of the spatial distribution of cortical cues able to guide spindle orientation. In vitro studies of cell division have revealed similar mechanisms. Recent data support the possibility that the orientation of cell division in mammalian cells is dominated by cell adhesion and the associated traction forces developed in interphase. Cell shape is a manifestation of these adhesive and tensional patterns. These patterns control the spatial distribution of cortical signals and thereby guide spindle orientation and daughter cell positioning. From these data, cell division appears to be a continuous transformation ensuring the maintenance of tissue mechanical integrity.     

Biometrics,biomathematics and the morphometric synthesis

At the core of contemporarymorphometrics—the quantitative study of biological shape variation—is a synthesis of two originally divergent methodological styles. One contributory tradition is the multivariate analysis of covariance matrices originally developed as biometrics and now dominant across a broad expanse of applied statistics. This approach, couched solely in the linear geometry of covariance structures, ignores biomathematical aspects of the original measurements. The other tributary emphasizes the direct visualization of changes in biological form. However, making objective the biological meaning of the features seen in those diagrams was always problematical; also, the representation of variation, as distinct from pairwise difference, proved infeasible. To combine these two variants of biomathematical modeling into a valid praxis for quantitative studies of biological shape was a goal earnestly sought though most of this century. That goal was finally achieved in the 1980s when techniques from mathematical statistics, multivariate biometrics, non-Euclidean geometry and computer graphics were combined in a coherent new system of tools for the complete regionalized quantitative analysis oflandmark points together with the biomedical images in which they are seen. In this morphometric synthesis, correspondence of landmarks (biologically labeled geometric points, like “bridge of the nose”) across specimens is taken as a biomathematical primitive. The shapes of configurations of landmarks are defined as equivalence classes with respect to the Euclidean similarity group and then represented as single points in David Kendall'sshape space, a Riemannian manifold with Procrustes distance as metric. All conventional multivariate strategies carry over to the study of shape variation and covariation when shapes are interpreted in the tangent space to the shape manifold at an average shape. For biomathematical interpretation of such analyses, one needs a basis for the tangent space compatible with the reality of local biotheoretical processes and explanations at many different geometric scales, and one needs graphics for visualizing average shape differences and other statistical contrasts there. Both of these needs are managed by thethin-plate spline, a deformation function that has an unusually helpful linear algebra. The spline also links the biometrics of landmarks to deformation analysis of the images from which the landmarks originally arose. This article reviews the history and principal tools of this synthesis in their biomathematical and biometrical context and demonstrates their usefulness in a study of focal neuroanatomical anomalies in schizophrenia.     

Morphomechanical programming of morphogenesis in cnidarian embryos

The factors governing the pattern formation process in the early morphogenesis of a marine colonial hydroid, Dynamena pumila, have been studied. Two different types of morphogenesis have been distinguished. Morphogenesis of the first type goes on via changes in cell shape and cell axis orientation, while morphogenesis of the second type is based upon the active coordinated cell movements associated with cell rearrangements. It was shown that morphogenesis of both types can be considered as cascades in which any event is a consequence of the previous one. The spatial structure of each developmental stage contains information about the direction and the initial conditions of further morphogenesis. So, an "epigenetic program" of morphogenesis gradually originates in the course of development and provides the stable reproduction of spatial structures. It is reasonable to consider the activity of epigenetic factors guiding Dynamena morphogenesis (geometry/topology of an embryo, heterogeneity of an embryo spatial structure, configuration of the field of mechanical stresses of the embryo surface) as "morphomechanical programming" of morphogenesis.     

Systematic spatial mapping of proteins at exocytic and endocytic structures

Vesicular secretion (exocytosis) involves the release and then compensatory recycling of vesicle components through endocytosis. This fundamental cellular process is controlled by the coordinated assembly and interactions of dozens of proteins at the plasma membrane. Understanding the molecular composition of individual exocytic and endocytic structures and their organization across the plasma membrane is critical to understanding the behavior and regulation of these two cellular processes. Here we develop a high-resolution and high-throughput fluorescence imaging–based approach for the unbiased mapping of 78 proteins at single exocytic vesicles and endocytic structures in neuroendocrine PC12 cells. This analysis uses two-color single-frame images to provide a systems-level map of the steady-state distributions of proteins at individual exocytic and endocytic structures in the cell. Along with this quantitative map, we find that both calcium-regulated exocytic vesicles (dense core vesicles) and endocytic structures (clathrin-coated structures) and the proteins associated with these structures exhibit a random spatial distribution in unstimulated neuroendocrine PC12 cells. This approach is broadly applicable for quantitatively mapping the molecular composition and spatial organization of discrete cellular processes with central molecular hubs.     

Pixel-based analysis of multiple images for the identification of changes: a novel approach applied to unravel proteome patterns [corrected] of 2-D electrophoresis gel images

A novel approach for revealing patterns of proteome variation among series of 2-DE gel images is presented. The approach utilises image alignment to ensure that each pixel represents the same information across all gels. Gel images are normalised, and background corrected, followed by unfolding of the images to 1-D pixel vectors and analysing pixel vectors by multivariate data modelling. Information resulting from the data analysis is refolded back to the image domain for visualisation and interpretation. The method is rapid and suitable for automatic routines applied after the gel alignment. The approach is compared with spot volume analysis to illustrate how this approach can solve persistent problems like mismatch of protein spots, erroneous missing values and failure to detect variation in overlapping proteins. The method may also detect variation in the border area of saturated proteins. The approach is given the name pixel-based analysis of multiple images for the identification of changes (PMC). The method can be used for multiple images in general. Effects of pretreatment of the images are discussed.     

A method for quantitative analysis of spatially variable physiological processes across leaf surfaces

Many physiological processes are spatially variable across leaf surfaces. While maps of photosynthesis, stomatal conductance, gene expression, water transport, and the production of reactive oxygen species (ROS) for individual leaves are readily obtained, analytical methods for quantifying spatial heterogeneity and combining information gathered from the same leaf but with different instruments are not widely used. We present a novel application of tools from the field of geographical imaging to the multivariate analysis of physiological images. Procedures for registration and resampling, cluster analysis, and classification provide a general framework for the analysis of spatially resolved physiological data. Two experiments were conducted to illustrate the utility of this approach. Quantitative analysis of images of chlorophyll fluorescence and the production of ROS following simultaneous exposure of soybean leaves to atmospheric O₃ and soybean mosaic virus revealed that areas of the leaf where the operating quantum efficiency of PSII was depressed also experienced an accumulation of ROS. This correlation suggests a causal relationship between oxidative stress and inhibition of photosynthesis. Overlaying maps of leaf surface temperature and chlorophyll fluorescence following a photoinhibition treatment indicated that areas with low operating quantum efficiency of PSII also experienced reduced stomatal conductance (high temperature). While each of these experiments explored the covariance of two processes by overlaying independent images gathered with different instruments, the same procedures can be used to analyze the covariance of information from multiple images. The application of tools from geographic image analysis to physiological processes occurring over small spatial scales will help reveal the mechanisms generating spatial variation across leaves.     

Developmental plasticity, morphological variation and evolvability: a multilevel analysis of morphometric integration in the shape of compound leaves

The structure of compound leaves provides flexibility for morphological change by variation in the shapes, sizes and arrangement of leaflets. Here, we conduct a multilevel analysis of shape variation in compound leaves to explore the developmental plasticity and evolutionary potential that are the basis of diversification in leaf shape. We use the methods of geometric morphometrics to study the shapes of individual leaflets and whole leaves in 20 taxa of Potentilla (sensu lato). A newly developed test based on the bootstrap approach suggests that uncertainty in the molecular phylogeny precludes firm conclusions whether there is a phylogenetic signal in the data on leaf shape. For variation among taxa, variation within taxa, as well as fluctuating asymmetry, there is evidence of strong morphological integration. The patterns of variation are similar across all three levels, suggesting that integration within taxa may act as a constraint on evolutionary change.     

Haltere morphology and campaniform sensilla arrangement across Diptera

One of the primary specializations of true flies (order Diptera) is the modification of the hind wings into club-shaped halteres. Halteres are complex mechanosensory structures that provide sensory feedback essential for stable flight control via an array of campaniform sensilla at the haltere base. The morphology of these sensilla has previously been described in a small number of dipteran species, but little is known about how they vary across fly taxa. Using a synoptic set of specimens representing 42 families from all of the major infraorders of Diptera, we used scanning electron microscopy to map the gross and fine structures of halteres, including sensillum shape and arrangement. We found that several features of haltere morphology correspond with dipteran phylogeny: Schizophora generally have smaller halteres with stereotyped and highly organized sensilla compared to nematoceran flies. We also found a previously undocumented high variation of haltere sensillum shape in nematoceran dipterans, as well as the absence of a dorsal sensillum field in multiple families. Overall, variation in haltere sensillar morphology across the dipteran phylogeny provides insight into the evolution of a highly specialized proprioceptive organ and a basis for future studies on haltere sensory function.     

Geometric morphometrics meets metacommunity ecology: environment and lineage distribution affects spatial variation in shape

Patterns of univariate trait variation across metacommunities are widely explored, as are searches for their underlying causes. Surprisingly, patterns of multivariate shape remain unknown, and the search for drivers of functional traits of communities often neglect the biogeographical distribution of phylogenetic clades. Our aim was to investigate multivariate shape distribution across metacommunities and to determine the main environmental drivers of shape beyond/taking into account the phylogenetic distribution of lineages. We obtained mean skull and mandible shape for 228 species of Neotropical sigmodontine rodents through geometric morphometrics (GM), and then calculated mean shapes for 1° × 1° cells across the Neotropics based on the incidence of sigmodontines. We investigated the effects of lineage distribution on mean trait variation by using phylogenetic fuzzy weighting to calculate principal coordinates of phylogenetic structure (PCPS). Effects of environmental variables on shape variation incorporating phylogenetic composition were realized through redundancy analysis. We found that the different distributions of major lineages throughout the Neotropics were responsible for much of the mean shape variation. The association of landscape features with tribal groupings (Oryzomyini with Amazonia and Phyllotini and Abrotrichini with the Andes) were standouts. Environmental variables and lineage distribution explain the same (i.e. shared) portion of shape variation, suggesting phylogenetic niche conservatism at the metacommunity level. Seasonality in temperature and land cover were the best environmental predictors of mean shape: larger tympanic bullae, incisive foramina, and check teeth are all associated with highly seasonal and less vegetated areas. Our new approach of using GM shape across metacommunities was demonstrably useful in understanding large‐scale biogeographical patterns of shape variation and identifying its underlying causes. The overlap between environmental variables and phylogenetic lineage distribution suggests that a process of niche conservatism is likely: the phenotype–environment correlation is mediated by the differential biogeographical distribution of the main clades.     

Predicting present and future intra-specific genetic structure through niche hindcasting across 24 millennia

Paleoclimatic reconstructions coupled with species distribution models and identification of extant spatial genetic structure have the potential to provide insights into the demographic events that shape the distribution of intra-specific genetic variation across time. Using the globeflower Trollius europaeus as a case-study, we combined (1) Amplified Fragment Length Polymorphisms, (2) suites of 1000-years stepwise hindcasted species distributions and (3) a model of diffusion through time over the last 24,000 years, to trace the spatial dynamics that most likely fits the species' current genetic structure. We show that the globeflower comprises four gene pools in Europe which, from the dry period preceding the Last Glacial Maximum, dispersed while tracking the conditions fitting its climatic niche. Among these four gene pools, two are predicted to experience drastic range retraction in the near future. Our interdisciplinary approach, applicable to virtually any taxon, is an advance in inferring how climate change impacts species' genetic structures.     

Spatial relationships between fibroblasts during the growth of rat-tail tendon

Sections of tendons from the base of the tail of rats were taken at eight time intervals from 18 days in utero until 244 days after birth and were examined in the electron microscope. For each time period, measurements were made of the relative area of fibroblasts, collagen and interstitial material, of the number of fibroblasts per unit area of tendon and of the average area of individual fibroblasts. The spatial arrangement of fibroblasts in the tendon sections was described quantitatively using the "nearest neighbor" method. Initially there was a rapid increase in the area of collagen accompanied by a decrease in the area occupied by fibroblasts but after 104 days of age these values changed very little. The numbers of fibroblasts per unit area decreased steadily from the embryo until 104 days whereas the average size of each cell increased to reach a maximum area at 40 days of age and then declined. At all time intervals cells were arranged in a regular, dispersed pattern across the tendon fascicles. Growth in width of the rat tail appears to involve the secretion of collagen and other intercellular material symmetrically around each fibroblast, so as to gradually separate the cells until a stage is reached at which cells are sufficiently far apart that there is little contact between adjacent cell processes. This may interfere with the integration of metabolic activity in the tissue. As a consequence, there is shrinkage of the cell bodies and a reduction in secretory activity so that, between 55 and 104 days of age, the tendon enters a period of terminal senescence.     

Nucleotide-dependent shape changes in the reverse direction motor, myosin VI

We have studied the shape of myosin VI, the actin minus-end directed motor, by negative stain and metal shadow electron microscopy. Single particle processing was used to make two-dimensional averages of the stain images, which greatly increases the clarity and allows detailed comparisons with crystal structures. A total of 169,964 particle images were obtained from two different constructs in six different states (four nucleotide states and with and without Ca²⁺). The shape of truncated apo myosin VI was very similar to the apo crystal structure, with the lever arm bent strongly backward and around the motor domain. In the full-length molecule, the C-terminal part of the tail has an additional bend taking it back across the motor domain, which may reflect a regulated state. Addition of ATP, ADP, or ATP-γS resulted in a large change, straightening the molecule from the bent shape and swinging the lever by ∼140°. Although these nucleotides would not be expected to produce the pre-powerstroke state, myosin VI in their presence was most similar to the truncated crystal structure with bound ADP-VO₄, which is thought to show the pre-powerstroke shape. The nucleotide data were therefore substantially different from expectation based on crystal structures. The full-length molecule was almost completely monomeric; only ∼1% were dimers, joined through the ends of the tail. Addition of calcium ions appeared to result in release of the second calmodulin light chain. In negatively stained molecules there was little indication of extended α-helical structure in the tail, but molecules viewed by metal shadowing had a tail ∼3× longer, 29 vs. 9 nm, part of which is likely to be a single α-helix.     

Detecting Nematode Features from Digital Images

Procedures for estimating and calibrating nematode features from digitial images are described and evaluated by illustration and mathematical formulae. Technical problems, such as capturing and cleaning raw images, standardizing the grey level range of images, and the detection of characteristics of the body habitus, presence or absence of stylet knobs, and tail and lip region shape are discussed. This study is the first of a series aimed at developing a set of automated methods to permit more rapid, objective characterizations of nematode features than is achievable by cumbersome conventional methods.