Genetic correlations between sides and heritability of asymmetry for nonmetric traits in rhesus macaques on Cayo Santiago

The use of nonmetric traits for estimation of biological distance is a long-standing practice in biological anthropology. Nonmetric traits can be scored using either the individual or the side of the individual as the unit of measure. If sides of the individual are genetically correlated the use of sides would produce redundant genetic information. For this reason, Korey (Am. J. Phys. Anthropol. 53:19-23, 1980) argues for the use of individuals as the unit of measure for nonmetric traits. Ossenberg (Am. J. Phys, Anthropol. 54:471-479, 1981), however, argues that bilateral occurrence of nonmetric traits indicates greater genetic liability for the trait and that therefore the sides are the more biologically correct unit of measure. Genetic correlations for 13 cranial nonmetric traits are estimated for a sample of rhesus macaque skeletons from Cayo Santiago. In addition, heritability of asymmetry is estimated for these 13 traits as a test of Ossenberg's contention that asymmetry is genetically influenced. Significant genetic correlations between sides support Korey's contention that nonmetric traits should be scored by individual. Only two asymmetry heritabilities were significantly different from zero, providing no significant support for Ossenberg's contention that asymmetry is genetically determined. Our results support the theory that asymmetry represents a measure of the ability of an organism to buffer stresses. Therefore, a measure of the heritability of asymmetry is a measure of the heritability of the ability to buffer stresses. This ability does not appear to be heritable in this sample.     

The incidence of bilateral nonmetric skeletal traits: A reanalysis of sampling procedures

For reporting the incidence of bilateral skeletal traits, the choice between sampling statistics depends upon more than their relative efficiencies. Of overreaching importance are the fundamentally different assumptions about the genetic significance of bilateral asymmetry represented by the two principal sampling approaches. Sampling by side is consistent with the premise that trait expression on each side reflects an additive component of genetic variation. Implicit in sampling by individual, by contrast, is the proposition that asymmetries of expression result chiefly from developmental noise. The pattern of age-regression indicated for many of these traits, suggesting a transient developmental role for unilateral expression, supports both this latter view and the thesis of stress asymmetry. Given this pattern, furthermore, incorporating unadjusted trait frequencies into divergence statistics would seem injudicious.     

The use of inverse sine transformations in the analysis of non-metric cranial data

In recent years osteologists have frequently used non-metric (dichotomous) cranial data to measure biological distance between skeletal samples of Homo sapiens. Applying methods used earlier by biologists, these workers begin by attempting to stabilize the variance of the measures used by transforming the observed trait frequencies using some type of inverse sine transformation. The frequently used Grewal-Smith transformation doesn't work well for small samples of the size often considered by osteologists. As a consequence the mean measure of divergence between populations determined by this method is strongly influenced by a bias which depends on sample size. This paper compares several transformations in terms of how close the actual variance of the transformed frequency corresponds to its nominal value. It is suggested that the traditional (Grewal-Smith) inverse sine transformation not be used, and several alternatives are considered.     

A study of non-metric (qualitative) variation in Gujarati crania

Three hundred and seventy adult skulls (284 crania of unknown sex, 58 males and 28 females) from Gujarat State of India were examined for the incidence of non-metric variants and compared with other populations to establish the distance between them. In general the Gujarati incidences are of similar order to those in other series. The mean measures of divergence between Gujarati and other populations were all statistically significant (P less than 0.001). The Gujarati differed most from Australian Aborigines, but only slightly from the Burma, Punjab and Egypt samples. From the same material side and sex dimorphism was also tested to ascertain that how far sides and sexes can be pooled in Indian sample for making comparison between populations. In Gujarati population out of 22 cranial variants only four show sex difference and in case of bilateral traits, none of the variant has shown significant (P less than 0.05) side to side difference.     

An argument for the use of total side frequencies of bilateral nonmetric skeletal traits in population distance analysis: The regression of symmetry on incidence

In population studies based on frequencies of bilateral nonmetric skeletal traits, the choice between sampling by individuals or by sides should depend less on the exigencies imposed by fragmentary remains than on fundamental assumptions about the biological meaning of symmetry/asymmetry. Though the latter has been interpreted in various ways, little attention has focused on the possibility that bilateral correlation is meaningful in quantifying genetic liability for a trait. Analysis of two independent mandibular features, mylohyoid bridge and suppressed third molar, in Indian and Eskimo population samples (total N ? 1,200) reveals a statistically significant pattern of increasing bilateral occurrence with increasing population incidence. This pattern is consistent with the theory that liability for a “quasi-continuous” variant is normally distributed with constant increment between thresholds on an underlying scale. According to theory, phenotypes with more pronounced expression (bilateral occurrence) have greater genetic potential than those with less pronounced expression (unilateral occurrence). Therefore, scoring traits in total left and right sides, by giving greater weight to bilaterally affected individuals, may provide a better estimate of the liability for the trait in the population. Viewed in a theoretical context broader than that of the sampling debate, this pattern of positive regression of symmetry on incidence means that prevalence of unilateral occurrence probably cannot be used to assess the relative strength of genetic versus nongenetic control of threshold variants.     

An entropy-based measure for QTL mapping using extreme samples of population

Quantitative trait locus (QTL) mapping can be accomplished through the method of selective genotyping, which is based on the differences of frequencies between an upper sample and a lower sample in population. However, amplifying the differences in marker allele frequencies in extreme samples may increase the probability for QTL mapping. Shannon entropy, which is a nonlinear function of allele frequencies, can be used to amplify the differences in marker allele frequencies. In this paper, we present a novel measure for linkage disequilibrium (LD) between a marker and single QTL, that is based on the comparison of the entropy and conditional entropy in a marker in extreme samples of population. This measure of LD between the marker and the trait locus can be used when the marker allele frequencies are known in the extreme samples of a population. We investigate the mapping performance in both analytic and simulation scenarios of a single QTL linked to a single marker. Our results show that the measure has very reasonable performance. In addition, a simulation study is performed on the basis of the haplotype frequencies of 10 SNPs of angiotensin-I converting enzyme (ACE) genes.     

Ancient teeth and modern human origins: an expanded comparison of African Plio-Pleistocene and recent world dental samples

Previous research by the first author revealed that, relative to other modern peoples, sub-Saharan Africans exhibit the highest frequencies of ancestral (or plesiomorphic) dental traits and, thus, appear to be least derived dentally from an ancestral hominin state. This determination, in conjunction with various other lines of dental morphological evidence, was interpreted to be supportive of an African origin for modern humans. The present investigation expands upon this work by using: 1) direct observations of fossil hominin teeth, rather than data gleaned from published sources, 2) a single morphological scoring system (the Arizona State University Dental Anthropology System) with consistent trait breakpoints, and 3) data from larger and more varied modern human comparative samples. As before, a multivariate distance statistic, the mean measure of divergence, was used to assess diachronic phenetic affinities among the Plio-Pleistocene hominins and modern humans. The present study also employed principal components analysis on dental trait frequencies across samples. Both methods yielded similar results, which support the previous findings; that is, of all modern human samples, sub-Saharan Africans again exhibit the closest phenetic similarity to various African Plio-Pleistocene hominins-through their shared prevalence of morphologically complex crown and root traits. The fact that sub-Saharan Africans express these apparently plesiomorphic characters, along with additional information on their affinity to other modern populations, evident intra-population heterogeneity, and a world-wide dental cline emanating from the sub-continent, provides further evidence that is consistent with an African origin model.     

Comparison of linkage-disequilibrium methods for localization of genes influencing quantitative traits in humans

Linkage disequilibrium has been used to help in the identification of genes predisposing to certain qualitative diseases. Although several linkage-disequilibrium tests have been developed for localization of genes influencing quantitative traits, these tests have not been thoroughly compared with one another. In this report we compare, under a variety of conditions, several different linkage-disequilibrium tests for identification of loci affecting quantitative traits. These tests use either single individuals or parent-child trios. When we compared tests with equal samples, we found that the truncated measured allele (TMA) test was the most powerful. The trait allele frequencies, the stringency of sample ascertainment, the number of marker alleles, and the linked genetic variance affected the power, but the presence of polygenes did not. When there were more than two trait alleles at a locus in the population, power to detect disequilibrium was greatly diminished. The presence of unlinked disequilibrium (D'*) increased the false-positive error rates of disequilibrium tests involving single individuals but did not affect the error rates of tests using family trios. The increase in error rates was affected by the stringency of selection, the trait allele frequency, and the linked genetic variance but not by polygenic factors. In an equilibrium population, the TMA test is most powerful, but, when adjusted for the presence of admixture, Allison test 3 becomes the most powerful whenever D'*>.15.     

Fluctuating asymmetry,fecundity and development time in Drosophila: is there an association under optimal and stress conditions?

A number of studies have reported a significant negative association between fluctuating asymmetry (FA) of bilateral morphological traits and individual fitness traits, but almost all of these are unreplicated and based on small sample sizes using single trait estimates of FA. We therefore tested if there was a relationship between the FA of five bilateral traits and fecundity and development time in Drosophila in a multiple replicated experimental design. Stressed treatments were included to increase the variability of traits and to test whether associations among traits were affected by changes in the environment. Significant positive relationships were found between the size of wing characters and mean fecundity for the 5‐day period and this relationship tended to be stronger in the stress treatments. No association was found between FA and mean fecundity for any of the traits measured. Similarly, a significant positive relationship was detected between wing trait size and development time but no association was detected between trait FA and development time. There were no differences between mean fecundity or development time of extreme asymmetry phenotypes compared with modal phenotypes. These results are discussed with reference to suggestions in the literature that FA can be used to estimate individual fitness.     

Female mate choice of male signals is unlikely to promote ecological adaptation in Enchenopa treehoppers (Hemiptera: Membracidae)

A key question in speciation research is how ecological and sexual divergence arise and interact. We tested the hypothesis that mate choice causes local adaptation and ecological divergence using the rationale that the performance~signal trait relationship should parallel the attractiveness~signal trait relationship. We used female fecundity as a measure of ecological performance. We used a species in the Enchenopa binotata treehopper complex, wherein speciation involves adaptation to novel environments and divergence in sexual communication. We used a full‐sibling, split‐family rearing design to estimate genetic correlations (r_G) between fecundity and signal traits, and compared those relationships against population‐level mate preferences for the signal traits. Animal model estimates for r_G between female fecundity and male signal traits overlapped zero—rejecting the hypothesis—but could reflect sample size limitations. The magnitude of r_G correlated with the strength of the mate preferences for the corresponding signal traits, especially for signal frequency, which has the strongest mate preference and the most divergence in the complex. However, signal frequencies favored by the population‐level mate preference are not associated with high fecundity. Therefore, mate preferences do not appear to have been selected to favor high‐performance genotypes. Our findings suggest that ecological and sexual divergence may arise separately, but reinforce each other, during speciation.     

CWM and Rao’s quadratic diversity: a unified framework for functional ecology

Assessing the effects of environmental constraints on community structure often relies on methods that consider changes in species functional traits in response to environmental processes. Various indices have been proposed to measure relevant aspects of community trait composition from different viewpoints and perspectives. Among these, the ‘community-weighted mean trait value’ (CWM) and the Rao coefficient have been widely used in ecological research for summarizing different facets of functional composition and diversity. Analyzing changes in functional diversity of bee communities along a post-fire successional gradient in southern Switzerland we show that these two measures may be used to describe two complementary aspects of community structure, such as the mean and the dispersion of functional traits within a given species assemblage. While CWM can be adequately used to summarize shifts in mean trait values within communities due to environmental selection for certain functional traits, the Rao coefficient can be effectively applied to analyze patterns of trait convergence or divergence compared to a random expectation.     

Fluctuating asymmetry in the otoliths of larval fish as an indicator of condition: conceptual and methodological aspects

Developmental instability, measured as fluctuating asymmetry (FA) in bilateral traits, has been used widely as an indicator of genetic or environmental stress in a variety of plant and animal taxa. FA arises as small deviations from perfect bilateral symmetry which reflect 'mistakes' in developmental processes resulting from the inability of the genotype to buffer itself effectively against environmental perturbations. Recently, it has been proposed that FA in the otoliths can be used as an indicator of condition in larval fish. This paper reviews the conceptual and methodological aspects of FA relevant to its potential use as a measure of well-being. Its simplicity makes FA an attractive tool to measure developmental precision and condition. However, there are several pitfalls, such as measuring error or potentially size dependence. Subtle asymmetries, i.e. between sides variation of a trait at the individual level, may not always be indicative of condition and should be interpreted with caution. The past-growth record of otoliths may provide a powerful means of studying the development of asymmetries at the level of the individual.     

Sexual selection and population divergence II. Divergence in different sexual traits and signal modalities in field crickets (Teleogryllus oceanicus)

Sexual selection can target many different types of traits. However, the relative influence of different sexually selected traits during evolutionary divergence is poorly understood. We used the field cricket Teleogryllus oceanicus to quantify and compare how five traits from each of three sexual signal modalities and components diverge among allopatric populations: male advertisement song, cuticular hydrocarbon (CHC) profiles and forewing morphology. Population divergence was unexpectedly consistent: we estimated the among‐population (genetic) variance‐covariance matrix, D , for all 15 traits, and D_max explained nearly two‐thirds of its variation. CHC and wing traits were most tightly integrated, whereas song varied more independently. We modeled the dependence of among‐population trait divergence on genetic distance estimated from neutral markers to test for signatures of selection versus neutral divergence. For all three sexual trait types, phenotypic variation among populations was largely explained by a neutral model of divergence. Our findings illustrate how phenotypic integration across different types of sexual traits might impose constraints on the evolution of mating isolation and divergence via sexual selection.     

AnthropMMD: An R package with a graphical user interface for the mean measure of divergence

The mean measure of divergence is a dissimilarity measure between groups of individuals described by dichotomous variables. It is well suited to datasets with many missing values, and it is generally used to compute distance matrices and represent phenograms. Although often used in biological anthropology and archaeozoology, this method suffers from a lack of implementation in common statistical software. A package for the R statistical software, AnthropMMD, is presented here. Offering a dynamic graphical user interface, it is the first one dedicated to Smith's mean measure of divergence. The package also provides facilities for graphical representations and the crucial step of trait selection, so that the entire analysis can be performed through the graphical user interface. Its use is demonstrated using an artificial dataset, and the impact of trait selection is discussed. Finally, AnthropMMD is compared to three other free tools available for calculating the mean measure of divergence, and is proven to be consistent with them.     

Microevolution of African American dental morphology

The African American (AA) gene pool is primarily the result of gene flow between two biologically disparate groups: West Africans (WA) and Americans of western European descent (EA). This research utilizes characteristics of dental morphology to trace genetic relationships among WA, western Europeans (EU), AA, and European Americans. Dental morphological traits are useful for this purpose because they are heritable, do not remodel during life (although they can be lost to wear or pathology), and can be compared equally among samples from past and present populations. The results of this research provide new information about human microevolution through time and space in a biocultural setting. The mean measure of divergence is used to analyze dental morphological data from 1,265 individuals in 25 samples grouped by ancestry and time. Three hypotheses associated with admixture in AA are tested. When compared with known history, results from dental morphological data are equivocal in documenting admixture in AA. Dental morphological traits do appear to reflect admixture in AA. However, changes in trait frequencies do not closely correspond with important cultural events and trends such as the institutionalized racism of the Civil War and Jim Crow era. Results are mixed concerning whether AA with greater admixture were more likely to take part in the Great Migration to southern urban centers and to the North.     

A Unified Method for Detecting Secondary Trait Associations with Rare Variants: Application to Sequence Data

Next-generation sequencing has made possible the detection of rare variant (RV) associations with quantitative traits (QT). Due to high sequencing cost, many studies can only sequence a modest number of selected samples with extreme QT. Therefore association testing in individual studies can be underpowered. Besides the primary trait, many clinically important secondary traits are often measured. It is highly beneficial if multiple studies can be jointly analyzed for detecting associations with commonly measured traits. However, analyzing secondary traits in selected samples can be biased if sample ascertainment is not properly modeled. Some methods exist for analyzing secondary traits in selected samples, where some burden tests can be implemented. However p-values can only be evaluated analytically via asymptotic approximations, which may not be accurate. Additionally, potentially more powerful sequence kernel association tests, variable selection-based methods, and burden tests that require permutations cannot be incorporated. To overcome these limitations, we developed a unified method for analyzing secondary trait associations with RVs (STAR) in selected samples, incorporating all RV tests. Statistical significance can be evaluated either through permutations or analytically. STAR makes it possible to apply more powerful RV tests to analyze secondary trait associations. It also enables jointly analyzing multiple cohorts ascertained under different study designs, which greatly boosts power. The performance of STAR and commonly used RV association tests were comprehensively evaluated using simulation studies. STAR was also implemented to analyze a dataset from the SardiNIA project where samples with extreme low-density lipoprotein levels were sequenced. A significant association between LDLR and systolic blood pressure was identified, which is supported by pharmacogenetic studies. In summary, for sequencing studies, STAR is an important tool for detecting secondary-trait RV associations.     

Directional asymmetry and the measurement of developmental instability

Three widely used methods of estimating fluctuating asymmetry may yield serious overestimates if directional asymmetry is present. When two sides of a bilateral trait grow at different rates, then the asymmetry variance (Var[l-r]) increases with size, even when developmental noise is nil. But the residual variance around a population's mean developmental trajectory is invariant with respect to size. Thus, it can be used as a measure of developmental instability. We introduce a measure of developmental instability, the residual variance (s²δ), obtainable from either a major axis regression, which is equivalent to a principal component analysis on l and r, or a general structural model. This residual variance can be estimated from directionally asymmetric or even antisymmetric traits. We present examples of developmental instability estimated from directionally asymmetric mandibles (house mouse) and leaves (soybean), and antisymmetric claws (fiddler crab).     

Who were the ancient Egyptians? Dental affinities among Neolithic through postdynastic peoples

Qualitative and quantitative methods are employed to describe and compare up to 36 dental morphological variants in 15 Neolithic through Roman-period Egyptian samples. Trait frequencies are determined, and phenetic affinities are calculated using the mean measure of divergence and Mahalanobis D2 statistics for discrete traits; the most important traits in generating this intersample variation are identified with correspondence analysis. Assuming that the samples are representative of the populations from which they derive, and that phenetic similarity provides an estimate of genetic relatedness, these affinities are suggestive of overall population continuity. That is, other than a few outliers exhibiting extreme frequencies of nine influential traits, the dental samples appear to be largely homogenous and can be characterized as having morphologically simple, mass-reduced teeth. These findings are contrasted with those resulting from previous skeletal and other studies, and are used to appraise the viability of five Egyptian peopling scenarios. Specifically, affinities among the 15 time-successive samples suggest that: 1) there may be a connection between Neolithic and subsequent predynastic Egyptians, 2) predynastic Badarian and Naqada peoples may be closely related, 3) the dynastic period is likely an indigenous continuation of the Naqada culture, 4) there is support for overall biological uniformity through the dynastic period, and 5) this uniformity may continue into postdynastic times.     

The phenetic distances of the living Druze from other human populations suggest a major genetic drift from the Western Eurasian ancestral category

The objectives were to determine the expression frequency and sexual dimorphism of 16 non-metric crown traits on the sample of permanent dentitions of the living Druze population (a Near Eastern genetic isolate) in Jordan, and to assess the biological affinity of this sample to 21 regional groups, and to the living general Jordanian population, based on these traits. Druze schoolchildren (46 males, 40 females; mean age = 16.0, sd = 0.5 years) were studied in 2011. The traits were classified using the Arizona State University dental anthropology system, counted with the individual count method, and dichotomized according to the criteria of Scott and Turner for the purpose of group comparisons. Fisher's exact test for dichotomized scores was used to assess sexual dimorphism in these traits. Smith's mean measure of divergence was used to measure all pairwise distance values among the groups. Sexual dimorphism was found in five traits (i.e., UI2 interruption grooves, 3-cusped UM2, UM1 Carabelli's tubercle/cusp, 4-cusp LM1, and LM2 Y-groove pattern). This study revealed that the dental pattern of living Druze, which is similar to that of the general Jordanian population, is sufficiently distinct from the Western Eurasian pattern and all other known dental patterns to form a distinct dental pattern for the regional group or subcategory to which these two populations belong. Moreover, the relatively large distance values of the living Druze and Jordanians from the other world groups considered, including the Western Eurasian groups, suggest a similar major genetic difference of these two populations from the Western Eurasian ancestry.     

Q(St) meets the G matrix: the dimensionality of adaptive divergence in multiple correlated quantitative traits.

The Q(ST)-F(ST) comparison has become an increasingly common method for inferring adaptive quantitative trait divergence among populations. For cases in which there is divergence in multiple traits, most studies have applied the method by performing multiple univariate Q(ST)-F(ST) comparisons. However, because traits are often genetically correlated, such univariate analyses are likely to paint a simplified picture of adaptive divergence. Here we show how the multivariate analogue of Q(ST), F(STq), which accounts for genetic correlations among traits, can be used to supply a more detailed picture of multitrait divergence. We apply the method to naturally occurring genetic variation for a suite of sexually selected display traits in Drosophila serrata. The analyses suggest the operation of divergent multivariate selection that has influenced multiple independent axes of genetic variance in a sex-specific manner. Finally, we show how a comparison of the components of F(STq), the average within and among population genetic variance-covariance matrices, G(W) and G(B), can be used as an additional test of the null expectation of neutral divergence, and allows for an investigation of whether natural populations have diverged along major or minor axes of genetic variance.