Human genetic adaptation related to cellular zinc homeostasis.

SLC30A9 encodes a ubiquitously zinc transporter (ZnT9) and has been consistently suggested as a candidate for positive selection in humans. However, no direct adaptive molecular phenotype has been demonstrated. Our results provide evidence for directional selection operating in two major complementary haplotypes in Africa and East Asia. These haplotypes are associated with differential gene expression but also differ in the Met50Val substitution (rs1047626) in ZnT9, which we show is found in homozygosis in the Denisovan genome and displays accompanying signatures suggestive of archaic introgression. Although we found no significant differences in systemic zinc content between individuals with different rs1047626 genotypes, we demonstrate that the expression of the derived isoform (ZnT9 50Val) in HEK293 cells shows a gain of function when compared with the ancestral (ZnT9 50Met) variant. Notably, the ZnT9 50Val variant was found associated with differences in zinc handling by the mitochondria and endoplasmic reticulum, with an impact on mitochondrial metabolism. Given the essential role of the mitochondria in skeletal muscle and since the derived allele at rs1047626 is known to be associated with greater susceptibility to several neuropsychiatric traits, we propose that adaptation to cold may have driven this selection event, while also impacting predisposition to neuropsychiatric disorders in modern humans.

Exploring Adaptive Phenotypes for the Human Calcium-Sensing Receptor Polymorphism R990G.

Rainforest hunter-gatherers from Southeast Asia are characterized by specific morphological features including a particularly dark skin color (D), short stature (S), woolly hair (W), and the presence of steatopygia (S)-fat accumulation localized in the hips (DSWS phenotype). Based on previous evidence in the Andamanese population, we first characterized signatures of adaptive natural selection around the calcium-sensing receptor gene in Southeast Asian rainforest groups presenting the DSWS phenotype and identified the R990G substitution (rs1042636) as a putative adaptive variant for experimental follow-up. Although the calcium-sensing receptor has a critical role in calcium homeostasis by directly regulating the parathyroid hormone secretion, it is expressed in different tissues and has been described to be involved in many biological functions. Previous works have also characterized the R990G substitution as an activating polymorphism of the calcium-sensing receptor associated with hypocalcemia. Therefore, we generated a knock-in mouse for this substitution and investigated organismal phenotypes that could have become adaptive in rainforest hunter-gatherers from Southeast Asia. Interestingly, we found that mouse homozygous for the derived allele show not only lower serum calcium concentration but also greater body weight and fat accumulation, probably because of enhanced preadipocyte differentiation and lipolysis impairment resulting from the calcium-sensing receptor activation mediated by R990G. We speculate that such differential features in humans could have facilitated the survival of hunter-gatherer groups during periods of nutritional stress in the challenging conditions of the Southeast Asian tropical rainforests.

Uncovering Signals of Positive Selection in Peruvian Populations from Three Ecological Regions.

Peru hosts extremely diverse ecosystems which can be broadly classified into the following three major ecoregions: the Pacific desert coast, the Andean highlands, and the Amazon rainforest. Since its initial peopling approximately 12,000 years ago, the populations inhabiting such ecoregions might have differentially adapted to their contrasting environmental pressures. Previous studies have described several candidate genes underlying adaptation to hypobaric hypoxia among Andean highlanders. However, the adaptive genetic diversity of coastal and rainforest populations has been less studied. Here, we gathered genome-wide single-nucleotide polymorphism-array data from 286 Peruvians living across the three ecoregions and analyzed signals of recent positive selection through population differentiation and haplotype-based selection scans. Among highland populations, we identify candidate genes related to cardiovascular function (TLL1, DUSP27, TBX5, PLXNA4, SGCD), to the Hypoxia-Inducible Factor pathway (TGFA, APIP), to skin pigmentation (MITF), as well as to glucose (GLIS3) and glycogen metabolism (PPP1R3C, GANC). In contrast, most signatures of adaptation in coastal and rainforest populations comprise candidate genes related to the immune system (including SIGLEC8, TRIM21, CD44, and ICAM1 in the coast; CBLB and PRDM1 in the rainforest; and BRD2, HLA-DOA, HLA-DPA1 regions in both), possibly as a result of strong pathogen-driven selection. This study identifies candidate genes related to human adaptation to the diverse environments of South America.

The Counteracting Effects of Demography on Functional Genomic Variation: The Roma Paradigm.

Demographic history plays a major role in shaping the distribution of genomic variation. Yet the interaction between different demographic forces and their effects in the genomes is not fully resolved in human populations. Here, we focus on the Roma population, the largest transnational ethnic minority in Europe. They have a South Asian origin and their demographic history is characterized by recent dispersals, multiple founder events, and extensive gene flow from non-Roma groups. Through the analyses of new high-coverage whole exome sequences and genome-wide array data for 89 Iberian Roma individuals together with forward simulations, we show that founder effects have reduced their genetic diversity and proportion of rare variants, gene flow has counteracted the increase in mutational load, runs of homozygosity show ancestry-specific patterns of accumulation of deleterious homozygotes, and selection signals primarily derive from preadmixture adaptation in the Roma population sources. The present study shows how two demographic forces, bottlenecks and admixture, act in opposite directions and have long-term balancing effects on the Roma genomes. Understanding how demography and gene flow shape the genome of an admixed population provides an opportunity to elucidate how genomic variation is modeled in human populations.

The shared genetic architecture of schizophrenia, bipolar disorder and lifespan.

Psychiatric disorders such as Schizophrenia (SCZ) and Bipolar Disorder (BD) represent an evolutionary paradox, as they exhibit strong negative effects on fitness, such as decreased fecundity and early mortality, yet they persist at a worldwide prevalence of approximately 1%. Molecular mechanisms affecting lifespan, which may be widely common among complex diseases with fitness effects, can be studied by the integrated analysis of data from genome-wide association studies (GWAS) of human longevity together with any disease of interest. Here, we report the first of such studies, focusing on the genetic overlap-pleiotropy-between two psychiatric disorders with shortened lifespan, SCZ and BD, and human parental lifespan (PLS) as a surrogate of life expectancy. Our results are twofold: first, we demonstrate extensive polygenic overlap between SCZ and PLS and to a lesser extent between BD and PLS. Second, we identified novel loci shared between PLS and SCZ (n = 39), and BD (n = 8). Whereas most of the identified SCZ (66%) and BD (62%) pleiotropic risk alleles were associated with reduced lifespan, we also detected some antagonistic protective alleles associated to shorter lifespans. In fact, top-associated SNPs with SCZ seems to explain longevity variance explained (LVE) better than many other life-threatening diseases, including Type 2 diabetes and most cancers, probably due to a high overlap with smoking-related pathways. Overall, our study provides evidence of a genetic burden driven through premature mortality among people with SCZ, which can have profound implications for understanding, and potentially treating, the mortality gap associated with this psychiatric disorder.

Contribution of Evolutionary Selected Immune Gene Polymorphism to Immune-Related Disorders: The Case of Lymphocyte Scavenger Receptors CD5 and CD6.

Pathogens are one of the main selective pressures that ancestral humans had to adapt to. Components of the immune response system have been preferential targets of natural selection in response to such pathogen-driven pressure. In turn, there is compelling evidence showing that positively selected immune gene variants conferring increased resistance to past or present infectious agents are today associated with increased risk for autoimmune or inflammatory disorders but decreased risk of cancer, the other side of the same coin. CD5 and CD6 are lymphocytic scavenger receptors at the interphase of the innate and adaptive immune responses since they are involved in both: (i) microbial-associated pattern recognition; and (ii) modulation of intracellular signals mediated by the clonotypic antigen-specific receptor present in T and B cells (TCR and BCR, respectively). Here, we review available information on CD5 and CD6 as targets of natural selection as well as on the role of CD5 and CD6 variation in autoimmunity and cancer.

Targeted resequencing reveals rare variants enrichment in multiple sclerosis susceptibility genes.

Although genome-wide association studies have identified a number of common variants associated with multiple sclerosis (MS) susceptibility, little is known about the relevance of rare variants. Here, we aimed to explore the role of rare variants in 14 MS risk genes (FCRL1, RGS1, TIMMDC1, HHEX, CXCR5, LTBR, TSFM, GALC, TRAF3, STAT3, TNFSF14, IFI30, CD40, and CYP24A1) by targeted resequencing in an Iberian population of 524 MS cases and 546 healthy controls. Four rare variants-enriched regions within CYP24A1, FCRL1, RGS1, and TRAF3 were identified as significantly associated with MS. Functional studies revealed significantly decreased regulator of G protein signaling 1 (RGS1) gene expression levels in peripheral blood mononuclear cells from MS patients with RGS1 rare variants compared to noncarriers, whereas no significant differences in gene expression were observed for CYP24A1, FCRL1, and TRAF3 between rare variants carriers and noncarriers. Immunophenotyping showed significant decrease in RGS1 expression in peripheral blood B lymphocytes from MS patients with RGS1 rare variants relative to noncarriers. Lastly, peripheral blood mononuclear cell from MS patients carrying RGS1 rare variants showed significantly lower induction of RGS1 gene expression by interferon-ß compared to MS patients lacking RGS1 variants. The presence of rare variants in RGS1 reinforce the ideas of high genetic heterogeneity and a role of rare variants in MS pathogenesis.

Properties of human disease genes and the role of genes linked to Mendelian disorders in complex disease aetiology.

Do genes presenting variation that has been linked to human disease have different biological properties than genes that have never been related to disease? What is the relationship between disease and fitness? Are the evolutionary pressures that affect genes linked to Mendelian diseases the same to those acting on genes whose variation contributes to complex disorders? The answers to these questions could shed light on the architecture of human genetic disorders and may have relevant implications when designing mapping strategies in future genetic studies. Here we show that, relative to non-disease genes, human disease (HD) genes have specific evolutionary profiles and protein network properties. Additionally, our results indicate that the mutation-selection balance renders an insufficient account of the evolutionary history of some HD genes and that adaptive selection could also contribute to shape their genetic architecture. Notably, several biological features of HD genes depend on the type of pathology (complex or Mendelian) with which they are related. For example, genes harbouring both causal variants for Mendelian disorders and risk factors for complex disease traits (Complex-Mendelian genes), tend to present higher functional relevance in the protein network and higher expression levels than genes associated only with complex disorders. Moreover, risk variants in Complex-Mendelian genes tend to present higher odds ratios than those on genes associated with the same complex disorders but with no link to Mendelian diseases. Taken together, our results suggest that genetic variation at genes linked to Mendelian disorders plays an important role in driving susceptibility to complex disease.

Signatures of Evolutionary Adaptation in Quantitative Trait Loci Influencing Trace Element Homeostasis in Liver.

Essential trace elements possess vital functions at molecular, cellular, and physiological levels in health and disease, and they are tightly regulated in the human body. In order to assess variability and potential adaptive evolution of trace element homeostasis, we quantified 18 trace elements in 150 liver samples, together with the expression levels of 90 genes and abundances of 40 proteins involved in their homeostasis. Additionally, we genotyped 169 single nucleotide polymorphism (SNPs) in the same sample set. We detected significant associations for 8 protein quantitative trait loci (pQTL), 10 expression quantitative trait loci (eQTLs), and 15 micronutrient quantitative trait loci (nutriQTL). Six of these exceeded the false discovery rate cutoff and were related to essential trace elements: 1) one pQTL for GPX2 (rs10133290); 2) two previously described eQTLs for HFE (rs12346) and SELO (rs4838862) expression; and 3) three nutriQTLs: The pathogenic C282Y mutation at HFE affecting iron (rs1800562), and two SNPs within several clustered metallothionein genes determining selenium concentration (rs1811322 and rs904773). Within the complete set of significant QTLs (which involved 30 SNPs and 20 gene regions), we identified 12 SNPs with extreme patterns of population differentiation (FST values in the top 5% percentile in at least one HapMap population pair) and significant evidence for selective sweeps involving QTLs at GPX1, SELENBP1, GPX3, SLC30A9, and SLC39A8. Overall, this detailed study of various molecular phenotypes illustrates the role of regulatory variants in explaining differences in trace element homeostasis among populations and in the human adaptive response to environmental pressures related to micronutrients.

Mendelian genes for Parkinson's disease contribute to the sporadic forms of the disease.

Parkinson's disease (PD) can be divided into familial (Mendelian) and sporadic forms. A number of causal genes have been discovered for the Mendelian form, which constitutes 10-20% of the total cases. Genome-wide association studies have successfully uncovered a number of susceptibility loci for sporadic cases but those only explain a small fraction (6-7%) of PD heritability. It has been observed that some genes that confer susceptibility to PD through common risk variants also contain rare causing mutations for the Mendelian forms of the disease. These results suggest a possible functional link between Mendelian and sporadic PD and led us to investigate the role that rare and low-frequency variants could have on the sporadic form. Through a targeting approach, we have resequenced at 49× coverage the exons and regulatory regions of 38 genes (including Mendelian and susceptibility PD genes) in 249 sporadic PD patients and 145 unrelated controls of European origin. Unlike susceptibility genes, Mendelian genes show a clear general enrichment of rare functional variants in PD cases, observed directly as well as with Tajima's D statistic and several collapsing methods. Our findings suggest that rare variation on PD Mendelian genes may have a role in the sporadic forms of the disease.

Detection of genomic rearrangements from targeted resequencing data in Parkinson's disease patients.

BACKGROUND: The analysis of coverage depth in next-generation sequencing data allows the detection of gene dose alterations. We explore the frequency of such structural events in a Spanish cohort of sporadic PD cases. METHODS: Gene dose alterations were detected with the eXome-Hidden Markov Model (XHMM) software from depth of coverage in resequencing data available for 38 Mendelian and other risk PD loci in 394 individuals (249 cases and 145 controls) and subsequently validated by quantitative PCR. RESULTS: We identified 10 PD patients with exon dosage alterations in PARK2, GBA-GBAP1, and DJ1. Additional functional variants, including 2 novel nonsense mutations (p.Arg1552Ter in LRRK2 and p.Trp90Ter in PINK1), were confirmed by Sanger sequencing. This combined approach disclosed the genetic cause of 12 PD cases. CONCLUSIONS: Gene dose alterations related to PD can be correctly identified from targeting resequencing data. This approach substantially improves the detection rate of cases with causal genetic alterations. © 2016 The Authors. Movement Disorders published by Wiley Periodicals, Inc. on behalf of International Parkinson and Movement Disorder Society.

Extreme population differences in the human zinc transporter ZIP4 (SLC39A4) are explained by positive selection in Sub-Saharan Africa.

Extreme differences in allele frequency between West Africans and Eurasians were observed for a leucine-to-valine substitution (Leu372Val) in the human intestinal zinc uptake transporter, ZIP4, yet no further evidence was found for a selective sweep around the ZIP4 gene (SLC39A4). By interrogating allele frequencies in more than 100 diverse human populations and resequencing Neanderthal DNA, we confirmed the ancestral state of this locus and found a strong geographical gradient for the derived allele (Val372), with near fixation in West Africa. In extensive coalescent simulations, we show that the extreme differences in allele frequency, yet absence of a classical sweep signature, can be explained by the effect of a local recombination hotspot, together with directional selection favoring the Val372 allele in Sub-Saharan Africans. The possible functional effect of the Leu372Val substitution, together with two pathological mutations at the same codon (Leu372Pro and Leu372Arg) that cause acrodermatitis enteropathica (a disease phenotype characterized by extreme zinc deficiency), was investigated by transient overexpression of human ZIP4 protein in HeLa cells. Both acrodermatitis mutations cause absence of the ZIP4 transporter cell surface expression and nearly absent zinc uptake, while the Val372 variant displayed significantly reduced surface protein expression, reduced basal levels of intracellular zinc, and reduced zinc uptake in comparison with the Leu372 variant. We speculate that reduced zinc uptake by the ZIP4-derived Val372 isoform may act by starving certain pathogens of zinc, and hence may have been advantageous in Sub-Saharan Africa. Moreover, these functional results may indicate differences in zinc homeostasis among modern human populations with possible relevance for disease risk.

Copy number variation analysis in the great apes reveals species-specific patterns of structural variation.

Copy number variants (CNVs) are increasingly acknowledged as an important source of evolutionary novelties in the human lineage. However, our understanding of their significance is still hindered by the lack of primate CNV data. We performed intraspecific comparative genomic hybridizations to identify loci harboring copy number variants in each of the four great apes: bonobos, chimpanzees, gorillas, and orangutans. For the first time, we could analyze differences in CNV location and frequency in these four species, and compare them with human CNVs and primate segmental duplication (SD) maps. In addition, for bonobo and gorilla, patterns of CNV and nucleotide diversity were studied in the same individuals. We show that CNVs have been subject to different selective pressures in different lineages. Evidence for purifying selection is stronger in gorilla CNVs overlapping genes, while positive selection appears to have driven the fixation of structural variants in the orangutan lineage. In contrast, chimpanzees and bonobos present high levels of common structural polymorphism, which is indicative of relaxed purifying selection together with the higher mutation rates induced by the known burst of segmental duplication in the ancestor of the African apes. Indeed, the impact of the duplication burst is noticeable by the fact that bonobo and chimpanzee share more CNVs with gorilla than expected. Finally, we identified a number of interesting genomic regions that present high-frequency CNVs in all great apes, while containing only very rare or even pathogenic structural variants in humans.

Evolutionary and functional evidence for positive selection at the human CD5 immune receptor gene.

CD5 is a lymphocyte surface coreceptor of still incompletely understood function. Currently available information indicates that CD5 participates not only in cell-to-cell immune interactions through still poorly defined endogenous ligands expressed on hemopoietic and nonhemopoietic cells but also in recognition of exogenous and highly conserved microbial structures such as fungal ß-glucans. Preceding single nucleotide polymorphism (SNP) data analysis provided evidence for a recent selective sweep in East Asia and suggested a nonsynonymous substitution at position 471 (A471V; rs2229177) of the cytoplasmatic region of the CD5 receptor as the most plausible target of selection. The present report further investigates the role of natural selection in the CD5 gene by a resequencing approach in 60 individuals representing populations from 3 different continents (20 Africans, 20 Europeans and 20 East Asians) and by functionally assaying the relevance of the A471V replacement on CD5 signaling. The high differentiation pattern found at the nonsynonymous A471V site together with the low diversity, most of the performed neutrality tests (Tajima's D, Fu and Li's F* and D*, and Fu's Fs) and the predominance of a major haplotype in East Asians strongly argue in favor of positive selection for the A471V site. Importantly, anti-CD5 monoclonal antibody cross-linking unveiled significant differences among A471V variants regarding the mitogen-activated protein kinase (MAPK) cascade activation on COS7 and on human peripheral blood mononuclear cells. Similar differences on MAPK activation and IL-8 cytokine release were also observed upon exposure of HEK293 cell transfectants expressing the A471V variants to Zymosan, a ß-glucan-rich fungal particle. Taken together, the results provide evidence for the hypothesis of an adaptive role of the A471V substitution to environmental challenges, most likely infectious pathogens, in East Asian populations.

Identifying signatures of positive selection in human populations from North Africa.

Because of its location, North Africa (NA) has witnessed continuous demographic movements with an impact on the genomes of present-day human populations. Genomic data describe a complex scenario with varying proportions of at least four main ancestry components: Maghrebi, Middle Eastern-, European-, and West-and-East-African-like. However, the footprint of positive selection in NA has not been studied. Here, we compile genome-wide genotyping data from 190 North Africans and individuals from surrounding populations, investigate for signatures of positive selection using allele frequencies and linkage disequilibrium-based methods and infer ancestry proportions to discern adaptive admixture from post-admixture selection events. Our results show private candidate genes for selection in NA involved in insulin processing (KIF5A), immune function (KIF5A, IL1RN, TLR3), and haemoglobin phenotypes (BCL11A). We also detect signatures of positive selection related to skin pigmentation (SLC24A5, KITLG), and immunity function (IL1R1, CD44, JAK1) shared with European populations and candidate genes associated with haemoglobin phenotypes (HPSE2, HBE1, HBG2), other immune-related (DOCK2) traits, and insulin processing (GLIS3) traits shared with West and East African populations. Finally, the SLC8A1 gene, which codifies for a sodium-calcium exchanger, was the only candidate identified under post-admixture selection in Western NA.

Polygenic risk scores enhance prediction of body mass index increase in individuals with a first episode of psychosis.

BACKGROUND: Individuals with a first episode of psychosis (FEP) show rapid weight gain during the first months of treatment, which is associated with a reduction in general physical health. Although genetics is assumed to be a significant contributor to weight gain, its exact role is unknown. METHODS: We assembled a population-based FEP cohort of 381 individuals that was split into a Training (n = 224) set and a Validation (n = 157) set to calculate the polygenic risk score (PRS) in a two-step process. In parallel, we obtained reference genome-wide association studies for body mass index (BMI) and schizophrenia (SCZ) to examine the pleiotropic landscape between the two traits. BMI PRSs were added to linear models that included sociodemographic and clinical variables to predict BMI increase (∆BMI) in the Validation set. RESULTS: The results confirmed considerable shared genetic susceptibility for the two traits involving 449 near-independent genomic loci. The inclusion of BMI PRSs significantly improved the prediction of ∆BMI at 12 months after the onset of antipsychotic treatment by 49.4% compared to a clinical model. In addition, we demonstrated that the PRS containing pleiotropic information between BMI and SCZ predicted ∆BMI better at 3 (12.2%) and 12 months (53.2%). CONCLUSIONS: We prove for the first time that genetic factors play a key role in determining ∆BMI during the FEP. This finding has important clinical implications for the early identification of individuals most vulnerable to weight gain and highlights the importance of examining genetic pleiotropy in the context of medically important comorbidities for predicting future outcomes.

Sleep Disturbances in At-Risk Mental States and First Episode of Psychosis: A Narrative Review on Interventions.

Sleep disturbances are a common yet often overlooked symptom of psychosis that can drastically affect the quality of life and well-being of those living with the condition. Sleep disorders are common in people diagnosed with schizophrenia and have significant negative effects on the clinical course of the illness and the functional outcomes and quality of life of patients. There is a limited number of studies addressing this question in first-episode psychosis (FEP). In this narrative review, we aimed to provide an overview of sleep disorders in populations with FEP and at-risk mental states (ARMS). The review was focused on the various treatments currently used for sleep disorders, including both non-pharmacological and pharmacological treatments. A total of 48 studies were included. We found that sleep disturbances are associated with attenuated psychotic symptoms and other psychopathological symptoms in ARMSs. The association of sleep disturbances with the transition to psychosis has been poorly investigated. Sleep disturbances have an impact on the quality of life and the psychopathological symptoms of people suffering from FEP. The non-pharmacological treatments include cognitive behavioral therapy for insomnia, bright light therapy, cognitive restructuring techniques, sleep restriction therapy, basic sleep hygiene education, and the provision of portable sleep trackers. Other treatments include antipsychotics in acute phases and melatonin. The early intervention in sleep disturbances may improve overall prognosis in emerging psychosis populations.

Understanding signatures of positive natural selection in human zinc transporter genes.

Zinc is an essential micronutrient with a tightly regulated systemic and cellular homeostasis. In humans, some zinc transporter genes (ZTGs) have been previously reported as candidates for strong geographically restricted selective sweeps. However, since zinc homeostasis is maintained by the joint action of 24 ZTGs, other more subtle modes of selection could have also facilitated human adaptation to zinc availability. Here, we studied whether the complete set of ZTGs are enriched for signals of positive selection in worldwide populations and population groups from South Asia. ZTGs showed higher levels of genetic differentiation between African and non-African populations than would be randomly expected, as well as other signals of polygenic selection outside Africa. Moreover, in several South Asian population groups, ZTGs were significantly enriched for SNPs with unusually extended haplotypes and displayed SNP genotype-environmental correlations when considering zinc deficiency levels in soil in that geographical area. Our study replicated some well-characterized targets for positive selection in East Asia and sub-Saharan Africa, and proposes new candidates for follow-up in South Asia (SLC39A5) and Africa (SLC39A7). Finally, we identified candidate variants for adaptation in ZTGs that could contribute to different disease susceptibilities and zinc-related human health traits.

Recent human evolution has shaped geographical differences in susceptibility to disease.

BACKGROUND: Searching for associations between genetic variants and complex diseases has been a very active area of research for over two decades. More than 51,000 potential associations have been studied and published, a figure that keeps increasing, especially with the recent explosion of array-based Genome-Wide Association Studies. Even if the number of true associations described so far is high, many of the putative risk variants detected so far have failed to be consistently replicated and are widely considered false positives. Here, we focus on the world-wide patterns of replicability of published association studies. RESULTS: We report three main findings. First, contrary to previous results, genes associated to complex diseases present lower degrees of genetic differentiation among human populations than average genome-wide levels. Second, also contrary to previous results, the differences in replicability of disease associated-loci between Europeans and East Asians are highly correlated with genetic differentiation between these populations. Finally, highly replicated genes present increased levels of high-frequency derived alleles in European and Asian populations when compared to African populations. CONCLUSIONS: Our findings highlight the heterogeneous nature of the genetic etiology of complex disease, confirm the importance of the recent evolutionary history of our species in current patterns of disease susceptibility and could cast doubts on the status as false positives of some associations that have failed to replicate across populations.