A candidate gene association study for growth performance in an improved Giant Freshwater Prawn (Macrobrachium rosenbergii) culture line

A candidate gene approach using type I single nucleotide polymorphism (SNP) markers can provide an effective method for detecting genes and gene regions that underlie phenotypic variation in adaptively significant traits. In the absence of available genomic data resources, transcriptomes were recently generated in Macrobrachium rosenbergii to identify candidate genes and markers potentially associated with growth. The characterisation of 47 candidate loci by ABI re-sequencing of four cultured and eight wild samples revealed 342 putative SNPs. Among these, 28 SNPs were selected in 23 growth-related candidate genes to genotype in 200 animals selected for improved growth performance in an experimental GFP culture line in Vietnam. The associations between SNP markers and individual growth performance were then examined. For additive and dominant effects, a total of three exonic SNPs in glycogen phosphorylase (additive), heat shock protein 90 (additive and dominant) and peroxidasin (additive), and a total of six intronic SNPs in ankyrin repeats-like protein (additive and dominant), rolling pebbles (dominant), transforming growth factor-β induced precursor (dominant), and UTP-glucose-1-phosphate uridylyltransferase 2 (dominant) genes showed significant associations with the estimated breeding values in the experimental animals (P =0.001−0.031). Individually, they explained 2.6−4.8 % of the genetic variance (R2=0.026−0.048). This is the first large set of SNP markers reported for M. rosenbergii and will be useful for confirmation of associations in other samples or culture lines as well as having applications in marker-assisted selection in future breeding programs.

Identity-by-Descent Mapping to Detect Rare Variants Conferring Susceptibility to Multiple Sclerosis

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. 'Population-based linkage analysis' (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10-6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority. © 2013 Lin et al.

Association study of genes related to bone formation and resorption and the extent of radiographic change in ankylosing spondylitis

Objective: To identify genetic associations with severity of radiographic damage in ankylosing spondylitis (AS). Method: We studied 1537 AS cases of European descent; all fulfilled the modified New York Criteria. Radiographic severity was assessed from digitised lateral radiographs of the cervical and lumbar spine using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS). A two-phase genotyping design was used. In phase 1, 498 single nucleotide polymorphisms (SNPs) were genotyped in 688 cases; these were selected to capture >90% of the common haplotypic variation in the exons, exon-intron boundaries, and 5 kb flanking DNA in the 5' and 3' UTR of 74 genes involved in anabolic or catabolic bone pathways. In phase 2, 15 SNPs exhibiting p<0.05 were genotyped in a further cohort of 830 AS cases; results were analysed both separately and in combination with the discovery phase data. Association was tested by contingency tables after separating the samples into 'mild' and 'severe' groups, defined as the bottom and top 40% by mSASSS, adjusted for gender and disease duration. Results: Experiment-wise association was observed with the SNP rs8092336 (combined OR 0.32, p=1.2×10-5), which lies within RANK (receptor activator of NF?B), a gene involved in osteoclastogenesis, and in the interaction between T cells and dendritic cells. Association was also found with the SNP rs1236913 in PTGS1 (prostaglandin-endoperoxide synthase 1, cyclooxygenase 1), giving an OR of 0.53 (p=2.6×10-3). There was no observed association between radiographic severity and HLA-B*27. Conclusions: These findings support roles for bone resorption and prostaglandins pathways in the osteoproliferative changes in AS.

The multiple sclerosis whole blood mRNA transcriptome and genetic associations indicate dysregulation of specific T cell pathways in pathogenesis

Multiple sclerosis (MS) is an autoimmune disease with a genetic component, caused at least in part by aberrant lymphocyte activity. The whole blood mRNA transcriptome was measured for 99 untreated MS patients: 43 primary progressive MS, 20 secondary progressive MS, 36 relapsing remitting MS and 45 age-matched healthy controls. The ANZgene Multiple Sclerosis Genetics Consortium genotyped more than 300 000 SNPs for 115 of these samples. Transcription from genes on translational regulation, oxidative phosphorylation, immune synapse and antigen presentation pathways was markedly increased in all forms of MS. Expression of genes tagging T cells was also upregulated (P < 10-12) in MS. A T cell gene signature predicts disease state with a concordance index of 0.79 with age and gender as co-variables, but the signature is not associated with clinical course or disability. The ANZgene genome wide association screen identified two novel regions with genome wide significance: one encoding the T cell co-stimulatory molecule, CD40; the other a region on chromosome 12q13-14. The CD40 haplotype associated with increased MS susceptibility has decreased gene expression in MS (P < 0.0007). The second MS susceptibility region includes 17 genes on 12q13-14 in tight linkage disequilibrium. Of these, only 13 are expressed in leukocytes, and of these the expression of one, FAM119B, is much lower in the susceptibility haplotype (P tdthomlt; 10-14). Overall, these data indicate dysregulation of T cells can be detected in the whole blood of untreated MS patients, and supports targeting of activated T cells in therapy for all forms of MS.

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same region.

A large-scale analysis of genetic variants within putative miRNA binding sites in prostate cancer

Prostate cancer is the second most common malignancy among men worldwide. Genome-wide association studies have identified 100 risk variants for prostate cancer, which can explain approximately 33% of the familial risk of the disease. We hypothesized that a comprehensive analysis of genetic variations found within the 3' untranslated region of genes predicted to affect miRNA binding (miRSNP) can identify additional prostate cancer risk variants. We investigated the association between 2,169 miRSNPs and prostate cancer risk in a large-scale analysis of 22,301 cases and 22,320 controls of European ancestry from 23 participating studies. Twenty-two miRSNPs were associated (P<2.3×10(-5)) with risk of prostate cancer, 10 of which were within 7 genes previously not mapped by GWAS studies. Further, using miRNA mimics and reporter gene assays, we showed that miR-3162-5p has specific affinity for the KLK3 rs1058205 miRSNP T-allele, whereas miR-370 has greater affinity for the VAMP8 rs1010 miRSNP A-allele, validating their functional role. SIGNIFICANCE Findings from this large association study suggest that a focus on miRSNPs, including functional evaluation, can identify candidate risk loci below currently accepted statistical levels of genome-wide significance. Studies of miRNAs and their interactions with SNPs could provide further insights into the mechanisms of prostate cancer risk.

Generalizability of established prostate cancer risk variants in men of African ancestry

Genome-wide association studies have identified more than 80 risk variants for prostate cancer, mainly in European or Asian populations. The generalizability of these variants in other racial/ethnic populations needs to be understood before the loci can be used widely in risk modeling. In our study, we examined 82 previously reported risk variants in 4,853 prostate cancer cases and 4,678 controls of African ancestry. We performed association testing for each variant using logistic regression adjusted for age, study and global ancestry. Of the 82 known risk variants, 68 (83%) had effects that were directionally consistent in their association with prostate cancer risk and 30 (37%) were significantly associated with risk at p < 0.05, with the most statistically significant variants being rs116041037 (p = 3.7 × 10(-26) ) and rs6983561 (p = 1.1 × 10(-16) ) at 8q24, as well as rs7210100 (p = 5.4 × 10(-8) ) at 17q21. By exploring each locus in search of better markers, the number of variants that captured risk in men of African ancestry (p < 0.05) increased from 30 (37%) to 44 (54%). An aggregate score comprised of these 44 markers was strongly associated with prostate cancer risk [per-allele odds ratio (OR) = 1.12, p = 7.3 × 10(-98) ]. In summary, the consistent directions of effects for the vast majority of variants in men of African ancestry indicate common functional alleles that are shared across populations. Further exploration of these susceptibility loci is needed to identify the underlying biologically relevant variants to improve prostate cancer risk modeling in populations of African ancestry.

A Phenomic Scan of the Norfolk Island Genetic Isolate Identifies a Major Pleiotropic Effect Locus Associated with Metabolic and Renal Disorder Markers

Multiphenotype genome-wide association studies (GWAS) may reveal pleiotropic genes, which would remain undetected using single phenotype analyses. Analysis of large pedigrees offers the added advantage of more accurately assessing trait heritability, which can help prioritise genetically influenced phenotypes for GWAS analysis. In this study we performed a principal component analysis (PCA), heritability (h2) estimation and pedigree-based GWAS of 37 cardiovascular disease -related phenotypes in 330 related individuals forming a large pedigree from the Norfolk Island genetic isolate. PCA revealed 13 components explaining >75% of the total variance. Nine components yielded statistically significant h2 values ranging from 0.22 to 0.54 (P<0.05). The most heritable component was loaded with 7 phenotypic measures reflecting metabolic and renal dysfunction. A GWAS of this composite phenotype revealed statistically significant associations for 3 adjacent SNPs on chromosome 1p22.2 (P<1x10-8). These SNPs form a 42kb haplotype block and explain 11% of the genetic variance for this renal function phenotype. Replication analysis of the tagging SNP (rs1396315) in an independent US cohort supports the association (P = 0.000011). Blood transcript analysis showed 35 genes were associated with rs1396315 (P<0.05). Gene set enrichment analysis of these genes revealed the most enriched pathway was purine metabolism (P = 0.0015). Overall, our findings provide convincing evidence for a major pleiotropic effect locus on chromosome 1p22.2 influencing risk of renal dysfunction via purine metabolism pathways in the Norfolk Island population. Further studies are now warranted to interrogate the functional relevance of this locus in terms of renal pathology and cardiovascular disease risk.

Genetic analysis for a shared biological basis between migraine and coronary artery disease

Objective: To apply genetic analysis of genome-wide association data to study the extent and nature of a shared biological basis between migraine and coronary artery disease (CAD). Methods: Four separate methods for cross-phenotype genetic analysis were applied on data from 2 large-scale genome-wide association studies of migraine (19,981 cases, 56,667 controls) and CAD (21,076 cases, 63,014 controls). The first 2 methods quantified the extent of overlapping risk variants and assessed the load of CAD risk loci in migraineurs. Genomic regions of shared risk were then identified by analysis of covariance patterns between the 2 phenotypes and by querying known genome-wide significant loci. Results: We found a significant overlap of genetic risk loci for migraine and CAD. When stratified by migraine subtype, this was limited to migraine without aura, and the overlap was protective in that patients with migraine had a lower load of CAD risk alleles than controls. Genes indicated by 16 shared risk loci point to mechanisms with potential roles in migraine pathogenesis and CAD, including endothelial dysfunction (PHACTR1) and insulin homeostasis (GIP). Conclusions: The results suggest that shared biological processes contribute to risk of migraine and CAD, but surprisingly this commonality is restricted to migraine without aura and the impact is in opposite directions. Understanding the mechanisms underlying these processes and their opposite relationship to migraine and CAD may improve our understanding of both disorders.

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.

Concordance of genetic risk across migraine subgroups: Impact on current and future genetic association studies

BACKGROUND There has been intensive debate whether migraine with aura (MA) and migraine without aura (MO) should be considered distinct subtypes or part of the same disease spectrum. There is also discussion to what extent migraine cases collected in specialised headache clinics differ from cases from population cohorts, and how female cases differ from male cases with respect to their migraine. To assess the genetic overlap between these migraine subgroups, we examined genome-wide association (GWA) results from analysis of 23,285 migraine cases and 95,425 population-matched controls. METHODS Detailed heterogeneity analysis of single-nucleotide polymorphism (SNP) effects (odds ratios) between migraine subgroups was performed for the 12 independent SNP loci significantly associated (p < 5 x 10(-8); thus surpassing the threshold for genome-wide significance) with migraine susceptibility. Overall genetic overlap was assessed using SNP effect concordance analysis (SECA) at over 23,000 independent SNPs. RESULTS: Significant heterogeneity of SNP effects (p het < 1.4 x 10(-3)) was observed between the MA and MO subgroups (for SNP rs9349379), and between the clinic- and population-based subgroups (for SNPs rs10915437, rs6790925 and rs6478241). However, for all 12 SNPs the risk-increasing allele was the same, and SECA found the majority of genome-wide SNP effects to be in the same direction across the subgroups. CONCLUSIONS Any differences in common genetic risk across these subgroups are outweighed by the similarities. Meta-analysis of additional migraine GWA datasets, regardless of their major subgroup composition, will identify new susceptibility loci for migraine.

DCAF4, a novel gene associated with leucocyte telomere length

Background Leucocyte telomere length (LTL), which is fashioned by multiple genes, has been linked to a host of human diseases, including sporadic melanoma. A number of genes associated with LTL have already been identified through genome-wide association studies. The main aim of this study was to establish whether DCAF4 (DDB1 and CUL4-associated factor 4) is associated with LTL. In addition, using ingenuity pathway analysis (IPA), we examined whether LTL-associated genes in the general population might partially explain the inherently longer LTL in patients with sporadic melanoma, the risk for which is increased with ultraviolet radiation (UVR). Results Genome-wide association (GWA) meta-analysis and de novo genotyping of 20 022 individuals revealed a novel association (p=6.4×10−10) between LTL and rs2535913, which lies within DCAF4. Notably, eQTL analysis showed that rs2535913 is associated with decline in DCAF4 expressions in both lymphoblastoid cells and sun-exposed skin (p=4.1×10−3 and 2×10−3, respectively). Moreover, IPA revealed that LTL-associated genes, derived from GWA meta-analysis (N=9190), are over-represented among genes engaged in melanoma pathways. Meeting increasingly stringent p value thresholds (p<0.05, <0.01, <0.005, <0.001) in the LTL-GWA meta-analysis, these genes were jointly over-represented for melanoma at p values ranging from 1.97×10−169 to 3.42×10−24. Conclusions We uncovered a new locus associated with LTL in the general population. We also provided preliminary findings that suggest a link of LTL through genetic mechanisms with UVR and melanoma in the general population.

Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine

Migraine is a debilitating neurological disorder affecting around 1 in 7 people worldwide, but its molecular mechanisms remain poorly understood. Some debate exists over whether migraine is a disease of vascular dysfunction, or a result of neuronal dysfunction with secondary vascular changes. Genome-wide association (GWA) studies have thus far identified 13 independent loci associated with migraine. To identify new susceptibility loci, we performed the largest genetic study of migraine to date, comprising 59,674 cases and 316,078 controls from 22 GWA studies. We identified 45 independent single nucleotide polymorphisms (SNPs) significantly associated with migraine risk (P < 5 x 10-8) that map to 38 distinct genomic loci, including 28 loci not previously reported and the first locus identified on chromosome X. Furthermore, a subset analysis for migraine without aura (MO) identified seven of the same loci as from the full sample, whereas no loci reached genome-wide significance in the migraine with aura (MA) subset. In subsequent computational analyzes, the identified loci showed enrichment for genes expressed in vascular and smooth muscle tissues, consistent with a predominant theory of migraine that highlights vascular etiologies.

Functional evaluation of genetic variants associated with endometriosis near GREB1

STUDY QUESTION: Do DNA variants in the growth regulation by estrogen in breast cancer 1 (GREB1) region regulate endometrial GREB1 expression and increase the risk of developing endometriosis in women? SUMMARY ANSWER: We identified new single nucleotide polymorphisms (SNPs) with strong association with endometriosis at the GREB1 locus although we did not detect altered GREB1 expression in endometriosis patients with defined genotypes. WHAT IS ALREADY KNOWN: Genome-wide association studies have identified the GREB1 region on chromosome 2p25.1 for increasing endometriosis risk. The differential expression of GREB1 has also been reported by others in association with endometriosis disease phenotype. STUDY DESIGN, SIZE, DURATION: Fine mapping studies comprehensively evaluated SNPs within the GREB1 region in a large-scale data set (>2500 cases and >4000 controls). Publicly available bioinformatics tools were employed to functionally annotate SNPs showing the strongest association signal with endometriosis risk. Endometrial GREB1 mRNA and protein expression was studied with respect to phases of the menstrual cycle (n = 2-45 per cycle stage) and expression quantitative trait loci (eQTL) analysis for significant SNPs were undertaken for GREB1 [mRNA (n = 94) and protein (n = 44) in endometrium]. PARTICIPANTS/MATERIALS, SETTING, METHODS: Participants in this study are females who provided blood and/or endometrial tissue samples in a hospital setting. The key SNPs were genotyped using Sequenom MassARRAY. The functional roles and regulatory annotations for identified SNPs are predicted by various publicly available bioinformatics tools. Endometrial GREB1 expression work employed qRT-PCR, western blotting and immunohistochemistry studies. MAIN RESULTS AND THE ROLE OF CHANCE: Fine mapping results identified a number of SNPs showing stronger association (0.004 < P < 0.032) with endometriosis risk than the original GWAS SNP (rs13394619) (P = 0.034). Some of these SNPs were predicted to have functional roles, for example, interaction with transcription factor motifs. The haplotype (a combination of alleles) formed by the risk alleles from two common SNPs showed significant association (P = 0.026) with endometriosis and epistasis analysis showed no evidence for interaction between the two SNPs, suggesting an additive effect of SNPs on endometriosis risk. In normal human endometrium, GREB1 protein expression was altered depending on the cycle stage (significantly different in late proliferative versus late secretory, P < 0.05) and cell type (glandular epithelium, not stromal cells). However, GREB1 expression in endometriosis cases versus controls and eQTL analyses did not reveal any significant changes. LIMITATIONS, REASONS FOR CAUTION: In silico prediction tools are generally based on cell lines different to our tissue and disease of interest. Functional annotations drawn from these analyses should be considered with this limitation in mind. We identified cell-specific and hormone-specific changes in GREB1 protein expression. The lack of a significant difference observed following our GREB1 expression studies may be the result of moderate power on mixed cell populations in the endometrial tissue samples. WIDER IMPLICATIONS OF THE FINDINGS: This study further implicates the GREB1 region on chromosome 2p25.1 and the GREB1 gene with involvement in endometriosis risk. More detailed functional studies are required to determine the role of the novel GREB1 transcripts in endometriosis pathophysiology. STUDY FUNDING/COMPETING INTERESTS: Funding for this work was provided by NHMRC Project Grants APP1012245, APP1026033, APP1049472 and APP1046880. There are no competing interests.

Inference of the genetic architecture underlying BMI and height with the use of 20,240 sibling pairs

Evidence that complex traits are highly polygenic has been presented by population-based genome-wide association studies (GWASs) through the identification of many significant variants, as well as by family-based de novo sequencing studies indicating that several traits have a large mutational target size. Here, using a third study design, we show results consistent with extreme polygenicity for body mass index (BMI) and height. On a sample of 20,240 siblings (from 9,570 nuclear families), we used a within-family method to obtain narrow-sense heritability estimates of 0.42 (SE = 0.17, p = 0.01) and 0.69 (SE = 0.14, p = 6 x 10(-)(7)) for BMI and height, respectively, after adjusting for covariates. The genomic inflation factors from locus-specific linkage analysis were 1.69 (SE = 0.21, p = 0.04) for BMI and 2.18 (SE = 0.21, p = 2 x 10(-10)) for height. This inflation is free of confounding and congruent with polygenicity, consistent with observations of ever-increasing genomic-inflation factors from GWASs with large sample sizes, implying that those signals are due to true genetic signals across the genome rather than population stratification. We also demonstrate that the distribution of the observed test statistics is consistent with both rare and common variants underlying a polygenic architecture and that previous reports of linkage signals in complex traits are probably a consequence of polygenic architecture rather than the segregation of variants with large effects. The convergent empirical evidence from GWASs, de novo studies, and within-family segregation implies that family-based sequencing studies for complex traits require very large sample sizes because the effects of causal variants are small on average.