Abstract

Knowledge about biological rhythms of diseases may not only help in understanding the pathophysiology of diseases but can also help health service policy makers and emergency department directors to allocate resources efficiently. Aneurysmal subarachnoid haemorrhage (SAH) has high rates of morbidity and mortality. The incidence of SAH has been attributed to patient-related factors such as characteristics of aneurysms, smoking, and hypertension. There are studies showing that the incidence of aneurysmal SAH appears to behave in periodic fashions over long time periods. However, there are inconsistencies in the literature regarding the impact of chronobiological factors such as circadian, seasonal, and lunar cycle factors on the occurrence of SAH. In this study, we focused on the analysis of a temporal pattern of SAH (infradian rhythms) with a novel approach using circular statistical methods. We aimed to see whether there is a circular pattern for the occurrence of SAH at all and if so, whether it can be related to known temporal patterns based on available literature. Our study did not support the notion that aneurysmal subarachnoid haemorrhages occur on any specific day in a cycle with specific lengths up to 365 days including specific weekdays, full moon, equinoxes, and solstices. Hence, we found no relationship between SAH incidence and timing. Study in larger populations using similar circular statistical methods is suggested.

1. Introduction

Biological rhythms are cyclic phenomena that help living organisms to adapt to different periodicities like solar or lunar rhythms. The biological rhythms are classified into (a) circadian (from the Latin circa dies, period of ∼24 hours), (b) ultradian (period <24 hours), and (c) infradian (period >24 hours) rhythms [1]. Several studies have pointed to health risks associated with desynchronized biological rhythms, including cancer, heart disease, digestive disorders, and menstrual irregularities [2]. Knowledge about biological rhythms of diseases may not only help in understanding the pathophysiology of diseases but can also help health service policy makers and emergency department directors to allocate resources efficiently.

Among all cerebrovascular pathologies, aneurysmal subarachnoid haemorrhage (SAH) represents 1–7% of all strokes [3] and has high rates of morbidity and mortality [4]. The incidence of SAH due to rupture of aneurysms may be attributed to genetic components, but it is likely to be largely affected by potential modifiable and/or predictable extrinsic factors [5].

Characteristics of the aneurysm (i.e., size [68], morphology [6, 9, 10], location [6, 10], and haemodynamic features [1113]) have been consistently shown to be related to the risks of aneurysm rupture and subarachnoid haemorrhage. Studies have also supported patient-related factors such as smoking, hypertension [14, 15], alcohol [14], and physical activity [15, 16] as contributors to aneurysm rupture risk. There are studies showing that the incidence of aneurysmal SAH appears to behave in periodic fashions over long time periods [17]. However, there are inconsistencies in the literature regarding the impact of circadian [15, 18, 19], seasonal [1830], and lunar cycle [3134] on the occurrence of SAH.

All previous studies have investigated specific periodicities, for example, weekly, monthly, seasonal, or annual cycles. In this study, we focused on the analysis of a temporal pattern (infradian rhythms—cycles longer than one day) of SAH with a novel approach using circular statistical methods. Firstly, we aim to see whether there is a circular pattern for the occurrence of SAH at all and if so, whether it can be related to known temporal patterns based on available literature.

2. Methods

We extracted dates of admission under the diagnosis of SAH to Royal Hobart Hospital in Tasmania (Australia) from 8/1/1986 to 9/5/2003 from previously collected data. The dates of incidences of SAH had been extracted from medical records, following the patient privacy and data protection regulations. Ethics approval was not required for this study as the study did not involve human participation or use of any patient’s personal data.

Firstly, we aimed to see whether there is a circular pattern (with cycle lengths from 1 to 365 days) for the occurrence of SAH at all. The statistical methods that deal with circular/directional data are referred to as circular statistical methods and require input data in angular format (degrees or radians). Dates of SAH incidences were converted to serial numbers (January 1, 1900, as serial number 1) using Microsoft Excel functions. The serial numbers (linear values) were then transformed to angular values measured on a scale of 0–359 degrees for each cycle length from 1 to 365 (days of a year). For example, for a four-day cycle, any serial number modulus 4 which equaled to two was transformed to 180 degrees. In order to investigate whether there is a circular pattern in our data, we needed statistical procedures to test the null hypothesis that points are spread uniformly around the circle without a preferred direction. Raleigh and Kuiper statistical tests are common circular statistical tests for uniformity. The readers are referred to related literature for further technical information about these tests [35]. Circular statistics analyses including Raleigh and Kuiper tests were performed using circstat module [36] in STATA (StataCorp. 2021. Stata Statistical Software: Release 17. College Station, TX: StataCorp LLC.).

Cycles whose both Raleigh and Kuiper test results were statistically significant () were selected. Raleigh and Kuiper test results with values <0.05 are considered statistically significant implying that the incidences of SAH are likely not to be distributed uniformly. Cycle lengths that were multiples of other significant cycle lengths were removed, as if a cycle length test result is statistically significant, all multiples of that cycle length are expected to be statistically significant as well. For example, if cycle length of seven was selected, all cycle lengths that were multiples of seven (14, 21, 28, …) were removed. Among the remaining cycle lengths, we focused on those whose lengths that were within the range of known cycle lengths such as weekly or lunar cycles.

PyEphem astronomy library for Python (version 4.1.4) was used to calculate the lunar phases, equinoxes, and solstices during the period of study [37].

Python codes (version 3.11.1) were used for manipulation of data. STATA MP 17 (StataCorp. 2021. Stata Statistical Software: Release 17. College Station, TX: StataCorp LLC.) was used for descriptive analysis and Pearson chi-square test.

3. Results

There were 453 days in the study period that the incidence of SAH was recorded. Cycle lengths of 4, 7, 31, 53, 83, 86, 165, 167, 225, 226, 227, and 230 had both statistically significant Raleigh and Kuiper tests and thus provided evidence to reject uniformity (Table 1). That simply meant that we could not rule out that occurrence of SAH could happen uniformly in cycles with any of the abovementioned cycle lengths. As the next step, we investigated each of those cycle lengths that have been considered in the literature as biological rhythms for humans [38]. As is seen in Table 1, three known cycles could fall within the 95% confidence interval range of the selected cycle lengths: weekly cycle (7 days), lunar cycle (29.53056 days long (29 d 12 h 44 m)), and 90 and 182 days (period between equinoxes and/or solstices). We included weekly cycle as it appears that our seven-day week, which is found in many ancient and modern civilizations including the three main monotheistic religions, may be an adaptation to an endogenous biologic rhythm rather than the rhythm being a societally impressed phenomenon [39].

Each of the abovementioned cycles was further studied using chi-square tests.

There were no statistically significant differences between the days of the week (Table 2). None of the SAH occurrences exactly coincided with full moon phases, equinoxes, or solstices during the period of study.

Due to the relatively small number of cases during the study period, in addition to the exact dates of full moon phases, equinoxes, and solstices, the number of SAH incidences was calculated in equal periods before and after each of these days. The length of periods was arbitrarily selected as 1/4th of the lunar phases and 1/6th of the times between the equinoxes or solstices. For example, a lunar month was divided into four equal phases.

The number of SAH incidences was not significantly different in any of the abovementioned periods compared with other periods statistically (Tables 35).

The calculated range of lunar phases was 29.27–29.80 days, the range of period between vernal and autumnal equinoxes was 178.8–186.4 days, and that between summer and winter solstices was 181.7–183.5 days.

4. Discussion

Circular analysis is commonly used in biological sciences, in assessment of directional (e.g., animal orientation) data to time-dependent data (e.g., annual cycles and circadian rhythms). The most common statistical exploration of circular data involves testing the null hypothesis that the data are uniformly distributed across all possible values around the circle using the Rayleigh test [40]. There is no consensus in the related published literature regarding the temporal pattern of SAH. For instance, while correlation between seasonal variation and the onset of SAH was limited to certain age groups in one study [18], its correlation was largely inconsistent within and between countries [1829]. Few studies reported higher incidence of SAH in winter/spring [18, 20, 2325], but other studies, even from the same countries, reported no seasonal variation in the SAH occurrence [15, 19, 20, 22, 25, 26, 41].

In addition, a potential association between lunar cycle and incidence of SAH still remains elusive following inconsistent findings both within and between countries [3133].

This generalizes further to stroke where there remain inconsistencies in the impact of chronobiological factors on different stroke types. Studies have shown that there is a relationship between moon phase and stroke, although this varied among classification of stroke [4246]. However, others have found no difference [47].

While there are studies reporting a weekly large peak on Mondays [48], other studies that have investigated a circaseptan variation in SAH onset did not find apparent weekly patterns [18, 24, 49].

Our study did not show any temporal patterns for occurrence of SAH for cycle length of 1 to 365. In other words, this study did not support the idea that SAH occurs more often on any particular day in different cycle lengths up to 365 days.

While inconsistencies between studies could be related to the different sample size, heterogeneous populations, or potential other extrinsic factors (i.e., environmental and geographic conditions), our results could be considered as another evidence that there is no actual temporal pattern (infradian rhythms) per se in occurrence of SAH. Studies with larger sample size using similar circular statistical methods could further clarify the matter.

All previous studies have investigated specific periodicities, for example, weekly, monthly, seasonal, or annual cycles. One of the novel aspects of this study is that we examined all possible cycle lengths between 1 and 365 days. The method used in the study has not been used previously in neurological disorders and can be considered as one of the advantages of this study.

5. Limitations

Our study is retrospective and relies upon hospital records as there is no official registry for all SAH cases in Tasmania. A possible concern about our findings is the reliability of the data source and missing cases of SAH. Royal Hobart Hospital is the only neurosurgical referral center in Tasmania. Additionally, the incidence of SAH in our study (6 per 100000 person-years) has not been statistically different from the reported estimate of incidences in Australia (mean 5.5, range 5.3–6.0 cases per 100,000 person-years). These pieces of information make it less likely that a significant number of SAH episodes occurring in Tasmania could have been excluded from our study.

This study explores only unimodal departure from uniformity; however, there might be more than one cluster around the circle or multimodal patterns of temporal distribution of SAHs.

We used days as unit of possible cycles as the recorded timing of occurrence of SAH (hours and minutes) was not available or reliable and was therefore not explored. Even though we explored arbitrary periods of time around the full moon, equinoxes, and solstices, we mainly tried to find out if there is any specific day in a cycle with specific length, when SAH occurs more frequently. This study has not specifically explored unimodal patterns of SAH occurrences with clusters longer than one day in different length of cycles or cycle lengths longer than 365 days, which could be considered as a limitation of the study.

6. Conclusion

Our study did not support the notion that aneurysmal subarachnoid haemorrhages occur on any specific day in a cycle with specific lengths up to 365 days including specific weekdays, full moon, equinoxes, and solstices. Hence, we found no relationship between SAH incidence and timing. The findings are consistent with some of the related published literature. Study in larger populations using similar circular statistical methods is suggested.

Data Availability

Data are available upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Open-access publishing was facilitated by the University of Sydney, as part of the Wiley-the University of Sydney agreement via the Council of Australian University Librarians.