Expression Levels Of Nearby Genes Research Articles

Mitochondrial DNA Haplogroup K Is Protective Against Autism Spectrum Disorder Risk in Populations of European Ancestry

Accumulative evidence indicates a critical role of mitochondrial function in autism spectrum disorders (ASD), implying that ASD risk may be linked to mitochondrial dysfunction due to DNA (mtDNA) variations. Although a few studies have explored the association between mtDNA variations and ASD, the role of mtDNA in ASD is still unclear. Here, we aimed to investigate whether mitochondrial DNA haplogroups are associated with the risk of ASD. Two European cohorts and an Ashkenazi Jewish (AJ) cohort were analyzed, including 2,062 ASD patients in comparison with 4,632 healthy controls. DNA samples were genotyped using Illumina HumanHap550/610 and Illumina 1M arrays, inclusive of mitochondrial markers. Mitochondrial DNA (mtDNA) haplogroups were identified from genotyping data using HaploGrep2. A mitochondrial genome imputation pipeline was established to detect mtDNA variants. We conducted a case-control study to investigate potential associations of mtDNA haplogroups and variants with the susceptibility to ASD. We observed that the ancient adaptive mtDNA haplogroup K was significantly associated with decreased risk of ASD by the investigation of 2 European cohorts including a total of 2,006 cases and 4,435 controls (odds ratio= 0.64, P=1.79× 10-5), and we replicated this association in an Ashkenazi Jewish (AJ) cohort including 56 cases and 197 controls (odds ratio= 0.35, P= 9.46× 10-3). Moreover, we demonstrate that the mtDNA variants rs28358571, rs28358584, and rs28358280 are significantly associated with ASD risk. Further expression quantitative trait loci (eQTLs) analysis indicated that the rs28358584 and rs28358280 genotypes are associated with expression levels of nearby genes in brain tissues, suggesting those mtDNA variants may confer risk for ASD via regulation of expression levels of genes encoded by the mitochondrial genome. This study helps to shed light on the contribution of mitochondria in ASD and provides new insights into the genetic mechanism underlying ASD, suggesting the potential involvement of mtDNA-encoded proteins in the development of ASD. Increasing evidence indicates that mitochondrial dysfunction may be linked to autism spectrum disorder (ASD). This study investigated potential associations of mitochondrial DNA (mtDNA) variants in 2 European and Ashkenazi Jewish cohorts including 2,062 individuals with ASD and 4,632 healthy controls. Researchers found that the ancient mtDNA haplogroup K was linked to a reduced risk of ASD in both European and Ashkenazi Jewish populations. Additionally, specific mtDNA variants were associated with ASD risk and were shown to influence the expression of nearby genes in the brain. These findings highlight the potential involvement of mtDNA in ASD development, offering new insights into the genetic mechanisms underlying the disorder.

DNA methylation clustered by predicted chromatin state reveals distinct association patterns with neuropathologic and biochemical measures in the AD brain

AbstractBackgroundAlzheimer’s disease (AD) is neuropathologically characterized by amyloid‐β (Aβ) plaques and tau neurofibrillary tangles often quantified by Thal phase and Braak stage, respectively. Aβ also frequently deposits in the cerebrovasculature with severity categorized by a cerebral amyloid angiopathy (CAA) score. These and related measures often show high variability within AD suggesting distinct underlying mechanisms, the identification of which may reveal important biological insights into AD pathophysiology. We hypothesize that, within the AD brain, neuropathology and levels of core AD‐related proteins are influenced by variations in DNA methylation (DNAm). To test this, we performed epigenome‐wide association studies (EWAS) using DNAm clusters from the temporal cortex (TCX) and cerebellum (CER) with AD‐related neuropathologic measures (Braak, Thal, and CAA) and brain biochemical levels of five proteins (apoE, Aβ40, Aβ42, tau, and p‐tau).MethodsDNAm from neuropathologically‐confirmed AD cases available through the Mayo Clinic brain bank was measured by reduced representation bisulfite sequencing (RRBS) from 471 brain samples, 200 of which had both TCX and CER measures. Neuropathology was scored by an experienced neuropathologist. TCX levels of five AD‐related proteins from three tissue fractions (buffer‐, detergent‐, and in‐soluble) were measured by ELISA (Liu 2020). CpG methylation (CpGm) was binned by a 15‐state chromatin model (Kundaje 2015) and averaged. Averaged CpGm clusters were tested for their association with each endophenotype through multi‐variable linear regression. The influence of CpGm clusters on expression levels of nearby genes was also investigated through multi‐variable regression.ResultsOur innovative binning method demonstrated expected and biologically relevant CpGm patterns. We identified multiple epigenome‐wide significant CpGm clusters unique to each AD‐related endophenotype and brain region. Comparisons between the TCX and CER further revealed brain region‐specific enrichment of distinct CpGm clusters. Significant CpGm clusters were further characterized by their influence on gene expression.ConclusionsAlthough all neuropathologic and biochemical endophenotypes tested are core to AD pathophysiology, our results suggest each has a distinct epigenetic architecture underlying their variability in the AD brain. By discovering brain region‐specific DNAm changes that contribute to each neuropathologic or biochemical outcome, we can identify core components of complex mechanisms that individually and jointly contribute to AD.

Novel Genetic Variants Associated with Chronic Kidney Disease Progression.

eGFR slope has been used as a surrogate outcome for progression of CKD. However, genetic markers associated with eGFR slope among patients with CKD were unknown. We aimed to identify genetic susceptibility loci associated with eGFR slope. A two-phase genome-wide association study identified single nucleotide polymorphisms (SNPs) in TPPP and FAT1-LINC02374 , and 22 of them were used to derive polygenic risk scores that mark the decline of eGFR by disrupting binding of nearby transcription factors. This work is the first to identify the impact of TPPP and FAT1-LINC02374 on CKD progression, providing predictive markers for the decline of eGFR in patients with CKD. The incidence of CKD is associated with genetic factors. However, genetic markers associated with the progression of CKD have not been fully elucidated. We conducted a genome-wide association study among 1738 patients with CKD, mainly from the KoreaN cohort study for Outcomes in patients With CKD. The outcome was eGFR slope. We performed a replication study for discovered single nucleotide polymorphisms (SNPs) with P <10 -6 in 2498 patients with CKD from the Chronic Renal Insufficiency Cohort study. Several expression quantitative trait loci (eQTL) studies, pathway enrichment analyses, exploration of epigenetic architecture, and predicting disruption of transcription factor (TF) binding sites explored potential biological implications of the loci. We developed and evaluated the effect of polygenic risk scores (PRS) on incident CKD outcomes. SNPs in two novel loci, TPPP and FAT1-LINC02374 , were replicated (rs59402340 in TPPP , Pdiscovery =7.11×10 -7 , PCRIC =8.13×10 -4 , Pmeta =7.23×10 -8 ; rs28629773 in FAT1-LINC02374 , Pdiscovery =6.08×10 -7 , PCRIC =4.33×10 -2 , Pmeta =1.87×10 -7 ). The eQTL studies revealed that the replicated SNPs regulated the expression level of nearby genes associated with kidney function. Furthermore, these SNPs were near gene enhancer regions and predicted to disrupt the binding of TFs. PRS based on the independently significant top 22 SNPs were significantly associated with CKD outcomes. This study demonstrates that SNP markers in the TPPP and FAT1-LINC02374 loci could be predictive markers for the decline of eGFR in patients with CKD.

Genome-wide evaluation of the effect of short tandem repeat variation on local DNA methylation.

Short tandem repeats (STRs) contribute significantly to genetic diversity in humans, including disease-causing variation. Although the effect of STR variation on gene expression has been extensively assessed, their impact on epigenetics has been poorly studied and limited to specific genomic regions. Here, we investigated the hypothesis that some STRs act as independent regulators of local DNA methylation in the human genome and modify risk of common human traits. To address these questions, we first analyzed two independent data sets comprising PCR-free whole-genome sequencing (WGS) and genome-wide DNA methylation levels derived from whole-blood samples in 245 (discovery cohort) and 484 individuals (replication cohort). Using genotypes for 131,635 polymorphic STRs derived from WGS using HipSTR, we identified 11,870 STRs that associated with DNA methylation levels (mSTRs) of 11,774 CpGs (Bonferroni P < 0.001) in our discovery cohort, with 90% successfully replicating in our second cohort. Subsequently, through fine-mapping using CAVIAR we defined 585 of these mSTRs as the likely causal variants underlying the observed associations (fm-mSTRs) and linked a fraction of these to previously reported genome-wide association study signals, providing insights into the mechanisms underlying complex human traits. Furthermore, by integrating gene expression data, we observed that 12.5% of the tested fm-mSTRs also modulate expression levels of nearby genes, reinforcing their regulatory potential. Overall, our findings expand the catalog of functional sequence variants that affect genome regulation, highlighting the importance of incorporating STRs in future genetic association analysis and epigenetics data for the interpretation of trait-associated variants.

Genome-wide Interaction Study with Smoking for Colorectal Cancer Risk Identifies Novel Genetic Loci Related to Tumor Suppression, Inflammation, and Immune Response.

Tobacco smoking is an established risk factor for colorectal cancer. However, genetically defined population subgroups may have increased susceptibility to smoking-related effects on colorectal cancer. A genome-wide interaction scan was performed including 33,756 colorectal cancer cases and 44,346 controls from three genetic consortia. Evidence of an interaction was observed between smoking status (ever vs. never smokers) and a locus on 3p12.1 (rs9880919, P = 4.58 × 10-8), with higher associated risk in subjects carrying the GG genotype [OR, 1.25; 95% confidence interval (CI), 1.20-1.30] compared with the other genotypes (OR <1.17 for GA and AA). Among ever smokers, we observed interactions between smoking intensity (increase in 10 cigarettes smoked per day) and two loci on 6p21.33 (rs4151657, P = 1.72 × 10-8) and 8q24.23 (rs7005722, P = 2.88 × 10-8). Subjects carrying the rs4151657 TT genotype showed higher risk (OR, 1.12; 95% CI, 1.09-1.16) compared with the other genotypes (OR <1.06 for TC and CC). Similarly, higher risk was observed among subjects carrying the rs7005722 AA genotype (OR, 1.17; 95% CI, 1.07-1.28) compared with the other genotypes (OR <1.13 for AC and CC). Functional annotation revealed that SNPs in 3p12.1 and 6p21.33 loci were located in regulatory regions, and were associated with expression levels of nearby genes. Genetic models predicting gene expression revealed that smoking parameters were associated with lower colorectal cancer risk with higher expression levels of CADM2 (3p12.1) and ATF6B (6p21.33). Our study identified novel genetic loci that may modulate the risk for colorectal cancer of smoking status and intensity, linked to tumor suppression and immune response. These findings can guide potential prevention treatments.

Genome-wide contribution of common short-tandem repeats to Parkinson's disease genetic risk.

Parkinson's disease is a complex neurodegenerative disorder with a strong genetic component, for which most known disease-associated variants are single nucleotide polymorphisms (SNPs) and small insertions and deletions (indels). DNA repetitive elements account for >50% of the human genome; however, little is known of their contribution to Parkinson's disease aetiology. While select short tandem repeats (STRs) within candidate genes have been studied in Parkinson's disease, their genome-wide contribution remains unknown. Here we present the first genome-wide association study of STRs in Parkinson's disease. Through a meta-analysis of 16 imputed genome-wide association study cohorts from the International Parkinson's Disease Genomic Consortium (IPDGC), totalling 39 087 individuals (16 642 cases and 22 445 controls of European ancestry), we identified 34 genome-wide significant STR loci (P < 5.34 × 10-6), with the strongest signal located in KANSL1 [chr17:44 205 351:[T]11, P = 3 × 10-39, odds ratio = 1.31 (95% confidence interval = 1.26-1.36)]. Conditional-joint analyses suggested that four significant STRs mapping nearby NDUFAF2, TRIML2, MIRNA-129-1 and NCOR1 were independent from known risk SNPs. Including STRs in heritability estimates increased the variance explained by SNPs alone. Gene expression analysis of STRs (eSTRs) in RNA sequencing data from 13 brain regions identified significant associations of STRs influencing the expression of multiple genes, including known Parkinson's disease genes. Further functional annotation of candidate STRs revealed that significant eSTRs within NUDFAF2 and ZSWIM7 overlap with regulatory features and are associated with change in the expression levels of nearby genes. Here, we show that STRs at known and novel candidate loci contribute to Parkinson's disease risk and have functional effects in disease-relevant tissues and pathways, supporting previously reported disease-associated genes and giving further evidence for their functional prioritization. These data represent a valuable resource for researchers currently dissecting Parkinson's disease risk loci.

Convergent lines of evidence support BIN1 as a risk gene of Alzheimer\u2019s disease

Genome-wide association studies (GWAS) have identified several susceptibility loci of Alzheimer’s disease (AD), which were mainly located in noncoding regions of the genome. Meanwhile, the putative biological mechanisms underlying AD susceptibility loci were still unclear. At present, identifying the functional variants of AD pathogenesis remains a major challenge. Herein, we first used summary data-based Mendelian randomization (SMR) with AD GWAS summary and expression quantitative trait loci (eQTL) data to identify variants who affects expression levels of nearby genes and contributed to the risk of AD. Using the SMR integrative analysis, we totally identified 14 SNPs significantly affected the expression level of 16 nearby genes in blood or brain tissues and contributed to the AD risk. Then, to confirm the results, we replicated the GWAS and eQTL results across multiple samples. Totally, four risk SNP (rs11682128, rs601945, rs3935067, and rs679515) were validated to be associated with AD and affected the expression level of nearby genes (BIN1, HLA-DRA, EPHA1-AS1, and CR1). Besides, our differential expression analysis showed that the BIN1 gene was significantly downregulated in the hippocampus (P = 2.0 × 10−3) and survived after multiple comparisons. These convergent lines of evidence suggest that the BIN1 gene identified by SMR has potential roles in the pathogenesis of AD. Further investigation of the roles of the BIN1 gene in the pathogenesis of AD is warranted.

Association of Variants in PLD1, 3p24.1, and 10q11.21 Regions With Hirschsprung’s Disease in Han Chinese Population

Background and Aims: Hirschsprung’s disease (HSCR) is a rare genetically heterogeneous congenital disorder. A recent study based on whole genome sequencing demonstrated that common variants at four novel loci, which contained two intronic variants on CASQ2 and PLD1, and intergenic variants located between SLC4A7 and EOMES at 3p24.1, and between LINC01518 and LOC283028 at 10q11.21, were associated with HSCR susceptibility. To validate these associations with HSCR susceptibility, we performed a case–control study in a Han Chinese sample set.Methods: We selected four previously identified single nucleotide polymorphisms (SNPs) for replication, along with tag SNPs to cover the four associated regions. In total, 61 SNPs were genotyped in 420 HSCR patients and 1,665 healthy controls from the Han Chinese population.Results: None of the 14 tag SNPs in the CASQ2 gene region, including the previously associated rs9428225, showed an association with HSCR. Among the 24 tag SNPs from the SLC4A7-EOMES region at 3p24.1, rs2642925 [odds ratio (OR) = 1.41, 95% confidence interval (95% CI) = 1.10–1.79; PAdditive = 0.007] and the previously associated SNP rs9851320 showed a suggestive association (OR = 1.22, 95% CI = 1.01–1.47; PAdditive = 0.042). A non-synonymous SNP, rs2287579, in PLD1 showed a suggestive association with HSCR susceptibility (OR = 1.71, 95% CI = 1.18–2.46; PAdditive = 0.004). Additionally, the previously associated PLD1 SNP rs12632766 showed a suggestive significance (OR = 1.20, 95% CI = 1.01–1.42, PAdditive = 0.038). In the LINC01518-LOC283028 region at 10q11.21, three SNPs meet the study-wide significance threshold. Rs17153309 was the most associated SNP (OR = 1.60, 95% CI = 1.34–1.90; PAdditive = 1.13 × 10–7). The previously associated SNP rs1414027 also showed significant association (OR = 1.43, 95% CI = 1.20–1.70, PAdditive = 3.92 × 10–5). Two associated SNPs at 10q11.21 (rs1414027 and rs624804) were expression quantitative trait loci in digestive tract tissues from GTEx databases.Conclusions: Our results confirmed that variants of the LINC01518-LOC283028 region were associated with HSCR in the Han Chinese population. Additionally, the susceptibility of SNPs in the LINC01518-LOC283028 region were associated with the expression levels of nearby genes. These results provide new insight into the pathogenesis of HSCR.

Identification of novel common variants associated with chronic pain using conditional false discovery rate analysis with major depressive disorder and assessment of pleiotropic effects of LRFN5

Chronic pain is a complex trait that is moderately heritable and genetically, as well as phenotypically, correlated with major depressive disorder (MDD). Use of the conditional false discovery rate (cFDR) approach, which leverages pleiotropy identified from existing GWAS outputs, has been successful in discovering novel associated variants in related phenotypes. Here, genome-wide association study outputs for both von Korff chronic pain grade and for MDD were used to identify variants meeting a cFDR threshold for each outcome phenotype separately, as well as a conjunctional cFDR (ccFDR) threshold for both phenotypes together. Using a moderately conservative threshold, we identified a total of 11 novel single nucleotide polymorphisms (SNPs), six of which were associated with chronic pain grade and nine of which were associated with MDD. Four SNPs on chromosome 14 were associated with both chronic pain grade and MDD. SNPs associated only with chronic pain grade were located within SLC16A7 on chromosome 12. SNPs associated only with MDD were located either in a gene-dense region on chromosome 1 harbouring LINC01360, LRRIQ3, FPGT and FPGT-TNNI3K, or within/close to LRFN5 on chromosome 14. The SNPs associated with both outcomes were also located within LRFN5. Several of the SNPs on chromosomes 1 and 14 were identified as being associated with expression levels of nearby genes in the brain and central nervous system. Overall, using the cFDR approach, we identified several novel genetic loci associated with chronic pain and we describe likely pleiotropic effects of a recently identified MDD locus on chronic pain.

Genome-wide DNA methylation and long-term ambient air pollution exposure in Korean adults

BackgroundAmbient air pollution is associated with numerous adverse health outcomes, but the underlying mechanisms are not well understood; epigenetic effects including altered DNA methylation could play a role. To evaluate associations of long-term air pollution exposure with DNA methylation in blood, we conducted an epigenome-wide association study in a Korean chronic obstructive pulmonary disease cohort (N = 100 including 60 cases) using Illumina’s Infinium HumanMethylation450K Beadchip. Annual average concentrations of particulate matter ≤ 10 μm in diameter (PM10) and nitrogen dioxide (NO2) were estimated at participants’ residential addresses using exposure prediction models. We used robust linear regression to identify differentially methylated probes (DMPs) and two different approaches, DMRcate and comb-p, to identify differentially methylated regions (DMRs).ResultsAfter multiple testing correction (false discovery rate < 0.05), there were 12 DMPs and 27 DMRs associated with PM10 and 45 DMPs and 57 DMRs related to NO2. DMP cg06992688 (OTUB2) and several DMRs were associated with both exposures. Eleven DMPs in relation to NO2 confirmed previous findings in Europeans; the remainder were novel. Methylation levels of 39 DMPs were associated with expression levels of nearby genes in a separate dataset of 3075 individuals. Enriched networks were related to outcomes associated with air pollution including cardiovascular and respiratory diseases as well as inflammatory and immune responses.ConclusionsThis study provides evidence that long-term ambient air pollution exposure impacts DNA methylation. The differential methylation signals can serve as potential air pollution biomarkers. These results may help better understand the influences of ambient air pollution on human health.

New Blood Pressure-Associated Loci Identified in Meta-Analyses of 475 000 Individuals.

Genome-wide association studies have recently identified >400 loci that harbor DNA sequence variants that influence blood pressure (BP). Our earlier studies identified and validated 56 single nucleotide variants (SNVs) associated with BP from meta-analyses of exome chip genotype data. An additional 100 variants yielded suggestive evidence of association. Here, we augment the sample with 140 886 European individuals from the UK Biobank, in whom 77 of the 100 suggestive SNVs were available for association analysis with systolic BP or diastolic BP or pulse pressure. We performed 2 meta-analyses, one in individuals of European, South Asian, African, and Hispanic descent (pan-ancestry, ≈475 000), and the other in the subset of individuals of European descent (≈423 000). Twenty-one SNVs were genome-wide significant (P<5×10-8) for BP, of which 4 are new BP loci: rs9678851 (missense, SLC4A1AP), rs7437940 (AFAP1), rs13303 (missense, STAB1), and rs1055144 (7p15.2). In addition, we identified a potentially independent novel BP-associated SNV, rs3416322 (missense, SYNPO2L) at a known locus, uncorrelated with the previously reported SNVs. Two SNVs are associated with expression levels of nearby genes, and SNVs at 3 loci are associated with other traits. One SNV with a minor allele frequency <0.01, (rs3025380 at DBH) was genome-wide significant. We report 4 novel loci associated with BP regulation, and 1 independent variant at an established BP locus. This analysis highlights several candidate genes with variation that alter protein function or gene expression for potential follow-up.

FIRE: functional inference of genetic variants that regulate gene expression

Interpreting genetic variation in noncoding regions of the genome is an important challenge for personal genome analysis. One mechanism by which noncoding single nucleotide variants (SNVs) influence downstream phenotypes is through the regulation of gene expression. Methods to predict whether or not individual SNVs are likely to regulate gene expression would aid interpretation of variants of unknown significance identified in whole-genome sequencing studies. We developed FIRE (Functional Inference of Regulators of Expression), a tool to score both noncoding and coding SNVs based on their potential to regulate the expression levels of nearby genes. FIRE consists of 23 random forests trained to recognize SNVs in cis-expression quantitative trait loci (cis-eQTLs) using a set of 92 genomic annotations as predictive features. FIRE scores discriminate cis-eQTL SNVs from non-eQTL SNVs in the training set with a cross-validated area under the receiver operating characteristic curve (AUC) of 0.807, and discriminate cis-eQTL SNVs shared across six populations of different ancestry from non-eQTL SNVs with an AUC of 0.939. FIRE scores are also predictive of cis-eQTL SNVs across a variety of tissue types. FIRE scores for genome-wide SNVs in hg19/GRCh37 are available for download at https://sites.google.com/site/fireregulatoryvariation/. nilah@stanford.edu. Supplementary data are available at Bioinformatics online.

Bayesian association scan reveals loci associated with human lifespan and linked biomarkers

The enormous variation in human lifespan is in part due to a myriad of sequence variants, only a few of which have been revealed to date. Since many life-shortening events are related to diseases, we developed a Mendelian randomization-based method combining 58 disease-related GWA studies to derive longevity priors for all HapMap SNPs. A Bayesian association scan, informed by these priors, for parental age of death in the UK Biobank study (n=116,279) revealed 16 independent SNPs with significant Bayes factor at a 5% false discovery rate (FDR). Eleven of them replicate (5% FDR) in five independent longevity studies combined; all but three are depleted of the life-shortening alleles in older Biobank participants. Further analysis revealed that brain expression levels of nearby genes (RBM6, SULT1A1 and CHRNA5) might be causally implicated in longevity. Gene expression and caloric restriction experiments in model organisms confirm the conserved role for RBM6 and SULT1A1 in modulating lifespan.

Identification of Susceptibility Loci and Genes for Colorectal Cancer Risk

Known genetic factors explain only a small fraction of genetic variation in colorectal cancer (CRC). We conducted a genome-wide association study to identify risk loci for CRC. This discovery stage included 8027 cases and 22,577 controls of East-Asian ancestry. Promising variants were evaluated in studies including as many as 11,044 cases and 12,047 controls. Tumor-adjacent normal tissues from 188 patients were analyzed to evaluate correlations of risk variants with expression levels of nearby genes. Potential functionality of risk variants were evaluated using public genomic and epigenomic databases. We identified 4 loci associated with CRC risk; P values for the most significant variant in each locus ranged from 3.92× 10(-8) to 1.24× 10(-12): 6p21.1 (rs4711689), 8q23.3 (rs2450115, rs6469656), 10q24.3 (rs4919687), and 12p13.3 (rs11064437). We also identified 2 risk variants at loci previously associated with CRC: 10q25.2 (rs10506868) and 20q13.3 (rs6061231). These risk variants, conferring an approximate 10%-18% increase in risk per allele, are located either inside or near protein-coding genes that include transcription factor EB (lysosome biogenesis and autophagy), eukaryotic translation initiation factor 3, subunit H (initiation of translation), cytochrome P450, family 17, subfamily A, polypeptide 1 (steroidogenesis), splA/ryanodine receptor domain and SOCS box containing 2 (proteasome degradation), and ribosomal protein S2 (ribosome biogenesis). Gene expressionanalyses showed a significant association (P < .05) for rs4711689 with transcription factor EB, rs6469656 with eukaryotic translation initiation factor 3, subunit H, rs11064437 with splA/ryanodine receptor domain and SOCS box containing 2, and rs6061231 with ribosomal protein S2. We identified susceptibility loci and genes associated with CRC risk, linking CRC predisposition to steroid hormone, protein synthesis and degradation, and autophagy pathways and providing added insight into the mechanism of CRC pathogenesis.

Association of Forced Vital Capacity with the Developmental Gene NCOR2.

BackgroundForced Vital Capacity (FVC) is an important predictor of all-cause mortality in the absence of chronic respiratory conditions. Epidemiological evidence highlights the role of early life factors on adult FVC, pointing to environmental exposures and genes affecting lung development as risk factors for low FVC later in life. Although highly heritable, a small number of genes have been found associated with FVC, and we aimed at identifying further genetic variants by focusing on lung development genes.MethodsPer-allele effects of 24,728 SNPs in 403 genes involved in lung development were tested in 7,749 adults from three studies (NFBC1966, ECRHS, EGEA). The most significant SNP for the top 25 genes was followed-up in 46,103 adults (CHARGE and SpiroMeta consortia) and 5,062 children (ALSPAC). Associations were considered replicated if the replication p-value survived Bonferroni correction (p<0.002; 0.05/25), with a nominal p-value considered as suggestive evidence. For SNPs with evidence of replication, effects on the expression levels of nearby genes in lung tissue were tested in 1,111 lung samples (Lung eQTL consortium), with further functional investigation performed using public epigenomic profiling data (ENCODE).ResultsNCOR2-rs12708369 showed strong replication in children (p = 0.0002), with replication unavailable in adults due to low imputation quality. This intronic variant is in a strong transcriptional enhancer element in lung fibroblasts, but its eQTL effects could not be tested due to low imputation quality in the eQTL dataset. SERPINE2-rs6754561 replicated at nominal level in both adults (p = 0.036) and children (p = 0.045), while WNT16-rs2707469 replicated at nominal level only in adults (p = 0.026). The eQTL analyses showed association of WNT16-rs2707469 with expression levels of the nearby gene CPED1. We found no statistically significant eQTL effects for SERPINE2-rs6754561.ConclusionsWe have identified a new gene, NCOR2, in the retinoic acid signalling pathway pointing to a role of vitamin A metabolism in the regulation of FVC. Our findings also support SERPINE2, a COPD gene with weak previous evidence of association with FVC, and suggest WNT16 as a further promising candidate.

Bacterial infection remodels the DNA methylation landscape of human dendritic cells

DNA methylation is an epigenetic mark thought to be robust to environmental perturbations on a short time scale. Here, we challenge that view by demonstrating that the infection of human dendritic cells (DCs) with a live pathogenic bacteria is associated with rapid and active demethylation at thousands of loci, independent of cell division. We performed an integrated analysis of data on genome-wide DNA methylation, histone mark patterns, chromatin accessibility, and gene expression, before and after infection. We found that infection-induced demethylation rarely occurs at promoter regions and instead localizes to distal enhancer elements, including those that regulate the activation of key immune transcription factors. Active demethylation is associated with extensive epigenetic remodeling, including the gain of histone activation marks and increased chromatin accessibility, and is strongly predictive of changes in the expression levels of nearby genes. Collectively, our observations show that active, rapid changes in DNA methylation in enhancers play a previously unappreciated role in regulating the transcriptional response to infection, even in nonproliferating cells.

Methylation QTLs are associated with coordinated changes in transcription factor binding, histone modifications, and gene expression levels.

DNA methylation is an important epigenetic regulator of gene expression. Recent studies have revealed widespread associations between genetic variation and methylation levels. However, the mechanistic links between genetic variation and methylation remain unclear. To begin addressing this gap, we collected methylation data at ∼300,000 loci in lymphoblastoid cell lines (LCLs) from 64 HapMap Yoruba individuals, and genome-wide bisulfite sequence data in ten of these individuals. We identified (at an FDR of 10%) 13,915 cis methylation QTLs (meQTLs)—i.e., CpG sites in which changes in DNA methylation are associated with genetic variation at proximal loci. We found that meQTLs are frequently associated with changes in methylation at multiple CpGs across regions of up to 3 kb. Interestingly, meQTLs are also frequently associated with variation in other properties of gene regulation, including histone modifications, DNase I accessibility, chromatin accessibility, and expression levels of nearby genes. These observations suggest that genetic variants may lead to coordinated molecular changes in all of these regulatory phenotypes. One plausible driver of coordinated changes in different regulatory mechanisms is variation in transcription factor (TF) binding. Indeed, we found that SNPs that change predicted TF binding affinities are significantly enriched for associations with DNA methylation at nearby CpGs.

Nucleosome positioning and histone modifications define relationships between regulatory elements and nearby gene expression in breast epithelial cells.

BackgroundThe precise nature of how cell type specific chromatin structures at enhancer sites affect gene expression is largely unknown. Here we identified cell type specific enhancers coupled with gene expression in two different types of breast epithelial cells, HMEC (normal breast epithelial cells) and MDAMB231 (triple negative breast cancer cell line).ResultsEnhancers were defined by modified neighboring histones [using chromatin immunoprecipitation followed by sequencing (ChIP-seq)] and nucleosome depletion [using formaldehyde-assisted isolation of regulatory elements followed by sequencing (FAIRE-seq)]. Histone modifications at enhancers were related to the expression levels of nearby genes up to 750 kb away. These expression levels were correlated with enhancer status (poised or active), defined by surrounding histone marks. Furthermore, about fifty percent of poised and active enhancers contained nucleosome-depleted regions. We also identified response element motifs enriched at these enhancer sites that revealed key transcription factors (e.g. TP63) likely involved in regulating breast epithelial enhancer-mediated gene expression. By utilizing expression data, potential target genes of more than 600 active enhancers were identified. These genes were involved in proteolysis, epidermis development, cell adhesion, mitosis, cell cycle, and DNA replication.ConclusionsThese findings facilitate the understanding of epigenetic regulation specifically, such as the relationships between regulatory elements and gene expression and generally, how breast epithelial cellular phenotypes are determined by cell type specific enhancers.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-331) contains supplementary material, which is available to authorized users.

Transcriptional regulation and prostate cancer risk loci.

1554 Background: While genome-wide association studies (GWAS) have identified numerous loci at which common single nucleotide polymorphisms (SNPs) are clearly associated with the risk of developing prostate cancer, the mechanism by which these variants influence disease is not clear. Based on the observation that regions of the genome with marks of transcriptional regulatory elements harbor SNPs that are more likely to be under negative selective pressure, we hypothesized that many prostate cancer associated SNPs, or their correlated proxies in linkage disequilibrium, function by altering regulation of nearby genes through alteration of transcription factor binding sites. Two predictions of this hypothesis are that there will be an enrichment of prostate cancer risk SNPs or their proxies in regulatory elements in prostate tissue and that some of these SNPs will correlate with the expression levels of nearby genes. Methods: To test the first prediction, we compared the number of prostate cancer risk SNPs (or their proxies) that overlap with DNase hypersensitive sites in the LNCaP prostate cancer cell line with randomly selected SNPs matched for allele frequency, distance from nearby genes, and local gene density. To test the second prediction, we asked if any of the known prostate cancer risk SNPs are associated with mRNA expression levels of nearby genes. Results: We found a four-fold enrichment of real prostate cancer risk SNPs overlapping with DNase hypersensitive sites versus the random set (p=0.03) at an r2 threshold of 0.95. In both tumor and benign prostate tissue, we found that a recently identified prostate cancer risk SNP is associated with expression levels of IRX4, a transcription factor previously implicated only in cardiac development. Conclusions: These data support the hypothesis that many GWAS-identified risk SNPs alter transcriptional regulation of nearby genes.

Comparison of class 2 transposable elements at superfamily resolution reveals conserved and distinct features in cereal grass genomes

BackgroundClass 2 transposable elements (TEs) are the predominant elements in and around plant genes where they generate significant allelic diversity. Using the complete sequences of four grasses, we have performed a novel comparative analysis of class 2 TEs. To ensure consistent comparative analyses, we re-annotated class 2 TEs in Brachypodium distachyon, Oryza sativa (rice), Sorghum bicolor and Zea mays and assigned them to one of the five cut-and-paste superfamilies found in plant genomes (Tc1/mariner, PIF/Harbinger, hAT, Mutator, CACTA). We have focused on noncoding elements because of their abundance, and compared superfamily copy number, size and genomic distribution as well as correlation with the level of nearby gene expression.ResultsOur comparison revealed both unique and conserved features. First, the average length or size distribution of elements in each superfamily is largely conserved, with the shortest always being Tc1/mariner elements, followed by PIF/Harbinger, hAT, Mutator and CACTA. This order also holds for the ratio of the copy numbers of noncoding to coding elements. Second, with the exception of CACTAs, noncoding TEs are enriched within and flanking genes, where they display conserved distribution patterns, having the highest peak in the promoter region. Finally, our analysis of microarray data revealed that genes associated with Tc1/mariner and PIF/Harbinger noncoding elements have significantly higher expression levels than genes without class 2 TEs. In contrast, genes with CACTA elements have significantly lower expression than genes without class 2 TEs.ConclusionsWe have achieved the most comprehensive annotation of class 2 TEs to date in these four grass genomes. Comparative analysis of this robust dataset led to the identification of several previously unknown features of each superfamily related to copy number, element size, genomic distribution and correlation with the expression levels of nearby genes. These results highlight the importance of distinguishing TE superfamilies when assessing their impact on gene and genome evolution.