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

Deficiency in Prader-Willi syndrome gene necdin leads to attenuated cardiac contractility

Prader-Willi syndrome (PWS) is a genetic disorder characterized by behavioral disturbances, hyperphagia, and intellectual disability. Several surveys indicate that PWS is also associated with cardiac abnormalities, possibly contributing to a high incidence of sudden death. However, the pathological mechanisms underlying cardiac dysfunction in PWS remain unclear. In this study, we found that deficiency in necdin, an intronless gene within PWS region, led to heart systolic and diastolic dysfunction in mice. Through yeast two-hybrid screening, we identified an interaction between necdin and non-muscle myosin regulatory light chain 12a/b (MYL12 A/B). We further showed that necdin stabilized MYL12 A/B via SGT1-heat shock protein 90 (HSP90) chaperone machinery. The zebrafish lacking the MYL12 A/B analog, MYL12.1, exhibited impaired heart function, while cardiac-specific overexpression of MYL12A normalized the heart dysfunction in necdin-deficient mice. Our findings revealed necdin dysfunction as a contributing factor to cardiomyopathy in PWS patients and emphasized the importance of HSP90 chaperone machinery and non-muscle myosin in heart fitness.

Genetic evidence for T-wave area from 12-lead electrocardiograms to monitor cardiovascular diseases in patients taking diabetes medications

Aims Many studies indicated use of diabetes medications can influence the electrocardiogram (ECG), which remains the simplest and fastest tool for assessing cardiac functions. However, few studies have explored the role of genetic factors in determining the relationship between the use of diabetes medications and ECG trace characteristics (ETC). Methods Genome-wide association studies (GWAS) were performed for 168 ETCs extracted from the 12-lead ECGs of 42,340 Europeans in the UK Biobank. The genetic correlations, causal relationships, and phenotypic relationships of these ETCs with medication usage, as well as the risk of cardiovascular diseases (CVDs), were estimated by linkage disequilibrium score regression (LDSC), Mendelian randomization (MR), and regression model, respectively. Results The GWAS identified 124 independent single nucleotide polymorphisms (SNPs) that were study-wise and genome-wide significantly associated with at least one ETC. Regression model and LDSC identified significant phenotypic and genetic correlations of T-wave area in lead aVR (aVR_T-area) with usage of diabetes medications (ATC code: A10 drugs, and metformin), and the risks of ischemic heart disease (IHD) and coronary atherosclerosis (CA). MR analyses support a putative causal effect of the use of diabetes medications on decreasing aVR_T-area, and on increasing risk of IHD and CA. ConclusionPatients taking diabetes medications are prone to have decreased aVR_T-area and an increased risk of IHD and CA. The aVR_T-area is therefore a potential ECG marker for pre-clinical prediction of IHD and CA in patients taking diabetes medications.

Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle

BACKGROUND

Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified. Our aim was to search for genomic regions associated with this disease.

METHODS

A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study. Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar's test.

RESULTS

Two independent genomic regions were identified as having the most significant association with BCHF: the arrestin domain-containing protein 3 gene ( ARRDC3), and the nuclear factor IA gene ( NFIA, mid- p-values, 1x10 -8 and 2x10 -7, respectively). Animals with two copies of risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates with either one or zero risk alleles at both genes (CI 95 = 3-17). Further, animals with two copies of risk alleles at both genes were 28-fold more likely to have BCHF than all others ( p-value = 1×10 -7, CI 95 = 4-206). A missense variant in ARRDC3 (C182Y) represents a potential functional variant since the C182 codon is conserved among all other jawed vertebrate species observed. A two-SNP test with markers in both genes showed 29% of 273 BCHF cases had homozygous risk genotypes in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This and other DNA tests may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here.

CONCLUSIONS

Although pathogenic roles for variants in the ARRDC3 and NFIA genes are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle

Background

Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified. Our aim was to search for genomic regions associated with this disease.

Methods

A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study. Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar's test.

Results

Two independent genomic regions were identified as having the most significant association with BCHF: the arrestin domain-containing protein 3 gene ( ARRDC3), and the nuclear factor IA gene ( NFIA, mid- p-values, 1x10 -8 and 2x10 -7, respectively). Animals with two copies of risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates with either one or zero risk alleles at both genes (CI 95 = 3-17). Further, animals with two copies of risk alleles at both genes were 28-fold more likely to have BCHF than all others ( p-value = 1×10 -7, CI 95 = 4-206). A missense variant in ARRDC3 (C182Y) represents a potential functional variant since the C182 codon is conserved among all other jawed vertebrate species observed. A two-SNP test with markers in both genes showed 29% of 273 BCHF cases had homozygous risk genotypes in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This and other DNA tests may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here.

Conclusions

Although pathogenic roles for variants in the ARRDC3 and NFIA genes are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Genetic variants and effect modifiers of QT interval prolongation in patients with sickle cell disease

BACKGROUND

Sickle cell disease (SCD) is a common inherited blood disorder among African Americans (AA), with premature mortality which has been associated with prolongation of the heart rate-corrected QT interval (QTc), a known risk factor for sudden cardiac death. Although numerous genetic variants have been identified as contributors to QT interval prolongation in the general population, their impact on SCD patients remains unclear. This study used an unweighted polygenic risk score (PRS) to validate the previously identified associations between SNPs and QTc interval in SCD patients, and to explore possible interactions with other factors that prolong QTc interval in AA individuals with SCD.

METHODS

In SCD patients, candidate genetic variants associated with the QTc interval were genotyped. To identify any risk SNPs that may be correlated with QTc interval prolongation, linear regression was employed, and an unweighted PRS was subsequently constructed. The effect of PRS on the QTc interval was evaluated using linear regression, while stratification analysis was used to assess the influence of serum alanine transaminase (ALT), a biomarker for liver disease, on the PRS effect. We also evaluated the PRS with the two subcomponents of QTc, the QRS and JTc intervals.

RESULTS

Out of 26 candidate SNPs, five risk SNPs were identified for QTc duration under the recessive model. For every unit increase in PRS, the QTc interval prolonged by 4.0 ms (95% CI: [2.0, 6.1]; p-value: <0.001) in the additive model and 9.4 ms in the recessive model (95% CI: [4.6, 14.1]; p-value: <0.001). Serum ALT showed a modification effect on PRS-QTc prolongation under the recessive model. In the normal ALT group, each PRS unit increased QTc interval by 11.7 ms (95% CI: [6.3, 17.1]; p-value: 2.60E-5), whereas this effect was not observed in the elevated ALT group (0.9 ms; 95% CI: [-7.0, 8.8]; p-value: 0.823).

CONCLUSION

Several candidate genetic variants are associated with QTc interval prolongation in SCD patients, and serum ALT acts as a modifying factor. The association of a CPS1 gene variant in both QTc and JTc duration adds to NOS1AP as evidence of involvement of the urea cycle and nitric oxide metabolism in cardiac repolarization in SCD. Larger replication studies are needed to confirm these findings and elucidate the underlying mechanisms.

Discovery of TBX20 as a Novel Gene Underlying Atrial Fibrillation

Atrial fibrillation (AF), the most prevalent type of sustained cardiac dysrhythmia globally, confers strikingly enhanced risks for cognitive dysfunction, stroke, chronic cardiac failure, and sudden cardiovascular demise. Aggregating studies underscore the crucial roles of inherited determinants in the occurrence and perpetuation of AF. However, due to conspicuous genetic heterogeneity, the inherited defects accounting for AF remain largely indefinite. Here, via whole-genome genotyping with genetic markers and a linkage assay in a family suffering from AF, a new AF-causative locus was located at human chromosome 7p14.2-p14.3, a ~4.89 cM (~4.43-Mb) interval between the markers D7S526 and D7S2250. An exome-wide sequencing assay unveiled that, at the defined locus, the mutation in the TBX20 gene, NM_001077653.2: c.695A>G; p.(His232Arg), was solely co-segregated with AF in the family. Additionally, a Sanger sequencing assay of TBX20 in another family suffering from AF uncovered a novel mutation, NM_001077653.2: c.862G>C; p.(Asp288His). Neither of the two mutations were observed in 600 unrelated control individuals. Functional investigations demonstrated that the two mutations both significantly reduced the transactivation of the target gene KCNH2 (a well-established AF-causing gene) and the ability to bind the promoter of KCNH2, while they had no effect on the nuclear distribution of TBX20. Conclusively, these findings reveal a new AF-causative locus at human chromosome 7p14.2-p14.3 and strongly indicate TBX20 as a novel AF-predisposing gene, shedding light on the mechanism underlying AF and suggesting clinical significance for the allele-specific treatment of AF patients.

Importance of Diversity in Precision Medicine: Generalizability of Genetic Associations Across Ancestry Groups Toward Better Identification of Disease Susceptibility Variants

Genome-wide association studies (GWAS) revolutionized our understanding of common genetic variation and its impact on common human disease and traits. Developed and adopted in the mid-2000s, GWAS led to searchable genotype-phenotype catalogs and genome-wide datasets available for further data mining and analysis for the eventual development of translational applications. The GWAS revolution was swift and specific, including almost exclusively populations of European descent, to the neglect of the majority of the world's genetic diversity. In this narrative review, we recount the GWAS landscape of the early years that established a genotype-phenotype catalog that is now universally understood to be inadequate for a complete understanding of complex human genetics. We then describe approaches taken to augment the genotype-phenotype catalog, including the study populations, collaborative consortia, and study design approaches aimed to generalize and then ultimately discover genome-wide associations in non-European descent populations. The collaborations and data resources established in the efforts to diversify genomic findings undoubtedly provide the foundations of the next chapters of genetic association studies with the advent of budget-friendly whole-genome sequencing.

Antithrombin, Protein C, and Protein S: Genome and Transcriptome-Wide Association Studies Identify 7 Novel Loci Regulating Plasma Levels

BACKGROUND

Antithrombin, PC (protein C), and PS (protein S) are circulating natural anticoagulant proteins that regulate hemostasis and of which partial deficiencies are causes of venous thromboembolism. Previous genetic association studies involving antithrombin, PC, and PS were limited by modest sample sizes or by being restricted to candidate genes. In the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, we meta-analyzed across ancestries the results from 10 genome-wide association studies of plasma levels of antithrombin, PC, PS free, and PS total.

METHODS

Study participants were of European and African ancestries, and genotype data were imputed to TOPMed, a dense multiancestry reference panel. Each of the 10 studies conducted a genome-wide association studies for each phenotype and summary results were meta-analyzed, stratified by ancestry. Analysis of antithrombin included 25 243 European ancestry and 2688 African ancestry participants, PC analysis included 16 597 European ancestry and 2688 African ancestry participants, PSF and PST analysis included 4113 and 6409 European ancestry participants. We also conducted transcriptome-wide association analyses and multiphenotype analysis to discover additional associations. Novel genome-wide association studies and transcriptome-wide association analyses findings were validated by in vitro functional experiments. Mendelian randomization was performed to assess the causal relationship between these proteins and cardiovascular outcomes.

RESULTS

Genome-wide association studies meta-analyses identified 4 newly associated loci: 3 with antithrombin levels (GCKR, BAZ1B, and HP-TXNL4B) and 1 with PS levels (ORM1-ORM2). transcriptome-wide association analyses identified 3 newly associated genes: 1 with antithrombin level (FCGRT), 1 with PC (GOLM2), and 1 with PS (MYL7). In addition, we replicated 7 independent loci reported in previous studies. Functional experiments provided evidence for the involvement of GCKR, SNX17, and HP genes in antithrombin regulation.

CONCLUSIONS

The use of larger sample sizes, diverse populations, and a denser imputation reference panel allowed the detection of 7 novel genomic loci associated with plasma antithrombin, PC, and PS levels.

Genetic architecture of heart failure with preserved versus reduced ejection fraction

Pharmacologic clinical trials for heart failure with preserved ejection fraction have been largely unsuccessful as compared to those for heart failure with reduced ejection fraction. Whether differences in the genetic underpinnings of these major heart failure subtypes may provide insights into the disparate outcomes of clinical trials remains unknown. We utilize a large, uniformly phenotyped, single cohort of heart failure sub-classified into heart failure with reduced and with preserved ejection fractions based on current clinical definitions, to conduct detailed genetic analyses of the two heart failure sub-types. We find different genetic architectures and distinct genetic association profiles between heart failure with reduced and with preserved ejection fraction suggesting differences in underlying pathobiology. The modest genetic discovery for heart failure with preserved ejection fraction (one locus) compared to heart failure with reduced ejection fraction (13 loci) despite comparable sample sizes indicates that clinically defined heart failure with preserved ejection fraction likely represents the amalgamation of several, distinct pathobiological entities. Development of consensus sub-phenotyping of heart failure with preserved ejection fraction is paramount to better dissect the underlying genetic signals and contributors to this highly prevalent condition.

A large genome-wide association study of QT interval length utilizing electronic health records

QT interval length is an important risk factor for adverse cardiovascular outcomes; however, the genetic architecture of QT interval remains incompletely understood. We conducted a genome-wide association study of 76,995 ancestrally diverse Kaiser Permanente Northern California members enrolled in the Genetic Epidemiology Research on Adult Health and Aging cohort using 448,517 longitudinal QT interval measurements, uncovering 9 novel variants, most replicating in 40,537 individuals in the UK Biobank and Population Architecture using Genomics and Epidemiology studies. A meta-analysis of all 3 cohorts (n = 117,532) uncovered an additional 19 novel variants. Conditional analysis identified 15 additional variants, 3 of which were novel. Little, if any, difference was seen when adjusting for putative QT interval lengthening medications genome-wide. Using multiple measurements in Genetic Epidemiology Research on Adult Health and Aging increased variance explained by 163%, and we show that the ≈6 measurements in Genetic Epidemiology Research on Adult Health and Aging was equivalent to a 2.4× increase in sample size of a design with a single measurement. The array heritability was estimated at ≈17%, approximately half of our estimate of 36% from family correlations. Heritability enrichment was estimated highest and most significant in cardiovascular tissue (enrichment 7.2, 95% CI = 5.7-8.7, P = 2.1e-10), and many of the novel variants included expression quantitative trait loci in heart and other relevant tissues. Comparing our results to other cardiac function traits, it appears that QT interval has a multifactorial genetic etiology.

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.

DeepLUCIA: predicting tissue-specific chromatin loops using Deep Learning-based Universal Chromatin Interaction Annotator

MOTIVATION

The importance of chromatin loops in gene regulation is broadly accepted. There are mainly two approaches to predict chromatin loops: transcription factor (TF) binding-dependent approach and genomic variation-based approach. However, neither of these approaches provides an adequate understanding of gene regulation in human tissues. To address this issue, we developed a deep learning-based chromatin loop prediction model called DeepLUCIA (Deep Learning-based Universal Chromatin Interaction Annotator).

RESULTS

Although DeepLUCIA does not use TF binding profile data which previous TF binding-dependent methods critically rely on, its prediction accuracies are comparable to those of the previous TF binding-dependent methods. More importantly, DeepLUCIA enables the tissue-specific chromatin loop predictions from tissue-specific epigenomes that cannot be handled by genomic variation-based approach. We demonstrated the utility of the DeepLUCIA by predicting several novel target genes of SNPs identified in genome-wide association studies targeting Brugada syndrome, COVID-19 severity, and age-related macular degeneration.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of genetic effects underlying type 2 diabetes in South Asian and European populations

South Asians are at high risk of developing type 2 diabetes (T2D). We carried out a genome-wide association meta-analysis with South Asian T2D cases (n = 16,677) and controls (n = 33,856), followed by combined analyses with Europeans (neff = 231,420). We identify 21 novel genetic loci for significant association with T2D (P = 4.7 × 10-8 to 5.2 × 10-12), to the best of our knowledge at the point of analysis. The loci are enriched for regulatory features, including DNA methylation and gene expression in relevant tissues, and highlight CHMP4B, PDHB, LRIG1 and other genes linked to adiposity and glucose metabolism. A polygenic risk score based on South Asian-derived summary statistics shows ~4-fold higher risk for T2D between the top and bottom quartile. Our results provide further insights into the genetic mechanisms underlying T2D, and highlight the opportunities for discovery from joint analysis of data from across ancestral populations.

Environmentally sensitive hotspots in the methylome of the early human embryo

In humans, DNA methylation marks inherited from gametes are largely erased following fertilisation, prior to construction of the embryonic methylome. Exploiting a natural experiment of seasonal variation including changes in diet and nutritional status in rural Gambia, we analysed three datasets covering two independent child cohorts and identified 259 CpGs showing consistent associations between season of conception (SoC) and DNA methylation. SoC effects were most apparent in early infancy, with evidence of attenuation by mid-childhood. SoC-associated CpGs were enriched for metastable epialleles, parent-of-origin-specific methylation and germline differentially methylated regions, supporting a periconceptional environmental influence. Many SoC-associated CpGs overlapped enhancers or sites of active transcription in H1 embryonic stem cells and fetal tissues. Half were influenced but not determined by measured genetic variants that were independent of SoC. Environmental 'hotspots' providing a record of environmental influence at periconception constitute a valuable resource for investigating epigenetic mechanisms linking early exposures to lifelong health and disease.

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.

HAND transcription factors cooperatively specify the aorta and pulmonary trunk

Congenital heart defects (CHDs) affecting the cardiac outflow tract (OFT) constitute a significant cause of morbidity and mortality. The OFT develops from migratory cell populations which include the cardiac neural crest cells (cNCCs) and secondary heart field (SHF) derived myocardium and endocardium. The related transcription factors HAND1 and HAND2 have been implicated in human CHDs involving the OFT. Although Hand1 is expressed within the OFT, Hand1 NCC-specific conditional knockout mice (H1CKOs) are viable. Here we show that these H1CKOs present a low penetrance of OFT phenotypes, whereas SHF-specific Hand1 ablation does not reveal any cardiac phenotypes. Further, HAND1 and HAND2 appear functionally redundant within the cNCCs, as a reduction/ablation of Hand2 on an NCC-specific H1CKO background causes pronounced OFT defects. Double conditional Hand1 and Hand2 NCC knockouts exhibit persistent truncus arteriosus (PTA) with 100% penetrance. NCC lineage-tracing and Sema3c in situ mRNA expression reveal that Sema3c-expressing cells are mis-localized, resulting in a malformed septal bridge within the OFTs of H1CKO;H2CKO embryos. Interestingly, Hand1 and Hand2 also genetically interact within the SHF, as SHF H1CKOs on a heterozygous Hand2 background exhibit Ventricular Septal Defects (VSDs) with incomplete penetrance. Previously, we identified a BMP, HAND2, and GATA-dependent Hand1 OFT enhancer sufficient to drive reporter gene expression within the nascent OFT and aorta. Using these transcription inputs as a probe, we identify a novel Hand2 OFT enhancer, suggesting that a conserved BMP-GATA dependent mechanism transcriptionally regulates both HAND factors. These findings support the hypothesis that HAND factors interpret BMP signaling within the cNCCs to cooperatively coordinate OFT morphogenesis.

Mog1 knockout causes cardiac hypertrophy and heart failure by downregulating tbx5-cryab-hspb2 signalling in zebrafish

AIMS

MOG1 is a small protein that can bind to small GTPase RAN and regulate transport of RNA and proteins between the cytoplasm and nucleus. However, the in vivo physiological role of mog1 in the heart needs to be fully defined.

METHODS

Mog1 knockout zebrafish was generated by TALEN. Echocardiography, histological analysis, and electrocardiograms were used to examine cardiac structure and function. RNA sequencing and real-time RT-PCR were used to elucidate the molecular mechanism and to analyse the gene expression. Isoproterenol was used to induce cardiac hypertrophy. Whole-mount in situ hybridization was used to observe cardiac morphogenesis.

RESULTS

Mog1 knockout zebrafish developed cardiac hypertrophy and heart failure (enlarged pericardium, increased nppa and nppb expression and ventricular wall thickness, and reduced ejection fraction), which was aggravated by isoproterenol. RNAseq and KEGG pathway analyses revealed the effect of mog1 knockout on the pathways of cardiac hypertrophy, dilatation and contraction. Mechanistic studies revealed that mog1 knockout decreased expression of tbx5, which reduced expression of cryab and hspb2, resulting in cardiac hypertrophy and heart failure. Overexpression of cryab, hspb2 and tbx5 rescued the cardiac oedema phenotype of mog1 KO zebrafish. Telemetry electrocardiogram monitoring showed QRS and QTc prolongation and a reduced heart rate in mog1 knockout zebrafish, which was associated with reduced scn1b expression. Moreover, mog1 knockout resulted in abnormal cardiac looping during embryogenesis because of the reduced expression of nkx2.5, gata4 and hand2.

CONCLUSION

Our data identified an important molecular determinant for cardiac hypertrophy and heart failure, and rhythm maintenance of the heart.

The Role of Tbx20 in Cardiovascular Development and Function

Tbx20 is a member of the Tbx1 subfamily of T-box-containing genes and is known to play a variety of fundamental roles in cardiovascular development and homeostasis as well as cardiac remodeling in response to pathophysiological stresses. Mutations in TBX20 are widely associated with the complex spectrum of congenital heart defects (CHDs) in humans, which includes defects in chamber septation, chamber growth, and valvulogenesis. In addition, genetic variants of TBX20 have been found to be associated with dilated cardiomyopathy and heart arrhythmia. This broad spectrum of cardiac morphogenetic and functional defects is likely due to its broad expression pattern in multiple cardiogenic cell lineages and its critical regulation of transcriptional networks during cardiac development. In this review, we summarize recent findings in our general understanding of the role of Tbx20 in regulating several important aspects of cardiac development and homeostasis and heart function.

Association of T66A polymorphism in CASQ2 with PR interval in a Chinese population

OBJECTIVE

The aim of this study was to explore the relationship between arrhythmia-associated or electrocardiogram (ECG)-associated common variants and PR interval, QRS duration, QT_corrected, and heart rate in a Chinese cohort.

METHODS

We studied the association between 26 single-nucleotide polymorphisms (SNPs) and digital ECG data from 379 unrelated Han Chinese individuals collected in an epidemiological survey in Beijing. All subjects were 45 years of age or older and were free of cardiovascular diseases and diabetes. The SNPs were genotyped in a multiplex panel using the Sequenom MassARRAY platform.

RESULTS

Missense variant T66A (Thr66Ala, rs4074536) of the CASQ2 gene, which was previously reported to be associated with QRS complex in European populations, was significantly associated with PR interval prolongation in our sample (p_adjusted = 0.006, beta_adjusted = 3.983 ms). A two-tailed t test showed that the CC genotype (n = 86) had a significantly longer PR interval (162.9 ± 19.4 ms) than the non-CC genotypes (n = 288, PR interval: 154.6 ± 20.9 ms), with a remarkable difference of 8.2 ms between the groups (p = 0.001). Interestingly, this association between T66A of CASQ2 and PR interval was more evident in females (p_adjusted = 0.007, beta_adjusted = 5.723 ms) than in males (p_adjusted = 0.177, beta_adjusted = 2.725 ms). In addition, rs3822714 in the HAND1 locus might be associated with QRS duration (p_adjusted = 0.034, beta_adjusted = -2.268 ms).

CONCLUSION

We identified a novel signal of an association between the CC genotype of T66A in CASQ2 and PR interval prolongation in a Chinese population, particularly in females. This association deserves further exploration given its possible effects on calcium handling in cardiac electrophysiology.

A phenome-wide association study (PheWAS) in the Population Architecture using Genomics and Epidemiology (PAGE) study reveals potential pleiotropy in African Americans

We performed a hypothesis-generating phenome-wide association study (PheWAS) to identify and characterize cross-phenotype associations, where one SNP is associated with two or more phenotypes, between thousands of genetic variants assayed on the Metabochip and hundreds of phenotypes in 5,897 African Americans as part of the Population Architecture using Genomics and Epidemiology (PAGE) I study. The PAGE I study was a National Human Genome Research Institute-funded collaboration of four study sites accessing diverse epidemiologic studies genotyped on the Metabochip, a custom genotyping chip that has dense coverage of regions in the genome previously associated with cardio-metabolic traits and outcomes in mostly European-descent populations. Here we focus on identifying novel phenome-genome relationships, where SNPs are associated with more than one phenotype. To do this, we performed a PheWAS, testing each SNP on the Metabochip for an association with up to 273 phenotypes in the participating PAGE I study sites. We identified 133 putative pleiotropic variants, defined as SNPs associated at an empirically derived p-value threshold of p<0.01 in two or more PAGE study sites for two or more phenotype classes. We further annotated these PheWAS-identified variants using publicly available functional data and local genetic ancestry. Amongst our novel findings is SPARC rs4958487, associated with increased glucose levels and hypertension. SPARC has been implicated in the pathogenesis of diabetes and is also known to have a potential role in fibrosis, a common consequence of multiple conditions including hypertension. The SPARC example and others highlight the potential that PheWAS approaches have in improving our understanding of complex disease architecture by identifying novel relationships between genetic variants and an array of common human phenotypes.

Sequence variants with large effects on cardiac electrophysiology and disease

Features of the QRS complex of the electrocardiogram, reflecting ventricular depolarisation, associate with various physiologic functions and several pathologic conditions. We test 32.5 million variants for association with ten measures of the QRS complex in 12 leads, using 405,732 electrocardiograms from 81,192 Icelanders. We identify 190 associations at 130 loci, the majority of which have not been reported before, including associations with 21 rare or low-frequency coding variants. Assessment of genes expressed in the heart yields an additional 13 rare QRS coding variants at 12 loci. We find 51 unreported associations between the QRS variants and echocardiographic traits and cardiovascular diseases, including atrial fibrillation, complete AV block, heart failure and supraventricular tachycardia. We demonstrate the advantage of in-depth analysis of the QRS complex in conjunction with other cardiovascular phenotypes to enhance our understanding of the genetic basis of myocardial mass, cardiac conduction and disease.

Aberrant morphology of the QRS complex in an electrocardiogram can be associated with cardiac morbidity and mortality. Here, the authors perform genome-wide association studies for ten measures of the QRS complex in 81,192 individuals and find 86 previously unreported loci that associate with at least one parameter.

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.

Variation in a Left Ventricle-Specific Hand1 Enhancer Impairs GATA Transcription Factor Binding and Disrupts Conduction System Development and Function

RATIONALE

The ventricular conduction system (VCS) rapidly propagates electrical impulses through the working myocardium of the ventricles to coordinate chamber contraction. GWAS (Genome-wide association studies) have associated nucleotide polymorphisms, most are located within regulatory intergenic or intronic sequences, with variation in VCS function. Two highly correlated polymorphisms (r²>0.99) associated with VCS functional variation (rs13165478 and rs13185595) occur 5' to the gene encoding the basic helix-loop-helix transcription factor HAND1 (heart- and neural crest derivatives-expressed protein 1).

OBJECTIVE

Here, we test the hypothesis that these polymorphisms influence HAND1 transcription thereby influencing VCS development and function.

METHODS AND RESULTS

We employed transgenic mouse models to identify an enhancer that is sufficient for left ventricle (LV) cis-regulatory activity. Two evolutionarily conserved GATA transcription factor cis-binding elements within this enhancer are bound by GATA4 and are necessary for cis-regulatory activity, as shown by in vitro DNA binding assays. CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9-mediated deletion of this enhancer dramatically reduces Hand1 expression solely within the LV but does not phenocopy previously published mouse models of cardiac Hand1 loss-of-function. Electrophysiological and morphological analyses reveals that mice homozygous for this deleted enhancer display a morphologically abnormal VCS and a conduction system phenotype consistent with right bundle branch block. Using 1000 Genomes Project data, we identify 3 additional single nucleotide polymorphisms (SNPs), located within the Hand1 LV enhancer, that compose a haplotype with rs13165478 and rs13185595. One of these SNPs, rs10054375, overlaps with a critical GATA cis-regulatory element within the Hand1 LV enhancer. This SNP, when tested in electrophoretic mobility shift assays, disrupts GATA4 DNA-binding. Modeling 2 of these SNPs in mice causes diminished Hand1 expression and mice present with abnormal VCS function.

CONCLUSIONS

Together, these findings reveal that SNP rs10054375, which is located within a necessary and sufficient LV-specific Hand1 enhancer, exhibits reduces GATA DNA-binding in electrophoretic mobility shift assay, and this enhancer in total, is required for VCS development and function in mice and perhaps humans.

GWAS of QRS duration identifies new loci specific to Hispanic/Latino populations

BACKGROUND

The electrocardiographically quantified QRS duration measures ventricular depolarization and conduction. QRS prolongation has been associated with poor heart failure prognosis and cardiovascular mortality, including sudden death. While previous genome-wide association studies (GWAS) have identified 32 QRS SNPs across 26 loci among European, African, and Asian-descent populations, the genetics of QRS among Hispanics/Latinos has not been previously explored.

METHODS

We performed a GWAS of QRS duration among Hispanic/Latino ancestry populations (n = 15,124) from four studies using 1000 Genomes imputed genotype data (adjusted for age, sex, global ancestry, clinical and study-specific covariates). Study-specific results were combined using fixed-effects, inverse variance-weighted meta-analysis.

RESULTS

We identified six loci associated with QRS (P<5x10-8), including two novel loci: MYOCD, a nuclear protein expressed in the heart, and SYT1, an integral membrane protein. The top SNP in the MYOCD locus, intronic SNP rs16946539, was found in Hispanics/Latinos with a minor allele frequency (MAF) of 0.04, but is monomorphic in European and African descent populations. The most significant QRS duration association was with intronic SNP rs3922344 (P = 1.19x10-24) in SCN5A/SCN10A. Three other previously identified loci, CDKN1A, VTI1A, and HAND1, also exceeded the GWAS significance threshold among Hispanics/Latinos. A total of 27 of 32 previously identified QRS duration SNPs were shown to generalize in Hispanics/Latinos.

CONCLUSIONS

Our QRS duration GWAS, the first in Hispanic/Latino populations, identified two new loci, underscoring the utility of extending large scale genomic studies to currently under-examined populations.

Hemodynamic and ECG responses to stress test in early adolescent athletes explain ethnicity-related cardiac differences

BACKGROUND

Ethnicity is an important determinant of athletes' cardiovascular adaptation. Black adolescent and adult athletes exhibit a left ventricular (LV) hypertrophy with a concentric remodelling higher than their Caucasian counterparts. Scant data, however, are available on race-related differences in hemodynamic response of adolescent athletes to exercise and its relation with heart remodelling. We evaluated if race-specific, sport-related structural and electrical remodelling in adolescent athletes of Caucasian and African ethnicity exclusively depends on race itself rather than on different cardiovascular responses to physical exercise.

METHODS

We examined 90 adolescent athletes, 60 Caucasian (WA) and 30 Black (BA). All participants underwent thorough clinical, echocardiographic and stress test evaluations.

RESULTS

BA had greater indexed LV mass (LVM/BSA) with increased relative wall thickness (RWT) implying a concentric remodelling. BA showed higher systolic blood pressure (SBP) compared to WA during the whole exercise test. ECG data showed that BA vs WA had a significant shorter QRS duration in each step considered with a significant greater QT dispersion. BA reached a higher relative pressure peak as compared to WA. RWT was strongly influenced by ethnicity and less by SBP at peak of exercise (PE), although LVM/BSA was significantly related to SBP at PE and just marginally to age and not significantly to race.

CONCLUSIONS

Black adolescent athletes showed higher SBP during all steps of exercise associated to a different trend. Ethnicity was the main determinant of RWT, suggesting that LV geometry is principally race-related rather than influenced by a different hemodynamic profile to physical activity.

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.

GWAS of the electrocardiographic QT interval in Hispanics/Latinos generalizes previously identified loci and identifies population-specific signals

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.