Search

The association of the cystic fibrosis transmembrane conductance regulator (CFTR) and epithelial sodium channel (ENaC) in the pathophysiology of cystic fibrosis (CF) is controversial. Previously, we demonstrated a close physical association between wild-type (WT) CFTR and WT ENaC. We have also shown that the F508del CFTR fails to associate with ENaC unless the mutant protein is rescued pharmacologically or by low temperature. In this study, we present the evidence for a direct physical association between WT CFTR and ENaC subunits carrying Liddle’s syndrome mutations. We show that all three ENaC subunits bearing Liddle’s syndrome mutations (both point mutations and the complete truncation of the carboxy terminus), could be coimmunoprecipitated with WT CFTR. The biochemical studies were complemented by fluorescence lifetime imaging microscopy (FLIM), a distance-dependent approach that monitors protein-protein interactions between fluorescently labeled molecules. Our measurements revealed significantly increased fluorescence resonance energy transfer between CFTR and all tested ENaC combinations as compared with controls (ECFP and EYFP cotransfected cells). Our findings are consistent with the notion that CFTR and ENaC are within reach of each other even in the setting of Liddle’s syndrome mutations, suggestive of a direct intermolecular interaction between these two proteins.

We describe the protocol for identifying COVID-19 severity specific cell types and their regulatory marker genes using single-cell transcriptomics data. We construct COVID-19 comorbid disease-associated gene list using multiple databases and literature resources. Next, we identify specific cell type where comorbid genes are upregulated. We further characterize the identified cell type using gene enrichment analysis. We detect upregulation of marker gene restricted to severe COVID-19 cell type and validate our findings using

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.

The advent of long-read sequencing offers a new assessment method of detecting genomic structural variation (SV) in numerous rare genetic diseases. For autism spectrum disorders (ASD) cases where pathogenic variants fail to be found in the protein-coding genic regions along chromosomes, we proposed a scalable workflow to characterize the risk factor of SVs impacting non-coding elements of the genome. We applied whole-genome sequencing on an Emirati family having three children with ASD using long and short-read sequencing technology. A series of analytical pipelines were established to identify a set of SVs with high sensitivity and specificity. At 15-fold coverage, we observed that long-read sequencing technology (987 variants) detected a significantly higher number of SVs when compared to variants detected using short-read technology (509 variants) (p-value &lt, 1.1020 × 10−57). Further comparison showed 97.9% of long-read sequencing variants were spanning within the 1–100 kb size range (p-value &lt, 9.080 × 10−67) and impacting over 5000 genes. Moreover, long-read variants detected 604 non-coding RNAs (p-value &lt, 9.02 × 10−9), comprising 58% microRNA, 31.9% lncRNA, and 9.1% snoRNA. Even at low coverage, long-read sequencing has shown to be a reliable technology in detecting SVs impacting complex elements of the genome.

In order to investigate the early phase of the amyloid formation by the short amyloidogenic octapeptide sequence (‘NFGAILSS’) derived from IAPP, we carried out a 100ns all-atom molecular dynamics (MD) simulations of systems that contain 27 peptides and over 30,000 water molecules. The large-scale calculations were performed for the wild type sequence and seven alanine-scanned sequences using AMBER 8.0 on RIKEN’s special purpose MD-GRAPE3 supercomputer, using the all-atom point charge force field ff99, which do not favor β-structures. Large peptide clusters (size 18–26 mers) were observed for all simulations, and our calculations indicated that isoleucine at position 5 played important role in the formation of β-rich clusters. In the oligomeric state, the wild type and the S7A sequences had the highest β-structure content (~14%), as calculated by DSSP, in line with experimental observations, whereas I5A and G3A had the highest helical content (~20%). Importantly, the β-structure preferences of wild type IAPP originate from its association into clusters and are not intrinsic to its sequence. Altogether, the results of this first large-scale, multi-peptide all-atom molecular dynamics simulation appear to provide insights into the mechanism of amyloidogenic and non-amyloidogenic oligomers that mainly corroborate previous experimental observations.