Single-cell reconstruction and mutation enrichment analysis identifies dysregulated cardiomyocyte and endothelial cells in congenital heart disease.

Congenital heart disease (CHD) is one of the most prevalent neonatal congenital anomalies. To catalog the putative candidate CHD risk genes, we collected 16,349 variants [single-nucleotide variants (SNVs) and Indels] impacting 8,308 genes in 3,166 CHD cases for a comprehensive meta-analysis. Using American College of Medical Genetics (ACMG) guidelines, we excluded the 0.1% of benign/likely benign variants and the resulting dataset consisted of 83% predicted loss of function variants and 17% missense variants. Seventeen percent were de novo variants. A stepwise analysis identified 90 variant-enriched CHD genes, of which six (GPATCH1, NYNRIN, TCLD2, CEP95, MAP3K19, and TTC36) were novel candidate CHD genes. Single-cell transcriptome cluster reconstruction analysis on six CHD tissues and four controls revealed upregulation of the top 10 frequently mutated genes primarily in cardiomyocytes. NOTCH1 (highest number of variants) and MYH6 (highest number of recurrent variants) expression was elevated in endocardial cells and cardiomyocytes, respectively, and 60% of these gene variants were associated with tetralogy of Fallot and coarctation of the aorta, respectively. Pseudobulk analysis using the single-cell transcriptome revealed significant (P < 0.05) upregulation of both NOTCH1 (endocardial cells) and MYH6 (cardiomyocytes) in the control heart data. We observed nine different subpopulations of CHD heart cardiomyocytes of which only four were observed in the control heart. This is the first comprehensive meta-analysis combining genomics and CHD single-cell transcriptomics, identifying the most frequently mutated CHD genes, and demonstrating CHD gene heterogeneity, suggesting that multiple genes contribute to the phenotypic heterogeneity of CHD. Cardiomyocytes and endocardial cells are identified as major CHD-related cell types.NEW & NOTEWORTHY Congential heart disease (CHD) is one of the most prevalent neonatal congenital anomalies. We present a comprehensive analysis combining genomics and CHD single-cell transcriptome. Our study identifies 90 potential candidate CHD risk genes of which 6 are novel. The risk genes have heterogenous expression suggestive of multiple genes contributing to the phenotypic heterogeneity of CHD. Cardiomyocytes and endocardial cells are identified as major CHD-related cell types.

Single-cell transcriptomics trajectory and molecular convergence of clinically relevant mutations in Brugada syndrome.

Brugada syndrome (BrS) is a rare, inherited arrhythmia with high risk of sudden cardiac death. To evaluate the molecular convergence of clinically relevant mutations and to identify developmental cardiac cell types that are associated with BrS etiology, we collected 733 mutations represented by 16 sodium, calcium, potassium channels, and regulatory and structural genes related to BrS. Among the clinically relevant mutations, 266 are unique singletons and 88 mutations are recurrent. We observed an over-representation of clinically relevant mutations (∼80%) in SCN5A gene and also identified several candidate genes, including GPD1L, TRPM4, and SCN10A. Furthermore, protein domain enrichment analysis revealed that a large proportion of the mutations impacted ion transport domains in multiple genes, including SCN5A, TRPM4, and SCN10A. A comparative protein domain analysis of SCN5A further established a significant (P = 0.04) enrichment of clinically relevant mutations within ion transport domain, including a significant (P = 0.02) mutation hotspot within 1321-1380 residue. The enrichment of clinically relevant mutations within SCN5A ion transport domain is stronger (P = 0.00003) among early onset of BrS. Our spatiotemporal cellular heart developmental (prenatal to adult) trajectory analysis applying single-cell transcriptome identified the most frequently BrS-mutated genes (SCN5A and GPD1L) are significantly upregulated in the prenatal cardiomyocytes. A more restrictive cellular expression trajectory is prominent in the adult heart ventricular cardiomyocytes compared to prenatal. Our study suggests that genomic and proteomic hotspots in BrS converge into ion transport pathway and cardiomyocyte as a major BrS-associated cell type that provides insight into the complex genetic etiology of BrS.NEW & NOTEWORTHY Brugada syndrome is a rare inherited arrhythmia with high risk of sudden cardiac death. We present the findings for a molecular convergence of clinically relevant mutations and identification of a single-cell transcriptome-derived cardiac cell types that are associated with the etiology of BrS. Our study suggests that genomic and proteomic hotspots in BrS converge into ion transport pathway and cardiomyocyte as a major BrS-associated cell type that provides insight into the complex genetic etiology of BrS.

SARS-CoV-2 may hijack GPCR signaling pathways to dysregulate lung ion and fluid transport.

The tropism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a virus responsible for the ongoing coronavirus disease 2019 (COVID-19) pandemic, toward the host cells is determined, at least in part, by the expression and distribution of its cell surface receptor, angiotensin-converting enzyme 2 (ACE2). The virus further exploits the host cellular machinery to gain access into the cells; its spike protein is cleaved by a host cell surface transmembrane serine protease 2 (TMPRSS2) shortly after binding ACE2, followed by its proteolytic activation at a furin cleavage site. The virus primarily targets the epithelium of the respiratory tract, which is covered by a tightly regulated airway surface liquid (ASL) layer that serves as a primary defense mechanism against respiratory pathogens. The volume and viscosity of this fluid layer is regulated and maintained by a coordinated function of different transport pathways in the respiratory epithelium. We argue that SARS-CoV-2 may potentially alter evolutionary conserved second-messenger signaling cascades via activation of G protein-coupled receptors (GPCRs) or by directly modulating G protein signaling. Such signaling may in turn adversely modulate transepithelial transport processes, especially those involving cystic fibrosis transmembrane conductance regulator (CFTR) and epithelial Na+ channel (ENaC), thereby shifting the delicate balance between anion secretion and sodium absorption, which controls homeostasis of this fluid layer. As a result, activation of the secretory pathways including CFTR-mediated Cl- transport may overwhelm the absorptive pathways, such as ENaC-dependent Na+ uptake, and initiate a pathophysiological cascade leading to lung edema, one of the most serious and potentially deadly clinical manifestations of COVID-19.