Open Access

Coincident frequencies and relative phases among brain activity and hormonal signals

  • Silvia Solís-Ortíz1,
  • Rafael G Campos2,
  • Julián Félix3 and
  • Octavio Obregón3Email author
Behavioral and Brain Functions20095:18

DOI: 10.1186/1744-9081-5-18

Received: 13 October 2008

Accepted: 14 March 2009

Published: 14 March 2009

Abstract

Background

Fourier transform is a basic tool for analyzing biological signals and is computed for a finite sequence of data sample. The electroencephalographic (EEG) signals analyzed with this method provide only information based on the frequency range, for short periods. In some cases, for long periods it can be useful to know whether EEG signals coincide or have a relative phase between them or with other biological signals. Some studies have evidenced that sex hormones and EEG signals show oscillations in their frequencies across a period of 28 days; so it seems of relevance to seek after possible patterns relating EEG signals and endogenous sex hormones, assumed as long time-periodic functions to determine their typical periods, frequencies and relative phases.

Methods

In this work we propose a method that can be used to analyze brain signals and hormonal levels and obtain frequencies and relative phases among them. This method involves the application of a discrete Fourier Transform on previously reported datasets of absolute power of brain signals delta, theta, alpha1, alpha2, beta1 and beta2 and the endogenous estrogen and progesterone levels along 28 days.

Results

Applying the proposed method to exemplary datasets and comparing each brain signal with both sex hormones signals, we found a characteristic profile of coincident periods and typical relative phases. For the corresponding coincident periods the progesterone seems to be essentially in phase with theta, alpha1, alpha2 and beta1, while delta and beta2 go oppositely. For the relevant coincident periods, the estrogen goes in phase with delta and theta and goes oppositely with alpha2.

Conclusion

Findings suggest that the procedure applied here provides a method to analyze typical frequencies, or periods and phases between signals with the same period. It generates specific patterns for brain signals and hormones and relations among them.

Background

Fourier transform is a basic tool for analyzing biological signals [13]. Mostly a fast Fourier transform is computed for finite sequence of data sample [4]. In the case of brain signals obtained from electroencephalographic (EEG) recording, the information can be extracted by spectral analysis, whereby amplitude characteristics of the frequency domain of the EEG signal can be assessed using a technique known as fast Fourier transform, specific frequency bands are identified and the signal amplitude within each of these wavebands calculated for short periods [5, 6]. Generally this analysis provides information, according with the frequency range of the well known brain signals delta (1.5 – 3.5 Hz), theta (3.3 – 7.5 Hz), alpha1 (7.5 – 9.5 Hz), alpha2 (9.5 – 12.5 Hz), beta1 (12.5 – 17.5 Hz), beta2 (17.5 – 30.0 Hz) and gamma waves (30–50 Hz), which have been widely used both in clinical and experimental designs [711]. However, the EEG signals analyzed with this method provide only information based in frequency range, for short periods. In some cases, for long periods it can be useful to know whether EEG signals coincide or have a relative phase between them or with other biological signals of scientific interest. The endogenous sex hormones are biological signals of interest in many aspects [12], and some studies have evidenced that these signals and EEG signals show oscillations in their frequencies across a period of 28 days [1319], so it seems of relevance to seek after possible patterns relating with EEG signals and endogenous sex hormones data. In the present paper, we propose one method that can be used to obtain coincidence and relative phases among EEG signals and hormonal levels for long periods. The procedure is established as follows: Assigning certain numerical value, i.e., the absolute power, to each EEG signal at certain sampling times, generates data that can be interpolated through a long period, yielding an absolute power function of time for each signal. A similar procedure is followed in case of hormonal values. A discrete Fourier transform is then performed to analyze these new functions, finding typical frequencies and their corresponding periods for each one of these signals and, also, relative phases for coincident periods between two or more signals. The technique is based on a discrete Fourier transform which gives accurate results even in the case of few point samples [20, 21]. We used exemplary dataset from real frequencies EEG recording and values of hormonal levels obtained along of 28 days to illustrate this method. In this way, we found typical frequencies and characteristics phases of the brain signals and of the progesterone and estrogen levels.

Methods

Brain signals and hormonal levels

We used a dataset of EEG recording of women and the variations of their estrogen and progesterone levels previously reported [17, 22, 23]. Briefly, 9 female volunteers aged between 20 and 34 (mean = 27.11; SD = 4.45) participated in the study. They were university students and they were selected on account of their regular cycles (28 days ± 2). None of them were under any type of medication nor were taking or have taken oral contraceptives. A total of 10 sessions of the EEG activity were recorded, 3 times a week every second day, except Sundays, at the same time of the day (between 16 and 19 hrs), at rest with eyes open. EEG recordings started at the subject own convenience, thus the day of menstruation of different subjects were randomly distributed over the sessions. Electrodes were placed according to the 10–20 International System al F3, F4, C3, C4, P3, P4, O1 and O2, referred to ipsilateral earlobes. EEG was recorded on a Grass model 16E polygraph set to pass frequencies between 1 and 35 Hz. Impedance was kept below 10 k 0 hms. EEG was captured in a PC computer using the program CAPTUSEN [24], by a 12-bits resolution analog/digital converter, at a sampling rate of 128 Hz. The stored EEG signals were segmented into non-overlapping 2.048s epochs and were carefully inspected offline to eliminate artifacts. Absolute power was obtained for each epoch using the program POTENCOR [25]. Absolute power values were averaged over frequency bins to obtain the following broad bands: delta (1.5–3.5 Hz), theta (3.5–7.5 Hz), alpha1 (7.5–9.5 Hz), alpha2 (9.5–12.5 Hz), beta1 (12.5–17.5 Hz) and beta2 (17.5–30 Hz). The resulting values were then averaged over all epochs for each period of the cycle and for each subject. Note that this time span has a total number of 10 epochs for each subject of lengths 2, 4, 1, 5, 1, 6, 1, 3, 1, and 3 days respectively. The absolute power of the brain signals, as the described EEG recording delivered them, was obtained during the days 1, 3, 7, 8, 13, 14, 20, 21, 24 and 25, corresponding to a cycle. Table 1 shows the average of these absolute powers for each one of the brain signals. Table 2 shows the values of progesterone and estrogen levels during 28 days in an arbitrary scale convenient for the Fourier analysis.
Table 1

Primary absolute average power brain signals, as the EEG delivered them, from the nine young women.

Brain signals (Hz) (absolute power)

Day

1

3

7

8

13

14

20

21

24

25

Delta

388.7

399.0

336.8

402.6

366.8

376.2

388.6

375.3

382.8

396.3

Theta

374.5

359.4

362.9

368.5

339.6

346.9

384.9

320.9

373.4

363.7

Alpha1

378.6

329.6

417.3

323.5

331.5

283.6

413.0

304.4

390.7

360.0

Alpha2

230.5

231.3

194.6

270.3

256.4

345.1

262.5

264.7

274.5

392.1

Beta1

137.0

154.8

129.4

137.5

126.8

134.2

135.6

145.7

160.4

130.5

Beta2

172.5

167.1

156.6

159.0

139.3

144.2

146.4

170.6

213.5

140.8

Menstrual cycle days are divided as follows: 1,3: Menstrual phase; 7,8: Follicular phase; 13,14: Ovulatory phase; 20,21: Early luteal phase; 24,25: Late luteal phase.

Table 2

Progesterone and estrogen levels rescaled in appropriate units suitable for the described Fourier analysis through 28 days.

Hormone/day

1

2

3

4

5

6

7

8

9

10

Progesterone

0.04

0.04

0.15

1.15

0.90

0.03

0.13

0.40

0.13

0.64

Estrogen

0.65

0.80

0.20

0.20

0.20

0.20

0.65

0.55

0.70

0.60

Hormone/day

11

12

13

14

15

16

17

18

19

20

Progesterone

0.35

0.15

0.50

0.15

0.05

1.35

0.85

2.15

3.00

5.00

Estrogen

0.75

1.00

2.40

3.10

2.30

1.90

1.10

0.50

0.85

0.65

Hormone/day

21

22

23

24

25

26

27

28

  

Progesterone

4.40

4.65

3.18

2.00

1.05

0.65

0.35

0.15

  

Estrogen

1.40

1.05

1.00

1.05

1.00

1.35

0.40

0.35

  

The discrete Fourier transform

We present here the method of analysis used to study the brain waves and sex hormone levels. Our discrete Fourier transform is presented in this section and proofs and examples are given elsewhere [20, 21, 26]. This nonstandard Hermite-Gaussian quadrature yields a truly discrete Fourier transform which can be used to compute the spectrum of a sampled nonperiodic (or periodic) signal and to break it down into their sinusoidal components of specific frequencies. Some of the numerical examples given in [20, 26] deal with singular or oscillatory signals for which our discrete Fourier transform performs well, whereas the usual discrete Fourier transform, i.e., the so-called Fast Fourier Transform, fails to yield the correct spectrum. Our discrete Fourier transform can also be used in the frequency analysis of transient signals, like those appearing in some voice recognition problems [27].

Let us denote by f ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaaaaa@2D3A@ (ω) the Fourier transform of the function f(t) defined by
f ^ ( ω ) = e i ω t f ( t ) d t , MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaacqGGOaakcqaHjpWDcqGGPaqkcqGH9aqpdaWdXaqaaiabdwgaLnaaCaaaleqabaGaemyAaKMaeqyYdCNaemiDaqhaaOGaemOzayMaeiikaGIaemiDaqNaeiykaKIaemizaqMaemiDaqhaleaacqGHsislcqGHEisPaeaacqGHEisPa0Gaey4kIipakiabcYcaSaaa@45FE@
(1)
and let H N (t) and t1, t2,...,t N be the N th Hermite polynomial and its N zeros, respectively. Then, if f(t) is a function satisfying certain integrability conditions, the following quadrature formula holds:
f ^ ( ω j ) = e i ω j t f ( t ) d t = k = 1 N F j k f ( t k ) + R j , j = 1 , 2 , ... , N . MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabeGaaaqaaiqbdAgaMzaajaGaeiikaGIaeqyYdC3aaSbaaSqaaiabdQgaQbqabaGccqGGPaqkcqGH9aqpdaWdXaqaaiabdwgaLnaaCaaaleqabaGaemyAaKMaeqyYdC3aaSbaaWqaaiabdQgaQbqabaWccqWG0baDaaGccqWGMbGzcqGGOaakcqWG0baDcqGGPaqkcqWGKbazcqWG0baDaSqaaiabgkHiTiabg6HiLcqaaiabg6HiLcqdcqGHRiI8aOGaeyypa0ZaaabCaeaacqWGgbGrdaWgaaWcbaGaemOAaOMaem4AaSgabeaakiabdAgaMjabcIcaOiabdsha0naaBaaaleaacqWGRbWAaeqaaOGaeiykaKIaey4kaSIaemOuai1aaSbaaSqaaiabdQgaQbqabaaabaGaem4AaSMaeyypa0JaeGymaedabaGaemOta4eaniabggHiLdGccqGGSaalaeaacqWGQbGAcqGH9aqpcqaIXaqmcqGGSaalcqaIYaGmcqGGSaalcqGGUaGlcqGGUaGlcqGGUaGlcqGGSaalcqWGobGtcqGGUaGlaaaaaa@6A88@
(2)
Here, ω j is again a zero of H N (t) and its numerical value corresponds to t j , i.e., ω j = t j , and
F j k = 2 π ( 1 ) j + k 2 N 1 ( N 1 ) ! N H N 1 ( ω j ) H N 1 ( t k ) l = 0 N 1 i l 2 l l ! H l ( ω j ) H l ( t k ) . MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aaSbaaSqaaiabdQgaQjabdUgaRbqabaGccqGH9aqpjuaGdaWcaaqaamaakaaabaGaeGOmaiJaeqiWdahabeaacqGGOaakcqGHsislcqaIXaqmcqGGPaqkdaahaaqabeaacqWGQbGAcqGHRaWkcqWGRbWAaaGaeGOmaiZaaWbaaeqabaGaemOta4KaeyOeI0IaeGymaedaaiabcIcaOiabd6eaojabgkHiTiabigdaXiabcMcaPiabcgcaHaqaaiabd6eaojabdIeainaaBaaabaGaemOta4KaeyOeI0IaeGymaedabeaacqGGOaakcqaHjpWDdaWgaaqaaiabdQgaQbqabaGaeiykaKIaemisaG0aaSbaaeaacqWGobGtcqGHsislcqaIXaqmaeqaaiabcIcaOiabdsha0naaBaaabaGaem4AaSgabeaacqGGPaqkaaGcdaaeWbqaaKqbaoaalaaabaGaemyAaK2aaWbaaeqabaGaemiBaWgaaaqaaiabikdaYmaaCaaabeqaaiabdYgaSbaacqWGSbaBcqGGHaqiaaGccqWGibasdaWgaaWcbaGaemiBaWgabeaakiabcIcaOiabeM8a3naaBaaaleaacqWGQbGAaeqaaOGaeiykaKIaemisaG0aaSbaaSqaaiabdYgaSbqabaGccqGGOaakcqWG0baDdaWgaaWcbaGaem4AaSgabeaakiabcMcaPaWcbaGaemiBaWMaeyypa0JaeGimaadabaGaemOta4KaeyOeI0IaeGymaedaniabggHiLdGccqGGUaGlaaa@798E@
(3)
The component R j stands for the error obtained when the integral transform is substituted by the sum of the right-hand side. Note that (2) can be written in matrix form as
f ^ = F f + R , MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaacqGH9aqpcqWGgbGrcqWGMbGzcqGHRaWkcqWGsbGucqGGSaalaaa@33E7@
(4)

where f ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaaaaa@2D3A@ and f are the vectors whose elements are given by the Fourier transform f ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaaaaa@2D3A@ (ω) and the function f(t) evaluated at the zeros of H N (t) respectively, F is the matrix whose elements are given by (3) and R is the residual vector whose elements are of order O MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWenfgDOvwBHrxAJfwnHbqeg0uy0HwzTfgDPnwy1aaceaGae8NdX=eaaa@3765@ (1/N). This upper bound for the error takes into account critical cases where the signal to be transformed may have singularities or a slow decay at infinity. Thus, F is a matrix representing the Fourier transform in a vector space of finite dimension.

Analysis of biological signals

Let f(t) be one of the signals studied here. Since the sex hormones levels can be represented by periodic functions of 28-day period and our aim is to find correlations between these and the brain signals, we assume that any signal f(t) is a 28-day periodic function which is sampled at M time values (days) 0 <τ1, τ2,...τ M < 28. The first part of the analysis of the signals (basically the calculation of the Fourier transform f ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaaaaa@2D3A@ (ω) of f(t)) was carried out according to the following steps:
  1. 1.

    Select a number N of zeros of a Hermite polynomial. It should be chosen great enough to reduce the residual vector R but small enough to compute the matrix F with significant numerical precision.

     
  2. 2.

    Compute the Fourier matrix F.

     
  3. 3.

    Interpolate the M values f(τ k ) with a trigonometric polynomial.

     
  4. 4.

    Shift and scale the interpolated function to the interval [-π, π] to yield a 2π-periodic signal f ˜ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaGaaaaa@2D39@ (t) (this step is not necessary, but convenient from the numerical point of view).

     
  5. 5.

    Evaluate the polynomial f ˜ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaGaaaaa@2D39@ (t) obtained in the previous step at the N zeros of Hermite to yield the vector f.

     
  6. 6.

    Compute the Fourier vector through the relation f ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaaaaa@2D3A@ = Ff.

     
  7. 7.

    Interpolate to a trigonometric polynomial f ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaaaaa@2D3A@ (ω) the values f ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaaaaa@2D3A@ k .

     

Basically, the sequence of analysis is: first, interpolate the samples to give a periodic function and second, break this function down into their sinusoidal components through our discrete Fourier transform.

The second part of the analysis was done according to the following. Since the Fourier transform f ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmOzayMbaKaaaaa@2D3A@ (ω) of the signal is composed by a real and an imaginary part, it is necessary to determine if one of them or both are important for the frequency analysis. Thus, we first compute the Fourier transforms of each signal and take apart the real and imaginary components. Concerning the brain signals, it comes out that the values of the real part of the transformed signals at ω ≠ 0 are all negligible compared with the great values they take at ω = 0. This behavior, which is typical of the Fourier Transform of a constant function, is common to all the real part of the transformed signals, whereas the imaginary part varies significantly among the different brain signals. Thus, the real part of a transformed signal is more related to the power delivered than the imaginary part, which is more related to the frequency dependence of the signal. Therefore, only the imaginary part of the transformed brain signals will be taken into account in our analysis. On the other hand, in the case of the Progesterone and Estrogen levels, the real part of the Fourier Transform comes out as important as the imaginary part, consequently both of them are taken into account.

Each signal was normalized to the height of its biggest peak. This is not necessary for Fourier analysis, but it is convenient to make uniform the data in order to compare them. Clearly, the coincident maxima define common frequencies for some signals at frequencies 1, 2, 3 and 4 units, or equivalently, at periods of 28, 14, 10 and 7 days (A frequency unit means that the signal will repeat itself each 28 days, a higher frequency means a shorter period). For coincident frequencies we will analyze only those frequencies where the value of the transformed signals are above 1/3. This percentile value gives a criterion for what "dominant" means and in this way we safely cut out all the contribution of the exposed technique to the data background, ensuring that we treat mainly the signal of interest.

There appear dominant frequencies (A word of caution is important. These frequencies correspond to the main sinusoidal components of a given signal and are related to the time needed to have that main component back again. The EEG frequency is of different nature). Some of them coincide, in particular, for a brain signal and one of the hormones. A pair of signals (or more) centered around a characteristic frequency of interest is taken into account in our study only if both of them have the value 1/3 at least, as already mentioned. This condition seems arbitrary, however the obtained results are plausible. For values smaller than 1/3, signals probably arise from noise generated during the discrete Fourier transform. The results correspond to trigonometric functions which main components present a novel and interesting relative phase. Therefore, signals with equal periods present, in general, a phase between them that can be estimated directly from a simple plot. In order to get the relative phase between two coincident periods, we consider only the values around the desired frequency. These windowed signals are then back-transformed to the time variable to get their relative phase by comparing the corresponding Fourier components. By this procedure, the EEG recording and progesterone and estrogen levels are characterized and related through the use of the Fourier transform we have selected.

Results

The transformed signals are plotted in Figure 1 against frequency in units of 2π/28 days-1. The frequency patterns for each one of the EEG-signals and the two hormone levels are plotted below for positive values of the frequencies. As it can be seen, each one of these eight signals presents a peculiar pattern, a characteristic spectrum, each one in different color and normalized to 1. In blue are plotted the estrogen signals; in red, the progesterone signals; in green, the brain signals. In (a), according to the selection criterion, the spectrum of the delta brain signal versus estrogen shows a coincidence in 4; versus progesterone, a coincidence in 2. In (b), the spectrum of the theta brain signal versus estrogen shows a coincidence in 4; versus progesterone, a coincidence in 3. In (c), the spectrum of the alpha1 brain signal versus estrogen shows non coincidence; versus progesterone, a coincidence in 2 and 3. In (d), the spectrum of the alpha2 brain signal versus estrogen shows a coincidence in 1; versus progesterone, a coincidence in 1, 2, and 3. In (e), the spectrum of the beta1 brain signal versus estrogen shows non coincidence; versus progesterone, a coincidence in 2 and 3. In (f), the spectrum of the beta2 brain signal versus estrogen shows non coincidence; versus progesterone, a coincidence in 3. The Fourier transform of each brain signal and hormone levels were plotted superimposed each other, to analyze relative phases between them. When one of the relevant brain signal periods essentially coincides with one of the hormone levels, the relative phase between the two corresponding time functions is displayed in Figure 2 and Figure 3. It should be remarked that the phases presented in each one of these figures correspond to the relevant trigonometric components of the whole signal, centered around the frequency of interest, for we have related coincident relevant frequencies of the two signals considered.
Figure 1

Intensity of the different brain signals superimposed to hormone level signals. The horizontal scale is 2π/28 days-1. The vertical scale has normalized arbitrary units. It is shown period coincidences and relative intensities in different points; when the peaks agree and both are bigger than 1/3 the coincidences were taken into account for this analysis, considering the smaller peaks as noise generated from the discrete Fourier transform.

Figure 2

Estrogen signals versus brain signals. Red dots and red continuous lines are used indistinctly to clarify the figures. Delta and theta are in phase with estrogen signal; alpha2, in relative phase of 10 days.

Figure 3

Progesterone signals versus the six brain signals. Different relative phases are presented between brain signals and progesterone signal.

In Figure 2(a), the estrogen levels and delta, with 7-day period, show a relative phase of 0 days. The estrogen levels and theta, with 7-day period, in (b), show a relative phase of 0 days. The estrogen levels and alpha2, with 28-day period, in (c), show a relative phase of 10 days. In Figure 3, the progesterone levels and delta, 14-day period, in (a), show a relative phase of 5 days. The progesterone levels and theta, 10-day period, in (b) show a relative phase of 0 days. The progesterone levels and alpha1, 14-day period, in (c), show a relative phase of 2.5 days. The progesterone levels and alpha2, 10-day period and 28-day period, in (d), show a relative phase of 0 days and 5 days respectively. The progesterone levels and beta1, 14-day period, in (e), show a relative phase of 2.5 days. The progesterone levels and beta2, 10-day period, in (f), show a relative phase of 5 days. The phases are directly taken from the figures.

Discussion

We have presented a proposal to analyze long-period biological signals. In the cases studied here, two or more biological signals present nearly coincident frequencies and, consequently, equal periods. We also obtained the relative phases between two or more of these coincident signals. The most relevant frequency components of each signal represent an important part of the total power of this signal. Since each one of these components represent an important part of the whole signal, this behavior seems to indicate the degree in which certain brain signal couple to a specific sex hormone during the 28-day period. In the cases presented here, the brain signals delta and theta seem to accompany the estrogen, i.e., they are in phase all the time, however, alpha2 acts against to the estrogen evolution, i.e., they are out of phase all the time. Theta, alpha1, and beta1 evolve essentially as the progesterone, whereas delta and beta2 go in the contrary sense, and alpha2 follows progesterone, presenting in one case a difference of phase.

It is worth to notice that the reliability of our results was tested by adding some noise to the measured data and then carrying out the same analysis used for the unmodified signals. As expected, the analysis of the altered signals comes to different conclusions, pointing to the robustness and reliability of the method under consideration in this paper. Also notice that the method followed here is more deterministic than statistical: once the averaged signals are given, the conclusions can be drawn directly.

The reliability of repeated EEG recordings has been previously studied along 28 days [28] and considering also nine months [7]. The EEG data utilized in the present study [17] was taken and processed under similar conditions. Our analysis for long time periodic signals (28 days), as mentioned, is based on these reliable data of hormone levels and brain signals. Since brain signals have been associated with functional abilities and difficulties for certain brain processes [2931] as well as cognitive processes, particularly, associated with women sex hormones [13, 14, 3234], the method presented here could be used to give additional information on the importance of sex hormones and brain signals on several functional abilities and cognitive processes. A careful study of the meaning of our results, in connection with these aspects, is beyond the scope of our present work.

Limitations

Before ending, we would like to emphasize some assumptions and procedures to which our analysis is circumscribed.

With respect to the Fourier analysis: 1) We have assumed that both the brain signals and sex hormones levels are described by functions that repeat themselves every 28 days. This is an important assumption but usually taken by granted, 2) the method requires a pre-analysis to determine if the real part or the imaginary part (or both) of the Fourie-transformed signal should be taken into account, 3) the method requires a criterion to decide which peak of the Fourie-transformed signal corresponds to a relevant part of the signal of interest.

Concerning the measured data: 1) The EEG data utilized in this study was obtained from subjects at rest and with open eyes. This is a particular condition to take into account, 2) finally, the EEG data correspond to subjects with ages in certain range.

Conclusion

Our findings suggest that the proposed method achieved in the dataset is useful to find a characteristic profile of coincident frequencies and relative phases among brain activity and hormonal signals. The approach and kind of analysis we introduce here could be generalized to study other biological signals and to search for possible similar patterns and relations to those shown in this work. The knowledge of the typical periods and phases among biological signals of clinical interest seems to provide us with an auxiliary tool of analysis of possible diagnostic use.

Declarations

Acknowledgements

We thank M. Corsi-Cabrera for comments during the development of this work. Partially supported by CONCYTEG grants 06-16-K117-142, and 06-16-K117-99, CONACYT grants 51306-F, 39941, 60645, and 52365 and Universidad de Guanajuato.

Authors’ Affiliations

(1)
Departamento de Ciencias Médicas, División de Ciencias de la Salud, Campus León, Universidad de Guanajuato
(2)
Facultad de Ciencias Físico-Matemáticas, Universidad Michoacana
(3)
Departamento de Física, División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato

References

  1. Guler I, Hardalac F, Muldur S: Determination of aorta failure with the application of FFT, AR and wavelet methods to Doppler technique. Comput Biol Med. 2001, 31: 229-238.View ArticlePubMedGoogle Scholar
  2. Harris CM: The Fourier analysis of biological transients. J Neurosc Methods. 1998, 83 (1): 15-34.View ArticleGoogle Scholar
  3. Lin Z, Chen JD: Advances in time-frequency analysis of biomedical signals. Crit Rev Biomed Eng. 1996, 24: 1-72.View ArticlePubMedGoogle Scholar
  4. Thompson WJ, Macdonald JR: Discrete and integral Fourier transforms: analytical examples. Proc Natl Acad Sci USA. 1993, 90: 6904-6908.PubMed CentralView ArticlePubMedGoogle Scholar
  5. Rippon G: Electroencephalography. Methods in Mind. Edited by: Senior C, Russel T, Gazzaniga MS. 2006, Cambridge MA: MIT Press, 237-262.Google Scholar
  6. Fisch BJ: Descriptors of EEG activity. Basic Principles of Digital and Analog EEG. Edited by: Fisch BJ. 1999, The Netherlands: Elsevier Science, 145-154. 3Google Scholar
  7. Corsi-Cabrera M, Galindo-Vilchis L, Del Río Portilla Y, Arce C, Ramos-Loyo J: Within-subjects reliability and inter-session stability of EEG power and coherent activity in women evaluated monthly over nine months. Clin Neurophysiol. 2007, 118: 9-21.View ArticlePubMedGoogle Scholar
  8. Ahmed SN, Unterman J, McCloskey B, Sinclair DB: Clinical significance of photic stimulation during routine EEGs of adult patients. Am J Electroneurodiagnostic Technol. 2006, 46: 356-362.PubMedGoogle Scholar
  9. Zani A, Mado A: The Cognitive Electrophysiology of Mind and Brain. Edited by: Zani A, Mado A. 2003, San Diego CA USA: Academic Press, 3-12.View ArticleGoogle Scholar
  10. Niedermeyer E: The normal EEG of the waking adult. Electroencephalography. Basic principles, clinical applications and related fields. Edited by: Nierdermeyer E, Lopes Da Silva F. 1993, Maryland USA: Williams & Wilkins, 131-152. 3Google Scholar
  11. Steriade M, Gloor P, Llinas RR, Lopes da Silva FH, Mesulam MM: Basic mechanisms of cerebral rhythmic activities. Electroencephalogr Clin Neurophysiol. 1990, 76: 481-508.View ArticlePubMedGoogle Scholar
  12. Majewska MD: Steroids and brain activity, essential dialogue between body and mind. Biochem Pharmacol. 1987, 36 (22): 3781-3788.View ArticlePubMedGoogle Scholar
  13. Walpurger V, Pietrowsky R, Kirschbaum C, Wolf OT: Effects of the menstrual cycle on auditory event-related potentials. Horm Behav. 2004, 46: 600-606.View ArticlePubMedGoogle Scholar
  14. Solís-Ortiz S, Guevara MA, Corsi-Cabrera M: Performance in a test demanding prefrontal functions is favored by early luteal phase progesterone. An electroencephalographic study. Psychoneuroendocrinology. 2004, 29: 1047-1057.View ArticlePubMedGoogle Scholar
  15. Smith MJ, Adams LF, Schmith PJ, Rubinow DR, Wassermann EM: Effects of ovarian hormones on human cortical excitability. Ann Neurol. 2002, 51: 599-603.View ArticlePubMedGoogle Scholar
  16. Krug R, Plihal W, Fehm HL, Born J: Selective influence of the menstrual cycle on perception of stimuli with reproductive significance: an event-related potential study. Psychophysiology. 2000, 37: 111-122.View ArticlePubMedGoogle Scholar
  17. Solís-Ortiz S, Corsi-Cabrera M, Ramos J, Guevara MA, Arce C: EEG oscillations during menstrual cycle. Int J Neurosci. 1994, 76: 279-292.View ArticlePubMedGoogle Scholar
  18. Fleck KM, Polich J: P300 and the menstrual cycle. Electroencephalogr Clin Neurophysiol. 1988, 71: 157-160.View ArticlePubMedGoogle Scholar
  19. Becker D, Creutzfeldt OD, Schwibbe M, Wuttke W: Changes in physiological, EEG and psychological parameters in women during the spontaneous menstrual cycle and following oral contraceptives. Psychoneuroendocrinology. 1982, 7: 75-90.View ArticlePubMedGoogle Scholar
  20. Campos RG: A quadrature formula for the Hankel transform. Num Algorithms. 1995, 9: 343-354.View ArticleGoogle Scholar
  21. Campos RG, Juárez LZ: A discretization of the continuous Fourier transform. Nuovo Cimento B. 1992, 100: 703-711.View ArticleGoogle Scholar
  22. Solís-Ortiz S, Corsi-Cabrera M: EEG pattern of anxiety along the menstrual cycle. Rev Mex Psicol. 2002, 12: 187-195.Google Scholar
  23. Marshall JC: Hormonal regulation of the menstrual cycle and mechanisms of ovulation. Endocrinology. Edited by: Leslie J, De Groot J, Jamenson L. 2001, Philadelphia: W.B. Saunders Company, 2073-2085.Google Scholar
  24. Guevera MA, Ramos J, Corsi-Cabrera M: CAPTUSEN: un sistema para la adquisición computarizada del EEG y potenciales relacionados a eventos. Rev Mex Psicol. 2000, 17: 77-88.Google Scholar
  25. Guevera MA, Ramos J, Hernández-González M, Zarabozo D, Corsi-Cabrera M: POTENCOR; a program to calculate power and correlation spectra of EEG signals. Comput Meth Prog Biomed. 2002, 72: 241-250.View ArticleGoogle Scholar
  26. Campos RG, Domínguez Mota F, Coronado E: Quadrature formulas for integrals transforms generated by orthogonal polynomials. arXiv:0805.2111v1. 2008, http://www.arXiv.orgGoogle Scholar
  27. Camarena-Ibarrola A, Chavez E: Using a new discretization of the Fourier transform to discriminate voiced from unvoiced speech. Proceedings of the Seventh Mexican International Conference on Computer Science. 2006, 127-134. doi: 10.1109/ENC.2006.36Google Scholar
  28. Corsi-Cabrera M, Solís-Ortiz S, Guevara MA: Stability of EEG inter and in-trahemisferic correlation in women. Electroenceph Clin Neurophysiol. 1997, 102: 248-255.View ArticlePubMedGoogle Scholar
  29. Steriade M: Grouping of brain rhythms in corticothalamic systems. Neuroscience. 2006, 137 (4): 1087-1106.View ArticlePubMedGoogle Scholar
  30. Klimesch W, Schack B, Sauseng P: The functional significance of theta and upper alpha oscillations. Exp Psychol. 2005, 52: 99-108.View ArticlePubMedGoogle Scholar
  31. Harmony T, Fernandez T, Silva J, Bernal J, Díaz-Comas L, Reyes A, Marosi E, Rodríguez M: EEG delta activity: An indicator of attention to internal processing during performance of mental tasks. Int J Psychophysiol. 1996, 24: 161-171.View ArticlePubMedGoogle Scholar
  32. Solís-Ortiz S, Corsi-Cabrera M: Sustained attention is favored by progesterone during eraly luteal phase and visuo-spatial memory by estrogens during ovulatory phase in young women. Psychoneuroendocrinology. 2008, 33: 989-998.View ArticlePubMedGoogle Scholar
  33. Hollander A, Hausmann M, Hamm JP, Corballis MC: Sex hormonal modulation of hemispheric asymmetries in the attentional blink. J Int Neuropsychol Soc. 2005, 11: 263-272.View ArticlePubMedGoogle Scholar
  34. Fernandez G, Weis S, Stoffel-Wagner B, Tendolkar I, Reuber M, Beyenburg S, Klaver P, Fell J, de Greiff A, Ruhlmann J, Reul J, Elger CE: Menstrual cycle dependent neural plasticity in the adult human brain is hormone, task, and region specific. J Neurosci. 2003, 23: 3790-3795.PubMedGoogle Scholar

Copyright

© Solís-Ortíz et al; licensee BioMed Central Ltd. 2009

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.