An intron endonuclease facilitates interference competition between coinfecting viruses.

Introns containing homing endonucleases are widespread in nature and have long been assumed to be selfish elements that provide no benefit to the host organism. These genetic elements are common in viruses, but whether they confer a selective advantage is unclear. In this work, we studied intron-encoded homing endonuclease gp210 in bacteriophage ΦPA3 and found that it contributes to viral competition by interfering with the replication of a coinfecting phage, ΦKZ. We show that gp210 targets a specific sequence in ΦKZ, which prevents the assembly of progeny viruses. This work demonstrates how a homing endonuclease can be deployed in interference competition among viruses and provide a relative fitness advantage. Given the ubiquity of homing endonucleases, this selective advantage likely has widespread evolutionary implications in diverse plasmid and viral competition as well as virus-host interactions.

Dynamic evolution of the heterochromatin sensing histone demethylase IBM1.

Heterochromatin is critical for maintaining genome stability, especially in flowering plants, where it relies on a feedback loop involving the H3K9 methyltransferase, KRYPTONITE (KYP), and the DNA methyltransferase CHROMOMETHYLASE3 (CMT3). The H3K9 demethylase INCREASED IN BONSAI METHYLATION 1 (IBM1) counteracts the detrimental consequences of KYP-CMT3 activity in transcribed genes. IBM1 expression in Arabidopsis is uniquely regulated by methylation of the 7th intron, allowing it to monitor global H3K9me2 levels. We show the methylated intron is prevalent across flowering plants and its underlying sequence exhibits dynamic evolution. We also find extensive genetic and expression variations in KYP, CMT3, and IBM1 across flowering plants. We identify Arabidopsis accessions resembling weak ibm1 mutants and Brassicaceae species with reduced IBM1 expression or deletions. Evolution towards reduced IBM1 activity in some flowering plants could explain the frequent natural occurrence of diminished or lost CMT3 activity and loss of gene body DNA methylation, as cmt3 mutants in A. thaliana mitigate the deleterious effects of IBM1.

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.

Reassessing the exon-foldon correspondence using frustration analysis.

Protein folding and evolution are intimately linked phenomena. Here, we revisit the concept of exons as potential protein folding modules across a set of 38 abundant and conserved protein families. Taking advantage of genomic exon-intron organization and extensive protein sequence data, we explore exon boundary conservation and assess the foldon-like behavior of exons using energy landscape theoretic measurements. We found deviations in the exon size distribution from exponential decay indicating selection in evolution. We show that when taken together there is a pronounced tendency to independent foldability for segments corresponding to the more conserved exons, supporting the idea of exon-foldon correspondence. While 45% of the families follow this general trend when analyzed individually, there are some families for which other stronger functional determinants, such as preserving frustrated active sites, may be acting. We further develop a systematic partitioning of protein domains using exon boundary hotspots, showing that minimal common exons correspond with uninterrupted alpha and/or beta elements for the majority of the families but not for all of them.

IFIH1 (MDA5) is required for innate immune detection of intron-containing RNA expressed from the HIV-1 provirus.

Intron-containing RNA expressed from the HIV-1 provirus activates type 1 interferon in primary human blood cells, including CD4+ T cells, macrophages, and dendritic cells. To identify the innate immune receptor required for detection of intron-containing RNA expressed from the HIV-1 provirus, a loss-of-function screen was performed with short hairpin RNA-expressing lentivectors targeting twenty-one candidate genes in human monocyte-derived dendritic cells. Among the candidate genes tested, only knockdown of XPO1 (CRM1), IFIH1 (MDA5), or MAVS prevented activation of the interferon-stimulated gene ISG15. The importance of IFIH1 protein was demonstrated by rescue of the knockdown with nontargetable IFIH1 coding sequence. Inhibition of HIV-1-induced ISG15 by the IFIH1-specific Nipah virus V protein, and by IFIH1-transdominant 2-CARD domain-deletion or phosphomimetic point mutations, indicates that IFIH1 (MDA5) filament formation, dephosphorylation, and association with MAVS are all required for innate immune activation in response to HIV-1 transduction. Since both IFIH1 (MDA5) and DDX58 (RIG-I) signal via MAVS, the specificity of HIV-1 RNA detection by IFIH1 was demonstrated by the fact that DDX58 knockdown had no effect on activation. RNA-Seq showed that IFIH1 knockdown in dendritic cells globally disrupted the induction of IFN-stimulated genes by HIV-1. Finally, specific enrichment of unspliced HIV-1 RNA by IFIH1 (MDA5), over two orders of magnitude, was revealed by formaldehyde cross-linking immunoprecipitation (f-CLIP). These results demonstrate that IFIH1 is the innate immune receptor for intron-containing RNA from the HIV-1 provirus and that IFIH1 potentially contributes to chronic inflammation in people living with HIV-1, even in the presence of effective antiretroviral therapy.

Antagonistic conflict between transposon-encoded introns and guide RNAs.

TnpB nucleases represent the evolutionary precursors to CRISPR-Cas12 and are widespread in all domains of life. IS605-family TnpB homologs function as programmable RNA-guided homing endonucleases in bacteria, driving transposon maintenance through DNA double-strand break-stimulated homologous recombination. In this work, we uncovered molecular mechanisms of the transposition life cycle of IS607-family elements that, notably, also encode group I introns. We identified specific features for a candidate "IStron" from Clostridium botulinum that allow the element to carefully control the relative levels of spliced products versus functional guide RNAs. Our results suggest that IStron transcripts evolved an ability to balance competing and mutually exclusive activities that promote selfish transposon spread while limiting adverse fitness costs on the host. Collectively, this work highlights molecular innovation in the multifunctional utility of transposon-encoded noncoding RNAs.

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.

Identification of the new allele HLA-C*04:01:01:186 in a Greek individual using next generation sequencing.

The newly discovered HLA-C*04:01:01:186 allele differs from HLA-C*04:01:01:01 by a single nucleotide substitution in intron 3.

Characterisation of the novel HLA-A*32:01:01:41 allele by next-generation sequencing.

HLA-A*32:01:01:41 differs from the HLA-A*32:01:01:01 allele by one nucleotide substitution in the intron 3.

Unique features of conventional and nonconventional introns in Euglena gracilis.

BACKGROUND: Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns. RESULTS: A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans. CONCLUSIONS: Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Genome-wide identification and evolution of the tubulin gene family in Camelina sativa.

BACKGROUND: Tubulins play crucial roles in numerous fundamental processes of plant development. In flowering plants, tubulins are grouped into α-, ß- and γ-subfamilies, while α- and ß-tubulins possess a large isotype diversity and gene number variations among different species. This circumstance leads to insufficient recognition of orthologous isotypes and significantly complicates extrapolation of obtained experimental results, and brings difficulties for the identification of particular tubulin isotype function. The aim of this research is to identify and characterize tubulins of an emerging biofuel crop Camelina sativa. RESULTS: We report comprehensive identification and characterization of tubulin gene family in C. sativa, including analyses of exon-intron organization, duplicated genes comparison, proper isotype designation, phylogenetic analysis, and expression patterns in different tissues. 17 α-, 34 ß- and 6 γ-tubulin genes were identified and assigned to a particular isotype. Recognition of orthologous tubulin isotypes was cross-referred, involving data of phylogeny, synteny analyses and genes allocation on reconstructed genomic blocks of Ancestral Crucifer Karyotype. An investigation of expression patterns of tubulin homeologs revealed the predominant role of N6 (A) and N7 (B) subgenomes in tubulin expression at various developmental stages, contrarily to general the dominance of transcripts of H7 (C) subgenome. CONCLUSIONS: For the first time a complete set of tubulin gene family members was identified and characterized for allohexaploid C. sativa species. The study demonstrates the comprehensive approach of precise inferring gene orthology. The applied technique allowed not only identifying C. sativa tubulin orthologs in model Arabidopsis species and tracking tubulin gene evolution, but also uncovered that A. thaliana is missing orthologs for several particular isotypes of α- and ß-tubulins.

Detection of the novel HLA-DRB5*01:01:01:07 allele by Pacific Biosciences HiFi sequencing in a Chinese individual.

HLA-DRB5*01:01:01:07 differs from HLA-DRB5*01:01:01:01 by two nucleotide changes in intron 1 and intron 2.

Genome and RNA sequencing were essential to reveal cryptic intronic variants associated to defective ATP6AP1 mRNA processing.

The diagnosis of Mendelian disorders has notably advanced with integration of whole exome and genome sequencing (WES and WGS) in clinical practice. However, challenges in variant interpretation and uncovered variants by WES still leave a substantial percentage of patients undiagnosed. In this context, integrating RNA sequencing (RNA-seq) improves diagnostic workflows, particularly for WES inconclusive cases. Additionally, functional studies are often necessary to elucidate the impact of prioritized variants on gene expression and protein function. Our study focused on three unrelated male patients (P1-P3) with ATP6AP1-CDG (congenital disorder of glycosylation), presenting with intellectual disability and varying degrees of hepatopathy, glycosylation defects, and an initially inconclusive diagnosis through WES. Subsequent RNA-seq was pivotal in identifying the underlying genetic causes in P1 and P2, detecting ATP6AP1 underexpression and aberrant splicing. Molecular studies in fibroblasts confirmed these findings and identified the rare intronic variants c.289-233C > T and c.289-289G > A in P1 and P2, respectively. Trio-WGS also revealed the variant c.289-289G > A in P3, which was a de novo change in both patients. Functional assays expressing the mutant alleles in HAP1 cells demonstrated the pathogenic impact of these variants by reproducing the splicing alterations observed in patients. Our study underscores the role of RNA-seq and WGS in enhancing diagnostic rates for genetic diseases such as CDG, providing new insights into ATP6AP1-CDG molecular bases by identifying the first two deep intronic variants in this X-linked gene. Additionally, our study highlights the need to integrate RNA-seq and WGS, followed by functional validation, in routine diagnostics for a comprehensive evaluation of patients with an unidentified molecular etiology.

Identification and characterization of a novel intronic splicing mutation in CSF1R-related leukoencephalopathy.

AIMS: Colony stimulating factor 1 receptor (CSF1R)-related leukoencephalopathy is a rapidly progressing neurodegenerative disease caused by CSF1R gene mutations. This study aimed to identify and investigate the effect of a novel intronic mutation (c.1754-3C>G) of CSF1R on splicing. METHODS: A novel intronic mutation was identified using whole-exome sequencing. To investigate the impact of this mutation, we employed various bioinformatics tools to analyze the transcription of the CSF1R gene and the three-dimensional structure of its encoded protein. Furthermore, reverse transcription polymerase chain reaction (RT-PCR) was performed to validate the findings. RESULTS: A novel mutation (c.1754-3C>G) in CSF1R was identified, which results in exon 13 skipping due to the disruption of the 3' splice site consensus sequence NYAG/G. This exon skipping event was further validated in the peripheral blood of the mutation carrier through RT-PCR and Sanger sequencing. Protein structure prediction indicated a disruption in the tyrosine kinase domain, with the truncated protein showing significant structural alterations. CONCLUSIONS: Our findings underscore the importance of intronic mis-splicing mutations in the diagnosis and management of CSF1R-related leukoencephalopathy.

Light regulates nuclear detainment of intron-retained transcripts through COP1-spliceosome to modulate photomorphogenesis.

Intron retention (IR) is the most common alternative splicing event in Arabidopsis. An increasing number of studies have demonstrated the major role of IR in gene expression regulation. The impacts of IR on plant growth and development and response to environments remain underexplored. Here, we found that IR functions directly in gene expression regulation on a genome-wide scale through the detainment of intron-retained transcripts (IRTs) in the nucleus. Nuclear-retained IRTs can be kept away from translation through this mechanism. COP1-dependent light modulation of the IRTs of light signaling genes, such as PIF4, RVE1, and ABA3, contribute to seedling morphological development in response to changing light conditions. Furthermore, light-induced IR changes are under the control of the spliceosome, and in part through COP1-dependent ubiquitination and degradation of DCS1, a plant-specific spliceosomal component. Our data suggest that light regulates the activity of the spliceosome and the consequent IRT nucleus detainment to modulate photomorphogenesis through COP1.

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.

Comparative analysis the chloroplast genomes of Celastrus (Celastraceae) species: Provide insights into molecular evolution, species identification and phylogenetic relationships.

BACKGROUND: The genus Celastrus is an important medicinal plant resource. The similarity of morphology and the lack of complete chloroplast genome analysis have significantly impeded the exploration of species identification, molecular evolution and phylogeny of Celastrus. PURPOSE: In order to resolve the phylogenic controversy of Celastrus species, the chloroplast genome comparative analysis was performed to provide genetic evidence. METHODS: In this study, we collected and sequenced ten chloroplast genomes of Celastrus species from China and downloaded three chloroplast genomes from the databases. The chloroplast genomes were compared and analyzed to explore their characteristics and evolution. Furthermore, the phylogenetic relationships of Celastrus species were inferred based on the whole chloroplast genomes and protein-coding genes. RESULTS: All the 13 Celastrus species chloroplast genomes showed a typical quadripartite structure with genome sizes ranging from 155,113 to 157,366 bp. The intron loss of the rps16 gene occurred in all the 13 Celastrus species. The GC content, gene sequence, repeat types and codon bias pattern were highly conserved. Ten highly variation regions were identified, which can be used as potential DNA markers in molecular identification of Celastrus species. Eight genes, including accD, atp4, ndhB, rpoC1, rbcL, rpl2, rpl20 and ycf1, were detected to experience positive selection. Phylogenetic analysis showed that Celastrus was a monophyletic group and Tripterygium was the closest sister-group. Noteworthy, C. gemmatus Loes. and C. orbiculatus Thunb. can be discriminated using the chloroplast genome as a super barcode. The comparative and phylogenetic analysis results proposed that C. tonkinensis Pitard. was the synonym of C. hindsii Benth. CONCLUSION: The comparative analysis of the Celastrus chloroplast genomes can provide comprehensive genetic evidence for molecular evolution, species identification and phylogenetic relationships.

Targeting the conserved active site of splicing machines with specific and selective small molecule modulators.

The self-splicing group II introns are bacterial and organellar ancestors of the nuclear spliceosome and retro-transposable elements of pharmacological and biotechnological importance. Integrating enzymatic, crystallographic, and simulation studies, we demonstrate how these introns recognize small molecules through their conserved active site. These RNA-binding small molecules selectively inhibit the two steps of splicing by adopting distinctive poses at different stages of catalysis, and by preventing crucial active site conformational changes that are essential for splicing progression. Our data exemplify the enormous power of RNA binders to mechanistically probe vital cellular pathways. Most importantly, by proving that the evolutionarily-conserved RNA core of splicing machines can recognize small molecules specifically, our work provides a solid basis for the rational design of splicing modulators not only against bacterial and organellar introns, but also against the human spliceosome, which is a validated drug target for the treatment of congenital diseases and cancers.

Red2Flpe-SCON: a versatile, multicolor strategy for generating mosaic conditional knockout mice.

Image-based lineage tracing enables tissue turnover kinetics and lineage potentials of different adult cell populations to be investigated. Previously, we reported a genetic mouse model system, Red2Onco, which ectopically expressed mutated oncogenes together with red fluorescent proteins (RFP). This system enabled the expansion kinetics and neighboring effects of oncogenic clones to be dissected. We now report Red2Flpe-SCON: a mosaic knockout system that uses multicolor reporters to label both mutant and wild-type cells. We develop the Red2Flpe mouse line for red clone-specific Flpe expression, as well as the FRT-based SCON (Short Conditional IntrON) method to facilitate tunable conditional mosaic knockouts in mice. We use the Red2Flpe-SCON method to study Sox2 mutant clonal analysis in the esophageal epithelium of adult mice which reveal that the stem cell gene, Sox2, is less essential for adult stem cell maintenance itself, but rather for stem cell proliferation and differentiation.

A systematic screen identifies Saf5 as a link between splicing and transcription in fission yeast.

Splicing is an important step of gene expression regulation in eukaryotes, as there are many mRNA precursors that can be alternatively spliced in different tissues, at different cell cycle phases or under different external stimuli. We have developed several integrated fluorescence-based in vivo splicing reporter constructs that allow the quantification of fission yeast splicing in vivo on intact cells, and we have compared their splicing efficiency in a wild type strain and in a prp2-1 (U2AF65) genetic background, showing a clear dependency between Prp2 and a consensus signal at 5' splicing site (5'SS). To isolate novel genes involved in regulated splicing, we have crossed the reporter showing more intron retention with the Schizosaccharomyces pombe knock out collection. Among the candidate genes involved in the regulation of splicing, we have detected strong splicing defects in two of the mutants -Δcwf12, a member of the NineTeen Complex (NTC) and Δsaf5, a methylosome subunit that acts together with the survival motor neuron (SMN) complex in small nuclear ribonucleoproteins (snRNP) biogenesis. We have identified that strains with mutations in cwf12 have inefficient splicing, mainly when the 5'SS differs from the consensus. However, although Δsaf5 cells also have some dependency on 5'SS sequence, we noticed that when one intron of a given pre-mRNA was affected, the rest of the introns of the same pre-mRNA had high probabilities of being also affected. This observation points Saf5 as a link between transcription rate and splicing.