Cochlea development shapes bat sensory system evolution

Sensory organs must develop alongside the skull within which they are largely encased, and this relationship can manifest as the skull constraining the organs, organs constraining the skull, or organs constraining one another in relative size. How this interplay between sensory organs and the developing skull plays out during the evolution of sensory diversity; however, remains unknown. Here, we examine the developmental sequence of the cochlea, the organ responsible for hearing and echolocation, in species with distinct diet and echolocation types within the ecologically diverse bat super-family Noctilionoidea. We found the size and shape of the cochlea largely correlates with skull size, with exceptions of Pteronotus parnellii, whose high duty cycle echolocation (nearly constant emission of sound pulses during their echolocation process allowing for detailed information gathering, also called constant frequency echolocation) corresponds to a larger cochlear and basal turn, and Monophyllus redmani, a small-bodied nectarivorous bat, for which interactions with other sensory organs restrict cochlea size. Our findings support the existence of developmental constraints, suggesting that both developmental and anatomical factors may act synergistically during the development of sensory systems in noctilionoid bats.

Variation in ontogenetic trajectories of limb dimensions in humans is attributable to both climatic effects and neutral evolution

Previous studies showed that there is variation in ontogenetic trajectories of human limb dimensions and proportions. However, little is known about the evolutionary significance of this variation. This study used a global sample of modern human immature long bone measurements and a multivariate linear mixed-effects model to study 1) whether the variation in ontogenetic trajectories of limb dimensions is consistent with ecogeographic predictions and 2) the effects of different evolutionary forces on the variation in ontogenetic trajectories. We found that genetic relatedness arising from neutral (nonselective) evolution, allometric variation associated with the change in size, and directional effects from climate all contributed to the variation in ontogenetic trajectories of all major long bone dimensions in modern humans. After accounting for the effects of neutral evolution and holding other effects considered in the current study constant, extreme temperatures have weak, positive associations with diaphyseal length and breadth measurements, while mean temperature shows negative associations with diaphyseal dimensions. The association with extreme temperatures fits the expectations of ecogeographic rules, while the association with mean temperature may explain the observed among-group variation in intralimb indices. The association with climate is present throughout ontogeny, suggesting an explanation of adaptation by natural selection as the most likely cause. On the other hand, genetic relatedness among groups, as structured by neutral evolutionary factors, is an important consideration when interpreting skeletal morphology, even for nonadult individuals.

Relating multivariate shapes to genescapes using phenotype-biological process associations for craniofacial shape

Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high-dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait – craniofacial shape – in a sample of n = 1145 mice from the Diversity Outbred (DO) experimental line. We use biological process Gene Ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (process MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application that provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex trait genetics.

Lewontin did not commit Lewontin's fallacy, his critics do: Why racial taxonomy is not useful for the scientific study of human variation

In 1972, R.C. Lewontin concluded that it follows from the fact that the large majority of human genetic variation (≈ 85%) is among individuals within local populations that racial taxonomy is unjustified. Three decades later, Edwards demonstrated that while the accuracy with which individuals may be assigned to groups is poor for a single locus, consideration of multi-locus data allows for highly accurate assignments. Edwards concluded that Lewontin's dismissal of racial taxonomy was unwarranted. Edwards misidentified the aim of Lewontin's critique, which was directed at the utility of racial classification and not at assigning individuals to groups using genetic data. Moreover, Edwards conflated distinct kinds of correlation when sketching out his argument. If we follow Edwards' argument to its natural terminus, it becomes clear that it is consideration of all of the correlation structure among local groups in human genetic data that renders racial taxonomy scientifically useless. Lewontin considers the correlation structure relevant to his analysis of racial taxonomy and does not make his eponymous misstep. Rather, critics of Lewontin who use racial taxonomies in their work are the primary offenders when it comes to committing Lewontin's fallacy.

Mammal Molar Size Ratios and the Inhibitory Cascade at the Intraspecific Scale

Mammalian molar crowns form a module in which measurements of size for individual teeth within a tooth row covary with one another. Molar crown size covariation is proposed to fit the inhibitory cascade model (ICM) or its variant the molar module component (MMC) model, but the inability of the former model to fit across biological scales is a concern in the few cases where it has been tested in Primates. The ICM has thus far failed to explain patterns of intraspecific variation, an intermediate biological scale, even though it explains patterns at both smaller organ-level and larger between-species biological scales. Studies of this topic in a much broader range of taxa are needed, but the properties of a sample appropriate for testing the ICM at the intraspecific level are unclear. Here, we assess intraspecific variation in relative molar sizes of the cotton mouse, Peromyscus gossypinus, to further test the ICM and to develop recommendations for appropriate sampling protocols in future intraspecific studies of molar size variation across Mammalia. To develop these recommendations, we model the sensitivity of estimates of molar ratios to sample size and simulate the use of composite molar rows when complete ones are unavailable. Similar to past studies on primates, our results show that intraspecific variance structure of molar ratios within the rodent P. gossypinus does not meet predictions of the ICM or MMC. When we extend these analyses to include the MMC, one model does not fit observed patterns of variation better than the other. Standing variation in molar size ratios is relatively constant across mammalian samples containing all three molars. In future studies, analyzing average ratio values will require relatively small minimum sample sizes of two or more complete molar rows. Even composite-based estimates from four or more specimens per tooth position can accurately estimate mean molar ratios. Analyzing variance structure will require relatively large sample sizes of at least 40-50 complete specimens, and composite molar rows cannot accurately reconstruct variance structure of ratios in a sample. Based on these results, we propose guidelines for intraspecific studies of molar size covariation. In particular, we note that the suitability of composite specimens for averaging mean molar ratios is promising for the inclusion of isolated molars and incomplete molar rows from the fossil record in future studies of the evolution of molar modules, as long as variance structure is not a key component of such studies.

Facial shape and allometry quantitative trait locus intervals in the Diversity Outbred mouse are enriched for known skeletal and facial development genes

The biology of how faces are built and come to differ from one another is complex. Discovering normal variants that contribute to differences in facial morphology is one key to untangling this complexity, with important implications for medicine and evolutionary biology. This study maps quantitative trait loci (QTL) for skeletal facial shape using Diversity Outbred (DO) mice. The DO is a randomly outcrossed population with high heterozygosity that captures the allelic diversity of eight inbred mouse lines from three subspecies. The study uses a sample of 1147 DO animals (the largest sample yet employed for a shape QTL study in mouse), each characterized by 22 three-dimensional landmarks, 56,885 autosomal and X-chromosome markers, and sex and age classifiers. We identified 37 facial shape QTL across 20 shape principal components (PCs) using a mixed effects regression that accounts for kinship among observations. The QTL include some previously identified intervals as well as new regions that expand the list of potential targets for future experimental study. Three QTL characterized shape associations with size (allometry). Median support interval size was 3.5 Mb. Narrowing additional analysis to QTL for the five largest magnitude shape PCs, we found significant overrepresentation of genes with known roles in growth, skeletal and facial development, and sensory organ development. For most intervals, one or more of these genes lies within 0.25 Mb of the QTL's peak. QTL effect sizes were small, with none explaining more than 0.5% of facial shape variation. Thus, our results are consistent with a model of facial diversity that is influenced by key genes in skeletal and facial development and, simultaneously, is highly polygenic.

Exerting an influence on evolution

Experiments on mice have shown that developmental processes are influencing the generation of phenotypic variation in a way that shapes evolution.

Variation in mouse pelvic morphology maps to locations enriched in Sox9 Class II and Pitx1 regulatory features

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high-resolution quantitative trait locus (QTL) mapping in the F₃₄ generation of an advanced intercross experiment (LG,SM-G₃₄ ) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome- or chromosome-wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation-seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Complexity, Genetic Causation, and Hereditarianism

Hereditarians have claimed that recent advances in psychological and psychiatric genetics support their contention that socially important aspects of behavior and cognition in individuals and groups are largely insensitive to environmental context. This has been countered by anti-hereditarians who (correctly) claim that the conclusion of genetic ineluctability is false. Anti-hereditarians, however, sometimes use problematic arguments based on complexity and the ignorance that comes with complexity and a demand for mechanistic, as opposed to variational, explanations for the ways in which genes affect phenotype. I argue here, as a committed anti-hereditarian, that the complexity gambit and the demand for mechanisms open anti-hereditarian arguments to counterattack from hereditarians. Refocusing the argument onto issues about when heritability, genotypic scores, and genome-wide association studies may be appropriately applied and reemphasizing the point that context matters are stronger measures to counter hereditarian claims.

Integration and the Developmental Genetics of Allometry

Allometry refers to the ways in which organismal shape is associated with size. It is a special case of integration, or the tendency for traits to covary, in that variation in size is ubiquitous and evolutionarily important. Allometric variation is so commonly observed that it is routinely removed from morphometric analyses or invoked as an explanation for evolutionary change. In this case, familiarity is mistaken for understanding because rarely do we know the mechanisms by which shape correlates with size or understand their significance. As with other forms of integration, allometric variation is generated by variation in developmental processes that affect multiple traits, resulting in patterns of covariation. Given this perspective, we can dissect the genetic and developmental determinants of allometric variation. Our work on the developmental and genetic basis for allometric variation in craniofacial shape in mice and humans has revealed that allometric variation is highly polygenic. Different measures of size are associated with distinct but overlapping patterns of allometric variation. These patterns converge in part on a common genetic basis. Finally, environmental modulation of size often generates variation along allometric trajectories, but the timing of genetic and environmental perturbations can produce deviations from allometric patterns when traits are differentially sensitive over developmental time. These results question the validity of viewing allometry as a singular phenomenon distinct from morphological integration more generally.

Nonlinear gene expression-phenotype relationships contribute to variation and clefting in the A/WySn mouse

BACKGROUND

Cleft lip and palate is one of the most common human birth defects, but the underlying etiology is poorly understood. The A/WySn mouse is a spontaneously occurring model of multigenic clefting in which 20% to 30% of individuals develop an orofacial cleft. Recent work has shown altered methylation at a specific retrotransposon insertion downstream of the Wnt9b locus in clefting animals, which results in decreased Wnt9b expression.

RESULTS

Using a newly developed protocol that allows us to measure morphology, gene expression, and DNA methylation in the same embryo, we relate gene expression in an individual embryo directly to its three-dimensional morphology for the first time. We find that methylation at the retrotransposon relates to Wnt9b expression and morphology. IAP methylation relates to shape of the nasal process in a manner consistent with clefting. Embryos with low IAP methylation exhibit increased among-individual variance in facial shape.

CONCLUSIONS

Methylation and gene expression relate nonlinearly to nasal process morphology. Individuals at one end of a continuum of phenotypic states display a clinical phenotype and increased phenotypic variation. Variable penetrance and expressivity in this model is likely determined both by among-individual variation in methylation and changes in phenotypic robustness along the underlying liability distribution for orofacial clefting.

Hybridization in human evolution: Insights from other organisms

During the late Pleistocene, isolated lineages of hominins exchanged genes thus influencing genomic variation in humans in both the past and present. However, the dynamics of this genetic exchange and associated phenotypic consequences through time remain poorly understood. Gene exchange across divergent lineages can result in myriad outcomes arising from these dynamics and the environmental conditions under which it occurs. Here we draw from our collective research across various organisms, illustrating some of the ways in which gene exchange can structure genomic/phenotypic diversity within/among species. We present a range of examples relevant to questions about the evolution of hominins. These examples are not meant to be exhaustive, but rather illustrative of the diverse evolutionary causes/consequences of hybridization, highlighting potential drivers of human evolution in the context of hybridization including: influences on adaptive evolution, climate change, developmental systems, sex-differences in behavior, Haldane's rule and the large X-effect, and transgressive phenotypic variation.

The developmental-genetics of canalization

Canalization, or robustness to genetic or environmental perturbations, is fundamental to complex organisms. While there is strong evidence for canalization as an evolved property that varies among genotypes, the developmental and genetic mechanisms that produce this phenomenon are very poorly understood. For evolutionary biology, understanding how canalization arises is important because, by modulating the phenotypic variation that arises in response to genetic differences, canalization is a determinant of evolvability. For genetics of disease in humans and for economically important traits in agriculture, this subject is important because canalization is a potentially significant cause of missing heritability that confounds genomic prediction of phenotypes. We review the major lines of thought on the developmental-genetic basis for canalization. These fall into two groups. One proposes specific evolved molecular mechanisms while the other deals with robustness or canalization as a more general feature of development. These explanations for canalization are not mutually exclusive and they overlap in several ways. General explanations for canalization are more likely to involve emergent features of development than specific molecular mechanisms. Disentangling these explanations is also complicated by differences in perspectives between genetics and developmental biology. Understanding canalization at a mechanistic level will require conceptual and methodological approaches that integrate quantitative genetics and developmental biology.

Developmental constraint through negative pleiotropy in the zygomatic arch

BACKGROUND

Previous analysis suggested that the relative contribution of individual bones to regional skull lengths differ between inbred mouse strains. If the negative correlation of adjacent bone lengths is associated with genetic variation in a heterogeneous population, it would be an example of negative pleiotropy, which occurs when a genetic factor leads to opposite effects in two phenotypes. Confirming negative pleiotropy and determining its basis may reveal important information about the maintenance of overall skull integration and developmental constraint on skull morphology.

RESULTS

We identified negative correlations between the lengths of the frontal and parietal bones in the midline cranial vault as well as the zygomatic bone and zygomatic process of the maxilla, which contribute to the zygomatic arch. Through gene association mapping of a large heterogeneous population of Diversity Outbred (DO) mice, we identified a quantitative trait locus on chromosome 17 driving the antagonistic contribution of these two zygomatic arch bones to total zygomatic arch length. Candidate genes in this region were identified and real-time PCR of the maxillary processes of DO founder strain embryos indicated differences in the RNA expression levels for two of the candidate genes, Camkmt and Six2.

CONCLUSIONS

A genomic region underlying negative pleiotropy of two zygomatic arch bones was identified, which provides a mechanism for antagonism in component bone lengths while constraining overall zygomatic arch length. This type of mechanism may have led to variation in the contribution of individual bones to the zygomatic arch noted across mammals. Given that similar genetic and developmental mechanisms may underlie negative correlations in other parts of the skull, these results provide an important step toward understanding the developmental basis of evolutionary variation and constraint in skull morphology.

Developmental nonlinearity drives phenotypic robustness

Robustness to perturbation is a fundamental feature of complex organisms. Mutations are the raw material for evolution, yet robustness to their effects is required for species survival. The mechanisms that produce robustness are poorly understood. Nonlinearities are a ubiquitous feature of development that may link variation in development to phenotypic robustness. Here, we manipulate the gene dosage of a signaling molecule, Fgf8, a critical regulator of vertebrate development. We demonstrate that variation in Fgf8 expression has a nonlinear relationship to phenotypic variation, predicting levels of robustness among genotypes. Differences in robustness are not due to gene expression variance or dysregulation, but emerge from the nonlinearity of the genotype–phenotype curve. In this instance, embedded features of development explain robustness differences. How such features vary in natural populations and relate to genetic variation are key questions for unraveling the origin and evolvability of this feature of organismal development.

Developmental processes often involve nonlinearities, but the consequences for translating genotype to phenotype are not well characterized. Here, Green et al. vary Fgf8 signaling across allelic series of mice and show that phenotypic robustness in craniofacial shape is explained by a nonlinear effect of Fgf8 expression.

Facial shape manifestations of growth faltering in Tanzanian children

Variation in the shape of the human face and in stature is determined by complex interactions between genetic and environmental influences. One such environmental influence is malnourishment, which can result in growth faltering, usually diagnosed by means of comparing an individual's stature with a set of age-appropriate standards. These standards for stature, however, are typically ascertained in groups where people are at low risk for growth faltering. Moreover, genetic differences among populations with respect to stature are well established, further complicating the generalizability of stature-based diagnostic tools. In a large sample of children aged 5-19 years, we obtained high-resolution genomic data, anthropometric measures and 3D facial images from individuals within and around the city of Mwanza, Tanzania. With genome-wide complex trait analysis, we partitioned genetic and environmental variance for growth outcomes and facial shape. We found that children with growth faltering have faces that look like those of older and taller children, in a direction opposite to the expected allometric trajectory, and in ways predicted by the environmental portion of covariance at the community and individual levels. The environmental variance for facial shape varied subtly but significantly among communities, whereas genetic differences were minimal. These results reveal that facial shape preserves information about exposure to undernourishment, with important implications for refining assessments of nutritional status in children and the developmental-genetics of craniofacial variation alike.

Body size and allometric variation in facial shape in children

OBJECTIVES

Morphological integration, or the tendency for covariation, is commonly seen in complex traits such as the human face. The effects of growth on shape, or allometry, represent a ubiquitous but poorly understood axis of integration. We address the question of to what extent age and measures of size converge on a single pattern of allometry for human facial shape.

METHODS

Our study is based on two large cross-sectional cohorts of children, one from Tanzania and the other from the United States (N = 7,173). We employ 3D facial imaging and geometric morphometrics to relate facial shape to age and anthropometric measures.

RESULTS

The two populations differ significantly in facial shape, but the magnitude of this difference is small relative to the variation within each group. Allometric variation for facial shape is similar in both populations, representing a small but significant proportion of total variation in facial shape. Different measures of size are associated with overlapping but statistically distinct aspects of shape variation. Only half of the size-related variation in facial shape can be explained by the first principal component of four size measures and age while the remainder associates distinctly with individual measures.

CONCLUSIONS

Allometric variation in the human face is complex and should not be regarded as a singular effect. This finding has important implications for how size is treated in studies of human facial shape and for the developmental basis for allometric variation more generally.

Ecogeography, genetics, and the evolution of human body form

Genetic resemblances among groups are non-randomly distributed in humans. This population structure may influence the correlations between traits and environmental drivers of natural selection thus complicating the interpretation of the fossil record when modern human variation is used as a referential model. In this paper, we examine the effects of population structure and natural selection on postcranial traits that reflect body size and shape with application to the more general issue of how climate - using latitude as a proxy - has influenced hominin morphological variation. We compare models that include terms reflecting population structure, ascertained from globally distributed microsatellite data, and latitude on postcranial phenotypes derived from skeletal dimensions taken from a large global sample of modern humans. We find that models with a population structure term fit better than a model of natural selection along a latitudinal cline in all cases. A model including both latitude and population structure terms is a good fit to distal limb element lengths and bi-iliac breadth, indicating that multiple evolutionary forces shaped these morphologies. In contrast, a model that included only a population structure term best explained femoral head diameter and the crural index. The results demonstrate that population structure is an important part of human postcranial variation, and that clinally distributed natural selection is not sufficient to explain among-group differentiation. The distribution of human body form is strongly influenced by the contingencies of modern human origins, which calls for new ways to approach problems in the evolution of human variation, past and present.

Genetic and environmental contributions to variation in baboon cranial morphology

The development, function, and integration of morphological characteristics are all hypothesized to influence the utility of traits for phylogenetic reconstruction by affecting the way in which morphological characteristics evolve. We use a baboon model to test the hypotheses about phenotypic and quantitative genetic variation of traits in the cranium that bear on a phenotype's propensity to evolve. We test the hypotheses that: 1) individual traits in different functionally and developmentally defined regions of the cranium are differentially environmentally, genetically, and phenotypically variable; 2) genetic covariance with other traits constrains traits in one region of the cranium more than those in others; 3) and regions of the cranium subject to different levels of mechanical strain differ in the magnitude of variation in individual traits. We find that the levels of environmental and genetic variation in individual traits are randomly distributed across regions of the cranium rather than being structured by developmental origin or degree of exposure to strain. Individual traits in the cranial vault tend to be more constrained by covariance with other traits than those in other regions. Traits in regions subject to high degrees of strain during mastication are not any more variable at any level than other traits. If these results are generalizable to other populations, they indicate that there is no reason to suppose that individual traits from any one part of the cranium are intrinsically less useful for reconstructing patterns of evolution than those from any other part.

Heritability of alveolar bone loss from periodontal disease in a baboon population: a pilot study

BACKGROUND

Reports from studies of twins, disease aggregation in families, animal models for periodontal disease, and various genetic-analysis studies have determined that genetics plays a role in the susceptibility to periodontal disease. The purpose of this pilot study is to evaluate the effect of genetics on periodontal disease by evaluating the heritability of alveolar bone loss in a captive baboon population.

METHODS

A collection of baboon skulls from a pedigreed colony (for which scientists and veterinarians maintain complete genealogic and veterinary records) was obtained from the Southwest National Primate Research Center, San Antonio, Texas and used in this pilot study. Measurements of alveolar bone loss were performed on 390 dry baboon skulls. A periodontal probe was used to measure alveolar bone loss. Maximum likelihood methods (designed to handle complex genealogies) were used to determine the heritability of alveolar bone loss. This software used known pedigrees in the captive baboon sample and tested the relationship between pairwise kinship and alveolar bone loss data to determine the heritability of alveolar bone loss from periodontal disease.

RESULTS

Genetic data were available for 347 of the 390 specimens. Using age and sex as covariates, genetic analysis indicated a heritability of 35% (standard error = 20%; P = 0.01). Although gender was not a significant factor in periodontal disease (P = 0.96), age was highly significantly associated with periodontal disease (P <0.0001).

CONCLUSIONS

In this pilot study, analysis of alveolar bone loss measurements from captive baboons indicates that bone loss increases with age and that a portion of periodontal disease risk may be caused by genetic variance. These findings provide evidence that periodontal disease is heritable in captive baboons and indicate that a larger, more-detailed study is warranted.

Fine-mapping of obesity-related quantitative trait loci in an F9/10 advanced intercross line

Obesity develops in response to a combination of environmental effects and multiple genes of small effect. Although there has been significant progress in characterizing genes in many pathways contributing to metabolic disease, knowledge about the relationships of these genes to each other and their joint effects upon obesity lags behind. The LG,SM advanced intercross line (AIL) model of obesity has been used to characterize over 70 loci involved in fatpad weight, body weight, and organ weights. Each of these quantitative trait loci (QTLs) encompasses large regions of the genome and require fine-mapping to isolate causative sequence changes and possible mechanisms of action as indicated by the genetic architecture. In this study we fine-map QTLs first identified in the F(2) and F(2/3) populations in the combined F(9/10) advanced intercross generations. We observed significantly narrowed QTL confidence regions, identified many single QTL that resolve into multiple QTL peaks, and identified new QTLs that may have been previously masked due to opposite gene effects at closely linked loci. We also present further characterization of the pleiotropic and epistatic interactions underlying these obesity-related traits.

Replication of long-bone length QTL in the F9-F10 LG,SM advanced intercross

Quantitative trait locus (QTL) mapping techniques are frequently used to identify genomic regions associated with variation in phenotypes of interest. However, the F(2) intercross and congenic strain populations usually employed have limited genetic resolution resulting in relatively large confidence intervals that greatly inhibit functional confirmation of statistical results. Here we use the increased resolution of the combined F(9) and F(10) generations (n = 1455) of the LG,SM advanced intercross to fine-map previously identified QTL associated with the lengths of the humerus, ulna, femur, and tibia. We detected 81 QTL affecting long-bone lengths. Of these, 49 were previously identified in the combined F(2)-F(3) population of this intercross, while 32 represent novel contributors to trait variance. Pleiotropy analysis suggests that most QTL affect three to four long bones or serially homologous limb segments. We also identified 72 epistatic interactions involving 38 QTL and 88 novel regions. This analysis shows that using later generations of an advanced intercross greatly facilitates fine-mapping of confidence intervals, resolving three F(2)-F(3) QTL into multiple linked loci and narrowing confidence intervals of other loci, as well as allowing identification of additional QTL. Further characterization of the biological bases of these QTL will help provide a better understanding of the genetics of small variations in long-bone length.

Comparison of Mandibular Phenotypic and Genetic Integration between Baboon and Mouse

In this study we compare patterns of mandibular integration between mice and baboons using both phenotypic and quantitative genetic data. Specifically, we test how well each species fits with the mosaic model of mandibular integration suggested by Atchley and Hall (Biol Rev Camb Philos Soc 66:101-157, 1991) based on developmental modules. We hypothesize that patterns of integration will be similar for mice and baboons and that both species will show strong integration within developmental modules and weaker integration between modules. Corresponding landmark data were collected from the hemi-mandibles of an advanced intercross mouse sample (N = 1239) and mandibles from a baboon sample of known pedigree from the Southwest Foundation for Biomedical Research (N = 430). We used four methods of analysis to quantify and compare the degree of mandibular integration between species including two methods based on a priori assumptions, and two a posteriori analyses. We found that patterns of integration are broadly similar for baboon and mouse mandibles, with both species displaying a modular pattern of integration. While there is a general trend of similarity in integration patterns between species, there were some marked differences. Mice are strongly correlated among distances within the coronoid process and the incisive alveolar region, whereas baboons are strongly integrated within the condylar process. We discuss the potential evolutionary implications of the similar patterns of integration between these species with an emphasis on the role of modularity.

Genetic variation in baboon craniofacial sexual dimorphism

Sexual dimorphism is a widespread phenomenon and contributes greatly to intraspecies variation. Despite a long history of active research, the genetic basis of dimorphism for complex traits remains unknown. Understanding the sex-specific differences in genetic architecture for cranial traits in a highly dimorphic species could identify possible mechanisms through which selection acts to produce dimorphism. Using distances calculated from three-dimensional landmark data from CT scans of 402 baboon skulls from a known genealogy, we estimated genetic variance parameters in both sexes to determine the presence of gene-by-sex (G x S) interactions and X-linked heritability. We hypothesize that traits exhibiting the greatest degree of sexual dimorphism (facial traits in baboons) will demonstrate either stronger G x S interactions or X-linked effects. We found G x S interactions and X-linked effects for a few measures that span the areas connecting the face to the neurocranium but for no traits restricted to the face. This finding suggests that facial traits will have a limited response to selection for further evolution of dimorphism in this population. We discuss the implications of our results with respect to the origins of cranial sexual dimorphism in this baboon sample, and how the genetic architecture of these traits affects their potential for future evolution.

Genetic architecture of adiposity and organ weight using combined generation QTL analysis

We present here a detailed study of the genetic contributions to adult body size and adiposity in the LG,SM advanced intercross line (AIL), an obesity model. This study represents a first step in fine-mapping obesity quantitative trait loci (QTLs) in an AIL. QTLs for adiposity in this model were previously isolated to chromosomes 1, 6, 7, 8, 9, 12, 13, and 18. This study focuses on heritable contributions and the genetic architecture of fatpad and organ weights. We analyzed both the F(2) and F(3) generations of the LG,SM AIL population single-nucleotide polymorphism (SNP) genotyped with a marker density of approximately 4 cM. We replicate 88% of the previously identified obesity QTLs and identify 13 new obesity QTLs. Nearly half of the single-trait QTLs were sex-specific. Several broad QTL regions were resolved into multiple, narrower peaks. The 113 single-trait QTLs for organs and body weight clustered into 27 pleiotropic loci. A large number of epistatic interactions are described which begin to elucidate potential interacting molecular networks. We present a relatively rapid means to obtain fine-mapping details from AILs using dense marker maps and consecutive generations. Analysis of the complex genetic architecture underlying fatpad and organ weights in this model may eventually help to elucidate not only heritable contributions to obesity but also common gene sets for obesity and its comorbidities.

Genetic variation in pleiotropy: differential epistasis as a source of variation in the allometric relationship between long bone lengths and body weight

Pleiotropy is an aspect of genetic architecture underlying the phenotypic covariance structure. The presence of genetic variation in pleiotropy is necessary for natural selection to shape patterns of covariation between traits. We examined the contribution of differential epistasis to variation in the intertrait relationship and the nature of this variation. Genetic variation in pleiotropy was revealed by mapping quantitative trait loci (QTLs) affecting the allometry of mouse limb and tail length relative to body weight in the mouse-inbred strain LG/J by SM/J intercross. These relationship QTLs (rQTLs) modify relationships between the traits affected by a common pleiotropic locus. We detected 11 rQTLs, mostly affecting allometry of multiple bones. We further identified epistatic interactions responsible for the observed allometric variation. Forty loci that interact epistatically with the detected rQTLs were identified. We demonstrate how these epistatic interactions differentially affect the body size variance and the covariance of traits with body size. We conclude that epistasis, by differentially affecting both the canalization and mean values of the traits of a pleiotropic domain, causes variation in the covariance structure. Variation in pleiotropy maintains evolvability of the genetic architecture, in particular the evolvability of its modular organization.

Were neandertal and modern human cranial differences produced by natural selection or genetic drift?

Most evolutionary explanations for cranial differences between Neandertals and modern humans emphasize adaptation by natural selection. Features of the crania of Neandertals could be adaptations to the glacial climate of Pleistocene Europe or to the high mechanical strains produced by habitually using the front teeth as tools, while those of modern humans could be adaptations for articulate speech production. A few researchers have proposed non-adaptive explanations. These stress that isolation between Neandertal and modern human populations would have lead to cranial diversification by genetic drift (chance changes in the frequencies of alleles at genetic loci contributing to variation in cranial morphology). Here we use a variety of statistical tests founded on explicit predictions from quantitative- and population-genetic theory to show that genetic drift can explain cranial differences between Neandertals and modern humans. These tests are based on thirty-seven standard cranial measurements from a sample of 2524 modern humans from 30 populations and 20 Neandertal fossils. As a further test, we compare our results for modern human cranial measurements with those for a genetic dataset consisting of 377 microsatellites typed for a sample of 1056 modern humans from 52 populations. We conclude that rather than requiring special adaptive accounts, Neandertal and modern human crania may simply represent two outcomes from a vast space of random evolutionary possibilities.

Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa

Equilibrium models of isolation by distance predict an increase in genetic differentiation with geographic distance. Here we find a linear relationship between genetic and geographic distance in a worldwide sample of human populations, with major deviations from the fitted line explicable by admixture or extreme isolation. A close relationship is shown to exist between the correlation of geographic distance and genetic differentiation (as measured by F(ST)) and the geographic pattern of heterozygosity across populations. Considering a worldwide set of geographic locations as possible sources of the human expansion, we find that heterozygosities in the globally distributed populations of the data set are best explained by an expansion originating in Africa and that no geographic origin outside of Africa accounts as well for the observed patterns of genetic diversity. Although the relationship between F(ST) and geographic distance has been interpreted in the past as the result of an equilibrium model of drift and dispersal, simulation shows that the geographic pattern of heterozygosities in this data set is consistent with a model of a serial founder effect starting at a single origin. Given this serial-founder scenario, the relationship between genetic and geographic distance allows us to derive bounds for the effects of drift and natural selection on human genetic variation.

Multivariate apportionment of global human craniometric diversity

Extensive research in human genetics on presumably neutral loci has shown that the overwhelming majority of human diversity is found among individuals within local populations. Previous apportionments of craniometric diversity are similar to these genetic apportionments, implying that interregionally differing selection pressures have played a limited role in producing contemporary human cranial diversity. Here we show that when cranial variation is independently partitioned using principal components analysis of Mosimann shape variables, some estimates of interregional craniometric differentiation are larger than those for most genetic loci. These estimates are similar to estimates for genetic loci where interregionally differing selection pressures are implicated in producing the observed patterns of variation within and among regions. These results suggest that differences among regions in at least some cranial features, particularly in the nasal region, are in part the product of interregionally differing selection pressures. Moreover, these results have implications for assessing the biological affinities of prehistoric samples that are temporally separated from contemporary human reference populations.

Detecting interregionally diversifying natural selection on modern human cranial form by using matched molecular and morphometric data

This comparison of morphological and neutral genetic variation in 10 human populations was designed to test a neutral hypothesis of cranial evolution in living and recent humans and to explain deviations from neutrality where detected. Overall, among-population differences in extant Homo sapiens cranial morphology are proportional to among-population differences in neutral molecular characteristics. For most of the populations studied, cranial morphology varies among regions in a manner consistent with neutral expectations. Removal of the effects of shared population history and structure by using the partial Mantel's test, however, does not remove the correlation between some aspects of cranial morphology and a measure of coldness of climate. The excess differentiation is most apparent in those population comparisons that involve a Siberian population living in an extremely cold environment. This finding suggests the action of natural selection, associated with regional variation in temperature, leading to among-population differentiation in excess of neutral expectations for some cranial dimensions. Those dimensions reflect the breadth of the skull, cranial vault size and shape, and aspects of nasal morphology. Although morphology for most of the world appears to vary among populations in accordance with neutral expectations in the context of population structure and history, morphology of the Siberian population appears to have undergone adaptation by natural selection.

Excavating Y-chromosome haplotype strata in Anatolia

Analysis of 89 biallelic polymorphisms in 523 Turkish Y chromosomes revealed 52 distinct haplotypes with considerable haplogroup substructure, as exemplified by their respective levels of accumulated diversity at ten short tandem repeat (STR) loci. The major components (haplogroups E3b, G, J, I, L, N, K2, and R1; 94.1%) are shared with European and neighboring Near Eastern populations and contrast with only a minor share of haplogroups related to Central Asian (C, Q and O; 3.4%), Indian (H, R2; 1.5%) and African (A, E3*, E3a; 1%) affinity. The expansion times for 20 haplogroup assemblages was estimated from associated STR diversity. This comprehensive characterization of Y-chromosome heritage addresses many multifaceted aspects of Anatolian prehistory, including: (1) the most frequent haplogroup, J, splits into two sub-clades, one of which (J2) shows decreasing variances with increasing latitude, compatible with a northward expansion; (2) haplogroups G1 and L show affinities with south Caucasus populations in their geographic distribution as well as STR motifs; (3) frequency of haplogroup I, which originated in Europe, declines with increasing longitude, indicating gene flow arriving from Europe; (4) conversely, haplogroup G2 radiates towards Europe; (5) haplogroup E3b3 displays a latitudinal correlation with decreasing frequency northward; (6) haplogroup R1b3 emanates from Turkey towards Southeast Europe and Caucasia and; (7) high resolution SNP analysis provides evidence of a detectable yet weak signal (<9%) of recent paternal gene flow from Central Asia. The variety of Turkish haplotypes is witness to Turkey being both an important source and recipient of gene flow.

Independent and linked migrants: determinants of African American interstate migration

"Through kinship and other links to destinations, many African American interstate migrants in the United States join other people in destination households. These ¿linked' migrants contrast to ¿independent migrants' who move as individuals or intact groups and set up their own households at the destination. Using U.S. Census Public Use Micro Sample data, this paper first shows that, in the 1985-90 period, about 45 percent of all Black interstate migrants were independent, compared to 38 percent who were linked to housing at the destination and 17 percent who moved into group quarters. Second, a multinomial logit model, incorporating individual and state-level variables, is specified that contrasts the determinants of independent and linked migration.... It is concluded that the understanding of Black migration must take into account a variety of factors beyond traditional labor market conditions, including links to the destination and individual housing circumstances."

The development of migration expectations: changes throughout the lifecourse

This research examines how the relative influence of factors that determine an individual's expectations of moving varies throughout the lifecourse. Factors representing personal characteristics, ties to origin, and ties to potential destinations were used to discriminate between expected migrants and nonmigrants in three population subgroups representing different stages of life. These subgroups included a general adult population, a preretirement population, and an elderly population. The results indicate that personal characteristics are most influential among preretirement persons, ties to origin are most influential among the general and elderly populations, and ties to potential destinations are influential among all population subgroups.

An AID analysis of factors associated with contemporary net migration in the United States

The authors use Automatic Interaction Detector (AID) analysis to examine variables associated with net U.S. migration by county for the period 1970-1975. "The AID analyses also reveal various ways in which several factors combine to characterize places of substantial growth or decline through migration." Government employment is identified as a major factor associated with net positive migration