Detection and characterization of constitutive replication origins defined by DNA polymerase epsilon.

BACKGROUND: Despite the process of DNA replication being mechanistically highly conserved, the location of origins of replication (ORI) may vary from one tissue to the next, or between rounds of replication in eukaryotes, suggesting flexibility in the choice of locations to initiate replication. Lists of human ORI therefore vary widely in number and location, and there are currently no methods available to compare them. Here, we propose a method of detection of ORI based on somatic mutation patterns generated by the mutator phenotype of damaged DNA polymerase epsilon (POLE). RESULTS: We report the genome-wide localization of constitutive ORI in POLE-mutated human tumors using whole genome sequencing data. Mutations accumulated after many rounds of replication of unsynchronized dividing cell populations in tumors allow to identify constitutive origins, which we show are shared with high fidelity between individuals and tumor types. Using a Smith-Waterman-like dynamic programming approach, we compared replication origin positions obtained from multiple different methods. The comparison allowed us to define a consensus set of replication origins, identified consistently by multiple ORI detection methods. Many DNA features co-localized with the consensus set of ORI, including chromatin loop anchors, G-quadruplexes, S/MARs, and CpGs. Among all features, the H2A.Z histone exhibited the most significant association. CONCLUSIONS: Our results show that mutation-based detection of replication origins is a viable approach to determining their location and associated sequence features.

Multiomics-Based Feature Extraction and Selection for the Prediction of Lung Cancer Survival.

Lung cancer is a global health challenge, hindered by delayed diagnosis and the disease's complex molecular landscape. Accurate patient survival prediction is critical, motivating the exploration of various -omics datasets using machine learning methods. Leveraging multi-omics data, this study seeks to enhance the accuracy of survival prediction by proposing new feature extraction techniques combined with unbiased feature selection. Two lung adenocarcinoma multi-omics datasets, originating from the TCGA and CPTAC-3 projects, were employed for this purpose, emphasizing gene expression, methylation, and mutations as the most relevant data sources that provide features for the survival prediction models. Additionally, gene set aggregation was shown to be the most effective feature extraction method for mutation and copy number variation data. Using the TCGA dataset, we identified 32 molecular features that allowed the construction of a 2-year survival prediction model with an AUC of 0.839. The selected features were additionally tested on an independent CPTAC-3 dataset, achieving an AUC of 0.815 in nested cross-validation, which confirmed the robustness of the identified features.

Multiomics to investigate the mechanisms contributing to repression of PTPRC and SOCS2 in pediatric T-ALL: Focus on miR-363-3p and promoter methylation.

T-cell acute lymphoblastic leukemia (T-ALL) is a heterogeneous and aggressive malignancy arising from T-cell precursors. MiRNAs are implicated in negative regulation of gene expression and when aberrantly expressed contribute to various cancer types, including T-ALL. Previously we demonstrated the oncogenic potential of miR-363-3p overexpression in a subgroup of T-ALL patients. Here, using combined proteomic and transcriptomic approaches, we show that miR-363-3p enhances cell growth of T-ALL in vitro via inhibition of PTPRC and SOCS2, which are implicated in repression of the JAK-STAT pathway. We propose that overexpression of miR-363-3p is a novel mechanism potentially contributing to overactivation of JAK-STAT pathway. Additionally, by combining the transcriptomic and methylation data of T-ALL patients, we show that promoter methylation may also contribute to downregulation of SOCS2 expression and thus potentially to JAK-STAT activation. In conclusion, we highlight aberrant miRNA expression and aberrant promoter methylation as mechanisms, alternative to mutations of JAK-STAT-related genes, which might lead to the upregulation of JAK-dependent signaling in T-ALL.

RNA-seq library preparation for comprehensive transcriptome analysis in cancer cells: The impact of insert size.

With long reads and high coverage, RNA-seq enables comprehensive transcriptome analysis of cancer cells, provided that optimal length of libraries (and their inserts) is assured, to avoid overlap of paired reads and consequent loss of sequencing data. We assessed TruSeq Stranded library preparation protocols (poly(A) enrichment-PA and rRNA depletion-RD) for the thoroughness of transcriptome analysis of a heterogeneous cancer, acute lymphoblastic leukemia. We applied 2x150PE sequencing, >150 M reads/sample on Illumina NovaSeq6000. We show that PA outperforms RD for the analysis of gene expression and structural aberrations. RD is more suitable for detection of various classes of RNAs, mutations or polymorphisms. We demonstrate that reduced RNA fragmentation time (generating longer inserts) positively affects detection of structural RNA changes, without introducing bias into gene expression analysis. We recommend this modification for all RNA-seq studies utilizing reads longer than 75 nt, aimed to go beyond gene expression analysis and to detect also structural changes.

Small RNA-Seq Reveals Similar miRNA Transcriptome in Children and Young Adults with T-ALL and Indicates miR-143-3p as Novel Candidate Tumor Suppressor in This Leukemia.

We aimed to identify miRNAs and pathways specifically deregulated in adolescent and young adult (AYA) T-ALL patients. Small RNA-seq showed no major differences between AYA and pediatric T-ALL, but it revealed downregulation of miR-143-3p in T-ALL patients. Prediction algorithms identified several known and putative oncogenes targeted by this miRNA, including KRAS, FGF1, and FGF9. Pathway analysis indicated signaling pathways related to cell growth and proliferation, including FGFR signaling and PI3K-AKT signaling, with the majority of genes overrepresented in these pathways being predicted targets of hsa-miR-143-3p. By luciferase reporter assays, we validated direct interactions of this miRNA with KRAS, FGF1 and FGF9. In cell proliferation assays, we showed reduction of cell growth upon miR-143-3p overexpression in two T-ALL cell lines. Our study is the first description of the miRNA transcriptome in AYA T-ALL patients and the first report on tumor suppressor potential of miR-143-3p in T-ALL. Downregulation of this miRNA in T-ALL patients might contribute to enhanced growth and viability of leukemic cells. We also discuss the potential role of miR-143-3p in FGFR signaling. Although this requires more extensive validation, it might be an interesting direction, since FGFR inhibition proved promising in preclinical studies in various cancers.

Discrimination between human populations using a small number of differentially methylated CpG sites: a preliminary study using lymphoblastoid cell lines and peripheral blood samples of European and Chinese origin.

BACKGROUND: Epigenetics is one of the factors shaping natural variability observed among human populations. A small proportion of heritable inter-population differences are observed in the context of both the genome-wide methylation level and the methylation status of individual CpG sites. It has been demonstrated that a limited number of carefully selected differentially methylated sites may allow discrimination between main human populations. However, most of the few published results have been performed exclusively on B-lymphocyte cell lines. RESULTS: The goal of our study was to identify a set of CpG sites sufficient to discriminate between populations of European and Chinese ancestry based on the difference in the DNA methylation profile not only in cell lines but also in primary cell samples. The preliminary selection of CpG sites differentially methylated in these two populations (pop-CpGs) was based on the analysis of two groups of commercially available ethnically-specific B-lymphocyte cell lines, performed using Illumina Infinium Human Methylation 450 BeadChip Array. A subset of 10 pop-CpGs characterized by the best differentiating criteria (|Mdiff| > 1, q < 0.05; lack of the confounding genomic features), and 10 additional CpGs in their immediate vicinity, were further tested using pyrosequencing technology in both B-lymphocyte cell lines and in the primary samples of the peripheral blood representing two analyzed populations. To assess the population-discriminating potential of the selected set of CpGs (further referred to as "composite pop (CEU-CHB)-CpG marker"), three classification methods were applied. The predictive ability of the composite 8-site pop (CEU-CHB)-CpG marker was assessed using 10-fold cross-validation method on two independent sets of samples. CONCLUSIONS: Our results showed that less than 10 pop-CpG sites may distinguish populations of European and Chinese ancestry; importantly, this small composite pop-CpG marker performs well in both lymphoblastoid cell lines and in non-homogenous blood samples regardless of a gender.

Effect of the unfolded protein response and oxidative stress on mutagenesis in CSF3R: a model for evolution of severe congenital neutropenia to myelodysplastic syndrome/acute myeloid leukemia.

Severe congenital neutropenia (SCN) is a rare blood disorder characterised by abnormally low levels of circulating neutrophils. The most common recurrent mutations that cause SCN involve neutrophil elastase (ELANE). The treatment of choice for SCN is the administration of granulocyte-colony stimulating factor (G-CSF), which increases the neutrophil number and improves the survival and quality of life. Long-term survival is however linked to the development of myelodysplastic syndrome/acute myeloid leukemia (MDS/AML). About 70% of MDS/AML patients acquire nonsense mutations affecting the cytoplasmic domain of CSF3R (the G-CSF receptor). About 70% of SCN patients with AML harbour additional mutations in RUNX1. We hypothesised that this coding region of CSF3R constitutes a hotspot vulnerable to mutations resulting from excessive oxidative stress or endoplasmic reticulum (ER) stress. We used the murine Ba/F3 cell line to measure the effect of induced oxidative or ER stress on the mutation rate in our hypothesised hotspot of the exogenous human CSF3R, the corresponding region in the endogenous Csf3r, and Runx1. Ba/F3 cells transduced with the cDNA for partial C-terminal of CSF3R fused in-frame with a green fluorescent protein (GFP) tag were subjected to stress-inducing treatment for 30 days (~51 doubling times). The amplicon-based targeted deep sequencing data for days 15 and 30 samples show that although there was increased mutagenesis observed in all the three genes of interest (partial CSF3R, Csf3r and Runx1), there were more mutations in the GFP region compared with the partial CSF3R region. Our findings also indicate that there is no correlation between the stress-inducing chemical treatments and mutagenesis in Ba/F3 cells. Our data suggest that oxidative or ER stress induction does not promote genomic instability, affecting partial C-terminal of the transduced CSF3R, the endogenous Csf3R and the endogenous Runx1 in Ba/F3 cells that could account for these targets to being mutational hotspots. We conclude that other mechanisms to acquire mutations of CSF3R that help drive the evolution of SCN to MDS/AML.

A mathematical model as a tool to identify microRNAs with highest impact on transcriptome changes.

BACKGROUND: Rapid changes in the expression of many messenger RNA (mRNA) species follow exposure of cells to ionizing radiation. One of the hypothetical mechanisms of this response may include microRNA (miRNA) regulation, since the amounts of miRNAs in cells also vary upon irradiation. To address this possibility, we designed experiments using cancer-derived cell lines transfected with luciferase reporter gene containing sequences targeted by different miRNA species in its 3'- untranslated region. We focus on the early time-course response (1 h past irradiation) to eliminate secondary mRNA expression waves. RESULTS: Experiments revealed that the irradiation-induced changes in the mRNA expression depend on the miRNAs which interact with mRNA. To identify the strongest interactions, we propose a mathematical model which predicts the mRNA fold expression changes, caused by perturbation of microRNA-mRNA interactions. Model was applied to experimental data including various cell lines, irradiation doses and observation times, both ours and literature-based. Comparison of modelled and experimental mRNA expression levels given miRNA level changes allows estimating how many and which miRNAs play a significant role in transcriptome response to stress conditions in different cell types. As an example, in the human melanoma cell line the comparison suggests that, globally, a major part of the irradiation-induced changes of mRNA expression can be explained by perturbed miRNA-mRNA interactions. A subset of about 30 out of a few hundred miRNAs expressed in these cells appears to account for the changes. These miRNAs play crucial roles in regulatory mechanisms observed after irradiation. In addition, these miRNAs have a higher average content of GC and a higher number of targeted transcripts, and many have been reported to play a role in the development of cancer. CONCLUSIONS: Our proposed mathematical modeling approach may be used to identify miRNAs which participate in responses of cells to ionizing radiation, and other stress factors such as extremes of temperature, exposure to toxins, and drugs.

Identification of Endogenous Control miRNAs for RT-qPCR in T-Cell Acute Lymphoblastic Leukemia.

Optimal endogenous controls enable reliable normalization of microRNA (miRNA) expression in reverse-transcription quantitative PCR (RT-qPCR). This is particularly important when miRNAs are considered as candidate diagnostic or prognostic biomarkers. Universal endogenous controls are lacking, thus candidate normalizers must be evaluated individually for each experiment. Here we present a strategy that we applied to the identification of optimal control miRNAs for RT-qPCR profiling of miRNA expression in T-cell acute lymphoblastic leukemia (T-ALL) and in normal cells of T-lineage. First, using NormFinder for an iterative analysis of miRNA stability in our miRNA-seq data, we established the number of control miRNAs to be used in RT-qPCR. Then, we identified optimal control miRNAs by a comprehensive analysis of miRNA stability in miRNA-seq data and in RT-qPCR by analysis of RT-qPCR amplification efficiency and expression across a variety of T-lineage samples and T-ALL cell line culture conditions. We then showed the utility of the combination of three miRNAs as endogenous normalizers (hsa-miR-16-5p, hsa-miR-25-3p, and hsa-let-7a-5p). These miRNAs might serve as first-line candidate endogenous controls for RT-qPCR analysis of miRNAs in different types of T-lineage samples: T-ALL patient samples, T-ALL cell lines, normal immature thymocytes, and mature T-lymphocytes. The strategy we present is universal and can be transferred to other RT-qPCR experiments.

Gene signature of the post-Chernobyl papillary thyroid cancer.

PURPOSE: Following the nuclear accidents in Chernobyl and later in Fukushima, the nuclear community has been faced with important issues concerning how to search for and diagnose biological consequences of low-dose internal radiation contamination. Although after the Chernobyl accident an increase in childhood papillary thyroid cancer (PTC) was observed, it is still not clear whether the molecular biology of PTCs associated with low-dose radiation exposure differs from that of sporadic PTC. METHODS: We investigated tissue samples from 65 children/young adults with PTC using DNA microarray (Affymetrix, Human Genome U133 2.0 Plus) with the aim of identifying molecular differences between radiation-induced (exposed to Chernobyl radiation, ECR) and sporadic PTC. All participants were resident in the same region so that confounding factors related to genetics or environment were minimized. RESULTS: There were small but significant differences in the gene expression profiles between ECR and non-ECR PTC (global test, p < 0.01), with 300 differently expressed probe sets (p < 0.001) corresponding to 239 genes. Multifactorial analysis of variance showed that besides radiation exposure history, the BRAF mutation exhibited independent effects on the PTC expression profile; the histological subset and patient age at diagnosis had negligible effects. Ten genes (PPME1, HDAC11, SOCS7, CIC, THRA, ERBB2, PPP1R9A, HDGF, RAD51AP1, and CDK1) from the 19 investigated with quantitative RT-PCR were confirmed as being associated with radiation exposure in an independent, validation set of samples. CONCLUSION: Significant, but subtle, differences in gene expression in the post-Chernobyl PTC are associated with previous low-dose radiation exposure.

Cross talk between cytokine and hyperthermia-induced pathways: identification of different subsets of NF-κB-dependent genes regulated by TNFα and heat shock.

Heat shock inhibits NF-κB signaling, yet the knowledge about its influence on the regulation of NF-κB-dependent genes is limited. Using genomic approaches, i.e., expression microarrays and ChIP-Seq, we aimed to establish a global picture for heat shock-mediated impact on the expression of genes regulated by TNFα cytokine. We found that 193 genes changed expression in human U-2 osteosarcoma cells stimulated with cytokine (including 77 genes with the κB motif in the proximal promoters). A large overlap between sets of genes modulated by cytokine or by heat shock was revealed (86 genes were similarly affected by both stimuli). Binding sites for heat shock-induced HSF1 were detected in regulatory regions of 1/3 of these genes. Furthermore, pre-treatment with heat shock affected the expression of 2/3 of cytokine-modulated genes. In the largest subset of co-affected genes, heat shock suppressed the cytokine-mediated activation (antagonistic effect, 83 genes), which genes were associated with the canonical functions of NF-κB signaling. However, subsets of co-activated and co-repressed genes were also revealed. Importantly, pre-treatment with heat shock resulted in the suppression of NF-κB binding in the promoters of the cytokine-upregulated genes, either antagonized or co-activated by both stimuli. In conclusion, we confirmed that heat shock inhibited activation of genes involved in the classical cytokine-mediated functions of NF-κB. On the other hand, genes involved in transcription regulation were over-represented in the subset of genes upregulated by both stimuli. This suggests the replacement of NF-κB-mediated regulation by heat shock-mediated regulation in the latter subset of genes, which may contribute to the robust response of cells to both stress conditions.

Inferring Bladder Cancer Evolution from Mucosal field Effects by Whole-Organ Spatial Mutational, Proteomic, and Metabolomic Mapping.

Multi-platform mutational, proteomic, and metabolomic spatial mapping was used on the whole-organ scale to identify the molecular evolution of bladder cancer from mucosal field effects. We identified complex proteomic and metabolomic dysregulations in microscopically normal areas of bladder mucosa adjacent to dysplasia and carcinoma in situ. The mutational landscape developed in a background of complex defects of protein homeostasis which included dysregulated nucleocytoplasmic transport, splicesome, ribosome biogenesis, and peroxisome. These changes were combined with altered urothelial differentiation which involved lipid metabolism and protein degradations controlled by PPAR. The complex alterations of proteome were accompanied by dysregulation of gluco-lipid energy-related metabolism. The analysis of mutational landscape identified three types of mutations based on their geographic distribution and variant allele frequencies. The most common were low frequency α mutations restricted to individual mucosal samples. The two other groups of mutations were associated with clonal expansion. The first of this group referred to as ß mutations occurred at low frequencies across the mucosa. The second of this group called γ mutations increased in frequency with disease progression. Modeling of the mutations revealed that carcinogenesis may span nearly 30 years and can be divided into dormant and progressive phases. The α mutations developed gradually in the dormant phase. The progressive phase lasted approximately five years and was signified by the advent of ß mutations, but it was driven by γ mutations which developed during the last 2-3 years of disease progression to invasive cancer. Our study indicates that the understanding of complex alterations involving mucosal microenvironment initiating bladder carcinogenesis can be inferred from the multi-platform whole-organ mapping.

Adaptive filtering of microarray gene expression data based on Gaussian mixture decomposition.

BACKGROUND: DNA microarrays are used for discovery of genes expressed differentially between various biological conditions. In microarray experiments the number of analyzed samples is often much lower than the number of genes (probe sets) which leads to many false discoveries. Multiple testing correction methods control the number of false discoveries but decrease the sensitivity of discovering differentially expressed genes. Concerning this problem, filtering methods for improving the power of detection of differentially expressed genes were proposed in earlier papers. These techniques are two-step procedures, where in the first step some pool of non-informative genes is removed and in the second step only the pool of the retained genes is used for searching for differentially expressed genes. RESULTS: A very important parameter to choose is the proportion between the sizes of the pools of removed and retained genes. A new method, which we propose, allow to determine close to optimal threshold values for sample means and sample variances for gene filtering. The method is adaptive and based on the decomposition of the histogram of gene expression means or variances into mixture of Gaussian components. CONCLUSIONS: By performing analyses of several publicly available datasets and simulated datasets we demonstrate that our adaptive method increases sensitivity of finding differentially expressed genes compared to previous methods of filtering microarray data based on using fixed threshold values.

Sources of high variance between probe signals in Affymetrix short oligonucleotide microarrays.

High density oligonucleotide microarrays present a big challenge for statistical data processing methods which aim to separate changes induced by experimental factors from those caused by artifacts and measurement inaccuracies. Despite huge advances in the field of microarray probe design methods, the signal variation between probes that target a single transcript is substantially larger than their between-replicate array variability, suggesting a large influence of various probe-specific effects that introduce bias to the data. In this work we present the influence of probe-related design variations on the expression intensities of individual probes, focusing on five potential sources of high probe signal variance: the GC composition of the probe, the distance between individual probe target sites, G-quadruplex formation in the probe sequence, the occurrence of sequence motifs complementary to the oligo(dT) primer, and the specificity of unrecognized alternative splicing probeset assignment. By focusing on two high quality microarray datasets based on two distinct array designs we show the extent of variance between probes that target a specific transcript providing guidelines for the future design of microarrays and data processing methods.

Correction: Kurpas et al. Genomic Analysis of SARS-CoV-2 Alpha, Beta and Delta Variants of Concern Uncovers Signatures of Neutral and Non-Neutral Evolution. Viruses 2022, 14, 2375.

Missing Funding [...].

Hematopoietic stem and progenitor cells confer cross-protective trained immunity in mouse models.

Recent studies suggest that infection reprograms hematopoietic stem and progenitor cells (HSPCs) to enhance innate immune responses upon secondary infectious challenge, a process called "trained immunity." However, the specificity and cell types responsible for this response remain poorly defined. We established a model of trained immunity in mice in response to Mycobacterium avium infection. scRNA-seq analysis revealed that HSPCs activate interferon gamma-response genes heterogeneously upon primary challenge, while rare cell populations expand. Macrophages derived from trained HSPCs demonstrated enhanced bacterial killing and metabolism, and a single dose of recombinant interferon gamma exposure was sufficient to induce similar training. Mice transplanted with influenza-trained HSPCs displayed enhanced immunity against M. avium challenge and vice versa, demonstrating cross protection against antigenically distinct pathogens. Together, these results indicate that heterogeneous responses to infection by HSPCs can lead to long-term production of bone marrow derived macrophages with enhanced function and confer cross-protection against alternative pathogens.

Genomic Analysis of SARS-CoV-2 Alpha, Beta and Delta Variants of Concern Uncovers Signatures of Neutral and Non-Neutral Evolution.

Due to the emergence of new variants of the SARS-CoV-2 coronavirus, the question of how the viral genomes evolved, leading to the formation of highly infectious strains, becomes particularly important. Three major emergent strains, Alpha, Beta and Delta, characterized by a significant number of missense mutations, provide a natural test field. We accumulated and aligned 4.7 million SARS-CoV-2 genomes from the GISAID database and carried out a comprehensive set of analyses. This collection covers the period until the end of October 2021, i.e., the beginnings of the Omicron variant. First, we explored combinatorial complexity of the genomic variants emerging and their timing, indicating very strong, albeit hidden, selection forces. Our analyses show that the mutations that define variants of concern did not arise gradually but rather co-evolved rapidly, leading to the emergence of the full variant strain. To explore in more detail the evolutionary forces at work, we developed time trajectories of mutations at all 29,903 sites of the SARS-CoV-2 genome, week by week, and stratified them into trends related to (i) point substitutions, (ii) deletions and (iii) non-sequenceable regions. We focused on classifying the genetic forces active at different ranges of the mutational spectrum. We observed the agreement of the lowest-frequency mutation spectrum with the Griffiths-Tavaré theory, under the Infinite Sites Model and neutrality. If we widen the frequency range, we observe the site frequency spectra much more consistently with the Tung-Durrett model assuming clone competition and selection. The coefficients of the fitting model indicate the possibility of selection acting to promote gradual growth slowdown, as observed in the history of the variants of concern. These results add up to a model of genomic evolution, which partly fits into the classical drift barrier ideas. Certain observations, such as mutation "bands" persistent over the epidemic history, suggest contribution of genetic forces different from mutation, drift and selection, including recombination or other genome transformations. In addition, we show that a "toy" mathematical model can qualitatively reproduce how new variants (clones) stem from rare advantageous driver mutations, and then acquire neutral or disadvantageous passenger mutations which gradually reduce their fitness so they can be then outcompeted by new variants due to other driver mutations.

Involvement of miRNAs in cellular responses to radiation.

PURPOSE: Exposure of living cells to ionizing radiation has different consequences, depending on the dose and cell type. Changes in gene expression at the level of transcription and translation, including those regulated by microRNAs (miRNAs), play a role in the intrinsic radiosensitivity of different cells and define their fate, survival or death. The aim of our work was to examine how ionizing radiation may influence the expression of genes regulated by different miRNAs and miRNA biogenesis. MATERIALS AND METHODS: The work was performed on cultured human melanoma Me45 cells, transiently transfected with plasmids containing Renilla luciferase reporter gene targeted by miRNAs Let-7, miR-21 or miR-24. The levels of reporter mRNAs and mRNAs coding for proteins participating in miRNA biogenesis were assayed at different time points in irradiated and non-irradiated cells using RT-qPCR, and reporter protein by luciferase activity assays. MiRNA-targeted motifs in mRNAs coding for proteins engaged in miRNA biogenesis were extracted from the miRTarBase database. RESULTS: Messenger RNA and protein levels of transfected luciferase genes fluctuated in time in patterns that depended on the type of miRNA regulation and changed upon irradiation of the cells. The average levels of reporter mRNAs were higher in irradiated cells, whereas the levels of proteins changed in either direction. Radiation also influenced the levels of miRNAs and the expression of genes engaged in their biogenesis suggesting that the changes in gene expression following ionizing radiation result mainly from these changes in expression of genes regulating miRNA biogenesis and the influence of miRNA on mRNA translation. CONCLUSIONS: Currently, the responses of cells to ionizing radiation are mainly ascribed to changes in their redox conditions and increased intracellular levels of ROS, but the experiments described here suggest that a further important factor is modulation of translation through changes in biogenesis and levels of miRNAs.

CRISPRi for specific inhibition of miRNA clusters and miRNAs with high sequence homology.

miRNAs form a class of noncoding RNAs, involved in post-transcriptional regulation of gene expression, broadly studied for their involvement in physiological and pathological context. Inhibition of mature miRNA transcripts, commonly used in miRNA loss-of-function experiments, may not be specific in case of miRNAs with high sequence homology, e.g. miRNAs from the same seed family. Phenotypic effects of miRNA repression might be biased by the repression of highly similar miRNAs. Another challenge is simultaneous inhibition of multiple miRNAs encoded within policistronic clusters, potentially co-regulating common biological processes. To elucidate roles of miRNA clusters and miRNAs with high sequence homology, it is of key importance to selectively repress only the miRNAs of interest. Targeting miRNAs on genomic level with CRISPR/dCas9-based methods is an attractive alternative to blocking mature miRNAs. Yet, so far no clear guidelines on the design of CRISPR inhibition (CRISPRi) experiments, specifically for miRNA repression, have been proposed. To address this need, here we propose a strategy for effective inhibition of miRNAs and miRNA clusters using CRISPRi. We provide clues on how to approach the challenges in using CRISPR/dCas in miRNA studies, which include prediction of miRNA transcription start sites (TSSs) and the design of single guide RNAs (sgRNAs). The strategy implements three TSS prediction online tools, dedicated specifically for miRNAs: miRStart, FANTOM 5 miRNA atlas, DIANA-miRGen, and CRISPOR tool for sgRNAs design; it includes testing and selection of optimal sgRNAs. We demonstrate that compared to siRNA/shRNA-based miRNA silencing, CRISPRi improves the repression specificity for miRNAs with highly similar sequence and contribute to higher uniformity of the effects of silencing the whole miRNA clusters. This strategy may be adapted for CRISPR-mediated activation (CRISPRa) of miRNA expression.