Heritabilities, proportions of heritabilities explained by GWAS findings, and implications of cross-phenotype effects on PR interval

Evolutionary Origins of Metabolic Reprogramming in Cancer

Metabolic changes that facilitate tumor growth are one of the hallmarks of cancer. These changes are not specific to tumors but also take place during the physiological growth of tissues. Indeed, the cellular and tissue mechanisms present in the tumor have their physiological counterpart in the repair of tissue lesions and wound healing. These molecular mechanisms have been acquired during metazoan evolution, first to eliminate the infection of the tissue injury, then to enter an effective regenerative phase. Cancer itself could be considered a phenomenon of antagonistic pleiotropy of the genes involved in effective tissue repair. Cancer and tissue repair are complex traits that share many intermediate phenotypes at the molecular, cellular, and tissue levels, and all of these are integrated within a Systems Biology structure. Complex traits are influenced by a multitude of common genes, each with a weak effect. This polygenic component of complex traits is mainly unknown and so makes up part of the missing heritability. Here, we try to integrate these different perspectives from the point of view of the metabolic changes observed in cancer.

Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Genome-wide association studies of cardiac electrical phenotypes

The genetic basis of cardiac electrical phenotypes has in the last 25 years been the subject of intense investigation. While in the first years, such efforts were dominated by the study of familial arrhythmia syndromes, in recent years, large consortia of investigators have successfully pursued genome-wide association studies (GWAS) for the identification of single-nucleotide polymorphisms that govern inter-individual variability in electrocardiographic parameters in the general population. We here provide a review of GWAS conducted on cardiac electrical phenotypes in the last 14 years and discuss the implications of these discoveries for our understanding of the genetic basis of disease susceptibility and variability in disease severity. Furthermore, we review functional follow-up studies that have been conducted on GWAS loci associated with cardiac electrical phenotypes and highlight the challenges and opportunities offered by such studies.

P-wave signal-averaged electrocardiography: Reference values, clinical correlates, and heritability in the Framingham Heart Study

BACKGROUND

P-wave signal-averaged electrocardiography (P-SAECG) quantifies atrial electrical activity. P-SAECG measures and their clinical correlates and heritability have had limited characterization in community-based cohorts.

OBJECTIVE

The purpose of this study was to (1) establish reference values; (2) identify clinical risk factors associated with P-SAECG; and (3) estimate genetic heritability for P-SAECG traits.

METHODS

We performed P-SAECG in 2 generations of Framingham Heart Study participants. We performed backward elimination regression models to assess associations of clinical factors with each SAECG trait (P-wave [PW] duration, root mean square voltage in terminal 40 ms [RMS40], terminal 30 ms RMS30, terminal 20 ms RMS20, RMS PW, and PW integral). We estimated the adjusted genetic heritability of P-SAECG measures using the Sequential Oligogenic Linkage Analysis Routines (SOLAR) program.

RESULTS

We included 4307 participants (age 55 ± 14 years; 56% female). The reference values were derived from 1752 participants without cardiovascular risk factors. Median (2.5th percentile; 97.5th percentile) total PW duration was 118 ms (93; 146) in women and 128 ms (104; 158) in men in the reference sample, and 121 ms (94; 151) in women and 129 ms (103; 159) in the entire study cohort (broad sample). In the broad sample, after adjusting for age and sex, total PW duration was positively associated with height, weight, prevalent heart failure, history of atrial fibrillation (AF), and atrioventricular node blockers, and negatively associated with smoking, waist circumference, heart rate, and diabetes. The estimated heritability of P-SAECG traits was moderate, ranging from 11.9% for RMS30 to 24.9% for PW integral.

CONCLUSION

P-SAECG traits are associated with multiple AF-related risk factors and are moderately heritable.

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.

Multi-Ethnic Genome-Wide Association Study of Decomposed Cardioelectric Phenotypes Illustrates Strategies to Identify and Characterize Evidence of Shared Genetic Effects for Complex Traits

BACKGROUND

We examined how expanding electrocardiographic trait genome-wide association studies to include ancestrally diverse populations, prioritize more precise phenotypic measures, and evaluate evidence for shared genetic effects enabled the detection and characterization of loci.

METHODS

We decomposed 10 seconds, 12-lead electrocardiograms from 34 668 multi-ethnic participants (15% Black; 30% Hispanic/Latino) into 6 contiguous, physiologically distinct (P wave, PR segment, QRS interval, ST segment, T wave, and TP segment) and 2 composite, conventional (PR interval and QT interval) interval scale traits and conducted multivariable-adjusted, trait-specific univariate genome-wide association studies using 1000-G imputed single-nucleotide polymorphisms. Evidence of shared genetic effects was evaluated by aggregating meta-analyzed univariate results across the 6 continuous electrocardiographic traits using the combined phenotype adaptive sum of powered scores test.

RESULTS

We identified 6 novels (CD36, PITX2, EMB, ZNF592, YPEL2, and BC043580) and 87 known loci (adaptive sum of powered score test P<5×10-9). Lead single-nucleotide polymorphism rs3211938 at CD36 was common in Blacks (minor allele frequency=10%), near monomorphic in European Americans, and had effects on the QT interval and TP segment that ranked among the largest reported to date for common variants. The other 5 novel loci were observed when evaluating the contiguous but not the composite electrocardiographic traits. Combined phenotype testing did not identify novel electrocardiographic loci unapparent using traditional univariate approaches, although this approach did assist with the characterization of known loci.

CONCLUSIONS

Despite including one-third as many participants as published electrocardiographic trait genome-wide association studies, our study identified 6 novel loci, emphasizing the importance of ancestral diversity and phenotype resolution in this era of ever-growing genome-wide association studies.

Familial Clustering of Cardiac Conduction Defects and Pacemaker Insertion

BACKGROUND

The etiopathogenesis of electrocardiographic bundle branch and atrioventricular blocks is not fully understood. We investigated familial clustering of cardiac conduction defects and pacemaker insertion in the FHS (Framingham Heart Study). Additionally, we assessed familial clustering of pacemaker insertion in the Danish general population.

METHODS

In FHS, we used multivariable-adjusted logistic regression models to investigate the association of parental atrioventricular block (PR interval, ≥0.2 s), complete bundle branch block (QRS, ≥0.12 s), or pacemaker insertion with the occurrence of cardiac conduction abnormalities in their offspring. The Danish nationwide administrative registries were interrogated to assess the relations of parental pacemaker insertion with offspring pacemaker insertion.

RESULTS

In FHS (n=371 cases with first-degree atrioventricular block, complete bundle branch block, or pacemaker insertion, and 1471 age- and sex-matched controls), individuals with at least 1 affected parent with a conduction defect had a 1.65-fold odds (odds ratio, 95% CI, 1.32-2.07) for manifesting an atrioventricular block and a 1.62-fold odds (95% CI, 1.08-2.42) for developing a complete bundle branch block. If at least 1 parent had any electrocardiographic conduction defect or pacemaker insertion, the offspring had a 1.62-fold odds (95% CI, 1.31-2.00) for experiencing any of these conditions. In Denmark (n=2 824 199 individuals; 5397 incident pacemaker implantations), individuals with at least 1 first-degree relative with history of pacemaker insertion had a multivariable-adjusted 1.68-fold (incidence rate ratio, 95% CI, 1.49-1.89) risk of undergoing a pacemaker insertion. If the affected relative was ≤45 years of age, the incidence rate ratio was markedly increased to 51.0 (95% CI, 32.7-79.9).

CONCLUSIONS

Cardiac conduction blocks and risk for pacemaker insertion cluster within families. A family history of conduction system disturbance or pacemaker insertion should trigger increased awareness of a similar propensity in other family members, especially so when the conduction system disease occurs at a younger age.

Age-related late-onset disease heritability patterns and implications for genome-wide association studies

Genome-wide association studies (GWASs) and other computational biology techniques are gradually discovering the causal gene variants that contribute to late-onset human diseases. After more than a decade of genome-wide association study efforts, these can account for only a fraction of the heritability implied by familial studies, the so-called "missing heritability" problem. Computer simulations of polygenic late-onset diseases (LODs) in an aging population have quantified the risk allele frequency decrease at older ages caused by individuals with higher polygenic risk scores (PRSs) becoming ill proportionately earlier. This effect is most prominent for diseases characterized by high cumulative incidence and high heritability, examples of which include Alzheimer's disease, coronary artery disease, cerebral stroke, and type 2 diabetes. The incidence rate for LODs grows exponentially for decades after early onset ages, guaranteeing that the cohorts used for GWASs overrepresent older individuals with lower PRSs, whose disease cases are disproportionately due to environmental causes such as old age itself. This mechanism explains the decline in clinical predictive power with age and the lower discovery power of familial studies of heritability and GWASs. It also explains the relatively constant-with-age heritability found for LODs of lower prevalence, exemplified by cancers.

Genome-wide association meta-analysis of 30,000 samples identifies seven novel loci for quantitative ECG traits

Genome-wide association studies (GWAS) of quantitative electrocardiographic (ECG) traits in large consortia have identified more than 130 loci associated with QT interval, QRS duration, PR interval, and heart rate (RR interval). In the current study, we meta-analyzed genome-wide association results from 30,000 mostly Dutch samples on four ECG traits: PR interval, QRS duration, QT interval, and RR interval. SNP genotype data was imputed using the Genome of the Netherlands reference panel encompassing 19 million SNPs, including millions of rare SNPs (minor allele frequency < 5%). In addition to many known loci, we identified seven novel locus-trait associations: KCND3, NR3C1, and PLN for PR interval, KCNE1, SGIP1, and NFKB1 for QT interval, and ATP2A2 for QRS duration, of which six were successfully replicated. At these seven loci, we performed conditional analyses and annotated significant SNPs (in exons and regulatory regions), demonstrating involvement of cardiac-related pathways and regulation of nearby genes.

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.

Pharmacogenomics study of thiazide diuretics and QT interval in multi-ethnic populations: the cohorts for heart and aging research in genomic epidemiology

Thiazide diuretics, commonly used antihypertensives, may cause QT interval (QT) prolongation, a risk factor for highly fatal and difficult to predict ventricular arrhythmias. We examined whether common single-nucleotide polymorphisms (SNPs) modified the association between thiazide use and QT or its component parts (QRS interval, JT interval) by performing ancestry-specific, trans-ethnic and cross-phenotype genome-wide analyses of European (66%), African American (15%) and Hispanic (19%) populations (N=78 199), leveraging longitudinal data, incorporating corrected standard errors to account for underestimation of interaction estimate variances and evaluating evidence for pathway enrichment. Although no loci achieved genome-wide significance (P<5 × 10-8), we found suggestive evidence (P<5 × 10-6) for SNPs modifying the thiazide-QT association at 22 loci, including ion transport loci (for example, NELL1, KCNQ3). The biologic plausibility of our suggestive results and simulations demonstrating modest power to detect interaction effects at genome-wide significant levels indicate that larger studies and innovative statistical methods are warranted in future efforts evaluating thiazide-SNP interactions.

A combined linkage, microarray and exome analysis suggests MAP3K11 as a candidate gene for left ventricular hypertrophy

BACKGROUND

Electrocardiographic measures of left ventricular hypertrophy (LVH) are used as predictors of cardiovascular risk. We combined linkage and association analyses to discover novel rare genetic variants involved in three such measures and two principal components derived from them.

METHODS

The study was conducted among participants from the Erasmus Rucphen Family Study (ERF), a Dutch family-based sample from the southwestern Netherlands. Variance components linkage analyses were performed using Merlin. Regions of interest (LOD > 1.9) were fine-mapped using microarray and exome sequence data.

RESULTS

We observed one significant LOD score for the second principal component on chromosome 15 (LOD score = 3.01) and 12 suggestive LOD scores. Several loci contained variants identified in GWAS for these traits; however, these did not explain the linkage peaks, nor did other common variants. Exome sequence data identified two associated variants after multiple testing corrections were applied.

CONCLUSIONS

We did not find common SNPs explaining these linkage signals. Exome sequencing uncovered a relatively rare variant in MAPK3K11 on chromosome 11 (MAF = 0.01) that helped account for the suggestive linkage peak observed for the first principal component. Conditional analysis revealed a drop in LOD from 2.01 to 0.88 for MAP3K11, suggesting that this variant may partially explain the linkage signal at this chromosomal location. MAP3K11 is related to the JNK pathway and is a pro-apoptotic kinase that plays an important role in the induction of cardiomyocyte apoptosis in various pathologies, including LVH.

Heritability in a SCN5A-mutation founder population with increased female susceptibility to non-nocturnal ventricular tachyarrhythmia and sudden cardiac death

BACKGROUND

Heritable cardiac-sodium channel dysfunction is associated with various arrhythmia syndromes, some predisposing to ventricular fibrillation. Phenotypic diversity among carriers of identical-by-descent mutations is often remarkable, suggesting influences of genetic modifiers.

OBJECTIVE

The purpose of this study was to identify a unique SCN5A-mutation founder population with mixed clinical phenotypes and sudden cardiac death, and to investigate the heritability of electromechanical traits besides the SCN5A-mutation effect.

METHODS

The 16-generation founder population segregating SCN5A c.4850_4852delTCT, p.(Phe1617del), was comprehensively phenotyped. Variance component analysis was used to evaluate the mutation's effects and assess heritability.

RESULTS

In 45 p.(Phe1617del) positives, the mutation associated strongly with QTc prolongation (472 ± 60 ms vs 423 ± 35 ms in 26 mutation negatives; P <.001; odds ratio for long-QT syndrome 22.4; 95% confidence interval 4.5-224.2; P <.001) and electromechanical window (EMW) negativity (-29 ± 47 ms vs 34 ± 26 ms; P <.001). Overlapping phenotypes including conduction delay and Brugada syndrome were noted in 19. Polymorphic ventricular tachyarrhythmias occurred mostly in the daytime, after arousal-evoked heart-rate acceleration and repolarization prolongation. Cox proportional hazards regression analysis revealed female gender as an independent risk factor for cardiac events (hazard ratio 5.1; 95% confidence interval 1.6-16.3; P = .006). p.(Phe1617del) was an important determinant of QTc_baseline, QTc_max, and EMW, explaining 18%, 28%, and 37%, respectively, of the trait's variance. Significant heritability was observed for PQ interval (P = .003) after accounting for the p.(Phe1617del) effect.

CONCLUSION

This SCN5A-p.(Phe1617del) founder population with phenotypic divergence and overlap reveals long-QT syndrome-related and arousal-evoked ventricular tachyarrhythmias with a female preponderance. Variance component analysis indicates additional genetic variance for PQ interval hidden in the genome, besides a dominant p.(Phe1617del) effect on QTc and EMW.

A Combined Linkage and Exome Sequencing Analysis for Electrocardiogram Parameters in the Erasmus Rucphen Family Study

Electrocardiogram (ECG) measurements play a key role in the diagnosis and prediction of cardiac arrhythmias and sudden cardiac death. ECG parameters, such as the PR, QRS, and QT intervals, are known to be heritable and genome-wide association studies of these phenotypes have been successful in identifying common variants; however, a large proportion of the genetic variability of these traits remains to be elucidated. The aim of this study was to discover loci potentially harboring rare variants utilizing variance component linkage analysis in 1547 individuals from a large family-based study, the Erasmus Rucphen Family Study (ERF). Linked regions were further explored using exome sequencing. Five suggestive linkage peaks were identified: two for QT interval (1q24, LOD = 2.63; 2q34, LOD = 2.05), one for QRS interval (1p35, LOD = 2.52) and two for PR interval (9p22, LOD = 2.20; 14q11, LOD = 2.29). Fine-mapping using exome sequence data identified a C > G missense variant (c.713C > G, p.Ser238Cys) in the FCRL2 gene associated with QT (rs74608430; P = 2.8 × 10^-4, minor allele frequency = 0.019). Heritability analysis demonstrated that the SNP explained 2.42% of the trait’s genetic variability in ERF (P = 0.02). Pathway analysis suggested that the gene is involved in cytosolic Ca²⁺ levels (P = 3.3 × 10^-3) and AMPK stimulated fatty acid oxidation in muscle (P = 4.1 × 10^-3). Look-ups in bioinformatics resources showed that expression of FCRL2 is associated with ARHGAP24 and SETBP1 expression. This finding was not replicated in the Rotterdam study. Combining the bioinformatics information with the association and linkage analyses, FCRL2 emerges as a strong candidate gene for QT interval.

Missing heritability of complex diseases: Enlightenment by genetic variants from intermediate phenotypes

Diseases of complex origin have a component of quantitative genetics that contributes to their susceptibility and phenotypic variability. However, after several studies, a major part of the genetic component of complex phenotypes has still not been found, a situation known as "missing heritability." Although there have been many hypotheses put forward to explain the reasons for the missing heritability, its definitive causes remain unknown. Complex diseases are caused by multiple intermediate phenotypes involved in their pathogenesis and, very often, each one of these intermediate phenotypes also has a component of quantitative inheritance. Here we propose that at least part of the missing heritability can be explained by the genetic component of intermediate phenotypes that is not detectable at the level of the main complex trait. At the same time, the identification of the genetic component of intermediate phenotypes provides an opportunity to identify part of the missing heritability of complex diseases.