Genetics of heart rate in heart failure patients (GenHRate)

Structural variation in humans and our primate kin in the era of telomere-to-telomere genomes and pangenomics

Structural variants (SVs) account for the majority of base pair differences both within and between primate species. However, our understanding of inter- and intra-species SV has been historically hampered by the quality of draft primate genomes and the absence of genome resources for key taxa. Recently, advances in long-read sequencing and genome assembly have begun to radically reshape our understanding of SVs. Two landmark achievements include the publication of a human telomere-to-telomere (T2T) genome as well as the development of the first human pangenome reference. In this review, we first look back to the major works laying the foundation for these projects. We then examine the ways in which T2T genome assemblies and pangenomes are transforming our understanding of and approach to primate SV. Finally, we discuss what the future of primate SV research may look like in the era of T2T genomes and pangenomics.

Using Omics to Identify Novel Therapeutic Targets in Heart Failure

Omics refers to the measurement and analysis of the totality of molecules or biological processes involved within an organism. Examples of omics data include genomics, transcriptomics, epigenomics, proteomics, metabolomics, and more. In this review, we present the available literature reporting omics data on heart failure that can inform the development of novel treatments or innovative treatment strategies for this disease. This includes polygenic risk scores to improve prediction of genomic data and the potential of multiomics to more efficiently identify potential treatment targets for further study. We also discuss the limitations of omic analyses and the barriers that must be overcome to maximize the utility of these types of studies. Finally, we address the current state of the field and future opportunities for using multiomics to better personalize heart failure treatment strategies.

Novel functional insights into ischemic stroke biology provided by the first genome-wide association study of stroke in indigenous Africans

BACKGROUND

African ancestry populations have the highest burden of stroke worldwide, yet the genetic basis of stroke in these populations is obscure. The Stroke Investigative Research and Educational Network (SIREN) is a multicenter study involving 16 sites in West Africa. We conducted the first-ever genome-wide association study (GWAS) of stroke in indigenous Africans.

METHODS

Cases were consecutively recruited consenting adults (aged > 18 years) with neuroimaging-confirmed ischemic stroke. Stroke-free controls were ascertained using a locally validated Questionnaire for Verifying Stroke-Free Status. DNA genotyping with the H3Africa array was performed, and following initial quality control, GWAS datasets were imputed into the NIH Trans-Omics for Precision Medicine (TOPMed) release2 from BioData Catalyst. Furthermore, we performed fine-mapping, trans-ethnic meta-analysis, and in silico functional characterization to identify likely causal variants with a functional interpretation.

RESULTS

We observed genome-wide significant (P-value < 5.0E-8) SNPs associations near AADACL2 and miRNA (MIR5186) genes in chromosome 3 after adjusting for hypertension, diabetes, dyslipidemia, and cardiac status in the base model as covariates. SNPs near the miRNA (MIR4458) gene in chromosome 5 were also associated with stroke (P-value < 1.0E-6). The putative genes near AADACL2, MIR5186, and MIR4458 genes were protective and novel. SNPs associations with stroke in chromosome 2 were more than 77 kb from the closest gene LINC01854 and SNPs in chromosome 7 were more than 116 kb to the closest gene LINC01446 (P-value < 1.0E-6). In addition, we observed SNPs in genes STXBP5-AS1 (chromosome 6), GALTN9 (chromosome 12), FANCA (chromosome 16), and DLGAP1 (chromosome 18) (P-value < 1.0E-6). Both genomic regions near genes AADACL2 and MIR4458 remained significant following fine mapping.

CONCLUSIONS

Our findings identify potential roles of regulatory miRNA, intergenic non-coding DNA, and intronic non-coding RNA in the biology of ischemic stroke. These findings reveal new molecular targets that promise to help close the current gaps in accurate African ancestry-based genetic stroke's risk prediction and development of new targeted interventions to prevent or treat stroke.

Whole-Exome sequencing analysis identified TMSB10/TRABD2A locus to be associated with carfilzomib-related cardiotoxicity among patients with multiple myeloma

Background

Proteasome inhibitor Carfilzomib (CFZ) is effective in treating patients with refractory or relapsed multiple myeloma (MM) but has been associated with cardiovascular adverse events (CVAE) such as hypertension, cardiomyopathy, and heart failure. This study aimed to investigate the contribution of germline genetic variants in protein-coding genes in CFZ-CVAE among MM patients using whole-exome sequencing (WES) analysis.

Methods

Exome-wide single-variant association analysis, gene-based analysis, and rare variant analyses were performed on 603,920 variants in 247 patients with MM who have been treated with CFZ and enrolled in the Oncology Research Information Exchange Network (ORIEN) at the Moffitt Cancer Center. Separate analyses were performed in European Americans and African Americans followed by a trans-ethnic meta-analysis.

Results

The most significant variant in the exome-wide single variant analysis was a missense variant rs7148 in the thymosin beta-10/TraB Domain Containing 2A (TMSB10/TRABD2A) locus. The effect allele of rs7148 was associated with a higher risk of CVAE [odds ratio (OR) = 9.3 with a 95% confidence interval of 3.9-22.3, p = 5.42*10^-7]. MM patients with rs7148 AG or AA genotype had a higher risk of CVAE (50%) than those with GG genotype (10%). rs7148 is an expression quantitative trait locus (eQTL) for TRABD2A and TMSB10. The gene-based analysis also showed TRABD2A as the most significant gene associated with CFZ-CVAE (p = 1.06*10^-6).

Conclusions

We identified a missense SNP rs7148 in the TMSB10/TRABD2A as associated with CFZ-CVAE in MM patients. More investigation is needed to understand the underlying mechanisms of these associations.

Validation of Polygenic Risk Scores for Coronary Heart Disease in a Middle Eastern Cohort Using Whole Genome Sequencing

BACKGROUND

Enthusiasm for using polygenic risk scores (PRSs) in clinical practice is tempered by concerns about their portability to diverse ancestry groups, thus motivating genome-wide association studies in non-European ancestry cohorts.

METHODS

We conducted a genome-wide association study for coronary heart disease in a Middle Eastern cohort using whole genome sequencing and assessed the performance of 6 PRSs developed with methods including LDpred (PGS000296), metaGRS (PGS000018), Pruning and Thresholding (PGS000337), and an EnsemblePRS we developed. Additionally, we evaluated the burden of rare variants in lipid genes in cases and controls. Whole genome sequencing at 30× coverage was performed in 1067 coronary heart disease cases (mean age=59 years; 70.3% males) and 6170 controls (mean age=40 years; 43.5% males).

RESULTS

The majority of PRSs performed well; odds ratio (OR) per 1 SD increase (OR1sd) was highest for PGS000337 (OR1sd=1.81, 95% CI [1.66-1.98], P=3.07×10-41). EnsemblePRS performed better than individual PRSs (OR1sd=1.8, 95% CI [1.66-1.96], P=5.89×10-44). The OR for the 10th decile versus the remaining deciles was >3.2 for PGS000337, PGS000296, PGS000018, and reached 4.58 for EnsemblePRS. Of 400 known genome-wide significant loci, 33 replicated at P<10-4. However, the 9p21 locus did not replicate. Six suggestive (P<10-5) new loci/genes with plausible biological function were identified (eg, CORO7, RBM47, PDE4D). The burden of rare functional variants in LDLR, APOB, PCSK9, and ANGPTL4 was greater in cases than controls.

CONCLUSIONS

Overall, we demonstrate that PRSs derived from European ancestry genome-wide association studies performed well in a Middle Eastern cohort, suggesting these could be used in the clinical setting while ancestry-specific PRSs are developed.

Genome-wide association analyses of common infections in a large practice-based biobank

INTRODUCTION

Infectious diseases are common causes of morbidity and mortality worldwide. Susceptibility to infection is highly heritable; however, little has been done to identify the genetic determinants underlying common infectious diseases. One GWAS was performed using 23andMe information about self-reported infections; we set out to confirm previous loci and identify new ones using medically diagnosed infections.

METHODS

We used the electronic health record (EHR)-based biobank at Vanderbilt and diagnosis codes to identify cases of 12 infectious diseases in white patients: urinary tract infection, pneumonia, chronic sinus infections, otitis media, candidiasis, streptococcal pharyngitis, herpes zoster, herpes labialis, hepatitis B, infectious mononucleosis, tuberculosis (TB) or a positive TB test, and hepatitis C. We selected controls from patients with no diagnosis code for the candidate disease and matched by year of birth, sex, and calendar year at first and last EHR visits. We conducted GWAS using SAIGE and transcriptome-wide analysis (TWAS) using S-PrediXcan. We also conducted phenome-wide association study to understand associations between identified genetic variants and clinical phenotypes.

RESULTS

We replicated three 23andMe loci (p ≤ 0.05): herpes zoster and rs7047299-A (p = 2.6 × 10^-3) and rs2808290-C (p = 9.6 × 10^-3;); otitis media and rs114947103-C (p = 0.04). We also identified 2 novel regions (p ≤ 5 × 10^-8): rs113235453-G for otitis media (p = 3.04 × 10^-8), and rs10422015-T for candidiasis (p = 3.11 × 10^-8). In TWAS, four gene-disease associations were significant: SLC30A9 for otitis media (p = 8.06 × 10^-7); LRP3 and WDR88 for candidiasis (p = 3.91 × 10^-7 and p = 1.95 × 10^-6); and AAMDC for hepatitis B (p = 1.51 × 10^-6).

CONCLUSION

We conducted GWAS and TWAS for 12 infectious diseases and identified novel genetic contributors to the susceptibility of infectious diseases.

Associations between the rs5498 (A > G) and rs281432 (C > G) polymorphisms of the ICAM1 gene and atherosclerotic cardiovascular disease risk, including hypercholesterolemia

Background

Atherosclerotic cardiovascular disease (ASCVD) originates from complex risk factors, including age, gender, dyslipidemia, obesity, race, genetic and genetic variation. ICAM1 gene polymorphisms are a significant risk factor for ASCVD. However, the impact of the rs5498 and rs281432 polymorphisms on the prevalence of hypercholesterolemia (HCL) has not been reported. Therefore, we determine the relationships between single nucleotide polymorphisms (SNPs), including rs5498 and rs281432 on Intercellular adhesion molecule 1 gene (ICAM1) and ASCVD susceptibility in patients with HCL.

Methods

The clinical characteristics of 278 participants were assessed, and classified to groups having HCL and without HCL. ICAM1 SNPs genotyping was performed by DNA sequencing, and ICAM1 expression was measured using real-time PCR.

Results

Positive dominant model rs5498 participants had twice the risk of HCL (95% confidence interval (CI): [1.24-3.23], P = 0.005). The frequency of the G allele in rs5498 was 1.69 times higher in participants with HCL than in controls (95% CI [1.15-2.47], P = 0.007). Participants with the rs5498 AG or GG variants and high ICAM1 mRNA expression (≥3.12) had 2.49 times the risk (95% CI [1.42-4.38], P = 0.001), and those with a high LDL-C concentration (≥3.36 mmol/L) had 2.09 times the risk (95% CI [1.19-3.66], P = 0.010) of developing ASCVD compared with those with low ICAM1 mRNA and LDL-C levels. Interestingly, participants carrying the rs5498 AG or GG variants who had tachycardia (resting heart rates (RHRs) >100 beats/min) had a 5.02-times higher risk than those with a lower RHR (95% CI [1.35-18.63], P = 0.016).

Conclusions

It may consider the G allele in ICAM1 rs5498 is associated with a higher risk of ASCVD in Thai people with HCL, and is also positively associated with ICAM1 mRNA expression, LDL-C concentration, and RHR.

Transcriptomic and Lipidomic Mapping of Macrophages in the Hub of Chronic Beta-Adrenergic-Stimulation Unravels Hypertrophy-, Proliferation-, and Lipid Metabolism-Related Genes as Novel Potential Markers of Early Hypertrophy or Heart Failure

Sympathetic nervous system overdrive with chronic release of catecholamines is the most important neurohormonal mechanism activated to maintain cardiac output in response to heart stress. Beta-adrenergic signaling behaves first as a compensatory pathway improving cardiac contractility and maladaptive remodeling but becomes dysfunctional leading to pathological hypertrophy and heart failure (HF). Cardiac remodeling is a complex inflammatory syndrome where macrophages play a determinant role. This study aimed at characterizing the temporal transcriptomic evolution of cardiac macrophages in mice subjected to beta-adrenergic-stimulation using RNA sequencing. Owing to a comprehensive bibliographic analysis and complementary lipidomic experiments, this study deciphers typical gene profiles in early compensated hypertrophy (ECH) versus late dilated remodeling related to HF. We uncover cardiac hypertrophy- and proliferation-related transcription programs typical of ECH or HF macrophages and identify lipid metabolism-associated and Na⁺ or K⁺ channel-related genes as markers of ECH and HF macrophages, respectively. In addition, our results substantiate the key time-dependent role of inflammatory, metabolic, and functional gene regulation in macrophages during beta-adrenergic dependent remodeling. This study provides important and novel knowledge to better understand the prevalent key role of resident macrophages in response to chronically activated beta-adrenergic signaling, an effective diagnostic and therapeutic target in failing hearts.

Haplotype-resolved diverse human genomes and integrated analysis of structural variation

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

Minor allele of GJA1 gene polymorphism is associated with higher heart rate during atrial fibrillation

Atrial fibrillation (AF) tachycardia causes heart failure and requires more attention. The genetic background of individual heart rate (HR) variations during AF are unclear. We hypothesized that HR-associated single nucleotide polymorphisms (SNPs) reported in Genome-Wide Association Studies (GWAS) are also associated with HR during AF. We enrolled patients with persistent AF (311 for screening and 146 for replication) who underwent AF ablation and were genotyped for the 21 h-associated SNPs reported in GWAS. The patients underwent 24-h Holter monitoring before AF ablation and electrophysiological study after AF ablation during sinus rhythm. Only the GJA1 SNP rs1015451 (T>C) was significantly associated with total HR (TT 110,643 ± 17,542 beats/day, TC 116,350 ± 19,060 beats/day, CC 122,163 ± 25,684 beats/day, P = 8.5 × 10⁻⁴). We also confirmed this significant association in the replication set. The intra-atrial conduction was faster in AF patients with the GJA1 minor allele than in those without it. Multivariate analysis revealed the presence of a GJA1 SNP rs1015451 additive model, female gender, lower left ventricular ejection fraction, and higher 1:1 atrioventricular nodal conduction were independently associated with higher HR during AF. The GJA1 SNP might be a new genetic marker for AF tachycardia.