A quantitative genetic analysis of localized morphology in mandibles of inbred mice using finite element scaling analysis.

We analyzed patterns of mandibular genetic and phenotypic morphological integration and the relationship of genealogy to interstrain molecular and morphological differences in ten inbred strains of mice. Positions of mandibular landmarks in two-dimensional space were used to construct a finite element mesh for each individual, then all individuals from the ten strains were compared to the average mandible from a standard strain (SEA/GnJ). Measures of size and shape associated with finite element scaling analysis were then used in a quantitative genetic analysis of mandibular variation. Significant genetic variation for mandibular size and shape was uncovered. Patterns of both genetic and phenotypic correlation for measures of landmark-specific sizes were consistent with models of morphological integration based on the developmental origin of parts of the mandible and on the effects of muscle attachment on mandibular morphology. Shape differences local to particular landmarks did not show these forms of morphological integration. Although interstrain distances based on local shape magnitudes were significantly correlated with genealogical relationship, distances based on local size differences were not. Even higher than the correlation of genealogy with distances based on local shape magnitude was the genealogical-molecular distance correlation. Patterns of morphometric mandibular variation corresponded to expected effects of epigenetic developmental processes. Also, when detailed shape differences were considered, morphology served as a rough guide to genealogy, although molecular distances showed a stronger relationship.     

The Measurement of Local Variation in Shape

Geometric morphometrics comprises tools for measuring and analyzing shape as captured by an entire set of landmark configurations. Many interesting questions in evolutionary, genetic, and developmental research, however, are only meaningful at a local level, where a focus on ??parts?? or ??traits?? takes priority over properties of wholes. To study variational properties of such traits, current approaches partition configurations into subsets of landmarks which are then studied separately. This approach is unable to fully capture both variational and spatial characteristics of these subsets because interpretability of shape differences is context-dependent. Landmarks omitted from a partition usually contain information about that partition??s shape. We present an interpolation-based approach that can be used to model shape differences at a local, infinitesimal level as a function of information available globally. This approach belongs in a large family of methods that see shape differences as continuous ??fields?? spanning an entire structure, for which landmarks serve as reference parameters rather than as data. We show, via analyses of simulated and real data, how interpolation models provide a more accurate representation of regional shapes than partitioned data. A key difference of this interpolation approach from current morphometric practice is that one must assume an explicit interpolation model, which in turn implies a particular kind of behavior of the regions between landmarks. This choice presents novel methodological challenges, but also an opportunity to incorporate and test biomechanical models that have sought to explain tissue-level processes underlying the generation of morphological shape.     

An Integrated Approach for Landmark-Based Resistant Shape Analysis in 3D

The study of shape changes in morphology has seen a significant renovation in the last 20 years, particularly as a consequence of the development of geometric morphometric methods based on Cartesian coordinates of points. In order to extract information about shape differences when Cartesian coordinates are used, it is necessary to establish a common reference frame or system for all specimens to be compared. Therefore, a central issue in coordinate-based methods is which criterion should be used to align these configurations of points, since shape differences highly depend on those alignments. This is usually accomplished by aligning the configurations in a way that the sum of squared distances between coordinates of homologous points (landmarks) is minimized: the least-squares superimposition method. However, it is widely recognized that this method has some limitations when shape differences are not homogeneous across landmarks. Here we present an integrated approach for the resistant shape comparison of 3D landmark sets. It includes a new ordinary resistant Procrustes superimposition and its corresponding generalized resistant Procrustes version. In addition, they are combined with existing resistant multivariate statistical techniques for depicting the results. We demonstrate, by using both simulated and real datasets, that resistant Procrustes better detects and measures localized shape variation whenever present in up to half but one of the landmarks. The resistant Procrustes results are highly concordant with a priori biological information, and might dramatically improve the quality of inferences on patterns of shape variation.     

Use of a natural hybrid zone for genomewide association mapping of craniofacial traits in the house mouse

The identification of the genes involved in morphological variation in nature is still a major challenge. Here, we explore a new approach: we combine 178 samples from a natural hybrid zone between two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus), and high coverage of the genome (~ 145K SNPs) to identify loci underlying craniofacial shape variation. Due to the long history of recombination in the hybrid zone, high mapping resolution is anticipated. The combination of genomes from subspecies allows the mapping of both, variation within subspecies and inter‐subspecific differences, thereby increasing the overall amount of causal genetic variation that can be detected. Skull and mandible shape were measured using 3D landmarks and geometric morphometrics. Using principal component axes as phenotypes, and a linear mixed model accounting for genetic relatedness in the mapping populations, we identified nine genomic regions associated with skull shape and 10 with mandible shape. High mapping resolution (median size of significant regions = 148 kb) enabled identification of single or few candidate genes in most cases. Some of the genes act as regulators or modifiers of signalling pathways relevant for morphological development and bone formation, including several with known craniofacial phenotypes in mice and humans. The significant associations combined explain 13% and 7% of the skull and mandible shape variation, respectively. In addition, a positive correlation was found between chromosomal length and proportion of variation explained. Our results suggest a complex genetic architecture for shape traits and support a polygenic model.     

Morphological divergence of lake and stream Phoxinus of Northern Italy and the Danube basin based on geometric morphometric analysis

Minnows of the genus Phoxinus are promising candidates to investigate adaptive divergence, as they inhabit both still and running waters of a variety of altitudes and climatic zones in Europe. We used landmark‐based geometric morphometric methods to quantify the level of morphological variability in Phoxinus populations from streams and lakes of Northern Italy and the Danube basin. We analyzed body shape differences of populations in the dorsal, lateral, and ventral planes, using a large array of landmarks and semilandmarks. As the species identification of Phoxinus on morphological characters is ambiguous, we used two mitochondrial genes to determine the genetic background of the samples and to ensure we are comparing homogenous groups. We have found significant body shape differences between habitats: Minnow populations inhabiting streams had a deeper body and caudal peduncle and more laterally inserted pectoral fins than minnows inhabiting lakes. We have also found significant body shape differences between genetic groups: Italian minnows had deeper bodies, deeper and shorter caudal peduncles, and a shorter and wider gape than both groups from the Danube. Our results show that the morphology of Phoxinus is highly influenced by habitat and that body shape variation between habitats was within the same range as between genetic groups. These morphological differences are possibly linked to different modes of swimming and foraging in the respective habitats and are likely results of phenotypic plasticity. However, differences in shape and interlandmark distances between the groups suggest that some (though few) morphometric characters might be useful for separating Phoxinus species.     

Analysis of quantitative trait locus effects on the size and shape of mandibular molars in mice

While >50 genes have been found to influence the development of teeth in mice, we still know very little about the genetic basis for the adaptive characteristics of teeth, such as size and shape. We applied interval mapping procedures to Procrustes size and shape data obtained from 10 morphological landmarks on the mandibular molar row of the F(2) progeny from a cross between the LG/J and SM/J strains of mice. This revealed many more QTL for molar shape (18) than for molar centroid size (3), although levels of dominance effects were comparable among QTL for size and shape. Comparisons of patterns of Procrustes additive and dominance shape effects and ordination of QTL effects by principal components analysis suggested that the effects of the shape QTL were dispersed among the three molars and thus that none of these molars represents a genetically distinct developmental structure. The results of an analysis of co-occurrence of QTL for molar shape, mandible shape, and cranial dimensions in these mice suggested that many of the QTL for molar shape may be the same as those affecting these other sets of characters, although in some cases this could be due to effects of closely linked genes.     

Estimation of Distances and Variances in Bookstein's Landmark Model

Planar random figures can be described by series of landmarks. Bookstein's model assumes that the landmarks result from independent individual fluctuations around fixed fictive landmark centres. The aim of this paper is to estimate the distances between the centers and the variances of fluctuations around them. Furthermore, two new shape variable estimators are suggested and compared with an estimator of Bookstein.     

A hierarchical Bayesian model for next-generation population genomics

The demography of populations and natural selection shape genetic variation across the genome and understanding the genomic consequences of these evolutionary processes is a fundamental aim of population genetics. We have developed a hierarchical Bayesian model to quantify genome-wide population structure and identify candidate genetic regions affected by selection. This model improves on existing methods by accounting for stochastic sampling of sequences inherent in next-generation sequencing (with pooled or indexed individual samples) and by incorporating genetic distances among haplotypes in measures of genetic differentiation. Using simulations we demonstrate that this model has a low false-positive rate for classifying neutral genetic regions as selected genes (i.e., Φ(ST) outliers), but can detect recent selective sweeps, particularly when genetic regions in multiple populations are affected by selection. Nonetheless, selection affecting just a single population was difficult to detect and resulted in a high false-negative rate under certain conditions. We applied the Bayesian model to two large sets of human population genetic data. We found evidence of widespread positive and balancing selection among worldwide human populations, including many genetic regions previously thought to be under selection. Additionally, we identified novel candidate genes for selection, several of which have been linked to human diseases. This model will facilitate the population genetic analysis of a wide range of organisms on the basis of next-generation sequence data.     

Geometric estimates of heritability in biological shape

The recently developed geometric morphometrics methods represent an important contribution of statistics and geometry to the study of biological shapes. We propose simple protocols using shape distances that incorporate geometric techniques into linear quantitative genetic models that should provide insights into the contribution of genetics to shape variation in organisms. The geometric approaches use Procrustes distances in a curved shape space and distances in tangent spaces within and among families to estimate shape heritability. We illustrate the protocols with an example of wing shape variation in the honeybee, Apis mellifera. The heritability of overall shape variation was small, but some localized components depicting shape changes on distal wing regions showed medium to large heritabilities. The genetic variance-covariance matrix of the geometric shape variables was significantly correlated with the phenotypic shape variance-covariance matrix. A comparison of the results of geometric methods with the traditional multivariate analysis of interlandmark distances indicated that even with a larger dimensionality, the interlandmark distances were not as rich in shape information as the landmark coordinates. Quantitative genetics studies of shape should greatly benefit from the application of geometric methods.     

Variation in guenon skulls (I): species divergence, ecological and genetic differences

Guenons are the most diverse clade of African monkeys. They have varied ecologies, include arboreal and terrestrial species, and can be found in nearly every region of sub-Saharan Africa. Species boundaries are often uncertain, with a variable number of species and subspecies mostly recognised on the basis of their geographic distribution and pelage. If guenon soft tissue patterns show high variability, the same does not seem to hold for skull morphology. Guenon skulls are traditionally considered relatively undifferentiated and homogeneous. However, patterns of variation in skulls have never been examined using a large number of specimens sampled across the breadth of species diversity. Thus, in the present study, skulls of adult guenons and two outgroup species are analysed using three-dimensional geometric morphometrics. Three-dimensional coordinates of 86 anatomical landmarks were measured on 1,315 adult specimens belonging to all living guenon species except Cercopithecus dryas. Species are well-discriminated using shape but the best discrimination occurs when species have either a long evolutionary history (e.g., Allenopithecus nigroviridis) or represent extremes of size variation (Miopithecus sp. and Erythrocebus patas). Interspecific phenetic relationships reflect size differences. Four main clusters are found that mainly correspond to four size groups: the smallest species (Miopithecus sp.), the largest species (E. patas plus the study outgroups), a group of medium-small arboreal guenons, and a group of medium-large arboreal and terrestrial guenons. Correlations between interspecific shape distances and interspecific differences in size are higher than between shape distances and genetic distances. However, if only the component of interspecific shape variation which is not correlated to evolutionary allometry is used in the comparison with genetic distances, correlations are up to 1.4 times larger than those including allometric shape. The smallest correlations are those between shape and ecological distances, which is consistent with the lack of clusters clearly reflecting broad ecological specialisations (e.g., arboreality versus terrestriality). Thus, size, which is generally considered more evolutionarily labile than shape, seems to have played a major role in the evolution of the guenons. The incongruence between interspecific shape differences and phylogeny might be explained by a large proportion of shape changes having occurred along allometric trajectories that tend to be conserved within this clade.     

Genetic linkage analysis and homology relationships of genes located on human chromosome 11q.

We have used DNA polymorphisms detected by probes for 11q to order 16 genes and to determine the genetic distances between them. Our map includes the genes for CD20, tyrosinase, progesterone receptor, stromelysin, collagenase, N-CAM, dopamine-D2 receptor, apolipoproteins AI-CIII-AIV, CD3-epsilon, -delta, and -gamma, porphobilinogen deaminase, thy-1, and ets-1. These genes have previously been sequenced as well as placed on the 11q cytogenetic map, which now makes them anchor points between the cytogenetic, genetic, and physical maps of this region. The ordering and distances between these genes are of immediate use in testing hypotheses of candidate genes for human genetic diseases associated with chromosome 11q. A comparison between our genetic map and similar maps from other species defines regions of homologous synteny that may be useful in mapping human genetic disease genes localized to the 11q region. Analysis of such homology provides additional bases for speculation of the evolutionary histories of gene families in this region.     

Principles and methods of geometric morphometrics

The basic concepts, notions and methods of geometric morphometrics (GM) are considered. This approach implies multivariate analysis of landmark coordinates located following certain rules on the surface of a morphological object. The aim of GM is to reveal differences between morphological objects by their shapes as such, the "size factor" being excluded. The GM is based on the concept of Kendall's space (KS) defined as a hypersphere with points distributed on its surface. These points are the shapes defined as aligned landmark configurations. KS is a non-Euclidian space, its metrics called Procrustes is defined by landmark configuration of a reference shape relative to which other shapes are aligned and compared. The differences among shapes are measured as Procrustes distances between respective points. For the linear methods of multivariate statistics to be applied to comparison of shapes, the respective points are projected onto the tangent plane (tangent space), the tangent point being defined by the reference. There are two principal methods of shape comparisons in GM: the Procrustes superimposition (a version of the least squares analysis) and thin-plate spline analysis. In the first case, Procrustes residuals are the outcome shape variables which remain after isometric alignment of the shapes being compared. Their summation over all landmarks yields Procrustes distances among these shapes. The Procrustes distances can be used in multivariate analyses just as the Euclidian distances. In the second case, the shapes are fitted to the references by stretching/compressing and shearing until complete identity of their landmark configurations. Eigenvectors of resulting bending energy matrix are defined as new shape variables, principal warps which yield another shape space with the origin defined by the reference. Projections of the shapes being compared onto principal warps yield partial warps, and their covariance matrix decomposition into eigenvectors yields relative warps which are similar to principal components (in particular, they are mutually orthogonal). Both partial and relative warps can be used in many multivariate statistic analyses as quantitative shape variables. Results of thin-plate spline analysis can be represented graphically by transformation grid which displays type, amount and localization of the shape differences. Basis rules of sample composition and landmark positioning to be used in GM are considered. At present, rigid (with minimal degrees of freedom) 2D morphological objects are most suitable for GM applications. It is important to recognize three type of real landmarks, and additionally semi-landmarks and "virtual" landmarks. Some procedures of thin-plate spline analysis are considered exemplified by some study cases, as well as applications of some standard multivariate methods to GM results. They make it possible to evaluate correlation between different shapes, as well as between a shape and some non-shape variables (linear measurements etc); to evaluate the differences among organisms by shape of a morphological structure; to identify landmarks which most accounted for both correlation and differences between the shapes. An annotated list of most popular softwares for GM is provided.     

The developmental trajectory of brain-scalp distance from birth through childhood: implications for functional neuroimaging

Measurements of human brain function in children are of increasing interest in cognitive neuroscience. Many techniques for brain mapping used in children, including functional near-infrared spectroscopy (fNIRS), electroencephalography (EEG), magnetoencephalography (MEG) and transcranial magnetic stimulation (TMS), use probes placed on or near the scalp. The distance between the scalp and the brain is a key variable for these techniques because optical, electrical and magnetic signals are attenuated by distance. However, little is known about how scalp-brain distance differs between different cortical regions in children or how it changes with development. We investigated scalp-brain distance in 71 children, from newborn to age 12 years, using structural T1-weighted MRI scans of the whole head. Three-dimensional reconstructions were created from the scalp surface to allow for accurate calculation of brain-scalp distance. Nine brain landmarks in different cortical regions were manually selected in each subject based on the published fNIRS literature. Significant effects were found for age, cortical region and hemisphere. Brain-scalp distances were lowest in young children, and increased with age to up to double the newborn distance. There were also dramatic differences between brain regions, with up to 50% differences between landmarks. In frontal and temporal regions, scalp-brain distances were significantly greater in the right hemisphere than in the left hemisphere. The largest contributors to developmental changes in brain-scalp distance were increases in the corticospinal fluid (CSF) and inner table of the cranium. These results have important implications for functional imaging studies of children: age and brain-region related differences in fNIRS signals could be due to the confounding factor of brain-scalp distance and not true differences in brain activity.     

GEOMETRIC MORPHOMETRICS OF DEVELOPMENTAL INSTABILITY: ANALYZING PATTERNS OF FLUCTUATING ASYMMETRY WITH PROCRUSTES METHODS

Although fluctuating asymmetry has become popular as a measure of developmental instability, few studies have examined its developmental basis. We propose an approach to investigate the role of development for morphological asymmetry by means of morphometric methods. Our approach combines geometric morphometrics with the two-way ANOVA customary for conventional analyses of fluctuating asymmetry and can discover localized features of shape variation by examining the patterns of covariance among landmarks. This approach extends the notion of form used in studies of fluctuating asymmetry from collections of distances between morphological landmarks to an explicitly geometric concept of shape characterized by the configuration of landmarks. We demonstrate this approach with a study of asymmetry in the wings of tsetse flies (Glossina palpalis gambiensis). The analysis revealed significant fluctuating and directional asymmetry for shape as well as ample shape variation among individuals and between the offspring of young and old females. The morphological landmarks differed markedly in their degree of variability but multivariate patterns of landmark covariation identified by principal component analysis were generally similar between fluctuating asymmetry (within-individual variability) and variation among individuals. Therefore there is no evidence that special developmental processes control fluctuating asymmetry. We relate some of the morphometric patterns to processes known to be involved in the development of fly wings.     

Morphological integration of soft-tissue facial morphology in Down Syndrome and siblings

Down syndrome (DS), resulting from trisomy of chromosome 21, is the most common live-born human aneuploidy. The phenotypic expression of trisomy 21 produces variable, though characteristic, facial morphology. Although certain facial features have been documented quantitatively and qualitatively as characteristic of DS (e.g., epicanthic folds, macroglossia, and hypertelorism), all of these traits occur in other craniofacial conditions with an underlying genetic cause. We hypothesize that the typical DS face is integrated differently than the face of non-DS siblings, and that the pattern of morphological integration unique to individuals with DS will yield information about underlying developmental associations between facial regions. We statistically compared morphological integration patterns of immature DS faces (N = 53) with those of non-DS siblings (N = 54), aged 6-12 years using 31 distances estimated from 3D coordinate data representing 17 anthropometric landmarks recorded on 3D digital photographic images. Facial features are affected differentially in DS, as evidenced by statistically significant differences in integration both within and between facial regions. Our results suggest a differential affect of trisomy on facial prominences during craniofacial development.     

Multiple capacitors for natural genetic variation in Drosophila melanogaster

Cryptic genetic variation (CGV) or a standing genetic variation that is not ordinarily expressed as a phenotype is released when the robustness of organisms is impaired under environmental or genetic perturbations. Evolutionary capacitors modulate the amount of genetic variation exposed to natural selection and hidden cryptically; they have a fundamental effect on the evolvability of traits on evolutionary timescales. In this study, I have demonstrated the effects of multiple genomic regions of Drosophila melanogaster on CGV in wing shape. I examined the effects of 61 genomic deficiencies on quantitative and qualitative natural genetic variation in the wing shape of D. melanogaster. I have identified 10 genomic deficiencies that do not encompass a known candidate evolutionary capacitor, Hsp90, exposing natural CGV differently depending on the location of the deficiencies in the genome. Furthermore, five genomic deficiencies uncovered qualitative CGV in wing morphology. These findings suggest that CGV in wing shape of wild‐type D. melanogaster is regulated by multiple capacitors with divergent functions. Future analysis of genes encompassed by these genomic regions would help elucidate novel capacitor genes and better understand the general features of capacitors regarding natural genetic variation.     

Genetic and morphological assessment of a vulnerable large catfish,Silonia silondia (Hamilton, 1822), in natural populations from India

Silonia silondia is a commercially important fish distributed in Asian countries, which is under threat due to overexploitation. This study focuses on the morphological analysis and genetic variation of S. silondia individuals, through truss network and sequencing of two mitochondrial regions, respectively, from six wild populations of the Ganga and Mahanadi river systems in India. A total of 38 haplotypes was observed by analysing combined mitochondrial genes (cytochrome b + ATPase 6/8) in 247 individuals of S. silondia collected from six populations. Average haplotype and nucleotide diversities were 0.8508 and 0.00231, respectively. Genetic structure analysis showed the predominant cause of genetic variation to be within populations. The two clades were observed among the haplotypes and time of divergence from their most probable ancestor was estimated to be around 0.3949 mya. Analysis of combined mitochondrial genes in six populations of S. silondia resulted into three management units or genetic stocks. The truss network analysis was carried out by interconnecting 12 landmarks from digital images of specimens to identify phenotypic stocks. Sixty-five truss morphometric variables were analysed for geometric shape variation which revealed morphological divergence in River Son specimens. The present study presents molecular markers and genetic diversity data which can be critical input for conservation and management of differentiated populations and future monitoring of the genetic bottleneck. The morphological shape analysis clearly shows that variation in the insertion of adipose fin is an important parameter influencing the morphological discrimination.     

Genetic divergence within frog species is greater in topographically more complex regions

Most global hotspots of biodiversity and endemism are in montane regions. One explanation is that montane regions have intrinsically higher speciation rates than lowland regions because complex mountain topography and climate variation facilitate genetic isolation among populations. Here, we ask from an intraspecific perspective whether frog species whose haplotypes are connected by topographically/climatically complex regions display strong genetic isolation (greater scaled genetic distances), compared with species whose haplotypes are connected by less complex regions. We analysed published DNA sequences of several frog species from tropical Central and South America for the mitochondrial cob, cox1 and 16S rRNA genes. Pairwise genetic distances among haplotypes within each species were scaled to the geographic distances between each pair of haplotypes. Topographic complexity was positively correlated with scaled genetic distances, and isolation‐by‐resistance was supported only in species from more topographically complex regions. This suggests that heterogeneous topographies increase landscape resistance, which in turn favours the appearance of isolation‐by‐resistance. Moreover, we found that the potential barriers that restrict gene flow within species are more closely related to factors associated with temperature and topography than to precipitation.     

Linkage map of Escherichia coli K-12, edition 8.

The linkage map of Escherichia coli K-12 depicts the arrangement of genes on the circular chromosome of this organism. The basic units of the map are minutes, determined by the time-of-entry of markers from Hfr into F- strains in interrupted-conjugation experiments. The time-of-entry distances have been refined over the years by determination of the frequency of cotransduction of loci in transduction experiments utilizing bacteriophage P1, which transduces segments of DNA approximately 2 min in length. In recent years, the relative positions of many genes have been determined even more precisely by physical techniques, including the mapping of restriction fragments and the sequencing of many small regions of the chromosome. On the whole, the agreement between results obtained by genetic and physical methods has been remarkably good considering the different levels of accuracy to be expected of the methods used. There are now few regions of the map whose length is still in some doubt. In some regions, genetic experiments utilizing different mutant strains give different map distances. In other regions, the genetic markers available have not been close enough to give accurate cotransduction data. The chromosome is now known to contain several inserted elements apparently derived from lambdoid phages and other sources. The nature of the region in which the termination of replication of the chromosome occurs is now known to be much more complex than the picture given in the previous map. The present map is based upon the published literature through June of 1988. There are now 1,403 loci placed on the linkage group, which may represent between one-third and one-half of the genes in this organism.     

Crossing over in chicken oogenesis: cytological and chiasma-based genetic maps of the chicken lampbrush chromosome 1

Chiasmata in diplotene bivalents are located at the points of physical exchange (crossing-over) between homologous chromosomes. We have studied chiasma distribution within chicken lampbrush chromosome 1 to estimate the crossing-over frequency between chromosome landmarks. The position of the centromere and chromosome region 1q3.3-1q3.6 on lampbrush chromosome 1 were determined by comparative physical mapping of the TTAGGG repeats in the chicken mitotic and lampbrush chromosomes. The comparison of the chiasma (=crossing over)-based genetic distances on chicken chromosome 1 with the genetic linkage map obtained in genetic experiments showed that current genetic distances estimated by the high-resolution genetic mapping of the East Lansing, Compton, and Wageningen chicken reference populations are 1.2-1.9 times longer than those based on chiasma counts. Conceivable reasons for this discrepancy are discussed.