The Search for and Functional Analysis of Genetic Variants in microRNA-Binding Sites using Massively Parallel Reporter Assay.

A large-scale search for the genetic variants with a bias in the representation of alleles in transcriptome data (AE SNPs) and the binding sites in microRNA 3'-UTRs was performed and their functional significance was assessed using massively parallel reporter assay (MPRA). Of the 629,559 associated "SNP-gene" pairs (eQTLs) discovered in the human liver tissue according to the GTEx Analysis V8 data, 4394 polymorphic positions in the 3'-UTRs of the genes, which represent the eQTLs for these genes were selected. The TargetScanHuman 7.0 algorithm and PolymiRTS database were searched for the potential microRNA-binding sites. Of the predicted microRNA sites affected by eQTL-SNPs, we selected 51 sites with the best evidence of functionality according to Ago2-CLIP-seq, CLEAR-CLIP, and eCLIP-seq for RNA-binding proteins. For MPRA, a library of the plasmids carrying the main and alternative alleles for each AE SNP (in total, 102 constructs) was created. Allele-specific expression for 6 SNPs was detected by transfection of the HepG2 cell line with the constructed plasmid library and sequencing of target DNA and RNA sequences using the Illumina (MiSeq) platform.

RatDEGdb: a knowledge base of differentially expressed genes in the rat as a model object in biomedical research.

The animal models used in biomedical research cover virtually every human disease. RatDEGdb, a knowledge base of the differentially expressed genes (DEGs) of the rat as a model object in biomedical research is a collection of published data on gene expression in rat strains simulating arterial hypertension, age-related diseases, psychopathological conditions and other human afflictions. The current release contains information on 25,101 DEGs representing 14,320 unique rat genes that change transcription levels in 21 tissues of 10 genetic rat strains used as models of 11 human diseases based on 45 original scientific papers. RatDEGdb is novel in that, unlike any other biomedical database, it offers the manually curated annotations of DEGs in model rats with the use of independent clinical data on equal changes in the expression of homologous genes revealed in people with pathologies. The rat DEGs put in RatDEGdb were annotated with equal changes in the expression of their human homologs in affected people. In its current release, RatDEGdb contains 94,873 such annotations for 321 human genes in 836 diseases based on 959 original scientific papers found in the current PubMed. RatDEGdb may be interesting first of all to human geneticists, molecular biologists, clinical physicians, genetic advisors as well as experts in biopharmaceutics, bioinformatics and personalized genomics. RatDEGdb is publicly available at https://www.sysbio.ru/RatDEGdb.

Social defeat stress in adult mice causes alterations in gene expression, alternative splicing, and the epigenetic landscape of H3K4me3 in the prefrontal cortex: An impact of early-life stress.

Chronic stress is the leading risk factor of a broad range of severe psychopathologies. Nonetheless, the molecular mechanisms triggering these pathological processes are not well understood. In our study, we investigated the effects of 15-day social defeat stress (SDS) on the genome-wide landscape of trimethylation at the 4th lysine residue of histone H3 (H3K4me3) and on the transcriptome in the prefrontal cortex of mice that were reared normally (group SDS) or subjected to maternal separation early in life (group MS+SDS). The mice with the history of stress early in life showed increased susceptibility to SDS in adulthood and demonstrated long-lasting genome-wide alterations in gene expression and splicing as well as in the H3K4me3 epigenetic landscape in the prefrontal cortex. Thus, the high-throughput techniques applied here allowed us to simultaneously detect, for the first time, genome-wide epigenetic and transcriptional changes in the murine prefrontal cortex that are associated with both chronic SDS and increased susceptibility to this stressor.

Data of correlation analysis between the density of H3K4me3 in promoters of genes and gene expression: Data from RNA-seq and ChIP-seq analyses of the murine prefrontal cortex.

H3K4me3 is typically found in the promoter region of genes and is a mark associated with an open chromatin state and active gene transcription. Nonetheless, the role of H3K4me3 in the regulation of transcription is still debated. To improve the understanding of the connection between H3K4me3 density in promoters and gene expression, we assessed the correlation between these two parameters. We utilized genome-wide high-throughput RNA sequencing (RNA-seq) data and H3K4me3-based chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq), carried out on the same samples of the prefrontal cortex from 10 male C57Bl6 mice with different stress experience [Social defeat stress in adult mice causes alterations in gene expression, alternative splicing, and the epigenetic landscape of H3K4me3 in the prefrontal cortex: an impact of early-life stress, 1]. In addition, we assessed the correlation between H3K4me3 density and gene expression in datasets of cell-specific genes. Altogether, the results are useful for the elucidation of H3K4me3 involvement in the regulation of transcription in the murine prefrontal cortex.

Evaluation of various RNA-seq approaches for identification of gene outrons in the flatworm Opisthorchis felineus.

The parasitic flatworm Opisthorchis felineus is one of the causative agents of opisthorchiasis in humans. Recently, we assembled the O. felineus genome, but the correct genome annotation by means of standard methods was hampered by the presence of spliced leader trans-splicing (SLTS). As a result of SLTS, the original 5'-end (outron) of the transcripts is replaced by a short spliced leader sequence donated from a specialized SL RNA. SLTS is involved in the RNA processing of more than half of O. felineus genes, making it hard to determine the structure of outrons and bona fide transcription start sites of the corresponding genes and operons, being based solely on mRNA-seq data. In the current study, we tested various experimental approaches for identifying the sequences of outrons in O. felineus using massive parallel sequencing. Two of them were developed by us for targeted sequencing of already processed branched outrons. One was based on sequence-specific reverse transcription from the SL intron toward the 5'-end of the Y-branched outron. The other used outron hybridization with an immobilized single-stranded DNA probe complementary to the SL intron. Additionally, two approaches to the sequencing of rRNA-depleted total RNA were used, allowing the identification of a wider range of transcripts compared to mRNAseq. One is based on the enzymatic elimination of overrepresented cDNAs, the other utilizes exonucleolytic degradation of uncapped RNA by Terminator enzyme. By using the outron-targeting methods, we were not able to obtain the enrichment of RNA preparations by processed outrons, which is most likely indicative of a rapid turnover of these trans-splicing intermediate products. Of the two rRNA depletion methods, a method based on the enzymatic normalization of cDNA (Zymo-Seq RiboFree) showed high efficiency. Compared to mRNA-seq, it provides an approximately twofold increase in the fraction of reads originating from outrons and introns. The results suggest that unprocessed nascent transcripts are the main source of outron sequences in the RNA pool of O. felineus.

Genomics and proteomics of the liver fluke Opisthorchis felineus.

The causative agent of opisthorchiasis, the liver fluke Opisthorchis felineus (Rivolta, 1884) is one of the helminths of humans and animals in Russia. Together with closely related species of trematodes O. viverrini (Poirier, 1886) and Clonorchis sinensis (Loos, 1907), O. felineus is a part of a triad of epidemiologically important trematodes in the family Opisthorchiidae. Adult O. felineus worms infest the hepatobiliary system of warm-blooded animals and might provoke the development of severe pathologies, including malignancy of bile duct epithelium. The high medical importance of O. felineus attracts the attention of researchers. This review briefly summarizes the data about O. felineus genomics and proteomics. The review provides a comparative analysis of the number of genes and sizes of nuclear genomes of a number of flatworms, the distribution of intron lengths, as well as results of synteny between the O. felineus, O. viverrini and C. sinensis genomes. Special attention is paid to a particular form of RNA processing known as trans-splicing, widely presented in the opisthorchiid genomes. We also provide the results of a comparative analysis of the xenobiotic metabolizing system between parasitic and free-living flatworms. Moreover, data on parasitic granulins, which are potential promoters of cholangiocyte neoplasia, are also presented. Data on the O. felineus genomics and proteomics provide first insights into the structural and functional organization of the genome of this parasitic flatworm with a complex life cycle as well as provide a significant contribution to our understanding of "host-parasite" interaction and evolution of this group of parasitic flatworms.

Strain-Specific Single-Nucleotide Polymorphisms in Hypertensive ISIAH Rats.

Single-nucleotide polymorphisms (SNPs) in the coding and regulatory regions of genes can affect transcription rate and translation efficiency, modify protein function, and, in some cases, cause the development of diseases. In the current study, the RNA-Seq approach has been used to discover strain-specific SNPs in ISIAH (inherited stress-induced arterial hypertension) rats, which are known as a model of stress-induced arterial hypertension. The comparison of the ISIAH SNPs with genome sequencing data available for another 42 rat strains and substrains, 11 of them known as hypertensive, showed a considerable genetic distance between the genotypes of ISIAH and all other rat strains and substrains. The study revealed 1849 novel SNPs specific for ISIAH rats and 158 SNPs present only in the genotypes of hypertensive rats. Amino acid substitutions with possible deleterious effect on protein function were detected. Several of them were found in the genes associated with hypertension. These SNPs may be considered as novel molecular targets for further studies aimed at assessing their potential in the therapy of stress-induced hypertension.

Analysis of the placental tissue transcriptome of normal and preeclampsia complicated pregnancies.

Preeclampsia is one of the most severe gestational complications which is one of the leading causes of maternal and perinatal morbidity and mortality. A growth in the incidence of severe and combined forms of the pathology has been observed in recent years. According to modern concepts, inadequate cytotrophoblast invasion into the spiral arteries of the uterus and development of the ischemia-reperfusion syndrome in the placental tissue play the leading role in the development of preeclampsia, which is characterized by multipleorgan failure. In this regard, our work was aimed at studying the patterns of placental tissue transcriptome that are specific to females with PE and with physiological pregnancy, as well as identifying the potential promising biomarkers and molecular mechanisms of this pathology. We have identified 63 genes whose expression proved to differ significantly in the placental tissue of females with PE and with physiological pregnancy. A cluster of differentially expressed genes (DEG) whose expression level is increased in patients with preeclampsia includes not only the known candidate genes that have been identified in many other genome-wide studies (e.g., LEP, BHLHB2, SIGLEC6, RDH13, BCL6), but also new genes (ANKRD37, SYDE1, CYBA, ITGB2, etc.), which can be considered as new biological markers of preeclampsia and are of further interest. The results of a functional annotation of DEG show that the development of preeclampsia may be related to a stress response, immune processes, the regulation of cell-cell interactions, intracellular signaling cascades, etc. In addition, the features of the differential gene expression depending on preeclampsia severity were revealed. We have found evidence of the important role of the molecular mechanisms responsible for the failure of immunological tolerance and initiation of the pro-inflammatory cascade in the development of severe preeclampsia. The results obtained elaborate the concept of the pathophysiology of preeclampsia and contain the information necessary to work out measures for targeted therapy of this disease. ;

Detection of target genes of FOXA transcription factors involved in proliferation control.

To reveal the mechanism of tumor-suppressing activity of FOXA proteins in liver, a search for potential target genes of these transcription factors involved in proliferation control was carried out. In the first step, we have used data from the literature concerning gene expression in mouse liver (high content of FOXA proteins) and kidney (FOXA expression is absent) obtained by hybridization on microchips. A search for FOXA binding sites in regulatory regions of forty differentially expressing genes involved in proliferation control was carried out using the computer method SITECON. Eleven genes containing clusters of potential FOXA sites incorporating 3-6-fold repeats of TTTG were revealed. The FOXA-specific interaction with such microsatellite sites was confirmed by gel-retardation technique using the GST-fused protein containing the DNA-binding domain of FOXA2. Six genes containing clusters of confirmed binding sites--Cul2, Cdc73, Ptk, Pdcd, Creb, and Ppp2r5d--were selected. The effect of hepatocarcinogen orthoaminoazotoluene (OAT), which lowers the FOXA activity, on expression of these genes was studied by the real-time PCR. OAT was shown to increase sharply the level of mRNA of the Cul2 and Cdc73 genes.

FOXA transcription factors determine the amplitude of glucocorticoid induction of tyrosine aminotransferase in mice.

o-Aminoazotoluene was more potent than 3'-methyl-4-dimethylaminoazobenzene in modulating glucocorticoid induction of tyrosine aminotransferase and DNA-binding activity of FOXA (HNF3) in 12-day-old ICR mice. In adult animals, induction of tyrosine aminotransferase and FOXA activity were modulated by o-aminoazotoluene, while 3'-methyl-4-dimethylaminoazobenzene was ineffective. Our results suggest that FOXA proteins determine glucocorticoid induction of tyrosine aminotransferase in mice (similarly to rats).

[An effect of the microsporidian Vairimorpha ephestiae (Microsporidia: Burenellidae) on activity and spectrum of nonspecific esterases in different tissues of the greater wax moth galleria mellonella (Lepidoptera: Pyralidae) larvae].

The effect of the microsporidian Vairimorpha ephestiae Matted (Microsporidia: Burenellidae) on nonspecific esterases was studied in hemolymph, fat body and midgut of the larvae of Galleria mellonella L. (Lepidoptera: Pyralidae). Esterase patterns were analyzed by the polyacrylamide gel electrophoresis, the total esterase activity was detected spectrophotometrically. The increase of total esterase activity was registered in hemolymph of inflected larvae. An overexpression of esterase isozyme in hemolymph was already detected at the 3rd day post infection. No changes in esterases pattern were observed in the fat body's homogenates of the G. mellonella larvae possessing the symptoms of microsporidiosis. The degradation of esterase isozymes and the decrease of total esterase activity in the pattern of the midgut homogenates of infected larvae were registered during parasite sporogony. The greatest esterase activity in hemolymph and midgut tissues was registered during vegetative reproduction of parasite, but the least level of esterase activity was observed during mass sporogony of microsporidia.