Novel loci associated with PR interval in a genome-wide association study of 10 African American cohorts

Advancing drug development for atrial fibrillation by prioritising findings from human genetic association studies

BACKGROUND

Drug development for atrial fibrillation (AF) has failed to yield new approved compounds. We sought to identify and prioritise potential druggable targets with support from human genetics, by integrating the available evidence with bioinformatics sources relevant for AF drug development.

METHODS

Genetic hits for AF and related traits were identified through structured search of MEDLINE. Genes derived from each paper were cross-referenced with the OpenTargets platform for drug interactions. Confirmation/validation was demonstrated through structured searches and review of evidence on MEDLINE and ClinialTrials.gov for each drug and its association with AF.

FINDINGS

613 unique drugs were identified, with 21 already included in AF Guidelines. Cardiovascular drugs from classes not currently used for AF (e.g. ranolazine and carperitide) and anti-inflammatory drugs (e.g. dexamethasone and mehylprednisolone) had evidence of potential benefit. Further targets were considered druggable but remain open for drug development.

INTERPRETATION

Our systematic approach, combining evidence from different bioinformatics platforms, identified drug repurposing opportunities and druggable targets for AF.

FUNDING

KK is supported by Barts Charity grant G-002089 and is mentored on the AFGen 2023-24 Fellowship funded by the AFGen NIH/NHLBI grant R01HL092577. RP is supported by the UCL BHF Research Accelerator AA/18/6/34223 and NIHR grant NIHR129463. AFS is supported by the BHF grants PG/18/5033837, PG/22/10989 and UCL BHF Accelerator AA/18/6/34223 as well as the UK Research and Innovation (UKRI) under the UK government's Horizon Europe funding guarantee EP/Z000211/1 and by the UKRI-NIHR grant MR/V033867/1 for the Multimorbidity Mechanism and Therapeutics Research Collaboration. AF is supported by UCL BHF Accelerator AA/18/6/34223. CF is supported by UCL BHF Accelerator AA/18/6/34223.

Differential genome-wide associated variants and enriched pathways of ECG parameters among people with versus without HIV

OBJECTIVES

People with HIV (PWH) are more likely to develop ECG abnormalities. Substantial evidence exists for genetic contribution to ECG parameters among general population. However, whether and how would host genome associate with ECG parameters among PWH is unclear. Our research aims to analyze and compare genetic variants, mapped genes, and enriched pathways of ECG parameters among PWH and HIV-negative controls.

DESIGN

A cross-sectional study.

METHOD

We performed a large original genome-wide association study (GWAS) of ECG parameters among PWH ( n  = 1730) and HIV-negative controls ( n =  3746). Genome-wide interaction analyses were also conducted.

RESULTS

A total of 18 novel variants were detected among PWH, six for PR interval including rs76345397 at ATL2 , 11 for QRS duration including rs10483994 at KCNK10 and rs2478830 at JCAD , and one for QTc interval (rs9815364). Among HIV-negative controls, we identified variants located at previously reported ECG-related genes ( SCN5A , CNOT1 ). Genetic variants had a significant interaction with HIV infection ( P  < 5 × 10 -8 ), implying that HIV infection and host genome might jointly influence ECG parameters. Mapped genes for PR interval and QRS duration among PWH were enriched in the biological process of viral genome replication and host response to virus, respectively, whereas enriched pathways for PR interval among HIV-negative controls were in the cellular component of voltage-gated sodium channel complex.

CONCLUSION

The present GWAS indicated a distinctive impact of host genome on quantitative ECG parameters among PWH. Different from HIV-negative controls, host genome might influence the cardiac electrical activity by interfering with HIV viral infection, production, and latency among PWH.

Complex interactions between p.His558Arg and linked variants in the sodium voltage-gated channel alpha subunit 5 (Na V 1.5)

Common genetic polymorphisms may modify the phenotypic outcome when co-occurring with a disease-causing variant, and therefore understanding their modulating role in health and disease is of great importance. The polymorphic p.His558Arg variant of the sodium voltage-gated channel alpha subunit 5 (Na _V 1.5) encoded by the SCN5A gene is a case in point, as several studies have shown it can modify the clinical phenotype in a number of cardiac diseases. To evaluate the genetic backgrounds associated with this modulating effect, we reanalysed previous electrophysiological findings regarding the p.His558Arg variant and further assessed its patterns of genetic diversity in human populations. The Na _V 1.5 p.His558Arg variant was found to be in linkage disequilibrium with six other polymorphic variants that previously were also associated with cardiac traits in GWAS analyses. On account of this, incongruent reports that Arg558 allele can compensate, aggravate or have no effect on Na _V 1.5, likely might have arose due to a role of p.His558Arg depending on the additional linked variants. Altogether, these results indicate a major influence of the epistatic interactions between SCN5A variants, revealing also that phenotypic severity may depend on the polymorphic background associated to each individual genome.

Atrial Remodeling in Atrial Fibrillation. Comorbidities and Markers of Disease Progression Predict Catheter Ablation Outcome

Atrial fibrillation is the most common supraventricular arrhythmia affecting an increasing proportion of the population in which mainstream therapy, i.e. catheter ablation, provides freedom from arrhythmia in only a limited number of patients. Understanding the mechanism is key in order to find more effective therapies and to improve patient selection. In this review, the structural and electrophysiological changes of the atrial musculature that constitute atrial remodeling in atrial fibrillaton and how risk factors and markers of disease progression can predict catheter ablation outcome will be discussed in detail.

A genome-wide association and polygenic risk score study on abnormal electrocardiogram in a Chinese population

Electrocardiography is a common and widely-performed medical examination based on the measurement and evaluation of electrocardiogram (ECG) to assess the up-to-date cardiac rhythms and thus suggest the health conditions of cardiovascular system and on a larger level the individual’s wellness. Abnormal ECG assessment from the detection of abnormal heart rhythms may have clinical implications including blood clots in formation, ongoing heart attack, coronary artery blockage, etc. Past genetic-phenotypic research focused primarily on the physical parameters of ECG but not the medical evaluation. To unbiasedly uncover the underlying links of genetic variants with normal vs. abnormal ECG assessment, a genome-wide association study (GWAS) is carried out in a 1006-participant cohort of Chinese population effectively genotyped for 243487 single nucleotide polymorphisms (SNPs). Both age and sex are influential factors, and six novel SNPs are identified for potential association with abnormal ECG. With the selected SNPs, a polygenic risk score (PRS) differentiates the case–control subgroups, and correlates well with increased risk of abnormal ECG. The findings are reproduced in an independent validation cohort. The derived PRS may function as a potential biomarker for prospectively screening the high-risk subgroup of heart issues in the Chinese population.

Pathophysiology and therapeutic relevance of PI3K(p110α) protein in atrial fibrillation: A non-interventional molecular therapy strategy

Genetically modified animal studies have revealed specific expression patterns and unequivocal roles of class I PI3K isoenzymes. PI3K(p110α), a catalytic subunit of class I PI3Ks is ubiquitously expressed and is well characterised in the cardiovascular system. Given that genetic inhibition of PI3K(p110α) causes lethal phenotype embryonically, the catalytic subunit is critically important in housekeeping and biological processes. A growing number of studies underpin crucial roles of PI3K(p110α) in cell survival, proliferation, hypertrophy and arrhythmogenesis. While the studies provide great insights, the precise mechanisms involved in PI3K(p110α) hypofunction and atrial fibrillation (AF) are not fully known. AF is a well recognised clinical problem with significant management limitations. In this translational review, we attempted a narration of PI3K(p110α) hypofunction in the molecular basis of AF pathophysiology. We sought to cautiously highlight the relevance of this molecule in the therapeutic approaches for AF management per se (i.e without conditions associate with cell proliferation, like cancer), and in mitigating effects of clinical risk factors in atrial substrate formation leading to AF progression. We also considered PI3K(p110α) in AF gene association, with the aim of identifying mechanistic links between the ever increasingly well-defined genetic loci (regions and genes) and AF. Such mechanisms will aid in identifying new drug targets for arrhythmogenic substrate and AF.

Multi-ancestry GWAS of the electrocardiographic PR interval identifies 202 loci underlying cardiac conduction

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N = 293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.

Allele-specific NKX2-5 binding underlies multiple genetic associations with human electrocardiographic traits

The cardiac transcription factor (TF) gene NKX2-5 has been associated with electrocardiographic (EKG) traits through genome-wide association studies (GWASs), but the extent to which differential binding of NKX2-5 at common regulatory variants contributes to these traits has not yet been studied. We analyzed transcriptomic and epigenomic data from induced pluripotent stem cell-derived cardiomyocytes from seven related individuals, and identified ~2,000 single-nucleotide variants associated with allele-specific effects (ASE-SNVs) on NKX2-5 binding. NKX2-5 ASE-SNVs were enriched for altered TF motifs, for heart-specific expression quantitative trait loci and for EKG GWAS signals. Using fine-mapping combined with epigenomic data from induced pluripotent stem cell-derived cardiomyocytes, we prioritized candidate causal variants for EKG traits, many of which were NKX2-5 ASE-SNVs. Experimentally characterizing two NKX2-5 ASE-SNVs (rs3807989 and rs590041) showed that they modulate the expression of target genes via differential protein binding in cardiac cells, indicating that they are functional variants underlying EKG GWAS signals. Our results show that differential NKX2-5 binding at numerous regulatory variants across the genome contributes to EKG phenotypes.

Genetic variants on chromosomes 7p31 and 12p12 are associated with abnormal atrial electrical activation in patients with early-onset lone atrial fibrillation

BACKGROUND

Abnormal P-wave morphology (PWM) has been associated with a history of atrial fibrillation (AF) in earlier studies. Although lone AF is believed to have substantial genetic basis, studies on associations between single nucleotide polymorphisms (SNP) linked to lone AF and PWM have not been reported. We aimed to assess whether SNPs previously associated with lone AF (rs2200733, rs13376333, rs3807989, and rs11047543) are also linked to P-wave abnormalities.

METHODS

Four SNPs were studied in 176 unrelated individuals with early-onset lone AF (age at onset <50 years), median age 38 years (19-63 years), 149 men. Using sinus rhythm ECG, orthogonal PWM was classified as Type 1-positive in leads X and Y and negative in lead Z, Type 2-positive in leads X and Y and biphasic (-/+) in lead Z, Type 3-positive in lead X and biphasic in lead Y (+/-), and the remaining as atypical.

RESULTS

Two SNPs were found to be significantly associated with altered P-wave morphology distribution: rs3807989 near the gene CAV1/CAV2 and rs11047543 near the gene SOX5. Both SNPs were associated with a higher risk of non-Type 1 P-wave morphology (rs3807989: OR = 4.8, 95% CI = 2.3-10.2, p < 0.001; rs11047543: OR = 4.7, 95% CI = 1.1-20.5, p = 0.04). No association was observed for rs2200733 and rs13376333.

CONCLUSION

In this study, the two variants rs3807989 and rs11047543, previously associated with PR interval and lone AF, were associated with altered P-wave morphology distribution in patients with early-onset lone AF. These findings suggest that common genetic variants may modify atrial conduction properties.

Functional characterization of SCN10A variants in several cases of sudden unexplained death

BACKGROUND

Multiple genome-wide association studies (GWAS) and targeted gene sequencing have identified common variants in SCN10A in cases of PR and QRS duration abnormalities, atrial fibrillation and Brugada syndrome. The New York City Office of Chief Medical Examiner has now also identified five SCN10A variants of uncertain significance in six separate cases within a cohort of 330 sudden unexplained death events. The gene product of SCN10A is the Nav1.8 sodium channel. The purpose of this study was to characterize effects of these variants on Nav1.8 channel function to provide better information for the reclassification of these variants.

METHODS AND RESULTS

Patch clamp studies were performed to assess effects of the variants on whole-cell Nav1.8 currents. We also performed RNA-seq analysis and immunofluorescence confocal microcopy to determine Nav1.8 expression in heart. We show that four of the five rare 'variants of unknown significance' (L388M, L867F, P1102S and V1518I) are associated with altered functional phenotypes. The R756W variant behaved similar to wild-type under our experimental conditions. We failed to detect Nav1.8 protein expression in immunofluorescence microscopy in rat heart. Furthermore, RNA-seq analysis failed to detect full-length SCN10A mRNA transcripts in human ventricle or mouse specialized cardiac conduction system, suggesting that the effect of Nav1.8 on cardiac function is likely to be extra-cardiac in origin.

CONCLUSIONS

We have demonstrated that four of five SCN10A variants of uncertain significance, identified in unexplained death, have deleterious effects on channel function. These data extend the genetic testing of SUD cases, but significantly more clinical evidence is needed to satisfy the criteria needed to associate these variants with the onset of SUD.

Role of MEIS1 in restless legs syndrome: From GWAS to functional studies in mice

MEIS1 is a transcription factor playing an important role in the development of several organs, including central and peripheral nervous systems. A genetic locus spanning the MEIS1 coding region has been associated with the risk of RLS in genome-wide association studies, with increasing evidence that MEIS1 is the causal RLS gene. The RLS-linked genetic signal has been mapped to an intronic regulatory element within MEIS1. This element plays a role in the ganglionic eminences of the developing forebrain, with the RLS risk allele related to a reduced activation of the enhancer. This suggests that the ganglionic eminences play an important role in the development of genetic susceptibility to RLS. In addition, rare variants within MEIS1 have been shown to contribute to the disease risk. These variants were identified first in RLS families and later found in further RLS cases by targeted sequencing. Some of these variants alone are sufficient to suppress MEIS1 function in neural development, providing further evidence of the importance of neurodevelopmental processes in the pathological mechanism of MEIS1 in RLS. Heterozygous Meis1 inactivation in mice causes hyperactivity at the onset of the inactive period, consistent with human RLS. In addition, these mice revealed an effect of MEIS1 on the dopaminergic system at both the spinal and supraspinal level. More studies are needed in human genetics to determine the exact role of MEIS1 variants in the risk of RLS, as well as in functional genetics and animal studies to further elucidate the pathological mechanism of MEIS1 in RLS.

SLMAP3 isoform modulates cardiac gene expression and function

The sarcolemmal membrane associated proteins (SLMAPs) belong to the super family of tail anchored membrane proteins which serve diverse roles in biology including cell growth, protein trafficking and ion channel regulation. Mutations in human SLMAP have been linked to Brugada syndrome with putative deficits in trafficking of the sodium channel (Na_v1.5) to the cell membrane resulting in aberrant electrical activity and heart function. Three main SLMAP isoforms (SLMAP1 (35 kDa), SLMAP2 (45 kDa), and SLMAP3 (91 kDa)) are expressed in myocardium but their precise role remains to be defined. Here we generated transgenic (Tg) mice with cardiac-specific expression of the SLMAP3 isoform during postnatal development which present with a significant decrease (20%) in fractional shortening and (11%) in cardiac output at 5 weeks of age. There was a lack of any notable cardiac remodeling (hypertrophy, fibrosis or fetal gene activation) in Tg hearts but the electrocardiogram indicated a significant increase (14%) in the PR interval and a decrease (43%) in the R amplitude. Western blot analysis indicated a selective and significant decrease (55%) in protein levels of Na_v1.5 while 45% drop in its transcript levels were detectable by qRT-PCR. Significant decreases in the protein and transcript levels of the calcium transport system of the sarcoplasmic reticulum (SERCA2a/PLN) were also evident in Tg hearts. These data reveal a novel role for SLMAP3 in the selective regulation of important ion transport proteins at the level of gene expression and suggest that it may be a unique target in cardiovascular function and disease.

PR interval genome-wide association meta-analysis identifies 50 loci associated with atrial and atrioventricular electrical activity

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.

Genomic approaches for the elucidation of genes and gene networks underlying cardiovascular traits

Genome-wide association studies have shed light on the association between natural genetic variation and cardiovascular traits. However, linking a cardiovascular trait associated locus to a candidate gene or set of candidate genes for prioritization for follow-up mechanistic studies is all but straightforward. Genomic technologies based on next-generation sequencing technology nowadays offer multiple opportunities to dissect gene regulatory networks underlying genetic cardiovascular trait associations, thereby aiding in the identification of candidate genes at unprecedented scale. RNA sequencing in particular becomes a powerful tool when combined with genotyping to identify loci that modulate transcript abundance, known as expression quantitative trait loci (eQTL), or loci modulating transcript splicing known as splicing quantitative trait loci (sQTL). Additionally, the allele-specific resolution of RNA-sequencing technology enables estimation of allelic imbalance, a state where the two alleles of a gene are expressed at a ratio differing from the expected 1:1 ratio. When multiple high-throughput approaches are combined with deep phenotyping in a single study, a comprehensive elucidation of the relationship between genotype and phenotype comes into view, an approach known as systems genetics. In this review, we cover key applications of systems genetics in the broad cardiovascular field.

Genome-wide association study of PR interval in Hispanics/Latinos identifies novel locus at ID2

OBJECTIVE

PR interval (PR) is a heritable electrocardiographic measure of atrial and atrioventricular nodal conduction. Changes in PR duration may be associated with atrial fibrillation, heart failure and all-cause mortality. Hispanic/Latino populations have high burdens of cardiovascular morbidity and mortality, are highly admixed and represent exceptional opportunities for novel locus identification. However, they remain chronically understudied. We present the first genome-wide association study (GWAS) of PR in 14 756 participants of Hispanic/Latino ancestry from three studies.

METHODS

Study-specific summary results of the association between 1000 Genomes Phase 1 imputed single-nucleotide polymorphisms (SNPs) and PR assumed an additive genetic model and were adjusted for global ancestry, study centre/region and clinical covariates. Results were combined using fixed-effects, inverse variance weighted meta-analysis. Sequential conditional analyses were used to identify independent signals. Replication of novel loci was performed in populations of Asian, African and European descent. ENCODE and RoadMap data were used to annotate results.

RESULTS

We identified a novel genome-wide association (P<5×10-8) with PR at ID2 (rs6730558), which replicated in Asian and European populations (P<0.017). Additionally, we generalised 10 previously identified PR loci to Hispanics/Latinos. Bioinformatics annotation provided evidence for regulatory function in cardiac tissue. Further, for six loci that generalised, the Hispanic/Latino index SNP was genome-wide significant and identical to (or in high linkage disequilibrium with) the previously identified GWAS lead SNP.

CONCLUSIONS

Our results suggest that genetic determinants of PR are consistent across race/ethnicity, but extending studies to admixed populations can identify novel associations, underscoring the importance of conducting genetic studies in diverse populations.

Common and Rare Coding Genetic Variation Underlying the Electrocardiographic PR Interval

BACKGROUND

Electrical conduction from the cardiac sinoatrial node to the ventricles is critical for normal heart function. Genome-wide association studies have identified more than a dozen common genetic loci that are associated with PR interval. However, it is unclear whether rare and low-frequency variants also contribute to PR interval heritability.

METHODS

We performed large-scale meta-analyses of the PR interval that included 83 367 participants of European ancestry and 9436 of African ancestry. We examined both common and rare variants associated with the PR interval.

RESULTS

We identified 31 genetic loci that were significantly associated with PR interval after Bonferroni correction (P<1.2×10-6), including 11 novel loci that have not been reported previously. Many of these loci are involved in heart morphogenesis. In gene-based analysis, we found that multiple rare variants at MYH6 (P=5.9×10-11) and SCN5A (P=1.1×10-7) were associated with PR interval. SCN5A locus also was implicated in the common variant analysis, whereas MYH6 was a novel locus.

CONCLUSIONS

We identified common variants at 11 novel loci and rare variants within 2 gene regions that were significantly associated with PR interval. Our findings provide novel insights to the current understanding of atrioventricular conduction, which is critical for cardiac activity and an important determinant of health.

Fifteen Genetic Loci Associated With the Electrocardiographic P Wave

BACKGROUND

The P wave on an ECG is a measure of atrial electric function, and its characteristics may serve as predictors for atrial arrhythmias. Increased mean P-wave duration and P-wave terminal force traditionally have been used as markers for left atrial enlargement, and both have been associated with increased risk of atrial fibrillation. Here, we explore the genetic basis of P-wave morphology through meta-analysis of genome-wide association study results for P-wave duration and P-wave terminal force from 12 cohort studies.

METHODS AND RESULTS

We included 44 456 individuals, of which 6778 (16%) were of African ancestry. Genotyping, imputation, and genome-wide association study were performed at each study site. Summary-level results were meta-analyzed centrally using inverse-variance weighting. In meta-analyses of P-wave duration, we identified 6 significant (P<5×10^-8) novel loci and replicated a prior association with SCN10A. We identified 3 loci at SCN5A, TBX5, and CAV1/CAV2 that were jointly associated with the PR interval, PR segment, and P-wave duration. We identified 6 novel loci in meta-analysis of P-wave terminal force. Four of the identified genetic loci were significantly associated with gene expression in 329 left atrial samples. Finally, we observed that some of the loci associated with the P wave were linked to overall atrial conduction, whereas others identified distinct phases of atrial conduction.

CONCLUSIONS

We have identified 6 novel genetic loci associated with P-wave duration and 6 novel loci associated with P-wave terminal force. Future studies of these loci may aid in identifying new targets for drugs that may modify atrial conduction or treat atrial arrhythmias.

PR Interval Associated Genes, Atrial Remodeling and Rhythm Outcome of Catheter Ablation of Atrial Fibrillation-A Gene-Based Analysis of GWAS Data

Background: PR interval prolongation has recently been shown to associate with advanced left atrial remodeling and atrial fibrillation (AF) recurrence after catheter ablation. While different genome-wide association studies (GWAS) have implicated 13 loci to associate with the PR interval as an AF endophenotype their subsequent associations with AF remodeling and response to catheter ablation are unknown. Here, we perform a gene-based analysis of GWAS data to test the hypothesis that PR interval candidate genes also associate with left atrial remodeling and arrhythmia recurrence following AF catheter ablation.

Methods and Results: Samples from 660 patients with paroxysmal (n = 370) or persistent AF (n = 290) undergoing AF catheter ablation were genotyped for ~1,000,000 SNPs. Gene-based association was investigated using VEGAS (versatile gene-based association study). Among the 13 candidate genes, SLC8A1, MEIS1, ITGA9, SCN5A, and SOX5 associated with the PR interval. Of those, ITGA9 and SOX5 were significantly associated with left atrial low voltage areas and left atrial diameter and subsequently with AF recurrence after radiofrequency catheter ablation.

Conclusion: This study suggests contributions of ITGA9 and SOX5 to AF remodeling expressed as PR interval prolongation, low voltage areas and left atrial dilatation and subsequently to response to catheter ablation. Future and larger studies are necessary to replicate and apply these findings with the aim of designing AF pathophysiology-based multi-locus risk scores.

GAREM1 regulates the PR interval on electrocardiograms

PR interval is the period from the onset of P wave to the start of the QRS complex on electrocardiograms. A recent genomewide association study (GWAS) suggested that GAREM1 was linked to the PR interval on electrocardiograms. This study was designed to validate this correlation using additional subjects and examined the function of Garem1 in a mouse model. We analyzed the association of rs17744182, a variant in the GAREM1 locus, with the PR interval in 5646 subjects who were recruited from 2 Korean replication sets, Yangpyeong (n = 2471) and Yonsei (n = 3175), and noted a significant genomewide association by meta-analysis (P = 2.39 × 10^-8). To confirm the function of Garem1 in mice, Garem1 siRNA was injected into mouse tail veins to reduce the expression of Garem1. Garem1 transcript levels declined by 53% in the atrium of the heart (P = 0.029), and Garem1-siRNA injected mice experienced a significant decrease in PR interval (43.27 ms vs. 44.89 ms in control, P = 0.007). We analyzed the expression pattern of Garem1 in the heart by immunohistology and observed specific expression of Garem1 in intracardiac ganglia. Garem1 was expressed in most neurons of the ganglion, including cholinergic and adrenergic cells. We have provided evidence that GAREM1 is involved in the PR interval of ECGs. These findings increase our understanding of the regulatory signals of heart rhythm through intracardiac ganglia of the autonomic nervous system and can be used to guide the development of a therapeutic target for heart conditions, such as atrial fibrillation.

Meis1: effects on motor phenotypes and the sensorimotor system in mice

MEIS1 encodes a developmental transcription factor and has been linked to restless legs syndrome (RLS) in genome-wide association studies. RLS is a movement disorder leading to severe sleep reduction and has a substantial impact on the quality of life of patients. In genome-wide association studies, MEIS1 has consistently been the gene with the highest effect size and functional studies suggest a disease-relevant downregulation. Therefore, haploinsufficiency of Meis1 could be the system with the most potential for modeling RLS in animals. We used heterozygous Meis1-knockout mice to study the effects of Meis1 haploinsufficiency on mouse behavioral and neurological phenotypes, and to relate the findings to human RLS. We exposed the Meis1-deficient mice to assays of motor, sensorimotor and cognitive ability, and assessed the effect of a dopaminergic receptor 2/3 agonist commonly used in the treatment of RLS. The mutant mice showed a pattern of circadian hyperactivity, which is compatible with human RLS. Moreover, we discovered a replicable prepulse inhibition (PPI) deficit in the Meis1-deficient animals. In addition, these mice were hyposensitive to the PPI-reducing effect of the dopaminergic receptor agonist, highlighting a role of Meis1 in the dopaminergic system. Other reported phenotypes include enhanced social recognition at an older age that was not related to alterations in adult olfactory bulb neurogenesis previously shown to be implicated in this behavior. In conclusion, the Meis1-deficient mice fulfill some of the hallmarks of an RLS animal model, and revealed the role of Meis1 in sensorimotor gating and in the dopaminergic systems modulating it.

Summary: Loss of Meis1 results in motor restlessness in mice, a phenotype resembling human restless legs syndrome, as well as altered sensorimotor gating and improved social discrimination memory.

Investigating the Genetic Architecture of the PR Interval Using Clinical Phenotypes

BACKGROUND

One potential use for the PR interval is as a biomarker of disease risk. We hypothesized that quantifying the shared genetic architectures of the PR interval and a set of clinical phenotypes would identify genetic mechanisms contributing to PR variability and identify diseases associated with a genetic predictor of PR variability.

METHODS AND RESULTS

We used ECG measurements from the ARIC study (Atherosclerosis Risk in Communities; n=6731 subjects) and 63 genetically modulated diseases from the eMERGE network (Electronic Medical Records and Genomics; n=12 978). We measured pairwise genetic correlations (rG) between PR phenotypes (PR interval, PR segment, P-wave duration) and each of the 63 phenotypes. The PR segment was genetically correlated with atrial fibrillation (rG=-0.88; P=0.0009). An analysis of metabolic phenotypes in ARIC also showed that the P wave was genetically correlated with waist circumference (rG=0.47; P=0.02). A genetically predicted PR interval phenotype based on 645 714 single-nucleotide polymorphisms was associated with atrial fibrillation (odds ratio=0.89 per SD change; 95% confidence interval, 0.83-0.95; P=0.0006). The differing pattern of associations among the PR phenotypes is consistent with analyses that show that the genetic correlation between the P wave and PR segment was not significantly different from 0 (rG=-0.03 [0.16]).

CONCLUSIONS

The genetic architecture of the PR interval comprises modulators of atrial fibrillation risk and obesity.

Genomic Contributors to Rhythm Outcome of Atrial Fibrillation Catheter Ablation - Pathway Enrichment Analysis of GWAS Data

BACKGROUND

Left atrial enlargement and persistent atrial fibrillation (AF) are well-known predictors for arrhythmia recurrence after AF catheter ablation (LRAF). In this study, by using pathway enrichment analysis of GWAS data, we tested the hypothesis that genetic pathways associated with these phenotypes are also associated with LRAF.

METHODS

Samples from 660 patients with paroxysmal (n = 370) or persistent AF (n = 290) undergoing de-novo AF catheter ablation were genotyped for ~1,000,000 SNPs. SNPs found to be significantly associated with left atrial diameter (LAD) or AF type were used for gene-based association tests in a systematic biological Knowledge-based mining system for Genome-wide Genetic studies (KGG). Associated genes were tested for pathway enrichment using WEB-based Gene SeT AnaLysis Toolkit (WebGestalt), the Gene Annotation Tool to Help Explain Relationships (GATHER) and the databases provided by Kyoto Encyclopedia of Genes and Genomes (KEGG). In a second step, the association of consistently enriched pathways and LRAF was tested.

RESULTS

By using sequential 7-day Holter ECGs, LRAF between 3 and 12 months was observed in 48% and was associated with LAD (B = 1.801, 95% CI 0.760-2.841, p = 1.0E-3) and persistent AF (OR = 2.1; 95% CI 1.567-2.931, p = 2.0E-6). WebGestalt (adj. p = 2.7E-22) and GATHER (adj. p = 5.2E-3) identified the calcium signaling pathway (hsa04020) as the only consistently enriched pathway for LAD, while the extracellular matrix (ECM) -receptor interaction pathway (hsa04512) was the only consistently enriched pathway for AF type (adj. p = 2.1E-15 in WebGestalt; adj. p = 9.3E-4 in GATHER). Both calcium signaling (adj. p = 2.2E-17 in WebGestalt; adj. p = 2.9E-2 in GATHER) and ECM-receptor interaction (adj. p = 1.2E-10 in WebGestalt; adj. p = 2.9E-2 in GATHER) were significantly associated with LRAF.

CONCLUSIONS

Calcium signaling and ECM-receptor interaction pathways are associated with LAD and AF type and, in turn, with LRAF. Future and larger studies are necessary to replicate and apply these findings.

Fine-mapping, novel loci identification, and SNP association transferability in a genome-wide association study of QRS duration in African Americans

The electrocardiographic QRS duration, a measure of ventricular depolarization and conduction, is associated with cardiovascular mortality. While single nucleotide polymorphisms (SNPs) associated with QRS duration have been identified at 22 loci in populations of European descent, the genetic architecture of QRS duration in non-European populations is largely unknown. We therefore performed a genome-wide association study (GWAS) meta-analysis of QRS duration in 13,031 African Americans from ten cohorts and a transethnic GWAS meta-analysis with additional results from populations of European descent. In the African American GWAS, a single genome-wide significant SNP association was identified (rs3922844, P = 4 × 10-14) in intron 16 of SCN5A, a voltage-gated cardiac sodium channel gene. The QRS-prolonging rs3922844 C allele was also associated with decreased SCN5A RNA expression in human atrial tissue (P = 1.1 × 10-4). High density genotyping revealed that the SCN5A association region in African Americans was confined to intron 16. Transethnic GWAS meta-analysis identified novel SNP associations on chromosome 18 in MYL12A (rs1662342, P = 4.9 × 10-8) and chromosome 1 near CD1E and SPTA1 (rs7547997, P = 7.9 × 10-9). The 22 QRS loci previously identified in populations of European descent were enriched for significant SNP associations with QRS duration in African Americans (P = 9.9 × 10-7), and index SNP associations in or near SCN5A, SCN10A, CDKN1A, NFIA, HAND1, TBX5 and SETBP1 replicated in African Americans. In summary, rs3922844 was associated with QRS duration and SCN5A expression, two novel QRS loci were identified using transethnic meta-analysis, and a significant proportion of QRS-SNP associations discovered in populations of European descent were transferable to African Americans when adequate power was achieved.

Gene regulatory networks in atrial fibrillation

Atrial fibrillation (AF) is the most frequent arrhythmogenic syndrome in humans. With an estimate incidence of 1%-2% in the general population, AF raises up to almost 10%-12% in 80+ years. Thus, AF represents nowadays a highly prevalent medical problem generating a large economic burden. At the electrophysiological level, distinct mechanisms have been elucidated. Yet, despite its prevalence, the genetic and molecular culprits of this pandemic cardiac electrophysiological abnormality have remained largely obscure. Molecular genetics of AF familiar cases have demonstrated that single nucleotide mutations in distinct genes encoding for ion channels underlie the onset of AF, albeit such alterations only explain a minor subset of patients with AF. In recent years, analyses by means of genome-wide association studies have unraveled a more complex picture of the etiology of AF, pointing out to distinct cardiac-enriched transcription factors, as well as to other regulatory genes. Furthermore a new layer of regulatory mechanisms have emerged, i.e., post-transcriptional regulation mediated by non-coding RNA, which have been demonstrated to exert pivotal roles in cardiac electrophysiology. In this manuscript, we aim to provide a comprehensive review of the genetic regulatory networks that if impaired exert electrophysiological abnormalities that contribute to the onset, and subsequently, on self-perpetuation of AF.

Common and rare variants in SCN10A modulate the risk of atrial fibrillation

BACKGROUND

Genome-wide association studies have shown that the common single nucleotide polymorphism rs6800541 located in SCN10A, encoding the voltage-gated Nav1.8 sodium channel, is associated with PR-interval prolongation and atrial fibrillation (AF). Single nucleotide polymorphism rs6800541 is in high linkage disequilibrium with the nonsynonymous variant in SCN10A, rs6795970 (V1073A, r(2)=0.933). We therefore sought to determine whether common and rare SCN10A variants are associated with early onset AF.

METHODS AND RESULTS

SCN10A was sequenced in 225 AF patients in whom there was no evidence of other cardiovascular disease or dysfunction (lone AF). In an association study of the rs6795970 single nucleotide polymorphism variant, we included 515 AF patients and 2 control cohorts of 730 individuals free of AF and 6161 randomly sampled individuals. Functional characterization of SCN10A variants was performed by whole-cell patch-clamping. In the lone AF cohort, 9 rare missense variants and 1 splice site donor variant were detected. Interestingly, AF patients were found to have higher G allele frequency of rs6795970, which encodes the alanine variant at position 1073 (described from here on as A1073, odds ratio =1.35 [1.16-1.54]; P=2.3×10(-5)). Both of the common variants, A1073 and P1092, induced a gain-of-channel function, whereas the rare missense variants, V94G and R1588Q, resulted in a loss-of-channel function.

CONCLUSIONS

The common variant A1073 is associated with increased susceptibility to AF. Both rare and common variants have effect on the function of the channel, indicating that these variants influence susceptibility to AF. Hence, our study suggests that SCN10A variations are involved in the genesis of AF.

Common and Rare Variants in SCN10A Modulate the Risk of Atrial Fibrillation

Background—

Genome-wide association studies have shown that the common single nucleotide polymorphism rs6800541 located in SCN10A , encoding the voltage-gated Na v 1.8 sodium channel, is associated with PR-interval prolongation and atrial fibrillation (AF). Single nucleotide polymorphism rs6800541 is in high linkage disequilibrium with the nonsynonymous variant in SCN10A , rs6795970 (V1073A, r 2 =0.933). We therefore sought to determine whether common and rare SCN10A variants are associated with early onset AF.

Methods and Results—

SCN10A was sequenced in 225 AF patients in whom there was no evidence of other cardiovascular disease or dysfunction (lone AF). In an association study of the rs6795970 single nucleotide polymorphism variant, we included 515 AF patients and 2 control cohorts of 730 individuals free of AF and 6161 randomly sampled individuals. Functional characterization of SCN10A variants was performed by whole-cell patch-clamping. In the lone AF cohort, 9 rare missense variants and 1 splice site donor variant were detected. Interestingly, AF patients were found to have higher G allele frequency of rs6795970, which encodes the alanine variant at position 1073 (described from here on as A1073, odds ratio =1.35 [1.16−1.54]; P =2.3×10 −5 ). Both of the common variants, A1073 and P1092, induced a gain-of-channel function, whereas the rare missense variants, V94G and R1588Q, resulted in a loss-of-channel function.

Conclusions—

The common variant A1073 is associated with increased susceptibility to AF. Both rare and common variants have effect on the function of the channel, indicating that these variants influence susceptibility to AF. Hence, our study suggests that SCN10A variations are involved in the genesis of AF.

Replication of SCN5A Associations with Electrocardio-graphic Traits in African Americans from Clinical and Epidemiologic Studies

The NAv1.5 sodium channel α subunit is the predominant α-subunit expressed in the heart and is associated with cardiac arrhythmias. We tested five previously identified SCN5A variants (rs7374138, rs7637849, rs7637849, rs7629265, and rs11129796) for an association with PR interval and QRS duration in two unique study populations: the Third National Health and Nutrition Examination Survey (NHANES III, n= 552) accessed by the Epidemiologic Architecture for Genes Linked to Environment (EAGLE) and a combined dataset (n= 455) from two biobanks linked to electronic medical records from Vanderbilt University (BioVU) and Northwestern University (NUgene) as part of the electronic Medical Records & Genomics (eMERGE) network. A meta-analysis including all three study populations (n~4,000) suggests that eight SCN5A associations were significant for both QRS duration and PR interval (p<5.0E-3) with little evidence for heterogeneity across the study populations. These results suggest that published SCN5A associations replicate across different study designs in a meta-analysis and represent an important first step in utility of multiple study designs for genetic studies and the identification/characterization of genetic variants associated with ECG traits in African-descent populations.

Genome-wide association studies in Africans and African Americans: expanding the framework of the genomics of human traits and disease

Genomic research is one of the tools for elucidating the pathogenesis of diseases of global health relevance and paving the research dimension to clinical and public health translation. Recent advances in genomic research and technologies have increased our understanding of human diseases, genes associated with these disorders, and the relevant mechanisms. Genome-wide association studies (GWAS) have proliferated since the first studies were published several years ago and have become an important tool in helping researchers comprehend human variation and the role genetic variants play in disease. However, the need to expand the diversity of populations in GWAS has become increasingly apparent as new knowledge is gained about genetic variation. Inclusion of diverse populations in genomic studies is critical to a more complete understanding of human variation and elucidation of the underpinnings of complex diseases. In this review, we summarize the available data on GWAS in recent African ancestry populations within the western hemisphere (i.e. African Americans and peoples of the Caribbean) and continental African populations. Furthermore, we highlight ways in which genomic studies in populations of recent African ancestry have led to advances in the areas of malaria, HIV, prostate cancer, and other diseases. Finally, we discuss the advantages of conducting GWAS in recent African ancestry populations in the context of addressing existing and emerging global health conditions.

Prolonged PR interval predicts clinical recurrence of atrial fibrillation after catheter ablation

BACKGROUND

A prolonged PR interval is known to be a poor prognostic factor in cardiovascular disease. The aim of this study was to investigate the association between PR interval and clinical outcome in patients undergoing radiofrequency catheter ablation (RFCA) of atrial fibrillation (AF).

METHODS AND RESULTS

We prospectively included 576 patients with AF (75.5% male, 57.8±11.6 years old, 68.8% paroxysmal AF) who underwent RFCA. We analyzed preprocedural sinus rhythm ECGs obtained in the absence of antiarrhythmic drug, and all enrolled patients were categorized into 4 groups based on the quartile values of the PR interval (166, 182, and 202 ms), and were analyzed according to the left atrium (LA) volume (CT; Computed tomography), LA voltage (NavX), and clinical outcome of AF ablation. Based on quartile value of PR interval, the highest quartile of PR interval (Q4; PR ≥202 ms) was oldest (P<0.001), and most likely to have persistent AF (P<0.001) and hypertension (P=0.013) compared with the other groups. However, there was no significant difference in LA conduction velocity and atrial effective refractory period. Q4 had the greatest LA dimension (P<0.001) and volume index (P<0.001), and lowest LA appendage-emptying velocity (P<0.032) and LA voltage (P<0.001) compared with the others. For 13.1±7.5 months, the classification based on the PR interval was a significant predictor of AF recurrence after RFCA of AF (HR=1.969, 95% CI 1.343 to 2.886, P=0.001).

CONCLUSIONS

The PR interval was closely associated with advanced LA remodeling due to AF, and had a noninvasive significant predictive value of clinical recurrence of AF after RFCA.

Genome-wide association study identifies three susceptibility loci for laryngeal squamous cell carcinoma in the Chinese population

To identify genetic markers for laryngeal squamous cell carcinoma (LSCC), we conducted a genome-wide association study (GWAS) on 993 individuals with LSCC (cases) and 1,995 cancer-free controls from Chinese populations. The most promising variants (association P < 1 × 10(-5)) were then replicated in 3 independent sets including 2,398 cases and 2,804 controls, among which we identified 3 new susceptibility loci at 11q12 (rs174549), 6p21 (rs2857595) and 12q24 (rs10492336). The minor alleles of each of these loci showed protective effects, with odds ratios (95% confidence intervals) of 0.73 (0.68-0.78; P = 1.00 × 10(-20)), 0.78 (0.72-0.84; P = 2.43 × 10(-15)) and 0.71 (0.65-0.77; P = 4.48 × 10(-14)), respectively. None of these variants showed an interaction with smoking or drinking. This is the first GWAS to our knowledge solely on LSCC, and the findings might advance understanding of the etiology of LSCC.

Genome-wide association study of electrocardiographic parameters identifies a new association for PR interval and confirms previously reported associations

Previous reports have described several associations of PR, QRS, QT and heart rate with genomic variations by genome-wide association studies (GWASs). In the present study, we examined the association of ∼2.5 million SNPs from 2994 Japanese healthy volunteers obtained from the JPDSC database with electrocardiographic parameters. We confirmed associations of PR interval, QRS duration and QT interval in individuals of Japanese ancestry with 11 of the 45 SNPs (6 of 20 for QT, 5 of 19 for PR and 0 of 6 for QRS) observed among individuals of European, African and Asian (Indian and Korean) ancestries. Those results indicate that many of the electrocardiographic associations with genes are shared by different ethnic groups including Japanese. Possible novel associations found in this study were validated by Korean data. As a result, we identified a novel association of SNP rs4952632[G] (maps near SLC8A1, sodium-calcium exchanger) (P = 7.595 × 10(-6)) with PR interval in Japanese individuals, and replication testing among Koreans confirmed the association of the same SNP with prolonged PR interval. Meta-analysis of the Japanese and Korean datasets demonstrated highly significant associations of SNP rs4952632[G] with a 2.325-ms (95% CI, 1.693-2.957 ms) longer PR interval per minor allele copy (P = 5.598 × 10(-13)). Cell-type-specific SLC8A1 knockout mice have demonstrated a regulatory role of sodium-calcium exchanger in automaticity and conduction in sinoatrial node, atrium and atrioventricular node. Our findings support a functional role of sodium-calcium exchanger in human atrial and atrioventricular nodal conduction as suggested by genetically modified mouse models.

Lack of association between rs3807989 in cav1 and atrial fibrillation

BACKGROUND

A novel gene Caveolin-1(CAV1) was identified to be susceptibility to PR interval and also associated with atrial fibrillation (AF) in two Genome wide associations studies (GWAS) studies in European ancestry. The purpose of this study was to determine the association of the SNPs in CAV1 gene of rs3807989 with AF in Chinese Han patients.

METHODS AND RESULTS

We attempted a replication in a cohort of 839 Chinese AF patients and 1215 healthy controls using melting temperature shift allele-specific genotyping analysis. One SNP in CAV1 (rs3807989) was genotyped. The final study cohort consisted of 839 AF patients and 1215 healthy controls. No significant association was detected between rs3807989 and AF in a Chinese Han population (allelic P-adj = 0.828 with OR = 1.02; genotypic P-adj = 0.815, 0.405, 0.760 with a dominant model, recessive model, and additive model). After logistic regression with multiple covariates, the association remained non-significant with adjusted P value 0.828. When the AF cases were further divided into lone AF (31.5%) and other types of AF (68.5%), no significant association was found between rs3807989 and lone AF (P-adj = 0.929 with OR = 0.990) and other types of AF (P-adj = 0.597 with OR = 1.060).

CONCLUSION

The SNP rs3807989 in CAV1 gene is not associated with AF or lone AF in our studies, which suggests that the SNP rs3807989 in CAV1 may not be a risk factor for AF in Chinese Han population.

Genetic determinants of P wave duration and PR segment

BACKGROUND

The PR interval on the ECG reflects atrial depolarization and atrioventricular nodal delay which can be partially differentiated by P wave duration and PR segment, respectively. Genome-wide association studies have identified several genetic loci for PR interval, but it remains to be determined whether this is driven by P wave duration, PR segment, or both.

METHODS AND RESULTS

We replicated 7 of the 9 known PR interval loci in 16 468 individuals of European ancestry. Four loci were unambiguously associated with PR segment, while the others were shared for P wave duration and PR segment. Next, we performed a genome-wide analysis on P wave duration and PR segment separately and identified 5 novel loci. Single-nucleotide polymorphisms in KCND3 (P=8.3×10(-11)) and FADS2 (P=2.7×10(-8)) were associated with P wave duration, whereas single-nucleotide polymorphisms near IL17D (P=2.3×10(-8)), in EFHA1 (P=3.3×10(-10)), and in LRCH1 (P=2.1×10(-8)) were associated with PR segment. Analysis on DNA elements indicated that genome-wide significant single-nucleotide polymorphisms were enriched at genomic regions suggesting active gene transcription in the human right atrium. Quantitative polymerase chain reaction showed that genes were significantly higher expressed in the right atrium and atrioventricular node compared with left ventricle (P=5.6×10(-6)).

CONCLUSIONS

Genetic associations of PR interval seem to be mainly driven by genetic determinants of the PR segment. Some of the PR interval associations are strengthened by a directional consistent effect of genetic determinants of P wave duration. Through genome-wide association we also identified genetic variants specifically associated with P wave duration which might be relevant for cardiac biology.

First-degree AV block-an entirely benign finding or a potentially curable cause of cardiac disease?

First-degree atrioventricular (AV) block is a delay within the AV conduction system and is defined as a prolongation of the PR interval beyond the upper limit of what is considered normal (generally 0.20 s). Up until recently, first-degree AV block was considered an entirely benign condition. In fact, some complain that it is a misnomer since there is only delay and no actual block in the AV conduction system (usually within the AV node). However, it has long been acknowledged that extreme forms of first-degree AV block (typically a PR interval exceeding 0.30 s) can cause symptoms due to inadequate timing of atrial and ventricular contractions, similar to the so-called pacemaker syndrome. Consequently, the current guidelines state that permanent pacemaker implantation is reasonable for first-degree AV block with symptoms similar to those of pacemaker syndrome or with hemodynamic compromise, but also stresses that there is little evidence to suggest that pacemakers improve survival in patients with isolated first-degree AV block. Recent reports suggest that it may be time to revisit the impact of first-degree AV block. Also, several findings in post hoc analyses of randomized device trials give important insights in possible treatment options. The present review aims to provide an update on the current knowledge concerning the impact of first-degree AV block and also to address the issue of pacing in patients with this condition.

Gene regulatory elements of the cardiac conduction system

The coordinated contraction of the heart relies on the generation and conduction of the electrical impulse. Aberrations of the function of the cardiac conduction system have been associated with various arrhythmogenic disorders and increased risk of sudden cardiac death. The genetics underlying conduction system function have been investigated using functional studies and genome-wide association studies. Both methods point towards the involvement of ion channel genes and the transcription factors that govern their activity. A large fraction of disease- and trait-associated sequence variants lie within non-coding sequences, enriched with epigenetic marks indicative of regulatory DNA. Although sequence conservation as a result of functional constraint has been a useful property to identify transcriptional enhancers, this identification process has been advanced through the development of techniques such as ChIP-seq and chromatin conformation capture technologies. The role of variation in gene regulatory elements in the cardiac conduction system has recently been demonstrated by studies on enhancers of SCN5A/SCN10A and TBX5. In both studies, a region harbouring a functionally implicated single-nucleotide polymorphism was shown to drive reproducible cardiac expression in a reporter gene assay. Furthermore, the risk variant of the allele abrogated enhancer function in both cases. Functional studies on regulatory DNA will likely receive a boost through recent developments in genome modification technologies.

Reference ranges of PR duration and P-wave indices in individuals free of cardiovascular disease: the Multi-Ethnic Study of Atherosclerosis (MESA)

In this brief report, we provide normal reference ranges for PR duration [unadjusted and heart rate adjusted] and P-wave indices [duration, amplitude and terminal force in V1] in individuals free of cardiovascular disease and its risk factors. We used automatically processed digital ECG data from 1252 US participants [mean age 59 (± 10) years, 738 women, 588 whites, 207 African-Americans, 217 Hispanics, 240 Chinese] from the Multi-Ethnic Study of Atherosclerosis [MESA]. In multivariable adjusted linear regression models with PR and each P-wave variable as a separate outcome, significant age, sex and race differences in these markers were observed. Subsequently, we report reference ranges for abnormal [2nd and 98th percentiles], borderline abnormal [5th and 95th percentiles] and mean [SD] values of PR and P-wave indices stratified by age [middle age (45-64 years) and seniors (65-84 years)], sex [men and women] and race [whites, African Americans, Hispanics and Chinese].