Pleiotropic contribution of rbfox1 to psychiatric and neurodevelopmental phenotypes in two zebrafish models.

RBFOX1 is a highly pleiotropic gene that contributes to several psychiatric and neurodevelopmental disorders. Both rare and common variants in RBFOX1 have been associated with several psychiatric conditions, but the mechanisms underlying the pleiotropic effects of RBFOX1 are not yet understood. Here we found that, in zebrafish, rbfox1 is expressed in spinal cord, mid- and hindbrain during developmental stages. In adults, expression is restricted to specific areas of the brain, including telencephalic and diencephalic regions with an important role in receiving and processing sensory information and in directing behaviour. To investigate the contribution of rbfox1 to behaviour, we used rbfox1sa15940, a zebrafish mutant line with TL background. We found that rbfox1sa15940 mutants present hyperactivity, thigmotaxis, decreased freezing behaviour and altered social behaviour. We repeated these behavioural tests in a second rbfox1 mutant line with a different genetic background (TU), rbfox1del19, and found that rbfox1 deficiency affects behaviour similarly in this line, although there were some differences. rbfox1del19 mutants present similar thigmotaxis, but stronger alterations in social behaviour and lower levels of hyperactivity than rbfox1sa15940 fish. Taken together, these results suggest that mutations in rbfox1 lead to multiple behavioural changes in zebrafish that might be modulated by environmental, epigenetic and genetic background effects, and that resemble phenotypic alterations present in Rbfox1-deficient mice and in patients with different psychiatric conditions. Our study, thus, highlights the evolutionary conservation of rbfox1 function in behaviour and paves the way to further investigate the mechanisms underlying rbfox1 pleiotropy on the onset of neurodevelopmental and psychiatric disorders.

BridgePRS leverages shared genetic effects across ancestries to increase polygenic risk score portability.

Here we present BridgePRS, a novel Bayesian polygenic risk score (PRS) method that leverages shared genetic effects across ancestries to increase PRS portability. We evaluate BridgePRS via simulations and real UK Biobank data across 19 traits in individuals of African, South Asian and East Asian ancestry, using both UK Biobank and Biobank Japan genome-wide association study summary statistics; out-of-cohort validation is performed in the Mount Sinai (New York) BioMe biobank. BridgePRS is compared with the leading alternative, PRS-CSx, and two other PRS methods. Simulations suggest that the performance of BridgePRS relative to PRS-CSx increases as uncertainty increases: with lower trait heritability, higher polygenicity and greater between-population genetic diversity; and when causal variants are not present in the data. In real data, BridgePRS has a 61% larger average R2 than PRS-CSx in out-of-cohort prediction of African ancestry samples in BioMe (P = 6 × 10-5). BridgePRS is a computationally efficient, user-friendly and powerful approach for PRS analyses in non-European ancestries.

Trumpet plots: visualizing the relationship between allele frequency and effect size in genetic association studies.

Recent advances in genome-wide association and sequencing studies have shown that the genetic architecture of complex traits and diseases involves a combination of rare and common genetic variants distributed throughout the genome. One way to better understand this architecture is to visualize genetic associations across a wide range of allele frequencies. However, there is currently no standardized or consistent graphical representation for effectively illustrating these results. Here we propose a standardized approach for visualizing the effect size of risk variants across the allele frequency spectrum. The proposed plots have a distinctive trumpet shape: with the majority of variants having high frequency and small effects, and a small number of variants having lower frequency and larger effects. To demonstrate the utility of trumpet plots in illustrating the relationship between the number of variants, their frequency, and the magnitude of their effects in shaping the genetic architecture of complex traits and diseases, we generated trumpet plots for more than one hundred traits in the UK Biobank. To facilitate their broader use, we developed an R package, 'TrumpetPlots' (available at the Comprehensive R Archive Network) and R Shiny application, 'Shiny Trumpets' (available at https://juditgg.shinyapps.io/shinytrumpets/) that allows users to explore these results and submit their own data.

BridgePRS : A powerful trans-ancestry Polygenic Risk Score method.

Polygenic Risk Scores (PRS) have huge potential to contribute to biomedical research and to a future of precision medicine, but to date their calculation relies largely on Europeanancestry GWAS data. This global bias makes most PRS substantially less accurate in individuals of non-European ancestry. Here we present BridgePRS , a novel Bayesian PRS method that leverages shared genetic effects across ancestries to increase the accuracy of PRS in non-European populations. The performance of BridgePRS is evaluated in simulated data and real UK Biobank (UKB) data across 19 traits in African, South Asian and East Asian ancestry individuals, using both UKB and Biobank Japan GWAS summary statistics. BridgePRS is compared to the leading alternative, PRS-CSx , and two single-ancestry PRS methods adapted for trans-ancestry prediction. PRS trained in the UK Biobank are then validated out-of-cohort in the independent Mount Sinai (New York) Bio Me Biobank. Simulations reveal that BridgePRS performance, relative to PRS-CSx , increases as uncertainty increases: with lower heritability, higher polygenicity, greater between-population genetic diversity, and when causal variants are not present in the data. Our simulation results are consistent with real data analyses in which BridgePRS has better predictive accuracy in African ancestry samples, especially in out-of-cohort prediction (into Bio Me ), which shows a 60% boost in mean R 2 compared to PRS-CSx ( P = 2 × 10 -6 ). BridgePRS performs the full PRS analysis pipeline, is computationally efficient, and is a powerful method for deriving PRS in diverse and under-represented ancestry populations.

Pleiotropic contribution of rbfox1 to psychiatric and neurodevelopmental phenotypes in a zebrafish model.

RBFOX1 is a highly pleiotropic gene that contributes to several psychiatric and neurodevelopmental disorders. Both rare and common variants in RBFOX1 have been associated with several psychiatric conditions, but the mechanisms underlying the pleiotropic effects of RBFOX1 are not yet understood. Here we found that, in zebrafish, rbfox1 is expressed in spinal cord, mid- and hindbrain during developmental stages. In adults, expression is restricted to specific areas of the brain, including telencephalic and diencephalic regions with an important role in receiving and processing sensory information and in directing behaviour. To investigate the effect of rbfox1 deficiency on behaviour, we used rbfox1sa15940, a rbfox1 loss-of-function line. We found that rbfox1sa15940 mutants present hyperactivity, thigmotaxis, decreased freezing behaviour and altered social behaviour. We repeated these behavioural tests in a second rbfox1 loss-of-function line with a different genetic background, rbfox1del19, and found that rbfox1 deficiency affects behaviour similarly in this line, although there were some differences. rbfox1del19 mutants present similar thigmotaxis, but stronger alterations in social behaviour and lower levels of hyperactivity than rbfox1sa15940 fish. Taken together, these results suggest that rbfox1 deficiency leads to multiple behavioural changes in zebrafish that might be modulated by environmental, epigenetic and genetic background effects, and that resemble phenotypic alterations present in Rbfox1-deficient mice and in patients with different psychiatric conditions. Our study thus highlights the evolutionary conservation of rbfox1 function in behaviour and paves the way to further investigate the mechanisms underlying rbfox1 pleiotropy on the onset of neurodevelopmental and psychiatric disorders.

PRSet: Pathway-based polygenic risk score analyses and software.

Polygenic risk scores (PRSs) have been among the leading advances in biomedicine in recent years. As a proxy of genetic liability, PRSs are utilised across multiple fields and applications. While numerous statistical and machine learning methods have been developed to optimise their predictive accuracy, these typically distil genetic liability to a single number based on aggregation of an individual's genome-wide risk alleles. This results in a key loss of information about an individual's genetic profile, which could be critical given the functional sub-structure of the genome and the heterogeneity of complex disease. In this manuscript, we introduce a 'pathway polygenic' paradigm of disease risk, in which multiple genetic liabilities underlie complex diseases, rather than a single genome-wide liability. We describe a method and accompanying software, PRSet, for computing and analysing pathway-based PRSs, in which polygenic scores are calculated across genomic pathways for each individual. We evaluate the potential of pathway PRSs in two distinct ways, creating two major sections: (1) In the first section, we benchmark PRSet as a pathway enrichment tool, evaluating its capacity to capture GWAS signal in pathways. We find that for target sample sizes of >10,000 individuals, pathway PRSs have similar power for evaluating pathway enrichment as leading methods MAGMA and LD score regression, with the distinct advantage of providing individual-level estimates of genetic liability for each pathway -opening up a range of pathway-based PRS applications, (2) In the second section, we evaluate the performance of pathway PRSs for disease stratification. We show that using a supervised disease stratification approach, pathway PRSs (computed by PRSet) outperform two standard genome-wide PRSs (computed by C+T and lassosum) for classifying disease subtypes in 20 of 21 scenarios tested. As the definition and functional annotation of pathways becomes increasingly refined, we expect pathway PRSs to offer key insights into the heterogeneity of complex disease and treatment response, to generate biologically tractable therapeutic targets from polygenic signal, and, ultimately, to provide a powerful path to precision medicine.

baz1b loss-of-function in zebrafish produces phenotypic alterations consistent with the domestication syndrome.

BAZ1B is a ubiquitously expressed nuclear protein with roles in chromatin remodeling, DNA replication and repair, and transcription. Reduced BAZ1B expression disrupts neuronal and neural crest development. Variation in the activity of BAZ1B has been proposed to underly morphological and behavioral aspects of domestication through disruption of neural crest development. Knockdown of baz1b in Xenopus embryos and Baz1b loss-of-function (LoF) in mice leads to craniofacial defects consistent with this hypothesis. We generated baz1b LoF zebrafish using CRISPR/Cas9 gene editing to test the hypothesis that baz1b regulates behavioral phenotypes associated with domestication in addition to craniofacial features. Zebrafish with baz1b LoF show mild underdevelopment at larval stages and distinctive craniofacial features later in life. Mutant zebrafish show reduced anxiety-associated phenotypes and an altered ontogeny of social behaviors. Thus, in zebrafish, developmental deficits in baz1b recapitulate both morphological and behavioral phenotypes associated with the domestication syndrome in other species.

Ankk1 Loss of Function Disrupts Dopaminergic Pathways in Zebrafish.

Ankyrin repeat and kinase domain containing 1 (ANKK1) is a member of the receptor-interacting protein serine/threonine kinase family, known to be involved in cell proliferation, differentiation and activation of transcription factors. Genetic variation within the ANKK1 locus is suggested to play a role in vulnerability to addictions. However, ANKK1 mechanism of action is still poorly understood. It has been suggested that ANKK1 may affect the development and/or functioning of dopaminergic pathways. To test this hypothesis, we generated a CRISPR-Cas9 loss of function ankk1 zebrafish line causing a 27 bp insertion that disrupts the ankk1 sequence introducing an early stop codon. We found that ankk1 transcript levels were significantly lower in ankk1 mutant (ankk127ins ) fish compared to their wild type (ankk1 +/+) siblings. In ankk1 +/+ adult zebrafish brain, ankk1 protein was detected in isocortex, hippocampus, basolateral amygdala, mesencephalon, and cerebellum, resembling the mammalian distribution pattern. In contrast, ankk1 protein was reduced in the brain of ankk127ins/27ins fish. Quantitative polymerase chain reaction analysis revealed an increase in expression of drd2b mRNA in ankk127ins at both larval and adult stages. In ankk1 +/+ adult zebrafish brain, drd2 protein was detected in cerebral cortex, cerebellum, hippocampus, and caudate homolog regions, resembling the pattern in humans. In contrast, drd2 expression was reduced in cortical regions of ankk127ins/27ins being predominantly found in the hindbrain. No differences in the number of cell bodies or axonal projections detected by anti-tyrosine hydroxylase immunostaining on 3 days post fertilization (dpf) larvae were found. Behavioral analysis revealed altered sensitivity to effects of both amisulpride and apomorphine on locomotion and startle habituation, consistent with a broad loss of both pre and post synaptic receptors. Ankk127ins mutants showed reduced sensitivity to the effect of the selective dopamine receptor antagonist amisulpride on locomotor responses to acoustic startle and were differentially sensitive to the effects of the non-selective dopamine agonist apomorphine on both locomotion and habituation. Taken together, our findings strengthen the hypothesis of a functional relationship between ANKK1 and DRD2, supporting a role for ANKK1 in the maintenance and/or functioning of dopaminergic pathways. Further work is needed to disentangle ANKK1's role at different developmental stages.

Behavioral Effects of Developmental Exposure to JWH-018 in Wild-Type and Disrupted in Schizophrenia 1 (disc1) Mutant Zebrafish.

Synthetic cannabinoids can cause acute adverse psychological effects, but the potential impact when exposure happens before birth is unknown. Use of synthetic cannabinoids during pregnancy may affect fetal brain development, and such effects could be moderated by the genetic makeup of an individual. Disrupted in schizophrenia 1 (DISC1) is a gene with important roles in neurodevelopment that has been associated with psychiatric disorders in pedigree analyses. Using zebrafish as a model, we investigated (1) the behavioral impact of developmental exposure to 3 µM 1-pentyl-3-(1-naphthoyl)-indole (JWH-018; a common psychoactive synthetic cannabinoid) and (2) whether disc1 moderates the effects of JWH-018. As altered anxiety responses are seen in several psychiatric disorders, we focused on zebrafish anxiety-like behavior. Zebrafish embryos were exposed to JWH-018 from one to six days post-fertilization. Anxiety-like behavior was assessed using forced light/dark and acoustic startle assays in larvae and novel tank diving in adults. Compared to controls, both acutely and developmentally exposed zebrafish larvae had impaired locomotion during the forced light/dark test, but anxiety levels and response to startle stimuli were unaltered. Adult zebrafish developmentally exposed to JWH-018 spent less time on the bottom of the tank, suggesting decreased anxiety. Loss-of-function in disc1 increased anxiety-like behavior in the tank diving assay but did not alter sensitivity to JWH-018. Results suggest developmental exposure to JWH-018 has a long-term behavioral impact in zebrafish, which is not moderated by disc1.

The effects of polygenic risk for psychiatric disorders and smoking behaviour on psychotic experiences in UK Biobank.

While psychotic experiences are core symptoms of mental health disorders like schizophrenia, they are also reported by 5-10% of the population. Both smoking behaviour and genetic risk for psychiatric disorders have been associated with psychotic experiences, but the interplay between these factors remains poorly understood. We tested whether smoking status, maternal smoking around birth, and number of packs smoked/year were associated with lifetime occurrence of three psychotic experiences phenotypes: delusions (n = 2067), hallucinations (n = 6689), and any psychotic experience (delusions or hallucinations; n = 7803) in 157,366 UK Biobank participants. We next calculated polygenic risk scores for schizophrenia (PRSSCZ), bipolar disorder (PRSBP), major depression (PRSDEP) and attention deficit hyperactivity disorder (PRSADHD) in 144,818 UK Biobank participants of European ancestry to assess whether association between smoking and psychotic experiences was attenuated after adjustment of diagnosis of psychiatric disorders and the PRSs. Finally, we investigated whether smoking exacerbates the effects of genetic predisposition on the psychotic phenotypes in gene-environment interaction models. Smoking status, maternal smoking, and number of packs smoked/year were associated with psychotic experiences (p < 1.77 × 10-5). Except for packs smoked/year, effects were attenuated but remained significant after adjustment for diagnosis of psychiatric disorders and PRSs (p < 1.99 × 10-3). Gene-environment interaction models showed the effects of PRSDEP and PRSADHD (but not PRSSCZ or PRSBP) on delusions (but not hallucinations) were significantly greater in current smokers compared to never smokers (p < 0.002). There were no significant gene-environment interactions for maternal smoking nor for number of packs smoked/year. Our results suggest that both genetic risk of psychiatric disorders and smoking status may have independent and synergistic effects on specific types of psychotic experiences.

Identification of slit3 as a locus affecting nicotine preference in zebrafish and human smoking behaviour.

To facilitate smoking genetics research we determined whether a screen of mutagenized zebrafish for nicotine preference could predict loci affecting smoking behaviour. From 30 screened F3 sibling groups, where each was derived from an individual ethyl-nitrosurea mutagenized F0 fish, two showed increased or decreased nicotine preference. Out of 25 inactivating mutations carried by the F3 fish, one in the slit3 gene segregated with increased nicotine preference in heterozygous individuals. Focussed SNP analysis of the human SLIT3 locus in cohorts from UK (n=863) and Finland (n=1715) identified two variants associated with cigarette consumption and likelihood of cessation. Characterisation of slit3 mutant larvae and adult fish revealed decreased sensitivity to the dopaminergic and serotonergic antagonist amisulpride, known to affect startle reflex that is correlated with addiction in humans, and increased htr1aa mRNA expression in mutant larvae. No effect on neuronal pathfinding was detected. These findings reveal a role for SLIT3 in development of pathways affecting responses to nicotine in zebrafish and smoking in humans.

Men and women differ in their perception of gender bias in research institutions.

There is extensive evidence of gender inequality in research leading to insufficient representation of women in leadership positions. Numbers revealing a gender gap in research are periodically reported by national and international institutions but data on perceptions of gender equality within the research community are scarce. In the present study, a questionnaire based on the British Athena Survey of Science, Engineering and Technology (ASSET 2016) was distributed among researchers working in Spain. Consistent with the original UK-based study, women in research perceived a greater degree of gender inequality than men. This difference was consistent from junior to senior positions, within public and private universities as well as research centres, and across all research disciplines. When responses were compared with the existing UK-based questionnaire, researchers in Spain felt that women and men are treated more equally in the workplace, yet they perceived their home departments to be less supportive regarding matters of gender equality. The results of this study provide clear evidence that men and women do not share the same perceptions of gender equality in science and that their differing perceptions are relatively consistent across two major European countries. The fact that men occupy the majority of senior positions while not perceiving the same inequality as women do, may be critical when it comes to ensuring the fair ascent of women to senior positions in an academic system. These data encourage the implementation of measures to ensure that both men and women are aware of gender biases in research.

Body composition in anorexia nervosa: Meta-analysis and meta-regression of cross-sectional and longitudinal studies.

OBJECTIVE: Clinically, anorexia nervosa (AN) presents with altered body composition. We quantified these alterations and evaluated their relationships with metabolites and hormones in patients with AN longitudinally. METHOD: In accordance with PRISMA guidelines, we conducted 94 meta-analyses on 62 samples published during 1996-2019, comparing up to 2,319 pretreatment, posttreatment, and weight-recovered female patients with AN with up to 1,879 controls. Primary outcomes were fat mass, fat-free mass, body fat percentage, and their regional distribution. Secondary outcomes were bone mineral density, metabolites, and hormones. Meta-regressions examined relationships among those measures and moderators. RESULTS: Pretreatment female patients with AN evidenced 50% lower fat mass (mean difference [MD]: -8.80 kg, 95% CI: -9.81, -7.79, Q = 1.01 × 10-63 ) and 4.98 kg (95% CI: -5.85, -4.12, Q = 1.99 × 10-28 ) lower fat-free mass, with fat mass preferentially stored in the trunk region during early weight restoration (4.2%, 95% CI: -2.1, -6.2, Q = 2.30 × 10-4 ). While the majority of traits returned to levels seen in healthy controls after weight restoration, fat-free mass (MD: -1.27 kg, 95% CI: -1.79, -0.75, Q = 5.49 × 10-6 ) and bone mineral density (MD: -0.10 kg, 95% CI: -0.18, -0.03, Q = 0.01) remained significantly altered. DISCUSSION: Body composition is markedly altered in AN, warranting research into these phenotypes as clinical risk or relapse predictors. Notably, the long-term altered levels of fat-free mass and bone mineral density suggest that these parameters should be investigated as potential AN trait markers. RESUMENOBJETIVO: Clínicamente, la anorexia nervosa (AN) se presenta con alteraciones en la composición corporal. Cuantificamos estas alteraciones y evaluamos longitudinalmente su relación con metabolitos y hormonas en pacientes con AN. MÉTODO: De acuerdo con las pautas PRISMA, realizamos 94 meta-análisis en 62 muestras publicadas entre 1996-2019, comparando hasta 2,319 pacientes mujeres en pre-tratamiento, post-tratamiento, y recuperadas en base al peso con hasta 1,879 controles. Las principales medidas fueron masa grasa, masa libre de grasa, porcentaje de grasa corporal y su distribución regional. Las medidas secundarias fueron densidad mineral ósea, metabolitos y hormonas. Las meta-regresiones examinaron las relaciones entre esas medidas y moderadores. RESULTADOS: Las pacientes femeninas con AN pre-tratamiento mostraron un 50% menos de masa grasa (MD: -8.80 kg, CI 95%: -9.81, -7.79, Q = 1.01 × 10-63 ) y 4.98 kg (CI 95%: -5.85, -4.12, Q = 1.99 × 10-28 ) menos de masa libre de grasa, con masa grasa preferentemente almacenada en la región del tronco durante la recuperación temprana del peso (4.2%, CI 95%: -2.1, -6.2, Q = 2.30 × 10-4 ). Aunque la mayoría de los rasgos regresaron a los niveles vistos en los controles sanos después de la restauración del peso, la masa libre de grasa (MD: -1.27 kg, CI 95%: -1.79, -0.75, Q = 5.49 × 10-6 ) y la densidad mineral ósea (MD: -0.10 kg, CI 95%: -0.18, -0.03, Q = 0.01) permanecieron significativamente alteradas. DISCUSIÓN: La composición corporal es marcadamente alterada en la AN, lo que garantiza la investigación en estos fenotipos como predictores de riesgo clínico o de recaída. Notablemente, la alteración a largo plazo de los niveles de masa libre de grasa y densidad mineral ósea sugieren que estos parámetros debe ser investigados como potenciales rasgos indicadores de AN.

Immunomodulatory effects of antipsychotic treatment on gene expression in first-episode psychosis.

Previous studies suggest immunological alterations in patients with first-episode psychosis (FEP). Some studies show that antipsychotic compounds may cause immunomodulatory effects. To evaluate the immunological changes and the possible immunomodulatory effects in FEP, we recruited patients with FEP (n = 67) and matched controls (n = 38), aged 18-40 years, from the catchment area of the Helsinki University Hospital and the City of Helsinki, Finland. Fasting peripheral blood samples were collected between 8 and 10 a.m. in 10 ml PAXgene tubes. We applied the NanoString nCounter in-solution hybridization technology to determine gene expression levels of 147 candidate genes reflecting activation of the immune system. Cases had higher gene expression levels of BDKRB1 and SPP1/osteopontin compared with controls. Of the individual medications used as monotherapy, risperidone was associated with a statistically significant upregulation of 11 immune system genes, including cytokines and cytokine receptors (SPP1, IL1R1, IL1R2), pattern recognition molecules (TLR1, TLR2 and TLR6, dectin-1/CLEC7A), molecules involved in apoptosis (FAS), and some other molecules with functions in immune activation (BDKRB1, IGF1R, CR1). In conclusion, risperidone possessed strong immunomodulatory properties affecting mainly innate immune response in FEP patients, whereas the observed effects of quetiapine and olanzapine were only marginal. Our results further emphasize the importance of understanding the immunomodulatory mechanisms of antipsychotic treatment, especially in terms of specific compounds, doses and duration of medication in patients with severe mental illness. Future studies should evaluate the response pre- and post-treatment, and the possible role of this inflammatory activation for the progression of psychiatric and metabolic symptoms.

Pharmacogenetics of antidepressant response: A polygenic approach.

BACKGROUND: Major depressive disorder (MDD) has a high personal and socio-economic burden and >60% of patients fail to achieve remission with the first antidepressant. The biological mechanisms behind antidepressant response are only partially known but genetic factors play a relevant role. A combined predictor across genetic variants may be useful to investigate this complex trait. METHODS: Polygenic risk scores (PRS) were used to estimate multi-allelic contribution to: 1) antidepressant efficacy; 2) its overlap with MDD and schizophrenia. We constructed PRS and tested whether these predicted symptom improvement or remission from the GENDEP study (n=736) to the STAR*D study (n=1409) and vice-versa, including the whole sample or only patients treated with escitalopram or citalopram. Using summary statistics from Psychiatric Genomics Consortium for MDD and schizophrenia, we tested whether PRS from these disorders predicted symptom improvement in GENDEP, STAR*D, and five further studies (n=3756). RESULTS: No significant prediction of antidepressant efficacy was obtained from PRS in GENDEP/STAR*D but this analysis might have been underpowered. There was no evidence of overlap in the genetics of antidepressant response with either MDD or schizophrenia, either in individual studies or a meta-analysis. Stratifying by antidepressant did not alter the results. DISCUSSION: We identified no significant predictive effect using PRS between pharmacogenetic studies. The genetic liability to MDD or schizophrenia did not predict response to antidepressants, suggesting differences between the genetic component of depression and treatment response. Larger or more homogeneous studies will be necessary to obtain a polygenic predictor of antidepressant response.

Validation and genotyping of multiple human polymorphic inversions mediated by inverted repeats reveals a high degree of recurrence.

In recent years different types of structural variants (SVs) have been discovered in the human genome and their functional impact has become increasingly clear. Inversions, however, are poorly characterized and more difficult to study, especially those mediated by inverted repeats or segmental duplications. Here, we describe the results of a simple and fast inverse PCR (iPCR) protocol for high-throughput genotyping of a wide variety of inversions using a small amount of DNA. In particular, we analyzed 22 inversions predicted in humans ranging from 5.1 kb to 226 kb and mediated by inverted repeat sequences of 1.6-24 kb. First, we validated 17 of the 22 inversions in a panel of nine HapMap individuals from different populations, and we genotyped them in 68 additional individuals of European origin, with correct genetic transmission in ∼ 12 mother-father-child trios. Global inversion minor allele frequency varied between 1% and 49% and inversion genotypes were consistent with Hardy-Weinberg equilibrium. By analyzing the nucleotide variation and the haplotypes in these regions, we found that only four inversions have linked tag-SNPs and that in many cases there are multiple shared SNPs between standard and inverted chromosomes, suggesting an unexpected high degree of inversion recurrence during human evolution. iPCR was also used to check 16 of these inversions in four chimpanzees and two gorillas, and 10 showed both orientations either within or between species, providing additional support for their multiple origin. Finally, we have identified several inversions that include genes in the inverted or breakpoint regions, and at least one disrupts a potential coding gene. Thus, these results represent a significant advance in our understanding of inversion polymorphism in human populations and challenge the common view of a single origin of inversions, with important implications for inversion analysis in SNP-based studies.