Common polygenic variation contributes to risk of migraine in the Norfolk Island population

Migraine has been defined as a common disabling primary headache disorder. Epidemiology studies have provided with the undeniable evidence of genetic components as active players in the development of the disease under a polygenic model in which multiple risk alleles exert modest individual effects. Our objective was to test the contribution of a polygenic effect to migraine risk in the Norfolk Island population using a panel of SNPs reported to be disease associated in published migraine GWAS. We also investigated whether individual SNPs were associated with gene expression levels measured in whole-blood. Polygenic scores were calculated in a total of 285 related individuals (74 cases, 211 controls) from the Norfolk Island using 51 SNPs previously reported to be associated with migraine in published GWAS. The association between polygenic score and migraine case-control status was tested using logistic regression. Results indicate that a migraine polygenic risk score was associated with migraine case-control status in this population (P=0.016). This supports the hypothesis that multiple SNPs with weak effects collectively contribute to migraine risk in this population. Amongst the SNPs included in the polygenic model, 4 were associated with the expression of the USMG5 gene, including rs171251 (P = 0.012). Results from this study provide evidence for a polygenic contribution to migraine risk in an isolated population and highlight specific SNPs that regulate the expression of USMG5, a gene critical for mitochondrial function.

Genetic studies of common rheumatological diseases

Editorial

Bahia grass pollen specific IgE is common in seasonal rhinitis patients but has limited cross-reactivity with Ryegrass

Background: Perennial Ryegrass is a major cause of rhinitis in spring and early summer. Bahia grass, Paspalum notatum, flowers late into summer and could account for allergic rhinitis at this time. We determined the frequency of serum immunoglobulin (Ig)E reactivity with Bahia grass in Ryegrass pollen allergic patients and investigated IgE cross-reactivity between Bahia and Ryegrass. Methods: Serum from 33 Ryegrass pollen allergic patients and 12 nonatopic donors were tested for IgE reactivity with Bahia and Ryegrass pollen extracts (PE) by enzyme-linked immunosorbent assay (ELISA), western blotting and inhibition ELISA. Allergen-specific antibodies from a pool of sera from allergic donors were affinity purified and tested for IgE cross-reactivity. Results: Seventy-eight per cent of the sera had IgE reactivity with Bahia grass, but more weakly than with Ryegrass. Antibodies eluted from the major Ryegrass pollen allergens, Lol p 1 and Lol p 5, showed IgE reactivity with allergens of Ryegrass and Canary but not Bahia or Bermuda grasses. Timothy, Canary and Ryegrass inhibited IgE reactivity with Ryegrass and Bahia grass, whereas Bahia, Johnson and Bermuda grass did not inhibit IgE reactivity with Ryegrass. Conclusions: The majority of Ryegrass allergic patients also showed serum IgE reactivity with Bahia grass PE. However, Bahia grass and Ryegrass had only limited IgE cross-reactivity indicating that Bahia grass should be considered in diagnosis and treatment of patients with hay fever late in' the grass pollen season.

A novel LEMD3 mutation common to patients with osteopoikilosis with and without melorheostosis

Recent studies have reported loss of function mutations in the LEMD3 gene, encoding an inner nuclear membrane protein that influences Smad signaling, as a cause of osteopoikilosis, Buschke-Ollendorff syndrome, and melorheostosis. We investigated LEMD3 in a three-generation family with osteopoikilosis from the Azores, an affected father and daughter from Ireland with osteopoikilosis (the daughter also had melorheostosis), and two other individuals from the UK with isolated melorheostosis. We found a novel C to T substitution at position 2032 bp (cDNA) in exon 8 of LEMD3, resulting in a premature stop codon at amino acid position 678. This mutation co-segregates with the osteopoikilosis phenotype in both the Azorean family and the Irish family. It was not detected in any of the six unaffected family members or in 342 healthy Caucasian individuals. No LEMD3 mutations were detected in the two patients with sporadic melorheostosis. The LEMD3 mutation reported was clearly the cause of osteopoikilosis in the two families but its relationship to melorheostosis in one of the family members is still unclear. Perhaps unsurprisingly in what is a segmental disease, we did not find LEMD3 mutations in peripheral-blood-derived DNA from the two other individuals with sporadic melorheostosis. The nature of the additional genetic and/or environmental influences required for the development of melorheostosis in those with osteopoikilosis requires further investigation.

Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis

To identify susceptibility loci for visceral leishmaniasis, we undertook genome-wide association studies in two populations: 989 cases and 1,089 controls from India and 357 cases in 308 Brazilian families (1,970 individuals). The HLA-DRB1-HLA-DQA1 locus was the only region to show strong evidence of association in both populations. Replication at this region was undertaken in a second Indian population comprising 941 cases and 990 controls, and combined analysis across the three cohorts for rs9271858 at this locus showed P combined = 2.76 × 10 -17 and odds ratio (OR) = 1.41, 95% confidence interval (CI) = 1.30-1.52. A conditional analysis provided evidence for multiple associations within the HLA-DRB1-HLA-DQA1 region, and a model in which risk differed between three groups of haplotypes better explained the signal and was significant in the Indian discovery and replication cohorts. In conclusion, the HLA-DRB1-HLA-DQA1 HLA class II region contributes to visceral leishmaniasis susceptibility in India and Brazil, suggesting shared genetic risk factors for visceral leishmaniasis that cross the epidemiological divides of geography and parasite species. © 2013 Nature America, Inc. All rights reserved.

Genome-wide association study to identify genetic determinants of severe asthma

Background: The genetic basis for developing asthma has been extensively studied. However, association studies to date have mostly focused on mild to moderate disease and genetic risk factors for severe asthma remain unclear. Objective: To identify common genetic variants affecting susceptibility to severe asthma. Methods: A genome-wide association study was undertaken in 933 European ancestry individuals with severe asthma based on Global Initiative for Asthma (GINA) criteria 3 or above and 3346 clean controls. After standard quality control measures, the association of 480 889 genotyped single nucleotide polymorphisms (SNPs) was tested. To improve the resolution of the association signals identified, non-genotyped SNPs were imputed in these regions using a dense reference panel of SNP genotypes from the 1000 Genomes Project. Then replication of SNPs of interest was undertaken in a further 231 cases and 1345 controls and a meta-analysis was performed to combine the results across studies. Results: An association was confirmed in subjects with severe asthma of loci previously identified for association with mild to moderate asthma. The strongest evidence was seen for the ORMDL3/GSDMB locus on chromosome 17q12-21 (rs4794820, p=1.03×10 (-8)following meta-analysis) meeting genome-wide significance. Strong evidence was also found for the IL1RL1/IL18R1 locus on 2q12 (rs9807989, p=5.59×10 (-8) following meta-analysis) just below this threshold. No novel loci for susceptibility to severe asthma met strict criteria for genome-wide significance. Conclusions: The largest genome-wide association study of severe asthma to date was carried out and strong evidence found for the association of two previously identified asthma susceptibility loci in patients with severe disease. A number of novel regions with suggestive evidence were also identified warranting further study.

Common variants at the MHC locus and at chromosome 16q24.1 predispose to Barrett's Esophagus

Barrett's esophagus is an increasingly common disease that is strongly associated with reflux of stomach acid and usually a hiatus hernia, and it strongly predisposes to esophageal adenocarcinoma (EAC), a tumor with a very poor prognosis. We report the first genome-wide association study on Barrett's esophagus, comprising 1,852 UK cases and 5,172 UK controls in the discovery stage and 5,986 cases and 12,825 controls in the replication stage. Variants at two loci were associated with disease risk: chromosome 6p21, rs9257809 (P combined = 4.09 × 10-9; odds ratio (OR) = 1.21, 95% confidence interval (CI) =1.13-1.28), within the major histocompatibility complex locus, and chromosome 16q24, rs9936833 (P combined = 2.74 × 10-10; OR = 1.14, 95% CI = 1.10-1.19), for which the closest protein-coding gene is FOXF1, which is implicated in esophageal development and structure. We found evidence that many common variants of small effect contribute to genetic susceptibility to Barrett's esophagus and that SNP alleles predisposing to obesity also increase risk for Barrett's esophagus. © 2012 Nature America, Inc. All rights reserved.

Common variants at 12q14 and 12q24 are associated with hippocampal volume

Aging is associated with reductions in hippocampal volume that are accelerated by Alzheimer's disease and vascular risk factors. Our genome-wide association study (GWAS) of dementia-free persons (n = 9,232) identified 46 SNPs at four loci with P values of <4.0 × 10 -7. In two additional samples (n = 2,318), associations were replicated at 12q14 within MSRB3-WIF1 (discovery and replication; rs17178006; P = 5.3 × 10 -11) and at 12q24 near HRK-FBXW8 (rs7294919; P = 2.9 × 10 -11). Remaining associations included one SNP at 2q24 within DPP4 (rs6741949; P = 2.9 × 10 -7) and nine SNPs at 9p33 within ASTN2 (rs7852872; P = 1.0 × 10 -7); along with the chromosome 12 associations, these loci were also associated with hippocampal volume (P < 0.05) in a third younger, more heterogeneous sample (n = 7,794). The SNP in ASTN2 also showed suggestive association with decline in cognition in a largely independent sample (n = 1,563). These associations implicate genes related to apoptosis (HRK), development (WIF1), oxidative stress (MSR3B), ubiquitination (FBXW8) and neuronal migration (ASTN2), as well as enzymes targeted by new diabetes medications (DPP4), indicating new genetic influences on hippocampal size and possibly the risk of cognitive decline and dementia.

Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes

Metformin is the most commonly used pharmacological therapy for type 2 diabetes. We report a genome-wide association study for glycemic response to metformin in 1,024 Scottish individuals with type 2 diabetes with replication in two cohorts including 1,783 Scottish individuals and 1,113 individuals from the UK Prospective Diabetes Study. In a combined meta-analysis, we identified a SNP, rs11212617, associated with treatment success (n = 3,920, P = 2.9 P×-9, odds ratio = 1.35, 95% CI 1.22-1.49) at a locus containing ATM, the ataxia telangiectasia mutated gene. In a rat hepatoma cell line, inhibition of ATM with KU-55933 attenuated the phosphorylation and activation of AMP-activated protein kinase in response to metformin. We conclude that ATM, a gene known to be involved in DNA repair and cell cycle control, plays a role in the effect of metformin upstream of AMP-activated protein kinase, and variation in this gene alters glycemic response to metformin. © 2011 Nature America, Inc. All rights reserved.

Dissection of the genetics of Parkinson's disease identifies an additional association 5' of SNCA and multiple associated haplotypes at 17q21

We performed a genome-wide association study (GWAS) in 1705 Parkinson's disease (PD) UK patients and 5175 UK controls, the largest sample size so far for a PD GWAS. Replication was attempted in an additional cohort of 1039 French PD cases and 1984 controls for the 27 regions showing the strongest evidence of association (P < 10 4). We replicated published associations in the 4q22/SNCA and 17q21/MAPT chromosome regions (P < 10 10) and found evidence for an additional independent association in 4q22/SNCA.A detailed analysis of the haplotype structure at 17q21 showed that there are three separate risk groups within this region. We found weak but consistent evidence of association for common variants located in three previously published associated regions (4p15/BST1, 4p16/GAK and 1q32/PARK16). We found no support for the previously reported SNP association in 12q12/LRRK2. We also found an association of the two SNPs in 4q22/SNCA with the age of onset of the disease. © The Author 2010. Published by Oxford University Press.

Genome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls

Copy number variants (CNVs) account for a major proportion of human genetic polymorphism and have been predicted to have an important role in genetic susceptibility to common disease. To address this we undertook a large, direct genome-wide study of association between CNVs and eight common human diseases. Using a purpose-designed array we typed 19,000 individuals into distinct copy-number classes at 3,432 polymorphic CNVs, including an estimated 50% of all common CNVs larger than 500 base pairs. We identified several biological artefacts that lead to false-positive associations, including systematic CNV differences between DNAs derived from blood and cell lines. Association testing and follow-up replication analyses confirmed three loci where CNVs were associated with diseaseIRGM for Crohns disease, HLA for Crohns disease, rheumatoid arthritis and type 1 diabetes, and TSPAN8 for type 2 diabetesalthough in each case the locus had previously been identified in single nucleotide polymorphism (SNP)-based studies, reflecting our observation that most common CNVs that are well-typed on our array are well tagged by SNPs and so have been indirectly explored through SNP studies. We conclude that common CNVs that can be typed on existing platforms are unlikely to contribute greatly to the genetic basis of common human diseases. © 2010 Macmillan Publishers Limited. All rights reserved.

Common variants in the region around Osterix are associated with bone mineral density and growth in childhood

Peak bone mass achieved in adolescence is a determinant of bone mass in later life. In order to identify genetic variants affecting bone mineral density (BMD), we performed a genome-wide association study of BMD and related traits in 1518 children from the Avon Longitudinal Study of Parents and Children (ALSPAC). We compared results with a scan of 134 adults with high or low hip BMD. We identified associations with BMD in an area of chromosome 12 containing the Osterix (SP7) locus, a transcription factor responsible for regulating osteoblast differentiation (ALSPAC: P = 5.8 × 10-4; Australia: P = 3.7 × 10-4). This region has previously shown evidence of association with adult hip and lumbar spine BMD in an Icelandic population, as well as nominal association in a UK population. A meta-analysis of these existing studies revealed strong association between SNPs in the Osterix region and adult lumbar spine BMD (P = 9.9 × 10-11). In light of these findings, we genotyped a further 3692 individuals from ALSPAC who had whole body BMD and confirmed the association in children as well (P = 5.4 × 10-5). Moreover, all SNPs were related to height in ALSPAC children, but not weight or body mass index, and when height was included as a covariate in the regression equation, the association with total body BMD was attenuated. We conclude that genetic variants in the region of Osterix are associated with BMD in children and adults probably through primary effects on growth.

Convergent evolution of marine mammals is associated with distinct substitutions in common genes

Phenotypic convergence is thought to be driven by parallel substitutions coupled with natural selection at the sequence level. Multiple independent evolutionary transitions of mammals to an aquatic environment offer an opportunity to test this thesis. Here, whole genome alignment of coding sequences identified widespread parallel amino acid substitutions in marine mammals; however, the majority of these changes were not unique to these animals. Conversely, we report that candidate aquatic adaptation genes, identified by signatures of likelihood convergence and/or elevated ratio of nonsynonymous to synonymous nucleotide substitution rate, are characterized by very few parallel substitutions and exhibit distinct sequence changes in each group. Moreover, no significant positive correlation was found between likelihood convergence and positive selection in all three marine lineages. These results suggest that convergence in protein coding genes associated with aquatic lifestyle is mainly characterized by independent substitutions and relaxed negative selection.

Rapid assessment of genotype-by-environment interactions and heritability for growth rate in aquaculture species using in vitro fertilisation and DNA tagging

Commercial environments may receive only a fraction of expected genetic gains for growth rate as predicted from the selection environment This fraction is the result of undesirable genotype-by-environment interactions (G x E) and measured by the genetic correlation (r(g)) of growth between environments. Rapid estimates of genetic correlation achieved in one generation are notoriously difficult to estimate with precision. A new design is proposed where genetic correlations can be estimated by utilising artificial mating from cryopreserved semen and unfertilised eggs stripped from a single female. We compare a traditional phenotype analysis of growth to a threshold model where only the largest fish are genotyped for sire identification. The threshold model was robust to differences in family mortality differing up to 30%. The design is unique as it negates potential re-ranking of families caused by an interaction between common maternal environmental effects and growing environment. The design is suitable for rapid assessment of G x E over one generation with a true 0.70 genetic correlation yielding standard errors as low as 0.07. Different design scenarios were tested for bias and accuracy with a range of heritability values, number of half-sib families created, number of progeny within each full-sib family, number of fish genotyped, number of fish stocked, differing family survival rates and at various simulated genetic correlation levels

Population structure, mortality and growth of Pinna nobilis Linnaeus, 1758 (Mollusca, Bivalvia) at different depths in Moraira bay (Alicante, Western Mediterranean)

An investigation to characterize the causes of Pinna nobilis population structure in Moraira bay (Western Mediterranean) was developed. Individuals of two areas of the same Posidonia meadow, located at different depths (A1, -13 and A2, -6 m), were inventoried, tagged, their positions accurately recorded and monitored from July 1997 to July 2002. On each area, different aspects of population demography were studied (i.e. spatial distribution, size structure, displacement evidences, mortality, growth and shell orientation). A comparison between both groups of individuals was carried out, finding important differences between them. In A1, the individuals were more aggregated and mean and maximum size were higher (A1, 10.3 and A2, 6 individuals/100 m(2); A1, x = 47.2 +/- 9.9; A2, x = 29.8 +/- 7.4 cm, P < 0.001, respectively). In A2, growth rate and mortality were higher, the latter concentrated on the largest individuals, in contrast to A1, where the smallest individuals had the higher mortality rate [A1, L = 56.03(1 - e(-0.17t)); A2, L = 37.59(1 - e(-0.40t)), P < 0.001; mean annual mortality A1: 32 dead individuals out of 135, 23.7% and A2: 16 dead individuals out of 36, 44.4%, and total mortality coefficients (z), z(A1(-30)) = 0.28, z(A1(31-45)) = 0.05, z(A1(46-)) = 0.08; z(A2(-30)) = 0.15, z(A2(31-45)) = 0.25]. A common shell orientation N-S, coincident with the maximum shore exposure, was observed in A2. Spatial distribution in both areas showed not enough evidence to discard a random distribution of the individuals, despite the greater aggregation on the deeper area (A1) (A1, chi(2) = 0.41, df = 3, P > 0.5, A2, chi(2)= 0.98, df = 2 and 0.3 < P < 0.5). The obtained results have demonstrated that the depth-related size segregation usually shown by P. nobilis is mainly caused by differences in mortality and growth among individuals located at different depths, rather than by the active displacement of individuals previously reported in the literature. Furthermore, dwarf individuals are observed in shallower levels and as a consequence, the relationship between size and age are not comparable even among groups of individuals inhabiting the same meadow at different depths. The final causes of the differences on mortality and growth are also discussed.