FEMALE SEXUAL FUNCTION OF HEALTHY WOMEN IN EASTERN CROATIA.

The aim was to evaluate sexual function of healthy women in Croatia and the possible impact of depression, anxiety, and sociodemographic factors. A total of 204 healthy women filled in a basic sociodemographic questionnaire, the Patient Health Questionnaire-9, Anxiety Disorder-7, and Female Sexual Function Index (FSFI). Almost half of study subjects (47.1%) reported at least some degree of sexual dysfunction defined as an FSFI score lower than 26.55. Study results suggest sexual dysfunction of women in Croatia as a still unrecognized problem. More room in research and in the public must be given to this issue.

Introduction

Sexual function is an integral component of every woman’s life. The impact of sexual health on overall life satisfaction and quality of life (QOL) is increasingly recognized. Female sexual function problems are multidimensional, associated with general health, affected by medication use or substance abuse, but also inseparable from relationship characteristics, cultural influences, socioeconomic status, and religious beliefs (). Sexual problems have a serious impact on mood, self-esteem, and QOL. They may lead to emotional distress, relationship problems, divorce, or they may even affect reproduction (). Sexual dysfunction (SD) is an important public health issue. According to epidemiological studies, sexual problems are widespread among women of all ages (-), ranging from 19% to 63% (). Female sexual dysfunction (FSD) is still largely unexplored, and actual prevalence figures vary widely due to differences in diagnostic criteria (). The International Consensus Development Conference on Female Sexual Dysfunctions classified FSD into four disorders () and developed the Female Sexual Function Index (FSFI), which assesses sexual function within six domains: desire, arousal, orgasm, pain, lubrication, and satisfaction (). Accurate estimates of the prevalence/incidence of FSD and identifying risk factors are important for prevention efforts and treatment planning. The prevalence of FSD in Europe and the United States has been investigated in many recent studies; however, data from Croatia on the prevalence and risk factors for this condition are rare. To the authors’ knowledge, only two research studies on female sexual problems in Croatia taking into account different domains of sexual function were conducted, one of them providing information on the prevalence of sexual problems and related risk factors in a community sample in the Croatian capital of Zagreb (), while the other investigated prevalence of sexual health difficulties among women aged 18-35 across all 21 counties of six major regions of Croatia (). There were no studies using the FSFI, but data relied solely on the participants’ statement that they had difficulties in different aspects of sexual functioning.

The aim of this study was to assess the level of sexual functioning and the prevalence of SD in healthy, sexually active women in eastern Croatia.

Materials and Methods

Participants

A cross-sectional study was conducted on 204 healthy female subjects aged 18-65 years (median age 42, IQR 32-51) from eastern Croatia as part of a larger study in which sexual function will be assessed in various clinical samples. Participants were gathered among acquaintances of the authors. Thus, the sample was a convenient one. Although a precise response rate was not determined, there were almost no declines to participate in the study. Women with severe somatic or psychiatric conditions, pregnant women, those under 18 and over 70 years of age, as well as those who declined to give their signed consent were excluded from the study. The study was approved by the Ethics Committees of the Osijek Faculty of Medicine and Osijek University Hospital Centre, and the study protocol conformed to the provisions of the Declaration of Helsinki. All study participants provided their written informed consent.

Questionnaire

The participants were presented a self-administered questionnaire consisting of four parts. In the first part, the participants were asked about their basic sociodemographic characteristics: age, current relationship status, type of relationship (marriage, extramarital relationship, same-sex partnership, and none of the above), place of residence (city, town, suburb, and village), education and parity.

The second part consisted of the Patient Health Questionnaire-9 (PHQ-9). PHQ-9 is a validated instrument with high sensitivity and specificity, consisting of nine questions that are answered on a Likert scale, each providing information on one of the nine criteria for depression as defined by Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV) (-).

The third part consisted of the Generalized Anxiety Disorder-7 (GAD-7) questionnaire. GAD-7 is a validated instrument developed to be used as an orientational evaluation of the anxiety symptoms. It consists of seven questions answered on a Likert scale (-).

The fourth part included the FSFI, a validated 19-item questionnaire that measures six domains of female sexual function (desire, arousal, lubrication, orgasm, satisfaction, and pain) () in four weeks before the examination. The answers vary from 0 (no sexual activity) or 1 to 5. Individual domain scores are defined as the sum of the scores of the questions pertaining to the domain, multiplied by a factor provided by the authors of the scale. Total FSFI is obtained by adding the six domain scores. Total score ranges from 2 to 36, with higher scores signifying better sexual function. A cutoff score of 26.55, as suggested by Wiegel et al., was used to determine SD ().

Since all the participants spoke Croatian, validated Croatian versions of the PHQ-9 and GAD-7 were used. Since, to the authors’ knowledge, Croatian version of the FSFI had not been previously prepared, a team of translators did translation using a modified back-translation method proposed by several authors (-). Two fluently bilingual physicians each translated the FSFI into Croatian. The two versions were compared and synthesized into a single translation by a team of three bilingual physicians. The aforementioned translation was then back-translated into English by a Croatian translator with master degree in English language and literature. The back-translation was then compared with the English original to assure linguistic validity by several fluently bilingual medical professionals and translators.

Statistical analysis

The sample size required to maintain a power of 80% and alpha error probability of 0.05 was calculated using the G*Power 3.1 software (, ). Kolmogorov-Smirnov test was used to test the normality of data distribution. Since all scalar variables deviated from normal distribution, nonparametric tests were used. Numerical data were expressed as median and interquartile range (IQR). To compare the means of two or more independent groups, Mann-Whitney and Kruskal-Wallis test was used, respectively. Correlation between variables was assessed with the Spearman’s rank-order correlation coefficient. The level of statistical significance was set at p<0.05. Multivariate analysis was performed to determine independent association of variables with the FSFI score. The analysis was conducted using the SPSS software (ver. 16.0, SPSS Inc., Chicago, IL, USA).

Results

Basic sociodemographic characteristics of the participants are shown in Table 1.

Table 1

Sociodemographic characteristics of study participants

	N=204
	n (%)
Current relationship status
In a relationship	172 (84.3)
Not in a relationship	32 (15.7)
Type of relationship
Marriage	126 (61.8)
Extramarital relationship	29 (14.2)
Same-sex partnership	0
None of the above	49 (24)
Place of residence
City	109 (53.4)
Small town	36 (17.6)
Suburb	22 (10.8)
Village	37 (18.1)
Level of education
Elementary school	11 (5.4)
High school	107 (52.5)
Bachelor’s degree	15 (6.9)
Four-year college degree, or higher	68 (33.3)
Number of children
0	60 (29.4)
1	41 (20.1)
2	71 (34.8)
3	18 (8.8)
4	8 (3.9)
5	3 (1.5)
6	3 (1.5)
Have you undergone serious surgery in your lifetime?
Yes	33 (16.2)
No	171 (83.8)

As shown in Table 1, 32 participants were not in a relationship in the month prior to the examination. Table 2 shows differences between the two groups divided according to relationship status. No differences were found in depression or anxiety score. As expected, significant differences were recorded for sexual function. Only the participants having reported a relationship were included in further analysis.

Table 2

Differences in PHQ-9, GAD-7, FSFI scores and age according to participant relationship status one month prior to the examination (N=204)

	Median (IQR)
	PHQ-9	GAD-7	FSFI	Age (yrs)
Not in a relationship (n=32)	7.5 (4.25-10)	6 (3.25-9)	5.2 (2.6-9.6)	49.5 (25-57.75)
In a relationship (n=172)	5 (2.25-9)	5 (2.25-8)	27.85 (22.7-31)	42 (33-50)
p value (Mann-Whitney test)	0.062	0.247	<0.001	0.565

PHQ-9 = Patient Health Questionnaire-9; GAD-7 = Generalized Anxiety Disorder-7; FSFI = Female Sexual Function Index; IQR = interquartile range

When considering the cut-off scores for PHQ-9, GAD-7 and FSFI, study participants were divided as follows: 110 (53.9%) subjects had no or very low level of anxiety, 56 (27.5%) had mild, 31 (15.2%) moderate and three (1.5%) subjects severe anxiety. Depressive symptoms were absent in 157 (77%) subjects, 34 (16.7%) had mild depressive symptoms, ten (4.9%) moderate and two (1%) subjects severe depressive symptoms.

In the group of women in a relationship, 91 (52.9%) had FSFI score over 26.55, indicating normal sexual function, whereas 81 (47.1%) showed some degree of sexual dysfunction. Sixteen (9.3%) women reported having had no sexual activity in the month prior to the examination.

There were no significant differences in age (p=0.4, Kruskal-Wallis test), PHQ-9 (p=0.874, Kruskal-Wallis test), GAD-7 (p=0.571, Kruskal-Wallis test) or FSFI score (p=0.063, Kruskal-Wallis test) according to the place of residence.

There were no significant differences in PHQ-9 (p=0.118, Kruskal-Wallis test), GAD-7 (p=0.578, Kruskal-Wallis test) or FSFI score (p=0.073, Kruskal-Wallis test) between the groups divided according to the level of education. However, there was a significant age difference according to the level of education, i.e. younger women had higher levels of education.

There was no significant difference in FSFI score according to the number of children (p=0.222, Kruskal-Wallis test), but there was a weak significant negative correlation between the number of children and FSFI score (Spearman’s rho=-0.166, p=0.029).

A significant difference was found in the level of anxiety according to the number of children (p=0.04, Kruskal-Wallis test), but there was no significant correlation between these variables (Spearman’s rho=0.086, p=0.222). A significant age difference was also found according to the number of children (p<0.001, Kruskal-Wallis test), with age positively correlating with the number of children (Spearman’s rho=0.596, p<0.001). There were no differences in PHQ-9 scores between the groups divided according to the number of children (p=0.351, Kruskal-Wallis test).

There were no differences in age (p=0.49, Kruskal-Wallis test), PHQ-9 (p=0.7, Kruskal-Wallis test) or GAD-7 (p=0.164, Kruskal-Wallis test) score related to the subjects’ having had sexual activity in the month prior to the examination or not.

A significant difference in the FSFI score was found according to the type of relationship (p=0.001, Kruskal-Wallis test), with married women (median 26.1, IQR 21.6-30.5) having lower scores than those in extramarital relationships (median 31, IQR 25.45-32.55) and those stating that they were in other types of relationship (median 30.15, IQR 25.2-32.75) (Table 2).

There was a significant age difference between the groups based on marital status (p<0.001, Kruskal-Wallis test). Married women were oldest (median 46, IQR 38-53), women in extramarital relationships were second oldest (median 32, IQR 25.5-42.5), and those in other kinds of relationship were youngest (median 26, IQR 20-36.5).

Median scores of the individual domains of the FSFI questionnaire, as well as median of total FSFI scores are presented in Table 3.

Table 3

FSFI domains and total FSFI scores (N=172)

	Median (IQR)
Desire	3.6 (3-3.6)
Arousal	4.2 (3.3-4.8)
Lubrication	4.95 (3.6-5.7)
Orgasm	4.8 (3.6-6)
Satisfaction	5.2 (4-6)
Pain	6 (4.4-6)
FSFI	27.85 (22.7-31)

FSFI = Female Sexual Function Index; IQR = interquartile range

Tables 4 and 5 illustrate correlations between scalar variables.

Table 4

Correlations (Spearman’s rho) of age, GAD-7, PHQ-9 and FSFI domain, and total scores of participants in relationship (N=172)

	GAD-7	PHQ-9	Desire	Arousal	Lubrication	Orgasm	Satisfaction	Pain	FSFI
Age	-0.023	-0.035	-0.349**	-0.36**	-0.215**	-0.241**	-0.205**	-0.056	-0.299**
GAD-7		0.827**	-0.193*	-0.224**	-0.094	-0.139	-0.19*	-0.151*	-0.173*
PHQ-9			-0.219**	-0.231**	-0.126	-0.219**	-0.257**	-0.22**	-0.23**
Desire				0.598**	0.419**	0.354**	0.456**	0.271**	0.609**
Arousal					0.645**	0.58**	0.547**	0.516**	0.771**
Lubrication						0.723**	0.667**	0.535**	0.865**
Orgasm							0.701**	0.458**	0.84**
Satisfaction								0.463**	0.838**
Pain									0.601**

PHQ-9 = Patient Health Questionnaire-9; GAD-7 = Generalized Anxiety Disorder-7; FSFI = Female Sexual Function Index; *p<0.05; **p<0.01

A stepwise multiple linear regression analysis was performed to determine how different variables affected the FSFI score. The number of children, GAD-7 and PHQ-9 scores were excluded from analysis. A model was created taking into account age and relationship status. The model explained 34.2% of variance as determined by the adjusted R Square (standard error (SE) 5.93). ANOVA suggested that the model had a significant explanatory power, F=50.069, p<0.001, df=2.

Both independent variables included in the model significantly affected the FSFI score, relationship status had unstandardized B 11.736, SE 1.236, t=9.493, p<0.001, while age had unstandardized B -0.090, SE 0.037, t=-2.419, p=0.017.

Discussion

Female sexual dysfunction is defined as a disorder of libido, arousal, orgasm, physical pleasure, and pain related to sexual intercourse (, ). SD has a significant impact on QOL and personal relationships, including aspects such as marriage satisfaction and social performance (). The purpose of our analysis was to provide insights into the female sexual function in eastern Croatia.

When using the cut-off score of 26.55 suggested by Wiegel et al. (), 47.1% of study participants showed some degree of SD (). This prevalence is similar to that found in similar studies conducted in Croatia, where 47% () and 37% () of women were found to have some difficulties. However, direct comparison of our results with other results from Croatia was not possible due to different measures used in those studies. Our results are also in line with the results of research in other areas. A study in the USA found the prevalence of FSD among women in the general population to be 43% (), and another study in England reports it to be 41% ().

Since our sample was a convenient one, with the participants gathered from the pool of the authors’ acquaintances, and since there were virtually no declines to participate in the study, we postulated that our results would be somewhat lower than those reported in the literature. Response rates in similar studies are most commonly below 50% (, ), and the women willing to participate are presumably more open, liberal and with better sexual function, whereas those who decline to participate presumably have more difficulties with their sexual function, and are therefore less likely to answer questions pertaining to their sexual life (). According to the Census of Population, Households, and Dwellings 2011, 86.28% of Croatian citizens were Catholic. Sexually restrictive religious morality can be responsible for reduction in the ability to focus on sexual feelings and sensations (). Such norms can block sexual expression by instilling feelings of guilt, at least in women.

Married women in our study had lower FSFI scores than both the women in extramarital relationships and those in other types of relationships, which is discordant with the results of both previous studies conducted in Croatia. The fact that partner communication and marriage significantly lower the odds of having sexual difficulties has been reported in many studies (, ), confirming the importance of relationship quality for female sexual health. The possible reason for the results in our sample might be age difference between the groups, with married women being oldest, and it is well documented that age leads to a decline in sexual function (, ). This was also confirmed in our study by the results of multiple linear regression analysis, which showed that age and relationship status explained 34.2% of variance with a significant explanatory power.

We found no significant differences in sexual function based on the number of children. This is in contrast with the findings of other studies, where women with small children were more likely to report sexual problems (, ).

Loss of interest, reduction in energy and inability to experience pleasure are characteristics of depression that may produce difficulties in sexual relationships. One study showed that as much as 72% of patients with unipolar depression reported loss of sexual interest (). Furthermore, comparative studies found higher levels of SD in depressed patients than in healthy controls (-). Although the specific type of sexual dysfunction varies in incidence, the most common complaint is loss of sexual desire, as was the case in our sample. In our sample, depression negatively correlated with every domain of sexual function, except for lubrication, which was expected and previously documented.

Aging affects a woman’s sexual function through a range of hormonal and physiological changes. There is a number of cross-sectional studies providing strong evidence that sexual function declines with age (, , ). On the other hand, SDs do not change or even decline with age (, -). This is interesting because one would expect an increase in SD given the evidence for sexual function decline with age. However, the key may be changes in sexually related distress, since the importance of sex decreases in older women ().

According to some studies, the frequency of SD in postmenopausal women was up to 81% (, ), which was accordant to the study by Peeyananjarassri et al., who found that 82.2% of postmenopausal women suffered from SD (). Our results are in line with those reported in the literature. There was a negative correlation between the FSFI scores and all its domains except for pain, and age.

A study performed in the USA found that approximately 23% of postmenopausal women reported problems with feeling sexual desire, 13% had arousal problems, and even 64% of postmenopausal women had not experienced orgasm during the intercourse. About 22% had no lubrication during sexual activity ().

Even in early psychodynamic formulations, SD anxiety was seen as one of the key factors leading to and promoting SD (). In our sample, the anxiety levels negatively correlated with sexual function, thus giving additional support to previous findings.

The study suffered from some limitations, as follows: our sample size was too small to represent total population of eastern Croatia, and the results presented here are preliminary results of a larger study. Much larger, multicenter studies are needed to provide better insight into the situation of sexual function of women in Croatia. The sample was a convenient one, with participants being gathered from the pool of the authors’ acquaintances. Although this could be considered a limitation, the response rate was far higher than in similar studies. Questions about the participants’ religious views, character, and additional questions about their relationships and sexual behavior beyond the scope of the FSFI might shed more light on the issue of the high prevalence of SD in Croatian women.

Conclusion

In this cross-sectional study, we questioned healthy women from eastern Croatia about their sexual function and other related parameters. Almost half of our subjects (47.1%) reported at least some degree of SD defined as an FSFI score lower than 26.55. Lower levels of sexual function were present in older women and those who showed higher levels of anxiety and depressive symptoms. Our results clearly suggest a great problem that is still unrecognized, i.e. SD of women in Croatia. This issue certainly requires more room in research and public sphere. Medical professionals should not neglect this aspect of a woman’s life, so as to help the woman increase her sense of wellbeing, self-esteem and QOL.