The Fat Mass and Obesity Associated Gene, FTO, Is Also Associated with Osteoporosis Phenotypes

Abstract

Obesity and osteoporosis are closely correlated genetically. FTO gene has been consistently identified to be associated with obesity phenotypes. A recent study reported that the mice lacking Fto could result in lower bone mineral density (BMD). Thus, we hypothesize that the FTO gene might be also important for osteoporosis phenotypes. To test for such a hypothesis, we performed an association analyses to investigate the relationship between SNPs in FTO and BMD at both hip and spine. A total of 141 SNPs were tested in two independent Chinese populations (818 and 809 unrelated Han subjects, respectively) and a Caucasian population (2,286 unrelated subjects). Combining the two Chinese samples, we identified 6 SNPs in FTO to be significantly associated with hip BMD after multiple testing adjustments, with the combined P values ranged from 4.99×10−4–1.47×10−4. These 6 SNPs are all located at the intron 8 of FTO and in high linkage disequilibrium. Each copy of the minor allele of each SNP was associated with increased hip BMD values with the effect size (beta) of ∼0.025 and ∼0.015 in the Chinese sample 1 and 2, respectively. However, none of these 6 SNPs showed significant association signal in the Caucasian sample, by presenting some extent of ethnic difference. Our findings, together with the prior biological evidence, suggest that the FTO gene might be a new candidate for BMD variation and osteoporosis in Chinese populations.

About the Authors:

Yan Guo

Affiliation: Key Laboratory of Biomedical Information Engineering of Ministry of Education, and Institute of Molecular Genetics, School of Life Science and Technology, Xi'an Jiaotong University, Xi'an, People's Republic of China

Hui Liu

Affiliation: Department of Articular Surgery, Xi'an Red Cross Hospital, Xi'an, People's Republic of China

Tie-Lin Yang

* E-mail: [email protected] (T-LY) (TY); [email protected] (H-WD) (HD)

Affiliation: Key Laboratory of Biomedical Information Engineering of Ministry of Education, and Institute of Molecular Genetics, School of Life Science and Technology, Xi'an Jiaotong University, Xi'an, People's Republic of China

Siyang M. Li

Affiliation: School of Medicine, University of Missouri - Kansas City, Kansas City, Missouri, United States of America

Siyuan K. Li

Affiliation: School of Medicine, University of Missouri - Kansas City, Kansas City, Missouri, United States of America

Qing Tian

Affiliation: School of Public Health and Tropical Medicine, Tulane University, New Orleans, Louisiana, United States of America

Yong-Jun Liu

Affiliation: School of Public Health and Tropical Medicine, Tulane University, New Orleans, Louisiana, United States of America

Hong-Wen Deng

* E-mail: [email protected] (T-LY) (TY); [email protected] (H-WD) (HD)

Affiliations School of Public Health and Tropical Medicine, Tulane University, New Orleans, Louisiana, United States of America, Center of Systematic Biomedical Research, University of Shanghai for Science and Technology, Shanghai, People's Republic of China

Introduction

Osteoporosis is a major public health problem with growing prevalence. Obesity, another common disease, has been demonstrated to be closely related with osteoporosis [1], [2], [3]. Adipocytes and osteoblasts arise from the same progenitor, bone marrow stromal cells, and can transdifferentiate into each other [4]. With aging, the composition of bone marrow shifts to favor the presence of adipocytes, osteoclast activity increases, and osteoblast function declines, leading to osteoporosis. Moreover, both osteoporosis and obesity have high genetic predisposition and the genetic correlation between them have been established across different ethnic groups [1], [3], [5]. Previous candidate gene and bivariate association studies have suggested some common genes both for obesity and osteoporosis, such as RANK [6], SP7 [7], and SOX6 [8].

Recently, several independent large-scale genome-wide association studies (GWAS) consistently identified a gene FTO (fat mass and obesity associated) to be associated with obesity-related traits and obesity risk [9], [10], [11]. The association has been further replicated in multiple studies in different populations [12], [13], [14], [15], resulting in much increasing attention on this gene. The FTO protein contains a double-stranded beta-helix fold homologous to those of Fe(II) and 2-oxoglutarate-dependent oxygenases, which might be involved in DNA demethylase [16]. Experimental animal studies provide direct functional evidence that FTO is a causal gene underlying obesity [17], [18]. Interestingly, a recent study found that the whole body Fto knockout mice displayed immediate postnatal growth retardation with shorter body length, lower body weight, and lower bone mineral density (BMD) than control mice [19]. This study reminded us that FTO might be a common genetic factor influencing not only obesity phenotypes, but also osteoporosis phenotypes, like BMD. To test for such hypothesis, we performed an association analyses to examine the relationship between the FTO gene and BMD. Our study was performed in three sample sets from two ethnicities, including two Chinese Han populations and a Caucasian population, in order to see whether the variants identified are common or ethnicity-specific.

Results

The basic characteristics of the study subjects are presented in Table 1. The major association results for SNPs in FTO with hip BMD are summarized in Table 2. Combining results from these two Chinese samples, 6 SNPs were identified to be significantly associated with hip BMD after multiple testing adjustments by FDR (Table 2). The most significant SNP was rs16952955, with the P values of 8.39×10−4, and 4.31×10−2 in the Chinese sample 1 and Chinese sample 2, respectively. After meta-analysis, the combined P value achieved a significant level of 1.47×10−4. Besides rs16952955, there were 5 additional SNPs (rs2540766, rs2540784, rs16952951, rs12447427, and rs2689247) showed significant association signals with hip BMD, with the combined P values ranged from 4.99×10−4–2.95×10−4 (Table 2). These 6 SNPs have a consistently protective effect on BMD, since each copy of the minor allele of each SNP was associated with increased hip BMD values with the effect size (beta) of ∼0.025 and ∼0.015 in the Chinese sample 1 and 2, respectively. All of these 6 SNPs are located at the intron 8 of FTO. We further characterized the LD for these SNPs using the regional association plot. As shown in Figure 1, these SNPs were in high LD with the top significant SNP rs16952955 (pairwise LD r2>0.9). For the Caucasian sample, however, no significant results were found for these 6 SNPs and the MAF of each SNP was quite different between Caucasian sample and Chinese samples (P<0.01 by χ2 test).

[Figure omitted. See PDF.]

Figure 1. Regional Association Plot for FTO on chromosome 16.

The color scheme of a white-to-black gradient reflects lower to higher linkage disequilibrium (LD, r2). The r2 is calculated between the top significant SNP rs16952955 and other SNPs. The scatter graph indicates the negative logarithm of P value for each SNP, which is based on the combined results in the two Chinese samples.

https://doi.org/10.1371/journal.pone.0027312.g001

[Figure omitted. See PDF.]

Table 1. Basic characteristics of the study subjects.

https://doi.org/10.1371/journal.pone.0027312.t001

[Figure omitted. See PDF.]

Table 2. Significant association results for SNPs in FTO with hip BMD.

https://doi.org/10.1371/journal.pone.0027312.t002

We further performed gender-specific association analyses for the above 6 SNPs with hip BMD in the Chinese samples. The association signals in each gender group were generally weaker than in the total sample, which might be largely due to the smaller sample sizes in each gender group. The associations were mainly driven by male subjects as reflected in the Chinese sample 1, with the P values ranged from 0.012–1.89×10−3 (Table 3).

[Figure omitted. See PDF.]

Table 3. Gender-specific association results for the six SNPs identified for hip BMD in the Chinese samples.

https://doi.org/10.1371/journal.pone.0027312.t003

For the spine BMD, 5 of above 6 SNPs achieved significant level in the Chinese sample 1 (rs16952951: P = 8.78×10−5; rs16952955: P = 1.47×10−4; rs2540784: P = 1.97×10−4; rs2689247: P = 2.11×10−4; rs12447427: P = 4.79×10−4; and rs2540766: P = 1.99×10−3). However, when considering the Chinese samples 1 and 2 together, no SNPs remained significant after FDR adjustments. For the Caucasian sample, we found another 4 SNPs which are also located at the intron 8 of FTO to be nominally associated with spine BMD (rs11076017: P = 5.23×10−3; rs1420318: P = 6.14×10−3; rs1876942: P = 7.21×10−3; and rs17226942: P = 9.23×10−3).

For the previously reported SNP rs9939609 identified by GWAS for obesity phenotypes [9], [10], [11], we only detected nominally significant association with spine BMD in the Chinese sample 1 (P = 0.037).

Discussion

The FTO gene has become the hotspot for researchers since it was reported as a novel obesity-susceptibility gene by a number of GWAS and follow-up replication studies. Given the evidence that obesity and osteoporosis share some common genetic determinations, we first performed an association study examining the SNPs in FTO for association with BMD. We identified a cluster of SNPs in FTO to be significantly associated with hip BMD in Chinese populations.

The human FTO gene is expressed in many tissues including mesenteric fat, pancreas, liver and adipose tissue, with the highest concentrations found in the hypothalamus [10], [20]. Functional studies have demonstrated the direct effect of FTO on obesity. For example, Fischer et al. have reported that the loss of Fto in mice leads to postnatal growth retardation and a significant reduction in adipose tissue and lean body mass [18]. Church et al. have shown that a mutation in the mouse Fto gene results in reduced fat mass, increased energy expenditure, and unchanged physical activity [17]. In respect to FTO's relevance to osteoporosis, a recent study by Gao et al. has found that Fto plays an essential role in postnatal growth. The mice lacking Fto completely displayed postnatal growth retardation manifested as reduced body weight and length, lower BMD [19]. Taking into account of this biological evidence and our statistical findings, we suggest that FTO might have potential role on BMD or osteoporosis.

It is necessary to examine the associations in different populations from different ethnicities, since the genomic variation is greater when compared across ethnicities. Our study identified consistent association for 6 SNPs in intron 8 of FTO with hip BMD in two Chinese populations. Unfortunately, such results were not replicated in the Caucasian population, which implied some extent of ethnic difference. Such ethnic difference could be explained from two aspects. First, the allele frequencies are quite different between the Chinese and Caucasian populations. The minor alleles of these 6 significant SNPs were much common in the Chinese than in the Caucasians, which may contribute to the overall effect. Second, the difference might be age-specific. The two Chinese populations were relatively younger than the Caucasians. Although we have included age as a covariate to adjust BMD, it could not eliminate the potential confounding effect of age on BMD variations thoroughly. Since animal studies have reported that Fto plays an important role in postnatal growth [19], this might suggest that the FTO gene might be a candidate genetic marker for peak bone mass acquisition. However, it is still too early to get the conclusion. We found several other SNPs which are also in intron 8 of FTO for nominally significant association with spine BMD in the Caucasian population, but not in the Chinese populations. Additional studies are needed to investigate the true effect of FTO on BMD in Caucasian populations.

The association of FTO with spine BMD was less significant than with hip BMD in our study, since we only observed significant associations between FTO and spine BMD in the Chinese sample 1. When considering the two Chinese samples together, no SNPs remained significant after multiple testing adjustments. Our results suggest that the effect of FTO on hip BMD might be stronger than spine BMD. However, it is still too soon to get such conclusion. The sample size of the Chinese samples was relatively small, which might decrease the statistical power to detect genetic associations. Increasing the sample size in further studies is needed to validate our results, and we are hopeful that publication of our findings will facilitate replication analyses by other groups.

In summary, our data provide novel evidence that a cluster of SNPs in FTO are associated BMD variations in Chinese populations. Considering the prior biological findings, we suggest that FTO might be a new candidate for osteoporosis. Further studies are warranted to explore the generality of our findings and elucidate the true functional variant.

Materials and Methods

Subjects

The study was approved by the Institutional Review Board or Research Administration of Xi'an Jiaotong University, Hunan Normal University, Creighton University and University of Missouri-Kansas City. Signed informed consent documents were obtained from all study participants before entering the study.

Chinese samples.

The Chinese sample 1 comprised 818 unrelated subjects, which were recruited from Southwest Chinese Han adults living in Changsha city and its neighboring areas. The Chinese sample 2 comprised 809 unrelated subjects drawing from Northwest Chinese Han adults in Xi'an City and its neighboring areas. For all the subjects, the exclusion criteria have been detailed in our earlier publication [21]. Briefly, subjects with chronic diseases and conditions that might potentially affect bone mass, structure, or metabolism were excluded from the study to minimize the influence of known environmental and therapeutic factors on bone variation. BMD at hip and spine was measured using Hologic 4500 W machines (Hologic Inc., Bedford, MA, USA) under the same strict protocols. The coefficient of variation (CV) values of the dual-energy X-ray absorptionmetry (DXA) measurements for spine and hip BMDs were approximately 1.01% and 1.33%, respectively.

Caucasian sample.

The Caucasian sample consisted of 2,286 unrelated adults. All of the subjects were US Caucasians of Northern European origin living in Midwestern area. The exclusion criteria were the same as with Chinese samples. BMD at hip and spine were measured using the same model Hologic 4500 W machines (Hologic Inc., Bedford, MA, USA) under the same strict protocols used in the Chinese sample. The CV values of the DXA measurements for spine and hip BMDs were approximately 1.98% and 1.87%, respectively.

Genotyping

Genomic DNA was extracted from peripheral blood leukocytes using standard protocols. For all the three samples, SNP genotyping was performed using Genome-Wide Human SNP Array 6.0 (Affymetrix, Santa Clara, CA, USA), according to the Affymetrix protocol. Only samples with a minimum call rate of 95% were included. Due to efforts of repeat experiments, all samples met this criteria and the final mean call rate reached a high level of 98.93% for Caucasian sample and 98.96% for Chinese sample. SNPs that deviated from Hardy-Weinberg equilibrium (HWE, P<0.0001) and had a minor allele frequency (MAF)<0.01 were discarded in each sample set. Thus, 141 SNPs in FTO were included for subsequent association analyses. The basic characteristics of these SNPs are summarized in Table S1.

Statistical analyses

Before association analyses, principal component analysis implemented in EIGENSTRAT [22] was used to adjust for potential population stratification that may lead to spurious association results. The first ten principal components emerging from the EIGENSTRAT analyses, along with sex, age, height, weight and BMI, were used as covariates to adjust the raw BMD values in each sample. The residues were used for association analyses. Linear regression implemented in PLINK [23] was used to test for association under the additive inheritance model.

Meta-analysis statistics were generated using METAL software package (http://genome.sph.umich.edu/wiki/METAL_Documentation) taking into account sample size and direction of effect. SNAP was used to characterize linkage disequilibrium (LD) and depict the regional association plot [24]. A raw P value of <0.05 in our study was considered nominally significant, which were further subjected to a false discovery rate (FDR) of Benjamini and Hochberg procedure [25] to account for multiple comparisons.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Properties of FTO SNPs tested in this study

https://doi.org/10.1371/journal.pone.0027312.s001

(DOC)

Author Contributions

Conceived and designed the experiments: YG. Performed the experiments: T-LY HL. Analyzed the data: T-LY QT. Contributed reagents/materials/analysis tools: Y-JL H-WD SML SKL. Wrote the paper: YG.

Citation: Guo Y, Liu H, Yang T-L, Li SM, Li SK, Tian Q, et al. (2011) The Fat Mass and Obesity Associated Gene, FTO, Is Also Associated with Osteoporosis Phenotypes. PLoS ONE6(11): e27312. https://doi.org/10.1371/journal.pone.0027312

References

1. Toth E, Ferenc V, Meszaros S, Csupor E, Horvath C (2005) [Effects of body mass index on bone mineral density in men]. Orv Hetil 146: 1489–1493.E. TothV. FerencS. MeszarosE. CsuporC. Horvath2005[Effects of body mass index on bone mineral density in men].Orv Hetil14614891493

2. Rosen CJ, Bouxsein ML (2006) Mechanisms of disease: is osteoporosis the obesity of bone? Nat Clin Pract Rheumatol 2: 35–43.CJ RosenML Bouxsein2006Mechanisms of disease: is osteoporosis the obesity of bone?Nat Clin Pract Rheumatol23543

3. Zhao LJ, Liu YJ, Liu PY, Hamilton J, Recker RR, et al. (2007) Relationship of obesity with osteoporosis. J Clin Endocrinol Metab 92: 1640–1646.LJ ZhaoYJ LiuPY LiuJ. HamiltonRR Recker2007Relationship of obesity with osteoporosis.J Clin Endocrinol Metab9216401646

4. Gimble JM, Robinson CE, Wu X, Kelly KA (1996) The function of adipocytes in the bone marrow stroma: an update. Bone 19: 421–428.JM GimbleCE RobinsonX. WuKA Kelly1996The function of adipocytes in the bone marrow stroma: an update.Bone19421428

5. Deng FY, Lei SF, Li MX, Jiang C, Dvornyk V, et al. (2006) Genetic determination and correlation of body mass index and bone mineral density at the spine and hip in Chinese Han ethnicity. Osteoporos Int 17: 119–124.FY DengSF LeiMX LiC. JiangV. Dvornyk2006Genetic determination and correlation of body mass index and bone mineral density at the spine and hip in Chinese Han ethnicity.Osteoporos Int17119124

6. Zhao LJ, Guo YF, Xiong DH, Xiao P, Recker RR, et al. (2006) Is a gene important for bone resorption a candidate for obesity? An association and linkage study on the RANK (receptor activator of nuclear factor-kappaB) gene in a large Caucasian sample. Hum Genet 120: 561–570.LJ ZhaoYF GuoDH XiongP. XiaoRR Recker2006Is a gene important for bone resorption a candidate for obesity? An association and linkage study on the RANK (receptor activator of nuclear factor-kappaB) gene in a large Caucasian sample.Hum Genet120561570

7. Zhao J, Bradfield JP, Li M, Zhang H, Mentch FD, et al. (2011) BMD-associated variation at the Osterix locus is correlated with childhood obesity in females. Obesity (Silver Spring) 19: 1311–1314.J. ZhaoJP BradfieldM. LiH. ZhangFD Mentch2011BMD-associated variation at the Osterix locus is correlated with childhood obesity in females.Obesity (Silver Spring)1913111314

8. Liu YZ, Pei YF, Liu JF, Yang F, Guo Y, et al. (2009) Powerful bivariate genome-wide association analyses suggest the SOX6 gene influencing both obesity and osteoporosis phenotypes in males. PLoS One 4: e6827.YZ LiuYF PeiJF LiuF. YangY. Guo2009Powerful bivariate genome-wide association analyses suggest the SOX6 gene influencing both obesity and osteoporosis phenotypes in males.PLoS One4e6827

9. Dina C, Meyre D, Gallina S, Durand E, Korner A, et al. (2007) Variation in FTO contributes to childhood obesity and severe adult obesity. Nat Genet 39: 724–726.C. DinaD. MeyreS. GallinaE. DurandA. Korner2007Variation in FTO contributes to childhood obesity and severe adult obesity.Nat Genet39724726

10. Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, et al. (2007) A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science 316: 889–894.TM FraylingNJ TimpsonMN WeedonE. ZegginiRM Freathy2007A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity.Science316889894

11. Scuteri A, Sanna S, Chen WM, Uda M, Albai G, et al. (2007) Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits. PLoS Genet 3: e115.A. ScuteriS. SannaWM ChenM. UdaG. Albai2007Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits.PLoS Genet3e115

12. Cha SW, Choi SM, Kim KS, Park BL, Kim JR, et al. (2008) Replication of genetic effects of FTO polymorphisms on BMI in a Korean population. Obesity (Silver Spring) 16: 2187–2189.SW ChaSM ChoiKS KimBL ParkJR Kim2008Replication of genetic effects of FTO polymorphisms on BMI in a Korean population.Obesity (Silver Spring)1621872189

13. Chang YC, Liu PH, Lee WJ, Chang TJ, Jiang YD, et al. (2008) Common variation in the fat mass and obesity-associated (FTO) gene confers risk of obesity and modulates BMI in the Chinese population. Diabetes 57: 2245–2252.YC ChangPH LiuWJ LeeTJ ChangYD Jiang2008Common variation in the fat mass and obesity-associated (FTO) gene confers risk of obesity and modulates BMI in the Chinese population.Diabetes5722452252

14. Hotta K, Nakata Y, Matsuo T, Kamohara S, Kotani K, et al. (2008) Variations in the FTO gene are associated with severe obesity in the Japanese. J Hum Genet 53: 546–553.K. HottaY. NakataT. MatsuoS. KamoharaK. Kotani2008Variations in the FTO gene are associated with severe obesity in the Japanese.J Hum Genet53546553

15. Hubacek JA, Bohuslavova R, Kuthanova L, Kubinova R, Peasey A, et al. (2008) The FTO gene and obesity in a large Eastern European population sample: the HAPIEE study. Obesity (Silver Spring) 16: 2764–2766.JA HubacekR. BohuslavovaL. KuthanovaR. KubinovaA. Peasey2008The FTO gene and obesity in a large Eastern European population sample: the HAPIEE study.Obesity (Silver Spring)1627642766

16. Gerken T, Girard CA, Tung YC, Webby CJ, Saudek V, et al. (2007) The obesity-associated FTO gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase. Science 318: 1469–1472.T. GerkenCA GirardYC TungCJ WebbyV. Saudek2007The obesity-associated FTO gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase.Science31814691472

17. Church C, Lee S, Bagg EA, McTaggart JS, Deacon R, et al. (2009) A mouse model for the metabolic effects of the human fat mass and obesity associated FTO gene. PLoS Genet 5: e1000599.C. ChurchS. LeeEA BaggJS McTaggartR. Deacon2009A mouse model for the metabolic effects of the human fat mass and obesity associated FTO gene.PLoS Genet5e1000599

18. Fischer J, Koch L, Emmerling C, Vierkotten J, Peters T, et al. (2009) Inactivation of the Fto gene protects from obesity. Nature 458: 894–898.J. FischerL. KochC. EmmerlingJ. VierkottenT. Peters2009Inactivation of the Fto gene protects from obesity.Nature458894898

19. Gao X, Shin YH, Li M, Wang F, Tong Q, et al. (2010) The fat mass and obesity associated gene FTO functions in the brain to regulate postnatal growth in mice. PLoS One 5: e14005.X. GaoYH ShinM. LiF. WangQ. Tong2010The fat mass and obesity associated gene FTO functions in the brain to regulate postnatal growth in mice.PLoS One5e14005

20. Stratigopoulos G, Padilla SL, LeDuc CA, Watson E, Hattersley AT, et al. (2008) Regulation of Fto/Ftm gene expression in mice and humans. Am J Physiol Regul Integr Comp Physiol 294: R1185–1196.G. StratigopoulosSL PadillaCA LeDucE. WatsonAT Hattersley2008Regulation of Fto/Ftm gene expression in mice and humans.Am J Physiol Regul Integr Comp Physiol294R11851196

21. Guo Y, Tan LJ, Lei SF, Yang TL, Chen XD, et al. (2010) Genome-wide association study identifies ALDH7A1 as a novel susceptibility gene for osteoporosis. PLoS Genet 6: e1000806.Y. GuoLJ TanSF LeiTL YangXD Chen2010Genome-wide association study identifies ALDH7A1 as a novel susceptibility gene for osteoporosis.PLoS Genet6e1000806

22. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, et al. (2006) Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 38: 904–909.AL PriceNJ PattersonRM PlengeME WeinblattNA Shadick2006Principal components analysis corrects for stratification in genome-wide association studies.Nat Genet38904909

23. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.S. PurcellB. NealeK. Todd-BrownL. ThomasMA Ferreira2007PLINK: a tool set for whole-genome association and population-based linkage analyses.Am J Hum Genet81559575

24. Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O'Donnell CJ, et al. (2008) SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics 24: 2938–2939.AD JohnsonRE HandsakerSL PulitMM NizzariCJ O'Donnell2008SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap.Bioinformatics2429382939

25. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I (2001) Controlling the false discovery rate in behavior genetics research. Behav Brain Res 125: 279–284.Y. BenjaminiD. DraiG. ElmerN. KafkafiI. Golani2001Controlling the false discovery rate in behavior genetics research.Behav Brain Res125279284