Melanocortin Derivatives Induced Vascularization and Neuroglial Proliferation in the Rat Brain under Conditions of Cerebral Ischemia.

Stroke remains the second leading cause of death worldwide. The development of new therapeutic agents focused on restoring vascular function and neuroprotection of viable tissues is required. In this study the neuroprotective activity of melanocortin-like ACTH(4-7)PGP and ACTH(6-9)PGP peptides was investigated in rat brain at 24 h after transient middle cerebral artery occlusion (tMCAO). The severity of ischemic damage, changes in the proliferative activity of neuroglial cells and vascularization of rat brain tissue were analyzed. The administration of peptides resulted in a significant increase in the volume density of neurons in the perifocal zone of infarction compared to rats subjected to ischemia and receiving saline. Immunohistochemical analysis of the proliferative activity of neuroglia cells using PCNA antibodies showed a significant increase in the number of proliferating cells in the penumbra and in the intact cerebral cortex of rats receiving peptide treatment. The effect of peptides on vascularization was examined using CD31 antibodies under tMCAO conditions, revealing a significant increase in the volume density of vessels and their sizes in the penumbra after administration of ACTH(4-7)PGP and ACTH(6-9)PGP. These findings confirm the neuroprotective effect of peptides due to the activation of neuroglia proliferation and the enhancement of collateral blood flow.

Circular RNAs Variously Participate in Coronary Atherogenesis.

Over the past decade, numerous studies have shown that circular RNAs (circRNAs) play a significant role in coronary artery atherogenesis and other cardiovascular diseases. They belong to the class of non-coding RNAs and arise as a result of non-canonical splicing of premature RNA, which results in the formation of closed single-stranded circRNA molecules that lack 5'-end caps and 3'-end poly(A) tails. circRNAs have broad post-transcriptional regulatory activity. Acting as a sponge for miRNAs, circRNAs compete with mRNAs for binding to miRNAs, acting as competing endogenous RNAs. Numerous circRNAs are involved in the circRNA-miRNA-mRNA regulatory axes associated with the pathogenesis of cardiomyopathy, chronic heart failure, hypertension, atherosclerosis, and coronary artery disease. Recent studies have shown that сirc_0001445, circ_0000345, circ_0093887, сircSmoc1-2, and circ_0003423 are involved in the pathogenesis of coronary artery disease (CAD) with an atheroprotective effect, while circ_0002984, circ_0029589, circ_0124644, circ_0091822, and circ_0050486 possess a proatherogenic effect. With their high resistance to endonucleases, circRNAs are promising diagnostic biomarkers and therapeutic targets. This review aims to provide updated information on the involvement of atherogenesis-related circRNAs in the pathogenesis of CAD. We also discuss the main modern approaches to detecting and studying circRNA-miRNA-mRNA interactions, as well as the prospects for using circRNAs as biomarkers and therapeutic targets for the treatment of cardiovascular diseases.

Differential Expression of Subsets of Genes Related to HDL Metabolism and Atherogenesis in the Peripheral Blood in Coronary Artery Disease.

Differential expression of genes (DEGs) in coronary artery disease (CAD) and the association between transcript level and high-density lipoprotein cholesterol (HDL-C) were studied with 76 male patients with CAD and 63 control patients. The transcript level of genes related to HDL metabolism (24 genes) and atherosclerosis-prone (41 genes) in RNA isolated from peripheral blood mononuclear cells was measured by real-time RT-PCR. Twenty-eight DEGs were identified. The expression of cholesterol transporters, ALB, APOA1, and LCAT was down-regulated, while the expression of AMN, APOE, LDLR, LPL, PLTP, PRKACA, and CETP was up-regulated. The systemic inflammation in CAD is evidenced by the up-regulation of IL1B, TLR8, CXCL5, and TNFRSF1A. For the controls, TLR8 and SOAT1 were negative predictors of the HDL-C level. For CAD patients, PRKACG, PRKCQ, and SREBF1 were positive predictors, while PRKACB, LCAT, and S100A8 were negative predictors. For CAD patients, the efficiency of reverse cholesterol transport is 73-79%, and intracellular free cholesterol seems to accumulate at hyperalphalipoproteinemia. Both atheroprotective (via S100A8) and proatherogenic (via SREBF1, LCAT, PRKACG, PRKACB, and PRKCQ) associations of gene expression with HDL-C determine HDL functionality in CAD patients. The selected key genes and involved pathways may represent HDL-specific targets for the diagnosis and treatment of CAD and atherosclerosis.

Evaluation of Density-Based Spatial Clustering for Identifying Genomic Loci Associated with Ischemic Stroke in Genome-Wide Data.

The genetic architecture of ischemic stroke (IS), which is one of the leading causes of death worldwide, is complex and underexplored. The traditional approach for associative gene mapping is genome-wide association studies (GWASs), testing individual single-nucleotide polymorphisms (SNPs) across the genomes of case and control groups. The purpose of this research is to develop an alternative approach in which groups of SNPs are examined rather than individual ones. We proposed, validated and applied to real data a new workflow consisting of three key stages: grouping SNPs in clusters, inferring the haplotypes in the clusters and testing haplotypes for the association with phenotype. To group SNPs, we applied the clustering algorithms DBSCAN and HDBSCAN to linkage disequilibrium (LD) matrices, representing pairwise r2 values between all genotyped SNPs. These clustering algorithms have never before been applied to genotype data as part of the workflow of associative studies. In total, 883,908 SNPs and insertion/deletion polymorphisms from people of European ancestry (4929 cases and 652 controls) were processed. The subsequent testing for frequencies of haplotypes restored in the clusters of SNPs revealed dozens of genes associated with IS and suggested the complex role that protocadherin molecules play in IS. The developed workflow was validated with the use of a simulated dataset of similar ancestry and the same sample sizes. The results of classic GWASs are also provided and discussed. The considered clustering algorithms can be applied to genotypic data to identify the genomic loci associated with different qualitative traits, using the workflow presented in this research.

The Study of the Association of Polymorphisms in LSP1, GPNMB, PDPN, TAGLN, TSPO, and TUBB6 Genes with the Risk and Outcome of Ischemic Stroke in the Russian Population.

To date, there has been great progress in understanding the genetic basis of ischemic stroke (IS); however, several aspects of the condition remain underexplored, including the influence of genetic factors on post-stroke outcomes and the identification of causative loci. We proposed that an analysis of the results obtained from animal models of brain ischemia could be helpful. To this end, we developed a bioinformatic approach for exploring single-nucleotide polymorphisms (SNPs) in human orthologs of rat genes expressed differentially after induced brain ischemia. Using this approach, we identified and analyzed 11 SNPs from 6 genes in 553 Russian individuals (331 patients with IS and 222 controls). We assessed the association of SNPs with the risk of IS and IS outcomes. We found that the SNPs rs858239 (GPNMB), rs907611 (LSP1), and rs494356 (TAGLN) were associated with different parameters of IS functional outcomes. In addition, the SNP rs1261025 (PDPN) was associated significantly with IS itself (p = 0.0188, recessive model). All these associations were demonstrated for the first time. Analysis of the literature suggests that they should be characterized as being inflammation related. This supports the pivotal role of inflammation in both the incidence of stroke and post-stroke outcomes. We believe the findings reported here will help with stroke prognosis in the future.

CAPN3 c.1746-20C>G variant is hypomorphic for LGMD R1 calpain 3-related.

The investigated intronic CAPN3 variant NM_000070.3:c.1746-20C>G occurs in the Central and Eastern Europe with a frequency of >1% and there are conflicting interpretations on its pathogenicity. We collected data on 14 patients carrying the CAPN3 c.1746-20C>G variant in trans position with another CAPN3 pathogenic/likely pathogenic variant. The patients compound heterozygous for the CAPN3 c.1746-20C>G variant presented a phenotype consistent with calpainopathy of mild/medium severity. This variant is most frequent in the North/West regions of Russia and may originate from that area. Molecular studies revealed that different splicing isoforms are produced in the muscle. We hypothesize that c.1746-20C>G is a hypomorphic variant with a reduction of RNA and protein expression and only individuals having a higher ratio of abnormal isoforms are affected. Reclassification of the CAPN3 variant c.1746-20C>G from variant with a conflicting interpretation of pathogenicity to hypomorphic variant explains many unidentified cases of limb girdle muscular dystrophy R1 calpain 3-related in Eastern and Central Europe.

Comparative Use of Contralateral and Sham-Operated Controls Reveals Traces of a Bilateral Genetic Response in the Rat Brain after Focal Stroke.

Ischemic stroke is a multifactorial disease with a complex etiology and global consequences. Model animals are widely used in stroke studies. Various controls, either brain samples from sham-operated (SO) animals or symmetrically located brain samples from the opposite (contralateral) hemisphere (CH), are often used to analyze the processes in the damaged (ipsilateral) hemisphere (IH) after focal stroke. However, previously, it was shown that focal ischemia can lead to metabolic and transcriptomic changes not only in the IH but also in the CH. Here, using a transient middle cerebral artery occlusion (tMCAO) model and genome-wide RNA sequencing, we identified 1941 overlapping differentially expressed genes (DEGs) with a cutoff value >1.5 and Padj < 0.05 that reflected the general transcriptome response of IH subcortical cells at 24 h after tMCAO using both SO and CH controls. Concomitantly, 861 genes were differentially expressed in IH vs. SO, whereas they were not vs. the CH control. Furthermore, they were associated with apoptosis, the cell cycle, and neurotransmitter responses. In turn, we identified 221 DEGs in IH vs. CH, which were non-DEGs vs. the SO control. Moreover, they were predominantly associated with immune-related response. We believe that both sets of non-overlapping genes recorded transcriptome changes in IH cells associated with transhemispheric differences after focal cerebral ischemia. Thus, the specific response of the CH transcriptome should be considered when using it as a control in studies of target brain regions in diseases that induce a global bilateral genetic response, such as stroke.

Pharmacotranscriptomics of peptide drugs with neuroprotective properties.

Here we present a review of studies on the effects of peptides with neuroprotective properties on gene transcription in nerve cells. The few published works in this area clearly demonstrate massive changes in cell transcriptomes induced by peptides under normal conditions and under conditions of experimental brain ischemia. These changes significantly affect signaling and metabolic pathways, affecting various body systems and confirming the multiple target actions of peptides. The importance of noncoding RNAs in the regulation of these processes is shown, and we discuss the prospects of research for determining the main mechanisms of peptide regulation, which is necessary for the further development of drugs with targeted neuroprotective effects.

Genome-wide sequence analyses of ethnic populations across Russia.

The Russian Federation is the largest and one of the most ethnically diverse countries in the world, however no centralized reference database of genetic variation exists to date. Such data are crucial for medical genetics and essential for studying population history. The Genome Russia Project aims at filling this gap by performing whole genome sequencing and analysis of peoples of the Russian Federation. Here we report the characterization of genome-wide variation of 264 healthy adults, including 60 newly sequenced samples. People of Russia carry known and novel genetic variants of adaptive, clinical and functional consequence that in many cases show allele frequency divergence from neighboring populations. Population genetics analyses revealed six phylogeographic partitions among indigenous ethnicities corresponding to their geographic locales. This study presents a characterization of population-specific genomic variation in Russia with results important for medical genetics and for understanding the dynamic population history of the world's largest country.

Antistress Action of Melanocortin Derivatives Associated with Correction of Gene Expression Patterns in the Hippocampus of Male Rats Following Acute Stress.

Natural melanocortins (MCs) have been used in the successful development of drugs with neuroprotective properties. Here, we studied the behavioral effects and molecular genetic mechanisms of two synthetic MC derivatives-ACTH(4-7)PGP (Semax) and ACTH(6-9)PGP under normal and acute restraint stress (ARS) conditions. Administration of Semax or ACTH(6-9)PGP (100 µg/kg) to rats 30 min before ARS attenuated ARS-induced behavioral alterations. Using high-throughput RNA sequencing (RNA-Seq), we identified 1359 differentially expressed genes (DEGs) in the hippocampus of vehicle-treated rats subjected to ARS, using a cutoff of >1.5 fold change and adjusted p-value (Padj) < 0.05, in samples collected 4.5 h after the ARS. Semax administration produced > 1500 DEGs, whereas ACTH(6-9)PGP administration led to <400 DEGs at 4.5 h after ARS. Nevertheless, ~250 overlapping DEGs were identified, and expression of these DEGs was changed unidirectionally by both peptides under ARS conditions. Modulation of the expression of genes associated with biogenesis, translation of RNA, DNA replication, and immune and nervous system function was produced by both peptides. Furthermore, both peptides upregulated the expression levels of many genes that displayed decreased expression after ARS, and vice versa, the MC peptides downregulated the expression levels of genes that were upregulated by ARS. Consequently, the antistress action of MC peptides may be associated with a correction of gene expression patterns that are disrupted during ARS.

Brain Protein Expression Profile Confirms the Protective Effect of the ACTH(4-7)PGP Peptide (Semax) in a Rat Model of Cerebral Ischemia-Reperfusion.

The Semax (Met-Glu-His-Phe-Pro-Gly-Pro) peptide is a synthetic melanocortin derivative that is used in the treatment of ischemic stroke. Previously, studies of the molecular mechanisms underlying the actions of Semax using models of cerebral ischemia in rats showed that the peptide enhanced the transcription of neurotrophins and their receptors and modulated the expression of genes involved in the immune response. A genome-wide RNA-Seq analysis revealed that, in the rat transient middle cerebral artery occlusion (tMCAO) model, Semax suppressed the expression of inflammatory genes and activated the expression of neurotransmitter genes. Here, we aimed to evaluate the effect of Semax in this model via the brain expression profiling of key proteins involved in inflammation and cell death processes (MMP-9, c-Fos, and JNK), as well as neuroprotection and recovery (CREB) in stroke. At 24 h after tMCAO, we observed the upregulation of active CREB in subcortical structures, including the focus of the ischemic damage; downregulation of MMP-9 and c-Fos in the adjacent frontoparietal cortex; and downregulation of active JNK in both tissues under the action of Semax. Moreover, a regulatory network was constructed. In conclusion, the suppression of inflammatory and cell death processes and the activation of recovery may contribute to the neuroprotective action of Semax at both the transcriptome and protein levels.

Reconstructing genetic history of Siberian and Northeastern European populations.

Siberia and Northwestern Russia are home to over 40 culturally and linguistically diverse indigenous ethnic groups, yet genetic variation and histories of peoples from this region are largely uncharacterized. We present deep whole-genome sequencing data (∼38×) from 28 individuals belonging to 14 distinct indigenous populations from that region. We combined these data sets with additional 32 modern-day and 46 ancient human genomes to reconstruct genetic histories of several indigenous Northern Eurasian populations. We found that Siberian and East Asian populations shared 38% of their ancestry with a 45,000-yr-old Ust'-Ishim individual who was previously believed to have no modern-day descendants. Western Siberians trace 57% of their ancestry to ancient North Eurasians, represented by the 24,000-yr-old Siberian Mal'ta boy MA-1. Eastern Siberian populations formed a distinct sublineage that separated from other East Asian populations ∼10,000 yr ago. In addition, we uncovered admixtures between Siberians and Eastern European hunter-gatherers from Samara, Karelia, Hungary, and Sweden (from 8000-6600 yr ago); Yamnaya people (5300-4700 yr ago); and modern-day Northeastern Europeans. Our results provide new insights into genetic histories of Siberian and Northeastern European populations and evidence of ancient gene flow from Siberia into Europe.

Genome-wide transcriptome analysis using RNA-Seq reveals a large number of differentially expressed genes in a transient MCAO rat model.

BACKGROUND: The transient middle cerebral artery occlusion (tMCAO) model is used for studying the molecular mechanisms of ischemic damage and neuroprotection. Numerous studies have demonstrated the role of individual genes and associated signaling pathways in the pathogenesis of ischemic stroke. Here, the tMCAO model was used to investigate the genome-wide response of the transcriptome of rat brain tissues to the damaging effect of ischemia and subsequent reperfusion. RESULTS: Magnetic resonance imaging and histological examination showed that the model of focal ischemia based on endovascular occlusion of the right middle cerebral artery for 90 min using a monofilament, followed by restoration of the blood flow, led to reproducible localization of ischemic damage in the subcortical structures of the brain. High-throughput RNA sequencing (RNA-Seq) revealed the presence of differentially expressed genes (DEGs) in subcortical structures of rat brains resulting from hemisphere damage by ischemia after tMCAO, as well as in the corresponding parts of the brains of sham-operated animals. Real-time reverse transcription polymerase chain reaction expression analysis of 20 genes confirmed the RNA-Seq results. We identified 469 and 1939 genes that exhibited changes in expression of > 1.5-fold at 4.5 and 24 h after tMCAO, respectively. Interestingly, we found 2741 and 752 DEGs under ischemia-reperfusion and sham-operation conditions at 24 h vs. 4.5 h after tMCAO, respectively. The activation of a large number of genes involved in inflammatory, immune and stress responses, apoptosis, ribosome function, DNA replication and other processes was observed in ischemia-reperfusion conditions. Simultaneously, massive down-regulation of the mRNA levels of genes involved in the functioning of neurotransmitter systems was recorded. A Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed that dozens of signaling pathways were associated with DEGs in ischemia-reperfusion conditions. CONCLUSIONS: The data obtained revealed a global profile of gene expression in the rat brain sub-cortex under tMCAO conditions that can be used to identify potential therapeutic targets in the development of new strategies for the prevention and treatment of ischemic stroke.

Developmental stage-specific expression of genes for sphingomyelin synthase in rat brain.

Sphingomyelin synthase genes (Sgms1 and Sgms2) encode the vital enzymes that participate in the processes of membrane transport, cell proliferation and apoptosis. We previously determined the exon-intron structure of Sgms1 and some features of its expression in human and rodent tissues. The circular RNAs (circRNAs) emerging from exons of the 5'-untranslated region (5'-UTR) of Sgms1 were determined. These circRNAs are represented at a high level in the adult brain. Here, we demonstrate that, in contrast to Sgms1, Sgms2 does not contain the multi-exon 5'-UTR but encodes circRNAs, which are composed of the coding region of the gene and are expressed at a low level. We present a study of the expression of sphingomyelin synthase genes in rat brain at embryonic days 7, 9, 13, 17 and 21 and in adult rat brain. In contrast to Sgms1, Sgms2 is expressed at a significantly low level in adult brain. In embryonic rat brain, the mRNA expression of sphingomyelin synthase genes is varied in a developmental stage-specific manner. The level of Sgms1 mRNAs, differing by 5'-UTR-in the formation of which alternative promoters can participate-changes significantly during the process of embryonic development. The expression of circRNAs of Sgms1 was significantly raised during rat embryonic brain development. We assume that the circRNAs are involved in the regulation of sphingomyelin synthase activity in rat brain in different developmental stages.

Circular RNAs-one of the enigmas of the brain.

Circular RNAs (circRNAs) provide a new and relatively unexplored class of noncoding RNAs that are predominantly found in mammalian cells. In this review, we present the latest data regarding the structural organization, possible mechanisms of synthesis, and functions of circRNAs. These transcripts were isolated as an RNA fraction that was resistant to RNase R treatment, which selectively destroys the linear forms of RNA molecules. circRNAs are encoded by orthologous genes in different organisms and show tissue- and organ-specific expression. Currently, the biogenesis and functional properties of circRNAs remain unclear; transcripts of this class, however, remain highly promising targets of research. Some of them have been ascribed the function of "molecular sponges" that can absorb microRNAs, RNA-binding proteins, and small nuclear RNAs. circRNAs are often formed from the RNA portions of protein-coding genes in the course of alternative splicing. Some features of the circRNAs of mammals were demonstrated using 11 circRNAs of the human sphingomyelin synthase 1 gene (SGMS1), which were discovered by us in the brain. These circRNAs consist mainly of portions of the multi-exon 5' untranslated region (5'UTR) of the SGMS1 gene and include one to five exons. The synthesis of circRNAs may be new, previously unknown, function of the multi-exon 5'UTR of genes. This feature is most clearly manifested in the brain, where the level of circRNAs is significantly higher.

Semax, an analog of ACTH(4-7), regulates expression of immune response genes during ischemic brain injury in rats.

Brain stroke continues to claim the lives of million people every year. To build the effective strategies for stroke treatment it is necessary to understand the neuroprotective mechanisms that are able to prevent the ischemic injury. Consisting of the ACTH(4-7) fragment and the tripeptide Pro-Gly-Pro (PGP), the synthetic peptide Semax effectively protects brain against ischemic stroke. However, the molecular mechanisms underlying its neuroprotection and participation of PGP in them are still needed to be clarified. To reveal biological processes and signaling pathways, which are affected by Semax and PGP, we performed the transcriptome analysis of cerebral cortex of rats with focal cerebral ischemia treated by these peptides. The genome-wide biochip data analysis detected the differentially expressed genes (DEGs) and bioinformatic web-tool Ingenuity iReport found DEGs associations with several biological processes and signaling pathways. The immune response is the process most markedly affected by the peptide: Semax enhances antigen presentation signaling pathway, intensifies the effect of ischemia on the interferon signaling pathways and affects the processes for synthesizing immunoglobulins. Semax significantly increased expression of the gene encoding the immunoglobulin heavy chain, highly affects on cytokine, stress response and ribosomal protein-encoding genes after occlusion. PGP treatment of rats with ischemia attenuates the immune activity and suppresses neurotransmission in the CNS. We suppose that neuroprotective mechanism of Semax is realized via the neuroimmune crosstalk, and the new properties of PGP were found under ischemia. Our results provided the basis for further proteomic investigations in the field of searching Semax neuroprotection mechanism.

A New Baltic Population-Specific Human Genetic Marker in the PMCA4 Gene.

OBJECTIVES: The PMCA gene family consists of 4 genes and at least 21 splice variants; among these, the Ca2+ ATPase 4 (PMCA4) gene encodes a plasma membrane protein abundantly expressed in several tissues, including the kidney, heart, and sperm. Knockout of PMCA4 causes infertility due to immotile sperm in mouse models. We therefore investigated variants in this gene for potential association with infertility in groups of Estonian (n = 191) and Latvian (n = 92) men with reduced sperm motility. METHODS: All exons, exon-intron boundaries, 5' and 3' untranslated regions, and the promoter region of the PMCA4 gene were analysed by direct sequencing for a group of Estonian infertile men. Genotyping of guanine and adenine alleles of rs147729934 was performed, using a custom-designed TaqMan® probe for a group of Latvian infertile men as well as additional groups from Latvia and several groups of people with proven ethnicity from the Baltic region. RESULTS: Although we did not identify any significant associations between variants in the gene and infertility, our results indicated that in all studied Latvian and Estonian groups the adenine allele of the variant rs147729934 was present at a higher frequency than expected. Analysis of additional samples indicated that the adenine allele of rs147729934 likely originated once in the modern-day Baltic or western Russia area, as the frequency of the minor adenine allele observed in this region is remarkably higher than that in the general European population. CONCLUSIONS: Our results revealed no significant difference in frequencies of genetic variants in PMCA4 gene between men with normal and those with reduced sperm motility. The adenine allele of the variant rs147729934 is potentially an informative tool for future population studies concerning ancient Baltic and Finno-Ugric history.

Second mutation in PARK2 is absent in patients with sporadic Parkinson's disease and heterozygous exonic deletions/duplications in parkin gene.

AIM OF THE STUDY: Mutations in PARK2 are one of the causes of Parkinson's disease (PD). Deletions and duplications/triplications of one exon or exon groups account for a large proportion of mutations in the gene. At the present time, it is still not fully clear whether heterozygous mutations cause the development of PD. Our study aimed at conducting screening for mutations in PARK2 in patients with a sporadic form of PD to clarify the role of PARK2 in the development of PD. MATERIALS AND METHODS: The cohort of 327 patients with PD was screened by quantitative real-time polimerase chain reaction (PCR) with subsequent Sanger sequencing. RESULTS: It was found that a sufficiently large proportion of these patients (21 patients, 6.4%) were carriers of heterozygous deletions or duplications in PARK2. Analysis of PARK2 exon rearrangement carriers for the presence of point mutations in PARK2 did not reveal any variants with pathogenic significance. CONCLUSIONS: Thus, our data indicate that heterozygous deletions and duplications can play an important role in the pathogenesis of PD and can be considered as dominant mutations with low penetrance.

Comparative analysis of sphingomyelin synthase 1 gene expression at the transcriptional and translational levels in human tissues.

Sphingomyelin synthase 1 (SMS1) catalyses the biosynthesis of sphingomyelin in eukaryotic cells. We have previously determined the structure of the SGMS1 gene encoding this enzyme and a number of its alternative transcripts. Here, we describe a study of the expression of the full-length SMS1 protein and the sum of the alternative transcripts encoding this protein in human tissues. The SMS1 protein and mRNA levels in tissues differed significantly and were not correlated, implying the active post-transcriptional regulation of SMS1 protein expression. The putative truncated isoforms of the SMS1 protein, which are encoded by a number of alternative transcripts, were not detected by immunoblotting and thus may be absent or present in only small amounts.

Circular RNA of the human sphingomyelin synthase 1 gene: Multiple splice variants, evolutionary conservatism and expression in different tissues.

The human sphingomyelin synthase 1 gene (SGMS1) encodes an essential enzyme that is involved in the synthesis of sphingomyelin and diacylglycerol from phosphatidylcholine and ceramide. Among the products of SGMS1, we found new transcripts, circular RNAs (circRNAs), that contain sequences of the gene's 5' untranslated region (5'UTR). Some of them include the gene's coding region and fragments of introns. An analysis of the abundance of circRNAs in human tissues showed that the largest transcripts were predominantly found in different parts of the brain. circRNAs of rat and mouse sphingomyelin synthase 1 orthologous genes were detected and are highly similar to the human SGMS1 gene transcripts. A quantitative analysis of the abundance of such transcripts also revealed their elevated amount in the brain. A computational analysis of sequences of human circRNAs showed their high potential of binding microRNAs (miRNAs), including the miRNAs that form complexes with Ago proteins and the mRNA of SGMS1. We assume that the circRNAs identified here participate in the regulation of the function of the SGMS1 gene in the brain.