Patients
This retrospective study was conducted to evaluate the changes in condylar position and morphology after mandibular reconstruction. Patients were eligible who had undergone mandibular reconstruction from January 2018 to December 2020 at No.1 unit of oral and maxillofacial surgery department, Nanjing Stomatological Hospital of the Medical School of Nanjing University (Nanjing, China). This study protocol was approved by the institutional ethics committee (2018NL-027(KS)).
The inclusion criteria were (1) mandibular defect with condylar preservation; (2) sequential CT documentation; (3) a follow-up period of at least 6 months; and (4) no denture for the occlusal defects caused by the surgery. Patients with preoperative temporomandibular joint (TMJ) dislocation, postoperative unstable occlusal relation, and inconsistent occlusal position during CT scan were excluded from the study.
In our grouping strategy, patients were allocated to the study group receiving CAS or to the control group with conventional surgery. All patients in both the study and control group underwent surgery performed by the same chief surgeon in NO.1 unit of the oral and maxillofacial surgery department.
Surgical procedures
Major procedures included tumor resection, bone flap harvest, and mandibular reconstruction. CAS was performed on all patients in the study group, including preoperative virtual surgery, real surgery, and postoperative accuracy analysis[14, 24, 25]. Patient-specific cutting and shaping guides and pre-bent surgical plates were designed and fabricated by UEG Medical (Shanghai, China). In the control group, all the operations were performed relying on the surgeon’s experience, including tumor resection, bone flap harvest and shaping, and grafts fixation. In situ fixation of the condylar head and fibula in both groups was accomplished using reconstructive titanium plates and (or) mini plates according to the flap shaping method.
Data Acquisition
Patients’ gender, age, pathology, type of defect, and resection range were documented. The mandibular defect was reported according to the Brown classification[26]. In the pre and postoperative accuracy analysis, the patient’s computed tomography (CT) data were used. CT scans were taken routinely by a spiral CT device (0.625mm slices) within 1 week before surgery (T0); 7 to 14 days after surgery (T1) to evaluate changes caused by surgery, and at least 6 months after surgery (T2) to assess long-term changes. Patients were instructed to bite their teeth together and stabilize their occlusion while the CT scans were taken. All patients were followed up for at least 6 months.
Parameters and measurement
The CT data were imported into Mimics software 20.0 (Materialise HQ, Leuven, Belgium) and Geomagic Studio 2013 (Geomagic®, Morrisville, North Carolina) for multiplanar reconstruction (MPR) and analysis. Various parameters were established to evaluate the condylar position and volume changes based on quantitative and qualitative measurements.
Initially, the Frankfort horizontal (FH) plane of all the CT images was adjusted perpendicular to the floor before the analysis. Then the condylar position was evaluated by two-dimensional images. The sagittal view, in which the condyle could be viewed with the widest transverse diameter, was selected as the reference view for measurement. The coronal plane image was selected by the same method.
Five parameters were established to evaluate the condylar position according to the methods proposed by Kamelchuk et al[27] and Tyan et al[28]. A sagittal plane image and an axial plane image were selected as the reference views for the measurement. In the sagittal view, point A and B were respectively defined as the most prominent anterior and posterior aspect of the condyle. Point C was defined as the most superior surface of the glenoid fossa. Line AA’, BB’, and CC’ were vertical to the line CA and CB, and FH plane respectively. The length of AA’, PP’, and CC’, representing the distance of anterior joint space (AJS), posterior joint space (PJS), and superior joint space (SJS) was measured (Fig. 1A). In the axial view, point F was defined as the most superior surface of the glenoid fossa. Point D and E were respectively the contact point on the medial and lateral surface of the condyle with line FD and FE, which were tangent to the medial and lateral surface of the condyle respectively. Line DD’ and EE’, vertical to the line FD and FE, represented the distance of medial joint space (MJS) and lateral joint space (LJS) (Fig 1B).
Four parameters were established to describe the morphological characteristics according to the methods proposed by Hoppenreijs et al[29] and Kinzinger et al[30]. In the sagittal view, Y-axis was tangent to posterior border of the condyle neck with the connect point G. X-axis was drawn perpendicular to Y-axis through the broadest point K of condyle. Top of condyle (point N) is related to midpoint (point M) on X-axis. The length of GK and MN represented the condylar width (CW) and condylar height (CH) respectively (Fig. 2A). In the axial view, the distance between the most prominent medial point Q and the most prominent lateral point R represented the condylar length (CL) (Fig. 2B). In terms of the condylar volume (CV), the contralateral condyle was dissected at the level of the sigmoid notch from the 3D mandible model, whereas the ipsilateral side was dissected at the level of the highest point of the remaining condylar section if the sigmoid notch was resected. The CV of all patients’ bilateral condyles was also calculated by Rhinoceros 5.0 (Robert McNeel & Associate, Seattle, USA). These four parameters were only measured in the CT images of T1 and T2.
All measurements were repeated 3 times, and the mean value was used for statistical analysis.
After the measurements, the condylar position was calculated respectively as ln(PJS/AJS) and ln(MJS/LJS) based on the method proposed by Pullinger and Hollender[31]. The condylar positions were divided into 3 categories in both sagittal and coronal views based on the results. In sagittal view, 1) concentric if ln(MJS/LJS) was at least -0.25 to no greater than 0.25; 2) medial if ln(MJS/LJS) was less than -0.25; or 3) lateral if ln(MJS/LJS) was greater than 0.25; In coronal view, 1) concentric if ln(PJS/AJS) was at least -0.25 to no greater than 0.25; 2) posterior if ln(PJS/AJS) was less than -0.25; or 3) anterior if ln(PJS/AJS) was greater than 0.25.
Finite element analysis
In this study, Finite element analysis was adopted to simulate the occlusal process and predict the deviation of mandible and condyle due to the imbalance caused by the reconstruction surgery.
CT scan data of a patient from the CAS group were served as the pattern for reconstructing a 3D model in Mimics software (Fig. 3). We excised the tumor on the left side of the mandible to obtain the defective mandibular and design the fibular grafts, titanium plates, and rivets in 3-Matic software (Materialise, Leuven, Belgium) according to the postoperative 3D model. The articular disc was obtained by offsetting the condylar surface by 2 mm[32, 33].
The constructed model above is meshed in 3-Matic, where the non-fluid assembly is used so that the titanium rivets and the mandible, bone graft, and titanium plate in contact with it are co-noded. All surfaces of the model were first meshed using triangular slices with side lengths set to 0.5. Linear tetrahedral (C3D4) was used to mesh the skull, articular disc, mandible, teeth, implants, and rivets.
The material properties of the mandible were set by gray values of CT images to obtain cortical bone and cancellous bone referring to previous literature[34]. The mean CT value and bone density of the cortical bone were set to 1600 Hu and 1.73 g/cm3 based on the reference value -1024 Hu and 0 g/cm3. The empirical formula proposed by Carter-Hayes and Zannoni[34] for density and elastic modulus was E=kρ3, k=4249 GPa (g/cm3) [21, 32, 35, 36]. All mechanical properties were set in Mimics software referring to the values in the previous studies (Fig. 4).
Table 1 Mechanical properties and Masticatory muscle forces set in Mimics software
Mechanical properties
|
Loading
|
Materials
|
Young's modulus (MPa)
|
Poisson's ratio
|
Muscles
|
Maximum force (N)
|
Bone grafts
|
7300
|
0.3
|
Superficial masseter
|
272.8
|
|
|
|
Deep anterior masseter
|
73.8
|
Titanium rivets and plates (Ti6Al4V)
|
110000
|
0.3
|
Deep posterior masseter
|
65.8
|
Anterior temporalis
|
308.0
|
Posterior temporalis
|
222.0
|
|
|
|
Medial pterygoid
|
240.0
|
Articular disc
|
44.1
|
0.4
|
Superior lateral pterygoid
|
38.0
|
Inferior lateral pterygoid
|
112.8
|
In this study, masticatory muscles were simulated as linear-elastic spring element- SpringA and added between the skull and mandible according to their anatomical location (Fig. 5). The stiffness values of springs were taken from the literature[37]: masseter muscle=16.35 N/mm, lateral pterygoid muscle=12 N/mm, medial pterygoid muscle=15 N/mm, and temporalis muscle=14 N/mm. The maximum muscle force was taken from the relevant research[38] (Table 1). The masseter muscle, medial pterygoid muscle, and temporalis muscle of the ipsilateral side were removed in this model due to the surgical resection of the corresponding mandible and soft tissue. The ligaments connected to the skull behind the mandibular condyles were also modeled as springs connected to the fixed point. A square was placed on the right molar to simulate food, and its elastic modulus and Poisson's ratio were set to 100 MPa and 0.4.
The contact between the condyle and the articular disc was set as frictionless contact in the model. The contact between the titanium rivets and the plates, the mandible, and the bone graft was bound. The contact was set as penalty contact both between the bone graft and the mandible and between the tooth and the simulative food. The top of the articular disc and the upper surface of the square were fixed in all degrees of freedom.
The final model was imported into ABAQUS (Dassault Systemes, Cedex, France) for subsequent finite element analysis.
Statistical analysis
Data were statistically analyzed using IBM SPSS Statistics 20.0 (IBM Corp, Armonk, NY). A systematic error was assessed using paired t-test, and the random error was calculated using the formula of Dahlberg. Continuous data were recorded as the mean with standard deviation (SD); comparisons were made using the two-sample t-test or one-way analysis of variance (ANOVA). Categorical data were recorded as the count with percentage; comparisons were made using the χ2test or Fisher’s exact test. All statistical tests were two-sided, and significance was defined as a P-value <0.05.