Next Article in Journal
Cumulative Exposure to Neighborhood Conditions and Substance Use Initiation among Low-Income Latinx and African American Adolescents
Next Article in Special Issue
Remotely Supervised Exercise during the COVID-19 Pandemic versus in-Person-Supervised Exercise in Achieving Long-Term Adherence to a Healthy Lifestyle
Previous Article in Journal
Comparison of the Mental Health Impact of COVID-19 on Vulnerable and Non-Vulnerable Groups: A Systematic Review and Meta-Analysis of Observational Studies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 18
Issue 20
10.3390/ijerph182010828
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association of COMT Polymorphisms with Multiple Physical Activity-Related Injuries among University Students in China

by
Shangmin Chen
1,2,
Weicong Cai
1,3,
Shiwei Duan
4,
Lijie Gao
5,
Wenda Yang
1,
Yang Gao
6,
Cunxian Jia
5,
Hongjuan Zhang
7 and
Liping Li
1,2,*
1
Injury Prevention Research Center, Shantou University Medical College, Shantou 515041, China
2
School of Public Health, Shantou University, Shantou 515041, China
3
Shenzhen Center for Chronic Disease Control, Shenzhen 518020, China
4
Medical Genetics Center, School of Medicine, Ningbo University, Ningbo 315211, China
5
Department of Epidemiology, Shandong University School of Public Health, Jinan 250012, China
6
Department of Sport and Physical Education, Hong Kong Baptist University, Kowloon Tong, Hong Kong 999077, China
7
School of Physical Education, Shanxi Normal University, Linfen 041000, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(20), 10828; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph182010828
Submission received: 12 September 2021 / Revised: 8 October 2021 / Accepted: 12 October 2021 / Published: 15 October 2021
(This article belongs to the Special Issue Sport, Physical Activity and Health at Time of COVID-19)
Browse Figures
Versions Notes

Abstract

:
The catechol-O-methyltransferase (COMT) is a candidate gene to provide promising evidence of psychiatric disorders, but there is a knowledge gap between the genetic factor and multiple physical activity-related injuries (PARIs). The aim of this study was to explore the contribution of COMT to the risk of PARIs among university students in the Chinese Han population. We can further search for the intrinsic risk factors for the occurrence of multiple physical activity injuries and provide a scientific basis for early screening and precise intervention for the high-risk group of college students with multiple PARIs. A 1:1 matched case-control study of 61 PARIs cases and 61 healthy controls were carried out. DNA samples of the participants were isolated from saliva and genotyped on eight SNPs of the COMT gene (rs9265, rs4680, rs6269, rs4818, rs4633, rs165655, rs165656, and rs165722) using the MALDI-TOF MS method. We found that rs6269 and rs4818 were significantly associated with PARIs, and rs6269-GG and rs4818-GG contributed to the reduced risk of PARIs. Further haplotype analysis showed a four-marker C-G-C-G haplotype (rs165722-rs6269-rs4633-rs4818) acted with a protective role in the development of PARIs (p = 0.037; OR: 0.474, 95% CI: 0.269 to 0.834). However, the interactions between club membership and rs6269 or rs4818 would significantly increase the risk of PARIs (both p < 0.001, OR: 5.121 and 4.977, respectively). This is the first study to find the contribution of COMT to PARIs occurrence, suggesting that the COMT polymorphisms and the gene–environment interactions may alter the risk of PARIs.

1. Introduction

Participation in physical activity (PA) is conducive to our individual well-being and helps to achieve and maintain optimal health levels in terms of physical, cognitive and psychological factors [1,2]. In light of these beneficial effects, practically all countries and regions in the world place great emphasis on PA promotion (a public health priority). More positively, there exist several recommendations involving the frequency and duration of PA participation [3,4,5,6]. With such a contemporary focus on a physically active lifestyle, however, a potentially increasing exposure to physical activity-related injury (PARI)—the major adverse consequence of PA participation—can be expected [7,8]. Earlier studies related to the PARI incidence have been primarily confirmed among sports-active populations. For example, basketball, football, and volleyball players in the National Collegiate Athletic Association have injury incidence densities of 6.54–7.97, 8.07–8.44, and 4.69–7.07 per 1000 athlete-exposure, respectively [9,10,11]. Also, our previous study noted that about 22.7% of Chinese university students had at least one PARI episode during the past 12-month period [12]. In fact, PARI could cause substantial medical costs and indirect harmful effects such as PA restriction, psychological fear, and physical discomfort [13,14]. Researchers strongly believe that most PARI could be effectively prevented given the targeted prophylactic interventions, on basis of the better understanding of the injury characteristics and etiology for the specific population [15]. According to previous research, the majority of injuries involved the lower extremities, were acute, newly occurring injuries, and occurred in non-contact situations [9,10,11,16]. Gender and grade differences contributed markedly to the occurrence and severity of PARI [17,18]. In addition, age, body mass index (BMI), PA levels, and family background were also associated with PARI risk [8,17,19]. At present, there is no direct research evidence to clarify the molecular genetic mechanism of PARI. However, since the psychological process involved in injury is mainly related to cognitive flexibility, especially executive function, aggressive behavior, and impulsive risk-taking behavior [17,20,21], by searching for these psychological variables closely related to injury, some genes that may be related to injury will be discovered.
The catechol-O-methyltransferase (COMT) gene is one of the well-studied genes involved in the metabolism of catecholamine neurotransmitters, including dopamine, adrenaline, and noradrenaline [22,23,24,25,26]. The COMT is located at “ 22q11.1–q11.2” and is approximately 27 Kbp in size, with 345 identified polymorphisms [27,28,29]. To date, the most widely studied COMT variant is rs4680 in exon 4, as rs4680 plays a key role in enzymatic activity and cognitive prefrontal cortical function [30]. Owing to the functional G to A polymorphism, rs4680 can result in substitution of valine (Val) to methionine (Met) in codon 108 (in soluble COMT) or 158 (in membrane-bound COMT), generating differences in the ability of the enzymes [24,28]. The enzymatic activity of the Val variant is nearly two to four times higher than the Met variant [30,31]. In addition to Val108/158Met, there are other COMT variants associated with the function or expression of regulatory enzymes [32,33]. Furthermore, the important role of COMT haplotype involving several polymorphic sites has also been identified. Notably, the haplotype including four COMT SNPs (rs4818 and rs4680—Val108/158Met—in exon 4, rs4633 in exon 3, and rs6269 in P1 promoter) can alter mRNA secondary structure and thereby mediate protein expression levels [34]. Compared to the low-function haplotype ACCGval, the high-function haplotype GCGGval expresses an 18- to 25-fold higher enzyme activity, in line with the increased protein levels [34].
Coupled with the previous studies, most genetic association studies of COMT focused on a variety of mental or psychological disorders [35,36], suicide behavior [27], and osteoporotic fracture [37]. No study investigating the association between the COMT and PARIs occurrence was found. Other than the similar influencing non-genetic factors of PARI stated earlier, genetic susceptibility may also play an important part in the occurrence of PARIs—especially in those who experienced PARI repeatedly in the same body part. Therefore, we hypothesized that COMT polymorphisms may alter the risk of PARIs. To test this hypothesis, we explored eight COMT SNPs (i.e., rs9265, rs4680, rs6269, rs4818, rs4633, rs165655, rs165656, and rs165722) through a 1:1 matched case-control study. Consideration of haplotypes are more powerful and effective tools to investigate the effect of genetic variations and gene-environment interaction is becoming increasingly important in the epidemiologic studies [24,38,39], and we also performed the haplotype analysis of COMT polymorphisms, and the interaction of the SNPs within COMT gene and environmental factors. The goal of our study was to find the contribution of COMT to the development of PARIs in the Chinese Han population of general university students.

2. Materials and Methods

2.1. Participants

The two-stage design was carried out in four universities in two Chinese cities, namely Shantou and Jinan. Chosen from 3390 students in grades 1 through 3 who participated in the baseline survey during March and April, 1421 students (41.9%) consented to be followed up and completed face-to-face interviews during April and May, 2017. No significant differences in basic demographics (i.e., gender, age, and year level) were found between the subjects in two stages (all p > 0.05). The interviewed participants who experienced multiple PARI episodes (at least three injury episodes) in the past 12 months were included in the case group (i.e., PARIs) and those who did not suffer from PARI event during the same period were recruited in the control group (i.e., non-PARI). This was a 1:1 case-control design that matched gender and year level. Given that one participant did not consent to participate in the genetic sub-study (i.e., saliva sample collection), our study finally enrolled 122 individuals identified during a face to face interview as either PARIs cases (n = 61, 20.02 ± 1.23 years) and non-PARI controls (n = 61, 19.89 ± 1.02 years) (Figure 1).
This study was in accordance with the Declaration of Helsinki, whose protocol was approved by Shantou University Medical College Ethics Committee (SUMC-2016-22). All participants (except for one subject) understood the purpose of our study and completed a written informed consent after our trained personnel verbally explained the study.

2.2. Data Collection

The baseline information including gender, age, grade, club membership, screen time (including cellphone and computer) were collected. Moreover, height and weight, and waist and hip circumferences were measured at the same time on the spot, which were calculated as BMI and waist-hip ratio (WHR), respectively.
Adapted from the reliability-validated Chinese version of Children’s Leisure Activities Study Survey (CLASS) [40] and the short version of the Minnesota Leisure Time Physical Activity Questionnaire (MLTAQ) [41], all subjects’ habitual participation in PA were estimated using a series of standardized questions. It has a good reliability (Cronbach’s α = 0.772) in the present study. Participants were asked whether they involved in any PA such as basketball, football, volleyball, and swimming per week in the past 12 months. If yes, they were required to provide information on the frequency and duration of this PA on both a weekday and a weekend, respectively. The weekly PA volume was then calculated.
The term PARI was defined as any injury suffered from during periods of involvement in physical education class, sports activities, or leisure time PA, like sprain, strain, fracture, etc. In addition, a countable PARI episode in this study must occur in the past 12 months and meet at least one of the four judgment criteria, which were fully described in our earlier report [17]. All students were requested to report and count their PARI experiences based on the predefined criteria. In order to ensure that the reported PARI met the aforementioned requirement, those injured participants were asked to provide information about the time of injury occurrence and what types of PA they involved when the injury occurred. Moreover, the injured students were also asked to describe other injury-related details such as causes, types, body parts, places, etc., which contributed to the validation of the PARI measure too. In this study, PARIs cases were defined as those who experienced at least three PARI episodes during the previous 12-month period.
A validated questionnaire was used, and participants were surveyed face-to-face, one-on-one by investigators and had relevant indicators (including height, weight, waist circumference, and hip circumference) measured on-site. In order to reduce bias, improve the degree of cooperation of the survey subjects and ensure the quality of the survey, before filling out the questionnaire, the investigator explained the purpose, significance and precautions of the study, answered relevant questions raised by the students during the filling process, and investigation using recall-promoting interrogation techniques. In order to test the reliability of the data collected by face-to-face interviews, the interviewed information were validated and reliability-tested against the data taken by telephone survey among 50 students with a 1-week interval after their completion of the earlier interviews (average kappa coefficient = 0.775 ± 0.262).

2.3. SNP Genotyping

We designed the primers using Assay Design 3.1 software (Sequenom Inc., San Diego, CA, USA), whose details were presented in Supplementary Table S1.
All consenting participants contributed a 2-mL saliva sample using DNA saliva self-collection tube (Beijing Thinkout Sci-Tech Co., Ltd., Beijing, China). DNA was extracted and purified from all saliva samples via a DNA purification kit (Beijing BioTeke Co., Ltd., Beijing, China). DNA fragments were amplified by the polymerase chain reaction (PCR), whose reaction volume was 5 μL containing 20 ng of genomic DNA template. The PCR product was digested with shrimp alkaline phosphatase (Sequenom Inc., San Diego, CA, USA). Following digestion, the amplified fragment was extended in a 9 μL reaction volume and resolved on a 4% agarose gel. All SNPs were genotyped using the high-throughput genotyping platform, conducted by matrix-assisted laser desorption ionization time of flight mass spectrometry (MALDI-TOF MS), and the data was collected via the TYPER 4.0 in the MassARRAY System (Sequenom Inc., San Diego, CA, USA). Furthermore, genotyping was repeated in 5% of all samples to verify the accuracy of the results, revealing that the proportion of concordance was greater than 99%.

2.4. Statistical Analysis

Descriptive statistics were used to evaluate the characteristics of the whole subjects. Continuous data that normally distributed were expressed as mean and standard deviation (SD), and the between-group differences were assessed by independent-sample t tests. Categorical variables were presented as number and percentage, and tested for differences by Pearson’s χ2 tests.
The power of the sample was calculated by the PS software for power and sample size calculation [42]. The obtained allele and genotype frequencies evaluated using the Hardy–Weinberg equilibrium (HWE) were tested by the χ2 tests for goodness-of-fit. Statistical significance for the between-group differences in allele, genotype and haplotype distribution was calculated by the χ2 tests. Whenever more than one fifth of the cells contained smaller than five subjects, a Fisher’s exact test was used. Odds ratios (ORs) and their corresponding 95% confidence intervals (95% CIs) were generated for associations between allele, genotypes and PARIs. We carried out dominant and recessive models as well, where crude and adjusted association were identified. The linkage disequilibrium (LD) between the SNPs was also calculated and haplotype blocks were identified via the CIs algorithm conducted in Haploview. Statistical analyses were performed with the PLINK 1.07, SPSS 23.0 (SPSS Inc. Chicago, IL, USA), and Haploview 4.0 (Broad Institute, Cambridge, MA, USA). We also applied the generalized multifactor dimensionality reduction (GMDR) using the GMDR 07 to identify the influence of the gene-environment interaction on the risk of PARIs.
All analyses with two-tailed p-value less than 0.05 were considered statistically significant.

3. Results

There is no difference in gender, club member, age, injury between people who agree to participate in the interview and those who refuse to participate (Table 1). Basic characteristics of the subjects with and without PARIs are presented in Table 2. A total of 60 males and 62 females were included in the present study, respectively. More students were in study year one (n = 56, 45.9%). The age of the whole subjects ranged from 18 to 23 years, with a mean age of 20.00 (SD: 1.127). There were no significant differences in the PARIs cases compared to the non-PARI controls in terms of age, BMI, WHR, and screen time (all p > 0.05). A significant difference was obtained in club membership between PARIs cases and non-PARI controls (χ2 = 30.500, p = 0.000). Also the cases had significantly longer duration in PA participation compared to the controls (1070.41 vs. 781.38 min per week on average, t = −2.846, p = 0.005).
The success genotyping rate of eight SNPs was between 93.2% and 99.4%. The genotype distribution of all SNPs within the COMT gene were in line with the HWE in the whole study populations (Supplementary Table S2). The frequencies of allele and genotype for eight COMT polymorphisms in cases and controls were summarized in the Table 3.
The genotype distribution of rs6269 differed significantly (χ2 = 6.047, p = 0.049). Although rare in both groups, rs6269-GG was significantly less found in cases (5.1% vs. 16.7%; OR: 0.200, 95% CI: 0.049 to 0.818). Regarding the allele frequency, the rs6269-G allele was significantly more common in controls (χ2 = 5.576, p = 0.018; 41.7% vs. 27.1%; OR: 0.521, 95% CI: 0.302 to 0.898). Only a moderate difference between cases and controls was observed in the genotype of rs4818 (χ2 = 5.788, p = 0.055), rs4818-GG homozygote was more frequent in the control group (14.8% vs. 5.2%; OR: 0.222, 95% CI: 0.054 to 0.923). In addition, more rs4818-G alleles were significantly detected in controls (χ2 = 5.385, p = 0.020; 41.0% vs. 26.7%; OR: 0.525, 95% CI: 0.304 to 0.908). There were no significant differences in either genotype or allele frequencies of the remaining COMT polymorphisms between cases and controls.
When applying a dominant or recessive model, we found that the frequency of G allele (i.e., GG + GA) carriers or GG carriers in rs6269 was lower in the case group (dominant model: 49.2% vs. 66.7%; χ2 = 3.745, p = 0.053; recessive model: 5.1% vs. 16.7%; χ2 = 4.101, p = 0.043). For the dominant model with no adjustments, a reduced risk of PARIs occurrence was found (OR: 0.483, 95% CI: 0.230 to 1.014). This relationship was diminished with the adjustment for gender, grade, and club membership (OR: 0.370; 95% CI: 0.150 to 0.914). After further adjustment in age, BMI, WHR, screen time, and PA participation, the results were further attenuated but remained significant (OR: 0.355; 95% CI: 0.137 to 0.918). There was no significant difference between cases and controls in the recessive model for rs6269 (AA + GA vs. GG). Similarly, the G allele (i.e., GG + GC) carriers or GG carriers of rs4818 in the dominant or recessive model were statistically underrepresented in the cases compared to the healthy control (dominant model: 48.3% vs. 67.2%; χ2 = 4.376, p = 0.036; recessive model: 5.2% vs. 14.8%; χ2 = 3.010, p = 0.083). A significant association with the lower risk of PARIs was found in a dominant model for rs4818-(GG + GC) vs. CC genotype without covariates (OR: 0.455, 95% CI: 0.217 to 0.956). After controlling various covariates, this association was still significant but weakened gradually. However, no significant difference was obtained in the recessive model for rs4818-(CC + GC) vs. GG genotype. Also, no dominant or recessive models of the remaining SNPs in two groups were found to be significantly correlated (Table 4).
Employing Haploview to calculate LD, COMT SNPs were found to be in strong LD with each other forming two haplotype modules (Figure 2). The first block was composed of rs165722, rs6269, rs4633, and rs4818, while the second block comprised rs9265 and rs165655. Only the C-G-C-G haplotype in block one was associated with a significantly reduced risk (OR: 0.474; 95% CI: 0.269 to 0.834). The higher frequency of the C-G-C-G haplotype was significantly found in the controls compared to the cases (39.3% vs. 26.8%; χ2 = 4.333, p = 0.037). No significant differences regarding the other COMT haplotypes in relation to PARIs could be detected (Table 5).
Investigation into the interaction between COMT SNPs (rs6269 and rs4818) and environmental factors (club membership and PA participation) on the risk of PARIs was performed by the GMDR approach, with gender, grade, age, BMI, WHR, and screen time as covariates. We identified a significantly important interaction model involving club membership and rs6269 based on 1000 times permutation test (cross-validation consistency, CVC: 9/10; test accuracy: 70.7%; p < 0.001), and the evidence of this interaction shown by the logistic regression model was 5.121 (95% CI: 1.938 to 13.532). In addition, club membership was also significantly interacted with rs4818 (CVC: 9/10; test accuracy: 71.5%; p < 0.001), with an effect of 4.977 (95% CI: 1.881 to 13.171). There were no significant interactions between PA participation and rs6269 or rs4818.

4. Discussion

The genetic studies into COMT could provide a promising contribution to the evidence of the pathophysiology, but there is a knowledge gap between the COMT gene and PARI (in particular PARIs occurrence). We conducted an association study involving 61 PARIs cases and 61 matched non-PARI controls. In this study, we found that rs6269 and rs4818 within the COMT gene were independently and significantly associated with PARIs, suggesting that both of them appeared to be involved in the development of PARIs in the Chinese Han population. Regardless of the adjustment for different covariates, the lower risk of PARIs occurrence was observed in the dominant inheritance model for both rs6269 and rs4818, whereas there was no significant association in the recessive model. These results further confirmed that rs6269-GG and rs4818-GG contributed to the reduced risk of PARIs. It is the first time these results of PARIs occurrence and COMT polymorphisms were found, and clearly more research is needed to clarify such an association more broadly.
Some researchers have proposed a “gene-brain-behavior” model, genes cause certain areas of the brain structure and function changes(e.g., hippocampus) [26,29], which affect the individual’s cognition, emotion and behavior [29,43]. Human behavior and the brain network that dominates behavior are formed by the ingenious combination of many genes through development and the environment. Prior research has demonstrated that there is a relationship between between the COMT rs4680 and the impact of COVID-19, which could plausibly be mediated by maladaptive anxiety-related behaviour [44]. Met allele is associated with poorer stress tolerance and the Val allele may confer an advantage in terms of threat perception, leading to better adherence to protective behaviour and minimizing the spread of the disease. The s4818 locus can explain more COMT activity variation, and the G allele at rs4818 locus means that the individual has a higher level of COMT enzyme activity and a lower level of prefrontal dopamine [45]. Studies in the field of cognition have found that subjects carrying the G allele at the rs4818 locus show higher problem-solving ability, but show lower efficiency in planning; The G allele at the rs6269 locus means high levels of COMT enzyme activity and low levels of prefrontal dopamine. Our research may contribute to early identification of individuals at risk in PARI, and suggests a mechanism that can be targeted in future studies aiming to prevent or limit negative outcomes. This may not only help us to reduce the injury during the COVID-19 pandemic, but also provide methods and ideas to explore the biological mechanisms of the impact of COVID-19.
When designing and interpreting association studies of phenotype and genotype, it is necessary to understand the haplotype structures in the human genome because haplotypes are more powerful and effective tools for investigating the potential effect of genetic variations on biological functions and disease susceptibility [38,46]. Many haplotype analysis studies suggested that the SNPs combination led to synergistic effects on enzymatic activity and protein function. [34] Moreover, the distribution of haplotypes varies between various settings and geographic locations and may affect the outcome [47]. In this study, we performed haplotype analysis based on the data for eight COMT polymorphisms generating two haplotype modules. Of these, only the inclusion of the C-G-C-G haplotype in the first block consisting of rs165722, rs6269, rs4633, and rs4818 exerted a significant association with the risk of PARIs (OR: 0.474; 95% CI: 0.269 to 0.834). Our result showed a significantly lower frequency of this haplotype in the PARIs group than the controls, suggesting that it acts as a protective role in the development of PARIs. This is the first study to explore the relationship between the COMT haplotype and PARIs occurrence, and the potential mechanism behind this protective effect of the low activity COMT haplotype on the development of PARIs is still unclear. Further work is needed to validate our findings and to build a more comprehensive view that better reflects the role the COMT haplotypes play in the development of PARIs.
The effects of a single allele might be small, but it might be overwhelmed by other factors [48]. More importantly, complex diseases are more likely to be affected by the gene–gene or gene–environment interactions [49]. For instance, earlier study assumed that the gene–environment interaction might relatively increase the risk for mental disorders with brain injury [48]. In the present study, using the GMDR method to detect the combined contribution of the gene and environment to PARIs occurrence, we found that the significant interactions between club membership and rs6269 or rs4818 played an important role in the development of PARIs. The results of the main effect of rs6269 and rs4818 showed that they were significantly associated with a lower risk of PARIs. Strikingly, when we considered the interactions of these two COMT SNPs with club membership, our results showed a dramatically increased risk of PARIs (OR: 5.121 and 4.977, respectively). The present findings provide a novel insight to the field of gene × environment that the significant joint contribution of club membership and rs6269 or rs4818 may be more responsible for the increased susceptibility of PARIs, helping to further note that the individual contribution of a specific polymorphism might be minor, but the gene–environment interactive effect could be rather significant. Moreover, we did not find significant interaction between PA participation and rs6269 or rs4818 on the PARIs occurrence although the main effects of COMT polymorphisms and PA levels on PARIs occurrence were found in the present study or previous reports [17]. Clearly, these findings warrant replication, and future research should take the gene–gene interaction into consideration as well.
This is the first study to explore the contribution of COMT to the risk of PARIs. But several limitations must be discussed when interpreting the present findings. First, the sample size was small, which may skew results and significance due to the lower statistical power. Thus, our study should be considered preliminary and hypothesis-generating. By contrast with the cross-sectional study design, we carried out a matched case-control study, whose cases and controls came from the same population. This enables us to better preclude the selection bias and obtain specific cause-effect relationship in the risk analysis. However, the participants in the present study were not randomly selected from the baseline survey. This may affect the generalizability of our findings somewhat. Therefore, our results need replication in future research, and further validation in larger sample sizes is undoubtedly necessary to confirm the genetic effect sizes and draw a firm conclusion from a good representative population. Second, although eight polymorphisms were included, we only analyzed the COMT gene in this study. We found that the gene–environment interaction could largely increase the risk of PARIs, but no gene–gene interaction was performed. Previous study noted that, if epistasis exists, the effect of one SNP may be affected by that of another SNP and thereby may impact the power to explore genetic association [50]. Thus, we should take the gene–gene interaction into account in future studies. Third, the information of PA participation was self-reported. The results of 1-week test–retest were validated to be reliable, but we could not fully preclude the possibility of the reporting bias. Hence, the PA measure may bias our gene-environment analyses. Finally, all participants were Han Chinese in order to control the effects of population stratification, but this means that it is not clear the findings can be applied to other races. Further validation in larger samples from other ethnic populations is thereby necessary in order to establish any firm conclusion.

5. Conclusions

In summary, this is the first study to find the association of COMT gene with the risk of PARIs. Our results noted that rs6269 and rs4818 were significantly associated with a lower risk of PARIs, especially the rs6269-GG and rs4818-GG. The haplotype of rs165722-C, rs6269-G, rs4633-C, and rs4818-G acted as a protective role in the development of PARIs. Moreover, the significant interactions of club membership with rs6269 or rs4818 could markedly increase the risk of PARIs.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijerph182010828/s1, Table S1. Basic information of primers for the SNPs in COMT gene; Table S2. Basic information of the eight SNPs within COMT gene.

Author Contributions

L.L., Y.G., C.J. conceived and designed the study, W.C., L.G. and W.Y. collected data and saliva samples, S.C., W.C., S.D. and H.Z. performed the statistical analyses, and W.C., S.C., S.D. and L.L. drafted the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China [grant number 31640038], the K.C. Wong Magna Fund in Ningbo University and the Guangdong Province Science and Technology Special Fund (Special Project and Task List) [grant number 200108155899596]. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

This study was in accordance with the Declaration of Helsinki, whose protocol was approved by Shantou University Medical College Ethics Committee (SUMC-2016-22).

Informed Consent Statement

All participants (except for one subject) understood the purpose of our study and completed a written informed consent after our trained personnel verbally explained the study.

Data Availability Statement

The datasets used and analyzed in the study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Garber, C.E.; Blissmer, B.; Deschenes, M.R.; Franklin, B.A.; Lamonte, M.J.; Lee, I.M.; Nieman, D.C.; Swain, D.P.; The American College of Sports Medicine. American College of Sports Medicine position stand. Quantity and quality of exercise for developing and maintaining cardiorespiratory, musculoskeletal, and neuromotor fitness in apparently healthy adults: Guidance for prescribing exercise. Med. Sci. Sports Exerc. 2011, 43, 1334–1359. [Google Scholar] [CrossRef]
  2. 2018 Physical Activity Guidelines Advisory Committee. 2018 Physical Activity Guidelines Advisory Committee Scientific Report; U.S. Department of Health and Human Services: Washington, DC, USA, 2018.
  3. Nielsen, R.O.; Bertelsen, M.L.; Ramskov, D.; Moller, M.; Hulme, A.; Theisen, D.; Finch, C.F.; Fortington, L.V.; Mansournia, M.A.; Parner, E.T. Time-to-event analysis for sports injury research part 2: Time-varying outcomes. Br. J. Sports Med. 2019, 53, 70–78. [Google Scholar] [CrossRef] [PubMed]
  4. Han, J.; Anson, J.; Waddington, G.; Adams, R.; Liu, Y. The Role of Ankle Proprioception for Balance Control in relation to Sports Performance and Injury. Biomed. Res. Int. 2015, 2015, 842804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Logerstedt, D.; Arundale, A.; Lynch, A.; Snyder-Mackler, L. A conceptual framework for a sports knee injury performance profile (SKIPP) and return to activity criteria (RTAC). Braz. J. Phys. Ther. 2015, 19, 340–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Weggemans, R.M.; Backx, F.J.G.; Borghouts, L.; Chinapaw, M.; Hopman, M.T.E.; Koster, A.; Kremers, S.; van Loon, L.J.C.; May, A.; Mosterd, A.; et al. The 2017 Dutch Physical Activity Guidelines. Int. J. Behav. Nutr. Phys. Act. 2018, 15, 58. [Google Scholar] [CrossRef]
  7. Verhagen, E.; Collard, D.; Paw, M.C.; van Mechelen, W. A prospective cohort study on physical activity and sports-related injuries in 10-12-year-old children. Br. J. Sports Med. 2009, 43, 1031–1035. [Google Scholar] [CrossRef] [Green Version]
  8. Bloemers, F.; Collard, D.; Mai, C.A.P.; Mechelen, W.V.; Twisk, J.; Verhagen, E. Physical inactivity is a risk factor for physical activity-related injuries in children. Br. J. Sports Med. 2012, 46, 669–674. [Google Scholar] [CrossRef]
  9. Zuckerman, S.L.; Wegner, A.M.; Roos, K.G.; Djoko, A.; Dompier, T.P.; Kerr, Z.Y. Injuries sustained in National Collegiate Athletic Association men’s and women’s basketball, 2009/2010-2014/2015. Br. J. Sports Med. 2018, 52, 261–268. [Google Scholar] [CrossRef]
  10. Roos, K.G.; Wasserman, E.B.; Dalton, S.L.; Gray, A.; Djoko, A.; Dompier, T.P.; Kerr, Z.Y. Epidemiology of 3825 injuries sustained in six seasons of National Collegiate Athletic Association men’s and women’s soccer (2009/2010–2014/2015). Br. J. Sports Med. 2016, 51, 1029–1034. [Google Scholar] [CrossRef]
  11. Baugh, C.M.; Weintraub, G.S.; Gregory, A.J.; Djoko, A.; Dompier, T.P.; Kerr, Z.Y. Descriptive Epidemiology of Injuries Sustained in National Collegiate Athletic Association Men’s and Women’s Volleyball, 2013–2014 to 2014–2015. Sports Health 2018, 10, 60–69. [Google Scholar] [CrossRef]
  12. Gao, Y.; Cai, W.; Gao, L.; Wang, J.; Liang, J.; Kwok, H.; Jia, C.; Li, L. Physical activity-related injuries among university students: A multicentre cross-sectional study in China. BMJ Open 2018, 8, e021845. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Collard, D.C.; Verhagen, E.A.; van Mechelen, W.; Heymans, M.W.; Chinapaw, M.J. Economic burden of physical activity-related injuries in Dutch children aged 10-12. Br. J. Sports Med. 2011, 45, 1058–1063. [Google Scholar] [CrossRef] [PubMed]
  14. Finch, C.; Cassell, E. The public health impact of injury during sport and active recreation. J. Sci. Med. Sport 2006, 9, 490–497. [Google Scholar] [CrossRef] [PubMed]
  15. Abernethy, L.; Bleakley, C. Strategies to prevent injury in adolescent sport: A systematic review. Br. J. Sports Med. 2007, 41, 627–638. [Google Scholar] [CrossRef] [PubMed]
  16. Goossens, L.; Verrelst, R.; Cardon, G.; De Clercq, D. Sports injuries in physical education teacher education students. Scand. J. Med. Sci. Sports 2014, 24, 683–691. [Google Scholar] [CrossRef]
  17. Cai, W.; Gao, Y.; Yang, W.; Cheng, F.; Tang, D.; Li, L. Physical Activity-Related Injury and Its Associated Factors among Middle School Students in Southern China. Int. J. Environ. Res. Public Health 2018, 15, 1244. [Google Scholar] [CrossRef] [Green Version]
  18. Cai, W.; Chen, S.; Li, L.; Yue, P.; Yu, X.; Gao, L.; Yang, W.; Jia, C.; Gao, Y. Gender-specific physical activity-related injuries and risk factors among university students in China: A multicentre population-based cross-sectional study. BMJ Open 2020, 10, e040865. [Google Scholar] [CrossRef]
  19. Pickett, W.; Molcho, M.; Simpson, K.; Janssen, I.; Kuntsche, E.; Mazur, J.; Harel, Y.; Boyce, W.F. Cross national study of injury and social determinants in adolescents. Inj. Prev. 2005, 11, 213–218. [Google Scholar] [CrossRef] [Green Version]
  20. Wiese-Bjornstal, D.M. Psychology and socioculture affect injury risk, response, and recovery in high-intensity athletes: A consensus statement. Scand. J. Med. Sci. Sports 2010, 20, 103–111. [Google Scholar] [CrossRef]
  21. Cai, W.; Tang, D.; Gao, Y.; Yang, W.; Chen, S.; Gao, L.; Jia, C.; Li, L. Characteristics and associated factors of physical activity-related injuries among university students in southern China. Sci. Rep. 2020, 10, 4174. [Google Scholar] [CrossRef]
  22. Perkovic, M.N.; Strac, D.S.; Tudor, L.; Konjevod, M.; Erjavec, G.N.; Pivac, N. Catechol-O-methyltransferase, Cognition and Alzheimer’s Disease. Curr. Alzheimer Res. 2018, 15, 408–419. [Google Scholar] [CrossRef]
  23. Bastos, P.; Gomes, T.; Ribeiro, L. Catechol-O-Methyltransferase (COMT): An Update on Its Role in Cancer, Neurological and Cardiovascular Diseases. Rev. Physiol. Biochem. Pharmacol. 2017, 173, 1–39. [Google Scholar] [CrossRef]
  24. Craddock, N.; Owen, M.J.; O’Donovan, M.C. The catechol-O-methyl transferase (COMT) gene as a candidate for psychiatric phenotypes: Evidence and lessons. Mol. Psychiatry 2006, 11, 446–458. [Google Scholar] [CrossRef] [Green Version]
  25. Crum, A.J.; Akinola, M.; Turnwald, B.P.; Kaptchuk, T.J.; Hall, K.T. Catechol-O-Methyltransferase moderates effect of stress mindset on affect and cognition. PLoS ONE 2018, 13, e0195883. [Google Scholar] [CrossRef] [Green Version]
  26. Smolka, M.N. Catechol-O-Methyltransferase val158met Genotype Affects Processing of Emotional Stimuli in the Amygdala and Prefrontal Cortex. J. Neurosci. 2005, 25, 836–842. [Google Scholar] [CrossRef] [Green Version]
  27. Raffaella, C.; Stefano, P.; Ina, G.; Hartmann, A.M.; Hans-Jürgen, M.; Diana, D.R.; Alessandro, S.; Dan, R. Catechol-o-methyltransferase gene modulation on suicidal behavior and personality traits: Review, meta-analysis and association study. J. Psychiatr. Res. 2011, 45, 309–321. [Google Scholar] [CrossRef]
  28. Harrison, P.J.; Tunbridge, E.M. Catechol-O-Methyltransferase (COMT): A Gene Contributing to Sex Differences in Brain Function, and to Sexual Dimorphism in the Predisposition to Psychiatric Disorders. Neuropsychopharmacology 2008, 33, 3037–3045. [Google Scholar] [CrossRef] [PubMed]
  29. van Rooij, S.J.; Stevens, J.S.; Ely, T.D.; Fani, N.; Smith, A.K.; Kerley, K.A.; Lori, A.; Ressler, K.J.; Jovanovic, T. Childhood Trauma and COMT Genotype Interact to Increase Hippocampal Activation in Resilient Individuals. Front Psychiatry 2016, 7, 156. [Google Scholar] [CrossRef] [Green Version]
  30. Haraldsson, H.M.; Ettinger, U.; Magnusdottir, B.B.; Sigmundsson, T.; Sigurdsson, E.; Ingason, A.; Petursson, H. Catechol-O-Methyltransferase Val158Met Polymorphism and Antisaccade Eye Movements in Schizophrenia. Schizophr. Bull. 2010, 36, 157–164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Chen, J.; Lipska, B.K.; Halim, N.; Ma, Q.D.; Matsumoto, M.; Melhem, S.; Kolachana, B.S.; Hyde, T.M.; Herman, M.M.; Apud, J.; et al. Functional analysis of genetic variation in catechol-O-methyltransferase (COMT): Effects on mRNA, protein, and enzyme activity in postmortem human brain. Am. J. Hum. Genet. 2004, 75, 807–821. [Google Scholar] [CrossRef] [Green Version]
  32. Hirata, H.; Hinoda, Y.; Okayama, N.; Suehiro, Y.; Kawamoto, K.; Kikuno, N.; Rabban, J.T.; Chen, L.-M.; Dahiya, R. COMT polymorphisms affecting protein expression are risk factors for endometrial cancer. Mol. Carcinog. 2008, 47, 768–774. [Google Scholar] [CrossRef] [PubMed]
  33. Diatchenko, L.; Slade, G.D.; Nackley, A.G.; Bhalang, K.; Sigurdsson, A.; Belfer, I.; Goldman, D.; Xu, K.; Shabalina, S.A.; Shagin, D.; et al. Genetic basis for individual variations in pain perception and the development of a chronic pain condition. Hum. Mol. Genet. 2005, 14, 135–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Nackley, A.G.; Shabalina, S.A.; Tchivileva, I.E.; Satterfield, K.; Korchynskyi, O.; Makarov, S.S.; Maixner, W.; Diatchenko, L. Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure. Science 2006, 314, 1930–1933. [Google Scholar] [CrossRef] [Green Version]
  35. Okochi, T.; Ikeda, M.; Kishi, T.; Kawashima, K.; Kinoshita, Y.; Kitajima, T.; Yamanouchi, Y.; Tomita, M.; Inada, T.; Ozaki, N.; et al. Meta-analysis of association between genetic variants in COMT and schizophrenia: An update. Schizophr. Res. 2009, 110, 140–148. [Google Scholar] [CrossRef] [PubMed]
  36. Chen, C.Y.; Lu, R.B.; Yeh, Y.W.; Shih, M.C.; Huang, S.Y. Association study of catechol-O-methyltransferase gene polymorphisms with schizophrenia and psychopathological symptoms in Han Chinese. Genes Brain Behav. 2011, 10, 316–324. [Google Scholar] [CrossRef] [PubMed]
  37. Cao, Y.; Wang, D.; Liu, B.; Yao, G.; Fu, Y.; Bi, Z. Catechol-O-methyltransferase (COMT) Val158Met polymorphism and risk of osteoporotic fracture. Mol. Biol. Rep. 2012, 39, 2975–2979. [Google Scholar] [CrossRef]
  38. Michaelovsky, E.; Gothelf, D.; Korostishevsky, M.; Frisch, A.; Burg, M.; Carmel, M.; Steinberg, T.; Inbar, D.; Apter, A.; Weizman, A. Association between a common haplotype in the COMT gene region and psychiatric disorders in individuals with 22q11.2DS. Int. J. Neuropsychopharmacol. 2008, 11, 351–363. [Google Scholar] [CrossRef] [Green Version]
  39. Gauderman, W.J. Sample size requirements for matched case-control studies of gene-environment interaction. Stat. Med. 2002, 21, 35–50. [Google Scholar] [CrossRef]
  40. Yan, L.H.; Chen, P.J.; Zhuang, J. Revision and reliability validity assessment of Children’s Leisure Activities Study Survey. Chin. J. Sch. Health 2011, 32, 268–270. [Google Scholar]
  41. Taylor, H.L.; Jacobs, D.R., Jr.; Schucker, B.; Knudsen, J.; Leon, A.S.; Debacker, G. A questionnaire for the assessment of leisure time physical activities. J. Chronic Dis. 1978, 31, 741–755. [Google Scholar] [CrossRef]
  42. Shieh, G. Power and sample size calculations for comparison of two regression lines with heterogeneous variances. PLoS ONE 2018, 13, e0207745. [Google Scholar] [CrossRef]
  43. Nogueira, N.; Bacelar, M.F.B.; Ferreira, B.P.; Parma, J.O.; Lage, G.M. Association between the catechol-O-methyltransferase (COMT) Val(158)Met polymorphism and motor behavior in healthy adults: A study review. Brain Res. Bull. 2019, 144, 223–232. [Google Scholar] [CrossRef] [PubMed]
  44. Rajkumar, R.P. Warriors, Worriers, and COVID-19: An Exploratory Study of the Catechol O-Methyltransferase Val158Met Polymorphism Across Populations. Cureus 2020, 12, e10103. [Google Scholar] [CrossRef] [PubMed]
  45. Roten, L.T.; Fenstad, M.H.; Forsmo, S.; Johnson, M.P.; Moses, E.K.; Austgulen, R.; Skorpen, F. A low COMT activity haplotype is associated with recurrent preeclampsia in a Norwegian population cohort (HUNT2). Mol. Hum. Reprod. 2011, 17, 439–446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Gabriel, S.B.; Schaffner, S.F.; Nguyen, H.; Moore, J.M.; Roy, J.; Blumenstiel, B.; Higgins, J.; DeFelice, M.; Lochner, A.; Faggart, M.; et al. The structure of haplotype blocks in the human genome. Science 2002, 296, 2225–2229. [Google Scholar] [CrossRef] [Green Version]
  47. Park, D.J.; Kim, S.H.; Nah, S.S.; Lee, J.H.; Kim, S.K.; Lee, Y.A.; Hong, S.J.; Kim, H.S.; Lee, H.S.; Kim, H.A.; et al. Association between catechol-O-methyl transferase gene polymorphisms and fibromyalgia in a Korean population: A case-control study. Eur. J. Pain 2016, 20, 1131–1139. [Google Scholar] [CrossRef]
  48. McAllister, T.W. Genetic factors in traumatic brain injury. Handb. Clin. Neurol. 2015, 128, 723–739. [Google Scholar] [CrossRef]
  49. Hou, T.T.; Lin, F.; Bai, S.; Cleves, M.A.; Xu, H.M.; Lou, X.Y. Generalized multifactor dimensionality reduction approaches to identification of genetic interactions underlying ordinal traits. Genet. Epidemiol. 2019, 43, 24–36. [Google Scholar] [CrossRef] [Green Version]
  50. Olesen, A.E.; Nielsen, L.M.; Feddersen, S.; Erlenwein, J.; Petzke, F.; Przemeck, M.; Christrup, L.L.; Drewes, A.M. Association Between Genetic Polymorphisms and Pain Sensitivity in Patients with Hip Osteoarthritis. Pain Pract. 2018, 18, 587–596. [Google Scholar] [CrossRef] [Green Version]
Ijerph 18 10828 g001 550
Figure 1. Flow diagram of participants included in the study.
Figure 1. Flow diagram of participants included in the study.
Ijerph 18 10828 g001
Ijerph 18 10828 g002 550
Figure 2. Linkage disequilibrium correlation coefficients of the COMT SNPs.
Figure 2. Linkage disequilibrium correlation coefficients of the COMT SNPs.
Ijerph 18 10828 g002
Table
Table 1. The differences between participants who agreed to participate in the interview and those who refused to participate.
Table 1. The differences between participants who agreed to participate in the interview and those who refused to participate.
CharacteristicsResponders (n = 1421)Non-Responders (n = 1969)χ2/tp-Value
Gender 0.9380.333
 Male377 (26.5%)552 (28.0)
 Female1 044 (73.5%)1417 (72.0)
Grade 30.890.001
 Year 1556 (39.1)620 (31.5)
 Year 2447 (31.5)605 (30.7)
 Year 3418 (29.4)744 (37.8)
Club member 3.1110.078
 Yes1065 (74.95)1527 (77.56)
 No356 (25.05)442 (22.44)
Age (Years)20.09 ± 1.1520.04 ± 1.380.5750.561
PA participation (min/wk)691.21 ± 632.64537.90 ± 702.4724.3820.001
PARI 5.3400.069
 01 132 (79.7)1517 (77.1)
 1–2227 (16.0)334 (17.0)
 ≥362 (4.4)118 (5.9)
Table
Table 2. Characteristics of physical activity-related injuries (PARIs) cases and non-PARI controls.
Table 2. Characteristics of physical activity-related injuries (PARIs) cases and non-PARI controls.
CharacteristicsControls (n = 61)Cases (n = 61)χ2/vsp-Value
Gender 0.0001.000
 Male30 (49.2)30 (49.2)
 Female31 (50.8)31 (50.8)
Grade 0.0001.000
 Year 128 (45.9)28 (45.9)
 Year 220 (32.8)20 (32.8)
 Year 313 (21.3)13 (21.3)
Club member 30.5000.001
 Yes10 (16.4)40 (65.6)
 No51 (83.6)21 (34.4)
Age (Years)19.89 ± 1.0220.02 ± 1.23−0.6410.523
BMI (kg/m2)20.24 ± 2.7721.12 ± 3.30−1.5930.114
WHR0.78 ± 0.060.79 ± 0.06−0.0860.932
Screen time (hours/day)7.16 ± 2.786.66 ± 2.810.9890.324
PA participation (min/wk)781.38 ± 575.871070.41 ± 545.52−2.8460.005
BMI, body mass index; WHR, waist-hip ratio; PA, physical activity. Figures in parentheses indicate percentages.
Table
Table 3. Allele and genotype distributions of the SNPs in the COMT gene in the PARIs cases and the non-PARI controls.
Table 3. Allele and genotype distributions of the SNPs in the COMT gene in the PARIs cases and the non-PARI controls.
SNPsControlsCasesχ2, p-ValueOR
rs9265 0.482, 0.786
CC17 (28.8)16 (27.6) 1.000 (ref.)
AC25 (42.4)28 (48.3) 0.875 (0.327–2.339)
AA17 (28.8)14 (24.1) 1.190 (0.499–2.840)
A%50.048.30.070, 0.7920.933 (0.559–1.558)
rs4680 3.607, 0.165
GG32 (56.1)23 (39.7) 1.000 (ref.)
AG24 (42.1)32 (55.2) 1.855 (0.874–3.939)
AA1 (1.8)3 (5.2) 4.174 (0.408–42.716)
A%22.832.82.835, 0.0921.649 (0.919–2.959)
rs6269 6.047, 0.049
AA20 (33.3)30 (50.8) 1.000 (ref.)
GA30 (50.0)26 (44.1) 0.578 (0.267–1.250)
GG10 (16.7)3 (5.1) 0.200 (0.049–0.818)
G%41.727.15.576, 0.0180.521 (0.302–0.898)
rs4818 5.788, 0.055
CC20 (32.8)30 (51.7) 1.000 (ref.)
GC32 (52.5)25 (43.1) 0.521 (0.241–1.126)
GG9 (14.8)3 (5.2) 0.222 (0.054–0.923)
G%41.026.75.385, 0.0200.525 (0.304–0.908)
rs4633 1.995, 0.369
CC39 (63.9)33 (54.1) 1.000 (ref.)
TC20 (32.8)23 (37.7) 1.359 (0.637–2.899)
TT2 (3.3)5 (8.2) 2.955 (0.538–16.239)
T%19.727.11.854, 0.1731.514 (0.832–2.756)
rs165655 0.081, 0.960
GG17 (27.9)18 (30.0) 1.000 (ref.)
AG34 (55.7)32 (53.3) 0.889 (0.291–2.018)
AA10 (16.4)10 (16.7) 0.944(0.315–2.834)
A%44.343.30.021, 0.8840.963 (0.579–1.600)
rs165656 1.395, 0.498
GG38 (63.3)32 (53.3) 1.000 (ref.)
CG19 (31.7)23 (38.3) 1.438 (0.667–3.099)
CC3 (5.0)5 (8.3) 1.979 (0.439–8.929)
C%20.827.51.455, 0.2281.441 (0.795–2.615)
rs165722 1.019, 0.601
CC37 (61.7)32 (55.2) 1.000 (ref.)
TC21 (35.0)22 (37.9) 1.211 (0.565–2.957)
TT2 (3.3)4 (6.9) 2.312 (0.397–13.469)
T%20.825.90.835, 0.3611.326 (0.723–2.429)
SNPs, single nucleotide polymorphisms; COMT, catechol-O-methyltransferase; A, Adenine; G, Guanine; C, Cytosine; T, Thymine. Figures in parentheses indicate percentages or 95% CIs.
Table
Table 4. Association of the dominant model and recessive model in the PARIs cases and the non-PARI controls.
Table 4. Association of the dominant model and recessive model in the PARIs cases and the non-PARI controls.
SNPsOR (95% CI) 1OR (95% CI) 2OR (95% CI) 3
rs9265
DOM1.062 (0.457–2.378)0.693 (0.267–1.800)0.719 (0.273–1.897)
REC0.786 (0.345–1.792)0.875 (0.332–2.302)0.850 (0.297–2.437)
rs4680
DOM1.948 (0.928–4.090)1.684 (0.726–3.908)1.938 (0.798–4.699)
REC3.055 (0.308–30.271)2.166 (0.168–27.946)0.967 (0.043–21.733)
rs6269
DOM0.483 (0.230–1.014)0.370 (0.150–0.914)0.355 (0.137–0.918)
REC0.268 (0.070–1.028)0.426 (0.096–1.887)0.351 (0.075–1.641)
rs4818
DOM0.455 (0.217–0.956)0.389 (0.161–0.938)0.364 (0.144–0.924)
REC0.315 (0.081–1.229)0.484 (0.107–2.182)0.385 (0.082–1.822)
rs4633
DOM1.504 (0.728–3.108)1.427 (0.611–3.334)1.523 (0.626–3.707)
REC2.634 (0.491–14.134)3.255 (0.522–20.286)2.418 (0.324–18.037)
rs165655
DOM0.902 (0.411–1.979)0.660 (0.262–1.663)0.711 (0.276–1.827)
REC1.020 (0.391–2.662)0.994 (0.314–3.144)0.900 (0.262–3.098)
rs165656
DOM1.511 (0.728–3.316)1.418 (0.598–3.362)1.532 (0.616–3.813)
REC1.727 (0.394–7.577)2.118 (0.411–10.906)1.545 (0.237–10.079)
rs165722
DOM1.307 (0.627–2.723)1.233 (0.525–2.899)1.320 (0.539–3.237)
REC2.148 (0.378–12.207)3.023 (0.465–19.631)2.214 (0.291–16.882)
1: Adjusted for none covariate; 2: Adjusted for gender, grade, and sports team membership; 3: Adjusted for gender, grade, sports team membership, age, body mass index, waist-hip ratio, screen time, and PA participation; SNPs, single nucleotide polymorphisms; DOM, dominant model; REC, recessive model.
Table
Table 5. Haplotype frequencies of the SNPs in the COMT gene in the PARIs cases and the non-PARI controls.
Table 5. Haplotype frequencies of the SNPs in the COMT gene in the PARIs cases and the non-PARI controls.
BlockHaplotypeControlsCasesχ2p-ValueOR (95% CI)
1 a
C-A-C-C0.3850.4521.1370.2861.490 (0.876–2.533)
C-G-C-G0.3930.2684.3330.0370.474 (0.269–0.834)
T-A-T-C0.1970.2701.8540.1731.466 (0.786–2.735)
2 b
C-G0.4990.5220.1310.7171.111 (0.664–1.860)
A-A0.4430.4410.0010.9800.992 (0.590–1.666)
A-G0.0590.0370.6210.4310.577 (0.164–2.026)
a: The bases of block one are listed in the following order: rs165722, rs6269, rs4633, and rs4818; b: The bases of block two are listed in the following order: rs9265 and rs165655; SNPs, single nucleotide polymorphisms; COMT, catechol-O-methyltransferase.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, S.; Cai, W.; Duan, S.; Gao, L.; Yang, W.; Gao, Y.; Jia, C.; Zhang, H.; Li, L. Association of COMT Polymorphisms with Multiple Physical Activity-Related Injuries among University Students in China. Int. J. Environ. Res. Public Health 2021, 18, 10828. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph182010828

AMA Style

Chen S, Cai W, Duan S, Gao L, Yang W, Gao Y, Jia C, Zhang H, Li L. Association of COMT Polymorphisms with Multiple Physical Activity-Related Injuries among University Students in China. International Journal of Environmental Research and Public Health. 2021; 18(20):10828. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph182010828

Chicago/Turabian Style

Chen, Shangmin, Weicong Cai, Shiwei Duan, Lijie Gao, Wenda Yang, Yang Gao, Cunxian Jia, Hongjuan Zhang, and Liping Li. 2021. "Association of COMT Polymorphisms with Multiple Physical Activity-Related Injuries among University Students in China" International Journal of Environmental Research and Public Health 18, no. 20: 10828. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph182010828

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