Gene-based pleiotropy across migraine with aura and migraine without aura patient groups

Introduction: It is unclear whether patients diagnosed according to International Classification of Headache Disorders criteria for migraine with aura (MA) and migraine without aura (MO) experience distinct disorders or whether their migraine subtypes are genetically related. Aim: Using a novel gene-based (statistical) approach, we aimed to identify individual genes and pathways associated both with MA and MO. Methods: Gene-based tests were performed using genome-wide association summary statistic results from the most recent International Headache Genetics Consortium study comparing 4505 MA cases with 34,813 controls and 4038 MO cases with 40,294 controls. After accounting for non-independence of gene-based test results, we examined the significance of the proportion of shared genes associated with MA and MO. Results: We found a significant overlap in genes associated with MA and MO. Of the total 1514 genes with a nominally significant gene-based p value (pgene-based ≤ 0.05) in the MA subgroup, 107 also produced pgene-based ≤ 0.05 in the MO subgroup. The proportion of overlapping genes is almost double the empirically derived null expectation, producing significant evidence of gene-based overlap (pleiotropy) (pbinomial-test = 1.5 × 10–4). Combining results across MA and MO, six genes produced genome-wide significant gene-based p values. Four of these genes (TRPM8, UFL1, FHL5 and LRP1) were located in close proximity to previously reported genome-wide significant SNPs for migraine, while two genes, TARBP2 and NPFF separated by just 259 bp on chromosome 12q13.13, represent a novel risk locus. The genes overlapping in both migraine types were enriched for functions related to inflammation, the cardiovascular system and connective tissue. Conclusions: Our results provide novel insight into the likely genes and biological mechanisms that underlie both MA and MO, and when combined with previous data, highlight the neuropeptide FF-amide peptide encoding gene (NPFF) as a novel candidate risk gene for both types of migraine.

CPAG: software for leveraging pleiotropy in GWAS to reveal similarity between human traits links plasma fatty acids and intestinal inflammation.

Meta-analyses of genome-wide association studies (GWAS) have demonstrated that the same genetic variants can be associated with multiple diseases and other complex traits. We present software called CPAG (Cross-Phenotype Analysis of GWAS) to look for similarities between 700 traits, build trees with informative clusters, and highlight underlying pathways. Clusters are consistent with pre-defined groups and literature-based validation but also reveal novel connections. We report similarity between plasma palmitoleic acid and Crohn's disease and find that specific fatty acids exacerbate enterocolitis in zebrafish. CPAG will become increasingly powerful as more genetic variants are uncovered, leading to a deeper understanding of complex traits. CPAG is freely available at www.sourceforge.net/projects/CPAG/.

Genetic pleiotropy between multiple sclerosis and schizophrenia but not bipolar disorder: differential involvement of immune-related gene loci

Converging evidence implicates immune abnormalities in schizophrenia (SCZ), and recent genome-wide association studies (GWAS) have identified immune-related single-nucleotide polymorphisms (SNPs) associated with SCZ. Using the conditional false discovery rate (FDR) approach, we evaluated pleiotropy in SNPs associated with SCZ (n=21 856) and multiple sclerosis (MS) (n=43 879), an inflammatory, demyelinating disease of the central nervous system. Because SCZ and bipolar disorder (BD) show substantial clinical and genetic overlap, we also investigated pleiotropy between BD (n=16 731) and MS. We found significant genetic overlap between SCZ and MS and identified 21 independent loci associated with SCZ, conditioned on association with MS. This enrichment was driven by the major histocompatibility complex (MHC). Importantly, we detected the involvement of the same human leukocyte antigen (HLA) alleles in both SCZ and MS, but with an opposite directionality of effect of associated HLA alleles (that is, MS risk alleles were associated with decreased SCZ risk). In contrast, we found no genetic overlap between BD and MS. Considered together, our findings demonstrate genetic pleiotropy between SCZ and MS and suggest that the MHC signals may differentiate SCZ from BD susceptibility.Molecular Psychiatry advance online publication, 28 January 2014; doi:10.1038/mp.2013.195.

Pleiotropy in microdeletion syndromes: neurologic and spermatogenic abnormalities in mice homozygous for the p6H deletion are likely due to dysfunction of a single gene.

Variability and complexity of phenotypes observed in microdeletion syndromes can be due to deletion of a single gene whose product participates in several aspects of development or can be due to the deletion of a number of tightly linked genes, each adding its own effect to the syndrome. The p6H deletion in mouse chromosome 7 presents a good model with which to address this question of multigene vs. single-gene pleiotropy. Mice homozygous for the p6H deletion are diluted in pigmentation, are smaller than their littermates, and manifest a nervous jerky-gait phenotype. Male homozygotes are sterile and exhibit profound abnormalities in spermiogenesis. By using N-ethyl-N-nitrosourea (EtNU) mutagenesis and a breeding protocol designed to recover recessive mutations expressed hemizygously opposite a large p-locus deletion, we have generated three noncomplementing mutations that map to the p6H deletion. Each of these EtNU-induced mutations has adverse effects on the size, nervous behavior, and progression of spermiogenesis that characterize p6H deletion homozygotes. Because EtNU is thought to induce primarily intragenic (point) mutations in mouse stem-cell spermatogonia, we propose that the trio of phenotypes (runtiness, nervous jerky gait, and male sterility) expressed in p6H deletion homozygotes is the result of deletion of a single highly pleiotropic gene. We also predict that a homologous single locus, quite possibly tightly linked and distal to the D15S12 (P) locus in human chromosome 15q11-q13, may be associated with similar developmental abnormalities in humans.

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn's disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn's disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.66 x 10-10, odds ratio (OR) = 0.74, 95% CI:0.68-0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10-4. OR = 0.86 (95% CI:0.79-0.93); rs744166, P = 2.6×10-5, OR = 0.84 (95% CI:0.77-0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10-5, OR = 0.65 (95% CI:0.54-0.79)), and further associations were detected for IL12B (rs10045431, P = 5.261025, OR = 0.83 (95% CI:0.76-0.91)), CDKAL1 (rs6908425, P = 1.1×10-4, OR = 0.82 (95% CI:0.74-0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10-5, OR = 1.92 (95% CI: 1.38-2.67)), and chr13q14 (rs3764147, P = 5.9×10-4, OR = 1.19 (95% CI: 1.08-1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn's disease and further highlight the involvement of common risk variants across multiple diseases.

Genetics of ankylosing spondylitis

Ankylosing spondylitis (AS) is a chronic inflammatory arthritis that affects the spine and sacroiliac joints. It causes significant disability and is associated with a number of other features including peripheral arthritis, anterior uveitis, psoriasis and inflammatory bowel disease (IBD). Significant progress has been made in the genetics of AS have in the last five years, leading to new treatments in trial, and major leaps in understanding of the aetiopathogenesis of the disease.

Shared common variants in prostate cancer and blood lipids

Background Epidemiological and clinical studies suggest comorbidity between prostate cancer (PCA) and cardiovascular disease (CVD) risk factors. However, the relationship between these two phenotypes is still not well understood. Here we sought to identify shared genetic loci between PCA and CVD risk factors. Methods We applied a genetic epidemiology method based on conjunction false discovery rate (FDR) that combines summary statistics from different genome-wide association studies (GWAS), and allows identification of genetic overlap between two phenotypes. We evaluated summary statistics from large, multi-centre GWA studies of PCA (n = 50 000) and CVD risk factors (n = 200 000) [triglycerides (TG), low-density lipoprotein (LDL) cholesterol and high-density lipoprotein (HDL) cholesterol, systolic blood pressure, body mass index, waist-hip ratio and type 2 diabetes (T2D)]. Enrichment of single nucleotide polymorphisms (SNPs) associated with PCA and CVD risk factors was assessed with conditional quantile-quantile plots and the Anderson-Darling test. Moreover, we pinpointed shared loci using conjunction FDR. Results We found the strongest enrichment of P-values in PCA was conditional on LDL and conditional on TG. In contrast, we found only weak enrichment conditional on HDL or conditional on the other traits investigated. Conjunction FDR identified altogether 17 loci; 10 loci were associated with PCA and LDL, 3 loci were associated with PCA and TG and additionally 4 loci were associated with PCA, LDL and TG jointly (conjunction FDR < 0.01). For T2D, we detected one locus adjacent to HNF1B. Conclusions We found polygenic overlap between PCA predisposition and blood lipids, in particular LDL and TG, and identified 17 pleiotropic gene loci between PCA and LDL, and PCA and TG, respectively. These findings provide novel pathobiological insights and may have implications for trials using targeting lipid-lowering agents in a prevention or cancer setting.

Mapping species differences for adventitious rooting in a Corymbia torelliana x Corymbia citriodora subspecies variegata hybrid

Quantitative trait loci (QTL) detection was carried out for adventitious rooting and associated propagation traits in a second-generation outbred Corymbia torelliana x Corymbia citriodora subspecies variegata hybrid family (n=186). The parental species of this cross are divergent in their capacity to develop roots adventitiously on stem cuttings and their propensity to form lignotubers. For the ten traits studied, there was one or two QTL detected, with some QTL explaining large amounts of phenotypic variation (e.g. 66% for one QTL for percentage rooting), suggesting that major effects influence rooting in this cross. Collocation of QTL for many strongly genetically correlated rooting traits to a single region on linkage group 12 suggested pleiotropy. A three locus model was most parsimonious for linkage group 12, however, as differences in QTL position and lower genetic correlations suggested separate loci for each of the traits of shoot production and root initiation. Species differences were thought to be the major source of phenotypic variation for some rooting rate and root quality traits because of the major QTL effects and up to 59-fold larger homospecific deviations (attributed to species differences) relative to heterospecific deviations (attributed to standing variation within species) evident at some QTL for these traits. A large homospecific/heterospecific ratio at major QTL suggested that the gene action evident in one cross may be indicative of gene action more broadly in hybrids between these species for some traits.

Genome-wide enrichment analysis between endometriosis and obesity-related traits reveals novel susceptibility loci

Endometriosis is a chronic inflammatory condition in women that results in pelvic pain and subfertility, and has been associated with decreased body mass index (BMI). Genetic variants contributing to the heritable component have started to emerge from genome-wide association studies (GWAS), although the majority remain unknown. Unexpectedly, we observed an intergenic locus on 7p15.2 that was genome-wide significantly associated with both endometriosis and fat distribution (waist-to-hip ratio adjusted for BMI; WHRadjBMI) in an independent meta-GWAS of European ancestry individuals. This led us to investigate the potential overlap in genetic variants underlying the aetiology of endometriosis, WHRadjBMI and BMI using GWAS data. Our analyses demonstrated significant enrichment of common variants between fat distribution and endometriosis (P = 3.7 x 10(-3)), which was stronger when we restricted the investigation to more severe (Stage B) cases (P = 4.5 x 10(-4)). However, no genetic enrichment was observed between endometriosis and BMI (P = 0.79). In addition to 7p15.2, we identify four more variants with statistically significant evidence of involvement in both endometriosis and WHRadjBMI (in/near KIFAP3, CAB39L, WNT4, GRB14); two of these, KIFAP3 and CAB39L, are novel associations for both traits. KIFAP3, WNT4 and 7p15.2 are associated with the WNT signalling pathway; formal pathway analysis confirmed a statistically significant (P = 6.41 x 10(-4)) overrepresentation of shared associations in developmental processes/WNT signalling between the two traits. Our results demonstrate an example of potential biological pleiotropy that was hitherto unknown, and represent an opportunity for functional follow-up of loci and further cross-phenotype comparisons to assess how fat distribution and endometriosis pathogenesis research fields can inform each other.

Quantitative trait loci detection and benefits from marker-assisted selection in dairy cattle

Knowing the chromosomal areas or actual genes affecting the traits under selection would add more information to be used in the selection decisions which would potentially lead to higher genetic response. The first objective of this study was to map quantitative trait loci (QTL) affecting economically important traits in the Finnish Ayrshire population. The second objective was to investigate the effects of using QTL information in marker-assisted selection (MAS) on the genetic response and the linkage disequilibrium between the different parts of the genome. Whole genome scans were carried out on a grand-daughter design with 12 half-sib families and a total of 493 sons. Twelve different traits were studied: milk yield, protein yield, protein content, fat yield, fat content, somatic cell score (SCS), mastitis treatments, other veterinary treatments, days open, fertility treatments, non-return rate, and calf mortality. The average spacing of the typed markers was 20 cM with 2 to 14 markers per chromosome. Associations between markers and traits were analyzed with multiple marker regression. Significance was determined by permutation and genome-wise P-values obtained by Bonferroni correction. The benefits from MAS were investigated by simulation: a conventional progeny testing scheme was compared to a scheme where QTL information was used within families to select among full-sibs in the male path. Two QTL on different chromosomes were modelled. The effects of different starting frequencies of the favourable alleles and different size of the QTL effects were evaluated. A large number of QTL, 48 in total, were detected at 5% or higher chromosome-wise significance. QTL for milk production were found on 8 chromosomes, for SCS on 6, for mastitis treatments on 1, for other veterinary treatments on 5, for days open on 7, for fertility treatments on 7, for calf mortality on 6, and for non-return rate on 2 chromosomes. In the simulation study the total genetic response was faster with MAS than with conventional selection and the advantage of MAS persisted over the studied generations. The rate of response and the difference between the selection schemes reflected clearly the changes in allele frequencies of the favourable QTL. The disequilibrium between the polygenes and QTL was always negative and it was larger with larger QTL size. The disequilibrium between the two QTL was larger with QTL of large effect and it was somewhat larger with MAS for scenarios with starting frequencies below 0.5 for QTL of moderate size and below 0.3 for large QTL. In conclusion, several QTL affecting economically important traits of dairy cattle were detected. Further studies are needed to verify these QTL, check their presence in the present breeding population, look for pleiotropy and fine map the most interesting QTL regions. The results of the simulation studies show that using MAS together with embryo transfer to pre-select young bulls within families is a useful approach to increase the genetic merit of the AI-bulls compared to conventional selection.

Abundant genetic overlap between blood lipids and immune-mediated diseases indicates shared molecular genetic mechanisms

Epidemiological studies suggest a relationship between blood lipids and immune-mediated diseases, but the nature of these associations is not well understood. We used genome-wide association studies (GWAS) to investigate shared single nucleotide polymorphisms (SNPs) between blood lipids and immune-mediated diseases. We analyzed data from GWAS (n~200,000 individuals), applying new False Discovery Rate (FDR) methods, to investigate genetic overlap between blood lipid levels [triglycerides (TG), low density lipoproteins (LDL), high density lipoproteins (HDL)] and a selection of archetypal immune-mediated diseases (Crohn's disease, ulcerative colitis, rheumatoid arthritis, type 1 diabetes, celiac disease, psoriasis and sarcoidosis). We found significant polygenic pleiotropy between the blood lipids and all the investigated immune-mediated diseases. We discovered several shared risk loci between the immune-mediated diseases and TG (n = 88), LDL (n = 87) and HDL (n = 52). Three-way analyses differentiated the pattern of pleiotropy among the immune-mediated diseases. The new pleiotropic loci increased the number of functional gene network nodes representing blood lipid loci by 40%. Pathway analyses implicated several novel shared mechanisms for immune pathogenesis and lipid biology, including glycosphingolipid synthesis (e.g. FUT2) and intestinal host-microbe interactions (e.g. ATG16L1). We demonstrate a shared genetic basis for blood lipids and immune-mediated diseases independent of environmental factors. Our findings provide novel mechanistic insights into dyslipidemia and immune-mediated diseases and may have implications for therapeutic trials involving lipid-lowering and anti-inflammatory agents.

Genetic Sharing with Cardiovascular Disease Risk Factors and Diabetes Reveals Novel Bone Mineral Density Loci

Bone Mineral Density (BMD) is a highly heritable trait, but genome-wide association studies have identified few genetic risk factors. Epidemiological studies suggest associations between BMD and several traits and diseases, but the nature of the suggestive comorbidity is still unknown. We used a novel genetic pleiotropy-informed conditional False Discovery Rate (FDR) method to identify single nucleotide polymorphisms (SNPs) associated with BMD by leveraging cardiovascular disease (CVD) associated disorders and metabolic traits. By conditioning on SNPs associated with the CVD-related phenotypes, type 1 diabetes, type 2 diabetes, systolic blood pressure, diastolic blood pressure, high density lipoprotein, low density lipoprotein, triglycerides and waist hip ratio, we identified 65 novel independent BMD loci (26 with femoral neck BMD and 47 with lumbar spine BMD) at conditional FDR < 0.01. Many of the loci were confirmed in genetic expression studies. Genes validated at the mRNA levels were characteristic for the osteoblast/osteocyte lineage, Wnt signaling pathway and bone metabolism. The results provide new insight into genetic mechanisms of variability in BMD, and a better understanding of the genetic underpinnings of clinical comorbidity.

Linking melanism to brain development: expression of a melanism-related gene in barn owl feather follicles covaries with sleep ontogeny.

BACKGROUND: Intra-specific variation in melanocyte pigmentation, common in the animal kingdom, has caught the eye of naturalists and biologists for centuries. In vertebrates, dark, eumelanin pigmentation is often genetically determined and associated with various behavioral and physiological traits, suggesting that the genes involved in melanism have far reaching pleiotropic effects. The mechanisms linking these traits remain poorly understood, and the potential involvement of developmental processes occurring in the brain early in life has not been investigated. We examined the ontogeny of rapid eye movement (REM) sleep, a state involved in brain development, in a wild population of barn owls (Tyto alba) exhibiting inter-individual variation in melanism and covarying traits. In addition to sleep, we measured melanistic feather spots and the expression of a gene in the feather follicles implicated in melanism (PCSK2). RESULTS: As in mammals, REM sleep declined with age across a period of brain development in owlets. In addition, inter-individual variation in REM sleep around this developmental trajectory was predicted by variation in PCSK2 expression in the feather follicles, with individuals expressing higher levels exhibiting a more precocial pattern characterized by less REM sleep. Finally, PCSK2 expression was positively correlated with feather spotting. CONCLUSIONS: We demonstrate that the pace of brain development, as reflected in age-related changes in REM sleep, covaries with the peripheral activation of the melanocortin system. Given its role in brain development, variation in nestling REM sleep may lead to variation in adult brain organization, and thereby contribute to the behavioral and physiological differences observed between adults expressing different degrees of melanism.

The role of Ken&Barbie in JAK/STAT signalling during development of Drosophila melanogaster

Cell-cell interactions during embryonic development are crucial in the co-ordination of growth, differentiation and maintenance of many different cell types. To achieve this co-ordination each cell must properly translate signals received from neighbouring cells, into spatially and temporally appropriate developmental responses. A surprisingly limited number of signal pathways are responsible for the differentiation of enormous variety of cell types. As a result, pathways are frequently 'reused' during development. Thus, in mammals the JAK/STAT pathway is required during early embryogenesis, mammary gland formation, hematopoiesis and, finally, plays a pivotal role in immune response. In the canonical way, the JAK/STAT pathway is represented by a transmembrane receptor associated with a Janus kinase (JAK), which upon stimulation by an extra-cellular ligand, phosphorylates itself, the receptor and, finally, the signal transducer and activator of transcription (STAT) molecules. Phosphorylated STATs dimerise and translocate to the nucleus where they activate transcription of target genes. The JAK/STAT pathway has been conserved throughout evolution, and all known components are present in the genome of Drosophila melanogaster. Besides hematopoietic and immunity functions, the pathway is also required during development for processes including embryonic segmentation, tracheal morphogenesis, posterior spiracle formation etc. This study describes Drosophila Ken&Barbie (Ken) as a selective regulator of JAK/STAT signalling. ken mutations identified in a screen for modulators of an eye overgrowth phenotype, caused by over-expression of the pathway ligand unpaired, also interact genetically with the pathway receptor domeless (dome) and the transcription factor stat92E. Over-expression of Ken can phenocopy developmental defects known to be caused by the loss of JAK/STAT signalling. These genetic interactions suggest that Ken may function as a negative regulator of the pathway. Ken has C-terminal Zn-finger domain, presumably for DNA binding, and N-terminal BTB/POZ domain, often found in transcriptional repressors. Using EGFP-fused construct expressed in vivo revealed nuclear accumulation of Ken. Therefore, it is proposed that Ken may act as a suppresser of STAT92E target genes. An in vitro assay, termed SELEX, determined that Ken specifically binds to a DNA sequence, with the essential for DNA recognition core overlapping that of STAT92E. This interesting observation suggests that not all STAT92E sites may also allow Ken binding. Strikingly, when effects of ectopic Ken on the expression of putative JAK/STAT pathway target genes were examined, only a subset of the genes tested, namely vvl, trh and kni, were down-regulated by Ken, whereas some others, such as eve and fj, appeared to be unresponsive. Further analysis of vvl, one of the genes susceptible to ectopic Ken, was undertaken. In the developing hindgut, expression of vvl is JAK/STAT pathway dependent, but remains repressed in the posterior spiracles, despite the stimulation of STAT92E by Upd in their primordia. Importantly, ken is also expressed in the developing posterior spiracles. Strikingly, up-regulation of vvl is observed in these tissues in ken mutant embryos. These imply that while ectopic Ken is sufficient to repress the expression of vvl in the hindgut, endogenous Ken is also necessary to prevent its activation in the posterior spiracles. It is therefore conceivable that ectopic vvl expression in the posterior spiracles of the ken mutants may be the result of de-repression of endogenous STAT92E activity. Another consequence of these observations is a fine balance that must exist between STAT92E and Ken activities. Apparently, endogenous level of Ken is sufficient to repress vvl, but not other, as yet unidentified, JAK/STAT pathway targets, whose presumable activation by STAT92E is required for posterior spiracle development as the embryos mutant for dome, the receptor of the pathway, show severe spiracle defects. These defects are also observed in the embryos mis-expressing Ken. Though it is possible that the posterior spiracle phenotype caused by higher levels of Ken results from a JAK/STAT pathway independent activity, it seems to be more likely that Ken acts in a dosage dependent manner, and extra Ken is able to further antagonise JAK/STAT pathway target genes. While STAT92E binding sites required for target gene expression have been poorly characterised, the existence of genome data allows the prediction of candidate STAT92E sites present in target genes promoters to be attempted. When a 6kb region containing the putative regulatory domains flanking the vvl locus are examined, only a single potential STAT92E binding site located 825bp upstream of the translational start can be detected. Strikingly, this site also includes a perfect Ken binding sequence. Such an in silico observation, though consistent with both Ken DNA binding assay in vitro and regulation of STAT92E target genes in vivo, however, requires further analysis. The JAK/STAT pathway is implicated in a variety of processes during embryonic and larval development as well as in imago. In each case, stimulation of the same transcription factor results in different developmental outcomes. While many potential mechanisms have been proposed and demonstrated to explain such pleiotropy, the present study indicates that Ken may represent another mechanism, with which signal transduction pathways are controlled. Ken selectively down-regulates a subset of potential target genes and so modifies the transcriptional profile generated by activated STAT92E - a mechanism, which may be partially responsible for differences in the morphogenetic processes elicited by JAK/STAT signalling during development.