Intron retention in PI-PLC γ1 mRNA as a key mechanism affecting MMP expression in human primary fibroblast-like synovial cells.

The intron retention (IR) is a phenomenon utilized by cells to allow diverse fates at the same mRNA, leading to a different pattern of synthesis of the same protein. In this study, we analyzed the modulation of phosphoinositide-specific phospholipase C (PI-PLC) enzymes by Harpagophytum procumbens extract (HPE) in synoviocytes from joins of osteoarthritis (OA) patients. In some samples, the PI-PLC γ1 isoform mature mRNA showed the IR and, in these synoviocytes, the HPE treatment increased the phenomenon. Moreover, we highlighted that as a consequence of IR, a lower amount of PI-PLC γ1 was produced. The decrease of PI-PLC γ1 was associated with the decrease of metalloprotease-3 (MMP-3), and MMP-13, and ADAMTS-5 after HPE treatment. The altered expression of MMPs is a hallmark of the onset and progression of OA, thus substances able to decrease their expression are very desirable. The interesting outcomes of this study are that 35% of analyzed synovial tissues showed the IR phenomenon in the PI-PLC γ1 mRNA and that the HPE treatment increased this phenomenon. For the first time, we found that the decrease of PI-PLC γ1 protein in synoviocytes interferes with MMP production, thus affecting the pathways involved in the MMP expression. This finding was validated by the silencing of PI-PLC γ1 in synoviocytes where the IR phenomenon was not present. Our results shed new light on the biochemical mechanisms involved in the degrading enzyme production in the joint of OA patients, suggesting a new therapeutic target and highlighting the importance of personalized medicine.

Innovative construction of the first reliable catalogue of bovine circular RNAs.

The aim of this study was to compare the circular transcriptome of divergent tissues in order to understand: i) the presence of circular RNAs (circRNAs) that are not exonic circRNAs, i.e. originated from backsplicing involving known exons and, ii) the origin of artificial circRNA (artif_circRNA), i.e. circRNA not generated in-vivo. CircRNA identification is mostly an in-silico process, and the analysis of data from the BovReg project (https://www.bovreg.eu/) provided an opportunity to explore new ways to identify reliable circRNAs. By considering 117 tissue samples, we characterized 23,926 exonic circRNAs, 337 circRNAs from 273 introns (191 ciRNAs, 146 intron circles), 108 circRNAs from small non-coding genes and nearly 36.6K circRNAs classified as other_circRNAs. Furthermore, for 63 of those samples we analysed in parallel data from total-RNAseq (ribosomal RNAs depleted prior to library preparation) with paired mRNAseq (library prepared with poly(A)-selected RNAs). The high number of circRNAs detected in mRNAseq, and the significant number of novel circRNAs, mainly other_circRNAs, led us to consider all circRNAs detected in mRNAseq as artificial. This study provided evidence of 189 false entries in the list of exonic circRNAs: 103 artif_circRNAs identified by total RNAseq/mRNAseq comparison using two circRNA tools, 26 probable artif_circRNAs, and 65 identified by deep annotation analysis. Extensive benchmarking was performed (including analyses with CIRI2 and CIRCexplorer-2) and confirmed 94% of the 23,737 reliable exonic circRNAs. Moreover, this study demonstrates the effectiveness of a panel of highly expressed exonic circRNAs (5-8%) in analysing the tissue specificity of the bovine circular transcriptome.

The spliceosome impacts morphogenesis in the human fungal pathogen Candida albicans.

At human body temperature, the fungal pathogen Candida albicans can transition from yeast to filamentous morphologies in response to host-relevant cues. Additionally, elevated temperatures encountered during febrile episodes can independently induce C. albicans filamentation. However, the underlying genetic pathways governing this developmental transition in response to elevated temperatures remain largely unexplored. Here, we conducted a functional genomic screen to unravel the genetic mechanisms orchestrating C. albicans filamentation specifically in response to elevated temperature, implicating 45% of genes associated with the spliceosome or pre-mRNA splicing in this process. Employing RNA-Seq to elucidate the relationship between mRNA splicing and filamentation, we identified greater levels of intron retention in filaments compared to yeast, which correlated with reduced expression of the affected genes. Intriguingly, homozygous deletion of a gene encoding a spliceosome component important for filamentation (PRP19) caused even greater levels of intron retention compared with wild type and displayed globally dysregulated gene expression. This suggests that intron retention is a mechanism for fine-tuning gene expression during filamentation, with perturbations of the spliceosome exacerbating this process and blocking filamentation. Overall, this study unveils a novel biological process governing C. albicans filamentation, providing new insights into the complex regulation of this key virulence trait.IMPORTANCEFungal pathogens such as Candida albicans can cause serious infections with high mortality rates in immunocompromised individuals. When C. albicans is grown at temperatures encountered during human febrile episodes, yeast cells undergo a transition to filamentous cells, and this process is key to its virulence. Here, we expanded our understanding of how C. albicans undergoes filamentation in response to elevated temperature and identified many genes involved in mRNA splicing that positively regulate filamentation. Through transcriptome analyses, we found that intron retention is a mechanism for fine-tuning gene expression in filaments, and perturbation of the spliceosome exacerbates intron retention and alters gene expression substantially, causing a block in filamentation. This work adds to the growing body of knowledge on the role of introns in fungi and provides new insights into the cellular processes that regulate a key virulence trait in C. albicans.

Evolution of Whirly1 in the angiosperms: sequence, splicing, and expression in a clade of early transitional mycoheterotrophic orchids.

The plastid-targeted transcription factor Whirly1 (WHY1) has been implicated in chloroplast biogenesis, plastid genome stability, and fungal defense response, which together represent characteristics of interest for the study of autotrophic losses across the angiosperms. While gene loss in the plastid and nuclear genomes has been well studied in mycoheterotrophic plants, the evolution of the molecular mechanisms impacting genome stability is completely unknown. Here, we characterize the evolution of WHY1 in four early transitional mycoheterotrophic orchid species in the genus Corallorhiza by synthesizing the results of phylogenetic, transcriptomic, and comparative genomic analyses with WHY1 genomic sequences sampled from 21 orders of angiosperms. We found an increased number of non-canonical WHY1 isoforms assembled from all but the greenest Corallorhiza species, including intron retention in some isoforms. Within Corallorhiza, phylotranscriptomic analyses revealed the presence of tissue-specific differential expression of WHY1 in only the most photosynthetically capable species and a coincident increase in the number of non-canonical WHY1 isoforms assembled from fully mycoheterotrophic species. Gene- and codon-level tests of WHY1 selective regimes did not infer significant signal of either relaxed selection or episodic diversifying selection in Corallorhiza but did so for relaxed selection in the late-stage full mycoheterotrophic orchids Epipogium aphyllum and Gastrodia elata. Additionally, nucleotide substitutions that most likely impact the function of WHY1, such as nonsense mutations, were only observed in late-stage mycoheterotrophs. We propose that our findings suggest that splicing and expression changes may precede the selective shifts we inferred for late-stage mycoheterotrophic species, which therefore does not support a primary role for WHY1 in the transition to mycoheterotrophy in the Orchidaceae. Taken together, this study provides the most comprehensive view of WHY1 evolution across the angiosperms to date.

SimSpliceEvol2: alternative splicing-aware simulation of biological sequence evolution and transcript phylogenies.

BACKGROUND: SimSpliceEvol is a tool for simulating the evolution of eukaryotic gene sequences that integrates exon-intron structure evolution as well as the evolution of the sets of transcripts produced from genes. It takes a guide gene tree as input and generates a gene sequence with its transcripts for each node of the tree, from the root to the leaves. However, the sets of transcripts simulated at different nodes of the guide gene tree lack evolutionary connections. Consequently, SimSpliceEvol is not suitable for evaluating methods for transcript phylogeny inference or gene phylogeny inference that rely on transcript conservation. RESULTS: Here, we introduce SimSpliceEvol2, which, compared to the first version, incorporates an explicit model of transcript evolution for simulating alternative transcripts along the branches of a guide gene tree, as well as the transcript phylogenies inferred. We offer a comprehensive software with a graphical user interface and an updated version of the web server, ensuring easy and user-friendly access to the tool. CONCLUSION: SimSpliceEvol2 generates synthetic datasets that are useful for evaluating methods and tools for spliced RNA sequence analysis, such as spliced alignment methods, methods for identifying conserved transcripts, and transcript phylogeny reconstruction methods. The web server is accessible at https://simspliceevol.cobius.usherbrooke.ca , where you can also download the standalone software. Comprehensive documentation for the software is available at the same address. For developers interested in the source code, which requires the installation of all prerequisites to run, it is provided at  https://github.com/UdeS-CoBIUS/SimSpliceEvol .

Methylation of the Selected TP53 Introns in Advanced-Stage Ovarian Carcinomas.

Objective: Advanced-stage ovarian cancer (OC) is among the most fatal female genital tract neoplasms worldwide. Although different genetic mechanisms have been shown to be involved in ovarian carcinogenesis, the role of TP53 introns methylation is still unresolved. We performed methylation analysis of introns 1, 3, and 4 of the TP53 to identify patterns in primary stage III OCs, corresponding metastases, and healthy tissues. Methods: The study involved samples of paraffin-embedded tissues obtained from 80 patients with stage III OCs, who underwent surgery at the Department of Gynecology and Gynecologic Oncology of the Military Institute of Medicine in Warsaw, Poland. Altogether, 40 serous-type G2/3 OCs and 40 endometrioid-type G2/3 OCs were included. From the same patient, metastatic and normal tissues were simultaneously analyzed. As a control group, 80 tissue samples were collected from patients after bariatric operations. Human ovarian cancer A2780 cell line was also investigated. Total genomic DNA was isolated from paraffin-embedded tissue blocks and the methylation analysis was performed by bisulfite DNA conversion, DNA amplification with specific primers, cloning, and DNA sequencing. Results: All of the samples of intron 1 of TP53 were un-methylated in OCs, metastatic tissues, and in healthy tissues from the same patient. Also, no methylation of TP53 intron 1 was detected in cells from the human A2780 ovarian cancer cell line and in all samples from control group. In all samples, introns 3 and 4 of the TP53 were methylated in primary tumors, metastatic tissue, and in healthy tissue from the same patient, in human A2780 ovarian cell line, and in DNA samples from healthy patients. None of the clinicopatholocal features was related to the TP53 introns methylation status. Conclusions: Our data on TP53 introns methylation sheds new light on the mechanism of p53 activity for a better understanding of cancer biology. The study suggests the existence of an additional regulation rule of TP53 activity that involves demethylation-methylation mechanisms. Methylation at introns 3 and 4 may also overall help in protecting TP53 against damage by viral restrictases or viral DNA integration.

Aberrant spliceosome activity via elevated intron retention and upregulation and phosphorylation of SF3B1 in chronic lymphocytic leukemia.

Splicing factor 3b subunit 1 (SF3B1) is the largest subunit and core component of the spliceosome. Inhibition of SF3B1 was associated with an increase in broad intron retention (IR) on most transcripts, suggesting that IR can be used as a marker of spliceosome inhibition in chronic lymphocytic leukemia (CLL) cells. Furthermore, we separately analyzed exonic and intronic mapped reads on annotated RNA-sequencing transcripts obtained from B cells (n = 98 CLL patients) and healthy volunteers (n = 9). We measured intron/exon ratio to use that as a surrogate for alternative RNA splicing (ARS) and found that 66% of CLL-B cell transcripts had significant IR elevation compared with normal B cells (NBCs) and that correlated with mRNA downregulation and low expression levels. Transcripts with the highest IR levels belonged to biological pathways associated with gene expression and RNA splicing. A >2-fold increase of active pSF3B1 was observed in CLL-B cells compared with NBCs. Additionally, when the CLL-B cells were treated with macrolides (pladienolide-B), a significant decrease in pSF3B1, but not total SF3B1 protein, was observed. These findings suggest that IR/ARS is increased in CLL, which is associated with SF3B1 phosphorylation and susceptibility to SF3B1 inhibitors. These data provide additional support to the relevance of ARS in carcinogenesis and evidence of pSF3B1 participation in this process.

Comprehensive analysis of insertion sequences within rRNA genes of CPR bacteria and biochemical characterization of a homing endonuclease encoded by these sequences.

The Candidate Phyla Radiation (CPR) represents an extensive bacterial clade comprising primarily uncultured lineages and is distinguished from other bacteria by a significant prevalence of insertion sequences (ISs) within their rRNA genes. However, our understanding of the taxonomic distribution and characteristics of these ISs remains limited. In this study, we used a comprehensive approach to systematically determine the nature of the rRNA ISs in CPR bacteria. The analysis of hundreds of rRNA gene sequences across 65 CPR phyla revealed that ISs are present in 48% of 16S rRNA genes and 82% of 23S rRNA genes, indicating a broad distribution across the CPR clade, with exceptions in the 16S and 23S rRNA genes of Candidatus (Ca.) Saccharibacteria and the 16S rRNA genes of Ca. Peregrinibacteria. Over half the ISs display a group-I-intron-like structure, whereas specific 16S rRNA gene ISs display features reminiscent of group II introns. The ISs frequently encode proteins with homing endonuclease (HE) domains, centered around the LAGLIDADG motif. The LAGLIDADG HE (LHE) proteins encoded by the rRNA ISs of CPR bacteria predominantly have a single-domain structure, deviating from the usual single- or double-domain configuration observed in typical prokaryotic LHEs. Experimental analysis of one LHE protein, I-ShaI from Ca. Shapirobacteria, confirmed that its endonuclease activity targets the DNA sequence of its insertion site, and chemical cross-linking experiments demonstrated its capacity to form homodimers. These results provide robust evidence supporting the hypothesis that the explosive proliferation of rRNA ISs in CPR bacteria was facilitated by mechanisms involving LHEs. IMPORTANCE: Insertion sequences (ISs) in rRNA genes are relatively limited and infrequent in most bacterial phyla. With a comprehensive bioinformatic analysis, we show that in CPR bacteria, these ISs occur in 48% of 16S rRNA genes and 82% of 23S rRNA genes. We also report the systematic and biochemical characterization of the LAGLIDADG homing endonucleases (LHEs) encoded by these ISs in the first such analysis of the CPR bacteria. This study significantly extends our understanding of the phylogenetic positions of rRNA ISs within CPR bacteria and the biochemical features of their LHEs.

Genes divided according to the relative position of the longest intron show increased representation in different KEGG pathways.

Despite the fact that introns mean an energy and time burden for eukaryotic cells, they play an irreplaceable role in the diversification and regulation of protein production. As a common feature of eukaryotic genomes, it has been reported that in protein-coding genes, the longest intron is usually one of the first introns. The goal of our work was to find a possible difference in the biological function of genes that fulfill this common feature compared to genes that do not. Data on the lengths of all introns in genes were extracted from the genomes of six vertebrates (human, mouse, koala, chicken, zebrafish and fugu) and two other model organisms (nematode worm and arabidopsis). We showed that more than 40% of protein-coding genes have the relative position of the longest intron located in the second or third tertile of all introns. Genes divided according to the relative position of the longest intron were found to be significantly increased in different KEGG pathways. Genes with the longest intron in the first tertile predominate in a range of pathways for amino acid and lipid metabolism, various signaling, cell junctions or ABC transporters. Genes with the longest intron in the second or third tertile show increased representation in pathways associated with the formation and function of the spliceosome and ribosomes. In the two groups of genes defined in this way, we further demonstrated the difference in the length of the longest introns and the distribution of their absolute positions. We also pointed out other characteristics, namely the positive correlation between the length of the longest intron and the sum of the lengths of all other introns in the gene and the preservation of the exact same absolute and relative position of the longest intron between orthologous genes.

From bedside to genetic analysis: New insights into pathophysiology of melanoma, basal cell carcinoma, and other cancers.

OBJECTIVE: Patients with myotonic muscular dystrophy (MMD) were observed to have numerous basal cell carcinoma (BCC) and abnormal dysplastic nevi (DN) on non-sun exposed skin. Simultaneously a large study published in the Journal of American Medical Association (JAMA) illustrated that patients with MMD have "overall" an increased risk for cancer development. Based on these findings, this author in 2010 postulated that dysregulation of RNA binding proteins (RBP), responsible for clinical manifestations of MMD, is also responsible for the development of BCC and melanoma. METHODS: To report new research elucidating the etiology of melanoma, BCC, MMD-induced cancers, and potentially other environmentally induced malignancies. RESULTS: Dysregulation of RBP induces aberrant mRNA splicing; recent data indicates that abnormal mRNA splicing not just plays a key role in the pathogenesis of melanoma but is a hallmark of essentially all human malignancies. CONCLUSION: The author's hypothesis is that ultraviolet (UV) radiation induces DNA damage in intronic regions of a variety of genes. Furthermore, these UV-induced abnormal DNA dimers, repeats and mutations interfere with normal mRNA splicing thus producing abnormal proteins. These abnormal proteins in turn activate oncogenic pathways such as hedgehog, MAP kinase, and WNT.

The complete mitochondrial genomes of five critical phytopathogenic Bipolaris species: features, evolution, and phylogeny.

In the present study, three mitogenomes from the Bipolaris genus (Bipolaris maydis, B. zeicola, and B. oryzae) were assembled and compared with the other two reported Bipolaris mitogenomes (B. oryzae and B. sorokiniana). The five mitogenomes were all circular DNA molecules, with lengths ranging from 106,403 bp to 135,790 bp. The mitogenomes of the five Bipolaris species mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). The PCG length, AT skew and GC skew showed large variability among the 13 PCGs in the five mitogenomes. Across the 13 core PCGs tested, nad6 had the least genetic distance among the 16 Pleosporales species we investigated, indicating that this gene was highly conserved. In addition, the Ka/Ks values for all 12 core PCGs (excluding rps3) were < 1, suggesting that these genes were subject to purifying selection. Comparative mitogenomic analyses indicate that introns were the main factor contributing to the size variation of Bipolaris mitogenomes. The introns of the cox1 gene experienced frequent gain/loss events in Pleosporales species. The gene arrangement and collinearity in the mitogenomes of the five Bipolaris species were almost highly conserved within the genus. Phylogenetic analysis based on combined mitochondrial gene datasets showed that the five Bipolaris species formed well-supported topologies. This study is the first report on the mitogenomes of B. maydis and B. zeicola, as well as the first comparison of mitogenomes among Bipolaris species. The findings of this study will further advance investigations into the population genetics, evolution, and genomics of Bipolaris species.

Design and Synthesis of Circular RNA Expression Vectors.

The recent success of the synthetic mRNA-based anti-COVID-19 vaccines has demonstrated the broad potential of the mRNA platform for applications in medicine, thanks to the combined efforts of a small community that has vastly improved key determinants such as design and formulation of synthetic mRNA during the past three decades. However, the cost of production and sensitivity to enzymatic degradation are still limiting the broader application of synthetic mRNA for therapeutic applications. The increased interest in mRNA-based technologies has spurred a renaissance for circular RNA (circRNA), as the lack of free 5' and 3' ends substantially increases resistance against enzymatic degradation in biological systems and does not require expensive cap analogs, as translation is controlled by an Internal Ribosome Entry Site (IRES) sequence. Thus, it can be expected that circRNA will play an important role for future mRNA therapeutics. Here we provide a detailed guide to the production of synthetic circRNA.

Intron lariat spliceosomes convert lariats to true circles: implications for intron transposition.

Rare, full-length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envisioned and tested a hypothesis for their formation using Saccharomyces cerevisiae, documenting full-length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full-length and processed circles. Postsplicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

A multi-tissue, splicing-based joint transcriptome-wide association study identifies susceptibility genes for breast cancer.

Splicing-based transcriptome-wide association studies (splicing-TWASs) of breast cancer have the potential to identify susceptibility genes. However, existing splicing-TWASs test the association of individual excised introns in breast tissue only and thus have limited power to detect susceptibility genes. In this study, we performed a multi-tissue joint splicing-TWAS that integrated splicing-TWAS signals of multiple excised introns in each gene across 11 tissues that are potentially relevant to breast cancer risk. We utilized summary statistics from a meta-analysis that combined genome-wide association study (GWAS) results of 424,650 women of European ancestry. Splicing-level prediction models were trained in GTEx (v.8) data. We identified 240 genes by the multi-tissue joint splicing-TWAS at the Bonferroni-corrected significance level; in the tissue-specific splicing-TWAS that combined TWAS signals of excised introns in genes in breast tissue only, we identified nine additional significant genes. Of these 249 genes, 88 genes in 62 loci have not been reported by previous TWASs, and 17 genes in seven loci are at least 1 Mb away from published GWAS index variants. By comparing the results of our splicing-TWASs with previous gene-expression-based TWASs that used the same summary statistics and expression prediction models trained in the same reference panel, we found that 110 genes in 70 loci that are identified only by the splicing-TWASs. Our results showed that for many genes, expression quantitative trait loci (eQTL) did not show a significant impact on breast cancer risk, whereas splicing quantitative trait loci (sQTL) showed a strong impact through intron excision events.

Comparative analysis of mitochondrial genomes in Ceratocystis fimbriata complex across diverse hosts.

The decline ofAcacia mangiumWilld. in Malaysia, especially in Sabah since 2010, is primarily due to Ceratocystiswilt and canker disease (CWCD) caused by theCeratocystis fimbriataEllis & Halst. complex. This study was aimed to investigate the mitochondrial genome architecture of two differentC. fimbriatacomplex isolates from Malaysia: one fromA. mangiumin Pahang (FRIM1162) and another fromEucalyptus pellitain Sarawak (FRIM1441). This research employed Next-Generation Sequencing (NGS) to contrast genomes from diverse hosts with nine additional mitochondrial sequences, identifying significant genetic diversity and mutational hotspots in the mitochondrial genome alignment. The mitochondrial genome-based phylogenetic analysis revealed a significant genetic relationship between the studied isolates and theC. fimbriatacomplex in the South American Subclade, indicating that theC. fimbriatacomplex discovered in Malaysia isC. manginecans. The comparative mitochondrial genome demonstrates the adaptability of the complex due to mobile genetic components and genomic rearrangements in the studiedfungal isolates. This research enhances our knowledge of the genetic diversity and evolutionary patterns within theC. fimbriatacomplex, aiding in a deeper understanding of fungal disease development and host adaption processes. The acquired insights are crucial for creating specific management strategies for CWCD, improving the overall understanding of fungal disease evolution and control.

Partial loss of desmin expression due to a leaky splice site variant in the human DES gene is associated with neuromuscular transmission defects.

Recessive desminopathies are rare and often present as severe early-onset myopathy. Here we report a milder phenotype in three unrelated patients from southern India (2 M, 1F) aged 16, 21, and 22 years, who presented with childhood-onset, gradually progressive, fatigable limb-girdle weakness, ptosis, speech and swallowing difficulties, without cardiac involvement. Serum creatine kinase was elevated, and repetitive nerve stimulation showed decrement in all. Clinical improvement was noted with pyridostigmine and salbutamol in two patients. All three patients had a homozygous substitution in intron 5: DES(NM_001927.4):c.1023+5G>A, predicted to cause a donor splice site defect. Muscle biopsy with ultrastructural analysis suggested myopathy with myofibrillar disarray, and immunohistochemistry showed partial loss of desmin with some residual staining, while western blot analysis showed reduced desmin. RT-PCR of patient muscle RNA revealed two transcripts: a reduced normal desmin transcript and a larger abnormal transcript suggesting leaky splicing at the intron 5 donor site. Sequencing of the PCR products confirmed the inclusion of intron 5 in the longer transcript, predicted to cause a premature stop codon. Thus, we provide evidence for a leaky splice site causing partial loss of desmin associated with a unique phenotypic presentation of a milder form of desmin-related recessive myopathy overlapping with congenital myasthenic syndrome.

TALEN-mediated intron editing of HSPCs enables transgene expression restricted to the myeloid lineage.

Gene therapy in hematopoietic stem and progenitor cells (HSPCs) shows great potential for the treatment of inborn metabolic diseases. Typical HSPC gene therapy approaches rely on constitutive promoters to express a therapeutic transgene, which is associated with multiple disadvantages. Here, we propose a novel promoterless intronic gene editing approach that triggers transgene expression only after cellular differentiation into the myeloid lineage. We integrated a splicing-competent eGFP cassette into the first intron of CD11b and observed expression of eGFP in the myeloid lineage but minimal to no expression in HSPCs or differentiated non-myeloid lineages. In vivo, edited HSPCs successfully engrafted in immunodeficient mice and displayed transgene expression in the myeloid compartment of multiple tissues. Using the same approach, we expressed alpha-L-iduronidase (IDUA), the defective enzyme in Mucopolysaccharidosis type I, and observed a 10-fold supraendogenous IDUA expression exclusively after myeloid differentiation. Edited cells efficiently populated bone marrow, blood, and spleen of immunodeficient mice, and retained the capacity to secrete IDUA ex vivo. Importantly, cells edited with the eGFP and IDUA transgenes were also found in the brain. This approach may unlock new therapeutic strategies for inborn metabolic and neurological diseases that require the delivery of therapeutics in brain.

Perturbation of the insomnia WDR90 genome-wide association studies locus pinpoints rs3752495 as a causal variant influencing distal expression of neighboring gene, PIG-Q.

Although genome-wide association studies (GWAS) have identified loci for sleep-related traits, they do not directly uncover the underlying causal variants and corresponding effector genes. The majority of such variants reside in non-coding regions and are therefore presumed to impact cis-regulatory elements. Our previously reported 'variant-to-gene mapping' effort in human induced pluripotent stem cell (iPSC)-derived neural progenitor cells (NPCs), combined with validation in both Drosophila and zebrafish, implicated phosphatidyl inositol glycan (PIG)-Q as a functionally relevant gene at the insomnia "WDR90" GWAS locus. However, importantly that effort did not characterize the corresponding underlying causal variant. Specifically, our previous 3D genomic datasets nominated a shortlist of three neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium within an intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. We sought to investigate the influence of these SNPs collectively and then individually on PIG-Q modulation to pinpoint the causal "regulatory" variant. Starting with gross level perturbation, deletion of the entire region in NPCs via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from individual luciferase reporter assays for each SNP in iPSCs revealed that the region with the rs3752495 risk allele (RA) induced a ~2.5-fold increase in luciferase expression. Importantly, rs3752495 also exhibited an allele-specific effect, with the RA increasing the luciferase expression by ~2-fold versus the non-RA. In conclusion, our variant-to-function approach and in vitro validation implicate rs3752495 as a causal insomnia variant embedded within WDR90 while modulating the expression of the distally located PIG-Q.

A concerted increase in readthrough and intron retention drives transposon expression during aging and senescence.

Aging and senescence are characterized by pervasive transcriptional dysfunction, including increased expression of transposons and introns. Our aim was to elucidate mechanisms behind this increased expression. Most transposons are found within genes and introns, with a large minority being close to genes. This raises the possibility that transcriptional readthrough and intron retention are responsible for age-related changes in transposon expression rather than expression of autonomous transposons. To test this, we compiled public RNA-seq datasets from aged human fibroblasts, replicative and drug-induced senescence in human cells, and RNA-seq from aging mice and senescent mouse cells. Indeed, our reanalysis revealed a correlation between transposons expression, intron retention, and transcriptional readthrough across samples and within samples. Both intron retention and readthrough increased with aging or cellular senescence and these transcriptional defects were more pronounced in human samples as compared to those of mice. In support of a causal connection between readthrough and transposon expression, analysis of models showing induced transcriptional readthrough confirmed that they also show elevated transposon expression. Taken together, our data suggest that elevated transposon reads during aging seen in various RNA-seq dataset are concomitant with multiple transcriptional defects. Intron retention and transcriptional readthrough are the most likely explanation for the expression of transposable elements that lack a functional promoter.

Identification of a novel ANK1 gene variant c.1504-9G>A and its mechanism of intron retention in hereditary spherocytosis.

Objective: The objective of this study was to pinpoint pathogenic genes and assess the mutagenic pathogenicity in two pediatric patients with hereditary spherocytosis. Methods: We utilized whole-exome sequencing (WES) for individual analysis (case 1) and family-based trio analysis (case 2). The significance of the intronic mutation was validated through a Minigene splicing assay and supported by subsequent in vitro experiments. Results: Both probands received a diagnosis of hereditary spherocytosis. WES identified a novel ANK1 c.1504-9G>A mutation in both patients, causing the retention of seven nucleotides at the 5' end of intron 13, as substantiated by the Minigene assay. This variant results in a premature stop codon and the production of a truncated protein. In vitro studies indicated a reduced expression of the ANK1 gene. Conclusion: The novel ANK1 c.1504-9G>A variant is established as the causative factor for hereditary spherocytosis, with the c.1504-9G site functioning as a splicing receptor.