Applying an evolutionary mismatch framework to understand disease susceptibility.

Noncommunicable diseases (NCDs) are on the rise worldwide. Obesity, cardiovascular disease, and type 2 diabetes are among a long list of "lifestyle" diseases that were rare throughout human history but are now common. The evolutionary mismatch hypothesis posits that humans evolved in environments that radically differ from those we currently experience; consequently, traits that were once advantageous may now be "mismatched" and disease causing. At the genetic level, this hypothesis predicts that loci with a history of selection will exhibit "genotype by environment" (GxE) interactions, with different health effects in "ancestral" versus "modern" environments. To identify such loci, we advocate for combining genomic tools in partnership with subsistence-level groups experiencing rapid lifestyle change. In these populations, comparisons of individuals falling on opposite extremes of the "matched" to "mismatched" spectrum are uniquely possible. More broadly, the work we propose will inform our understanding of environmental and genetic risk factors for NCDs across diverse ancestries and cultures.

Glucosylceramides impact cellulose deposition and cellulose synthase complex motility in Arabidopsis.

Cellulose is an abundant component of plant cell wall matrices, and this para-crystalline polysaccharide is synthesized at the plasma membrane by motile Cellulose Synthase Complexes (CSCs). However, the factors that control CSC activity and motility are not fully resolved. In a targeted chemical screen, we identified the alkylated nojirimycin analog N-Dodecyl Deoxynojirimycin (ND-DNJ) as a small molecule that severely impacts Arabidopsis seedling growth. Previous work suggests that ND-DNJ-related compounds inhibit the biosynthesis of glucosylceramides (GlcCers), a class of glycosphingolipid associated with plant membranes. Our work uncovered major changes in the sphingolipidome of plants treated with ND-DNJ, including reductions in GlcCer abundance and altered acyl chain length distributions. Crystalline cellulose content was also reduced in ND-DNJ-treated plants as well as plants treated with the known GlcCer biosynthesis inhibitor N-[2-hydroxy-1-(4-morpholinylmethyl)-2-phenyl ethyl]-decanamide (PDMP) or plants containing a genetic disruption in GLUCOSYLCERAMIDE SYNTHASE (GCS), the enzyme responsible for sphingolipid glucosylation that results in GlcCer synthesis. Live-cell imaging revealed that CSC speed distributions were reduced upon treatment with ND-DNJ or PDMP, further suggesting an important relationship between glycosylated sphingolipid composition and CSC motility across the plasma membrane. These results indicate that multiple interventions compromising GlcCer biosynthesis disrupt cellulose deposition and CSC motility, suggesting that GlcCers regulate cellulose biosynthesis in plants.

The effect of increased fruit and vegetable consumption on blood pressure and lipids: a pooled analysis of six randomised controlled fruit and vegetable intervention trials.

BACKGROUND: Increasing fruit and vegetable (FV) consumption is associated with reduced cardiovascular disease risk in observational studies but with little evidence from randomised controlled trials (RCTs). The impact of concurrent pharmacological therapy is unknown. OBJECTIVE: To pool data from six RCTs to examine the effect of increasing FV intake on blood pressure (BP) and lipid profile, also exploring whether effects differed by medication use. DESIGN: Across trials, dietary intake was assessed by diet diaries or histories, lipids by routine biochemical methods and BP by automated monitors. Linear regression provided an estimate of the change in lipid profile or BP associated with a one portion increase in self-reported daily FV intake, with interaction terms fitted for medication use. RESULTS: The pooled sample included a total of 554 participants (308 males and 246 females). Meta-analysis of regression coefficients revealed no significant change in either systolic or diastolic BP per portion FV increase, although there was significant heterogeneity across trials for systolic BP (I2 = 73%). Neither adjusting for change in body mass index, nor analysis according to use of anti-hypertensive medication altered the relationship. There was no significant change in lipid profile per portion FV increase, although there was a significant reduction in total cholesterol among those not on lipid-lowering therapy (P < 0.05 after Bonferroni correction). CONCLUSION: Pooled analysis of six individual FV trials showed no impact of increasing intake on BP or lipids, but there was a total cholesterol-lowering effect in those not on lipid-lowering therapy.

A novel in vitro high-content imaging assay for the prediction of drug-induced lung toxicity.

The development of inhaled drugs for respiratory diseases is frequently impacted by lung pathology in non-clinical safety studies. To enable design of novel candidate drugs with the right safety profile, predictive in vitro lung toxicity assays are required that can be applied during drug discovery for early hazard identification and mitigation. Here, we describe a novel high-content imaging-based screening assay that allows for quantification of the tight junction protein occludin in A549 cells, as a model for lung epithelial barrier integrity. We assessed a set of compounds with a known lung safety profile, defined by clinical safety or non-clinical in vivo toxicology data, and were able to correctly identify 9 of 10 compounds with a respiratory safety risk and 9 of 9 compounds without a respiratory safety risk (90% sensitivity, 100% specificity). The assay was sensitive at relevant compound concentrations to influence medicinal chemistry optimization programs and, with an accessible cell model in a 96-well plate format, short protocol and application of automated imaging analysis algorithms, this assay can be readily integrated in routine discovery safety screening to identify and mitigate respiratory toxicity early during drug discovery. Interestingly, when we applied physiologically-based pharmacokinetic (PBPK) modelling to predict epithelial lining fluid exposures of the respiratory tract after inhalation, we found a robust correlation between in vitro occludin assay data and lung pathology in vivo, suggesting the assay can inform translational risk assessment for inhaled small molecules.

CALCIUM-DEPENDENT PROTEIN KINASE32 regulates cellulose biosynthesis through post-translational modification of cellulose synthase.

Cellulose is an essential component of plant cell walls and an economically important source of food, paper, textiles, and biofuel. Despite its economic and biological significance, the regulation of cellulose biosynthesis is poorly understood. Phosphorylation and dephosphorylation of cellulose synthases (CESAs) were shown to impact the direction and velocity of cellulose synthase complexes (CSCs). However, the protein kinases that phosphorylate CESAs are largely unknown. We conducted research in Arabidopsis thaliana to reveal protein kinases that phosphorylate CESAs. In this study, we used yeast two-hybrid, protein biochemistry, genetics, and live-cell imaging to reveal the role of calcium-dependent protein kinase32 (CPK32) in the regulation of cellulose biosynthesis in A. thaliana. We identified CPK32 using CESA3 as a bait in a yeast two-hybrid assay. We showed that CPK32 phosphorylates CESA3 while it interacts with both CESA1 and CESA3. Overexpressing functionally defective CPK32 variant and phospho-dead mutation of CESA3 led to decreased motility of CSCs and reduced crystalline cellulose content in etiolated seedlings. Deregulation of CPKs impacted the stability of CSCs. We uncovered a new function of CPKs that regulates cellulose biosynthesis and a novel mechanism by which phosphorylation regulates the stability of CSCs.

Cellulose synthesis complexes are homo-oligomeric and hetero-oligomeric in Physcomitrium patens.

The common ancestor of seed plants and mosses contained homo-oligomeric cellulose synthesis complexes (CSCs) composed of identical subunits encoded by a single CELLULOSE SYNTHASE (CESA) gene. Seed plants use different CESA isoforms for primary and secondary cell wall deposition. Both primary and secondary CESAs form hetero-oligomeric CSCs that assemble and function in planta only when all the required isoforms are present. The moss Physcomitrium (Physcomitrella) patens has seven CESA genes that can be grouped into two functionally and phylogenetically distinct classes. Previously, we showed that PpCESA3 and/or PpCESA8 (class A) together with PpCESA6 and/or PpCESA7 (class B) form obligate hetero-oligomeric complexes required for normal secondary cell wall deposition. Here, we show that gametophore morphogenesis requires a member of class A, PpCESA5, and is sustained in the absence of other PpCESA isoforms. PpCESA5 also differs from the other class A PpCESAs as it is able to self-interact and does not co-immunoprecipitate with other PpCESA isoforms. These results are consistent with the hypothesis that homo-oligomeric CSCs containing only PpCESA5 subunits synthesize cellulose required for gametophore morphogenesis. Analysis of mutant phenotypes also revealed that, like secondary cell wall deposition, normal protonemal tip growth requires class B isoforms (PpCESA4 or PpCESA10), along with a class A partner (PpCESA3, PpCESA5, or PpCESA8). Thus, P. patens contains both homo-oligomeric and hetero-oligomeric CSCs.

Phalangeal curvature in a chimpanzee raised like a human: Implications for inferring arboreality in fossil hominins.

Arboreal primates such as chimpanzees exhibit pronounced curvature in their hand and foot phalanges, which is assumed to develop throughout life in response to mechanical loads produced by grasping and hanging from branches. Intriguingly, ancient fossil hominins also exhibit substantial phalangeal curvature, which, too, has been interpreted as a direct result of habitual arboreality during life. Here, we describe the phalangeal curvature of a chimpanzee who was raised during the 1930s in New York City to live much like a human, including by having very few opportunities to engage in arboreal activities. We show that the degree of hand and foot phalangeal curvature in this individual is indistinguishable from that of wild chimpanzees and distinct from humans. Thus, rather than being a direct effect of mechanical loads produced by lifetime arboreal activities, phalangeal curvature appears to be shaped largely by genetic factors. An important implication of this finding is that phalangeal curvature among fossil hominins is evidently best interpreted as a primitive trait inherited from an arboreal ancestral species rather than proof of engagement in arboreal activities during life.

Geographical analysis of the distribution of certified genetic counselors in the United States.

The number of certified genetic counselors (CGCs) in the genetic counseling workforce has increased over the past few decades as the number of training programs increases and CGCs expand into new patient-facing and non-patient-facing roles. Few studies have explored the distribution of CGCs across the United States. We sought to identify the U.S. geographical regions with the highest number of CGCs and those regions where the physical presence of CGCs is sparser. Deidentified city, state, and ZIP code information for each CGC in the United States were obtained from the American Board of Genetic Counseling (ABGC) database. A countrywide analysis of the distribution of CGCs was completed using geographic information system (GIS) mapping software. The data were organized into U.S. metropolitan or micropolitan statistical areas, if applicable, and analyzed by CGC per capita. We included a total of 4,554 data points (92.2%) in the analysis. Results showed there is one CGC for every 71,842 people nationwide. Of 3,141 total counties (or county equivalents) in the United States, 535 counties had at least one CGC (17.0%). The majority (98.7%) of CGCs live or work within metropolitan statistical areas (MSAs), which are defined by this study as geographical areas with greater than 50,000 people. Of the MSAs with a CGC, approximately half have more than one CGC per 100,000 people. These results are consistent with the overall distribution of the U.S. population. We believe that the MSAs with the most CGCs per capita are due to associations with specific institutions, that is, genetic counseling training programs, health system headquarters, or genetic laboratories. Although the present study cannot draw definite conclusions regarding direct patient care services provided by CGCs, it does provide a snapshot of current CGC distribution within the country. Knowing the distribution of CGCs provides a tool to conduct further workforce analyses to determine the number of CGCs needed to serve the U.S. population.

Signaling at Physical Barriers during Pollen-Pistil Interactions.

In angiosperms, double fertilization requires pollen tubes to transport non-motile sperm to distant egg cells housed in a specialized female structure known as the pistil, mediating the ultimate fusion between male and female gametes. During this journey, the pollen tube encounters numerous physical barriers that must be mechanically circumvented, including the penetration of the stigmatic papillae, style, transmitting tract, and synergid cells as well as the ultimate fusion of sperm cells to the egg or central cell. Additionally, the pollen tube must maintain structural integrity in these compact environments, while responding to positional guidance cues that lead the pollen tube to its destination. Here, we discuss the nature of these physical barriers as well as efforts to genetically and cellularly identify the factors that allow pollen tubes to successfully, specifically, and quickly circumnavigate them.

Functional analysis tools for post-translational modification: a post-translational modification database for analysis of proteins and metabolic pathways.

Post-translational modifications (PTMs) are critical regulators of protein function, and nearly 200 different types of PTM have been identified. Advances in high-resolution mass spectrometry have led to the identification of an unprecedented number of PTM sites in numerous organisms, potentially facilitating a more complete understanding of how PTMs regulate cellular behavior. While databases have been created to house the resulting data, most of these resources focus on individual types of PTM, do not consider quantitative PTM analyses or do not provide tools for the visualization and analysis of PTM data. Here, we describe the Functional Analysis Tools for Post-Translational Modifications (FAT-PTM) database (https://bioinformatics.cse.unr.edu/fat-ptm/), which currently supports eight different types of PTM and over 49 000 PTM sites identified in large-scale proteomic surveys of the model organism Arabidopsis thaliana. The FAT-PTM database currently supports tools to visualize protein-centric PTM networks, quantitative phosphorylation site data from over 10 different quantitative phosphoproteomic studies, PTM information displayed in protein-centric metabolic pathways and groups of proteins that are co-modified by multiple PTMs. Overall, the FAT-PTM database provides users with a robust platform to share and visualize experimentally supported PTM data, develop hypotheses related to target proteins or identify emergent patterns in PTM data for signaling and metabolic pathways.

Convergent evolution of hetero-oligomeric cellulose synthesis complexes in mosses and seed plants.

In seed plants, cellulose is synthesized by rosette-shaped cellulose synthesis complexes (CSCs) that are obligate hetero-oligomeric, comprising three non-interchangeable cellulose synthase (CESA) isoforms. The moss Physcomitrella patens has rosette CSCs and seven CESAs, but its common ancestor with seed plants had rosette CSCs and a single CESA gene. Therefore, if P. patens CSCs are hetero-oligomeric, then CSCs of this type evolved convergently in mosses and seed plants. Previous gene knockout and promoter swap experiments showed that PpCESAs from class A (PpCESA3 and PpCESA8) and class B (PpCESA6 and PpCESA7) have non-redundant functions in secondary cell wall cellulose deposition in leaf midribs, whereas the two members of each class are redundant. Based on these observations, we proposed the hypothesis that the secondary class A and class B PpCESAs associate to form hetero-oligomeric CSCs. Here we show that transcription of secondary class A PpCESAs is reduced when secondary class B PpCESAs are knocked out and vice versa, as expected for genes encoding isoforms that occupy distinct positions within the same CSC. The class A and class B isoforms co-accumulate in developing gametophores and co-immunoprecipitate, suggesting that they interact to form a complex in planta. Finally, secondary PpCESAs interact with each other, whereas three of four fail to self-interact when expressed in two different heterologous systems. These results are consistent with the hypothesis that obligate hetero-oligomeric CSCs evolved independently in mosses and seed plants and we propose the constructive neutral evolution hypothesis as a plausible explanation for convergent evolution of hetero-oligomeric CSCs.

Associations between phytohormones and cellulose biosynthesis in land plants.

BACKGROUND: Phytohormones are small molecules that regulate virtually every aspect of plant growth and development, from basic cellular processes, such as cell expansion and division, to whole plant environmental responses. While the phytohormone levels and distribution thus tell the plant how to adjust itself, the corresponding growth alterations are actuated by cell wall modification/synthesis and internal turgor. Plant cell walls are complex polysaccharide-rich extracellular matrixes that surround all plant cells. Among the cell wall components, cellulose is typically the major polysaccharide, and is the load-bearing structure of the walls. Hence, the cell wall distribution of cellulose, which is synthesized by large Cellulose Synthase protein complexes at the cell surface, directs plant growth. SCOPE: Here, we review the relationships between key phytohormone classes and cellulose deposition in plant systems. We present the core signalling pathways associated with each phytohormone and discuss the current understanding of how these signalling pathways impact cellulose biosynthesis with a particular focus on transcriptional and post-translational regulation. Because cortical microtubules underlying the plasma membrane significantly impact the trajectories of Cellulose Synthase Complexes, we also discuss the current understanding of how phytohormone signalling impacts the cortical microtubule array. CONCLUSION: Given the importance of cellulose deposition and phytohormone signalling in plant growth and development, one would expect that there is substantial cross-talk between these processes; however, mechanisms for many of these relationships remain unclear and should be considered as the target of future studies.

Secular decline in limb bone strength among South African Africans during the 19th and 20th centuries.

OBJECTIVES: South African Africans have been reported to have experienced negative or null secular trends in stature and other measures of skeletal structure across the 19th and 20th centuries, presumably due to poor living conditions during a time of intensifying racial discrimination. Here, we investigate whether any secular trend is apparent in limb bone strength during the same period. MATERIALS AND METHODS: Cadaver-derived skeletons (n = 221) were analyzed from female and male South African Africans who were born between 1839 and 1970, lived in and around Johannesburg, and died between 1925 and 1991 when they were 17-90 years of age. For each skeleton, a humerus and femur were scanned using computed tomography, and mid-diaphyseal cross-sectional geometric properties were calculated and scaled according to body size. RESULTS: In general linear mixed models accounting for sex, age at death, and skeletal element, year of birth was a significant (p < .05) negative predictor of size-standardized mid-diaphyseal cortical area (a proxy for resistance to axial loading) and polar moment of area (a proxy for resistance to bending and torsion), indicating a temporal trend toward diminishing limb bone strength. No significant interactions were detected between year of birth and age at death, suggesting that the decline in limb bone strength was mainly due to changes in skeletal maturation rather than severity of age-related bone loss. DISCUSSION: Limb bone strength is thus potentially another feature of the skeletal biology of South African Africans that was compromised by poor living conditions during the 19th and 20th centuries.

BRASSINOSTEROID INSENSITIVE2 negatively regulates cellulose synthesis in Arabidopsis by phosphorylating cellulose synthase 1.

The deposition of cellulose is a defining aspect of plant growth and development, but regulation of this process is poorly understood. Here, we demonstrate that the protein kinase BRASSINOSTEROID INSENSITIVE2 (BIN2), a key negative regulator of brassinosteroid (BR) signaling, can phosphorylate Arabidopsis cellulose synthase A1 (CESA1), a subunit of the primary cell wall cellulose synthase complex, and thereby negatively regulate cellulose biosynthesis. Accordingly, point mutations of the BIN2-mediated CESA1 phosphorylation site abolished BIN2-dependent regulation of cellulose synthase activity. Hence, we have uncovered a mechanism for how BR signaling can modulate cellulose synthesis in plants.

Knee osteoarthritis has doubled in prevalence since the mid-20th century.

Knee osteoarthritis (OA) is believed to be highly prevalent today because of recent increases in life expectancy and body mass index (BMI), but this assumption has not been tested using long-term historical or evolutionary data. We analyzed long-term trends in knee OA prevalence in the United States using cadaver-derived skeletons of people aged ≥50 y whose BMI at death was documented and who lived during the early industrial era (1800s to early 1900s; n = 1,581) and the modern postindustrial era (late 1900s to early 2000s; n = 819). Knee OA among individuals estimated to be ≥50 y old was also assessed in archeologically derived skeletons of prehistoric hunter-gatherers and early farmers (6000-300 B.P.; n = 176). OA was diagnosed based on the presence of eburnation (polish from bone-on-bone contact). Overall, knee OA prevalence was found to be 16% among the postindustrial sample but only 6% and 8% among the early industrial and prehistoric samples, respectively. After controlling for age, BMI, and other variables, knee OA prevalence was 2.1-fold higher (95% confidence interval, 1.5-3.1) in the postindustrial sample than in the early industrial sample. Our results indicate that increases in longevity and BMI are insufficient to explain the approximate doubling of knee OA prevalence that has occurred in the United States since the mid-20th century. Knee OA is thus more preventable than is commonly assumed, but prevention will require research on additional independent risk factors that either arose or have become amplified in the postindustrial era.

Secular trends in cranial size and shape among black South Africans over the late 19th and 20th centuries.

BACKGROUND: Previous studies of secular change in cranial size among black South Africans have produced conflicting results. AIM: We re-examined cranial size change in this population during the 19th and 20th century by evaluating its relationship with individual year-of-birth, and the significance of trends among eight decennial cohorts. SUBJECTS AND METHODS: This study is based on 102 male and 89 female adults born between 1865 and 1959. Linear regressions were employed to evaluate possible relationships between year-of-birth and cranial dimensions; ANOVAs were used to evaluate the significance of long-term trends among decennial cohorts. RESULTS: No analysis revealed a secular change in cranial length in either sex; however, the ANOVA for cranial length in the combined sex sample was significant. There is no secular trend in female cranial breadth, but males display a negative trend in this dimension. This results in a secular trend for increased male dolichocephaly. CONCLUSIONS: The factors that underlie the negative secular trend in male cranial breadth and the absence of a secular trend in overall cranial size in this population are unclear. Nevertheless, these observations accord with findings related to stature and long bone strength in this population and are consistent with observations for other sub-Saharan African populations.

Interspecies Bombolitins Exhibit Structural Diversity upon Membrane Binding, Leading to Cell Specificity.

Bombolitins, a class of peptides produced by bees of the genus Bombus, target and disrupt cellular membranes, leading to lysis. Antimicrobial peptides exhibit various mechanisms of action resulting from the interplay between peptide structure, lipid composition, and cellular target membrane selectivity. Herein, two bombolitins displaying significant amino-acid-sequence similarity, BII and BL6, were assessed for antimicrobial activity as well as correlated dodecylphosphocholine (DPC) micelle binding and membrane-induced peptide conformational changes. Infrared and circular dichroism spectroscopies were used to assess the structure-function relationship of each bombolitin, and the results indicate that BII forms a rigid and helically ordered secondary structure upon binding to DPC micelles, whereas BL6 largely lacks secondary structural order. Moreover, the binding affinity of each peptide to DPC micelles was determined, revealing that BL6 displayed a difference in binding affinity by over two orders of magnitude. Further investigations into the growth-inhibitory activity of the two bombolitins were performed against Escherichia coli and Saccharomyces cerevisiae. Interestingly, BII specifically targeted S. cerevisiae, whereas BL6 more effectively inhibited E. coli growth. Overall, the antimicrobial selectivity and specificity of BII and BL6 are largely dependent on the primary as well as secondary structural content of the peptides and the membrane composition.

Knee osteoarthritis risk in non-industrial societies undergoing an energy balance transition: evidence from the indigenous Tarahumara of Mexico.

Non-industrial societies with low energy balance levels are expected to be less vulnerable than industrial societies to diseases associated with obesity including knee osteoarthritis. However, as non-industrial societies undergo rapid lifestyle changes that promote positive energy balance, individuals whose metabolisms are adapted to energetic scarcity are encountering greater energy abundance, increasing their propensity to accumulate abdominal adipose tissue and thus potentially their sensitivity to obesity-related diseases. OBJECTIVES: Here, we propose that knee osteoarthritis is one such disease for which susceptibility is amplified by this energy balance transition. METHODS: Support for our hypothesis comes from comparisons of knee radiographs, knee pain and anthropometry among men aged ≥40 years in two populations: Tarahumara subsistence farmers in Mexico undergoing the energy balance transition and urban Americans from Framingham, Massachusetts. RESULTS: We show that despite having markedly lower obesity levels than the Americans, the Tarahumara appear predisposed to accrue greater abdominal adiposity (ie, larger abdomens) for a given body weight, and are more vulnerable to radiographic and symptomatic knee osteoarthritis at lower levels of body mass index. Also, proportionate increases in abdomen size in the two groups are associated with greater increases in radiographic knee osteoarthritis risk among the Tarahumara than the Americans, implying that the abdominal adipose tissue of the Tarahumara is a more potent stimulus for knee degeneration. CONCLUSIONS: Heightened vulnerability to knee osteoarthritis among non-industrial societies experiencing rapid lifestyle changes is a concern that warrants further investigation since such groups represent a large but understudied fraction of the global population.

A Putative Protein O-Fucosyltransferase Facilitates Pollen Tube Penetration through the Stigma-Style Interface.

During pollen-pistil interactions in angiosperms, the male gametophyte (pollen) germinates to produce a pollen tube. To fertilize ovules located within the female pistil, the pollen tube must physically penetrate specialized tissues. Whereas the process of pollen tube penetration through the pistil has been anatomically well described, the genetic regulation remains poorly understood. In this study, we identify a novel Arabidopsis (Arabidopsis thaliana) gene, O-FUCOSYLTRANSFERASE1 (AtOFT1), which plays a key role in pollen tube penetration through the stigma-style interface. Semi-in vivo growth assays demonstrate that oft1 mutant pollen tubes have a reduced ability to penetrate the stigma-style interface, leading to a nearly 2,000-fold decrease in oft1 pollen transmission efficiency and a 5- to 10-fold decreased seed set. We also demonstrate that AtOFT1 is localized to the Golgi apparatus, indicating its potential role in cellular glycosylation events. Finally, we demonstrate that AtOFT1 and other similar Arabidopsis genes represent a novel clade of sequences related to metazoan protein O-fucosyltransferases and that mutation of residues that are important for O-fucosyltransferase activity compromises AtOFT1 function in vivo. The results of this study elucidate a physiological function for AtOFT1 in pollen tube penetration through the stigma-style interface and highlight the potential importance of protein O-glycosylation events in pollen-pistil interactions.

Experimental evidence that physical activity affects the multivariate associations among muscle attachments (entheses).

The morphology of entheses (muscle/tendon attachment sites) on bones is routinely used in paleontological and bioarcheological studies to infer the physical activity patterns of ancient vertebrate species including hominins. However, such inferences have often been disputed owing to limitations of the quantitative methods commonly employed and a lack of experimental evidence demonstrating direct effects of physical activity on entheseal morphology. Recently, we introduced a new and improved method of quantifying and analyzing entheseal morphology that involves repeatable three-dimensional measurements combined with multivariate statistics focused on associations among multiple entheses. Here, to assess the validity of our method for investigating variation in entheseal morphology related to physical activity patterns, we analyzed femora of growing turkeys that were experimentally exercised for 10 weeks on either an inclined or declined treadmill or served as controls (N=15 individuals, 5 per group). Our multivariate approach identified certain patterns involving three different entheses (associated with the gluteus primus, medial gastrocnemius, vastus medialis and adductor magnus muscles) that clearly differentiated controls from runners. Importantly, these differences were not observable when comparing groups within each of the three entheseal structures separately. Body mass was not correlated with the resulting multivariate patterns. These results provide the first experimental evidence that variation in physical activity patterns has a direct influence on entheseal morphology. Moreover, our findings highlight the promise of our newly developed quantitative methods for analyzing the morphology of entheses to reconstruct the behavior of extinct vertebrate species based on their skeletal remains.