We performed genome-wide meta-analysis of lipid traits on three samples of Mexican and Mexican American ancestry comprising 4,383 individuals, and followed up significant and highly suggestive associations in three additional Hispanic samples comprising 7,876 individuals. Genome-wide significant signals were observed in or near CELSR2, ZNF259/APOA5, KANK2/DOCK6 and NCAN/MAU2 for total cholesterol, LPL, ABCA1, ZNF259/APOA5, LIPC and CETP for HDL cholesterol, CELSR2, APOB and NCAN/MAU2 for LDL cholesterol, and GCKR, TRIB1, ZNF259/APOA5 and NCAN/MAU2 for triglycerides. Linkage disequilibrium and conditional analyses indicate that signals observed at ABCA1 and LIPC for HDL cholesterol and NCAN/MAU2 for triglycerides are independent of previously reported lead SNP associations. Analyses of lead SNPs from the European Global Lipids Genetics Consortium (GLGC) dataset in our Hispanic samples show remarkable concordance of direction of effects as well as strong correlation in effect sizes. A meta-analysis of the European GLGC and our Hispanic datasets identified five novel regions reaching genome-wide significance: two for total cholesterol (FN1 and SAMM50), two for HDL cholesterol (LOC100996634 and COPB1) and one for LDL cholesterol (LINC00324/CTC1/PFAS). The top meta-analysis signals were found to be enriched for SNPs associated with gene expression in a tissue-specific fashion, suggesting an enrichment of tissue-specific function in lipid-associated loci.

The SORBS1 gene plays an important role in insulin signaling. We aimed to examine whether common single-nucleotide polymorphisms (SNPs) of SORBS1 are associated with prevalence and incidence of diabetes, age at onset of diabetes, and the related traits of glucose homeostasis. A total of 1135 siblings from 492 ethnic Chinese families were recruited at baseline, and 630 were followed up for 5.19 ± 0.96 years. Nine SNPs including rs7081076, rs2281939, rs3818540, rs2274490, rs61739184, rs726176, rs2296966, rs17849148, and rs3193970 were genotyped and examined. To deal with correlated data of subjects within the same families, the generalized estimating equations approach was applied throughout all association analyses. The GG genotype of rs2281939 was associated with a higher risk of diabetes at baseline, an earlier onset of diabetes, and higher steady-state plasma glucose levels in the modified insulin suppression test. The minor allele T of rs2296966 was associated with higher prevalence and incidence of diabetes, an earlier onset of diabetes, and higher 2-h glucose during oral glucose tolerance test. These two SNPs revealed independent associations with age of diabetes onset as well as risk of diabetes at baseline. These findings supported that SORBS1 gene participates in the pathogenesis of diabetes.

Viral infections associate with disease risk and select families of viruses encode miRNAs that control an efficient viral cycle. The association of viral miRNA expression with disease in a large human population has not been previously explored. We sequenced plasma RNA from 40 participants of the Framingham Heart Study (FHS, Offspring Cohort, Visit 8) and identified 3 viral miRNAs from 3 different human Herpesviridae. These miRNAs were mostly related to viral latency and have not been previously detected in human plasma. Viral miRNA expression was then screened in the plasma of 2763 participants of the remaining cohort utilizing high-throughput RT-qPCR. All 3 viral miRNAs associated with combinations of inflammatory or prothrombotic circulating biomarkers (sTNFRII, IL-6, sICAM1, OPG, P-selectin) but did not associate with hypertension, coronary heart disease or cancer. Using a large observational population, we demonstrate that the presence of select viral miRNAs in the human circulation associate with inflammatory biomarkers and possibly immune response, but fail to associate with overt disease. This study greatly extends smaller singular observations of viral miRNAs in the human circulation and suggests that select viral miRNAs, such as those for latency, may not impact disease manifestation.

The genetic basis of supraventricular and ventricular ectopy (SVE, VE) remains largely uncharacterized, despite established genetic mechanisms of arrhythmogenesis. To identify novel genetic variants associated with SVE/VE in ancestrally diverse human populations, we conducted a genome-wide association study of electrocardiographically identified SVE and VE in five cohorts including approximately 43,000 participants of African, European and Hispanic/Latino ancestry. In thirteen ancestry-stratified subgroups, we tested multivariable-adjusted associations of SVE and VE with single nucleotide polymorphism (SNP) dosage. We combined subgroup-specific association estimates in inverse variance-weighted, fixed-effects and Bayesian meta-analyses. We also combined fixed-effects meta-analytic t-test statistics for SVE and VE in multi-trait SNP association analyses. No loci reached genome-wide significance in trans-ethnic meta-analyses. However, we found genome-wide significant SNPs intronic to an apoptosis-enhancing gene previously associated with QRS interval duration (FAF1; lead SNP rs7545860; effect allele frequency = 0.02; P = 2.0 × 10-8) in multi-trait analysis among European ancestry participants and near a locus encoding calcium-dependent glycoproteins (DSC3; lead SNP rs8086068; effect allele frequency = 0.17) in meta-analysis of SVE (P = 4.0 × 10-8) and multi-trait analysis (P = 2.9 × 10-9) among African ancestry participants. The novel findings suggest several mechanisms by which genetic variation may predispose to ectopy in humans and highlight the potential value of leveraging pleiotropy in future studies of ectopy-related phenotypes.

Hypertension prevalence varies between ethnic groups, possibly due to differences in genetic, environmental, and cultural determinants. Hispanic/Latino Americans are a diverse and understudied population. We performed a genome-wide association study (GWAS) of blood pressure (BP) traits in 12,278 participants from the Hispanics Community Health Study/Study of Latinos (HCHS/SOL). In the discovery phase we identified eight previously unreported BP loci. In the replication stage, we tested these loci in the 1982 Pelotas Birth Cohort Study of admixed Southern Brazilians, the COGENT-BP study of African descent, women of European descent from the Women Health Initiative (WHI), and a sample of European descent from the UK Biobank. No loci met the Bonferroni-adjusted level of statistical significance (0.0024). Two loci had marginal evidence of replication: rs78701042 (NGF) with diastolic BP (P = 0.008 in the 1982 Pelotas Birth Cohort Study), and rs7315692 (SLC5A8) with systolic BP (P = 0.007 in European ancestry replication). We investigated whether previously reported loci associated with BP in studies of European, African, and Asian ancestry generalize to Hispanics/Latinos. Overall, 26% of the known associations in studies of individuals of European and Chinese ancestries generalized, while only a single association previously discovered in a people of African descent generalized.

QT interval prolongation is a heritable risk factor for ventricular arrhythmias and can predispose to sudden death. Most genome-wide association studies (GWAS) of QT were performed in European ancestral populations, leaving other groups uncharacterized. Herein we present the first QT GWAS of Hispanic/Latinos using data on 15,997 participants from four studies. Study-specific summary results of the association between 1000 Genomes Project (1000G) imputed SNPs and electrocardiographically measured QT were combined using fixed-effects meta-analysis. We identified 41 genome-wide significant SNPs that mapped to 13 previously identified QT loci. Conditional analyses distinguished six secondary signals at NOS1AP (n = 2), ATP1B1 (n = 2), SCN5A (n = 1), and KCNQ1 (n = 1). Comparison of linkage disequilibrium patterns between the 13 lead SNPs and six secondary signals with previously reported index SNPs in 1000G super populations suggested that the SCN5A and KCNE1 lead SNPs were potentially novel and population-specific. Finally, of the 42 suggestively associated loci, AJAP1 was suggestively associated with QT in a prior East Asian GWAS; in contrast BVES and CAP2 murine knockouts caused cardiac conduction defects. Our results indicate that whereas the same loci influence QT across populations, population-specific variation exists, motivating future trans-ethnic and ancestrally diverse QT GWAS.

It is unclear whether genetic markers interact with risk factors to influence atrial fibrillation (AF) risk. We performed genome-wide interaction analyses between genetic variants and age, sex, hypertension, and body mass index in the AFGen Consortium. Study-specific results were combined using meta-analysis (88,383 individuals of European descent, including 7,292 with AF). Variants with nominal interaction associations in the discovery analysis were tested for association in four independent studies (131,441 individuals, including 5,722 with AF). In the discovery analysis, the AF risk associated with the minor rs6817105 allele (at the PITX2 locus) was greater among subjects ≤ 65 years of age than among those > 65 years (interaction p-value = 4.0 × 10-5). The interaction p-value exceeded genome-wide significance in combined discovery and replication analyses (interaction p-value = 1.7 × 10-8). We observed one genome-wide significant interaction with body mass index and several suggestive interactions with age, sex, and body mass index in the discovery analysis. However, none was replicated in the independent sample. Our findings suggest that the pathogenesis of AF may differ according to age in individuals of European descent, but we did not observe evidence of statistically significant genetic interactions with sex, body mass index, or hypertension on AF risk.

Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to link SomaScan 1.3K proteomes to genomes for 2812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We correctly linked 90-95% of proteomes to their correct genome and for 95-99% we identify the 1% most likely links. The linking accuracy in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. We also linked proteomes-to-proteomes and used the proteome only to determine features such as sex, ancestry, and first-degree relatives. When serial proteomes are available, the linking algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

Five sequence variants in SLC16A11 (rs117767867, rs13342692, rs13342232, rs75418188, and rs75493593), which occur in two non-reference haplotypes, were recently shown to be associated with diabetes in Mexicans from the SIGMA consortium. We aimed to determine whether these previous findings would replicate in the HCHS/SOL Mexican origin group and whether genotypic effects were similar in other HCHS/SOL groups. We analyzed these five variants in 2492 diabetes cases and 5236 controls from the Hispanic Community Health Study/Study of Latinos (HCHS/SOL), which includes U.S. participants from six diverse background groups (Mainland groups: Mexican, Central American, and South American; and Caribbean groups: Puerto Rican, Cuban, and Dominican). We estimated the SNP-diabetes association in the six groups and in the combined sample. We found that the risk alleles occur in two non-reference haplotypes in HCHS/SOL, as in the SIGMA Mexicans. The haplotype frequencies were very similar between SIGMA Mexicans and the HCHS/SOL Mainland groups, but different in the Caribbean groups. The SLC16A11 sequence variants were significantly associated with risk for diabetes in the Mexican origin group (P = 0.025), replicating the SIGMA findings. However, these variants were not significantly associated with diabetes in a combined analysis of all groups, although the power to detect such effects was 85% (assuming homogeneity of effects among the groups). Additional analyses performed separately in each of the five non-Mexican origin groups were not significant. We also analyzed (1) exclusion of young controls and, (2) SNP by BMI interactions, but neither was significant in the HCHS/SOL data. The previously reported effects of SLC16A11 variants on diabetes in Mexican samples was replicated in a large Mexican-American sample, but these effects were not significant in five non-Mexican Hispanic/Latino groups sampled from U.S. populations. Lack of replication in the HCHS/SOL non-Mexicans, and in the entire HCHS/SOL sample combined may represent underlying genetic heterogeneity. These results indicate a need for future genetic research to consider heterogeneity of the Hispanic/Latino population in the assessment of disease risk, but add to the evidence suggesting SLC16A11 as a potential therapeutic target for type 2 diabetes.

Previous research has shown that genes play a substantial role in determining a person's susceptibility to age-related hearing impairment. The existing studies on this subject have different results, which may be caused by difficulties in determining the phenotype or the limited number of participants involved. Here, we have gathered the largest sample to date (discovery n = 9,675; replication n = 10,963; validation n = 356,141), and examined phenotypes that represented low/mid and high frequency hearing loss on the pure tone audiogram. We identified 7 loci that were either replicated and/or validated, of which 5 loci are novel in hearing. Especially the ILDR1 gene is a high profile candidate, as it contains our top SNP, is a known hearing loss gene, has been linked to age-related hearing impairment before, and in addition is preferentially expressed within hair cells of the inner ear. By verifying all previously published SNPs, we can present a paper that combines all new and existing findings to date, giving a complete overview of the genetic architecture of age-related hearing impairment. This is of importance as age-related hearing impairment is highly prevalent in our ageing society and represents a large socio-economic burden.

Prostate cancer (PCa) is a major cause of cancer death among men. The histopathological examination of post-surgical prostate specimens and manual annotation of PCa not only allow for detailed assessment of disease characteristics and extent, but also supply the ground truth for developing of computer-aided diagnosis (CAD) systems for PCa detection before definitive treatment. As manual cancer annotation is tedious and subjective, there have been a number of publications describing methods for automating the procedure via the analysis of digitized whole-slide images (WSIs). However, these studies have focused only on the analysis of WSIs stained with hematoxylin and eosin (H&E), even though there is additional information that could be obtained from immunohistochemical (IHC) staining. In this work, we propose a framework for automating the annotation of PCa that is based on automated colorimetric analysis of both H&E and IHC WSIs stained with a triple-antibody cocktail against high-molecular weight cytokeratin (HMWCK), p63, and α-methylacyl CoA racemase (AMACR). The analysis outputs were then used to train a regression model to estimate the distribution of cancerous epithelium within slides. The approach yielded an AUC of 0.951, sensitivity of 87.1%, and specificity of 90.7% as compared to slide-level annotations, and generalized well to cancers of all grades.

Rheumatoid arthritis (RA) is a risk factor for atherosclerotic cardiovascular diseases (CVD). Given the critical roles of the immune system and inflammatory signals in the pathogenesis of CVD, we hypothesized that interrogation of CVD-related proteins using integrative genomics might provide new insights into the pathophysiology of RA. We utilized two-sample Mendelian randomization (MR) for causal inference between circulating protein levels and RA by incorporating genetic variants, followed by colocalization to characterize the causal associations. Genetic variants from three sources were obtained: those associated with 71 CVD-related proteins measured in nearly 7000 Framingham Heart Study participants, a published genome-wide association study (GWAS) of RA (19 234 cases, 61 565 controls), and GWAS of rheumatoid factor (RF) levels from the UK Biobank (n = 30 565). We identified the soluble receptor for advanced glycation end products (sRAGE), a critical inflammatory pathway protein, as putatively causal and protective for both RA (odds ratio per 1-standard deviation increment in inverse-rank normalized sRAGE level = 0.364; 95% confidence interval 0.342-0.385; P = 6.40 × 10-241) and RF levels (β [change in RF level per sRAGE increment] = - 1.318; SE = 0.434; P = 0.002). Using an integrative genomic approach, we highlight the AGER/RAGE axis as a putatively causal and promising therapeutic target for RA.

Despite the prognostic value of arterial stiffness (AS) and pulsatile hemodynamics (PH) for cardiovascular morbidity and mortality, epigenetic modifications that contribute to AS/PH remain unknown. To gain a better understanding of the link between epigenetics (DNA methylation) and AS/PH, we examined the relationship of eight measures of AS/PH with CpG sites and co-methylated regions using multi-ancestry participants from Trans-Omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA) with sample sizes ranging from 438 to 874. Epigenome-wide association analysis identified one genome-wide significant CpG (cg20711926-CYP1B1) associated with aortic augmentation index (AIx). Follow-up analyses, including gene set enrichment analysis, expression quantitative trait methylation analysis, and functional enrichment analysis on differentially methylated positions and regions, further prioritized three CpGs and their annotated genes (cg23800023-ETS1, cg08426368-TGFB3, and cg17350632-HLA-DPB1) for AIx. Among these, ETS1 and TGFB3 have been previously prioritized as candidate genes. Furthermore, both ETS1 and HLA-DPB1 have significant tissue correlations between Whole Blood and Aorta in GTEx, which suggests ETS1 and HLA-DPB1 could be potential biomarkers in understanding pathophysiology of AS/PH. Overall, our findings support the possible role of epigenetic regulation via DNA methylation of specific genes associated with AIx as well as identifying potential targets for regulation of AS/PH.

Atrial fibrillation (AF) is a heritable disease that affects more than thirty million individuals worldwide. Extensive efforts have been devoted to the study of genetic determinants of AF. The objective of our study is to examine the effect of gene-gene interaction on AF susceptibility. We performed a large-scale association analysis of gene-gene interactions with AF in 8,173 AF cases, and 65,237 AF-free referents collected from 15 studies for discovery. We examined putative interactions between genome-wide SNPs and 17 known AF-related SNPs. The top interactions were then tested for association in an independent cohort for replication, which included more than 2,363 AF cases and 114,746 AF-free referents. One interaction, between rs7164883 at the HCN4 locus and rs4980345 at the SLC28A1 locus, was found to be significantly associated with AF in the discovery cohorts (interaction OR = 1.44, 95% CI: 1.27-1.65, P = 4.3 × 10-8). Eight additional gene-gene interactions were also marginally significant (P < 5 × 10-7). However, none of the top interactions were replicated. In summary, we did not find significant interactions that were associated with AF susceptibility. Future increases in sample size and denser genotyping might facilitate the identification of gene-gene interactions associated with AF.

Sixty genetic loci associated with abdominal obesity, measured by waist circumference (WC) and waist-hip ratio (WHR), have been previously identified, primarily from studies conducted in European-ancestry populations. We conducted a meta-analysis of associations of abdominal obesity with approximately 2.5 million single nucleotide polymorphisms (SNPs) among 53,052 (for WC) and 48,312 (for WHR) individuals of Asian descent, and replicated 33 selected SNPs among 3,762 to 17,110 additional individuals. We identified four novel loci near the EFEMP1, ADAMTSL3 , CNPY2, and GNAS genes that were associated with WC after adjustment for body mass index (BMI); two loci near the NID2 and HLA-DRB5 genes associated with WHR after adjustment for BMI, and three loci near the CEP120, TSC22D2, and SLC22A2 genes associated with WC without adjustment for BMI. Functional enrichment analyses revealed enrichment of corticotropin-releasing hormone signaling, GNRH signaling, and/or CDK5 signaling pathways for those newly-identified loci. Our study provides additional insight on genetic contribution to abdominal obesity.

Genome-wide association studies have identified numerous variants associated with lipid levels; yet, the majority are located in non-coding regions with unclear mechanisms. In the Insulin Resistance Atherosclerosis Family Study (IRASFS), heritability estimates suggest a strong genetic basis: low-density lipoprotein (LDL, h2 = 0.50), high-density lipoprotein (HDL, h2 = 0.57), total cholesterol (TC, h2 = 0.53), and triglyceride (TG, h2 = 0.42) levels. Exome sequencing of 1,205 Mexican Americans (90 pedigrees) from the IRASFS identified 548,889 variants and association and linkage analyses with lipid levels were performed. One genome-wide significant signal was detected in APOA5 with TG (rs651821, PTG = 3.67 × 10-10, LODTG = 2.36, MAF = 14.2%). In addition, two correlated SNPs (r2 = 1.0) rs189547099 (PTG = 6.31 × 10-08, LODTG = 3.13, MAF = 0.50%) and chr4:157997598 (PTG = 6.31 × 10-08, LODTG = 3.13, MAF = 0.50%) reached exome-wide significance (P < 9.11 × 10-08). rs189547099 is an intronic SNP in FNIP2 and SNP chr4:157997598 is intronic in GLRB. Linkage analysis revealed 46 SNPs with a LOD > 3 with the strongest signal at rs1141070 (LODLDL = 4.30, PLDL = 0.33, MAF = 21.6%) in DFFB. A total of 53 nominally associated variants (P < 5.00 × 10-05, MAF ≥ 1.0%) were selected for replication in six Mexican-American cohorts (N = 3,280). The strongest signal observed was a synonymous variant (rs1160983, PLDL = 4.44 × 10-17, MAF = 2.7%) in TOMM40. Beyond primary findings, previously reported lipid loci were fine-mapped using exome sequencing in IRASFS. These results support that exome sequencing complements and extends insights into the genetics of lipid levels.

Targeted delivery of lysosomal enzymes to the endocytic compartment of human cells represents a transformative technology for treating a large family of lysosomal storage diseases (LSDs). Gaucher disease is one of the most common types of LSDs caused by mutations to the lysosomal β-glucocerebrosidase (GBA). Here, we describe a genetic strategy to produce engineered exosomes loaded with GBA in two different spatial configurations for targeted delivery to the endocytic compartment of recipient cells. By fusing human GBA to an exosome-anchoring protein: vesicular stomatitis virus glycoprotein (VSVG), we demonstrate that the chimeric proteins were successfully integrated into exosomes which were secreted as extracellular vesicles (EVs) by producer cells. Isolation and molecular characterization of EVs confirmed that the fusion proteins were loaded onto exosomes without altering their surface markers, particle size or distribution. Further, enzyme-loaded exosomes/EVs added to cultured medium were taken up by recipient cells. Further, the endocytosed exosomes/EVs targeted to endocytic compartments exhibited a significant increase in GBA activity. Together, we have developed a novel method for targeting and delivery of lysosomal enzymes to their natural location: the endocytic compartment of recipient cells. Since exosomes/EVs have an intrinsic ability to cross the blood-brain-barrier, our technology may provide a new approach to treat severe types of LSDs, including Gaucher disease with neurological complications.

Gene function can be described using various measures. We integrated association studies of three types of omics data to provide insights into the pathophysiology of subclinical coronary disease and myocardial infarction (MI). Using multivariable regression models, we associated: (1) single nucleotide polymorphism, (2) DNA methylation, and (3) gene expression with coronary artery calcification (CAC) scores and MI. Among 3106 participants of the Framingham Heart Study, 65 (2.1%) had prevalent MI and 60 (1.9%) had incident MI, median CAC value was 67.8 [IQR 10.8, 274.9], and 1403 (45.2%) had CAC scores > 0 (prevalent CAC). Prevalent CAC was associated with AHRR (linked to smoking) and EXOC3 (affecting platelet function and promoting hemostasis). CAC score was associated with VWA1 (extracellular matrix protein associated with cartilage structure in endomysium). For prevalent MI we identified FYTTD1 (down-regulated in familial hypercholesterolemia) and PINK1 (linked to cardiac tissue homeostasis and ischemia-reperfusion injury). Incident MI was associated with IRX3 (enhancing browning of white adipose tissue) and STXBP3 (controlling trafficking of glucose transporter type 4 to plasma). Using an integrative trans-omics approach, we identified both putatively novel and known candidate genes associated with CAC and MI. Replication of findings is warranted.

Expression quantitative trait methylation (eQTM) analysis identifies DNA CpG sites at which methylation is associated with gene expression. The present study describes an eQTM resource of CpG-transcript pairs derived from whole blood DNA methylation and RNA sequencing gene expression data in 2115 Framingham Heart Study participants. We identified 70,047 significant cis CpG-transcript pairs at p < 1E-7 where the top most significant eGenes (i.e., gene transcripts associated with a CpG) were enriched in biological pathways related to cell signaling, and for 1208 clinical traits (enrichment false discovery rate [FDR] ≤ 0.05). We also identified 246,667 significant trans CpG-transcript pairs at p < 1E-14 where the top most significant eGenes were enriched in biological pathways related to activation of the immune response, and for 1191 clinical traits (enrichment FDR ≤ 0.05). Independent and external replication of the top 1000 significant cis and trans CpG-transcript pairs was completed in the Women's Health Initiative and Jackson Heart Study cohorts. Using significant cis CpG-transcript pairs, we identified significant mediation of the association between CpG sites and cardiometabolic traits through gene expression and identified shared genetic regulation between CpGs and transcripts associated with cardiometabolic traits. In conclusion, we developed a robust and powerful resource of whole blood eQTM CpG-transcript pairs that can help inform future functional studies that seek to understand the molecular basis of disease.

Atrial fibrillation (AF) is the most common cardiac arrhythmia, but little is known about the molecular mechanisms associated with AF arrhythmogenesis. DNA methylation is an important epigenetic mechanism that regulates gene expression and downstream biological processes. We hypothesize that DNA methylation might play an important role in the susceptibility to develop AF. A total of 2,639 participants from the Offspring Cohort of Framingham Heart Study were enrolled in the current study. These participants included 183 participants with prevalent AF and 220 with incident AF during up to 9 years follow up. Genome-wide methylation was profiled using the Illumina Infinium HumanMethylation450 BeadChip on blood-derived DNA collected during the eighth examination cycle (2005-2008). Two CpG sites were significantly associated with prevalent AF, and five CpGs were associated with incident AF after correction for multiple testing (FDR < 0.05). Fourteen previously reported genome-wide significant AF-related SNP were each associated with at least one CpG site; the most significant association was rs6490029 at the CUX2 locus and cg10833066 (P = 9.5 × 10-279). In summary, we performed genome-wide methylation profiling in a community-based cohort and identified seven methylation signatures associated with AF. Our study suggests that DNA methylation might play an important role in AF arrhythmogenesis.