key: cord-0003760-c4vpg8sw authors: Zhang, Jianying; Sun, Xiangzhao; Jia, Sen; Jiang, Xin; Deng, Tiange; Liu, Ping; Hu, Kaijin title: The role of lateral pterygoid muscle in the traumatic temporomandibular joint ankylosis: A gene chip based analysis date: 2019-03-22 journal: Mol Med Rep DOI: 10.3892/mmr.2019.10078 sha: 2a2a9ae443755df9e1d930c75c2023f0b1780a51 doc_id: 3760 cord_uid: c4vpg8sw Traumatic temporomandibular joint ankylosis (TMJA) is a common disease and disorder of the temporomandibular joint (TMJ); however, its pathogenesis has yet to be completely elucidated. In the authors' previous studies, the lateral pterygoid muscle (LPM) was confirmed to exert a function in distraction osteogenesis (DO) during the healing of a condylar fracture, which resulted in the formation of excess bone. The aim of the present study was to investigate alterations in the expression of any associated genes via an Affymetrix GeneChip method. The traumatic TMJA model was fabricated by a condylar fracture in the TMJ area of sheep with either a dissected LPM (LPD) or normal (LPN). The untreated sheep served as a control. At 4- and 12 weeks post-surgery, the condylar zone was isolated to perform the gene chip analysis, which was performed according to a standard Affymetrix protocol. The validated genes were further evaluated by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). The gene chip analysis indicated that the LPN gene expression pattern was similar compared with the DO process, while LPD was similar to that of normal bone fracture healing. The validated genes were collagen type II α1 chain, C-type lectin domain family 3 member A, interleukin 1A, cartilage oligomeric matrix protein, chondromodulin (LECT1), calcitonin receptor (CALCR), transforming growth factor (TGF)-β1, Fos proto-oncogene (FOS), bone γ-carboxyglutamate protein and bone morphogenic protein (BMP)7, among which, BMP7, LECT1, CALCR and FOS were confirmed by RT-qPCR. In conclusion, the present study demonstrated that LPM exerts a DO effect during the pathogenesis of traumatic TMJA, which may provide a novel target for preventing TMJA. The temporomandibular joint (TMJ) is crucial for normal mandibular movement, which is associated with numerous vital physiological processes, including chewing, talking and functions that require the opening of the mouth (1) . TMJ ankylosis (TMJA) restricts the movement of the mandible, which results in severe health problems in patients (2) . TMJA can result from various factors including trauma, infection, congenital issues and arthritis, among these, trauma in the TMJ region is the most common cause, such as a condylar fracture (3) . There are several hypotheses concerning the pathogenesis of traumatic TMJA. For example, it has been suggested that the abnormal healing of the fractured condyle is the principal reason for TMJA in the bone (4) . Specifically, injury in the TMJ region may result in haematoma, fibrosis and the abnormal formation of excessive bone (5) . A previous study hypothesized that the intra-articular bone fragment may help to increase the extent of ankylosis (6) . Based on the anatomy of the lateral pterygoid muscle (LPM), it can be inferred that the contraction of the LPM will generate sufficient distraction force to the condyle neck during the opening of the mouth and confirm the 'distraction osteogenesis (DO) effect of LPM' hypothesis, which is that the LPM may exert a DO effect during the healing of the fractured condyle (7) . In addition, the LPM exerts an important function in the reconstruction of the sagittal fractured condyle (8) and a previous study has reported that TMJA will be significantly alleviated if the LPM is dissected (9) . Furthermore, the study of Dai et al (10) also suggested that local injection of botulinum toxin A would temporarily block impulses from the LPM's and may effectively prevent the occurrence of traumatic TMJA. Consequently, blocking the LPM may serve as a promising therapeutic strategy for the prevention of traumatic TMJA. Unfortunately, the associated molecular mechanisms of traumatic TMJA and the effects of LPM have yet to be completely elucidated, which has limited the development of suitable therapeutic strategies. In the present study, a gene chip technique was used to systematically investigate any alterations in the expression of whole mRNA during the development of traumatic TMJA caused by a condylar fracture. In addition, the influence of LPM in the gene expression profile was also screened, in an attempt to identify the key molecules associated with this process. Animal model establishment and examination. Small-tailed Han sheep were provided by the Laboratory Animal Center of the Fourth Military Medical University (Xi'an, China) and all animal experiments were approved by the animal welfare ethics committee of the School of Stomatology, the Fourth Military Medical University. A total of 6 male sheep (aged 6 months, 25-35 kg) were kept under 20˚C, 40-60% humidity with a circadian rhythm (12:12 h light:dark cycle). The animals were fed with a standard diet and had free access to water. The traumatic TMJA model was constructed according to the protocol outlined in previous publications (9, 11, 12) . Briefly, the bilateral TMJ capsule was exposed from an incision to the anterior tragus. The joint capsule was incised and the condyle was exposed, and underwent a sagittal fracture by piezosurgery. Meanwhile, the outside quarter of the joint disc was dissected and a '#' shaped groove was created by piezosurgery on the fossa surface. In addition, the right side LPM was dissected (LPD) while the left side was normal (LPN) i.e. without dissection. The tissue was tightly sutured apart from the capsule and penicillin (1 mg pure penicillin G sodium per sheep for 3 days) was used to prevent infection. Sheep who did not receive surgery i.e. untreated normal sheep (UTR) were used as a blank control. The surgical procedures are illustrated in Fig. S1 . A total of 6 sheep were used and randomly divided into 2 groups of 3, which comprised the 4-and 12 weeks evaluation groups. In each group, one sheep was randomly selected as the UTR and the remaining sheep were subjected to surgery. For each animal in the group of surgical sheep, the right side of the TMJ was the LPD group while the left side was the LPN group. Gene chip RNA preparation. At either 4 or 12 weeks (8), the sheep were sacrificed using an overdose injection of pentobarbital sodium, and the bone tissues within the condylar fractured zones were separated and immediately frozen in liquid nitrogen (-195.79˚C ). Then the samples were subjected to total RNA extraction, purification and quality testing according to a standard Affymetrix protocol by Shanghai Biotechnology Corporation (Shanghai, China) as reported elsewhere (13) . The sample quality test is provided in Table SI . Gene chip analysis. The prepared RNA samples were subjected to gene chip hybridization, scanning and normalization according to the standard Affymetrix protocol. The differential genes were screened by a threshold method and the genes with a fold-change >2 were considered as differential genes (e.g. the ratio >2 or ratio <0.5). In addition, the data was further analyzed by cluster heat map (Cluster 3.0; http://bonsai. hgc.jp/~mdehoon/software/cluster/software.htm), scatter plot (R 3.5.0; https://www.r-project.org/), Database for Annotation, Visualization and Integrated Discovery (DAVID) Bioinformatics Resources 6.8 Platform (https://david.ncifcrf.gov/) and Gene Ontology (GO) term enrichment (Bioconductor 3.8 in R 3.5.0; http://bioconductor.org/), according to the respective manufacturer's protocol. The abbreviations for the biological processes included in the GO term enrichment analysis were listed in Table SII . Genes were identified as significantly differentially expressed by GO term enrichment analysis (P<0.05). Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) RNA preparation. The samples were frozen in liquid nitrogen for 48 h prior to being ground in a cold mortar. The total RNA was extracted by miRCURY RNA Isolation kit-tissue (Exiqon; Qiagen, Inc., Valencia, CA, USA) according to the manufacturer's protocol. The RNA concentration was measured by NanoDrop 2000 Spectrophotometer (Thermo Fisher Scientific, Inc., Waltham, MA, USA) and 1 µg total RNA was used for RT by the PrimeScript RT reagent kit (Takara Bio, Inc., Otsu, Japan) according to the manufacturer's protocol. Briefly, RT was conducted at 37˚C for 15 min, followed by 5 sec at 85˚C to terminate RT. Samples were then stored at 4˚C. RT-qPCR. The cDNA was amplified by SYBR Premix Ex Taq II (Takara Bio, Inc., Otsu, Japan) in order to perform the RT-qPCR procedure according to the manufacturer's protocol. qPCR was performed using an Applied Biosystems 7500 Fast Real-Time PCR System (Applied Biosystems; Thermo Fisher Scientific, Inc.) under the following conditions: 95˚C for 30 sec, followed by 40 cycles of 95˚C for 5 sec and 60˚C for 30 sec, and finally the melt curve stage (95˚C for 15 sec, 60˚C for 1 min and 95˚C for 15 sec). The associated primers were synthesized by Takara Bio, Inc., and RPL19 was used as endogenous reference, which are listed in Table SIII . The relative gene expression was calculated by the 2 -ΔΔCq method (14) . Three independent experiments were repeated. Statistical analysis. The quantitative data are presented as the mean ± standard deviation of at least three experimental repeats. The one-way analysis of variance with Student-Newman-Keuls post hoc test was used to compare differences between groups within a certain gene (GraphPad Prism 5; GraphPad Software, Inc., La Jolla, CA, USA). P<0.05 was considered to indicate a statistically significant difference. Cluster heat map. The overall gene expression was classified by heat map analysis and several regions of interest (ROI; the most different areas) were selected from the complete map (Fig. 1A ). The cluster similarity was similar for the results obtained at 12 and 4 weeks in the same treatment groups in the majority of clusters ( Fig. 1) , suggesting that the gene expression pattern was similar in animals in the same treatment groups. Of note, the gene expression in LPD was up-or downregulated in the opposite direction to LPN, at 4 and 12 weeks (Fig. 1B ). In addition, the genes associated with osteogenesis, including interleukin (IL)1A, collagen type II α1 chain (COL2A1), cartilage oligomeric matrix protein (COMP) and C-type lectin domain family 3 member A (CLEC3A) were included in these differential zones ( Fig. 2A-C) . The cluster heat map can directly display the general information about the expression levels of genes and the ROI can be easily selected from the differential genes zone. The data revealed that the cluster similarity was nearly identical between 4 and 12 weeks in each of the LPN, LPD and UTR groups, indicating that the different treatments can distinctively alter the gene expression patterns. There are certain specific muscle contraction associated genes including myosin light chain 2, myosin heavy chain (MYH)2, troponin C1 (TNNC1), MYH6 and actin α1 that were downregulated in LPD, suggesting that the LPM may serve an important role in the pathogenesis of traumatic TMJA. As bony ankylosis was the principal focus of the present study, the genes associated with osteogenesis were of particular interest. Based on the heat map, the IL1A, COL2A1, COMP and CLEC3A were preliminarily selected as potential critical genes. Scatter plot. The scatter plot was used to illustrate the genes expression distribution between two different groups. At 4 weeks, it was evident that the majority of genes in LPN were significantly upregulated, irrespective of whether it was being compared with LPD or UTR ( Fig. 3A and B) . However, the LPD and UTR exhibited near-identical expression patterns (Fig. 3C) . Additionally, at 12 weeks, the LPN group exhibited an increased number of downregulated genes when compared with the LPD group (Fig. 3D ) and there were more upregulated genes in the LPN group compared with the UTR group (Fig. 3E ). Whereas the LPD group had more upregulated genes compared with UTR at the 12 weeks stage (Fig. 3F) . The scatter plot can give direct insight into the number of up-and downregulated genes between two groups. Generally, in the present study, the LPN group exhibited a higher expression of active genes in the early stage, whereas more genes were activated in the LPD group at the later stage. Although the exact number of altered genes cannot be ascertained from the scatter plot, it corroborated the authors' hypothesis that genes associated with osteogenesis are dominant. During the process of DO, rapid osteogenesis is occurring in the first month, which may limit the degree to which the mouth can open and decrease the force of LPM (15) . As a result, the DO effect of LPM is weakened and returns to a stable state. Therefore, most genes are upregulated at an early stage and are then subsequently downregulated. Conversely, the dissection of LPM makes it similar to the normal bone fracture healing process and the majority of the osteogenesis associated genes are activated in a relatively late stage. However, the exact differential genes involved in this process require further analysis. DAVID and GO enrichment analysis. In order to determine the specific differential genes in differing treatments, DAVID gene classification and GO enrichment analysis were performed. The DAVID Bioinformatics Resources can be used to categorize the differential genes through different filters, which is extensively used for systematical analysis of large gene lists (16) (17) (18) . The GO term enrichment analysis can group the differential genes by different biological processes and determine which genes are significant (19, 20) . The genes were categorized by DAVID Bioinformatics Resources and more than 2-fold altered genes were considered different. Then the representative specific osteogenesis associated genes were listed. In LPN vs LPD at 4 weeks, the COL2A1, COMP, CLEC3A, calpain 3 and fibromodulin genes were upregulated up to >3-fold (Table I) . From the GO term enrichment data, the bone morphogenesis and endochondral bone growth were the most notable, in which COMP and (Table II) . It has previously been reported that the bone formed by DO is formed via intramembranous ossification (21) , in which COL2A1 is involved (22) . Taken together, these data indicate that the LPN group may be similar to the DO process and the genes associated with intramembranous ossification are activated at early stages. In addition, muscle contraction was also significantly upregulated in the LPN compared with the LPD group (P<0.05; Table II ), which may be attributed to the function of LPM. In LPN vs. UTR at 4 weeks, not only COMP and CLEC3A were upregulated but cluster of differentiation (CD)180, Toll-like receptor 7, histidine rich calcium binding protein and ryanodine receptor 1 (RYR1) genes were also significantly increased, while the expression of S100 calcium binding Fold protein A8 (S100A8), adiponectin (ADIPOQ) and BMP3 were downregulated (Table I ). In the GO-term enrichment data, the nuclear factor-κB transcription factor activity and ossification in bone mature and cellular response to chemical stimulus were significantly different, which also contained the S100A8, CD180, RYR1 and ADIPOQ genes (P<0.05; Table II ). In LPD vs. UTR at 4 weeks, the upregulated genes were FOS, matrix metalloproteinase 12 and IL1 receptor antagonist while the CLEC3A, COMP, COL2A1, BMP3, LECT1 and CALCR were markedly downregulated (Table I ). In the GO-term enrichment data, the ossification, cartilage development and response to external stimulus were more evident (Table II) . The COMP and COL2A genes were also revealed to be important genes in endochondral bone growth, ossification and cartilage development (Table II) . Since the cartilage development and endochondral ossification are the main processes during normal bone fracture healing (23, 24) , it may be inferred that LPD is similar to the normal healing of condylar fractures, which repair slowly over 4 weeks and therefore this explains why genes associated with osteogenesis including COMP, COL2A1, and BMP3 were downregulated. In LPN vs. LPD at 12 weeks, the ADIPOQ, myosin light chain, phosphorylatable, fast skeletal muscle, TNNC1 and calsequestrin 1 were upregulated while the CD180, FOS, LECT1 and IL1A were downregulated (Table III) . In this stage, the response to stress, the immune response, cellular response to chemical stimulus and response to external stimulus were significant, among which the FOS and IL1A exhibited multiple functions in the regulation of this process (Table IV) . Lewinson et al (25) reported that the expression of FOS is in line with the DO process, which is upregulated in the early stages and attenuated when the bone trabecula is formed. Consequently, the downregulation of FOS in LPN may suggest that the DO process is entering a stable stage. In LPN vs. UTR at 12 weeks, the osteogenesis associated genes of COMP, serpin family E member 1, CLEC3A, LECT1 were upregulated, while LEPR, BMP7, ADIPOQ, S100A9 and endothelin receptor type B were downregulated (Table III) . The genes of response to external stimulus, endochondral bone growth, cartilage development and bone development were significant, among which COMP was involved in multiple bone morphogenesis processes (Table IV) . This indicates that although the DO effect at 12 weeks is at a stable stage, the bone . LPD RS CHI3L1, TNFAIP3, MB21D1, IL1B, NR4A2, LYZ2, S100A9, SOD2, NOS2, CCL5, TLR2, MX2, EIF2AK2, DTX3L, IL1A, FAS, IER3, BPI, CD27, LAP, GBP5, RSAD2, CXCL16, IL1RN, ADIPOQ, CCL20, PARP9, DERL3, XDH, KMO, MT2A, LYZ,PTAFR, PLAC8, GADD45B, AQP9, MYF6, LY86, MYH2,TUBA4A, DOCK2, HERC5, FGF10, STC2, PML, SLC7A11,VCAM1, PRKCB, FOS, TCAP, TRIB1, NR4A3, ATF3, CCRL2,NR4A1, PTGS2, CASQ1, S100A8, S100A12, CD2, IFIT2, CYSLTR1, ACOD1, CCL7, CYBB, ISG15, DDX58, RFC3, ITGB2, CXCR4, PARPBP, CXCL8, MX1, TLR7, TLR8 IR TNFAIP3, MB21D1, IL1B, S100A9, NOS2, CCL5, TLR2, MX2, IL1A, FAS, GBP5, RSAD2, CXCL16 COMP, COL2A1 BM COMP, HAS2, FREM1, SFRP2, COL2A1 RS CHI3L1, TNFAIP3, IL1B, DUSP10, CAMP, NR4A2, GATA6, LYZ2, S100A9, SOD2, NOS2, CCL5, ITIH4, TLR2, EDN1, CXCL13, IL1A, FAS, IER3, LYZ1, BPI, CD27, LAP, CXCR6, CXCL16, IL1RN, ADIPOQ, HAS2, DERL3, XDH, CCL21, KMO, MT2A, CATHL2, LYZ,C7, PGLYRP1, SFRP1, HGF, SERPINE1, TNFSF8, LY86,LTF, MYH2, EGLN3, FABP4, DOCK2, TBX3, FOXF1,SELP, STC2, SFRP2, STAR, SLC7A11, FOS, CLEC7A, OLR1, TRIB1, NR4A3, DEFB128, LGMN, PLAUR, ATF3, C5AR1, HP, NR4A1, PTGS2, SGK1, S100A8, S100A12,CTSS, EDNRB, ACOD1, CCL7, CYBB, IFNA21, CAPN3, ITGB2, CXCR4, MYOCD, CD163, EYA4, LIPA, TLR10, CXCL8, CD14 Biological process abbreviations are defined in Table SII . BP, biological process; LPN, the lateral pterygoid muscle was normal; LPD, the lateral pterygoid muscle was dissected; UTR, Sham surgery without treatment. morphogenesis activity in LPN remains increased compared with in normal sheep. In LPD vs. UTR at 12 weeks, the osteogenesis associated genes COL2A1, IL1A, CLEC3A, FOS, CD180, COMP and LECT1 were upregulated and CALCR and ADIPOQ were downregulated (Table III) . The GO-term enrichment also suggested that endochondral bone growth, ossification and bone morphogenesis were significant processes (Table IV) . COMP and COL2A1 were revealed to be crucial genes in these biological processes (Table IV) . The upregulation of genes associated with osteogenesis indicates that during this stage, the LPD bone metabolism is more active, which exerts similarities to a normal bone fracture. In conclusion, through the gene chip analysis, the COL2A1, CLEC3A, IL1A, COMP, LECT1, CALCR, FOS and BMP7 genes were preliminarily selected as potentially critical genes during the pathogenesis of traumatic TMJA. However, due to the precision limitations of the gene chip technique, it should be taken into consideration that these validated genes may only be partially activated, or not the key vital genes during the establishment of traumatic TMJA. Therefore, the screened genes still require further analysis. The genes of interest i.e. BMP7, CALCR, FOS and LECT1 were confirmed by RT-qPCR. In LPN vs. LPD at 4 weeks, BMP7 (~1.5-fold), CALCR (~2.5-fold) and LECT1 (>5-fold) were significantly upregulated while FOS (~0.5-fold) was downregulated (P<0.05; Fig. 4A ). In LPN vs. UTR at 4 weeks, BMP7 (~0.75-fold) was significantly reduced while the FOS (~5-fold) and LECT1 (~5-fold) were upregulated (P<0.05; Fig. 4A ). The expression of CALCR exhibited no statistically significant difference between the LPN and UTR group (Fig. 4A) . However, for LPD vs. UTR at 4 weeks, BMP (<0.5-fold), CALCR (~0.2-fold) and LECT1 (~0.1-fold) were decreased while the FOS (~10-fold) was significantly improved (P<0.05; Fig. 4A ). In LPN vs. LPD at 12 weeks, the CALCR (<2-fold) was upregulated whereas the FOS (~0.1-fold) and LECT1 (~0.5-fold) were significantly inhibited (P<0.05; Fig. 4B ). The BMP7 results exhibited no statistically significant difference in the LPN group compared with the UTR (Fig. 4B) . In LPN vs. UTR at 12 weeks, BMP7 (~0.5-fold) was significantly downregulated whereas FOS (~1.5-fold) and LECT1 (>20-fold) were upregulated (P<0.05; Fig. 4B ). There was no difference in the expression of CALCR in LPN vs. UTR at 12 weeks (Fig. 4B) . However, in LPD vs. UTR at 12 weeks, BMP7 (<0.5-fold) was downregulated while the expression of FOS (>10-fold) and LECT1 (>40-fold) were markedly upregulated. In conclusion, the expression of BMP7 and LECT1 confirm the result from the chip data while the CALCR and FOS are slightly different. It is not uncommon that differences between Affymetrix GeneChip data and RT-qPCR measurements are detected. This may be due to the non-specific hybridization of genes and the fact that these processes rely on different algorithms (26) . Nevertheless, the gene chip data is considered accurate (27) and a previously study has also suggested that the gene chip data should be coupled with RT-qPCR analysis (28) . In conclusion, the present study demonstrated that normal LPM can significantly increase the expression of osteogenesis associated genes, which is similar to DO gene expression profiles. However, the dissection of LPM can significantly delay the expression of these genes, which is similar to the results observed in normal bone fracture healing. Numerous genes including COL2A1, CLEC3A, IL1A, COMP, LECT1, CALCR, FOS and BMP7 may serve vital roles during the pathogenesis of traumatic TMJA, and warrant further investigation in future studies. Temporomandibular joint disorders' impact on pain, function, and disability Surgical treatment of posttraumatic ankylosis of the TMJ with different pathogenic mechanisms Pathogenesis of post-traumatic ankylosis of the temporomandibular joint: A critical review Analysis of temporomandibular joint ankylosis caused by condylar fracture in adults A protocol for management of temporomandibular joint ankylosis The effect of an intra-articular bone fragment in the genesis of temporomandibular joint ankylosis A new hypothesis of mechanisms of traumatic ankylosis of temporomandibular joint The role of the lateral pterygoid muscle in the sagittal fracture of mandibular condyle (SFMC) healing process Influence of the lateral pterygoid muscle on traumatic temporomandibular joint bony ankylosis Injection of botulinum toxin A in lateral pterygoid muscle as a novel method for prevention of traumatic temporomandibular joint ankylosis Surgical induction of temporomandibular joint ankylosis: An animal model A sheep model for temporomandibular joint ankylosis Explaining differences in saturation levels for Affymetrix GeneChip arrays Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method Lateral pterygoid muscle volume and migraine in patients with temporomandibular disorders Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources The DAVID Gene Functional Classifcation Tool: A novel biological module-centric algorithm to functionally analyze large gene lists DAVID Bioinformatics Resources: Expanded annotation database and novel algorithms to better extract biology from large gene lists GOrilla: A tool for discovery and visualization of enriched GO terms in ranked gene lists Ontologizer 2.0-a multifunctional tool for GO term enrichment analysis and data exploration Biology of bone formation, fracture healing, and distraction osteogenesis The role of chondrocytes in intramembranous and endochondral ossification during distraction osteogenesis in the rabbit Forward deviation of the mandibular condyle enhances endochondral ossification of condylar cartilage indicated by increased expression of type X collagen Effect of alendronate on endochondral ossification in mandibular condyles of growing rats Stimulation of Fos-and Jun-related genes during distraction osteogenesis Probe set filtering increases correlation between Affymetrix GeneChip and qRT-PCR expression measurements Direct comparison of GeneChip and SAGE on the quantitative accuracy in transcript profiling analysis Coupling multiplex RT-PCR to a gene chip assay for sensitive and semiquantitative detection of severe acute respiratory syndrome-coronavirus The authors would like to thank Dr Yang Xue (Department of Oral Surgery, School of Stomatology, The Fourth Military Medical University, Xi'an, China), for her help during the experiment design and data analysis. The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. KH designed the study and provided financial support. JZ and XS participated in all the animal experiments as well as the gene chip analysis. SJ was involved in interpreting the data and drafted the manuscript. XJ, TD and PL contributed to the animal surgery and data analysis. All authors reviewed and approved the manuscript for submission. All the animal experiments were approved by the animal welfare ethics committee of the School of Stomatology, the Fourth Military Medical University (Xi'an, China) in accordance with the relevant guidelines and regulations. Not applicable. The authors declare they have no competing interests.