Genetic correlations between cartilage regeneration and degeneration reveal an inverse relationship.

OBJECTIVE: The etiology of osteoarthritis (OA) is unknown, however, there appears to be a significant contribution from genetics. We have identified recombinant inbred strains of mice derived from LG/J (large) and SM/J (small) strains that vary significantly in their ability to repair articular cartilage and susceptibility to post-traumatic OA due to their genetic composition. Here, we report cartilage repair phenotypes in the same strains of mice in which OA susceptibility was analyzed previously, and determine the genetic correlations between phenotypes. DESIGN: We used 12 recombinant inbred strains, including the parental strains, to test three phenotypes: ear-wound healing (n = 263), knee articular cartilage repair (n = 131), and post-traumatic OA (n = 53) induced by the surgical destabilization of the medial meniscus (DMM). Genetic correlations between various traits were calculated as Pearson's correlation coefficients of strain means. RESULTS: We found a significant positive correlation between ear-wound healing and articular cartilage regeneration (r = 0.71; P = 0.005). We observed a strong inverse correlation between articular cartilage regeneration and susceptibility to OA based on maximum (r = -0.54; P = 0.036) and summed Osteoarthritis Research Society International (OARSI) scores (r = -0.56; P = 0.028). Synovitis was not significantly correlated with articular cartilage regeneration but was significantly positively correlated with maximum (r = 0.63; P = 0.014) and summed (r = 0.70; P = 0.005) OARSI scores. Ectopic calcification was significantly positively correlated with articular cartilage regeneration (r = 0.59; P = 0.021). CONCLUSIONS: Using recombinant inbred strains, our study allows, for the first time, the measurement of genetic correlations of regeneration phenotypes with degeneration phenotypes, characteristic of OA (cartilage degeneration, synovitis). We demonstrate that OA is positively correlated with synovitis and inversely correlated with the ability to repair cartilage. These results suggest an addition to the risk paradigm for OA from a focus on degeneration to regeneration.

RNA-seq analysis of chondrocyte transcriptome reveals genetic heterogeneity in LG/J and SM/J murine strains.

OBJECTIVE: To investigate the transcriptomic differences in chondrocytes obtained from LG/J (large, healer) and SM/J (small, non-healer) murine strains in an attempt to discern the molecular pathways implicated in cartilage regeneration and susceptibility to osteoarthritis (OA). DESIGN: We performed RNA-sequencing on chondrocytes derived from LG/J (n = 16) and SM/J (n = 16) mice. We validated the expression of candidate genes and compared single nucleotide polymorphisms (SNPs) between the two mouse strains. We also examined gene expression of positional candidates for ear pinna regeneration and long bone length quantitative trait loci (QTLs) that display differences in cartilaginous expression. RESULTS: We observed a distinct genetic heterogeneity between cells derived from LG/J and SM/J mouse strains. We found that gene ontologies representing cell development, cartilage condensation, and regulation of cell differentiation were enriched in LG/J chondrocytes. In contrast, gene ontologies enriched in the SM/J chondrocytes were mainly related to inflammation and degeneration. Moreover, SNP analysis revealed that multiple validated genes vary in sequence between LG/J and SM/J in coding and highly conserved noncoding regions. Finally, we showed that most QTLs have 20-30% of their positional candidates displaying differential expression between the two mouse strains. CONCLUSIONS: While the enrichment of pathways related to cell differentiation, cartilage development and cartilage condensation infers superior healing potential of LG/J strain, the enrichment of pathways related to cytokine production, immune cell activation and inflammation entails greater susceptibility of SM/J strain to OA. These data provide novel insights into chondrocyte transcriptome and aid in identification of the quantitative trait genes and molecular differences underlying the phenotypic differences associated with individual QTLs.

Rate of evolutionary change in cranial morphology of the marsupial genus Monodelphis is constrained by the availability of additive genetic variation.

We tested the hypothesis that the rate of marsupial cranial evolution is dependent on the distribution of genetic variation in multivariate space. To do so, we carried out a genetic analysis of cranial morphological variation in laboratory strains of Monodelphis domestica and used estimates of genetic covariation to analyse the morphological diversification of the Monodelphis brevicaudata species group. We found that within-species genetic variation is concentrated in only a few axes of the morphospace and that this strong genetic covariation influenced the rate of morphological diversification of the brevicaudata group, with between-species divergence occurring fastest when occurring along the genetic line of least resistance. Accounting for the geometric distribution of genetic variation also increased our ability to detect the selective regimen underlying species diversification, with several instances of selection only being detected when genetic covariances were taken into account. Therefore, this work directly links patterns of genetic covariation among traits to macroevolutionary patterns of morphological divergence. Our findings also suggest that the limited distribution of Monodelphis species in morphospace is the result of a complex interplay between the limited dimensionality of available genetic variation and strong stabilizing selection along two major axes of genetic variation.

Fine-mapping quantitative trait loci affecting murine external ear tissue regeneration in the LG/J by SM/J advanced intercross line.

External ear hole closure in LG/J mice represents a model of regenerative response. It is accompanied by the formation of a blastema-like structure and the re-growth of multiple tissues, including cartilage. The ability to regenerate tissue is heritable. An F34 advanced intercross line of mice (Wustl:LG,SM-G34) was generated to identify genomic loci involved in ear hole closure over a 30-day healing period. We mapped 19 quantitative trait loci (QTL) for ear hole closure. Individual gene effects are relatively small (0.08 mm), and most loci have co-dominant effects with phenotypically intermediate heterozygotes. QTL support regions were limited to a median size of 2 Mb containing a median of 19 genes. Positional candidate genes were evaluated using differential transcript expression between LG/J and SM/J healing tissue, function analysis and bioinformatic analysis of single-nucleotide polymorphisms in and around positional candidate genes of interest. Analysis of the set of 34 positional candidate genes and those displaying expression differences revealed over-representation of genes involved in cell cycle regulation/DNA damage, cell migration and adhesion, developmentally related genes and metabolism. This indicates that the healing phenotype in LG/J mice involves multiple physiological mechanisms.

Relationship of age and body mass index to the expression of obesity and osteoarthritis-related genes in human meniscus.

OBJECTIVE: Aging and obesity contribute to the initiation and progression of osteoarthritis with little information on their relation to gene expression in joint tissues, particularly the meniscus. Here, we test the hypothesis that patient age and body mass index (BMI) correlate with the expression of osteoarthritis- and obesity-related gene signatures in the meniscus. DESIGN: Meniscus was obtained from patients (N=68) undergoing arthroscopic partial meniscectomy. The mRNA expression of 24 osteoarthritis-related and 4 obesity-related genes in meniscus was assessed by quantitative real-time PCR. The relationship between gene expression and patient age and BMI was analyzed using Spearman's rank-order correlation. Hierarchical cluster dendrogram and heat map were generated to study inter-gene associations. RESULTS: Age was negatively correlated (P<0.05) with the expression of MMP-1 (r=-0.447), NFκB2 (r=-0.361), NFκBIA (r=-0.312), IκBA (r=-0.308), IL-8 (r=-0.305), ADAMTS-4 (r=-0.294), APLN (apelin) (r=-0.250) and IL-6 (r=-0.244). Similarly, BMI was negatively correlated with the expression of APLN (r=-0.328), ACAN (r=-0.268) and MMP-1 (r=-0.261). After adjusting for the correlation between age and BMI (r=0.310; P=0.008), the only independent effect of BMI on gene expression was for APLN (r=-0.272). However, age had an independent effect on the expression on ADAMTS-4 (r=-0.253), MMP-1 (r=-0.399), IL-8 (r=-0.327), COL1A1 (r=-0.287), NFκBIA (r=-0.278), NFκB2 (r=-0.312) and IκBA (r=-0.299). The gene correlation analysis identified four clusters of potentially relevant genes: transcription factors, matrix-degrading enzymes, cytokines and chemokines, and obesity genes. CONCLUSION: Age and BMI were negatively correlated with several osteoarthritis- and obesity-related genes. Although the bulk of these changes appeared to be driven by age, expression of APLN was related to BMI. Inter-gene correlation analysis implicated a common role for strongly correlated genes. Although age-related variations in gene expression appear to be more relevant than obesity-related differences for the role of the meniscus in osteoarthritis development, further investigation into the role of APLN in meniscus and joint health is warranted.

Cartilage and bone changes during development of post-traumatic osteoarthritis in selected LGXSM recombinant inbred mice.

INTRODUCTION: Little evidence is available on the natural course of osteoarthritis (OA) development and the genes that protect and predispose individuals to it. This study was designed to compare strain-dependent development of OA and its association with tissue regeneration in mice. Two recombinant inbred lines LGXSM-6 and LGXSM-33 generated from LG/J and SM/J intercross were used. Previous studies indicated that LGXSM-6 can regenerate both articular cartilage and ear hole punch while LGXSM-33 cannot. METHODS: Transection of the medial meniscotibial ligament was performed on 10-week-old male mice to induce OA. Cartilage damage was analyzed by histology and bone morphology was evaluated using micro-computed tomography (CT). Ear punches were performed and evaluated by measurement of residual hole diameter. RESULTS: Cartilage analysis showed that LGXSM-33 developed a significantly higher grade of OA than LGXSM-6. Bone analysis showed that LGXSM-33 had substantial subchondral bone and trabecular bone thickening 8 weeks post-surgery, while LGXSM-6 showed bone loss over time. We also confirmed that LGXSM-6 can heal ear tissues significantly better than LGXSM-33. CONCLUSIONS: OA was found to be negatively correlated with the degree of tissue regeneration. LGXSM-33, a poor healer of ear tissues (and articular cartilage), developed more OA compared to LGXSM-6, which had better regenerative ability for ear tissues and articular cartilage. The phenotypic differences observed here are due to genetic differences further suggesting that similar sets of physiological processes and gene variants may mediate variation in OA development and tissue regeneration.

Healing quantitative trait loci in a combined cross analysis using related mouse strain crosses.

Inbred mouse strains MRL and LG share the ability to fully heal ear hole punches with the full range of appropriate tissues without scarring. They also share a common ancestry, MRL being formed from a multi-strain cross with two final backcrosses to LG before being inbred by brother-sister mating. Many gene-mapping studies for healing ability have been performed using these two strains, resulting in the location of about 20 quantitative trait loci (QTLs). Here, we combine two of these crosses (N = 638), MRL/lpr × C57BL/6NTac and LG/J × SM/J, in a single combined cross analysis to increase the mapping power, decrease QTL support intervals, separate multiple QTLs and establish allelic states at individual QTL. The combined cross analysis located 11 QTLs, 6 affecting only one cross (5 LG × SM and 1 MRL × B6) and 5 affecting both crosses, approximately the number of common QTLs expected given strain SNP similarity. Amongst the five QTLs mapped in both crosses, three had significantly different genetic effects, additive in one cross and over or underdominant in the other. It is possible that allelic states at these three loci are different in SM and B6 because they lead to differences in dominance interactions with the LG and MRL alleles. QTL support intervals are 40% smaller in the combined cross analysis than in either of the single crosses. Combined cross analysis was successful in enhancing the interpretation of earlier QTL results for these strains.

Genotype-dependent responses to levels of sibling competition over maternal resources in mice.

Research on phenotypic plasticity has often focused on how a given genotype responds to the changing physical environments such as temperature or diet. However, for many species the social environment has an equally important role because of competition for resources. During early development, the level of competition for limited (maternally provided) resources will often depend critically on the number of siblings. Therefore, competition among siblings should drive the evolution of genes that allow flexible responses to realized levels of competition and maternal resource availability. However, it is unknown whether genetically based differences between individuals exist in their response to the social environment that affect their future development. Using a quantitative trait locus approach in an experimental population of mice we demonstrate that effects of sibling number on body weight depend on individual genotype at seven loci, over and above the general negative litter size effect. Overall, these litter size-by-genotype interactions considerably modified the degree to which increasing litter size caused reduced weight. For example at one locus this effect leads to a 7% difference in body weight at week 7 between individuals experiencing the extremes of the normal range of litter sizes in our population (five to nine litter mates). The observed interaction between genotype and the competitive environment can produce differences in body weight that are similar in magnitude to the main effect of litter size on weight. Our results show that different genotypes respond to the social environment differentially and that interaction effects of genotype with litter size can be as important as genotype-independent effects of litter size.

The effect of a population bottleneck on the evolution of genetic variance/covariance structure.

It is well known that standard population genetic theory predicts decreased additive genetic variance (V(a) ) following a population bottleneck and that theoretical models including interallelic and intergenic interactions indicate such loss may be avoided. However, few empirical data from multicellular model systems are available, especially regarding variance/covariance (V/CV) relationships. Here, we compare the V/CV structure of seventeen traits related to body size and composition between control (60 mating pairs/generation) and bottlenecked (2 mating pairs/generation; average F = 0.39) strains of mice. Although results for individual traits vary considerably, multivariate analysis indicates that V(a) in the bottlenecked populations is greater than expected. Traits with patterns and amounts of epistasis predictive of enhanced V(a) also show the largest deviations from additive expectations. Finally, the correlation structure of weekly weights is not significantly different between control and experimental lines but correlations between necropsy traits do differ, especially those involving the heart, kidney and tail length.

Neuroanatomic and behavioral traits for autistic disorders in age-specific restricted index selection mice.

The pathogenesis of neurodevelopmental disorders such as autism is believed to be influenced by interactions between genetic and environmental factors, and appropriate animal models are needed to assess the influence of such factors on relevant neurodevelopmental phenotypes. A set of inbred mouse strains (Atchley strains) including A12 (E+L0) and A22 (E-L0) were generated by age-specific restricted index selection from a baseline random-bred ICR mouse population obtained from Harlan Sprague-Dawley [Atchley et al. (1997) Genetics 146(2):629-640; Indianapolis, IN, USA). As compared with the A22 strain, A12 mice had significantly increased early (P0-P10) body weight gain with minimal changes in late (P28-P56) body weight gain. We found that these strains also differed in brain weight, brain volume, cell proliferation, and FGF-2 levels in certain brain regions. Specifically, brain weight and volume were significantly greater in A12 mice than that in A22 mice at P10 and P28. Quantitative analysis of bromodeoxyuridine (BrdU) labeling of proliferating cells showed that the number of BrdU-positive cells in the A12 strain were significantly greater in the frontal cortex and lesser in the dentate gyrus than that in the A22 strain at P28. Western blot revealed that fibroblast growth factors-2 (FGF-2), but not brain-derived neurotrophic factor (BDNF), expression was significantly increased in the frontal cortex of A12 strain at P28. Also, A12 mice exhibited decreased intra-strain social interaction and increased repetitive stereotyped behaviors at P28. Our study suggests that A12 mice may partially mimic the anatomic and behavioral traits of patients with neurodevelopmental disorders such as autism spectrum disorders, and therefore may yield insights into the developmental mechanisms involved in their pathogenesis.

Reply to Lawler: feeding competition, cooperation, and the causes of primate sociality.

This is a reply to Richard Lawler's commentary on our previous work [Lawler, 2011; this issue] in which he develops a set of operational models to test socioecological theories of the evolutionary importance of feeding competition. We strongly agree that we need to critically re-evaluate the basic assumptions of all models of primate sociality, and to verify the explanatory power of alternative models. We also feel Lawler's commentary provides an important opportunity to broaden the debate concerning the fundamental roles of cooperation, competition, and aggression in understanding primate social systems. Lawler provides a number of suggestions as to how models developed in primate socioecology might be tested. We agree with these suggestions, make further suggestions, and call for specific operational definitions so that researchers might begin to develop and test various methodologies. However, we also call for testing alternative theories. Current socioecological theory is based on the assumption that competition and positive selection is always in operation and has driven the evolution of living organisms. We believe that this "explanation of choice" often is treated as an assumed truth to which data are forced to fit, rather than being seen as a theory to be tested. Furthermore, we agree with Weiss and Buchanan [2009. The Mermaid's Tale: Four Billion Years of Cooperation in the Making of Living Things] that on ecological and developmental scales, where organisms actually live out their lives, cooperation may play a more fundamental role than competition.

Sex-specific quantitative trait loci linked to autoresuscitation failure in SWR/J mice.

Autoresuscitation (AR) is a highly conserved response among mammals, which allows survival from transient extreme hypoxia. During hypoxia, bradycardia, and hypoxic gasping develop after a brief period of hyperactivity. Normally, AR occurs if oxygen is restored during the gasping period where an initial heart rate increase is rapidly followed resumption or eupneic breathing. Humans and other mammals can survive multiple immediately repeated AR. A defective AR capacity has been implicated in Sudden Infant Death Syndrome. We had reported earlier that inbred strains of mice such as BALB/cJ could survive a characteristic number of immediately repeated AR trials, but that SWR/J mice failed to AR from a single hypoxic episode. We now report that strains closely related to SWR/J, FVB/N and SJL/J exhibit partial resuscitation defects relative to BALB/cJ or other mouse strains, establishing a genetic basis for variation in AR failure. The AR trial phenotype of BALB/cJ x SWR/J intercross F(1) and F(2) mice was consistent with BALB/cJ dominance and a discrete number of loci. Genome-wide mapping conducted with 60 intercross F(2) animals linked two loci to the number of AR trials survived, including one sex-specific locus with male expression, consistent with the observed 50% male bias for Sudden Infant Death Syndrome in humans. A locus carried on SWR/J chromosome 10 seems to be particularly important in AR failure and was confirmed in a partial consomic line. These results establish a genetic basis for AR failure phenotype in mice, with relevance to Sudden Infant Death Syndrome.

Epistasis and the evolutionary dynamics of measured genotypic values during simulated serial bottlenecks.

The evolutionary effects of epistasis have been primarily explored analytically and most empirical studies have utilized yeast, viral and bacterial populations. Empirical analyses in multi-cellular organisms are rare because of experimental constraints. Here, we report the results of a genome-wide scan for two-way epistasis in 16 traits related to body size and composition in F(2) mice from the LG/J by SM/J intercross. We analyze two-locus genotypic values at quantitative trait loci (QTL), which provides an especially detailed view of epistatic architectures, to evaluate their predicted evolutionary consequences via Monte Carlo simulations. Epistatic profiles vary, but all traits show complicated genetic architectures which are largely hidden in single locus QTL scans. On average, detected epistatic effects are comparable in size to marginal effects. Simulations demonstrate an expected preservation, and often inflation, of heritable variance across several generations of small effective population size for many identified epistatic pairs over a range of starting allele frequencies.

A search for quantitative trait loci exhibiting imprinting effects on mouse mandible size and shape.

Genomic imprinting refers to the pattern of monoallelic parent-of-origin-dependent gene expression where one of the two alleles at a locus is expressed and the other silenced. Although some genes in mice are known to be imprinted, the true scope of imprinting and its impact on the genetic architecture of a wide range of morphometric traits is mostly unknown. We therefore searched for quantitative trait loci (QTL) exhibiting imprinting effects on mandible size and shape traits in a large F(3) population of mice originating from an intercross of the LG/J (Large) and SM/J (Small) inbred strains. We discovered a total of 51 QTL affecting mandible size and shape, 6 of which exhibited differences between reciprocal heterozygotes, the usual signature of imprinting effects. However, our analysis showed that only one of these QTL (affecting mandible size) exhibited a pattern consistent with true imprinting effects, whereas reciprocal heterozygote differences in the other five all were due to maternal genetic effects. We concluded that genomic imprinting has a negligible effect on these specific morphometric traits, and that maternal genetic effects may account for many of the previously reported instances of apparent genomic imprinting.

Quantitative trait loci linked to thalamus and cortex gray matter volumes in BXD recombinant inbred mice.

To investigate whether there are separate or shared genetic influences on the development of the thalamus and cerebral cortex, we identified quantitative trait loci (QTLs) for relevant structural volumes in BXD recombinant inbred (RI) strains of mice. In 34 BXD RI strains and two parental strains (C57BL/6J and DBA/2J), we measured the volumes of the entire thalamus and cortex gray matter using point counting and Cavalieri's rule. Heritability was calculated using analysis of variance (ANOVA), and QTL analysis was carried out using WebQTL (http://www.genenetwork.org). The heritability of thalamus volume was 36%, and three suggestive QTLs for thalamus volume were identified on chromosomes 10, 11 and 16. The heritability of cortical gray matter was 43%, and four suggestive QTLs for cortex gray matter volume were identified on chromosomes 2, 8, 16 and 19. The genetic correlation between thalamus and cortex gray matter volumes was 0.64. Also, a single QTL on chromosome 16 (D16Mit100) was identified for thalamus volume, cortex gray matter volume and Morris water maze search-time preference (r=0.71). These results suggest that there are separate and shared genetic influences on the development of the thalamus and cerebral cortex.

Genetic architecture of Arabidopsis thaliana response to infection by Pseudomonas syringae.

Plant pathogens can severely reduce host yield and fitness. Thus, investigating the genetic basis of plant response to pathogens is important to further understand plant-pathogen coevolution and to improve crop production. The interaction between Arabidopsis thaliana and Pseudomonas syringae is an important model for studying the genetic basis of plant-pathogen interactions. Studies in this model have led to the discovery of many genes that differentiate a resistant from a susceptible plant. However, little is known about the genetic basis of quantitative variation in response to P. syringae. In this study, we investigate the genetic basis of three aspects of A. thaliana's response to P. syringae: symptom severity, bacterial population size and fruit production using a quantitative trait loci (QTL) analysis. We found two QTL for symptom severity and two for fruit production (possible candidate genes for observed QTL are discussed). We also found significant two-locus epistatic effect on symptom severity and fruit production. Although bacterial population size and symptom severity were strongly phenotypically correlated, we did not detect any QTL for bacterial population size. Despite the detected genetic variation observed for susceptibility, we found only a weak overall relationship between susceptibility traits and fitness, suggesting that these traits may not respond to selection.

An epistatic genetic basis for fluctuating asymmetry of tooth size and shape in mice.

Although there typically is little additive genetic variation for fluctuating asymmetry (FA), or variation in nondirectional differences between left and right sides of bilateral characters, several investigators have hypothesized that FA may have an epistatic genetic basis. We tested this hypothesis by conducting a whole genome scan of FA of size and shape of the mandibular molars in house mice from an F2 intercross population generated from crossing the Large (LG/J) and Small (SM/J) inbred strains. Although no individual genes (QTLs=quantitative trait loci) on any of the 19 autosomes significantly affected FA for centroid size, and only two affected shape FA, a number of pairwise combinations of QTLs exhibited significant epistasis for FA in both molar size and shape. The QTLs involved in these interactions differed for FA in molar size versus FA in molar shape, but their epistatic contributions to the total variance was nearly the same (about 20%) for FA in both molar characters. It was noted that the genetic architecture of FA in the molar characters, consisting of little or no additive genetic variance but an abundance of epistatic genetic variance, is consistent with that of other typical fitness components such as litter size.

Epistasis affecting litter size in mice.

Litter size is an important reproductive trait as it makes a major contribution to fitness. Generally, traits closely related to fitness show low heritability perhaps because of the corrosive effects of directional natural selection on the additive genetic variance. Nonetheless, low heritability does not imply, necessarily, a complete absence of genetic variation because genetic interactions (epistasis and dominance) contribute to variation in traits displaying strong heterosis in crosses, such as litter size. In our study, we investigated the genetic architecture of litter size in 166 females from an F2 intercross of the SM/J and LG/J inbred mouse strains. Litter size had a low heritability (h2 = 12%) and a low repeatability (r = 33%). Using interval-mapping methods, we located two quantitative trait loci (QTL) affecting litter size at locations D7Mit21 + 0 cM and D12Mit6 + 8 cM, on chromosomes 7 and 12 respectively. These QTL accounted for 12.6% of the variance in litter size. In a two-way genome-wide epistasis scan we found eight QTL interacting epistatically involving chromosomes 2, 4, 5, 11, 14, 15 and 18. Taken together, the QTL and their interactions explain nearly 49% (39.5% adjusted multiple r2) of the phenotypic variation for litter size in this cross, an increase of 36% over the direct effects of the QTL. This indicates the importance of epistasis as a component of the genetic architecture of litter size and fitness in our intercross population.

Cranial evolution in sakis (Pithecia, Platyrrhini). II: Evolutionary processes and morphological integration.

Patterns of interspecific differentiation in saki monkeys (Pithecia) are quantitatively described and possible evolutionary processes producing them are examined. The comparison of species correlation matrices to expected patterns of morphological integration reveal significant and similar patterns of development-based cranial integration among species. Aspects of the facial region are more heavily influenced by general size variation than features of the neural region. The comparison of pooled within- and between-groups V/CV matrices suggests that genetic drift might be a sufficient explanation for saki cranial evolution. Differential natural selection gradients are also reconstructed because selection may also have caused population differentiation through evolutionary time. These gradients illustrate the inherent multivariate nature of selection, being a consequence of the interaction between existing morphological integration (correlation) among traits and the action of natural selection. Yet, our attempt to interpret selection gradients in terms of their functional significance did not result in any clear association between selection and function. Perhaps this is also an indication that morphological evolution in sakis was mostly neutral.