Raleighplots, or Rayleigh plots (also called circlegrams [1] and closely related to circular histograms, [2] phasor diagrams, [3] and wind roses [4] ), are statistical graphics that serve as graphical representations for a Raleigh test that map a mean vector to a circular plot. Raleigh plots have many applications in the field of chronobiology, such as in studying butterfly migration patterns or protein and gene expression, and in other fields such as geology, cognitive psychology, and physics.
Raleigh plots was first introduced by Lord Rayleigh. The concept of Raleigh plots evolved from Raleigh tests, also introduced by Lord Rayleigh in 1880. The Rayleigh test is a popular statistical test used to measure the concentration of data points around a circle, identifying any unimodal bias in the distribution. [5] Rayleigh plots emerged from this analysis as a means to illustrate the nature of the distribution.
In a Raleigh plot, each individual is assigned a unit vector with a corresponding angle. These unit vectors are averaged together in a Raleigh plot into the mean vector. The length of the mean vector is determined by r (or R), the mean resultant length. r is a measure of concentration, ranging in value between 0 and 1. If the individual angles of the unit vectors are tightly clustered, then the r value will be closer to 1 while if they are widely scattered, then r will be closer to zero. [6] The mean vector is positioned in the center of a circle. Dashes along the circumference of this circle denote desired values. Examples include angles from magnetic north (zero degrees) going clockwise (e.g., 90 degrees from magnetic north, or eastward); times of day, which can also be described in zeitgeber time and circadian time; and phase. Dots on the circumference are usually used to indicate individual unit vectors and their respective angle in regard to the values being measured. Raleigh plots can also use more than one mean vector, particularly if one wants to display the mean vector for different tested groups in the study or to compare mean vectors between groups. [7]
The example to the left is a Raleigh plot that has a high r value. Blue and yellow dots indicate individuals from different groups being tested, and the position of the dots indicate in which angle from magnetic north each tested individual is traveling. Due to the high overall concentration of individuals going at an angle between zero degrees and ninety degrees, the mean vector is much longer. Compare the figure to the left with the figure below:
In this second example (to the right), the Raleigh plot has a low r value. Both yellow and blue dots are spread along the circumference of the circle, indicating many individuals are traveling at different angles. The largest cluster of individuals, a group traveling between 180 degrees and 270 degrees, causes the mean vector to be pointed at an angle in that direction. Notably, due to the variability in the direction within this group, the mean vector is much smaller.
Raleigh plots have been used in chronobiology studies on the biological clocks behind monarch butterfly migration patterns. They are particularly relevant for studying sun compass orientation in migrating butterflies.
In butterfly migration studies, the Raleigh plot maps the orientation of the butterfly when allowed to fly, where the circumference is marked as a compass, with north (N) at the top position. Given the plotted data points, a mean r vector is drawn to indicate the mean orientation of the butterflies in a particular condition. [8] [9]
In his studies on the neurobiology of butterfly migration, Steven M. Reppert observes the oriented flight behavior of monarch butterflies. Reppert explains how Raleigh plots are used to handle butterfly orientation data and as tools for the data analysis. [10]
In a 2012 study by Reppert and colleagues on the sufficiency of an antenna for proper time compensation and sun compass orientation in the monarch butterfly, Raleigh plots were used to present the mean flight orientation of butterflies subjected to different study conditions. Along the edge of the circle, degrees 0º (magnetic north) through 360º are shown, the orientation of each butterfly is marked with a dot, and a mean vector is drawn to represent the mean flight orientation recorded. [7] In a 2016 annual review on the neurobiology of monarch butterfly migration, Reppert, Guerra, and Merlin also use Raleigh plots to present butterfly orientation data. The plots were used in the study of the time-compensated sun compass in monarch butterflies. [10] A 2018 review by Reppert and de Roode on the mechanisms of monarch butterfly migration also used Raleigh plots, or circlegrams, to represent butterfly orientation data. Each dot indicates the orientation in which a butterfly individual flew continuously for 5 minutes or longer, and the vector points in the mean direction (some degrees from north) with a magnitude proportional to the mean orientation. [1]
Raleigh plots can be used to visualize circadian rhythms in protein or gene expression, and how their phases are affected by other variables or induced conditions.
Jennifer Mohawk, then a postdoc at the University of Texas Southwestern Medical Center, used multiple Raleigh plots to illustrate PER2::LUC expression in her 2013 paper "Methamphetamine and Dopamine Receptor D1 Regulate Entrainment of Murine Circadian Oscillators." [11] Specifically, Mohawk investigated how injections of methamphetamine and D1 antagonist SCH-23390 would shift the peak time of PER2 expression in the liver, lung, pituitary gland, and salivary gland. In these plots, the Raleigh plot can be interpreted as a 24 hour clock with CT 0 at the top of the circle and CT 12 at the bottom of the circle. Each arrow represents the average peak phase of PER2::LUC expression of each group. The strength of the phase clustering is symbolized by the length of the arrow, meaning stronger clustering or closer data points resulting in longer arrows. The individual data points are plotted on the outside of the circle and their unique color and shape resemble the different groups of conditions. At ZT 7 is a pink box that shows the timing of the methamphetamine injection. Mohawk and collaborators compared the angle of the vector, or the mean phase, between the different groups in order to determine if methamphetamine injections induced statistically significant phase changes of PER2::LUC expression within the different glands and organs. [11]
Similarly, Tsedey Mekbib, then a PhD student at the Morehouse School of Medicine, utilized Raleigh plots to depict how the knockdown of SIAH2 impacted the rhythmic expression of all other genes in both males and females. After profiling the entire transcriptome via RNA sequence on the liver at frequent time intervals, the expression of peak timing for all rhythmically expressed genes was plotted on Raleigh plots for each group. These Raleigh plots contain vectors that represent the average peak phase. However, instead of differing the length of the vectors to illustrate the variability in the data points, Mekbib and the other collaborators added a +/- 95% confidence interval that is represented by red range along the circle. In addition to the typical Raleigh plot, the left half of the circle is shaded darker to better visualize the night phase occurring between CT 12 and CT 24. [12]
Raleigh plots can be used to visualize the circadian rhythms in protein or gene expression in the suprachiasmatic nucleus (SCN) specifically, and how they might influence other peripheral tissues.
Raleigh plots were used by Elizabeth Maywood, an English researcher at the MRC Laboratory of Molecular Biology, to visualize how pacemaking activity and synchrony between host SCN cells lacking vasoactive intestinal peptide (VIP) can be restored with a wild-type SCN graft. These plots show vectors that represent the phase of the host SCN cells, measured by PER2::LUC expression. Each plot has a value representing the mean vector length with time points where the cell phases are closer in sync having a value closer to 1. Maywood and collaborators showed that VIP-null host SCN cells synchrony deteriorated over time based on the mean vector length of the Raleigh plots decreasing, and concluded that paracrine signaling from an introduced wild-type SCN graft is sufficient to restore the synchrony between SCN cell pacemakers based on the mean vector length increasing after the graft was introduced. [13]
Mariko Izumo, while working at the University of Texas Southwestern Medical Center, used Raleigh plots to assess the effect of knocking out BMAL1 in the SCN on the circadian expression of PER2, measured using PER2::LUC expression, in the SCN and in peripheral tissues. Izumo and others found that knocking out BMAL1 led to the desynchronization and dampening of PER2 expression in peripheral tissues with Raleigh plots showing different mean phases of the rhythm of PER2 expression. They also show that light/dark cycles and feeding can restore synchrony in peripheral tissues with Raleigh plots showing the same mean phases of the rhythm of PER2 expression. [14]
Similarly, Yongli Shan, also while working at the University of Texas Southwestern Medical Center, used Raleigh plots to show that BMAL1 knockout AVP SCN neurons and VIP SCN neurons show a loss in rhythm in expression of PER2. BMAL1 knockout neurons have data points around the circular diagram and a small mean vector r, while wild-type neurons have a larger mean vector r and data points closer together. Additionally, Shan and others show that intercellular connections to the rest of the SCN was sufficient in restoring rhythmicity in BMAL1 knockout VIP neurons but not AVP neurons with an increase in the mean vector r in VIP neurons but not in AVP neurons. [15]
Raleigh plots or variations on Raleigh plots are used in fields beyond chronobiology. While Raleigh plots visualize a mean vector for data, the variational plots that are closely related to Raleigh plots may visualize histogram data in spokes on a circular chart.
In geology, circular histogram plots or rose diagrams can be used to characterize tectonic plate movements. For example, they may be used to visualize the frequency and direction of fault line motion. [16] In meteorology and climate studies, wind roses are used to present data on the direction, duration, and speed of winds that occurred at a given location. For example, in wind roses released by the Midwestern Regional Climate Center, the length of a spoke in a particular direction, representing a histogram bin, is proportional to the duration of time for which the wind was blowing in that direction, with different colors to show wind speed categories. [4]
Additionally, Raleigh plots can be used in cognitive psychology. Joëlle Provasi, professor at École pratique des hautes études, used them to explain the response of children with or without a lesion in their cerebellum due to surgically removed tumor to a rhythmic stimulus. Provasi and others show a Raleigh plot depicting responses that are close to the stimulus as data points around the top at 0 degrees with a mean vector value close to 1, and a Raleigh plot depicting responses that are irregular with dots spread around the circular plot and a mean vector value closer to 0. [17]
The suprachiasmatic nucleus or nuclei (SCN) is a small region of the brain in the hypothalamus, situated directly above the optic chiasm. The SCN is the principal circadian pacemaker in mammals, responsible for generating circadian rhythms. Reception of light inputs from photosensitive retinal ganglion cells allow the SCN to coordinate the subordinate cellular clocks of the body and entrain to the environment. The neuronal and hormonal activities it generates regulate many different body functions in an approximately 24-hour cycle.
Cryptochromes are a class of flavoproteins found in plants and animals that are sensitive to blue light. They are involved in the circadian rhythms and the sensing of magnetic fields in a number of species. The name cryptochrome was proposed as a portmanteau combining the chromatic nature of the photoreceptor, and the cryptogamic organisms on which many blue-light studies were carried out.
CLOCK is a gene encoding a basic helix-loop-helix-PAS transcription factor that is known to affect both the persistence and period of circadian rhythms.
Period (per) is a gene located on the X chromosome of Drosophila melanogaster. Oscillations in levels of both per transcript and its corresponding protein PER have a period of approximately 24 hours and together play a central role in the molecular mechanism of the Drosophila biological clock driving circadian rhythms in eclosion and locomotor activity. Mutations in the per gene can shorten (perS), lengthen (perL), and even abolish (per0) the period of the circadian rhythm.
Neuronal PAS domain protein 2 (NPAS2) also known as member of PAS protein 4 (MOP4) is a transcription factor protein that in humans is encoded by the NPAS2 gene. NPAS2 is paralogous to CLOCK, and both are key proteins involved in the maintenance of circadian rhythms in mammals. In the brain, NPAS2 functions as a generator and maintainer of mammalian circadian rhythms. More specifically, NPAS2 is an activator of transcription and translation of core clock and clock-controlled genes through its role in a negative feedback loop in the suprachiasmatic nucleus (SCN), the brain region responsible for the control of circadian rhythms.
The PER3 gene encodes the period circadian protein homolog 3 protein in humans. PER3 is a paralog to the PER1 and PER2 genes. It is a circadian gene associated with delayed sleep phase syndrome in humans.
PER2 is a protein in mammals encoded by the PER2 gene. PER2 is noted for its major role in circadian rhythms.
Aryl hydrocarbon receptor nuclear translocator-like 2, also known as Arntl2, Mop9, Bmal2, or Clif, is a gene.
Period circadian protein homolog 1 is a protein in humans that is encoded by the PER1 gene.
Basic helix-loop-helix ARNT-like protein 1 or aryl hydrocarbon receptor nuclear translocator-like protein 1 (ARNTL), or brain and muscle ARNT-like 1 is a protein that in humans is encoded by the BMAL1 gene on chromosome 11, region p15.3. It's also known as MOP3, and, less commonly, bHLHe5, BMAL, BMAL1C, JAP3, PASD3, and TIC.
In molecular biology, an oscillating gene is a gene that is expressed in a rhythmic pattern or in periodic cycles. Oscillating genes are usually circadian and can be identified by periodic changes in the state of an organism. Circadian rhythms, controlled by oscillating genes, have a period of approximately 24 hours. For example, plant leaves opening and closing at different times of the day or the sleep-wake schedule of animals can all include circadian rhythms. Other periods are also possible, such as 29.5 days resulting from circalunar rhythms or 12.4 hours resulting from circatidal rhythms. Oscillating genes include both core clock component genes and output genes. A core clock component gene is a gene necessary for to the pacemaker. However, an output oscillating gene, such as the AVP gene, is rhythmic but not necessary to the pacemaker.
Steven M. Reppert is an American neuroscientist known for his contributions to the fields of chronobiology and neuroethology. His research has focused primarily on the physiological, cellular, and molecular basis of circadian rhythms in mammals and more recently on the navigational mechanisms of migratory monarch butterflies. He was the Higgins Family Professor of Neuroscience at the University of Massachusetts Medical School from 2001 to 2017, and from 2001 to 2013 was the founding chair of the Department of Neurobiology. Reppert stepped down as chair in 2014. He is currently distinguished professor emeritus of neurobiology.
Monarch butterfly migration is the phenomenon, mainly across North America, where the subspecies Danaus plexippus plexippus migrates each summer and autumn to and from overwintering sites on the West Coast of California or mountainous sites in Central Mexico. Other subspecies perform minor migrations or none at all. This massive movement of butterflies has been called "one of the most spectacular natural phenomena in the world".
Ueli Schibler is a Swiss biologist, chronobiologist and a professor at the University of Geneva. His research has contributed significantly to the field of chronobiology and the understanding of circadian clocks in the body. Several of his studies have demonstrated strong evidence for the existence of robust, self-sustaining circadian clocks in the peripheral tissues.
Hitoshi Okamura is a Japanese scientist who specializes in chronobiology. He is currently a professor of Systems Biology at Kyoto University Graduate School of Pharmaceutical Sciences and the Research Director of the Japan Science Technology Institute, CREST. Okamura's research group cloned mammalian Period genes, visualized clock oscillation at the single cell level in the central clock of the SCN, and proposed a time-signal neuronal pathway to the adrenal gland. He received a Medal of Honor with Purple Ribbon in 2007 for his research and was awarded Aschoff's Ruler for his work on circadian rhythms in rodents. His lab recently revealed the effects of m6A mRNA methylation on the circadian clock, neuronal communications in jet lag, and the role of dysregulated clocks in salt-induced hypertension.
In the field of chronobiology, the dual circadian oscillator model refers to a model of entrainment initially proposed by Colin Pittendrigh and Serge Daan. The dual oscillator model suggests the presence of two coupled circadian oscillators: E (evening) and M (morning). The E oscillator is responsible for entraining the organism’s evening activity to dusk cues when the daylight fades, while the M oscillator is responsible for entraining the organism’s morning activity to dawn cues, when daylight increases. The E and M oscillators operate in an antiphase relationship. As the timing of the sun's position fluctuates over the course of the year, the oscillators' periods adjust accordingly. Other oscillators, including seasonal oscillators, have been found to work in conjunction with circadian oscillators in order to time different behaviors in organisms such as fruit flies.
Sato Honma is a Japanese chronobiologist who researches the biological mechanisms of circadian rhythms. She mainly collaborates with Ken-Ichi Honma on publications, and both of their primary research focuses are the human circadian clock under temporal isolation and the mammalian suprachiasmatic nucleus (SCN), its components, and associates. Honma is a retired professor at the Hokkaido University School of Medicine in Sapporo, Japan. She received her Ph.D. in physiology from Hokkaido University. She taught physiology at the School of Medicine and then at the Research and Education Center for Brain Science at Hokkaido University. She is currently the director at the Center for Sleep and Circadian Rhythm Disorders at Sapporo Hanazono Hospital and works as a somnologist.
Martha Ulbrick Gillette is a chronobiologist and neurobiologist with research focusing on the effects of circadian clocks on integrative brain functions metabolism and the molecular mechanisms involved in signaling pathways. She is a fellow of the American Association for the Advancement of Science.
Ken-Ichi Honma is a Japanese chronobiologist who researches the biological mechanisms underlying circadian rhythms. After graduating from Hokkaido University School of Medicine, he practiced clinical psychiatry before beginning his research. His recent research efforts are centered around photic and non-photic entrainment, the structure of circadian clocks, and the ontogeny of circadian clocks. He often collaborates with his wife, Sato Honma, in work involving the mammalian suprachiasmatic nucleus (SCN), its components, and associated topics.
Christine Merlin is a French chronobiologist and an associate professor of biology at Texas A&M University. Merlin's research focuses on the underlying genetics of the monarch butterfly circadian clock and explores how circadian rhythms modulate monarch behavior and navigation.