Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci.

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWASs) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N = 52) and nonsmokers (N = 171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and correcting for multiple testing using a two-stage procedure. We found >2 million significant meQTL variants (padj < 0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects, and five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTL variants for 958 unique eGenes (padj < 0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTN1 and ITIH4 colocalized across all data types (GWAS, meQTL, and eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Smoking-informed methylation and expression QTLs in human brain and colocalization with smoking-associated genetic loci.

Smoking is a leading cause of preventable morbidity and mortality. Smoking is heritable, and genome-wide association studies (GWAS) of smoking behaviors have identified hundreds of significant loci. Most GWAS-identified variants are noncoding with unknown neurobiological effects. We used genome-wide genotype, DNA methylation, and RNA sequencing data in postmortem human nucleus accumbens (NAc) to identify cis-methylation/expression quantitative trait loci (meQTLs/eQTLs), investigate variant-by-cigarette smoking interactions across the genome, and overlay QTL evidence at smoking GWAS-identified loci to evaluate their regulatory potential. Active smokers (N=52) and nonsmokers (N=171) were defined based on cotinine biomarker levels and next-of-kin reporting. We simultaneously tested variant and variant-by-smoking interaction effects on methylation and expression, separately, adjusting for biological and technical covariates and using a two-stage multiple testing approach with eigenMT and Bonferroni corrections. We found >2 million significant meQTL variants (padj<0.05) corresponding to 41,695 unique CpGs. Results were largely driven by main effects; five meQTLs, mapping to NUDT12, FAM53B, RNF39, and ADRA1B, showed a significant interaction with smoking. We found 57,683 significant eQTLs for 958 unique eGenes (padj<0.05) and no smoking interactions. Colocalization analyses identified loci with smoking-associated GWAS variants that overlapped meQTLs/eQTLs, suggesting that these heritable factors may influence smoking behaviors through functional effects on methylation/expression. One locus containing MUSTIN1 and ITIH4 colocalized across all data types (GWAS + meQTL + eQTL). In this first genome-wide meQTL map in the human NAc, the enriched overlap with smoking GWAS-identified genetic loci provides evidence that gene regulation in the brain helps explain the neurobiology of smoking behaviors.

Sleep apnea and pesticide exposure in a study of US farmers.

INTRODUCTION: Carbamate and organophosphate pesticides inhibit acetylcholinesterase, and poisoning leads to respiratory depression. Thus, involvement in sleep apnea is plausible, but no data exist at lower levels of exposure. Other pesticides could impact sleep apnea by different mechanisms but have not been studied. Our study examines the associations between pesticide exposure and sleep apnea among pesticide applicators from a US farming population. PARTICIPANTS AND METHODS: We analyzed data from 1569 male pesticide applicators, mostly farmers, from an asthma case-control study nested within the prospective Agricultural Health Study. On questionnaires, participants reported use of specific pesticides and physician diagnosis plus prescribed treatments for sleep apnea. We used multivariable logistic regression to estimate associations between ever use of 63 pesticides and sleep apnea (234 cases, 1335 noncases). RESULTS: The most notable association was for carbofuran, a carbamate (100 exposed cases, odds ratio 1.83, 95% confidence interval 1.34-2.51, P=.0002). Carbofuran use began before reported onset of sleep apnea in all cases. DISCUSSION: This study adds to the known adverse health outcomes of exposure to carbofuran, a pesticide canceled in the United States in 2009 for most agricultural purposes but persists in the environment and remains in use in some other countries. CONCLUSIONS: We conducted the first epidemiological study investigating the association of pesticide exposure and sleep apnea. Our results in a male agricultural population suggests that exposure to carbofuran is positively associated with sleep apnea.

Transcriptome analysis of adult and fetal trabecular meshwork, cornea, and ciliary body tissues by RNA sequencing.

PURPOSE: To characterize the transcriptional landscape of human adult and fetal trabecular meshwork (TM), cornea, and ciliary body (CB) tissues, and to evaluate the expression level of candidate genes selected from genetic association studies of primary-open angle glaucoma, central corneal thickness, intraocular pressure, vertical cup to disc ratio, and optic nerve parameters. METHODS: Deep RNA sequencing was performed on the selected human tissues. Transcriptome analyses were performed to 1) characterize the total number of expressed genes, 2) identify the most highly expressed genes, 3) estimate the number of novel transcripts, and 4) evaluate the expression of candidate genes in each tissue. Finally, a differential gene expression analysis was conducted to compare the adult and fetal ocular tissues. RESULTS: There was an average of 12,362 protein coding genes and 3725 novel transcripts expressed in each tissue. The top most expressed genes in each tissue included SPARC (fetal cornea and TM), APOD (adult TM), CLU (adult cornea), and PTGDS (adult and fetal CB). Twenty-nine candidate genes selected from genetic association studies primarily showed high expression levels in the trabecular meshwork and cornea. Comparison of adult and fetal samples identified 2012 and 1261 differentially expressed protein-coding genes within the cornea and trabecular meshwork, respectively. CONCLUSIONS: This study has provided an unbiased glimpse into the transcriptome of three essential anterior ocular tissues, resulting in the development of several novel hypotheses. These data can be used in the future to better guide ocular research questions.

House Dust Endotoxin Levels Are Associated with Adult Asthma in a U.S. Farming Population.

RATIONALE: Endotoxin initiates a proinflammatory response from the innate immune system. Studies in children suggest that endotoxin exposure from house dust may be an important risk factor for asthma, but few studies have been conducted in adult populations. OBJECTIVES: To investigate the association of house dust endotoxin levels with asthma and related phenotypes (wheeze, atopy, and pulmonary function) in a large U.S. farming population. METHODS: Dust was collected from the bedrooms (n = 2,485) of participants enrolled in a case-control study of current asthma (927 cases) nested within the Agricultural Health Study. Dust endotoxin was measured by Limulus amebocyte lysate assay. Outcomes were measured by questionnaire, spirometry, and blood draw. We evaluated associations using linear and logistic regression. MEASUREMENTS AND MAIN RESULTS: Endotoxin was significantly associated with current asthma (odds ratio [OR], 1.30; 95% confidence interval [CI], 1.14-1.47), and this relationship was modified by early-life farm exposure (born on a farm: OR, 1.18; 95% CI, 1.02-1.37; not born on a farm: OR, 1.67; 95% CI, 1.26-2.20; Interaction P = 0.05). Significant positive associations were seen with both atopic and nonatopic asthma. Endotoxin was not related to either atopy or wheeze. Higher endotoxin was related to lower FEV1/FVC in asthma cases only (Interaction P = 0.01). For asthma, there was suggestive evidence of a gene-by-environment interaction for the CD14 variant rs2569190 (Interaction P = 0.16) but not for the TLR4 variants rs4986790 and rs4986791. CONCLUSIONS: House dust endotoxin was associated with current atopic and nonatopic asthma in a U.S. farming population. The degree of the association with asthma depended on early-life farm exposures. Furthermore, endotoxin was associated with lower pulmonary function in patients with asthma.

The Genomic Basis of Postponed Senescence in Drosophila melanogaster.

Natural populations harbor considerable genetic variation for lifespan. While evolutionary theory provides general explanations for the existence of this variation, our knowledge of the genes harboring naturally occurring polymorphisms affecting lifespan is limited. Here, we assessed the genetic divergence between five Drosophila melanogaster lines selected for postponed senescence for over 170 generations (O lines) and five lines from the same base population maintained at a two week generation interval for over 850 generations (B lines). On average, O lines live 70% longer than B lines, are more productive at all ages, and have delayed senescence for other traits than reproduction. We performed population sequencing of pools of individuals from all B and O lines and identified 6,394 genetically divergent variants in or near 1,928 genes at a false discovery rate of 0.068. A 2.6 Mb region at the tip of the X chromosome contained many variants fixed for alternative alleles in the two populations, suggestive of a hard selective sweep. We also assessed genome wide gene expression of O and B lines at one and five weeks of age using RNA sequencing and identified genes with significant (false discovery rate < 0.05) effects on gene expression with age, population and the age by population interaction, separately for each sex. We identified transcripts that exhibited the transcriptional signature of postponed senescence and integrated the gene expression and genetic divergence data to identify 98 (175) top candidate genes in females (males) affecting postponed senescence and increased lifespan. While several of these genes have been previously associated with Drosophila lifespan, most are novel and constitute a rich resource for future functional validation.

Discovery and functional annotation of SIX6 variants in primary open-angle glaucoma.

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10(-11)). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10(-10)) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.