Several aspects of sleep behaviour such as timing, duration and quality have been demonstrated to be heritable. To identify common variants that influence sleep traits in the population, we conducted a genome-wide association study of 6 sleep phenotypes assessed by questionnaire in a sample of 2,323 individuals from the Australian Twin Registry. Genotyping was performed on the Illumina 317K, 370K and 610K arrays and the common Single Nucleotide Polymorphisms between platforms were used to impute non-genotyped SNPs. We tested for association with more than 2,000,000 common polymorphisms across the genome. While no SNPs reached the genome-wide significance threshold, we identified a number of associations in plausible candidate genes. Most notably, a group of SNPs in the 3rd intron of the CACNA1C gene ranked as most significant in the analysis of sleep latency (p = 1.3 × 10−6). We attempted to replicate this association in an independent sample from the Chronogen Consortium (n = 2,034), but found no evidence of association (p = 0.73). We have identified several other associations that await replication in an independent sample. Our study had good power to detect common single nucleotide polymorphisms that explain more than 2% of the phenotypic variance in self-report sleep phenotypes at a genome-wide significant level. No such variants were detected.
insomnia; genetics; mood; sleep; circadian
The pubertal height growth spurt is a distinctive feature of childhood growth reflecting both the central onset of puberty and local growth factors. Although little is known about the underlying genetics, growth variability during puberty correlates with adult risks for hormone-dependent cancer and adverse cardiometabolic health. The only gene so far associated with pubertal height growth, LIN28B, pleiotropically influences childhood growth, puberty and cancer progression, pointing to shared underlying mechanisms. To discover genetic loci influencing pubertal height and growth and to place them in context of overall growth and maturation, we performed genome-wide association meta-analyses in 18 737 European samples utilizing longitudinally collected height measurements. We found significant associations (P < 1.67 × 10−8) at 10 loci, including LIN28B. Five loci associated with pubertal timing, all impacting multiple aspects of growth. In particular, a novel variant correlated with expression of MAPK3, and associated both with increased prepubertal growth and earlier menarche. Another variant near ADCY3-POMC associated with increased body mass index, reduced pubertal growth and earlier puberty. Whereas epidemiological correlations suggest that early puberty marks a pathway from rapid prepubertal growth to reduced final height and adult obesity, our study shows that individual loci associating with pubertal growth have variable longitudinal growth patterns that may differ from epidemiological observations. Overall, this study uncovers part of the complex genetic architecture linking pubertal height growth, the timing of puberty and childhood obesity and provides new information to pinpoint processes linking these traits.
Partner and localizer of BRCA2 (PALB2) interacts with BRCA2 to enable double strand break repair through homologous recombination. Similar to BRCA2, germline mutations in PALB2 have been shown to predispose to Fanconi anaemia as well as pancreatic and breast cancer. The PALB2/BRCA2 protein interaction, as well as the increased melanoma risk observed in families harbouring BRCA2 mutations, makes PALB2 a candidate for melanoma susceptibility. In order to assess PALB2 as a melanoma predisposition gene, we sequenced the entire protein-coding sequence of PALB2 in probands from 182 melanoma families lacking pathogenic mutations in known high penetrance melanoma susceptibility genes: CDKN2A, CDK4, and BAP1. In addition, we interrogated whole-genome and exome data from another 19 kindreds with a strong family history of melanoma for deleterious mutations in PALB2. Here we report a rare known deleterious PALB2 mutation (rs118203998) causing a premature truncation of the protein (p.Y1183X) in an individual who had developed four different cancer types, including melanoma. Three other family members affected with melanoma did not carry the variant. Overall our data do not support a case for PALB2 being associated with melanoma predisposition.
Delta opioid receptors are implicated in a variety of psychiatric and neurological disorders. These receptors play a key role in the reinforcing properties of drugs of abuse, and polymorphisms in OPRD1 (the gene encoding delta opioid receptors) are associated with drug addiction. Delta opioid receptors are also involved in protecting neurons against hypoxic and ischemic stress. Here, we first examined a large sample of 738 elderly participants with neuroimaging and genetic data from the Alzheimer’s Disease Neuroimaging Initiative. We hypothesized that common variants in OPRD1 would be associated with differences in brain structure, particularly in regions relevant to addictive and neurodegenerative disorders. One very common variant (rs678849) predicted differences in regional brain volumes. We replicated the association of this single-nucleotide polymorphism with regional tissue volumes in a large sample of young participants in the Queensland Twin Imaging study. Although the same allele was associated with reduced volumes in both cohorts, the brain regions affected differed between the two samples. In healthy elderly, exploratory analyses suggested that the genotype associated with reduced brain volumes in both cohorts may also predict cerebrospinal fluid levels of neurodegenerative biomarkers, but this requires confirmation. If opiate receptor genetic variants are related to individual differences in brain structure, genotyping of these variants may be helpful when designing clinical trials targeting delta opioid receptors to treat neurological disorders.
neuroimaging; genetics; neurodegeneration; drug addiction; opiates
Elevated resting heart rate is associated with greater risk of cardiovascular disease and mortality. In a 2-stage meta-analysis of genome-wide association studies in up to 181,171 individuals, we identified 14 new loci associated with heart rate and confirmed associations with all 7 previously established loci. Experimental downregulation of gene expression in Drosophila melanogaster and Danio rerio identified 20 genes at 11 loci that are relevant for heart rate regulation and highlight a role for genes involved in signal transmission, embryonic cardiac development and the pathophysiology of dilated cardiomyopathy, congenital heart failure and/or sudden cardiac death. In addition, genetic susceptibility to increased heart rate is associated with altered cardiac conduction and reduced risk of sick sinus syndrome, and both heart rate–increasing and heart rate–decreasing variants associate with risk of atrial fibrillation. Our findings provide fresh insights into the mechanisms regulating heart rate and identify new therapeutic targets.
Deficits in lentiform nucleus volume and morphometry are implicated in a number of genetically influenced disorders, including Parkinson’s disease, schizophrenia, and ADHD. Here we performed genome-wide searches to discover common genetic variants associated with differences in lentiform nucleus volume in human populations. We assessed structural MRI scans of the brain in two large genotyped samples: the Alzheimer’s Disease Neuroimaging Initiative (ADNI; N=706) and the Queensland Twin Imaging Study (QTIM; N=639). Statistics of association from each cohort were combined meta-analytically using a fixed-effects model to boost power and to reduce the prevalence of false positive findings. We identified a number of associations in and around the flavin-containing monooxygenase (FMO) gene cluster. The most highly associated SNP, rs1795240, was located in the FMO3 gene; after meta-analysis, it showed genome-wide significant evidence of association with lentiform nucleus volume (PMA=4.79×10−8). This commonly-carried genetic variant accounted for 2.68 % and 0.84 % of the trait variability in the ADNI and QTIM samples, respectively, even though the QTIM sample was on average 50 years younger. Pathway enrichment analysis revealed significant contributions of this gene to the cytochrome P450 pathway, which is involved in metabolizing numerous therapeutic drugs for pain, seizures, mania, depression, anxiety, and psychosis. The genetic variants we identified provide replicated, genome-wide significant evidence for the FMO gene cluster’s involvement in lentiform nucleus volume differences in human populations.
Basal ganglia; Genome-wide association study (GWAS); MRI; Replication; Morphometry; Drug metabolism
Identification of single nucleotide polymorphisms (SNPs) associated with development of advanced colorectal adenomas.
Discovery Phase: 1,406 Caucasian patients (139 advanced adenoma cases and 1,267 controls) from the Adenoma Prevention with Celecoxib (APC) trial were included in a genome-wide association study (GWAS) to identify variants associated with post-polypectomy disease recurrence. Genome-wide significance was defined as false discovery rate < 0.05, unadjusted p=7.4×10−7. Validation Phase: Results were further evaluated using 4,175 familial colorectal adenoma or CRC cases and 5,036 controls from patients of European ancestry (COloRectal Gene Identification consortium, Scotland, Australia and VQ58).
Our study identified eight SNPs associated with advanced adenoma risk in the APC trial (rs2837156, rs7278863, rs2837237, rs2837241, rs2837254, rs741864 at 21q22.2, and rs1381392 and rs17651822 at 3p24.1, at p<10–7 level with odds ratio – OR>2). Five variants in strong pairwise linkage disequilbrium (rs7278863, rs2837237, rs741864, rs741864 and rs2837241, r2=0.8–1) are in or near the coding region for the tight junction adhesion protein, IGSF5. An additional variant associated with advanced adenomas, rs1535989 (minor allele frequency 0.11; OR 2.09; 95% confidence interval 1.50–2.91), also predicted CRC development in a validation analysis (p=0.019) using a series of adenoma cases or CRC (CORGI study) and 3 sets of CRC cases and controls (Scotland, VQ58 and Australia, N=9,211).
Our results suggest that common polymorphisms contribute to the risk of developing advanced adenomas and might also contribute to the risk of developing CRC. The variant at rs1535989 may identify patients whose risk for neoplasia warrants increased colonoscopic surveillance.
Colorectal adenomas; colorectal cancer screening; genetic predisposition
Despite the important role DNA methylation plays in transcriptional regulation, the transgenerational inheritance of DNA methylation is not well understood. The genetic heritability of DNA methylation has been estimated using twin pairs, although concern has been expressed whether the underlying assumption of equal common environmental effects are applicable due to intrauterine differences between monozygotic and dizygotic twins. We estimate the heritability of DNA methylation on peripheral blood leukocytes using Illumina HumanMethylation450 array using a family based sample of 614 people from 117 families, allowing comparison both within and across generations.
The correlations from the various available relative pairs indicate that on average the similarity in DNA methylation between relatives is predominantly due to genetic effects with any common environmental or zygotic effects being limited. The average heritability of DNA methylation measured at probes with no known SNPs is estimated as 0.187. The ten most heritable methylation probes were investigated with a genome-wide association study, all showing highly statistically significant cis mQTLs. Further investigation of one of these cis mQTL, found in the MHC region of chromosome 6, showed the most significantly associated SNP was also associated with over 200 other DNA methylation probes in this region and the gene expression level of 9 genes.
The majority of transgenerational similarity in DNA methylation is attributable to genetic effects, and approximately 20% of individual differences in DNA methylation in the population are caused by DNA sequence variation that is not located within CpG sites.
Imaging genetics aims to discover how variants in the human genome influence brain measures derived from images. Genome-wide association scans (GWAS) can screen the genome for common differences in our DNA that relate to brain measures. In small samples, GWAS has low power as individual gene effects are weak and one must also correct for multiple comparisons across the genome and the image. Here we extend recent work on genetic clustering of images, to analyze surface-based models of anatomy using GWAS. We performed spherical harmonic analysis of hippocampal surfaces, automatically extracted from brain MRI scans of 1254 subjects. We clustered hippocampal surface regions with common genetic influences by examining genetic correlations (rg) between the normalized deformation values at all pairs of surface points. Using genetic correlations to cluster surface measures, we were able to boost effect sizes for genetic associations, compared to clustering with traditional phenotypic correlations using Pearson's r.
heritability; GWAS; clustering; hippocampus; 3D surfaces; imaging genetics
To mine possibly hidden causal single nucleotide polymorphisms (SNPs) in the etiology of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome-wide association study (GWAS) dataset with 1804 melanoma cases and 1,026 cancer-free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in other two GWAS datasets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL dataset (P = 0.013 and 0.004, respectively), but failed validation in the Australia dataset. Genotype-phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression levels of PSMB9. Further experiments revealed that the promoter SNP rs2071480, which is in high LD with rs1351383 and rs2127675, involved in influencing transcription factor binding and gene expression. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.
GWAS; Cell cycle; PSMB9; Polymorphism; melanoma
Endometriosis, defined as estrogen-dependent lesions containing endometrial glands and stroma outside the uterus, is a chronic and often painful gynecological condition that affects 6% to 10% of reproductive age women. Endometriosis has estimated annual costs of US $12 419 per woman (approximately €9579), comprising one-third of the direct health care costs with two-thirds attributed to loss of productivity. Decreased quality of life is the most important predictor of direct health care and total costs. It has been estimated that there is a mean delay of 6.7 years between onset of symptoms and a surgical diagnosis of endometriosis, and each affected woman loses on average 10.8 hours of work weekly, mainly owing to reduced effectiveness while working. To encourage and facilitate research into this debilitating disease, a consensus workshop to define future directions for endometriosis research was held as part of the 11th World Congress on Endometriosis in September 2011 in Montpellier, France. The objective of this workshop was to review and update the endometriosis research priorities consensus statement developed following the 10th World Congress on Endometriosis in 2008.1 A total of 56 recommendations for research have been developed, grouped under 6 subheadings: (1) diagnosis, (2) classification and prognosis, (3) clinical trials, treatment, and outcomes, (4) epidemiology, (5) pathophysiology, and (6) research policy. By producing this consensus international research priorities statement, it is the hope of the workshop participants that researchers will be encouraged to develop new interdisciplinary research proposals that will attract increased funding support for work on endometriosis.
endometriosis; research directions; international workshop; consensus report
Disordered gambling is a moderately heritable trait, but the underlying genetic basis is largely unknown. We performed a genome-wide association study (GWAS) for disordered gambling using a quantitative factor score in 1,312 twins from 894 Australian families. Association was conducted for 2,381,914 single nucleotide polymorphisms (SNPs) using the family-based association test in Merlin followed by gene and pathway enrichment analyses. Although no SNP reached genome-wide significance, six achieved P-values < 1 × 10−5 with variants in three genes (MT1X, ATXN1 and VLDLR) implicated in disordered gambling. Secondary case-control analyses found two SNPs on chromosome 9 (rs1106076 and rs12305135 near VLDLR) and rs10812227 near FZD10 on chromosome 12 to be significantly associated with lifetime DSM-IV pathological gambling and SOGS classified probable pathological gambling status. Furthermore, several addiction-related pathways were enriched for SNPs associated with disordered gambling. Finally, gene-based analysis of 24 candidate genes for dopamine agonist induced gambling in individuals with Parkinson’s disease suggested an enrichment of SNPs associated with disordered gambling. We report the first GWAS of disordered gambling. While further replication is required, the identification of susceptibility loci and biological pathways will be important in characterizing the biological mechanisms that underpin disordered gambling.
association; disordered gambling; genomewide; MERLIN; quantitative
Several studies have identified genes associated with alcohol use disorders, but the variation in each of these genes explains only a small portion of the genetic vulnerability. The goal of the present study was to perform a genome-wide association study (GWAS) in extended families from the Collaborative Study on the Genetics of Alcoholism (COGA) to identify novel genes affecting risk for alcohol dependence. To maximize the power of the extended family design we used a quantitative endophenotype, measured in all individuals: number of alcohol dependence symptoms endorsed (symptom count). Secondary analyses were performed to determine if the single nucleotide polymorphisms (SNPs) associated with symptom count were also associated with the dichotomous phenotype, DSM-IV alcohol dependence. This family-based GWAS identified SNPs in C15orf53 that are strongly associated with DSM-IV alcohol (p=4.5×10−8, inflation corrected p=9.4×10−7). Results with DSM-IV alcohol dependence in the regions of interest support our findings with symptom count, though the associations were less significant. Attempted replications of the most promising association results were conducted in two independent samples: non-overlapping subjects from the Study of Addiction: Genes and Environment (SAGE) and the Australian twin-family study of alcohol use disorders (OZALC). Nominal association of C15orf53 with symptom count was observed in SAGE. The variant that showed strongest association with symptom count, rs12912251 and its highly correlated variants (D′=1, r2≥ 0.95), has previously been associated with risk for bipolar disorder.
DSM-IV alcohol dependence symptoms; Family-based GWAS; C15orf53; Quantitative traits
Inter-individual variation in mean leukocyte telomere length (LTL) is associated with cancer and several age-associated diseases. Here, in a genome-wide meta-analysis of 37,684 individuals with replication of selected variants in a further 10,739 individuals, we identified seven loci, including five novel loci, associated with mean LTL (P<5x10−8). Five of the loci contain genes (TERC, TERT, NAF1, OBFC1, RTEL1) that are known to be involved in telomere biology. Lead SNPs at two loci (TERC and TERT) associate with several cancers and other diseases, including idiopathic pulmonary fibrosis. Moreover, a genetic risk score analysis combining lead variants at all seven loci in 22,233 coronary artery disease cases and 64,762 controls showed an association of the alleles associated with shorter LTL with increased risk of CAD (21% (95% CI: 5–35%) per standard deviation in LTL, p=0.014). Our findings support a causal role of telomere length variation in some age-related diseases.
Identifying the downstream effects of disease-associated single nucleotide polymorphisms (SNPs) is challenging: the causal gene is often unknown or it is unclear how the SNP affects the causal gene, making it difficult to design experiments that reveal functional consequences. To help overcome this problem, we performed the largest expression quantitative trait locus (eQTL) meta-analysis so far reported in non-transformed peripheral blood samples of 5,311 individuals, with replication in 2,775 individuals. We identified and replicated trans-eQTLs for 233 SNPs (reflecting 103 independent loci) that were previously associated with complex traits at genome-wide significance. Although we did not study specific patient cohorts, we identified trait-associated SNPs that affect multiple trans-genes that are known to be markedly altered in patients: for example, systemic lupus erythematosus (SLE) SNP rs49170141 altered C1QB and five type 1 interferon response genes, both hallmarks of SLE2-4. Subsequent ChIP-seq data analysis on these trans-genes implicated transcription factor IKZF1 as the causal gene at this locus, with DeepSAGE RNA-sequencing revealing that rs4917014 strongly alters 3’ UTR levels of IKZF1. Variants associated with cholesterol metabolism and type 1 diabetes showed similar phenomena, indicating that large-scale eQTL mapping provides insight into the downstream effects of many trait-associated variants.
Approaches exploiting extremes of the trait distribution may reveal novel loci for common traits, but it is unknown whether such loci are generalizable to the general population. In a genome-wide search for loci associated with upper vs. lower 5th percentiles of body mass index, height and waist-hip ratio, as well as clinical classes of obesity including up to 263,407 European individuals, we identified four new loci (IGFBP4, H6PD, RSRC1, PPP2R2A) influencing height detected in the tails and seven new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3, ZZZ3) for clinical classes of obesity. Further, we show that there is large overlap in terms of genetic structure and distribution of variants between traits based on extremes and the general population and little etiologic heterogeneity between obesity subgroups.
Alcohol consumption has multiple biochemical consequences. Only a few of these are useful as diagnostic markers but many reflect potentially harmful or beneficial effects of alcohol. Average consumption of two to four drinks per day is associated with lower overall or cardiovascular mortality risk than either lower or higher intake. We have analysed the dose-response relationships between reported alcohol consumption and 17 biomarkers, with emphasis on intake of up to three drinks per day.
Biochemical tests were performed on serum from 8396 study participants (3750 men and 4646 women, aged 51 ± 13 years, range 18–93) who had provided information on alcohol consumption in the week preceding blood collection.
GGT, ALT, AST, CDT, urate, ferritin and bilirubin showed little or no change with alcohol consumption below two to three drinks per day, but increased with higher intake. HDL-C and albumin showed increasing results, and insulin showed decreasing results, across the entire range of alcohol use. Biphasic responses, where subjects reporting one to two drinks per day had lower results than those reporting either more or less alcohol use, occurred for triglycerides, glucose, C-reactive protein, alkaline phosphatase and butyrylcholinesterase. Increasing alcohol use was associated with decreasing LDL-C in younger women, but higher LDL-C in older men.
Some markers show threshold relationships with alcohol, others show continuous ones, and a third group show biphasic or U-shaped relationships. Overall the biochemical sequelae of low-to-moderate alcohol use are consistent with the epidemiological evidence on morbidity and mortality.
Alcohol; Biomarkers; Dose-response curve; Population study
Genome-wide association studies (GWASs) have mainly focused on top significant single nucleotide polymorphisms (SNPs), most of which did not have clear biological functions but were just surrogates for unknown causal variants. Studying SNPs with modest association and putative functions in biologically plausible pathways has become one complementary approach to GWASs. To unravel the key roles of mitogen-activated protein kinase (MAPK) pathways in cutaneous melanoma (CM) risk, we re-evaluated the associations between 47 818 SNPs in 280 MAPK genes and CM risk using our published GWAS dataset with 1804 CM cases and 1026 controls. We initially found 105 SNPs with P ≤ 0.001, more than expected by chance, 26 of which were predicted to be putatively functional SNPs. The risk associations with 16 SNPs around DUSP14 (rs1051849) and a previous reported melanoma locus MAFF/PLA2G6 (proxy SNP rs4608623) were replicated in the GenoMEL dataset (P < 0.01) but failed in the Australian dataset. Meta-analysis showed that rs1051849 in the 3ʹ untranslated regions of DUSP14 was associated with a reduced risk of melanoma (odds ratio = 0.89, 95% confidence interval: 0.82–0.96, P = 0.003, false discovery rate = 0.056). Further genotype–phenotype correlation analysis using the 90 HapMap lymphoblastoid cell lines from Caucasians showed significant correlations between two SNPs (rs1051849 and rs4608623) and messenger RNA expression levels of DUSP14 and MAFF (P = 0.025 and P = 0.010, respectively). Gene-based tests also revealed significant SNPs were over-represented in MAFF, PLA2G6, DUSP14 and other 16 genes. Our results suggest that functional SNPs in MAPK pathways may contribute to CM risk. Further studies are warranted to validate our findings.
Estrogen receptor (ER)-negative tumors represent 20–30% of all breast cancers, with a higher proportion occurring in younger women and women of African ancestry1. The etiology2 and clinical behavior3 of ER-negative tumors are different from those of tumors expressing ER (ER positive), including differences in genetic predisposition4. To identify susceptibility loci specific to ER-negative disease, we combined in a meta-analysis 3 genome-wide association studies of 4,193 ER-negative breast cancer cases and 35,194 controls with a series of 40 follow-up studies (6,514 cases and 41,455 controls), genotyped using a custom Illumina array, iCOGS, developed by the Collaborative Oncological Gene-environment Study (COGS). SNPs at four loci, 1q32.1 (MDM4, P = 2.1 × 10−12 and LGR6, P = 1.4 × 10−8), 2p24.1 (P = 4.6 × 10−8) and 16q12.2 (FTO, P = 4.0 × 10−8), were associated with ER-negative but not ER-positive breast cancer (P > 0.05). These findings provide further evidence for distinct etiological pathways associated with invasive ER-positive and ER-negative breast cancers.
The BRAF (V600E) mutation in colorectal cancers that are microsatellite stable (MSS) confers a poor patient prognosis, whereas BRAF mutant microsatellite-unstable (MSI) colorectal cancers have an excellent prognosis. BRAF wild type cancers are typically MSS and display chromosomal instability (CIN). CIN has not been extensively studied on a genome-wide basis in relation to BRAF mutational status in colorectal cancer. BRAF mutant/MSS (BRAFmut/MSS) cancers (n = 33) and BRAF mutant/MSI (BRAFmut/MSI) cancers (n = 30) were compared for presence of copy number aberrations (CNAs) indicative of CIN, with BRAF wild type/MSS (BRAFwt/MSS) cancers (n = 18) using Illumina CytoSNP-12 arrays. BRAFmut/MSS and BRAFwt/MSS cancers showed comparable numbers of CNAs/cancer at 32.8 and 29.8 respectively. However, there were differences in patterns of CNA length between MSS cohorts, with BRAFmut/MSS cancers having significantly greater proportions of focal CNAs compared to BRAFwt/MSS cancers (p<0.0001); whereas whole chromosomal arm CNAs were more common in BRAFwt/MSS cancers (p<0.0001). This related to a reduced average CNA length in BRAFmut/MSS compared to BRAFwt/MSS cancers (20.7 Mb vs 33.4 Mb;p<0.0001); and a smaller average percent of CIN affected genomes in BRAFmut/MSS compared to BRAFwt/MSS cancers (23.9% vs 34.9% respectively). BRAFmut/MSI cancers were confirmed to have low CNA rates (5.4/cancer) and minimal CIN-affected genomes (average of 4.5%) compared to MSS cohorts (p<0.0001). BRAFmut/MSS cancers had more frequent deletion CNAs compared to BRAFwt/MSS cancers on 6p and 17q at loci not typically correlated with colorectal cancer, and greater amplification CNAs on 8q and 18q compared to BRAFwt/MSS cancers. These results indicate that comparable rates of CIN occur between MSS subgroups, however significant differences in their patterns of instability exist, with BRAFmut/MSS cancers showing a ‘focal pattern’ and BRAFwt/MSS cancers having a ‘whole arm pattern’ of CIN. This and the genomic loci more frequently affected in BRAFmut/MSS cancers provides further evidence of the biological distinctions of this important cancer subgroup.
Neuronal nicotinic acetylcholine receptor (nAChR) genes (CHRNA5/CHRNA3/CHRNB4) have been reproducibly associated with nicotine dependence, smoking behaviors, and lung cancer risk. Of the few reports that have focused on early smoking behaviors, association results have been mixed. This meta-analysis examines early smoking phenotypes and SNPs in the gene cluster to determine: (1) whether the most robust association signal in this region (rs16969968) for other smoking behaviors is also associated with early behaviors, and/or (2) if additional statistically independent signals are important in early smoking. We focused on two phenotypes: age of tobacco initiation (AOI) and age of first regular tobacco use (AOS). This study included 56,034 subjects (41 groups) spanning nine countries and evaluated five SNPs including rs1948, rs16969968, rs578776, rs588765, and rs684513. Each dataset was analyzed using a centrally generated script. Meta-analyses were conducted from summary statistics. AOS yielded significant associations with SNPs rs578776 (beta = 0.02, P = 0.004), rs1948 (beta = 0.023, P = 0.018), and rs684513 (beta = 0.032, P = 0.017), indicating protective effects. There were no significant associations for the AOI phenotype. Importantly, rs16969968, the most replicated signal in this region for nicotine dependence, cigarettes per day, and cotinine levels, was not associated with AOI (P = 0.59) or AOS (P = 0.92). These results provide important insight into the complexity of smoking behavior phenotypes, and suggest that association signals in the CHRNA5/A3/B4 gene cluster affecting early smoking behaviors may be different from those affecting the mature nicotine dependence phenotype.
CHRNA5; CHRNA3; CHRNB4; meta-analysis; nicotine; smoke
The genetic contribution to liability for opioid dependence is well-established; identification of the responsible genes has proved challenging.
To examine association of 1430 candidate gene single-nucleotide polymorphisms (SNPs) with heroin dependence, reporting here only the 71 SNPs in the chromosome 11 gene cluster (NCAM1, TTC12, ANKK1, DRD2) that include the strongest observed associations.
Case-control genetic association study that included two control groups (lacking an established optimal control group).
Semi-structured psychiatric interviews
Australian cases (N=1459) ascertained from opioid replacement therapy (ORT) clinics, neighborhood controls (N=531) ascertained from economically disadvantaged areas near opioid replacement therapy clinics, and unrelated Australian Twin Registry (ATR) controls (N=1495) not dependent on alcohol or illicit drugs selected from a twin and family sample.
Main Outcome Measure
Lifetime heroin dependence
Comparison of cases with Australian Twin Registry controls found minimal evidence of association for all chromosome 11 cluster SNPs (p≥.01); a similar comparison to neighborhood controls revealed greater differences (p≥1.8 × 10−4). Comparing cases (N=1459) with the subgroup of neighborhood controls not dependent on illicit drugs (N=340), three SNPs were significantly associated (correcting for multiple testing): ANKK1 SNP rs877138 [most strongly associated; odds ratio 1.59; 95%CI (1.32–1.92); p=9.7 × 10−7], ANKK1 SNP rs4938013 and TTC12 SNP rs7130431. A similar pattern of association was observed when comparing illicit drug-dependent (N=191) and non-dependent (N=340) neighborhood controls, suggesting that liability likely extends to non-opioid illicit drug dependence. Aggregate heroin dependence risk associated with two SNPs, rs877138 and rs4492854 (located in NCAM1), varied more than 4-fold (p= 2.74 × 10−9 for the risk-associated linear trend).
Our results provide further evidence of association for chromosome 11 gene cluster SNPs with substance dependence, including extension of liability to illicit drug dependence. Our findings highlight the necessity of considering drug exposure history when selecting control groups for genetic investigations of illicit drug dependence.
Borderline personality disorder (BPD) and substance use disorders frequently cooccur; their dual presence predicts poor prognosis. The genetic underpinnings of BPD have not been well-characterized and could offer insight into comorbidity. The current report focuses on the association of Neurexin 3 (NRXN3) single nucleotide polymorphisms (SNPs) with BPD symptoms in heroin dependent cases and controls.
The sample of the Comorbidity and Trauma Study, a genetic association study of heroin dependence, consists of Australian heroin dependent cases ascertained from opioid replacement therapy clinics and controls ascertained in nearby economically-disadvantaged neighborhoods. The assessment included a screening instrument for BPD, used previously in Australian population surveys. Genotypic and BPD phenotypic data were available for 1439 cases and 507 controls. We examined the association of 1430 candidate gene SNPs with BPD phenotypes.
One or more NRXN3 SNPs were nominally associated with all BPD phenotypes; however, none met the conservative significance threshold we employed to correct for multiple testing. The most strongly associated SNPs included rs10144398 with identity disturbance (p=4.9 × 10−5) and rs10151731 with affective instability (p=8.8 × 10−5). The strongest association with screening positive for BPD was found for the NRXN3 SNP, rs10083466 (p=.0013). Neither the correlation of BPD phenotypes nor the linkage disequilibrium relationships of the SNPs account for the number of observed associations involving NRXN3 SNPs.
Our findings provide intriguing preliminary evidence for the association of NRXN3 with BPD phenotypes. The strongest associations were found for traits (i.e., affective instability; identity disturbance) also observed with other disorders.
borderline personality disorder; NRXN3; genetic association study; heroin dependence
People meeting diagnostic criteria for anxiety or depressive disorders tend to score high on the personality scale of neuroticism. Studying this personality dimension can give insights into the aetiology of these important psychiatric disorders.
To undertake a comprehensive genome-wide linkage study of neuroticism, using large study samples that have been measured multiple times. To compare the results between countries for replication and across time within countries for consistency.
Genome wide linkage scan.
Twin individuals and their family members from Australia (AU) and the Netherlands (NL).
19,635 sibling pairs completed self-report questionnaires for neuroticism up to five times over a period of up to 22 years. 5,069 sibling pairs were genotyped with microsatellite markers.
Non-parametric linkage analyses were conducted in Merlin-Regress for the mean neuroticism scores averaged across time. Additional analyses were conducted for the time specific measures of neuroticism from each country to investigate consistency of linkage results.
Three chromosomal regions exceeded empirically-derived thresholds for suggestive linkage using mean neuroticism scores: 10p 5 cM (NL), 14q 103 cM (NL) and 18q 117 cM (AU & NL combined), but only 14q retains significance after correction for multiple testing. These regions all showed evidence for linkage in individual time-specific measures of neuroticism and one (18q) showed some evidence for replication between countries. Linkage intervals for these regions all overlap with regions identified in other studies of neuroticism or related traits and/or in studies of anxiety in mice.
Our results demonstrate the value of the availability of multiple measures over time and add to the optimism reported in recent reviews for replication of linkage regions for neuroticism. These regions are likely to harbour causal variants for neuroticism and its related psychiatric disorders and can inform prioritisation of results from genome-wide association studies.