Genetic risk factors for endometriosis near estrogen receptor 1 and coexpression of genes in this region in endometrium.

The etiology and pathogenesis of endometriosis are complex with both genetic and environmental factors contributing to disease risk. Genome-wide association studies (GWAS) have identified multiple signals in the estrogen receptor 1 (ESR1) region associated with endometriosis and other reproductive traits and diseases. In addition, candidate gene association studies identified signals in the ESR1 region associated with endometriosis risk suggesting genetic regulation of genes in this region may be important for reproductive health. This study aimed to investigate hormonal and genetic regulation of genes in the ESR1 region in human endometrium. Changes in serum oestradiol and progesterone concentrations and expression of hormone receptors ESR1 and progesterone receptor (PGR) were assessed in endometrial samples from 135 women collected at various stages of the menstrual cycle. Correlation between hormone concentrations, receptor expression and expression of genes in the ESR1 locus was investigated. The effect of endometriosis risk variants on expression of genes in the region was analyzed to identify gene targets. Hormone concentrations and receptor expression varied significantly across the menstrual cycle. Expression of genes in the ESR1 region correlated with progesterone concentration; however, they were more strongly correlated with expression of ESR1 and PGR suggesting coregulation of genes. There was no evidence that endometriosis risk variants directly regulated expression of genes in the region. Limited sample size and cellular heterogeneity in endometrial tissue may impact the ability to detect significant genetic effects on gene expression. Effects of these variants should be validated in a larger dataset and in relevant individual cell types.

Prevalence and incidence of endometriosis in Australian women: a data linkage cohort study.

OBJECTIVE: To estimate the prevalence and incidence of endometriosis among Australian women. DESIGN: Population-based cohort study linked to administrative health records. SETTING: Secondary analysis of seven surveys collected between 2000 and 2018 from a population-based cohort study. POPULATION: A total of 13 508 Australian women, born 1973-78, from a prospective cohort study of 14 247 women conducted between 1996 and 2018. METHODS: During 2000 and 2018, self-reported longitudinal survey data were linked to three administrative health databases to separately identify women with clinically confirmed or suspected endometriosis across the multiple data sources. MAIN OUTCOME MEASURES: Prevalence and incidence of clinically confirmed endometriosis in the cohort were first estimated using national hospital data. Data were then combined with other administrative health databases and the survey data to capture all clinically confirmed and suspected diagnoses of endometriosis. RESULTS: The cumulative prevalence of clinically confirmed endometriosis was 6.0% (95% CI 5.8-6.2%) by age 40-44 years. The cumulative prevalence increased to 11.4% (95% CI 11.1-11.7%) when adding diagnoses of clinically suspected endometriosis. Age-specific incidence estimates peaked to 6 per 1000 person-years at age 30-34 years. CONCLUSIONS: Among 13 508 Australian women followed for 20 years, one in nine women had clinically confirmed or suspected endometriosis by the age of 44, with most diagnosed during their early thirties. Endometriosis is a significant public health issue requiring increased surveillance, clinical awareness and management. Efforts to expand knowledge on the aetiology of the disease and optimal methods for disease management are crucial to women's health. TWEETABLE ABSTRACT: In a national study of 13 508 Australian women, one in nine women were diagnosed with endometriosis by age 44.

Generation of immortalized human endometrial stromal cell lines with different endometriosis risk genotypes.

Endometriotic lesions are composed in part of endometrial-like stromal cells, however, there is a shortage of immortalized human endometrial stromal cultures available for research. As genetic factors play a role in endometriosis risk, it is important that genotype is also incorporated into analysis of pathological mechanisms. Human telomerase reverse transcriptase (hTERT) immortalization (using Lenti-hTERT-green fluorescent protein virus) took place following genotype selection; 13 patients homozygous for either the risk or non-risk 'other' allele for one or more important endometriosis risk single nucleotide polymorphism on chromosome 1p36.12 (rs3820282, rs56318008, rs55938609, rs12037376, rs7521902 or rs12061255). Short tandem repeat DNA profiling validated that donor tissue matched that of the immortalized cell lines and confirmed that cultures were genetically novel. Expression of morphological markers (vimentin and cytokeratin) and key genes of interest (telomerase, estrogen and progesterone receptors and LINC00339) were examined and functional assays for cell proliferation, steroid hormone and inflammatory responses were performed for 7/13 cultures. All endometrial stromal cell lines maintained their fibroblast-like morphology (vimentin-positive) and homozygous endometriosis-risk genotype following introduction of hTERT. Furthermore, the new stromal cultures demonstrated positive and diverse responses to hormones (proliferation and decidualisation changes) and inflammation (dose-dependent response), while maintaining hormone receptor expression. In conclusion, we successfully developed a range of human endometrial stromal cell lines that carry important endometriosis-risk alleles. The wider implications of this approach go beyond advancing endometriosis research; these cell lines will be valuable tools for multiple endometrial pathologies offering a level of genetic and phenotypic diversity not previously available.

A direct test of the diathesis-stress model for depression.

The diathesis-stress theory for depression states that the effects of stress on the depression risk are dependent on the diathesis or vulnerability, implying multiplicative interactive effects on the liability scale. We used polygenic risk scores for major depressive disorder (MDD) calculated from the results of the most recent analysis from the Psychiatric Genomics Consortium as a direct measure of the vulnerability for depression in a sample of 5221 individuals from 3083 families. In the same we also had measures of stressful life events and social support and a depression symptom score, as well as DSM-IV MDD diagnoses for most individuals. In order to estimate the variance in depression explained by the genetic vulnerability, the stressors and their interactions, we fitted linear mixed models controlling for relatedness for the whole sample as well as stratified by sex. We show a significant interaction of the polygenic risk scores with personal life events (0.12% of variance explained, P-value=0.0076) contributing positively to the risk of depression. Additionally, our results suggest possible differences in the aetiology of depression between women and men. In conclusion, our findings point to an extra risk for individuals with combined vulnerability and high number of reported personal life events beyond what would be expected from the additive contributions of these factors to the liability for depression, supporting the multiplicative diathesis-stress model for this disease.

Genome-wide association study identifies a novel locus for cannabis dependence.

Despite moderate heritability, only one study has identified genome-wide significant loci for cannabis-related phenotypes. We conducted meta-analyses of genome-wide association study data on 2080 cannabis-dependent cases and 6435 cannabis-exposed controls of European descent. A cluster of correlated single-nucleotide polymorphisms (SNPs) in a novel region on chromosome 10 was genome-wide significant (lowest P=1.3E-8). Among the SNPs, rs1409568 showed enrichment for H3K4me1 and H3K427ac marks, suggesting its role as an enhancer in addiction-relevant brain regions, such as the dorsolateral prefrontal cortex and the angular and cingulate gyri. This SNP is also predicted to modify binding scores for several transcription factors. We found modest evidence for replication for rs1409568 in an independent cohort of African American (896 cases and 1591 controls; P=0.03) but not European American (EA; 781 cases and 1905 controls) participants. The combined meta-analysis (3757 cases and 9931 controls) indicated trend-level significance for rs1409568 (P=2.85E-7). No genome-wide significant loci emerged for cannabis dependence criterion count (n=8050). There was also evidence that the minor allele of rs1409568 was associated with a 2.1% increase in right hippocampal volume in an independent sample of 430 EA college students (fwe-P=0.008). The identification and characterization of genome-wide significant loci for cannabis dependence is among the first steps toward understanding the biological contributions to the etiology of this psychiatric disorder, which appears to be rising in some developed nations.

The CODATwins Project: The Current Status and Recent Findings of COllaborative Project of Development of Anthropometrical Measures in Twins.

The COllaborative project of Development of Anthropometrical measures in Twins (CODATwins) project is a large international collaborative effort to analyze individual-level phenotype data from twins in multiple cohorts from different environments. The main objective is to study factors that modify genetic and environmental variation of height, body mass index (BMI, kg/m2) and size at birth, and additionally to address other research questions such as long-term consequences of birth size. The project started in 2013 and is open to all twin projects in the world having height and weight measures on twins with information on zygosity. Thus far, 54 twin projects from 24 countries have provided individual-level data. The CODATwins database includes 489,981 twin individuals (228,635 complete twin pairs). Since many twin cohorts have collected longitudinal data, there is a total of 1,049,785 height and weight observations. For many cohorts, we also have information on birth weight and length, own smoking behavior and own or parental education. We found that the heritability estimates of height and BMI systematically changed from infancy to old age. Remarkably, only minor differences in the heritability estimates were found across cultural-geographic regions, measurement time and birth cohort for height and BMI. In addition to genetic epidemiological studies, we looked at associations of height and BMI with education, birth weight and smoking status. Within-family analyses examined differences within same-sex and opposite-sex dizygotic twins in birth size and later development. The CODATwins project demonstrates the feasibility and value of international collaboration to address gene-by-exposure interactions that require large sample sizes and address the effects of different exposures across time, geographical regions and socioeconomic status.

A polymorphism in the OPRM1 3'-untranslated region is associated with methadone efficacy in treating opioid dependence.

The µ-opioid receptor (MOR) is the primary target of methadone and buprenorphine. The primary neuronal transcript of the OPRM1 gene, MOR-1, contains a ~13 kb 3' untranslated region with five common haplotypes in European-Americans. We analyzed the effects of these haplotypes on the percentage of opioid positive urine tests in European-Americans (n=582) during a 24-week, randomized, open-label trial of methadone or buprenorphine/naloxone (Suboxone) for the treatment of opioid dependence. A single haplotype, tagged by rs10485058, was significantly associated with patient urinalysis data in the methadone treatment group. Methadone patients with the A/A genotype at rs10485058 were less likely to have opioid-positive urine drug screens than those in the combined A/G and G/G genotypes group (relative risk=0.76, 95% confidence intervals=0.73-0.80, P=0.0064). Genotype at rs10485058 also predicted self-reported relapse rates in an independent population of Australian patients of European descent (n=1215) who were receiving opioid substitution therapy (P=0.003). In silico analysis predicted that miR-95-3p would interact with the G, but not the A allele of rs10485058. Luciferase assays indicated miR-95-3p decreased reporter activity of constructs containing the G, but not the A allele of rs10485058, suggesting a potential mechanism for the observed pharmacogenetic effect. These findings suggest that selection of a medication for opioid dependence based on rs10485058 genotype might improve outcomes in this ethnic group.

Common polygenic risk for autism spectrum disorder (ASD) is associated with cognitive ability in the general population.

Cognitive impairment is common among individuals diagnosed with autism spectrum disorder (ASD) and attention-deficit hyperactivity disorder (ADHD). It has been suggested that some aspects of intelligence are preserved or even superior in people with ASD compared with controls, but consistent evidence is lacking. Few studies have examined the genetic overlap between cognitive ability and ASD/ADHD. The aim of this study was to examine the polygenic overlap between ASD/ADHD and cognitive ability in individuals from the general population. Polygenic risk for ADHD and ASD was calculated from genome-wide association studies of ASD and ADHD conducted by the Psychiatric Genetics Consortium. Risk scores were created in three independent cohorts: Generation Scotland Scottish Family Health Study (GS:SFHS) (n=9863), the Lothian Birth Cohorts 1936 and 1921 (n=1522), and the Brisbane Adolescent Twin Sample (BATS) (n=921). We report that polygenic risk for ASD is positively correlated with general cognitive ability (beta=0.07, P=6 × 10(-7), r(2)=0.003), logical memory and verbal intelligence in GS:SFHS. This was replicated in BATS as a positive association with full-scale intelligent quotient (IQ) (beta=0.07, P=0.03, r(2)=0.005). We did not find consistent evidence that polygenic risk for ADHD was associated with cognitive function; however, a negative correlation with IQ at age 11 years (beta=-0.08, Z=-3.3, P=0.001) was observed in the Lothian Birth Cohorts. These findings are in individuals from the general population, suggesting that the relationship between genetic risk for ASD and intelligence is partly independent of clinical state. These data suggest that common genetic variation relevant for ASD influences general cognitive ability.

Evidence of CNIH3 involvement in opioid dependence.

Opioid dependence, a severe addictive disorder and major societal problem, has been demonstrated to be moderately heritable. We conducted a genome-wide association study in Comorbidity and Trauma Study data comparing opioid-dependent daily injectors (N=1167) with opioid misusers who never progressed to daily injection (N=161). The strongest associations, observed for CNIH3 single-nucleotide polymorphisms (SNPs), were confirmed in two independent samples, the Yale-Penn genetic studies of opioid, cocaine and alcohol dependence and the Study of Addiction: Genetics and Environment, which both contain non-dependent opioid misusers and opioid-dependent individuals. Meta-analyses found five genome-wide significant CNIH3 SNPs. The A allele of rs10799590, the most highly associated SNP, was robustly protective (P=4.30E-9; odds ratio 0.64 (95% confidence interval 0.55-0.74)). Epigenetic annotation predicts that this SNP is functional in fetal brain. Neuroimaging data from the Duke Neurogenetics Study (N=312) provide evidence of this SNP's in vivo functionality; rs10799590 A allele carriers displayed significantly greater right amygdala habituation to threat-related facial expressions, a phenotype associated with resilience to psychopathology. Computational genetic analyses of physical dependence on morphine across 23 mouse strains yielded significant correlations for haplotypes in CNIH3 and functionally related genes. These convergent findings support CNIH3 involvement in the pathophysiology of opioid dependence, complementing prior studies implicating the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) glutamate system.

Genome-wide autozygosity is associated with lower general cognitive ability.

Inbreeding depression refers to lower fitness among offspring of genetic relatives. This reduced fitness is caused by the inheritance of two identical chromosomal segments (autozygosity) across the genome, which may expose the effects of (partially) recessive deleterious mutations. Even among outbred populations, autozygosity can occur to varying degrees due to cryptic relatedness between parents. Using dense genome-wide single-nucleotide polymorphism (SNP) data, we examined the degree to which autozygosity associated with measured cognitive ability in an unselected sample of 4854 participants of European ancestry. We used runs of homozygosity-multiple homozygous SNPs in a row-to estimate autozygous tracts across the genome. We found that increased levels of autozygosity predicted lower general cognitive ability, and estimate a drop of 0.6 s.d. among the offspring of first cousins (P=0.003-0.02 depending on the model). This effect came predominantly from long and rare autozygous tracts, which theory predicts as more likely to be deleterious than short and common tracts. Association mapping of autozygous tracts did not reveal any specific regions that were predictive beyond chance after correcting for multiple testing genome wide. The observed effect size is consistent with studies of cognitive decline among offspring of known consanguineous relationships. These findings suggest a role for multiple recessive or partially recessive alleles in general cognitive ability, and that alleles decreasing general cognitive ability have been selected against over evolutionary time.

Twenty-eight loci that influence serum urate levels: analysis of association with gout.

OBJECTIVES: Twenty-eight genetic loci are associated with serum urate levels in Europeans. Evidence for association with gout at most loci is absent, equivocal or not replicated. Our aim was to test the loci for association with gout meeting the American College of Rheumatology gout classification criteria in New Zealand European and Polynesian case-control sample sets. METHODS: 648 European cases and 1550 controls, and 888 Polynesian (Ma¯ori and Pacific) cases and 1095 controls were genotyped. Association with gout was tested by logistic regression adjusting for age and sex. Power was adequate (>0.7) to detect effects of OR>1.3. RESULTS: We focused on 24 loci without previous consistent evidence for association with gout. In Europeans, we detected association at seven loci, one of which was the first report of association with gout (IGF1R). In Polynesian, association was detected at three loci. Meta-analysis revealed association at eight loci-two had not previously been associated with gout (PDZK1 and MAF). In participants with higher Polynesian ancestry, there was association in an opposing direction to Europeans at PRKAG2 and HLF (HLF is the first report of association with gout). There was obvious inconsistency of gout association at four loci (GCKR, INHBC, SLC22A11, SLC16A9) that display very similar effects on urate levels. CONCLUSIONS: We provide the first evidence for association with gout at four loci (IGF1R, PDZK1, MAF, HLF). Understanding why there is lack of correlation between urate and gout effect sizes will be important in understanding the aetiology of gout.

Examining non-syndromic autosomal recessive intellectual disability (NS-ARID) genes for an enriched association with intelligence differences.

Two themes are emerging regarding the molecular genetic aetiology of intelligence. The first is that intelligence is influenced by many variants and those that are tagged by common single nucleotide polymorphisms account for around 30% of the phenotypic variation. The second, in line with other polygenic traits such as height and schizophrenia, is that these variants are not randomly distributed across the genome but cluster in genes that work together. Less clear is whether the very low range of cognitive ability (intellectual disability) is simply one end of the normal distribution describing individual differences in cognitive ability across a population. Here, we examined 40 genes with a known association with non-syndromic autosomal recessive intellectual disability (NS-ARID) to determine if they are enriched for common variants associated with the normal range of intelligence differences. The current study used the 3511 individuals of the Cognitive Ageing Genetics in England and Scotland (CAGES) consortium. In addition, a text mining analysis was used to identify gene sets biologically related to the NS-ARID set. Gene-based tests indicated that genes implicated in NS-ARID were not significantly enriched for quantitative trait loci (QTL) associated with intelligence. These findings suggest that genes in which mutations can have a large and deleterious effect on intelligence are not associated with variation across the range of intelligence differences.

Genetic predisposition to schizophrenia associated with increased use of cannabis.

Cannabis is the most commonly used illicit drug worldwide. With debate surrounding the legalization and control of use, investigating its health risks has become a pressing area of research. One established association is that between cannabis use and schizophrenia, a debilitating psychiatric disorder affecting ~1% of the population over their lifetime. Although considerable evidence implicates cannabis use as a component cause of schizophrenia, it remains unclear whether this is entirely due to cannabis directly raising risk of psychosis, or whether the same genes that increases psychosis risk may also increase risk of cannabis use. In a sample of 2082 healthy individuals, we show an association between an individual's burden of schizophrenia risk alleles and use of cannabis. This was significant both for comparing those who have ever versus never used cannabis (P=2.6 × 10(-4)), and for quantity of use within users (P=3.0 × 10(-3)). Although directly predicting only a small amount of the variance in cannabis use, these findings suggest that part of the association between schizophrenia and cannabis is due to a shared genetic aetiology. This form of gene-environment correlation is an important consideration when calculating the impact of environmental risk factors, including cannabis use.

Genome-wide association study on detailed profiles of smoking behavior and nicotine dependence in a twin sample.

Smoking is a major risk factor for several somatic diseases and is also emerging as a causal factor for neuropsychiatric disorders. Genome-wide association (GWA) and candidate gene studies for smoking behavior and nicotine dependence (ND) have disclosed too few predisposing variants to account for the high estimated heritability. Previous large-scale GWA studies have had very limited phenotypic definitions of relevance to smoking-related behavior, which has likely impeded the discovery of genetic effects. We performed GWA analyses on 1114 adult twins ascertained for ever smoking from the population-based Finnish Twin Cohort study. The availability of 17 smoking-related phenotypes allowed us to comprehensively portray the dimensions of smoking behavior, clustered into the domains of smoking initiation, amount smoked and ND. Our results highlight a locus on 16p12.3, with several single-nucleotide polymorphisms (SNPs) in the vicinity of CLEC19A showing association (P<1 × 10(-6)) with smoking quantity. Interestingly, CLEC19A is located close to a previously reported attention-deficit hyperactivity disorder (ADHD) linkage locus and an evident link between ADHD and smoking has been established. Intriguing preliminary association (P<1 × 10(-5)) was detected between DSM-IV (Diagnostic and Statistical Manual of Mental Disorders, 4th edition) ND diagnosis and several SNPs in ERBB4, coding for a Neuregulin receptor, on 2q33. The association between ERBB4 and DSM-IV ND diagnosis was replicated in an independent Australian sample. Recently, a significant increase in ErbB4 and Neuregulin 3 (Nrg3) expression was revealed following chronic nicotine exposure and withdrawal in mice and an association between NRG3 SNPs and smoking cessation success was detected in a clinical trial. ERBB4 has previously been associated with schizophrenia; further, it is located within an established schizophrenia linkage locus and within a linkage locus for a smoker phenotype identified in this sample. In conclusion, we disclose novel tentative evidence for the involvement of ERBB4 in ND, suggesting the involvement of the Neuregulin/ErbB signalling pathway in addictions and providing a plausible link between the high co-morbidity of schizophrenia and ND.

Childhood intelligence is heritable, highly polygenic and associated with FNBP1L.

Intelligence in childhood, as measured by psychometric cognitive tests, is a strong predictor of many important life outcomes, including educational attainment, income, health and lifespan. Results from twin, family and adoption studies are consistent with general intelligence being highly heritable and genetically stable throughout the life course. No robustly associated genetic loci or variants for childhood intelligence have been reported. Here, we report the first genome-wide association study (GWAS) on childhood intelligence (age range 6-18 years) from 17,989 individuals in six discovery and three replication samples. Although no individual single-nucleotide polymorphisms (SNPs) were detected with genome-wide significance, we show that the aggregate effects of common SNPs explain 22-46% of phenotypic variation in childhood intelligence in the three largest cohorts (P=3.9 × 10(-15), 0.014 and 0.028). FNBP1L, previously reported to be the most significantly associated gene for adult intelligence, was also significantly associated with childhood intelligence (P=0.003). Polygenic prediction analyses resulted in a significant correlation between predictor and outcome in all replication cohorts. The proportion of childhood intelligence explained by the predictor reached 1.2% (P=6 × 10(-5)), 3.5% (P=10(-3)) and 0.5% (P=6 × 10(-5)) in three independent validation cohorts. Given the sample sizes, these genetic prediction results are consistent with expectations if the genetic architecture of childhood intelligence is like that of body mass index or height. Our study provides molecular support for the heritability and polygenic nature of childhood intelligence. Larger sample sizes will be required to detect individual variants with genome-wide significance.

The future for genetic studies in reproduction.

Genetic factors contribute to risk of many common diseases affecting reproduction and fertility. In recent years, methods for genome-wide association studies (GWAS) have revolutionized gene discovery for common traits and diseases. Results of GWAS are documented in the Catalog of Published Genome-Wide Association Studies at the National Human Genome Research Institute and report over 70 publications for 32 traits and diseases associated with reproduction. These include endometriosis, uterine fibroids, age at menarche and age at menopause. Results that pass appropriate stringent levels of significance are generally well replicated in independent studies. Examples of genetic variation affecting twinning rate, infertility, endometriosis and age at menarche demonstrate that the spectrum of disease-related variants for reproductive traits is similar to most other common diseases. GWAS 'hits' provide novel insights into biological pathways and the translational value of these studies lies in discovery of novel gene targets for biomarkers, drug development and greater understanding of environmental factors contributing to disease risk. Results also show that genetic data can help define sub-types of disease and co-morbidity with other traits and diseases. To date, many studies on reproductive traits have used relatively small samples. Future genetic marker studies in large samples with detailed phenotypic and clinical information will yield new insights into disease risk, disease classification and co-morbidity for many diseases associated with reproduction and infertility.

A genome-wide association study of alcohol-dependence symptom counts in extended pedigrees identifies C15orf53.

Several studies have identified genes associated with alcohol-use disorders (AUDs), 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 to identify novel genes affecting risk for alcohol dependence (AD). 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 (SC)). Secondary analyses were performed to determine if the single nucleotide polymorphisms (SNPs) associated with SC were also associated with the dichotomous phenotype, DSM-IV AD. This family-based GWAS identified SNPs in C15orf53 that are strongly associated with DSM-IV alcohol-dependence symptom counts (P=4.5 × 10(-8), inflation-corrected P=9.4 × 10(-7)). Results with DSM-IV AD in the regions of interest support our findings with SC, although the associations were less significant. Attempted replications of the most promising association results were conducted in two independent samples: nonoverlapping subjects from the Study of Addiction: Genes and Environment (SAGE) and the Australian Twin Family Study of AUDs (OZALC). Nominal association of C15orf53 with SC was observed in SAGE. The variant that showed strongest association with SC, rs12912251 and its highly correlated variants (D'=1, r(2) 0.95), have previously been associated with risk for bipolar disorder.

A new regulatory variant in the interleukin-6 receptor gene associates with asthma risk.

The main genetic determinant of soluble interleukin 6 receptor (sIL-6R) levels is the missense variant rs2228145 that maps to the cleavage site of IL-6R. For each Ala allele, sIL-6R serum levels increase by ≈ 20 ng ml(-1) and asthma risk by 1.09-fold. However, this variant does not explain the total heritability for sIL-6R levels. Additional independent variants in IL6R may therefore contribute to variation in sIL-6R levels and influence asthma risk. We imputed 471 variants in IL6R and tested these for association with sIL-6R serum levels in 360 individuals. An intronic variant (rs12083537) was associated with sIL-6R levels independently of rs4129267 (P=0.0005), a proxy single-nucleotide polymorphism for rs2228145. A significant and consistent association for rs12083537 was observed in a replication panel of 354 individuals (P=0.033). Each rs12083537:A allele increased sIL-6R serum levels by 2.4 ng ml(-1). Analysis of mRNA levels in two cohorts did not identify significant associations between rs12083537 and IL6R transcription levels. On the other hand, results from 16,705 asthmatics and 30,809 controls showed that the rs12083537:A allele increased asthma risk by 1.04-fold (P=0.0419). Genetic risk scores based on IL6R regulatory variants may prove useful in explaining variation in clinical response to tocilizumab, an anti-IL-6R monoclonal antibody.

Genome-wide association study of major depressive disorder: new results, meta-analysis, and lessons learned.

Major depressive disorder (MDD) is a common complex disorder with a partly genetic etiology. We conducted a genome-wide association study of the MDD2000+ sample (2431 cases, 3673 screened controls and >1 M imputed single-nucleotide polymorphisms (SNPs)). No SNPs achieved genome-wide significance either in the MDD2000+ study, or in meta-analysis with two other studies totaling 5763 cases and 6901 controls. These results imply that common variants of intermediate or large effect do not have main effects in the genetic architecture of MDD. Suggestive but notable results were (a) gene-based tests suggesting roles for adenylate cyclase 3 (ADCY3, 2p23.3) and galanin (GAL, 11q13.3); published functional evidence relates both of these to MDD and serotonergic signaling; (b) support for the bipolar disorder risk variant SNP rs1006737 in CACNA1C (P=0.020, odds ratio=1.10); and (c) lack of support for rs2251219, a SNP identified in a meta-analysis of affective disorder studies (P=0.51). We estimate that sample sizes 1.8- to 2.4-fold greater are needed for association studies of MDD compared with those for schizophrenia to detect variants that explain the same proportion of total variance in liability. Larger study cohorts characterized for genetic and environmental risk factors accumulated prospectively are likely to be needed to dissect more fully the etiology of MDD.

Genome-wide association analysis of coffee drinking suggests association with CYP1A1/CYP1A2 and NRCAM.

Coffee consumption is a model for addictive behavior. We performed a meta-analysis of genome-wide association studies (GWASs) on coffee intake from 8 Caucasian cohorts (N=18 176) and sought replication of our top findings in a further 7929 individuals. We also performed a gene expression analysis treating different cell lines with caffeine. Genome-wide significant association was observed for two single-nucleotide polymorphisms (SNPs) in the 15q24 region. The two SNPs rs2470893 and rs2472297 (P-values=1.6 × 10(-11) and 2.7 × 10(-11)), which were also in strong linkage disequilibrium (r(2)=0.7) with each other, lie in the 23-kb long commonly shared 5' flanking region between CYP1A1 and CYP1A2 genes. CYP1A1 was found to be downregulated in lymphoblastoid cell lines treated with caffeine. CYP1A1 is known to metabolize polycyclic aromatic hydrocarbons, which are important constituents of coffee, whereas CYP1A2 is involved in the primary metabolism of caffeine. Significant evidence of association was also detected at rs382140 (P-value=3.9 × 10(-09)) near NRCAM-a gene implicated in vulnerability to addiction, and at another independent hit rs6495122 (P-value=7.1 × 10(-09))-an SNP associated with blood pressure-in the 15q24 region near the gene ULK3, in the meta-analysis of discovery and replication cohorts. Our results from GWASs and expression analysis also strongly implicate CAB39L in coffee drinking. Pathway analysis of differentially expressed genes revealed significantly enriched ubiquitin proteasome (P-value=2.2 × 10(-05)) and Parkinson's disease pathways (P-value=3.6 × 10(-05)).