Characterization of the novel heterozygous SCN5A genetic variant Y739D associated with Brugada syndrome.

Genetic variants in SCN5A gene were identified in patients with various arrhythmogenic conditions including Brugada syndrome. Despite significant progress of last decades in studying the molecular mechanism of arrhythmia-associated SCN5A mutations, the understanding of relationship between genetics, electrophysiological consequences and clinical phenotype is lacking. We have found a novel genetic variant Y739D in the SCN5A-encoded sodium channel Nav1.5 of a male patient with Brugada syndrome (BrS). The objective of the study was to characterize the biophysical properties of Nav1.5-Y739D and provide possible explanation of the phenotype observed in the patient. The WT and Y739D channels were heterologously expressed in the HEK-293T cells and the whole-cell sodium currents were recorded. Substitution Y739D reduced the sodium current density by 47 ± 2% at -20 mV, positively shifted voltage-dependent activation, accelerated both fast and slow inactivation, and decelerated recovery from the slow inactivation. The Y739D loss-of-function phenotype likely causes the BrS manifestation. In the hNav1.5 homology models, which are based on the cryo-EM structure of rat Nav1.5 channel, Y739 in the extracellular loop IIS1-S2 forms H-bonds with K1381 and E1435 and pi-cation contacts with K1397 (all in loop IIIS5-P1). In contrast, Y739D accepts H-bonds from K1397 and Y1434. Substantially different contacts of Y739 and Y739D with loop IIIS5-P1 would differently transmit allosteric signals from VSD-II to the fast-inactivation gate at the N-end of helix IIIS5 and slow-inactivation gate at the C-end of helix IIIP1. This may underlie the atomic mechanism of the Y739D channel dysfunction.

The Activity of KIF14, Mieap, and EZR in a New Type of the Invasive Component, Torpedo-Like Structures, Predetermines the Metastatic Potential of Breast Cancer.

Intratumor morphological heterogeneity reflects patterns of invasive growth and is an indicator of the metastatic potential of breast cancer. In this study, we used this heterogeneity to identify molecules associated with breast cancer invasion and metastasis. The gene expression microarray data were used to identify genes differentially expressed between solid, trabecular, and other morphological arrangements of tumor cells. Immunohistochemistry was applied to evaluate the association of the selected proteins with metastasis. RNA-sequencing was performed to analyze the molecular makeup of metastatic tumor cells. High frequency of metastases and decreased metastasis-free survival were detected in patients either with positive expression of KIF14 or Mieap or negative expression of EZR at the tips of the torpedo-like structures in breast cancers. KIF14- and Mieap-positive and EZR-negative cells were mainly detected in the torpedo-like structures of the same breast tumors; however, their transcriptomic features differed. KIF14-positive cells showed a significant upregulation of genes involved in ether lipid metabolism. Mieap-positive cells were enriched in genes involved in mitophagy. EZR-negative cells displayed upregulated genes associated with phagocytosis and the chemokine-mediated signaling pathway. In conclusion, the positive expression of KIF14 and Mieap and negative expression of EZR at the tips of the torpedo-like structures are associated with breast cancer metastasis.

Systematic dissection of biases in whole-exome and whole-genome sequencing reveals major determinants of coding sequence coverage.

Advantages and diagnostic effectiveness of the two most widely used resequencing approaches, whole exome (WES) and whole genome (WGS) sequencing, are often debated. WES dominated large-scale resequencing projects because of lower cost and easier data storage and processing. Rapid development of 3rd generation sequencing methods and novel exome sequencing kits predicate the need for a robust statistical framework allowing informative and easy performance comparison of the emerging methods. In our study we developed a set of statistical tools to systematically assess coverage of coding regions provided by several modern WES platforms, as well as PCR-free WGS. We identified a substantial problem in most previously published comparisons which did not account for mappability limitations of short reads. Using regression analysis and simple machine learning, as well as several novel metrics of coverage evenness, we analyzed the contribution from the major determinants of CDS coverage. Contrary to a common view, most of the observed bias in modern WES stems from mappability limitations of short reads and exome probe design rather than sequence composition. We also identified the ~ 500 kb region of human exome that could not be effectively characterized using short read technology and should receive special attention during variant analysis. Using our novel metrics of sequencing coverage, we identified main determinants of WES and WGS performance. Overall, our study points out avenues for improvement of enrichment-based methods and development of novel approaches that would maximize variant discovery at optimal cost.

Characterization of a novel SCN5A genetic variant A1294G associated with mixed clinical phenotype.

Mutations in gene SCN5A, which encodes cardiac voltage-gated sodium channel Nav1.5, are associated with multiple clinical phenotypes. Here we describe a novel A1294G genetic variant detected in a male patient with combined clinical phenotype including atrioventricular II block, Brugada-like ECG, septal fibrosis, right ventricular dilatation and decreased left ventricular contractility. Residue A1294 is located in the IIIS3-S4 extracellular loop, in proximity to several residues whose mutations are associated with sodium channelopathies. The wild-type channel Nav1.5 and mutant Nav1.5-A1294G were expressed in the CHO-K1 and HEK293T cells and whole-cell sodium currents were recorded using the patch-clamp method. The A1294G channels demonstrated a negative shift of steady-state inactivation, accelerated fast and slow inactivation and decelerated recovery from intermediate inactivation. Our study reveals biophysical mechanism of the Nav1.5-A1294G dysfunction, which may underlie the combined phenotypic manifestation observed in the patient.

Atomic Mechanisms of Timothy Syndrome-Associated Mutations in Calcium Channel Cav1.2.

Timothy syndrome (TS) is a very rare multisystem disorder almost exclusively associated with mutations G402S and G406R in helix IS6 of Cav1.2. Recently, mutations R518C/H in helix IIS0 of the voltage sensing domain II (VSD-II) were described as a cause of cardiac-only TS. The three mutations are known to decelerate voltage-dependent inactivation (VDI). Here, we report a case of cardiac-only TS caused by mutation R518C. To explore possible impact of the three mutations on interdomain contacts, we modeled channel Cav1.2 using as templates Class Ia and Class II cryo-EM structures of presumably inactivated channel Cav1.1. In both models, R518 and several other residues in VSD-II donated H-bonds to the IS6-linked α1-interaction domain (AID). We further employed steered Monte Carlo energy minimizations to move helices S4-S5, S5, and S6 from the inactivated-state positions to those seen in the X-ray structures of the open and closed NavAb channel. In the open-state models, positions of AID and VSD-II were similar to those in Cav1.1. In the closed-state models, AID moved along the ß subunit (Cavß) toward the pore axis and shifted AID-bound VSD-II. In all the models R518 retained strong contacts with AID. Our calculations suggest that conformational changes in VSD-II upon its deactivation would shift AID along Cavß toward the pore axis. The AID-linked IS6 would bend at flexible G402 and G406, facilitating the activation gate closure. Mutations R518C/H weakened the IIS0-AID contacts and would retard the AID shift. Mutations G406R and G402S stabilized the open state and would resist the pore closure. Several Cav1.2 mutations associated with long QT syndromes are consistent with this proposition. Our results provide a mechanistic rationale for the VDI deceleration caused by TS-associated mutations and suggest targets for further studies of calcium channelopathies.

The exon junction complex factor Y14 is dynamic in the nucleus of the beetle Tribolium castaneum during late oogenesis.

BACKGROUND: The oocyte chromosomes of the red flour beetle, Tribolium castaneum, are gathered into a knot, forming a karyosphere at the diplotene stage of meiotic prophase. Chromatin rearrangement, which is a characteristic feature of oocyte maturation, is well documented. The T. castaneum karyosphere is surrounded by a complex extrachromosomal structure termed the karyosphere capsule. The capsule contains the vast majority of oocyte RNA. We have previously shown using a BrUTP assay that oocyte chromosomes in T. castaneum maintain residual transcription up to the very end of oocyte maturation. Karyosphere transcription requires evidently not only transcription factors but also mRNA processing factors, including the components of the exon junction complex with its core component, the splicing factor Y14. We employed a gene engineering approach with injection of mRNA derived from the Myc-tagged Y14 plasmid-based construct in order to monitor the newly synthesized fusion protein in the oocyte nuclei. RESULTS: Our preliminary data have been presented as a brief correspondence elsewhere. Here, we provide a full-length article including immunoelectron-microscopy localization data on Y14-Myc distribution in the nucleus of previtellogenic and vitellogenic oocytes. The injections of the fusion protein Y14-Myc mRNA into the oocytes showed a dynamic pattern of the protein distribution. At the previtellogenic stage, there are two main locations for the protein: SC35 domains (the analogues of interchromatin granule clusters or nuclear speckles) and the karyosphere capsule. At the vitellogenic stage, SC35 domains were devoid of labels, and Y14-Myc was found in the perichromatin region of the karyosphere, presumably at the places of residual transcription. We show that karyosphere formation is accompanied by the movement of a nuclear protein while the residual transcription occurs during genome inactivation. CONCLUSIONS: Our data indicate that the karyosphere capsule, being a destination site for a protein involved in mRNA splicing and export, is not only a specializes part of nuclear matrix separating the karyosphere from the products of chromosome activity, as believed previously, but represents a special nuclear compartment involved in the processes of gene expression in the case the karyosphere retains residual transcription activity.

Chicken rRNA Gene Cluster Structure.

Ribosomal RNA (rRNA) genes, whose activity results in nucleolus formation, constitute an extremely important part of genome. Despite the extensive exploration into avian genomes, no complete description of avian rRNA gene primary structure has been offered so far. We publish a complete chicken rRNA gene cluster sequence here, including 5'ETS (1836 bp), 18S rRNA gene (1823 bp), ITS1 (2530 bp), 5.8S rRNA gene (157 bp), ITS2 (733 bp), 28S rRNA gene (4441 bp) and 3'ETS (343 bp). The rRNA gene cluster sequence of 11863 bp was assembled from raw reads and deposited to GenBank under KT445934 accession number. The assembly was validated through in situ fluorescent hybridization analysis on chicken metaphase chromosomes using computed and synthesized specific probes, as well as through the reference assembly against de novo assembled rRNA gene cluster sequence using sequenced fragments of BAC-clone containing chicken NOR (nucleolus organizer region). The results have confirmed the chicken rRNA gene cluster validity.