Logo Medical Science Monitor

Call: +1.631.470.9640
Mon - Fri 10:00 am - 02:00 pm EST

Contact Us

Logo Medical Science Monitor Logo Medical Science Monitor Logo Medical Science Monitor

16 March 2021: Clinical Research  

Associations Between Cytochrome P450 (CYP) Gene Single-Nucleotide Polymorphisms and Second-to-Fourth Digit Ratio in Chinese University Students

Jing Zhang12ABCDEF, Mengyi Yang1B, Pengfei Luan1B, Wei Jia1B, Qiujun Liu1B, Zhanbing Ma1A, Jie Dang1A, Hong Lu1A, Qian Ma1B, Yanfeng Wang1B, Chunlan Mu1B, Zhenghao Huo13AG*

DOI: 10.12659/MSM.930591

Med Sci Monit 2021; 27:e930591

0 Comments

Abstract

BACKGROUND: Cytochrome P450 (CYP) genes are necessary for the production or metabolism of fetal sex hormones during pregnancy. The second-to-fourth digit ratio (2D: 4D) is formed in the early stage of human fetal development and considered an indicator reflecting prenatal sex steroids levels. We explored the association between 2D: 4D and single-nucleotide polymorphisms (SNPs) of CYP.

MATERIAL AND METHODS: Correlation analysis between 2D: 4D and 8 SNPs, rs2687133 (CPY3A7), rs7173655 (CYP11A1), rs1004467, rs17115149, and rs2486758 (CYP17A1), and rs4646, rs2255192, rs4275794 (CYP19A1), was performed using data from 426 female and 412 male Chinese university students. SNP genotyping was conducted using PCR. Digit lengths were photographed and measured by image processing software.

RESULTS: rs2486758 (CYP17A1) correlated with left hand 2D: 4D in men (P=0.026), and rs1004467 (CYP17A1) correlated with right hand 2D: 4D in men (P=0.008) and the whole population (P=0.032). In men, allele G rs1004467 decreased right hand 2D: 4D, while allele C of rs2486758 increased left hand 2D: 4D. In women, left hand 2D: 4D was higher in genotypes with allele A of SNP rs4646 (CYP19A1) under the dominant genetic model; female DR-L was higher in genotypes with allele T of rs17115149 (CYP11A1). SNPs rs2687133 (CYP3A7) and rs1004467 (CYP17A1) were significantly correlated with right hand 2D: 4D (P=0.0107).

CONCLUSIONS: SNPs rs1004467 and rs2486758 of CYP17A1 are significant in the relationship between 2D: 4D and CYP gene polymorphisms under different conditions. SNP interactions between CYP genes probably impact 2D: 4D. The correlation between 2D: 4D and some sex hormone-related diseases may be due to the effect of CYP variants on the 2 phenotypes.

Keywords: Fetal Development, Genetics, Gonadal Steroid Hormones, Polymorphism, Single Nucleotide, Alleles, aromatase, Asians, Case-Control Studies, Cholesterol Side-Chain Cleavage Enzyme, Cytochrome P-450 CYP3A, Cytochrome P-450 Enzyme System, Fingers, Gene Frequency, Genetic Association Studies, Genotype, Nucleotides, Steroid 17-alpha-Hydroxylase, Students, Universities, young adult

Background

The ratio of the second-to-fourth digit length (2D: 4D) in human hands is a sex dimorphism trait that tends to be higher in females than in males and is believed to reflect prenatal sex steroid levels [1]. This ratio is possibly determined by the balance of prenatal testosterone and estrogen levels in a narrow window during the fetal period and remains stable after birth [2]. Thus, 2D: 4D is commonly used as a noninvasive retrospective index of prenatal exposure to sex steroids, especially androgen and estrogen. Studies have shown that 2D: 4D is strongly associated with health conditions and diseases of sex dimorphisms or those affected by sex hormones. Therefore, 2D: 4D might be an indicator of human conditions such as body composition [3], personality traits [4], and physical performance [5], and of disease susceptibilities, including cancer [6], congenital adrenal hyperplasia [7], infertility [8], and mental disorders [9].

Genetic studies on 2D: 4D in twins have shown that 2D: 4D is highly heritable (about 80%) [10]. Gene polymorphisms and the 2D: 4D ratio are significant in studies on factors that can affect the upper limbs [11]. Therefore, 3 genome-wide association studies (GWAS) have been conducted to identify the single-nucleotide polymorphisms (SNPs) associated with 2D: 4D. In 2010, Medland et al studied a sample of 1382 twins and their siblings from Australia, and the SNP rs314277 of gene LIN28B was found to be strongly associated with 2D: 4D [12]. In 2013, Lawrance-Owen et al conducted a 2D: 4D GWAS in the UK and found that the SMOC1 polymorphism was correlated with 2D: 4D [13]. In 2018, Warrington et al performed the largest 2D: 4D genome-wide association meta-analysis, in which 11 loci were identified (GLIS1, rs4927012; EFNA1, rs11581730; LDAH, rs340600; OLA1, rs12474669; FLI1, rs10790969; HOXD11/12, rs847158; GLI3, rs77640775; SMOC1, rs2332175; SALL1, rs6499762; TOX3, rs1080014; and SALL3, rs4799176), and were believed to explain only 3.8% of the 2D: 4D discrepancy [14]. It is worth noting that none of these SNPs were located in the sex chromosome, and an Asian population was not involved in the studies. Later, other 2D: 4D correlation studies targeting hormone-related gene polymorphisms were conducted, with the majority mainly focusing on estrogen and androgen receptor genes [15,16]. However, 2D: 4D-associated SNPs of enzyme-related genes involved in hormone synthesizing and metabolism have been less frequently studied.

In humans, cytochrome P450s (CYPs) constitute one of the most abundant and diverse enzyme superfamilies, whose members catalyze the oxidative transformation of many organic substrates [17]. The function of CYP genes is critical for many metabolic processes, particularly for sex steroid conversion [18]. The altered expression of CYP genes involved in steroid metabolism may affect sex hormone levels and further increase the risk of some sex-related diseases, such as breast cancer [19], prostate cancer [20], infertility [21], polycystic ovary syndrome [22], and congenital adrenal hyperplasia [23]. During pregnancy, fetal sex hormones affect embryo development in many ways. The expression of specific CYP genes in the fetal liver, gonad, and placenta is indispensable across the development to synthesize sex hormones including progesterone, estrogen, and androgen [24]. First, CYP11A1, an enzyme initiating the steroidogenesis to converse cholesterol to pregnenolone, is pivotal for fetal survival in the early phase of pregnancy. Also, about 30% of the estrogen produced during pregnancy is yielded from the pregnenolone conversion by placental CYP17A1. In humans, dehydroepiandrosterone (DHEA) and DHEA sulfate produced by the fetal adrenal glands mainly undergo 16-hydroxylation by CYP3A7 in the liver, and are then metabolized by the placental enzyme aromatase (CYP19A1) [25]. Additionally, CYP19A1 is a rate-limiting enzyme involved in catalyzing the conversion of androgen to estrogen [26]. There is growing evidence to support that the analysis of sex hormone balance in a developing fetus might be helpful for predicting adult diseases during the lifespan [27].

In particular, some SNPs of CYP genes may affect the sex steroids, leading to physiological or pathological changes in humans. The C allele of rs45446698 in CYP3A7 causes the overexpression of the fetal CYP3A7 gene in adults and further impacts the levels of circulating endogenous sex hormones [28]. The SNP rs17115149 (C>A) of the CYP17A1 gene was strongly correlated with testosterone levels in men with infertility [29]. The SNPs in the CYP11A1 and CYP17 genes are significantly associated with testosterone plasma concentrations of patients with polycystic ovary syndrome, and SNPs rs4441215 and rs936306 of the CYP19A1 gene are associated with serum estradiol levels [30]. Further, SNPs in the CYP genes CYP3A7, CYP11A1, CYP17A1, and CYP19A1 have been shown to be associated with malignancy and androgenic enzyme expression in the fetal liver [31,32]. Therefore, it is proposed that CYP SNPs may have effects on prenatal sex hormones.

Herein, we conducted a study on 4 CYP genes involved in fetal or placental sex steroid synthesis and metabolism, namely CYP3A7, CYP11A1, CYP17A1, and CYP19A1. The association between these CYP gene polymorphisms and 2D: 4D were evaluated in a population of Chinese university students. Our study might elucidate the underlying mechanism of genetic influence on the formation of 2D: 4D. Also, the study might indirectly support 2D: 4D as a useful predictor for many hormone-related diseases, especially for those with significant sex differences in prevalence, incidence, and clinical manifestation.

Material and Methods

STUDY SUBJECTS:

The participants were 426 women aged 19.6±1.3 years and 412 men aged 19.2±1.2 years. All participants were recruited from the Ningxia Medical University in Yinchuan city of the Ningxia Hui Autonomous Region in northwestern China. All participants were physically and mentally healthy. The study design was approved by the Ningxia Medical University Ethics Committee (No.2019-057). Informed written consent was obtained from each participant. All the study samples are stored securely. The DNA samples, research data, and any information that can be traced to individual participants will not be provided to any other party.

SNP GENOTYPING:

Four CYP genes were chosen based on the review of related literature, namely CYP3A7 (critical for sex steroid synthesis during embryogenesis and expressed mainly in fetal liver), and CYP11A1, CYP17A1, and CYP19A1 (the genes of 3 rate-limiting enzymes catalyzing sequential steps in steroidogenesis). SNPs for the CYP genes were selected using the 1000 Genomes Browser database and Haploview 4.2 software, per the following criteria: (1) Functional SNPs from regions including 1000bp upstream and downstream of target genes were systematically screened using the online tools ENSEMBL (https://asia.ensembl.org) and SNPinfo Web Server (https://snpinfo.niehs.nih.gov/); (2) an r2 cut-off (0.8) and minor allele frequency cut-off (0.05) for the Han Chinese in Beijing population of the SNPs were defined; (3) the SNPs have been studied or reported in other diseases. Finally, we selected 8 SNPs, CPY3A7 (rs2687133 C>T), CYP11A1 (rs7173655 C>T), CYP17A1 (rs1004467 A>G, rs17115149 C>T, rs2486758 T>C), and CYP19A1 (rs4646 C>A, rs2255192 C>T, rs4275794 T>C).

Genomic DNA was extracted from peripheral blood samples using DNA extraction kits, according to the manufacturer’s instructions (DP304 TIANamp Genomic DNA Kit, TIANGEN Biotech Co. Ltd., Beijing, China). Genotyping of each variant was performed using the multiplex polymerase chain reaction (PCR) technique developed by Novogene Biotech Co., Ltd. (Beijing, China). The PCR program was performed in 2 rounds: round 1 consisted of 1 cycle at 95°C for 15 min; 4 cycles at 94°C for 30 s, 60°C for 10 min, and 72°C for 30 s; and 24 cycles at 94°C for 30s, 60°C for 1 min, and 72°C for 30 s; and round 2 consisted of 1 cycle at 95°C for 15 min; 5 cycles at 94°C for 30 s, 60°C for 4 min, and 72°C for 30 s; and 10 cycles at 94°C for 30 s, 65°C for 1 min, and 72°C for 30 s. After PCR, the amplified products were mixed in a centrifuge tube and shaken overnight. The PCR products were then purified using a CA2 absorption column. The purified PCR product was loaded and run on the Illumina X-10 sequencer. Illumina RTA software was used to perform quality control and raw data analysis.

DATA COLLECTION:

The participants put their 2 hands flat against the same white background with the palms facing up and fingers straight. A digital camera was used to collect the images of both hands. The photos were not used when the fine structures of each digit were not clearly visible. The Image-Pro Plus 6.0 (Media Cybernetics, Inc., USA) was used to measure the digit lengths of both hands (from the fingertip to the middle point of the most proximal crease to the palm). The measurement of each digit was performed 3 times by 3 different individuals who were well trained in the procedure, and the mean value was used for further analysis. The length ratios of both hands between the second and fourth digits were calculated.

STATISTICAL ANALYSES:

The Statistical Package for the Social Sciences (SPSS) 22.0 software (IBM, Armonk, NY, USA) was used to conduct the statistical analyses. The intra-class correlation coefficients (r1) were used for evaluating the digit length measurement credibility. The difference comparisons were performed using the t test for 2 groups and ANOVA for more than 2 groups. The Benjamini-Hochberg method was applied for multiple comparisons to control the false discovery rate. The χ2 test and Fisher’s exact test were used to evaluate the genotype and allele frequency distributions between groups. The statistical tests were 2-tailed and P<0.05 indicated statistical significance. The mean±standard deviation was used to represent all 2D: 4D data.

The Hardy-Weinberg equilibrium, minor allele frequency, and linkage disequilibrium were analyzed using Haploview software (Broad Institute of Harvard and MIT, MA, USA). The haplotype block was defined using the method of the Four Gamete Rule in Haploview.

Snpstats (https://www.snpstats.net/start.htm) was used to analyze the differences in haplotype, genotype, and allele frequencies between sexes. The linear regression analysis in the Plink package (v.1.9) was used to verify the associations between 2D: 4D and CYP SNPs. We assessed the interactions among SNPs of CYP genes using the generalized multifactor dimensionality reduction (GMDR) method [33]. The box-shaped scatter composite figure plotted by R was used to visualize the 2D: 4D genotype association characteristics.

Results

DESCRIPTION AND DISTRIBUTION OF CYP SNPS IN OUR POPULATION:

A total of 8 SNPs separately located in CPY3A7, CYP11A1, CYP17A1, and CYP19A1 were selected. The information and locations of the associated SNPs are listed in Table 1. The genotype frequencies of all the SNPs met the requirement of Hardy-Weinberg equilibrium in our population (P>0.05). The results indicated that the genotype and allele frequencies of all SNPs were not significantly different between men and women. A good linkage disequilibrium was found in the 2 groups of the following SNPs: 3 SNPs (rs1004467, rs17115149, rs2486758) of CYP17A1 and 3 SNPs (rs4646, rs2255192, rs4275794) of CYP19A1 (Figure 1). The haplotypes in the blocks were established and analyzed for differences between sexes, but no statistically significant differences were observed.

2D: 4D CHARACTERIZATIONS IN OUR POPULATION:

Remeasurement reliability of all digit lengths were acceptable for both hands in the intra-class correlation coefficient test (digit length of left hand: F=44.620, P<0.001; digit length of right hand: F=44.484, P<0.001). The significant difference was observed between different sexes in the left hand (L)2D: 4D (t=−2.951, P=0.003) and the right hand (R)2D: 4D) (t=−2.353, P=0.019). The L2D: 4D and R2D: 4D were both higher in women than in men, which was consistent with previous reports. However, the R2D: 4D minus L2D: 4D (DR-L) was not different between men and women (Table 2).

CORRELATION OF SINGLE CYP SNPS TO 2D: 4D IN THE STUDY POPULATION:

The L2D: 4D, R2D: 4D, and DR-L of CYP SNP genotypes were analyzed to explore the relationship between the CYP gene polymorphisms and the prenatal sex hormones. From the analytical results, the R2D: 4D was found to be significantly different in the 3 genotypes of SNP rs1004467 (F=4.113, P=0.017) (Table 3). The existence of minor allele G decreased the R2D: 4D, suggesting an increased fetal androgen or decreased fetal estrogen level. However, no genotype difference was observed in L2D: 4D and DR-L.

The associations between CYP SNPs and 2D: 4D by sex were then assessed, and only SNPs of CYP17A1 showed a relationship (Figure 2). For men, the L2D: 4D showed a statistically significant difference among different genotypes of SNP rs2486758 (F=4.771, P=0.009) and the genotypes with minor allele C had a higher L2D: 4D. The R2D: 4D was significantly different among different genotypes of SNPs rs1004467 (F=5.997, P=0.003) and rs2486758 (F=4.973, P=0.007). The minor allele G of rs1004467 was associated with a decrease in the R2D: 4D, while the minor allele C of rs2486758 was associated with an increase in the R2D: 4D. For women, a statistically significant difference in the L2D: 4D was observed in SNP rs4646 (F=3.142, P=0.044), and the L2D: 4D was higher in genotypes with minor allele A. Also, the DR-L in women was higher in genotypes containing the minor allele T of rs17115149 (F=3.260, P=0.039).

Further, the linear regression model was used to verify the causal relationship between CYP SNPs and 2D: 4D. As shown in Table 4, rs2486758 was correlated with the L2D: 4D (P=0.026) in men, and rs1004467 was correlated with the R2D: 4D in men (P=0.008) and the whole population (P=0.032). However, no relationship was found between CYP SNPs and the 2D: 4D in women. Additionally, there was no correlation between CYP SNPs and DR-L for either men or women.

CORRELATION OF CYP SNPS WITH 2D: 4D UNDER DIFFERENT GENETIC MODELS AND HAPLOTYPES:

The digit ratios of dominant and recessive models of SNP genotypes were compared. For the whole population, L2D: 4D was significantly lower in C/C than in A/C and A/A genotypes under the dominant model of rs4646 (t=−2.218, P=0.027), R2D: 4D was significantly higher in A/A than in A/G and G/G genotypes under the dominant model of rs1004467 (t=2.618, P=0.009), and DR-L was significantly different between the genotypes of SNP rs2687133 under the dominant model (t=2.067, P=0.039). In men, there were significant differences in L2D: 4D in the dominant model of rs2486758 (t=2.831, P=0.005), recessive models of rs1004467 (t=2.109, P=0.036) and rs2486758(t=−2.028, P=0.043), in R2D: 4D in dominant models of rs1004467 (t=2.878, P=0.004) and rs2486758 (t=−3.132, P=0.002), and the recessive model of rs1004467 (t=2.644, P=0.009). In women, the optimal models were the recessive model of rs17115149 (t=2.128, P=0.034) and the dominant model of rs4646 (t=2.419, P=0.016) for L2D: 4D, but there were no statistically significant genetic models for R2D: 4D and DR-L.

Next, 2D: 4D in different haplotypes of CYP17A1 and CYP19A1 were analyzed. As shown by the data summarized in Table 5, only rs17115149-rs1004467-rs2486758 in CYP17A1 showed significant differences in male R2D: 4D values (F=4.382, P=0.005). Multiple comparisons were performed, and the results indicated that R2D: 4D was significantly lower in the haplotype G-G-T.

IMPLICATIONS OF SNP-SNP INTERACTIONS WITHIN CYP GENES IN 2D: 4D:

Next, we used GMDR analysis to detect the potential high-dimension interactions among the CYP SNPs, with the interaction confirmed by cross-validation consistency. Key variables were defined as L2D: 4D, R2D: 4D, and DR-L for each participant, with SNP variables ranging from 0 to 2, indicating the number of minor alleles in an individual participant (0=0 minor allele, 1=1 minor allele, 2=2 minor alleles). Considering the effect on the R2D: 4D, CYP17A1 SNP rs1004667 had the highest testing balanced accuracy and CV consistency among the 8 SNPs, but without statistical significance. On testing, SNPs rs2687133 (CYP3A7) and rs1004467 (CYP17A1) had the highest accuracy in a 2-way interaction (P=0.0107) (Table 6, Figure 3).

Discussion

The current study was designed to help us understand the genetic polymorphisms of CYP genes and the association of these polymorphisms with the prenatal sex hormone balance in a northwestern Chinese population. Our results support the hypothesis that polymorphisms in CYP genes are correlated with the 2D: 4D; that is, variants of CYP genes might affect the synthesizing and metabolizing of sex hormones in utero and further result in the difference of digit ratio formation during early embryo development.

Genetic polymorphisms in CYP genes might lead to changes in CYP gene expression and thereby affect steroid synthetization and endocrine regulation [18]. Studies have found that the allele frequencies of many CYP SNPs vary between different ethnic groups [34] and diseases [35], but differences due to sex have been rarely investigated. For example, the CYP17 34 T>C allele distribution was found to be different between male (n=756) and female (n=915) healthy control participants in a study in an Austrian population, in which men had a statistically significant higher mutant C allele frequency than did women [36]. Our present study analyzed the genotype and allele frequencies of 8 SNPs of CYP genes in a population of northwestern China. However, the results indicated that there were no statistically significant differences in the genotype and allele frequencies of any SNPs between the men and women included in the study.

Research in the last decades has shown that unbalanced exposure to prenatal testosterone (PT) or prenatal estrogen (PE) may cause differences in 2D: 4D. High PT or low PE might lead to a low 2D: 4D, while low PT or low PE might result in a high 2D: 4D [37]. Many studies have demonstrated that women have a higher 2D: 4D than men, which is consistent with our findings. In addition, the DR-L or DL-R of 2D: 4D can reveal 2D: 4D asymmetry between the left and right hands [38]. Several reports have suggested that PT influences R2D: 4D and DR-L more than it influences L2D: 4D, and the DR-L of 2D: 4D is negatively correlated with prenatal exposure and sensitivity to testosterone [39]. However, in our present study, we did not observe marked differences in the range of DL-R between the sexes. Also, a study found that the number of TA(n) repeats in the promoter region of the ESR1 gene was significantly correlated with 2D: 4D in male left hands [40], but rs9340799 of the ESR1 gene might be related to the 2D: 4D of the right hands in school-aged boys [41]. Our research showed that some CYP gene polymorphisms were associated with 2D: 4D in the left or right hand, while others were associated with both left and right hands, indicating that the association between CYP gene polymorphisms and 2D: 4D had no left or right hand preference. Although the underlying mechanisms for this phenomenon are unclear, we assume that some SNP variants probably affect the CYP gene functions weakly and result in the weak effect on the associations between CYP SNPs and 2D: 4D.

CYP3A7, the most abundant CYP3A isoform in the fetus, is not typically expressed in adults. The SNP rs2687133 is an intron variant of CYP3A7. In the present study, the DR-L of T/C and C/C was higher than that of T/T genotypes in rs2687133, indicating that the minor allele C was possibly related to the prenatal sex hormone level. The SNP rs2687133 has been studied in research on drug metabolizing enzyme function and congenital adrenal hyperplasia, but no associations were observed [42]. Additional research may be required to confirm the relationship between rs2687133 and other sex steroid-related phenotypes and diseases. CYP11A1, the initial rate-limiting enzyme for steroidogenesis, is expressed mainly in the adrenal glands and gonads and is involved in the hypothalamus-pituitary-adrenal regulatory axis [43]. Because of its physiological importance, many studies have focused on the single-nucleotide polymorphism of CYP11A1. Women carrying the variant allele for SNPs rs7173655 have a higher risk of endometrial cancer than do women with 2 wild-type alleles in CYP11A1 [44]. Further, rs7173655 in CYP11A1 showed a statistically significant association with acrylamide intake for breast cancer risk [45]. Although 2D: 4D has been reported to be linked with cancers, especially breast cancer [46], no significant relationship between rs7173655 and 2D: 4D was observed in our study. We believe that the correlations between other SNPs in CYP11A1 and 2D: 4D are worthy of further study.

Three SNPs, namely rs1004467, rs17115149, and rs2486758 of CYP17A1, a gene encoding a key cytochrome P450 enzyme in the steroidogenic pathway, were analyzed to explore their roles in digit ratio formations. Our present data indicated that the R2D: 4D differed markedly in different genotypes and genetic models (dominant and recessive) of SNP rs1004467. In addition, we verified an association between rs1004467 and R2D: 4D by generalized linear analysis in men and the entire study population. The existence of the G allele significantly decreased the R2D: 4D, suggesting an increased PT or decreased PE. As previously reported, men carrying the GG genotype for SNP rs1004467 have significantly elevated concentrations of dihydrotestosterone and testosterone, which probably leads to a higher susceptibility for the development of prostate cancer [47]. Interestingly, it was found that low 2D: 4D is strongly associated with prostate cancer [6]. Additionally, the rs1004467 minor G allele of rs1004467 is statistically significantly correlated with a decreased risk of hypertension [48], and 2D: 4D is related to blood pressure. In terms of SNP rs2486758, which is mapped to the intergenic section near the 5′ UTR of the CYP17A1 gene, L2D: 4D and R2D: 4D showed significant correlations with SNP rs2486758 in men under different statistical or genetic models, and the minor allele C of rs2486758 increased the 2D: 4D values. The C allele of rs2486758 has been reported to be linked with a high serum estrogen level, prostate cancer, and coronary artery disease [49], which are also related to 2D: 4D [6,50]. In addition, rs17115149 showed a significant association with the development of prostate cancer in Korean men [51] and with testosterone levels of men with infertility [29]. In our present study, the DR-L in women was higher in minor allele T genotypes of rs17115149, and a significant difference of L2D: 4D was observed between TT and TC/CC in women. However, we did not find any association in men and failed to verify an association by linear analysis, indicating that a larger sample size and a more diverse ethnic population might be needed. It has been implied that rs1004467 and rs2486578 of CYP17A1 are important SNP loci for explaining digit ratio development and elucidating the fetal origin and organizational effect of disorders including prostate cancer, hypertension, and coronary artery disease, especially for men.

The product of the CYP19A1 gene, aromatase, is the key enzyme participating in converting androgen to estrogen. In human placental tissues, expression levels of CYP19A1 are elevated dramatically during pregnancy. A previous study on the relationship between 6 SNPs of CYP19A1 and 2D: 4D indicated that an association between rs4775936 and DR-L was found only in women [52]. In our study, we analyzed 3 SNPs which were also of clinical significance but we found no relationship between 2D: 4D and rs2255192 or rs4275794 of CYP19A1. However, we found rs4646 was correlated with L2D: 4D in women, under the dominant model. The L2D: 4D of genotype A/C and A/A was higher than that of genotype C/C, suggesting the A allele was associated with low PT or high PE. Researchers discovered that the rs4646 C/C genotype is correlated with higher estrogen levels in patients with breast cancer [53]. The homozygous minor allele (AA) of rs4646 is significantly related to the improved therapy outcomes of patients with early breast cancer [54]. It was found that a high 2D: 4D is linked with elevated breast cancer risk and bad prognosis [6]. Also, rs4775936 has been linked to both 2D: 4D and breast cancer in the study by Zhanbing et al [52]. For this reason, SNPs of CYP19A1 are valuable for identifying the 2D: 4D variances in women.

Heredity and GWAS studies have demonstrated that digit ratio is a polygenic genetic trait; that is, a single gene seldomly contributes to its development. The formation of 2D: 4D is influenced by multiple genes and environmental factors. Hence, the combined function and interaction among CYP genes involved in the synthesizing and metabolizing of sex steroids were studied to elucidate the relationship between CYP SNPs and 2D: 4D. The haplotype of rs17115149-rs1004467-rs2486758 in CYP17A1 showed remarkable differences compared with other haplotypes in R2D: 4D values in men, suggesting that multiple SNPs of CYP17A1 and the effect of their interaction might alter the digit ratios more obviously. As we know, progesterone is a common substrate for both CYP3A7 and CYP17A1 in producing downstream sex steroids [55]. It was found that the SNPs rs2687133 (CYP3A7) and rs1004467 (CYP17A1) have a strong interaction with R2D: 4D, which better explains the formation of 2D: 4D under the SNP-SNP interaction model between CYP genes.

Our study has the following limitations: 1) the SNPs of CYP genes selected in our research were not numerous enough to cover all potential function locations, and therefore more SNPs, as well as their interactions, might be required in future studies; 2) the source of all participants was limited to the Ningxia Medical University, and the small sample size of this study probably resulted in weak statistical significance; 3) the factor of a limited ethnic group might affect the results of SNP distributions and digit ratios; and 4) the study lacks functional validation. Therefore, an association study with more CYP SNP loci, a larger sample size, more ethic groups with different genetic backgrounds, and molecular function validation might be conducted in future.

Conclusions

Our research demonstrated that the SNPs rs1004467 and rs2486758 of the CYP17A1 gene are of great significance in the evaluation of the relationship between 2D: 4D and CYP gene polymorphisms under different conditions. SNP interactions between CYP genes also probably impact the association with 2D: 4D. Further, our study results suggest that the correlation between 2D: 4D and some sex hormone-related diseases is possibly due to the common effect of CYP variants on the 2 different phenotypes.

Figures

Linkage disequilibrium (LD) map of single-nucleotide polymorphisms (SNPs) in genes CYP17A1 and CYP19A1. A good LD was found among rs1004467, rs2486758, and rs17115149 of CYP17A1(rs1004467 – rs2486758: D’=0.9993, r2=0.1153; rs1004467 – rs17115149: D’=0.6885, r2=0.0202; rs2486758 – rs17115149: D’=0.9964, r2=0.0203); and rs2255192, rs4275794, and rs4646 of CYP19A1(rs2255192 – rs4275794: D’=0.9883, r2=0.9728; rs2255192 – rs4646: D’=0.9992, r2=0.0974; rs4275794 – rs4646: D’=0.9992, r2=0.0978). The haplotypes in the blocks were established based on the LD values.Figure 1. Linkage disequilibrium (LD) map of single-nucleotide polymorphisms (SNPs) in genes CYP17A1 and CYP19A1. A good LD was found among rs1004467, rs2486758, and rs17115149 of CYP17A1(rs1004467 – rs2486758: D’=0.9993, r2=0.1153; rs1004467 – rs17115149: D’=0.6885, r2=0.0202; rs2486758 – rs17115149: D’=0.9964, r2=0.0203); and rs2255192, rs4275794, and rs4646 of CYP19A1(rs2255192 – rs4275794: D’=0.9883, r2=0.9728; rs2255192 – rs4646: D’=0.9992, r2=0.0974; rs4275794 – rs4646: D’=0.9992, r2=0.0978). The haplotypes in the blocks were established based on the LD values. Second-to-fourth digit ratios (2D: 4D) of different sexes distributed in different genotypes of the CYP17A1 single-nucleotide polymorphisms (SNPs). (A–C) Left hand 2D: 4D of rs17115149, rs1004467, and rs2486758 genotypes; (D–F) Right hand 2D: 4D of rs17115149, rs1004467, and rs2486758 genotypes. The black scatter points are individual 2D: 4D observations (including 400 men and 404 women). The box along with the upper and lower vertical lines represent the quartile range of 2D: 4D. Different colors of boxes represent different phenotypes. The horizontal line inside the box represents the mean value of the 2D: 4D. Remeasurement reliability of all digit lengths among 3 operators were acceptable for both hands in the intra-class correlation coefficient test.Figure 2. Second-to-fourth digit ratios (2D: 4D) of different sexes distributed in different genotypes of the CYP17A1 single-nucleotide polymorphisms (SNPs). (A–C) Left hand 2D: 4D of rs17115149, rs1004467, and rs2486758 genotypes; (D–F) Right hand 2D: 4D of rs17115149, rs1004467, and rs2486758 genotypes. The black scatter points are individual 2D: 4D observations (including 400 men and 404 women). The box along with the upper and lower vertical lines represent the quartile range of 2D: 4D. Different colors of boxes represent different phenotypes. The horizontal line inside the box represents the mean value of the 2D: 4D. Remeasurement reliability of all digit lengths among 3 operators were acceptable for both hands in the intra-class correlation coefficient test. Single-nucleotide polymorphism (SNP) interaction graphs analyzed by generalized multifactor dimensionality reduction for second-to-fourth digit ratio (2D: 4D). (A, B) Optimal models as determined by generalized multifactor dimensionality reduction (GMDR) for CYP SNPs (0 = 0 minor allele; 1=1 minor allele; 2=2 minor alleles). The numbers within each small square represent positive scores (left) and negative scores (right). For each square, dark blu cells indicate the high effect on the 2D: 4D, and light blue cells represent the low effect on the 2D: 4D. The test was adjusted for sex.Figure 3. Single-nucleotide polymorphism (SNP) interaction graphs analyzed by generalized multifactor dimensionality reduction for second-to-fourth digit ratio (2D: 4D). (A, B) Optimal models as determined by generalized multifactor dimensionality reduction (GMDR) for CYP SNPs (0 = 0 minor allele; 1=1 minor allele; 2=2 minor alleles). The numbers within each small square represent positive scores (left) and negative scores (right). For each square, dark blu cells indicate the high effect on the 2D: 4D, and light blue cells represent the low effect on the 2D: 4D. The test was adjusted for sex.

References

1. Manning J, Cook C, Crewther B, Digit ratio (2D:4D) and testosterone supplementation: Early Hum Dev, 2019; 139; 104843

2. Zheng Z, Cohn MJ, Developmental basis of sexually dimorphic digit ratios: Proc Natl Acad Sci USA, 2011; 108(39); 16289-94

3. Wu L, Yao R, Zhang Y, The association between digit ratio (2D:4D) and overweight or obesity among Chinese children and adolescents:A cross-sectional study: Early Hum Dev, 2019; 136; 14-20

4. Nitschke JP, Bartz JA, Lower digit ratio and higher endogenous testosterone are associated with lower empathic accuracy: Horm Behav, 2020; 119; 104648

5. Eklund E, Ekström L, Thörngren JO, Ericsson M, Digit ratio (2D:4D) and physical performance in female olympic athletes: Front Endocrinol (Lausanne), 2020; 11; 292

6. Bunevicius A, The association of digit ratio (2D:4D) with cancer:A systematic review and meta-analysis: Dis Markers, 2018; 2018; 7698193

7. Richards G, Browne WV, Aydin E, Digit ratio (2D:4D) and congenital adrenal hyperplasia (CAH):Systematic literature review and meta-analysis: Horm Behav, 2020; 126; 104867

8. Vélez MP, Arbuckle TE, Monnier P, Fraser WD, Female digit length ratio (2D:4D) and time-to-pregnancy: Hum Reprod, 2016; 31(9); 2128-34

9. Tegin C, Canan F, El-Mallakh RS: J Affect Disord, 2019; 259; 27-30

10. Myers L, Van’t Westeinde A, Kuja-Halkola R, 2D:4D Ratio in neurodevelopmental disorders:A twin study: J Autism Dev Disord, 2018; 48(9); 3244-52

11. Slama N, Warner M, Mocarelli P: Early Hum Dev, 2019; 131; 45-50

12. Medland SE, Zayats T, Glaser B, A variant in LIN28B is associated with 2D:4D finger-length ratio, a putative retrospective biomarker of prenatal testosterone exposure: Am J Hum Genet, 2010; 86(4); 519-25

13. Lawrance-Owen AJ, Bargary G, Bosten JM, Genetic association suggests that SMOC1 mediates between prenatal sex hormones and digit ratio: Hum Genet, 2013; 132(4); 415-21

14. Warrington NM, Shevroja E, Hemani G, Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero: Hum Mol Genet, 2018; 27(11); 2025-38

15. Zhang K, Yang X, Zhang M, Revisiting the relationships of 2D:4D with androgen receptor (AR) gene and current testosterone levels:Replication study and meta-analyses: J Neurosci Res, 2020; 98(2); 353-70

16. Vaillancourt KL, Dinsdale NL, Hurd PL, Estrogen receptor 1 promoter polymorphism and digit ratio in men: Am J Hum Biol, 2012; 24(5); 682-89

17. Munro AW, McLean KJ, Grant JL, Makris TM, Structure and function of the cytochrome P450 peroxygenase enzymes: Biochem Soc Trans, 2018; 46(1); 183-96

18. Manikandan P, Nagini S, Cytochrome P450 structure, function and clinical significance:A review: Current Drug Targets, 2018; 19(1); 38-54

19. Luo B, Chen C, Wu X, Cytochrome P450 2U1 is a novel independent prognostic biomarker in breast cancer patients: Front Oncol, 2020; 10; 1379

20. Maksymchuk OV, Kashuba VI, Altered expression of cytochrome P450 enzymes involved in metabolism of androgens and vitamin D in the prostate as a risk factor for prostate cancer: Pharmacol Rep, 2020; 72(5); 1161-72

21. Salehi Z, Gholizadeh L, Vaziri H, Madani AH, Analysis of GSTM1, GSTT1, and CYP1A1 in idiopathic male infertility: Reprod Sci (Thousand Oaks, Calif), 2012; 19(1); 81-85

22. Munawar Lone N, Babar S, Sultan S, Association of the and gene polymorphisms in women with polycystic ovary syndrome from Punjab, Pakistan: Gynecol Endocrinol, 2020 [Online ahead of print]

23. Wang C, Pallan PS, Zhang W, Functional analysis of human cytochrome P450 21A2 variants involved in congenital adrenal hyperplasia: J Biol Chem, 2017; 292(26); 10767-78

24. Hill M, Paskova A, Kanceva R, Steroid profiling in pregnancy:A focus on the human fetus: J Steroid Biochem Mol Biol, 2014; 139; 201-22

25. Li H, Lampe JN, Neonatal cytochrome P450 CYP3A7:A comprehensive review of its role in development, disease, and xenobiotic metabolism: Arch Biochem Biophys, 2019; 673; 108078

26. Reinen J, Vermeulen NPE, Biotransformation of endocrine disrupting compounds by selected phase I and phase II enzymes – formation of estrogenic and chemically reactive metabolites by cytochromes P450 and sulfotransferases: Curr Med Chem, 2015; 22(4); 500-27

27. Burton GJ, Fowden AL, Thornburg KL, Placental origins of chronic disease: Physiol Rev, 2016; 96(4); 1509-65

28. Johnson N, De Ieso P, Migliorini G, Cytochrome P450 allele CYP3A7*1C associates with adverse outcomes in chronic lymphocytic leukemia, breast, and lung cancer: Cancer Res, 2016; 76(6); 1485-93

29. Park JH, Lee J, Kim CH, Lee S, The polymorphism (−600 C>A) of CpG methylation site at the promoter region of CYP17A1 and its association of male infertility and testosterone levels: Gene, 2014; 534(1); 107-12

30. Hosono S, Ito H, Oze I, Polymorphisms in CYP19A1, HSD17B1 and HSD17B2 genes and serum sex hormone level among postmenopausal Japanese women: Maturitas, 2015; 82(4); 394-401

31. Bethke L, Webb E, Sellick G, Polymorphisms in the cytochrome P450 genes CYP1A2, CYP1B1, CYP3A4, CYP3A5, CYP11A1, CYP17A1, CYP19A1 and colorectal cancer risk: BMC Cancer, 2007; 7; 123

32. O’Shaughnessy PJ, Monteiro A, Bhattacharya S, Steroidogenic enzyme expression in the human fetal liver and potential role in the endocrinology of pregnancy: Mol Hum Reprod, 2013; 19(3); 177-87

33. Xu HM, Xu LF, Hou TT, GMDR:Versatile software for detecting gene-gene and gene-environment interactions underlying complex traits: Curr Genomics, 2016; 17(5); 396-402

34. Hizbullah , Ahmed S, Noor Mumtaz M, Genetic variations in drug-metabolizing enzyme CYP2C9 among major ethnic groups of Pakistani population: Gene, 2020; 746; 144659

35. Al-Eitan LN, Rababa’h DM, Alghamdi MA, Khasawneh RH, Association of CYP gene polymorphisms with breast cancer risk and prognostic factors in the Jordanian population: BMC Med Genet, 2019; 20(1); 148

36. Bentz EK, Hefler LA, Kaufmann U, A polymorphism of the CYP17 gene related to sex steroid metabolism is associated with female-to-male but not male-to-female transsexualism: Fertil Steril, 2008; 90(1); 56-59

37. Manning JT, Fink B, Sexual dimorphism in the ontogeny of second (2D) and fourth (4D) digit lengths, and digit ratio (2D:4D): Am J Hum Biol, 2018; 30(4); e23138

38. Swaddle JP, Digit ratio:A pointer to fertility, behavior, and health: Heredity, 2002; 89(5); 403

39. Coyne SM, Manning JT, Ringer L, Bailey L, Directional asymmetry (right–left differences) in digit ratio (2D:4D) predict indirect aggression in women: Pers Individ Differ, 2007; 43(4); 865-72

40. Vaillancourt KL, Dinsdale NL, Hurd PL, Estrogen receptor 1 promoter polymorphism and digit ratio in men: Am J Hum Biol, 2012; 24(5); 682-89

41. Nishimura Y, Moriya K, Kobayashi S, Association between ESR1 polymorphisms and second to fourth digit ratio in school-aged children in the Hokkaido Study: Steroids, 2019; 141; 55-62

42. Nebesio TD, Renbarger JL, Nabhan ZM, Differential effects of hydrocortisone, prednisone, and dexamethasone on hormonal and pharmacokinetic profiles:A pilot study in children with congenital adrenal hyperplasia: Int J Pediatr Endocrinol, 2016; 2016; 17

43. Hu MC, Hsu HJ, Guo IC, Chung BC, Function of Cyp11a1 in animal models: Mol Cell Endocrinol, 2004; 215(1–2); 95-100

44. Terry K, McGrath M, Lee IM, Genetic variation in CYP11A1 and StAR in relation to endometrial cancer risk: Gynecol Oncol, 2010; 117(2); 255-59

45. Hogervorst JGF, van den Brandt PA, Godschalk RWL, Interaction between dietary acrylamide intake and genetic variants for estrogen receptor-positive breast cancer risk: Eur J Nutr, 2019; 58(3); 1033-45

46. Hong L, Zhan-Bing M, Zhi-Yun S, Digit ratio (2D:4D) in Chinese women with breast cancer: Am J Hum Biol, 2014; 26(4); 562-46

47. Poniah P, Mohamed Z, Apalasamy YD, Genetic polymorphisms in the androgen metabolism pathway and risk of prostate cancer in low incidence Malaysian ethnic groups: Int J Clin Exp Med, 2015; 10(8); 8

48. Zhang N, Jia J, Ding Q, Common variant rs11191548 near the CYP17A1 gene is associated with hypertension and the serum 25(OH) D levels in Han Chinese: J Hum Genet, 2018; 63(6); 731-37

49. Crucitta S, Del Re M, Paolieri F, CYP17A1 polymorphism c.-362T>C predicts clinical outcome in metastatic castration-resistance prostate cancer patients treated with abiraterone: Cancer Chemother Pharmacol, 2020; 86(4); 527-33

50. Wang L, Huo Z, Lu H, Digit ratio (2D:4D) and coronary artery disease in north Chinese women: Early Hum Dev, 2018; 116; 64-67

51. Myung S, Han J, Lee Y, Association between cytochrome CYP17A1, CYP3A4, and CYP3A43 polymorphisms and prostate cancer risk and aggressiveness in a Korean study population: Asian J Androl, 2015; 17(2); 285-91

52. Zhanbing M, Jie D, Chunyue B, Association of CYP19A1 single-nucleotide polymorphism with digit ratio (2D:4D) in a sample of men and women from Ningxia (China): Early Hum Dev, 2019; 132; 58-65

53. Armamento-Villareal R, Shah VO, The rs4646 and rs12592697 polymorphisms in CYP19A1 are associated with disease progression among patients with breast cancer from different racial/ethnic backgrounds: Front Genet, 2016; 7; 211

54. Johansson H, Aristarco V, Gandini S, Prognostic impact of genetic variants of CYP19A1 and UGT2B17 in a randomized trial for endocrine-responsive postmenopausal breast cancer: Pharmacogenomics J, 2020; 20(1); 19-26

55. Morel Y, Roucher F, Plotton I, Evolution of steroids during pregnancy:Maternal, placental and fetal synthesis: Ann Endocrinol (Paris), 2016; 77(2); 82-89

Figures

Figure 1. Linkage disequilibrium (LD) map of single-nucleotide polymorphisms (SNPs) in genes CYP17A1 and CYP19A1. A good LD was found among rs1004467, rs2486758, and rs17115149 of CYP17A1(rs1004467 – rs2486758: D’=0.9993, r2=0.1153; rs1004467 – rs17115149: D’=0.6885, r2=0.0202; rs2486758 – rs17115149: D’=0.9964, r2=0.0203); and rs2255192, rs4275794, and rs4646 of CYP19A1(rs2255192 – rs4275794: D’=0.9883, r2=0.9728; rs2255192 – rs4646: D’=0.9992, r2=0.0974; rs4275794 – rs4646: D’=0.9992, r2=0.0978). The haplotypes in the blocks were established based on the LD values.Figure 2. Second-to-fourth digit ratios (2D: 4D) of different sexes distributed in different genotypes of the CYP17A1 single-nucleotide polymorphisms (SNPs). (A–C) Left hand 2D: 4D of rs17115149, rs1004467, and rs2486758 genotypes; (D–F) Right hand 2D: 4D of rs17115149, rs1004467, and rs2486758 genotypes. The black scatter points are individual 2D: 4D observations (including 400 men and 404 women). The box along with the upper and lower vertical lines represent the quartile range of 2D: 4D. Different colors of boxes represent different phenotypes. The horizontal line inside the box represents the mean value of the 2D: 4D. Remeasurement reliability of all digit lengths among 3 operators were acceptable for both hands in the intra-class correlation coefficient test.Figure 3. Single-nucleotide polymorphism (SNP) interaction graphs analyzed by generalized multifactor dimensionality reduction for second-to-fourth digit ratio (2D: 4D). (A, B) Optimal models as determined by generalized multifactor dimensionality reduction (GMDR) for CYP SNPs (0 = 0 minor allele; 1=1 minor allele; 2=2 minor alleles). The numbers within each small square represent positive scores (left) and negative scores (right). For each square, dark blu cells indicate the high effect on the 2D: 4D, and light blue cells represent the low effect on the 2D: 4D. The test was adjusted for sex.

In Press

21 Feb 2024 : Clinical Research  

Potential Value of HSP90α in Prognosis of Triple-Negative Breast Cancer

Med Sci Monit In Press; DOI: 10.12659/MSM.943049  

22 Feb 2024 : Review article  

Differentiation of Native Vertebral Osteomyelitis: A Comprehensive Review of Imaging Techniques and Future ...

Med Sci Monit In Press; DOI: 10.12659/MSM.943168  

23 Feb 2024 : Clinical Research  

A Study of 60 Patients with Low Back Pain to Compare Outcomes Following Magnetotherapy, Ultrasound, Laser, ...

Med Sci Monit In Press; DOI: 10.12659/MSM.943732  

26 Feb 2024 : Clinical Research  

Predictive Value of Combined HbA1c and Neutrophil-to-Lymphocyte Ratio for Diabetic Peripheral Neuropathy in...

Med Sci Monit In Press; DOI: 10.12659/MSM.942509  

Most Viewed Current Articles

17 Jan 2024 : Review article  

Vaccination Guidelines for Pregnant Women: Addressing COVID-19 and the Omicron Variant

DOI :10.12659/MSM.942799

Med Sci Monit 2024; 30:e942799

0:00

16 May 2023 : Clinical Research  

Electrophysiological Testing for an Auditory Processing Disorder and Reading Performance in 54 School Stude...

DOI :10.12659/MSM.940387

Med Sci Monit 2023; 29:e940387

0:00

14 Dec 2022 : Clinical Research  

Prevalence and Variability of Allergen-Specific Immunoglobulin E in Patients with Elevated Tryptase Levels

DOI :10.12659/MSM.937990

Med Sci Monit 2022; 28:e937990

0:00

01 Jan 2022 : Editorial  

Editorial: Current Status of Oral Antiviral Drug Treatments for SARS-CoV-2 Infection in Non-Hospitalized Pa...

DOI :10.12659/MSM.935952

Med Sci Monit 2022; 28:e935952

0:00

Your Privacy

We use cookies to ensure the functionality of our website, to personalize content and advertising, to provide social media features, and to analyze our traffic. If you allow us to do so, we also inform our social media, advertising and analysis partners about your use of our website, You can decise for yourself which categories you you want to deny or allow. Please note that based on your settings not all functionalities of the site are available. View our privacy policy.

Medical Science Monitor eISSN: 1643-3750
Medical Science Monitor eISSN: 1643-3750