Next Article in Journal
Sleep Difficulties in Swiss Elite Athletes
Next Article in Special Issue
The Morphology of the Femur Influences the Fracture Risk during Stumbling and Falls on the Hip—A Computational Biomechanical Study
Previous Article in Journal
Screening for Breast Cancer: A Comparative Review of Guidelines
Previous Article in Special Issue
What Is New about the Semimembranosus Distal Tendon? Ultrasound, Anatomical, and Histological Study with Clinical and Therapeutic Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of 3D Footprint Morphology of Knee-Related Muscle Attachments Based on CT Data Reconstruction: A Feasibility Study

1
Institute for Anatomy, Rostock University Medical Center, Gertrudenstraße 9, 18057 Rostock, Germany
2
Institute of Molecular and Cellular Anatomy, University of Ulm, Albert-Einstein-Allee 11, 89081 Ulm, Germany
3
Department of Orthopaedics, Rostock University Medical Center, Doberaner Straße 142, 18055 Rostock, Germany
4
Department of Radiotherapy and Radiation Oncology, Rostock University Medical Center, Südring 75, 18059 Rostock, Germany
*
Author to whom correspondence should be addressed.
Life 2024, 14(6), 778; https://doi.org/10.3390/life14060778
Submission received: 23 May 2024 / Revised: 13 June 2024 / Accepted: 14 June 2024 / Published: 19 June 2024
(This article belongs to the Special Issue Topographic and Functional Anatomy of Musculoskeletal System)

Abstract

:
A three-dimensional (3D) understanding of muscle attachment footprints became increasingly relevant for musculoskeletal modeling. The established method to project attachments as points ignores patient-specific individuality. Research focuses on investigating certain muscle groups rather than comprehensively studying all muscles spanning a joint. Therefore, we present a reliable method to study several muscle attachments in order to reconstruct the attachment sites in 3D based on CT imaging for future applications in musculoskeletal modeling. For the present feasibility study, 23 knee-related muscle attachments were CT-scanned postmortem from four nonadipose male specimens. For this, the specific muscle attachments were dissected and marked with a barium sulfate containing paint (60 g BaSO4 in 30 mL water and 10 mL acrylic paint). Subsequently, bone geometries and muscle attachments were reconstructed and evaluated from CT datasets. Bone morphology and footprint variations were studied. Exemplarily, variations were high for pes anserinus insertions (mean 56%) and the origins of M. biceps femoris (mean 54%). In contrast, the origins of the vastus muscles as well as the insertion of the Achilles tendon showed low variation (mean 9% and 13%, respectively). Most attachment sites showed variation exceeding the individuality of bone morphology. In summary, the present data were consistent with the few published studies of specific muscle footprints. Our data shed light on the high variability of muscle attachments, which need to be addressed when studying muscle forces and movements through musculoskeletal modeling. This is the first step to achieving a more profound understanding of muscle morphology to be utilized in numerical simulations.

1. Introduction

The patellofemoral joint is known to bear some of the highest loads in the human body. It is estimated to carry up to 2.5 to 2.8× the body weight during easy walking compared to, e.g., 1.2× in the foot [1,2]. Furthermore, the patellofemoral load varies considerably during different activities, ranging from 0.6× body weight during walking to more than 8× body weight during a single-leg squat [3]. As such, diffuse knee pain syndromes and knee-related sports injuries are common. Multicausal anterior knee pain (AKP) is associated with muscle or tendon dysfunctions [4]. It is especially found in young athletes. Nearly 40% of adolescent athletes independent of primary sport played, age or sex reported AKP upon provocation [5]. AKP accounts for 20% of female sports injuries compared to 7.5% of male injuries [6]. Despite this, treatment and etiology of this disease are still insufficiently understood. Numerical simulations can help in predicting potential sources for dysfunctions but for reliable subject-specific models and meaningful conclusions accurate musculoskeletal geometry is essential [7,8,9]. Furthermore, musculoskeletal models are frequently used to support preoperative planning to evaluate joint dynamics and to simulate physiological or pathological movements [8,9,10]. Biomechanical engineers and clinicians aim to increase their efforts in building a diverse database of muscle attachment details in relation to various ages, body types and lifestyles of patients.
Qualitative or descriptive studies of specific muscle attachments and ligaments are invaluable and have progressed understanding of individuality in muscle anatomy [11,12,13]. However, attachment sites on the bone surfaces are often studied exclusively two-dimensionally based on macroscopic images, plain radiographs or magnetic resonance imaging (MRI) [14]. Three-dimensional (3D) characterization of footprints allows for more thorough morphological analysis [15,16]. The study of bone curvature and shape derived from 3D image data allows for reliable estimations of some attachment surfaces on bones, but it lacks single muscle resolution on conjoined attachment sites [17,18]. Furthermore, in statistical shape modeling, the error increases as the attachment size is reduced [19]. Additionally, Carbone et al. [7] introduced a valuable database (TLEM 2.0) which contains data on muscle attachment sites for musculoskeletal modeling based on the dissection of one human specimen. In this regard, Andreassen et al. [20] presented a male and female dataset based on medical images that included muscle geometries and attachment areas. However, these datasets provide no individual differences in the attachment areas.
It is also possible to deploy probabilistic approaches. For this, Fukuda et al. [21] proposed a probabilistic model, where they investigated the differences in muscle attachment areas for the hip muscles based on the dissection of eight human specimens. The individual muscle attachment areas of the hip region were defined using an optical tracker. Based on a probabilistic model with these measurements, specific muscle attachment areas for simulation models can be estimated and mapped using CT-derived models. Although this is a promising approach, there were also outliers in the measurement that had to be removed manually. Herteleer et al. [22] showed the variation of muscle attachment areas in the clavicle. However, they did not analyze the knee-relevant muscle areas or the centroids of the attachment areas. By using a 3D-digitizer and microscribes, thorough information of upper extremity attachments was acquired [23,24]. The hamstring origins were also studied in detail this way [25]. However, as such a set-up is expensive and highly sensitive to calibration, barium sulfate can be injected to visualize ligaments in computed tomography (CT) images [26,27]. As a radiopaque paint, it offers an inexpensive alternative to mark attachments and reconstruct them during segmentation [28]. This allows for high accuracy by avoiding the distance between points of measurement. Moreover, the gold standard in musculoskeletal modeling is to offer corresponding points of functional muscle subunits [7,29] or to use generalized models [30]. This approach ignores the anatomical variation in location or morphology of individual attachment areas [21] and leads to inaccurate predictions from imprecise modeling [8] as it can affect the muscle moment arms [31,32], muscle force calculation, as well as joint dynamics [8,9,33]. Davico et al. [9] investigated the personalization of musculoskeletal models and showed that the musculoskeletal anatomy and muscle activation patterns, particularly, had a considerable influence on joint force calculation.
Muscle attachments related to the patellofemoral and tibiofemoral sites are rarely studied extensively. Instead, investigators often focused on single muscles, evaluated muscle attachments in a simplified manner or did not focus on the individual differences in attachment areas [7,14,25]. In the present feasibility study, we introduced a reliable method to study several muscle attachments in order to reconstruct the attachment sites in 3D based on CT imaging for future applications in musculoskeletal modeling. We aimed to quantify the individual differences in muscle attachment areas for knee-related muscles.

2. Materials and Methods

2.1. Specimen Preparation

Four nonpaired, formalin-fixed human legs were used to prepare knee-related muscle attachment areas in four steps. This was approved by the local ethical committee of the University of Rostock (registration number A 2016-00083). Specimen preparation was performed by the authors assisted by medical prosectors. Attachment areas were freed of surrounding tissues by scissors and pulled off the bones. The outlines of the bony attachments were marked with a surgical marker and painted with a barium sulfate containing paint (60 g BaSO4 in 30 mL water and 10 mL acrylic paint; Figure 1A,B).
The muscles origins (O) and insertions (I) were dissected step-wise (plantaris and popliteus muscles were not evaluated) and a CT scan (Brilliance CT Big Bore, Philips, Amsterdam, The Netherlands) was performed after each dissection with a slice thickness of 1 mm (varying between 716 to 1063 slices per scan) covering the following structures: CT1: M. semimembranosus (O&I), M. sartorius (O&I), M. gastrocnemius (O), M. rectus femoris (O); CT2: M. biceps femoris (O&I), M. gracilis (O&I), M. vastus medialis (O); CT3: M. semitendinosus (O&I), M. vastus lateralis (O); CT4: Patellar tendon (O&I), M. vastus intermedius (O), M. quadriceps tendon (I) and Achilles tendon (I).
In the case of multipart insertions (e.g., M. semimembranosus), attachments into ligaments or joint capsules were excluded. Possibly conjoined tendons were carefully separated along the fibers until the surface of the bone was reached (e.g., hamstring origins); the quadriceps tendon, however, was taken as one. The patella is tediously surrounded, nevertheless, quadriceps insertion and patellar tendon origin were separated based on the direction of the fibers when attached to the bony surface. A small area of less than 5 mm in width along the patella was left unpainted in between for distinction during segmentation.

2.2. Bone Reconstruction

The bone geometry and the adjacent highly radiopaque paint were reconstructed from CT datasets in AMIRA® v.5.4.1 (v5.4.1, Zuse Institute Berlin, Berlin, Germany; Thermo Fisher Scientific, Waltham, MA, USA). The structures were segmented layer-by-layer using the image segmentation editor with an intensity threshold based on Hounsfield units (HU) (Figure 1C–E). As established, we used a primary HU value for bone of 250 up to 3000 [34]. Strongly calcified tendons were manually excluded (except in case of the quadriceps insertion). The surfaces of identical femur, pelvis, tibia and fibula were generated repeatedly, whenever an attachment was marked on it, to decrease inaccuracies of the labeling at the paint/bone contact area. The patella was labelled twice, with and without a mark. The foot was labelled only with a mark after CT4 since the maximal scanning length craniocaudally did not allow for a simultaneous scan from pelvis to foot. When segmentation was completed, 3D surfaces of the bones and the paint were reconstructed automatically using triangulated surfaces.

2.3. Generation of Coherent Models

These surface models were imported into GEOMAGIC studio v.13 (v2013, 3D Systems, Rock Hill, SC, USA) via a STL interface in ASCII mode. All bone surfaces were transformed based on the donor’s anatomy during CT1 and a 3D comparison was performed to evaluate the segmentation process. The different surfaces of the identical bones were merged, corrected (removal of holes and sharp edges) and smoothed, while the surfaces of the unmarked patella and the foot were only corrected and smoothed. Hence, one coherent 3D surface model for each leg was constructed (Figure 1F). Volume’s center of gravity, volume and surface area (SA) for every bone were calculated. Additionally, femoral head diameter (FHD), shaft diameter (SD), length of the mechanical axis (MA) and length of the transepicondylar axis (TEA) of the femur were measured. All data are expressed in the defined coordinate frame of the femur of each leg, with the center of the TEA as the origin, using a standardized coordinate system according to the International Society of Biomechanics definition [35]. Briefly, the y-axis is defined by the TEA midpoint and femoral head center point towards cranially, the z-axis lies perpendicular to the y-axis in the plane with the origin and femoral head center point, and the x-axis lies perpendicular to both other axes pointing anteriorly.

2.4. Projection of Attachments on Bone Surfaces

The surfaces of the paint marks were aligned using the transformation matrices of the bones calculated previously for the construction of coherent models. In case of overlap between the attachment paint marks (e.g., hamstring origins, origins around the linea aspera, pes anserinus superficialis insertions), the later applied mark was slightly corrected. The contact area of the paint with the surface of the bones was manually traced on the bone in GEOMAGIC studio and generated as a separate surface (Figure 1G–I).

2.5. Muscle Data Collection

The footprint and centroid of each attachment area were calculated in GEOMAGIC studio. For use in future musculoskeletal modeling, a breakthrough point of the centroid was generated as the intersection point with the attachment surface of a line perpendicular to a best-fit plane of the attachment area through the centroid. In this manner, we created attachment points on the surface of the respective bones. Some muscle attachments needed further division for meaningful point generation. Linear attachments were divided three times in length from the most proximal to the most distal point. The vastus intermedius was separated six times: in half by a plane passing through the most proximal point, the most distal point, and the center point between the lateral and medial edges, and three times equally in length.Vastus medialis and vastus lateralis were separated in a superior part (polygonal shaped) and an inferior part (linear shaped, separated three times) by the anterior and distal margin of the trochanter minor, respectively. The fibular insertion of biceps femoris was separated by a plane through the coordinate origin, the centroid and the apex of the distal crest of the fibular head. If a single breakthrough point was still not determinable for the attachment unit, it was separated along the thinnest part. The tibial insertion of the biceps femoris of Specimen 3 was naturally separated into two areas, which were combined for evaluation.

2.6. Data Analysis

Results are expressed as mean with the standard deviation (STD) of the population and the corresponding coefficient of variation (CV). The deviations from the reference model of CT1 were calculated using GEOMAGIC studio. Furthermore, volume’s center of gravity, volume and surface area (SA) for every bone were calculated.

3. Results

3.1. Specimen Characterization

The four legs (three left, one right) were from nonadipose male donors with an age between 68 and 78 years and a height of 170 to 180 cm. Causes of death were acute myocardial infarct, heart failure and multiorgan failure after metastatic cancer. Previous diseases were of the spectrum of cardiovascular diseases, chronic obstructive pulmonary disease and diabetes mellitus. Musculoskeletal diseases were not known or detected; however, some tendon calcifications were observed (Supplementary Figure S1). Muscle or bone quality were not evaluated.

3.2. Segmentation Deviation Analysis

We calculated the deviations (mm) of each subsequent CT reconstruction from the reference surface of CT1 (Table 1). All segmented femoral surfaces had a mean deviation from the reference surface of 0.48 ± 0.14 mm, the pelvic surfaces of 0.40 ± 0.10 mm, the tibiofibular surfaces of 0.54 ± 0.15 mm and the patellar surfaces of 0.70 ± 0.15 mm; altogether, there was a deviation from the reference surfaces of 0.50 ± 0.16 mm. The highest deviation after segmentation was seen underneath the painted muscle attachment areas (Figure 2A).

3.3. Morphological Data

Volumes, SA and volume’s centers of gravity of the bones were calculated (Table 2). Volumes and SA had a mean CV of 15.94 ± 1.69% and 11.42 ± 1.55%, respectively. We measured femoral morphology (Table 3), evaluated muscle attachment areas and centroids (Table 4) and generated breakthrough points for each attachment site (Supplementary Table S1). Femoral morphological parameters varied around 7.38 ± 2.01% between individuals. On the contrary, the mean surface areas of the attachments showed large variations (Figure 2B); the mean CV for all attachment site areas was 33.12 ± 20.82%. Attachment areas with high variation (CV > 50%) were rectus femoris caput reflexum O (58.94%), semitendinosus I (53.03%) and sartorius I (80.05%). The insertions forming the pes anserinus superficialis (semitendinosus, sartorius, gracilis) varied from moderately to strongly (53.03%, 80.05%, 33.90%, respectively), which might indicate some kind of dependence from one another. The attachment areas with relatively low variation (CV < 15%) are the origins of semimembranosus (10.32%), vastus lateralis (3.29%) and vastus intermedius (7.74%) as well as the Achilles tendon I (12.94%).

4. Discussion

The focus of our present study was to establish an accessible and inexpensive method to evaluate muscle attachment sites to construct a database of individual attachment morphologies. Furthermore, the aim was to quantify individual differences of attachment areas of knee-related muscles. In our feasibility study, we showed a high diversity of muscle attachment sites in size and location, mostly exceeding the individual variations in bone morphology. We analyzed the segmentation deviation and found a low deviation of 0.50 ± 0.16 mm. The bones around the dense cortex sometimes shared HU with the paint. Labeling within the contact area in these cases is difficult and prone to minor errors but does not alter the size and location of attachments. In a previous study [34], the 3D segmentation process was verified in an inter-laboratory study. Using different software packages and algorithms, the participating groups independently extracted the 3D geometry of a single human femur from CT data. Four different segmentation software packages were used: AMIRA® (FEI Visualization Sciences Group, Hillsboro, OR, USA), Mimics® (Materialise N.V., Leuven, Belgium), YaDiv (Welfenlab, Leibniz Universität Hannover, Hannover, Germany), and Fiji Life-Line. Stereolithography files were imported into GEOMAGIC studio v.2013 (Raindrop Geomagic, Triangle Park, NC, USA). They found no crucial differences compared with an optical scan of the original surface [34].
We investigated the muscle attachment sites of formalin-fixed human leg specimens. Herein, we analyzed the locations and areal shapes of the lower extremity muscles spanning the knee joint. We identified the variance of locations and shapes of the attachments of the different specimens. Therefore, the present study provides important insights into individual differences in musculoskeletal geometry based on methodology to identify the muscle attachment areas using a CT scan. A geometry dataset is provided for musculoskeletal modeling, similar to [7], but with different specimens and varying muscle attachment areas. Based on medical imaging data, we considered a more inter-individual anatomical variability in size and shape of the attachment footprints. The implementation of these attachment sites into a musculoskeletal multibody model requires the definition of an origin and insertion point. For this purpose, we calculated breakthrough points.
In addition, the influence of the varying muscle attachment areas on the predicted muscle forces can be investigated in future studies by means of inverse dynamics as these parameters have a considerable effect in model predictions [8,36]. To our knowledge, there are rarely sensitivity analyses that address this issue since there is no reliable numeric information for the areal shape and size of each of the individual muscle subunits. Fukuda et al. [21] investigated the differences in muscle attachment site for the hip muscles using an optical tracker. Contrarily, our feasibility study investigated the variation in the muscle attachment area of knee-related muscles in different specimens, providing the possibility to perform sensitivity studies. Future studies could investigate the impact on musculoskeletal model prediction, e.g., quadriceps force or tibiofemoral contact force, by perturbation of the origin and insertion points within the respective muscle attachment area [36].
Nevertheless, there are some limitations in the present study. We did not include muscle volumes and subcutaneous fat. Furthermore, we used a limited sample size, which did not cover all anatomical variations of human legs. In this context, female geometries should especially be integrated into future studies as presented by Andreassen et al. [20]. Moreover, muscle architecture can also be changed besides the musculoskeletal geometry [37]. The technique of marking the muscle attachment sites with barium paint and the use of a CT scanner are time-consuming as each muscle was outlined manually with a surgical marker and each leg was scanned and segmented four times. Nevertheless, given the inexpensive nature of the radiopaque paint, this method is a very cost efficient and easy-to-implement technique.

Variations of Muscle Attachments

Literature about attachment dimensions is sparsely found. Some studies described specific anatomical components but a coherent analysis of the diverse muscular aspects affecting knee function is missing. The variation of footprints is of great relevance for the musculoskeletal geometry, especially muscle moment arm as one of the most sensitive parameters [7]. Further studies are needed to investigate the practical consequences of our presented results in musculoskeletal modeling. Particularly, the errors we might expect when we incorporate our data into a musculoskeletal model of the lower extremity. Depending on the attachment points, muscle moment arms and muscle–tendon lines of action differ, thus changing the predicted muscle forces. Consequently, differences in moment arms lead to inaccuracies in the prediction of muscle and joint contact forces [8,9,33]. Such joint contact forces contribute to the progression of joint arthritis [38]. It was also shown that considering subject-specific bone structures and attachment areas leads to more accurate joint dynamics predictions [9].
Despite variations in methodology, studies investigating specific muscle attachments are largely consistent with our results [28,39,40]. The hamstring origin complex is a particularly well-studied attachment site due to the high prevalence of hamstring injuries in sports [28]. The origins of biceps femoris caput longum and semitendinosus usually attach as a strong conjoined tendon, therefore being particularly sensitive to preparation methods. Compared to the study by Philippon et al. (2015), our mean overall hamstring origin was larger with 1336.52 mm2 vs. 991.70 mm2, due to a greater conjoined tendon footprint of 864.74 mm2 vs. 567.00 mm2, while the origins of the semimembranosus were alike (471.78 mm2 vs. 412.40 mm2 by Philippon) [25].
The semimembranosus insertional attachment is less investigated despite the relevance of the tendon for the posteromedial corner of the knee and related injuries. The insertion tendon is a complex structure with slips inserted into the medial collateral ligament, the oblique popliteal ligament and a groove on the posteromedial aspect of the tibia named pes anserinus profundus. The arrangement of these insertion slips was previously analyzed [12,14]. Fukuda et al. [21] quantified the differences in the muscle attachment areas for hip-related muscles for eight human specimens using an optical tracker but did not provide insertion points using this dataset. To our knowledge, no study quantified the insertional footprint of the knee-related muscle considering its polygonal shape. A similar complex structure is the pes anserinus superficialis including the insertions of semitendinosus, sartorius and gracilis. The distal tendons of these muscles are regularly used as autografts for ligament reconstructions [41], but morphological quantification data of the attachment area is largely missing. In our study, the sartorius insertion showed high variation, something allured to be due to a highly varying inferior part [13,42].
The quadriceps group, the main antagonist of the hamstring group and extensor of the knee, and the gastrocnemius with the Achilles tendon are related to knee dysfunctions and common sports injuries [4,43], but attachment morphologies are unstudied. Ryan et al. (2014) measured the widths and lengths of the rectus femoris heads attachments [44]. Given the substantial variations seen in these sites, comparing the study (n = 6) with ours (n = 4) is difficult. The Achilles tendon was studied in descriptive ways [11,40] and reported with a similar mean area and variation to ours (18% vs. 13%) [40].
In summary, our results are largely consistent with the few previously conducted studies. Since the number of studies is small and the methods differ, more data is needed to assess physiological footprint dimensions sufficiently. However, we present the first step in establishing a thorough database of individual attachment morphologies. Our feasibility study had some limitations: the method was recently established and may be more prone to error than 3D-digitizer and the number of specimens evaluated in this study was small. For a more thorough analysis, however, a larger cohort is needed.

5. Conclusions

Depending on morphology and location, we found large variations regarding attachment dimensions within our limited number of human specimens. We assume that lifestyle or genetics notably influence the shape and size of muscle attachment sites. Future studies should, therefore, include environmental factors and patient’s lifestyles and gender in their investigations of muscle attachments and we would recommend increasing the impact of individuality in biomedical engineering research. Nevertheless, to have representative results for the young population that is more concerned with related pathologies, studies with similar elaborateness but nonterminal methodology need to be conducted. New, more advanced imaging methods like spectral CT that are able to characterize musculoskeletal pathologies noninvasively may be able to evaluate muscle attachment areas in this population in the future [45]. Since we processed the 3D data for musculoskeletal modeling, we want to progress the application of the muscle attachment data in respective musculoskeletal simulations. This study provides a technique for an anatomical description of the muscle attachment sites for knee-related muscles that can be used in future numerical studies. This may help to determine variability of the anatomical positions of muscle attachment sites, which is useful for orthopedic surgeons and computational biomechanics.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life14060778/s1, Figure S1: Pathologies of specimen. [A] Calcification of the hamstring tendon (and sacrotuberous ligament, excluded from segmentation) and [A’] on the patella. [B] Patella dysplasia with medial trochlear hypoplasia. [C] Broken acetabulum rim, possibly postmortem; Table S1: Generated breakthrough points for muscle attachments and subunits (O = origin, I = insertion).

Author Contributions

Conceptualization, R.B. and A.W.; methodology, R.B., A.W., M.K., G.H. and A.-M.N.; validation, A.-M.N., M.K., R.B. and A.W.; formal analysis, A.-M.N. and M.K.; investigation, A.-M.N. and M.K.; resources, R.B.; data curation, A.-M.N., M.K., G.H., R.B. and A.W.; writing—original draft preparation, A.-M.N., M.K., R.B. and A.W.; writing—review and editing, A.-M.N., M.K., R.B., G.H. and A.W.; visualization, A.-M.N.; supervision, R.B. and A.W.; project administration, R.B. and A.W.; funding acquisition, R.B. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by an internal research program of the Rostock University Medical Center.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the Medical Faculty, University of Rostock (registration number A 2016-00083).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets used during the current study are available from the first author on reasonable request.

Acknowledgments

We gratefully acknowledge Marko Schulze and Laura Hiepe for their excellent support in preparation of the specimen as well as the technical staff of the CT unit in the Department of Radiotherapy and Radiation Oncology for their eager participation.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. D’Lima, D.D.; Fregly, B.J.; Patil, S.; Steklov, N.; Colwell, C.W. Knee Joint Forces: Prediction, Measurement, and Significance. Proc. Inst. Mech. Eng. H 2012, 226, 95–102. [Google Scholar] [CrossRef]
  2. Chan, C.W.; Rudins, A. Foot Biomechanics during Walking and Running. Mayo Clin. Proc. 1994, 69, 448–461. [Google Scholar] [CrossRef] [PubMed]
  3. Song, K.; Scattone Silva, R.; Hullfish, T.J.; Silbernagel, K.G.; Baxter, J.R. Patellofemoral Joint Loading Progression Across 35 Weightbearing Rehabilitation Exercises and Activities of Daily Living. Am. J. Sports Med. 2023, 51, 2110–2119. [Google Scholar] [CrossRef]
  4. Hiemstra, L.A.; Kerslake, S.; Irving, C. Anterior Knee Pain in the Athlete. Clin. Sports Med. 2014, 33, 437–459. [Google Scholar] [CrossRef] [PubMed]
  5. Harris, M.; Edwards, S.; Rio, E.; Cook, J.; Cencini, S.; Hannington, M.C.; Bonello, C.; Docking, S. Nearly 40% of Adolescent Athletes Report Anterior Knee Pain Regardless of Maturation Status, Age, Sex or Sport Played. Phys. Ther. Sport 2021, 51, 29–35. [Google Scholar] [CrossRef] [PubMed]
  6. DeHaven, K.E.; Lintner, D.M. Athletic Injuries: Comparison by Age, Sport, and Gender. Am. J. Sports Med. 1986, 14, 218–224. [Google Scholar] [CrossRef] [PubMed]
  7. Carbone, V.; Fluit, R.; Pellikaan, P.; van der Krogt, M.M.; Janssen, D.; Damsgaard, M.; Vigneron, L.; Feilkas, T.; Koopman, H.F.J.M.; Verdonschot, N. TLEM 2.0—A Comprehensive Musculoskeletal Geometry Dataset for Subject-Specific Modeling of Lower Extremity. J. Biomech. 2015, 48, 734–741. [Google Scholar] [CrossRef]
  8. Gerus, P.; Sartori, M.; Besier, T.F.; Fregly, B.J.; Delp, S.L.; Banks, S.A.; Pandy, M.G.; D’Lima, D.D.; Lloyd, D.G. Subject-Specific Knee Joint Geometry Improves Predictions of Medial Tibiofemoral Contact Forces. J. Biomech. 2013, 46, 2778–2786. [Google Scholar] [CrossRef]
  9. Davico, G.; Lloyd, D.G.; Carty, C.P.; Killen, B.A.; Devaprakash, D.; Pizzolato, C. Multi-Level Personalization of Neuromusculoskeletal Models to Estimate Physiologically Plausible Knee Joint Contact Forces in Children. Biomech. Model. Mechanobiol. 2022, 21, 1873–1886. [Google Scholar] [CrossRef]
  10. Putame, G.; Terzini, M.; Rivera, F.; Kebbach, M.; Bader, R.; Bignardi, C. Kinematics and Kinetics Comparison of Ultra-Congruent versus Medial-Pivot Designs for Total Knee Arthroplasty by Multibody Analysis. Sci. Rep. 2022, 12, 3052. [Google Scholar] [CrossRef]
  11. Ballal, M.S.; Walker, C.R.; Molloy, A.P. The Anatomical Footprint of the Achilles Tendon: A Cadaveric Study. Bone Jt. J. 2014, 96-B, 1344–1348. [Google Scholar] [CrossRef]
  12. Benninger, B.; Delamarter, T. Distal Semimembranosus Muscle-Tendon-Unit Review: Morphology, Accurate Terminology, and Clinical Relevance. Folia Morphol. 2013, 72, 1–9. [Google Scholar] [CrossRef]
  13. Dziedzic, D.; Bogacka, U.; Ciszek, B. Anatomy of Sartorius Muscle. Folia Morphol. 2014, 73, 359–362. [Google Scholar] [CrossRef]
  14. De Maeseneer, M.; Shahabpour, M.; Lenchik, L.; Milants, A.; De Ridder, F.; De Mey, J.; Cattrysse, E. Distal Insertions of the Semimembranosus Tendon: MR Imaging with Anatomic Correlation. Skelet. Radiol. 2014, 43, 781–791. [Google Scholar] [CrossRef]
  15. Lin, J.; Zhang, S.; Xin, E.; Liang, M.; Yang, L.; Chen, J. Anterior Cruciate Ligament Femoral Footprint Is Oblong-Ovate, Triangular, or Two-Tears Shaped in Healthy Young Adults: Three-Dimensional MRI Analysis. Knee Surg. Sports Traumatol. Arthrosc. 2023, 31, 5514–5523. [Google Scholar] [CrossRef]
  16. Westermann, R.W.; Sybrowsky, C.; Ramme, A.J.; Amendola, A.; Wolf, B.R. Three-Dimensional Characterization of the Femoral Footprint of the Posterior Cruciate Ligament. Arthroscopy 2013, 29, 1811–1816. [Google Scholar] [CrossRef]
  17. Karakostis, F.A. Statistical Protocol for Analyzing 3D Muscle Attachment Sites Based on the “Validated Entheses-based Reconstruction of Activity” (VERA) Approach. Int. J. Osteoarchaeol. 2023, 33, 461–474. [Google Scholar] [CrossRef]
  18. Karakostis, F.A.; Lorenzo, C. Morphometric Patterns among the 3D Surface Areas of Human Hand Entheses. Am. J. Phys. Anthropol. 2016, 160, 694–707. [Google Scholar] [CrossRef]
  19. Salhi, A.; Burdin, V.; Mutsvangwa, T.; Sivarasu, S.; Brochard, S.; Borotikar, B. Subject-Specific Shoulder Muscle Attachment Region Prediction Using Statistical Shape Models: A Validity Study. In Proceedings of the 2017 39th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Jeju, Republic of Korea, 11–15 July 2017; pp. 1640–1643. [Google Scholar]
  20. Andreassen, T.E.; Hume, D.R.; Hamilton, L.D.; Walker, K.E.; Higinbotham, S.E.; Shelburne, K.B. Three Dimensional Lower Extremity Musculoskeletal Geometry of the Visible Human Female and Male. Sci. Data 2023, 10, 34. [Google Scholar] [CrossRef]
  21. Fukuda, N.; Otake, Y.; Takao, M.; Yokota, F.; Ogawa, T.; Uemura, K.; Nakaya, R.; Tamura, K.; Grupp, R.B.; Farvardin, A.; et al. Estimation of Attachment Regions of Hip Muscles in CT Image Using Muscle Attachment Probabilistic Atlas Constructed from Measurements in Eight Cadavers. Int. J. CARS 2017, 12, 733–742. [Google Scholar] [CrossRef]
  22. Herteleer, M.; Vancleef, S.; Herijgers, P.; Duflou, J.; Jonkers, I.; Vander Sloten, J.; Nijs, S. Variation of the Clavicle’s Muscle Insertion Footprints—A Cadaveric Study. Sci. Rep. 2019, 9, 16293. [Google Scholar] [CrossRef]
  23. Kamineni, S.; Bachoura, A.; Behrens, W.; Kamineni, E.; Deane, A. Distal Insertional Footprint of the Brachialis Muscle: 3D Morphometric Study. Anat. Res. Int. 2015, 2015, 786508. [Google Scholar] [CrossRef]
  24. Walton, C.; Li, Z.; Pennings, A.; Agur, A.; Elmaraghy, A. A 3-Dimensional Anatomic Study of the Distal Biceps Tendon: Implications for Surgical Repair and Reconstruction. Orthop. J. Sports Med. 2015, 3, 2325967115585113. [Google Scholar] [CrossRef]
  25. Philippon, M.J.; Ferro, F.P.; Campbell, K.J.; Michalski, M.P.; Goldsmith, M.T.; Devitt, B.M.; Wijdicks, C.A.; LaPrade, R.F. A Qualitative and Quantitative Analysis of the Attachment Sites of the Proximal Hamstrings. Knee Surg. Sports Traumatol. Arthrosc. 2015, 23, 2554–2561. [Google Scholar] [CrossRef]
  26. Bonnin, M.P.; de Kok, A.; Verstraete, M.; Van Hoof, T.; Van Der Straten, C.; Saffarini, M.; Victor, J. Popliteus Impingement after TKA May Occur with Well-Sized Prostheses. Knee Surg. Sports Traumatol. Arthrosc. 2017, 25, 1720–1730. [Google Scholar] [CrossRef]
  27. Van Hoof, T.; Cromheecke, M.; Tampere, T.; D’herde, K.; Victor, J.; Verdonk, P.C.M. The Posterior Cruciate Ligament: A Study on Its Bony and Soft Tissue Anatomy Using Novel 3D CT Technology. Knee Surg. Sports Traumatol. Arthrosc. 2013, 21, 1005–1010. [Google Scholar] [CrossRef]
  28. Neuschwander, T.B.; Benke, M.T.; Gerhardt, M.B. Anatomic Description of the Origin of the Proximal Hamstring. Arthroscopy 2015, 31, 1518–1521. [Google Scholar] [CrossRef]
  29. Klein Horsman, M.D.; Koopman, H.F.J.M.; van der Helm, F.C.T.; Prosé, L.P.; Veeger, H.E.J. Morphological Muscle and Joint Parameters for Musculoskeletal Modelling of the Lower Extremity. Clin. Biomech. 2007, 22, 239–247. [Google Scholar] [CrossRef]
  30. Arnold, E.M.; Ward, S.R.; Lieber, R.L.; Delp, S.L. A Model of the Lower Limb for Analysis of Human Movement. Ann. Biomed. Eng. 2010, 38, 269–279. [Google Scholar] [CrossRef]
  31. Kaya Keles, C.S.; Ates, F. How Mechanics of Individual Muscle-Tendon Units Define Knee and Ankle Joint Function in Health and Cerebral Palsy—A Narrative Review. Front. Bioeng. Biotechnol. 2023, 11, 1287385. [Google Scholar] [CrossRef]
  32. Scheys, L.; Desloovere, K.; Suetens, P.; Jonkers, I. Level of Subject-Specific Detail in Musculoskeletal Models Affects Hip Moment Arm Length Calculation during Gait in Pediatric Subjects with Increased Femoral Anteversion. J. Biomech. 2011, 44, 1346–1353. [Google Scholar] [CrossRef]
  33. Barzan, M.; Modenese, L.; Carty, C.P.; Maine, S.; Stockton, C.A.; Sancisi, N.; Lewis, A.; Grant, J.; Lloyd, D.G.; Brito Da Luz, S. Development and Validation of Subject-Specific Pediatric Multibody Knee Kinematic Models with Ligamentous Constraints. J. Biomech. 2019, 93, 194–203. [Google Scholar] [CrossRef]
  34. Soodmand, E.; Kluess, D.; Varady, P.A.; Cichon, R.; Schwarze, M.; Gehweiler, D.; Niemeyer, F.; Pahr, D.; Woiczinski, M. Interlaboratory Comparison of Femur Surface Reconstruction from CT Data Compared to Reference Optical 3D Scan. Biomed. Eng. Online 2018, 17, 29. [Google Scholar] [CrossRef]
  35. Wu, G.; Siegler, S.; Allard, P.; Kirtley, C.; Leardini, A.; Rosenbaum, D.; Whittle, M.; D’Lima, D.D.; Cristofolini, L.; Witte, H.; et al. ISB Recommendation on Definitions of Joint Coordinate System of Various Joints for the Reporting of Human Joint Motion--Part I: Ankle, Hip, and Spine. International Society of Biomechanics. J. Biomech. 2002, 35, 543–548. [Google Scholar] [CrossRef]
  36. Valente, G.; Pitto, L.; Testi, D.; Seth, A.; Delp, S.L.; Stagni, R.; Viceconti, M.; Taddei, F. Are Subject-Specific Musculoskeletal Models Robust to the Uncertainties in Parameter Identification? PLoS ONE 2014, 9, e112625. [Google Scholar] [CrossRef]
  37. Turcotte, C.M.; Rabey, K.N.; Green, D.J.; McFarlin, S.C. Muscle Attachment Sites and Behavioral Reconstruction: An Experimental Test of Muscle-bone Structural Response to Habitual Activity. Am. J. Biol. Anthropol. 2022, 177, 63–82. [Google Scholar] [CrossRef]
  38. Andriacchi, T.P.; Mündermann, A.; Smith, R.L.; Alexander, E.J.; Dyrby, C.O.; Koo, S. A Framework for the in Vivo Pathomechanics of Osteoarthritis at the Knee. Ann. Biomed. Eng. 2004, 32, 447–457. [Google Scholar] [CrossRef]
  39. Branch, E.A.; Anz, A.W. Distal Insertions of the Biceps Femoris: A Quantitative Analysis. Orthop. J. Sports Med. 2015, 3, 2325967115602255. [Google Scholar] [CrossRef]
  40. Lohrer, H.; Arentz, S.; Nauck, T.; Dorn-Lange, N.V.; Konerding, M.A. The Achilles Tendon Insertion Is Crescent-Shaped: An in Vitro Anatomic Investigation. Clin. Orthop. Relat. Res. 2008, 466, 2230–2237. [Google Scholar] [CrossRef]
  41. Lee, J.-H.; Kim, K.-J.; Jeong, Y.-G.; Lee, N.S.; Han, S.Y.; Lee, C.G.; Kim, K.-Y.; Han, S.-H. Pes Anserinus and Anserine Bursa: Anatomical Study. Anat. Cell Biol. 2014, 47, 127–131. [Google Scholar] [CrossRef]
  42. Olewnik, Ł.; Gonera, B.; Podgórski, M.; Polguj, M.; Jezierski, H.; Topol, M. A Proposal for a New Classification of Pes Anserinus Morphology. Knee Surg. Sports Traumatol. Arthrosc. 2019, 27, 2984–2993. [Google Scholar] [CrossRef] [PubMed]
  43. Tenforde, A.S.; Yin, A.; Hunt, K.J. Foot and Ankle Injuries in Runners. Phys. Med. Rehabil. Clin. N. Am. 2016, 27, 121–137. [Google Scholar] [CrossRef] [PubMed]
  44. Ryan, J.M.; Harris, J.D.; Graham, W.C.; Virk, S.S.; Ellis, T.J. Origin of the Direct and Reflected Head of the Rectus Femoris: An Anatomic Study. Arthroscopy 2014, 30, 796–802. [Google Scholar] [CrossRef] [PubMed]
  45. Sedaghat, S.; Langguth, P.; Larsen, N.; Campbell, G.; Both, M.; Jansen, O. Diagnostic Accuracy of Dual-Layer Spectral CT Using Electron Density Images to Detect Post-Traumatic Prevertebral Hematoma of the Cervical Spine. Rofo 2021, 193, 1445–1450. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The process of reconstruction of bones and attachment markings is as follows: (A) Contour of the patellar tendon insertion (PTI) footprint at the anterior tibia after dissection by a surgical marker. (B) Painted footprint of the PTI with radiopaque paint. (C) Radiopaque markings of attachments during CT. (D) Rough surfaces after the segmentation process. (E) Smoothed and corrected bone surfaces with unedited paint surfaces. (F) Coherent model with all bones in position after fixation. (GI) Process of footprint reconstruction. (G) Tracing of contact area through a transparent paint surface on the bone. (H) Fully traced contact area on a bone. (I) Generated footprint surface on the bone.
Figure 1. The process of reconstruction of bones and attachment markings is as follows: (A) Contour of the patellar tendon insertion (PTI) footprint at the anterior tibia after dissection by a surgical marker. (B) Painted footprint of the PTI with radiopaque paint. (C) Radiopaque markings of attachments during CT. (D) Rough surfaces after the segmentation process. (E) Smoothed and corrected bone surfaces with unedited paint surfaces. (F) Coherent model with all bones in position after fixation. (GI) Process of footprint reconstruction. (G) Tracing of contact area through a transparent paint surface on the bone. (H) Fully traced contact area on a bone. (I) Generated footprint surface on the bone.
Life 14 00778 g001
Figure 2. Exemplary deviation analysis of a femur and variation of muscle attachment footprints. (A) Deviation of the segmented surface of CT4 to the reference surface of CT1. Underneath the paint for the vastus intermedius footprint the surfaces deviate strongest because paint and bone surfaces interblend on the CT images. Scale expressed in mm. (B) Coefficient of variation (CV) of each muscle in relation to the mean CV of the bones’ volumes, surface areas and the mean CV of all muscles.
Figure 2. Exemplary deviation analysis of a femur and variation of muscle attachment footprints. (A) Deviation of the segmented surface of CT4 to the reference surface of CT1. Underneath the paint for the vastus intermedius footprint the surfaces deviate strongest because paint and bone surfaces interblend on the CT images. Scale expressed in mm. (B) Coefficient of variation (CV) of each muscle in relation to the mean CV of the bones’ volumes, surface areas and the mean CV of all muscles.
Life 14 00778 g002
Table 1. Deviations of the segmented surfaces from the references surface after CT1 expressed in mm.
Table 1. Deviations of the segmented surfaces from the references surface after CT1 expressed in mm.
SpecimenFemur CT2Femur CT3Femur CT4Pelvis CT2Pelvis CT3Tibia CT2Tibia CT3Tibia CT4Patella CT4
10.380.380.610.520.600.630.340.680.60
20.400.450.710.400.390.470.450.490.92
30.370.340.660.310.320.690.710.740.74
40.340.390.640.320.340.350.300.650.54
Mean0.370.390.660.390.410.530.450.640.70
STD0.020.040.040.080.110.130.160.090.15
Table 2. Morphological bone parameters.
Table 2. Morphological bone parameters.
BonesVolumes [mm3] Surface Area [mm2] Volume’s Center of Gravity [mm]
Mean ± STDCV [%]Mean ± STDCV [%]Mean x ± STDMean y ± STDMean z ± STD
Pelvis3.93 × 105 ± 5.32 × 10413.547.11 × 104 ± 6.21 × 1038.7416.92 ± 6.27454.79 ± 40.497.12 ± 4.32
Femur6.45 × 105 ± 1.11 × 10517.237.20 × 104 ± 8.77 × 10312.175.75 ± 3.10198.85 ± 24.3216.72 ± 2.24
Patella2.41 × 104 ± 4.18 × 10317.384.76 × 103 ± 5.64 × 10211.8549.73 ± 2.126.49 ± 3.724.23 ± 0.52
Tibia/Fibula3.58 × 105 ± 6.21 × 10417.336.71 × 104 ± 8.99 × 10313.4013.31 ± 9.16167.38 ± 12.128.77 ± 5.64
Foot2.45 × 105 ± 3.49 × 10414.246.03 × 104 ± 6.58 × 10310.9244.91 ± 34.53454.51 ± 42.0828.18 ± 21.02
Table 3. Established morphological dimensions and axis of the femur.
Table 3. Established morphological dimensions and axis of the femur.
SpecimenFemoral Morphological Parameter [mm]
MATEAFHDSD
1430.4097.8050.6036.60
2403.7095.0054.7040.40
3442.1088.5051.3034.90
4345.3085.3047.1030.90
Mean405.3891.6550.9335.70
STD37.374.982.703.41
CV [%]9.225.445.309.56
Table 4. Surface area, coefficient of variation and mean centroid of the muscle attachment footprints.
Table 4. Surface area, coefficient of variation and mean centroid of the muscle attachment footprints.
Muscle AttachmentSurface Area [mm2]CV [%]Centroid [mm]
Mean ± STDMean x ± STDMean y ± STDMean z ± STD
Gastrocnemius Caput Mediale O630.12 ± 182.7129.005.41 ± 1.7224.21 ± 3.9619.97 ± 3.51
Gastrocnemius Caput Laterale O788.24 ± 293.7737.272.89 ± 2.1419.10 ± 8.4135.02 ± 3.09
Achilles Tendon I1000.21 ± 129.4712.94114.80 ± 41.42417.36 ± 41.4727.23 ± 17.40
Rectus Femoris Caput Rectum O282.42 ± 73.2125.9240.90 ± 4.76441.16 ± 39.094.05 ± 2.94
Rectus Femoris Caput Reflexum O366.77 ± 216.1758.944.90 ± 2.71441.88 ± 40.2318.03 ± 1.32
Biceps Femoris Caput Longum O260.42 ± 162.7762.5062.22 ± 8.87381.46 ± 41.807.38 ± 7.11
Biceps Femoris Caput Breve O582.20 ± 259.9944.662.71 ± 2.26144.29 ± 17.8219.04 ± 0.56
Biceps Femoris Tibial I103.86 ± 82.8779.798.13 ± 5.2446.70 ± 4.9043.34 ± 3.49
Biceps Femoris Fibular I401.88 ± 94.1623.4322.57 ± 6.1952.71 ± 6.3748.95 ± 4.33
Semimembranosus O471.78 ± 48.6910.3249.83 ± 9.45370.54 ± 41.355.23 ± 4.93
Semimembranosus I1048.10 ± 304.7729.0816.88 ± 3.2155.22 ± 1.8621.42 ± 0.73
Semitendinosus O604.32 ± 155.2525.6964.19 ± 10.99360.41 ± 41.7110.73 ± 10.21
Semitendinosus I226.50 ± 120.1153.0311.50 ± 5.6197.74 ± 10.527.08 ± 4.02
Sartorius O115.94 ± 19.9217.1862.48 ± 5.84472.53 ± 39.859.45 ± 8.04
Sartorius I313.18 ± 250.6980.0514.28 ± 3.3486.16 ± 18.825.06 ± 4.82
Gracilis O132.23 ± 46.1334.8816.69 ± 12.57360.09 ± 34.5379.06 ± 4.38
Gracilis I67.18 ± 22.7833.9014.80 ± 6.4185.25 ± 5.924.76 ± 3.73
Vastus Lateralis O2493.56 ± 82.113.295.53 ± 3.11238.36 ± 54.3319.55 ± 8.51
Vastus Medialis O3410.55 ± 541.0715.861.54 ± 1.36301.81 ± 35.9450.51 ± 4.47
Vastus Intermedius O10,968.48 ± 849.017.7424.07 ± 4.15224.31 ± 22.8131.96 ± 2.28
Quadriceps Tendon I1480.95 ± 271.9418.3656.66 ± 2.4210.52 ± 6.664.08 ± 0.65
Patellar Tendon O992.16 ± 305.2830.7750.81 ± 1.9515.20 ± 6.475.20 ± 1.82
Patellar Tendon I733.14 ± 198.9227.1324.11 ± 5.8673.44 ± 8.9416.76 ± 6.60
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Neumann, A.-M.; Kebbach, M.; Bader, R.; Hildebrandt, G.; Wree, A. Evaluation of 3D Footprint Morphology of Knee-Related Muscle Attachments Based on CT Data Reconstruction: A Feasibility Study. Life 2024, 14, 778. https://doi.org/10.3390/life14060778

AMA Style

Neumann A-M, Kebbach M, Bader R, Hildebrandt G, Wree A. Evaluation of 3D Footprint Morphology of Knee-Related Muscle Attachments Based on CT Data Reconstruction: A Feasibility Study. Life. 2024; 14(6):778. https://doi.org/10.3390/life14060778

Chicago/Turabian Style

Neumann, Anne-Marie, Maeruan Kebbach, Rainer Bader, Guido Hildebrandt, and Andreas Wree. 2024. "Evaluation of 3D Footprint Morphology of Knee-Related Muscle Attachments Based on CT Data Reconstruction: A Feasibility Study" Life 14, no. 6: 778. https://doi.org/10.3390/life14060778

APA Style

Neumann, A. -M., Kebbach, M., Bader, R., Hildebrandt, G., & Wree, A. (2024). Evaluation of 3D Footprint Morphology of Knee-Related Muscle Attachments Based on CT Data Reconstruction: A Feasibility Study. Life, 14(6), 778. https://doi.org/10.3390/life14060778

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop