Effect of continuous light on diurnal rhythms in Cyanothece sp. ATCC 51142
© Elvitigala et al; licensee BioMed Central Ltd. 2009
Received: 23 December 2008
Accepted: 15 May 2009
Published: 15 May 2009
Life on earth is strongly affected by alternating day and night cycles. Accordingly, many organisms have evolved an internal timekeeping system with a period of approximately 24 hours. Cyanobacteria are the only known prokaryotes with robust rhythms under control of a central clock. Numerous studies have been conducted to elucidate components of the circadian clock and to identify circadian-controlled genes. However, the complex interactions between endogenous circadian rhythms and external cues are currently not well understood, and a direct and mathematical based comparison between light-mediated and circadian-controlled gene expression is still outstanding. Therefore, we combined and analyzed data from two independent microarray experiments, previously performed under alternating light-dark and continuous light conditions in Cyanothece sp. ATCC 51142, and sought to classify light responsive and circadian controlled genes.
Fourier Score-based methods together with random permutations and False Discovery Rates were used to identify genes with oscillatory expression patterns, and an angular distance based criterion was applied to recognize transient behaviors in gene expression under constant light conditions. Compared to previously reported mathematical approaches, the combination of these methods also facilitated the detection of modified amplitudes and phase-shifts of gene expression. Our analysis showed that the majority of diurnally regulated genes, essentially those genes that are maximally expressed during the middle of the light and dark period, are in fact light responsive. In contrast, most of the circadian controlled genes are up-regulated during the beginning of the dark or subjective dark, and are greatly enriched for genes associated with energy metabolism. Many of the circadian controlled and light responsive genes are found in gene clusters within the Cyanothece sp. ATCC 51142 genome. Interestingly, in addition to cyclic expression patterns with a period of 24 hours, we also found several genes that oscillate with an ultradian period of 12 hours, a novel finding among cyanobacteria.
We demonstrate that a combination of different analytical methods significantly improved the identification of cyclic and transient gene expression in Cyanothece sp. ATCC 51142. Our analyses provide an adaptable and novel analytical tool to study gene expression in a variety of organisms with diurnal, circadian and ultradian behavior.
Temporal periods of light have provided a significant abiotic selection pressure and greatly influenced evolution on earth. Accordingly, diurnal changes in various biological activities have been observed in different organisms, ranging from bacteria to mammals [1–4]. Central for the capacity to anticipate such environmental changes is an internal clock that controls circadian gene expression and orchestrates cell physiology in synchrony with a day and night cycle [5–8].
Cyanobacteria are the only prokaryotes with a circadian clock . These organisms thrive in environments as diverse as fresh and salt waters, glaciers, deserts, and hot springs, where they play a major role in global carbon sequestration and oxygen production. In addition, some cyanobacteria, such as Cyanothece sp. ATCC 51142 (Cyanothece 51142 hereafter) are able to fix atmospheric nitrogen and contribute to a great extent to the marine nitrogen cycle . Cyanothece 51142 is a unicellular cyanobacterium that separates oxygenic photosynthesis and nitrogen fixation temporally, thereby circumventing the irreversible inactivation of the oxygen sensitive nitrogenase enzyme . In Cyanothece 51142, photosynthesis occurs during the day and nitrogen fixation during the night of a diurnal cycle .
Previous work has focused on the identification of different clock components such as kaiA, kaiB and kaiC, and several detailed mathematical models describing the interactions between oscillator components have been published [13–15]. Most of these studies are based on modeling transcriptional and translational feedback loops as well as post-transcriptional modifications of clock genes. In other work, global DNA microarray analysis of Synechocystis sp. PCC 6803, a non-diazotrophic strain, under constant light conditions determined the extent of circadian control at the transcriptional level, leading to the prediction of 2–7% of cyclic genes . In contrast, microarray analyses of Cyanothece 51142 under alternating light and dark conditions revealed that ~30% of the genes in the genome are significantly diurnally regulated , whereas only 10% of these genes are circadian controlled . However, a detailed analytical approach to distinguish between diurnal and circadian regulated genes in Cyanothece 51142 is still outstanding.
Organisms such as Cyanothece 51142 perform different cellular processes at distinct phases of the diurnal cycle and thus depend on tightly regulated expression programs. On the other hand, the capability to adjust to periodic changes in their environment, in particular to alterations of light intensity and quantity, provides an important basis for the ecological success of these organisms. Light is an essential abiotic factor, not only because cyanobacteria obtain their entire energy through photosynthesis, but also for the activation and regulation of many central metabolic processes. In order to understand how Cyanothece 51142 balances gene expression under nitrogen-fixing conditions in consecutive light and dark cycles with altered physiological requirements under constant light conditions, we analyzed data sets from two different microarray experiments, to elucidate details beyond of what has been found in the former studies [16–18].
For this study we developed various advanced mathematical tools to identify rhythmic and transient patterns of gene expression. We demonstrate that a combination of different analytical approaches substantially improves the distinction between transient and cyclic gene expression. Our analysis revealed that a high percentage of the previously observed strongly diurnally regulated genes  are significantly affected under constant light conditions. Interestingly, the majority of circadian controlled genes in Cyanothece 51142 are maximally expressed during the subjective dark and largely comprised of genes related to energy metabolism. Furthermore, we found a significant number of genes that show ultradian rhythms in their expression and oscillate with a period of 12 hours. This study provides a novel and adaptable analytical tool for comparison of data from different microarray experiments to study gene expression in a variety of organisms with diurnal, circadian and ultradian behavior.
Results and Discussion
Microarray data sets
We have analyzed two independent Cyanothece 51142 microarray data sets that were performed over two consecutive diurnal periods with a sampling rate of every four hours and a shift in sampling time of one hour between the experiments. The studies were conducted using Agilent http://www.agilent.com custom-made two-channel microarrays. The work of Stöckel et al.  focused on global gene expression in Cyanothece ATCC 51142 under alternating 12 hour light and dark cycles, while Toepel et al.  investigated changes in gene expression under a 12 hour light and 12 hour dark period followed by a constant light period of 24 hours. In both experiments Cyanothece 51142 cells were grown under nitrogen-fixing conditions.
Data analyses and identification of cyclic expressed genes
In Stöckel et al.  diurnal regulated genes were classified using a correlation network and a fold-change cutoff of 1.3, in which genes with a period other than 24 hours were excluded, since they did not correlate with the main oscillatory network. In Toepel et al. [17, 18] cyclic genes were identified using a 2-fold cutoff for maximal changes in gene expression. Although the main oscillatory behaviors are detectable using these methods, alterations in gene expression such as changes of amplitudes and phase-shifts were not observed, especially because the altered light conditions were applied only for a short period of time. Therefore we proposed a combination of an angular distance based criteria to characterize such alterations in gene expression.
The two main frequencies of 12 and 24 hours were used for further Fourier score calculations. For the determination of cyclic behavior, a Fourier score cutoff of 6.5 corresponding to a False Discovery Rate of 2% under alternating light and dark conditions and of 3% under continuous light was used. Based on this analysis, 659 genes that oscillate with a period of 24 hours and 9 genes that cycle with a period of 12 hours were identified as common in both microarray datasets. An additional 1551 genes that oscillate with a period of 24 hours and 83 with a period of 12 hours were detected only in the Stöckel et.al data set. All of these genes were used for further analyses (see below).
Classification of light responsive genes
Average distance from different light and dark periods in each of the four classified gene groups.
Stöckel et al. data set
Toepel et al. data set
Based on the Fourier score and angular distance based methods, four major groups of genes could be identified. The genes in group 1 continued to oscillate under continuous light and were classified as circadian controlled genes, whereas genes from group 2 were light responsive and showed oscillating expression only under alternating light and dark cycles. Both groups represent genes that cycle with a 24 hour period. In contrast, genes from group 3 and 4 oscillate with an ultradian period of about 12 hours, with genes from group 3 being circadian controlled and the genes from group 4 light responsive (Table 1, [see Additional files 1, 2]).
Taken together, the combination of the angular distance and Fourier Score based methods lent an additional level of confidence to the identification of cyclic expressed genes in Cyanothece 51142. These analyses uncovered that most of the previously identified diurnal genes are indeed light responsive.
Diurnal regulated genes belonging to different functional categories
Considering the central role of photosynthesis in the cellular metabolism of Cyanothece 51142 and other cyanobacteria, a great impact of light on gene expression was anticipated. We found that the expression of the majority of regulatory genes is light responsive, presumably the result of comprehensive changes that take place in order to adjust to altered light conditions [see Additional file 1]. In addition, the expression of many genes related to transcription and translation is light regulated and likely accounts for a higher turnover of proteins under continuous light. Our analyses also revealed that several genes associated with photosynthesis, especially with photosystem II (PSII) such as psbE, psbF, and psbJ are circadian controlled [see Additional file 1]. The PsbJ protein in Synechocystis sp. PCC 6803 has been shown to control the amount of functionally assembled PSII complexes in the thylakoid membrane . Moreover, psbE and psbF which encode the α and β subunits of cytochrome b559 are required to stabilize the reaction center of PSII . Interestingly, a number of subunits of the ATP-synthase complex, including AtpF and AtpH, which are encoded by multicopy genes revealed that one gene copy is light responsive while the other gene copy was found to be circadian controlled [see Additional file 1]. The gene atpH has been identified as clock controlled in previous studies of Synechocystis sp. PCC 6803 .
A more detailed analysis of different metabolic processes revealed that nitrate and sulfate assimilation are largely under circadian control [see Additional file 1]. Furthermore, the genes hoxEFUYH encoding the bidirectional hydrogenase show a significant cyclic behavior in their expression [see Additional files 3ab. In comparison, a recent microarray analysis performed in Cyanothece 51142 over 24 hours under alternating 6 hours light/dark conditions  led the authors to the conclusion that the expression of the hox genes under their conditions are not circadian. However, an alternative interpretation of their data (Figure 3B in ) suggests that the hoxE gene does in fact show an oscillatory expression with a period of 24 hours, which would support our findings. This is also supported by previous studies in Synechococcus sp. PCC 7942 and Synechocystis sp. PCC 6803 which uncovered a strong circadian component in the expression of different hox-genes [5, 22]. In total, our study revealed that 9 out of 17 genes involved in hydrogen metabolism in Cyanothece 51142 were expressed at the same level during the dark and subjective dark and therefore classified as circadian controlled. In contrast, the expression of hupS and hupL encoding subunits of the uptake hydrogenase are severely affected during the subjective dark, with hupS showing an 8-time lower expression level (see Additional files 3C,D). Earlier studies of different Nostoc species uncovered that the expression of hupL is substantially stimulated by supplemented hydrogen [23, 24]. Accordingly, the reduced transcript abundance of hupS and hupL is presumably coupled to a decline in the activity of the nitrogenase enzyme, which has been observed during the subjective dark period under continuous light conditions . Furthermore, even though the genes involved in nitrogen fixation are organized in a single cluster consisting of two adjacent regulons on opposite strands in the genome , different expression profiles for many of those genes were observed under continuous light. The genes nifB, nifS, nifE, nifN, nifX and nifW, which are involved in biogenesis and assembly of the Mo-Fe cofactor, revealed a more than two fold increase in their expression levels  in addition to a shift in phase of about four hours under continuous light conditions [see Additional files 4ab vs. 4cd].
In fact, several light responsive and circadian controlled genes are found in clusters of three or more genes throughout the Cyanothece 51142 genome [see Additional files 1, 5 and 6]. Such clustering of similarly regulated genes into different regions of the genome might provide an opportunity for the cell to control the gene expression more efficiently. On the other hand, additional diurnal and circadian controlled genes can be identified which might not have been detected due to low amplitudes of oscillations. Thus, there is a high probability that an unclassified gene within a group of genes assigned to one of the two categories belongs to the same class.
In this work we compared two independent microarray experiments to study changes in gene expression under alternating and constant light conditions in Cyanothece 51142. Two mathematical approaches were applied to differentiate transient behaviors from cyclic gene expression. Fourier score based methods can successfully be used for identification of cyclic behaviors but they fail to detect some transient behaviors such as shifts of the expression peaks or smaller amplitudes under constant light conditions. Therefore, an angular distance based method was applied which improved the detection of such changes. These methods are valid for distinguishing transient from cyclic behaviors even for the limited number of time points in one of the different light treatments. This approach provides an adaptable tool for future studies of diurnal, circadian and ultradian rhythms in various organisms. Furthermore, the combined analysis of two different microarray data sets enabled the extraction of novel details about diurnal and circadian regulation in Cyanothece 51142 which would not be obtained from the individual experiments alone.
For the Stöckel et al. microarray data set , Cyanothece 51142 cultures were grown in ASP2 medium without NaNO3  at 30°C with air bubbling in alternating 12 hour light and dark cycles with 50 μmol photons/m2*s of white light. After 7 days, 150 mL samples were collected every 4 hours for 2 days, starting with 1 hour into the dark period (time point D1). In total, 12 samples were collected.
For the Toepel et al. data set , Cyanothece 51142 cultures were grown in an airlift bioreactor (6L BioFlo 3000; New Brunswick Scientific, Edison, NJ) in ASP2 medium without NO3 at 30°C in 12 hour light and 12 hour dark cycles. The culture was illuminated by two light-emitting-diode panels using orange (640 nm) and red (680 nm) light, yielding an intensity of 100 μmol photons/m2*s inside the bioreactor. The cultures were grown for 5–6 days under alternating light and dark conditions prior to collecting samples for every 4 hours over a time period of 2 days, starting at the transition from dark to light (time point L0) for 24 hours followed by 24 hours of continuous light.
Algorithms and data processing
Preliminary Data Processing
In previous studies [16–18], the data were processed and analyzed differently, which made a direct comparison of the results difficult. Therefore, in this communication both data sets have been combined and analyzed uniformly. The raw data were normalized using the LOWESS normalization algorithm to avoid a systematic intensity based bias which is commonly observed with two-channel microarrays  and preliminary data processing of both data sets was performed according to .
Identification of Cyclic Genes using Fourier Score
where ω = 2πf.
Initial Fast Fourier Transform calculations indicated existence of oscillations with 24 h and 12 h periods. Therefore Fourier score calculations were performed using those frequencies.
If a given signal consists of a dominant cyclic component of the corresponding frequency, a larger Fourier score is expected. The significance of the Fourier score was quantified by comparing the Fourier score of the original signal with scores of large collection of random signals. These random signals were obtained by using different permutations of the original signal. The original gene expressions were scaled to have a unit standard deviation, which enabled direct comparisons of Fourier scores from different genes.
Identification of transient behaviors using angular distance
where x1 and x2 are the vectors from two different 12 hour periods. D1,2 can have any value between 0 and 2, with 0 representing vectors with the same direction and 2 representing vectors with the opposite direction.
We thank Wenxue Wang and all members in the Pakrasi Lab for collegial discussions. This work is part of a Membrane Biology EMSL Scientific Grand Challenge project at the W. R. Wiley Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the U.S. Department of Energy's Office of Biological and Environmental Research (BER) program located at Pacific Northwest National Laboratory. PNNL is operated or the Department of Energy by Battelle. The project is also partially supported by the National Science Foundation FIBR program under grant number 0425749 as well as the DOE-BES program.
- Iwasaki H, Kondo T: Circadian timing mechanism in the prokaryotic clock system of cyanobacteria. J Biol Rhythms. 2004, 19: 436-444. 10.1177/0748730404269060.View ArticlePubMedGoogle Scholar
- Johnson CH, Mori T, Xu Y: A cyanobacterial circadian clockwork. Curr Biol. 2008, 18: 816-825. 10.1016/j.cub.2008.07.012.View ArticleGoogle Scholar
- Dong G, Golden SS: How a cyanobacterium tells time. Curr Opin Microbiol. 2008, 11: 541-546. 10.1016/j.mib.2008.10.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Harmer SL, Panda S, Kay SA: Molecular bases of circadian rhythms. Annu Rev Cell Dev Biol. 2001, 17: 215-253. 10.1146/annurev.cellbio.17.1.215.View ArticlePubMedGoogle Scholar
- Kucho K, Okamoto K, Tsuchiya Y, Nomura S, Nango M, Kanehisa M, Ishiura M: Global analysis of circadian expression in the cyanobacterium Synechocystis sp. strain PCC 6803. J Bacteriol. 2005, 187: 2190-2199. 10.1128/JB.187.6.2190-2199.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Schaffer R, Landgraf J, Accerbi M, Simon V, Larson M, Wisman E: Microarray analysis of diurnal and circadian-regulated genes in Arabidopsis. Plant Cell. 2001, 13: 113-123. 10.1105/tpc.13.1.113.PubMed CentralView ArticlePubMedGoogle Scholar
- Harmer SL, Hogenesch JB, Straume M, Chang HS, Han B, Zhu T, Wang X, Kreps JA, Kay SA: Orchestrated transcription of key pathways in Arabidopsis by the circadian clock. Science. 2000, 290: 2110-2113. 10.1126/science.290.5499.2110.View ArticlePubMedGoogle Scholar
- Panda S, Antoch MP, Miller BH, Su AI, Schook AB, Straume M, Schultz PG, Kay SA, Takahashi JS, Hogenesch JB: Coordinated transcription of key pathways in the mouse by the circadian clock. Cell. 2002, 109: 307-320. 10.1016/S0092-8674(02)00722-5.View ArticlePubMedGoogle Scholar
- Kondo T, Strayer CA, Kulkarni RD, Taylor W, Ishiura M, Golden SS, Johnson CH: Circadian rhythms in prokaryotes: luciferase as a reporter of circadian gene expression in cyanobacteria. Proc Natl Acad Sci USA. 1993, 90: 5672-5676. 10.1073/pnas.90.12.5672.PubMed CentralView ArticlePubMedGoogle Scholar
- Zehr JP, Waterbury JB, Turner PJ, Montoya JP, Omoregie E, Steward GF, Hansen A, Karl DM: Unicellular cyanobacteria fix N2 in the subtropical North Pacific Ocean. Nature. 2001, 412: 635-638. 10.1038/35088063.View ArticlePubMedGoogle Scholar
- Reddy KJ, Haskell JB, Sherman DM, Sherman LA: Unicellular, aerobic nitrogenfixing cyanobacteria of the genus Cyanothece. J Bacteriol. 1993, 175: 1284-1292.PubMed CentralPubMedGoogle Scholar
- Sherman LA, Meunier P, Colon-Lopez MS: Diurnal rhythms in metabolism: A day in the life of a unicellular, diazotrophic cyanobacterium. Photosynth Res. 1998, 58: 25-42. 10.1023/A:1006137605802.View ArticleGoogle Scholar
- Mehra A, Hong CI, Shi M, Loros JJ, Dunlap JC, Ruoff P: Circadian rhythmicity by autocatalysis. PLoS Computational Biology. 2006, 2: e96-10.1371/journal.pcbi.0020096.PubMed CentralView ArticlePubMedGoogle Scholar
- Miyoshi F, Nakayama Y, Kaizu K, Iwasaki H, Tomita M: A Mathematical Model for the Kai-Protein-Based Chemical Oscillator and Clock Gene Expression Rhythms in Cyanobacteria. J Biol Rhythms. 2007, 22: 69-80. 10.1177/0748730406295749.View ArticlePubMedGoogle Scholar
- Rust MJ, Markson JS, Lane WS, Fisher DS, O'Shea EK: Ordered Phosphorylation Governs Oscillation of a Three-Protein Circadian Clock. Science. 2007, 318: 809-812. 10.1126/science.1148596.PubMed CentralView ArticlePubMedGoogle Scholar
- Stöckel J, Welsh EA, Liberton M, Kunnvakkam R, Aurora R, Pakrasi HB: Global transcriptomic analysis of Cyanothece 51142 reveals robust diurnal oscillation of central metabolic processes. Proc Natl Acad Sci USA. 2008, 105: 6156-6161. 10.1073/pnas.0711068105.PubMed CentralView ArticlePubMedGoogle Scholar
- Toepel J, Welsh E, Summerfield TC, Pakrasi HB, Sherman LA: Differential transcriptional analysis of the cyanobacterium Cyanothece sp. strain ATCC 51142 during light-dark and continuous-light growth. J Bacteriol. 2008, 190: 3904-3913. 10.1128/JB.00206-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Toepel J, McDermott J, Summerfield TC, Sherman LA: Transcriptional analysis of the unicellular diazotrophic cyanobacterium Cyanothece sp. ATCC 51142 grown under short day/night cycles. J Phycology. 2009, DOI: 10.1111/j.1529-8817.2009.00674.x.,Google Scholar
- Duhamel P, Vetterli M: Fast Fourier Transforms: A Tutorial Review and a State of the Art. Signal Processing. 1990, 19: 259-299. 10.1016/0165-1684(90)90158-U.View ArticleGoogle Scholar
- Lind LK, Shukla VK, Nyhus KJ, Pakrasi HB: Genetic and immunological analyses of the cyanobacterium Synechocystis sp. PCC 6803 show that the protein encoded by the psbJ gene regulates the number of photosystem II centers in thylakoid membranes. J Biol Chem. 1993, 268: 1575-1579.PubMedGoogle Scholar
- Pakrasi HB, De Ciechi P, Whitmarsh J: Site directed mutagenesis of the heme axial ligands of cytochrome b559 affects the stability of the photosystem II complex. EMBO J. 1991, 10: 1619-1627.PubMed CentralPubMedGoogle Scholar
- Schmitz O, Boison G, Bothe H: Quantitative analysis of expression of two circadian clock-controlled gene clusters coding for the bidirectional hydrogenase in the cyanobacterium Synechococcus sp. PCC7942. Mol Microbiol. 2001, 41: 1409-1417. 10.1046/j.1365-2958.2001.02612.x.View ArticlePubMedGoogle Scholar
- Tamagnini P, Leitao E, Oliveira P, Ferreira D, Pinto F, Harris DJ, Heidorn T, Lindblad P: Cyanobacterial hydrogenases: diversity, regulation and applications. FEMS Microbiol Rev. 2007, 31: 692-720. 10.1111/j.1574-6976.2007.00085.x.View ArticlePubMedGoogle Scholar
- Axelsson R, Lindblad P: Transcriptional regulation of Nostoc hydrogenases: effects of oxygen, hydrogen, and nickel. Appl Environ Microbiol. 2002, 68: 444-447. 10.1128/AEM.68.1.444-447.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Welsh EA, Liberton M, Stöckel J, Loh T, Elvitigala T, Wang C, Wollam A, Fulton RS, Clifton SW, Jacobs JM, Aurora R, Ghosh BK, Sherman LA, Smith RD, Wilson RK, Pakrasi HB: The genome of Cyanothece 51142, a unicellular diazotrophic cyanobacterium important in the marine nitrogen cycle. Proc Natl Acad Sci USA. 1142, 105: 15094-15099. 10.1073/pnas.0805418105.View ArticleGoogle Scholar
- Lloyd D: Circadian and ultradian clock-controlled rhythms in unicellular microorganisms. Adv Microb Physiol. 1998, 39: 291-338. 10.1016/S0065-2911(08)60019-3.View ArticlePubMedGoogle Scholar
- Min H, Guo H, Xiong J: Rhythmic gene expression in a purple photosynthetic bacterium, Rhodobacter sphaeroides. FEBS Lett. 2005, 579: 808-812. 10.1016/j.febslet.2005.01.003.View ArticlePubMedGoogle Scholar
- Van Praag E, Degli Agosti R, Bachofen R: Rhythmic activity of uptake hydrogenase in the prokaryote Rhodospirillum rubrum. J Biol Rhythms. 2000, 15: 218-224. 10.1177/074873040001500303.View ArticlePubMedGoogle Scholar
- Kippert F, Hunt P: Ultradian clocks in eukaryotic microbes: from behavioural observation to functional genomics. Bioessays. 2000, 22: 16-22. 10.1002/(SICI)1521-1878(200001)22:1<16::AID-BIES5>3.0.CO;2-1.View ArticlePubMedGoogle Scholar
- Schibler U, Naef F: Cellular oscillators: rhythmic gene expression and metabolism. Curr Opin Cell Biol. 2005, 17: 223-229. 10.1016/j.ceb.2005.01.007.View ArticlePubMedGoogle Scholar
- Quackenbush J: Microarray data normalization and transformation. Nature Genetics. 2002, 32: 496-501. 10.1038/ng1032.View ArticlePubMedGoogle Scholar
- de Lichtenberg U, Jensen J, Fausbøll A, Jensen TS, Bork P, Brunak S: Comparison of computational methods for the identification of cell cycle-regulated genes. Bioinformatics. 2005, 21: 1164-1171. 10.1093/bioinformatics/bti093.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B. 1995, 57: 289-300.Google Scholar