Splice modulation strategy applied to deep intronic variants in COL7A1 causing recessive dystrophic epidermolysis bullosa.

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare and most often severe genetic disease characterized by recurrent blistering and erosions of the skin and mucous membranes after minor trauma, leading to major local and systemic complications. The disease is caused by loss-of-function variants in COL7A1 encoding type VII collagen (C7), the main component of anchoring fibrils, which form attachment structures stabilizing the cutaneous basement membrane zone. Alterations in C7 protein structure and/or expression lead to abnormal, rare or absent anchoring fibrils resulting in loss of dermal-epidermal adherence and skin blistering. To date, more than 1,200 distinct COL7A1 deleterious variants have been reported and 19% are splice variants. Here, we describe two RDEB patients for whom we identified two pathogenic deep intronic pathogenic variants in COL7A1. One of these variants (c.7795-97C > G) promotes the inclusion of a pseudoexon between exons 104 and 105 in the COL7A1 transcript, while the other causes partial or complete retention of intron 51. We used antisense oligonucleotide (ASO) mediated exon skipping to correct these aberrant splicing events in vitro. This led to increased normal mRNA splicing above 94% and restoration of C7 protein expression at a level (up to 56%) that should be sufficient to reverse the phenotype. This first report of exon skipping applied to counteract deep intronic variants in COL7A1 represents a promising therapeutic strategy for personalized medicine directed at patients with intronic variants at a distance of consensus splice sites.

MIPs: multi-locus intron polymorphisms in species identification and population genomics.

The study of species groups in which the presence of interspecific hybridization or introgression phenomena is known or suspected involves analysing shared bi-parentally inherited molecular markers. Current methods are based on different categories of markers among which the classical microsatellites or the more recent genome wide approaches for the analyses of thousands of SNPs or hundreds of microhaplotypes through high throughput sequencing. Our approach utilizes intron-targeted amplicon sequencing to characterise multi-locus intron polymorphisms (MIPs) and assess genetic diversity. These highly variable intron regions, combined with inter-specific transferable loci, serve as powerful multiple-SNP markers potentially suitable for various applications, from species and hybrid identification to population comparisons, without prior species knowledge. We developed the first panel of MIPs highly transferable across fish genomes, effectively distinguishing between species, even those closely related, and populations with different structures. MIPs offer versatile, hypervariable nuclear markers and promise to be especially useful when multiple nuclear loci must be genotyped across different species, such as for the monitoring of interspecific hybridization. Moreover, the relatively long sequences obtained ease the development of single-locus PCR-based diagnostic markers. This method, here demonstrated in teleost fishes, can be readily applied to other taxa, unlocking a new source of genetic variation.

Evolution of retrocopies in the context of HUSH silencing.

Retrotransposition is one of the main factors responsible for gene duplication and thus genome evolution. However, the sequences that undergo this process are not only an excellent source of biological diversity, but in certain cases also pose a threat to the integrity of the DNA. One of the mechanisms that protects against the incorporation of mobile elements is the HUSH complex, which is responsible for silencing long, intronless, transcriptionally active transposed sequences that are rich in adenine on the sense strand. In this study, broad sets of human and porcine retrocopies were analysed with respect to the above factors, taking into account evolution of these molecules. Analysis of expression pattern, genomic structure, transcript length, and nucleotide substitution frequency showed the strong relationship between the expression level and exon length as well as the protective nature of introns. The results of the studies also showed that there is no direct correlation between the expression level and adenine content. However, protein-coding retrocopies, which have a lower adenine content, have a significantly higher expression level than the adenine-rich non-coding but expressed retrocopies. Therefore, although the mechanism of HUSH silencing may be an important part of the regulation of retrocopy expression, it is one component of a more complex molecular network that remains to be elucidated.

A RORE-dependent Intronic Enhancer in the IL-7 Receptor-α Locus Controls Glucose Metabolism via Vγ4+ γδT17 Cells.

The IL-7R regulates the homeostasis, activation, and distribution of T cells in peripheral tissues. Although several transcriptional enhancers that regulate IL-7Rα expression in αß T cells have been identified, enhancers active in γδ T cells remain unknown. In this article, we discovered an evolutionarily conserved noncoding sequence (CNS) in intron 2 of the IL-7Rα-chain (IL-7Rα) locus and named this region CNS9. CNS9 contained a conserved retinoic acid receptor-related orphan receptor (ROR)-responsive element (RORE) and exerted RORγt-dependent enhancer activity in vitro. Mice harboring point mutations in the RORE in CNS9 (CNS9-RORmut) showed reduced IL-7Rα expression in IL-17-producing Vγ4+ γδ T cells. In addition, the cell number and IL-17A production of Vγ4+ γδ T cells were reduced in the adipose tissue of CNS9-RORmut mice. Consistent with the reduction in IL-17A, CNS9-RORmut mice exhibited decreased IL-33 expression in the adipose tissue, resulting in fewer regulatory T cells and glucose intolerance. The CNS9-ROR motif was partially responsible for IL-7Rα expression in RORγt+ regulatory T cells, whereas IL-7Rα expression was unaffected in RORγt-expressing Vγ2+ γδ T cells, Th17 cells, type 3 innate lymphoid cells, and invariant NKT cells. Our results indicate that CNS9 is a RORΕ-dependent, Vγ4+ γδ T cell-specific IL-7Rα enhancer that plays a critical role in adipose tissue homeostasis via regulatory T cells, suggesting that the evolutionarily conserved RORΕ in IL-7Rα intron 2 may influence the incidence of type 2 diabetes.

Investigating Splice Defects in USH2A Using Targeted Long-Read Sequencing.

Biallelic variants in USH2A are associated with retinitis pigmentosa (RP) and Type 2 Usher Syndrome (USH2), leading to impaired vision and, additionally, hearing loss in the latter. Although the introduction of next-generation sequencing into clinical diagnostics has led to a significant uplift in molecular diagnostic rates, many patients remain molecularly unsolved. It is thought that non-coding variants or variants of uncertain significance contribute significantly to this diagnostic gap. This study aims to demonstrate the clinical utility of the reverse transcription-polymerase chain reaction (RT-PCR)-Oxford Nanopore Technology (ONT) sequencing of USH2A mRNA transcripts from nasal epithelial cells to determine the splice-altering effect of candidate variants. Five affected individuals with USH2 or non-syndromic RP who had undergone whole genome sequencing were recruited for further investigation. All individuals had uncertain genotypes in USH2A, including deep intronic rare variants, c.8682-654C>G, c.9055+389G>A, and c.9959-2971C>T; a synonymous variant of uncertain significance, c.2139C>T; p.(Gly713=); and a predicted loss of function duplication spanning an intron/exon boundary, c.3812-3_3837dup p.(Met1280Ter). In silico assessment using SpliceAI provided splice-altering predictions for all candidate variants which were investigated using ONT sequencing. All predictions were found to be accurate; however, in the case of c.3812-3_3837dup, the outcome was a complex cryptic splicing pattern with predominant in-frame exon 18 skipping and a low level of exon 18 inclusion leading to the predicted stop gain. This study detected and functionally characterised simple and complex mis-splicing patterns in USH2A arising from previously unknown deep intronic variants and previously reported variants of uncertain significance, confirming the pathogenicity of the variants.

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.

A novel splice variant in intron 10 of PEX6 is associated with Zellweger Syndrome in a Chinese neonate.

BACKGROUND: Zellweger Syndrome (ZS), or cerebrohepatorenal syndrome, is a rare disorder due to PEX gene mutations affecting peroxisome function. While PEX6 coding mutations are known to cause ZS, the impact of noncoding mutations is less clear. METHODS: A Chinese neonate and his family were subjected to whole exome sequencing (WES) and bioinformatics to assess variant pathogenicity. A minigene assay was also performed for detailed splicing variant analysis. RESULTS: WES identified compound heterozygous PEX6 variants: c.315G>A (p. Trp105Ter) and c.2095-3 T>G. Minigene assays indicated that the latter variant led to abnormal mRNA splicing and the loss of exon 11 in PEX6 expression, potentially causing nonsense-mediated mRNA decay (NMD) or truncated protein structure. CONCLUSION: The study suggests that PEX6: c.2095-3 T>G might be a genetic contributor to the patient's condition, broadening the known mutation spectrum of PEX6. These insights lay groundwork for potential gene therapy for such variants.

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.

The Effect of Alternative Splicing Sites on Mirtron Formation and Arm Selection of Precursor microRNAs.

Mirtrons represent a subclass of microRNAs (miRNAs) that rely on the splicing machinery for their maturation. However, the molecular details of this Drosha-independent processing are still not fully understood; as an example, the Microprocessor complex cannot process the mirtronic pre-miRNA from the transcript even if splice site mutations are present. To investigate the influence of alternative splicing sites on mirtron formation, we generated Enhanced Green Fluorescent Protein (EGFP) reporters containing artificial introns to compare the processing of canonical miRNAs and mirtrons. Although mutations of both splice sites generated a complex pattern of alternative transcripts, mirtron formation was always severely affected as opposed to the normal processing of the canonical hsa-mir-33b miRNA. However, we also detected that while its formation was also hindered, the mirtron-derived hsa-mir-877-3p miRNA was less affected by certain mutations than the hsa-mir-877-5p species. By knocking down Drosha, we showed that this phenomenon is not dependent on Microprocessor activity but rather points toward the potential stability difference between the miRNAs from the different arms. Our results indicate that when the major splice sites are mutated, mirtron formation cannot be rescued by nearby alternative splice sites, and stability differences between 5p and 3p species should also be considered for functional studies of mirtrons.

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.

Rapid functional activation of horizontally transferred eukaryotic intron-containing genes in the bacterial recipient.

Horizontal gene transfer has occurred across all domains of life and contributed substantially to the evolution of both prokaryotes and eukaryotes. Previous studies suggest that many horizontally transferred eukaryotic genes conferred selective advantages to bacterial recipients, but how these eukaryotic genes evolved into functional bacterial genes remained unclear, particularly how bacteria overcome the expressional barrier posed by eukaryotic introns. Here, we first confirmed that the presence of intron would inactivate the horizontally transferred gene in Escherichia coli even if this gene could be efficiently transcribed. Subsequent large-scale genetic screens for activation of gene function revealed that activation events could rapidly occur within several days of selective cultivation. Molecular analysis of activation events uncovered two distinct mechanisms how bacteria overcome the intron barrier: (i) intron was partially deleted and the resulting stop codon-removed mutation led to one intact foreign protein or (ii) intron was intactly retained but it mediated the translation initiation and the interaction of two split small proteins (derived from coding sequences up- and downstream of intron, respectively) to restore gene function. Our findings underscore the likelihood that horizontally transferred eukaryotic intron-containing genes could rapidly acquire functionality if they confer a selective advantage to the prokaryotic recipient.

Functional Characterization of the Human BRCA1 ∆11 Splicing Isoforms in Yeast.

BRCA1, a crucial tumor suppressor gene, has several splicing isoforms, including Δ9-11, Δ11, and Δ11q, which lack exon 11, coding for significant portions of the protein. These isoforms are naturally present in both normal and cancerous cells, exhibiting altered activity compared to the full-length BRCA1. Despite this, the impact on cancer risk of the germline intronic variants promoting the exclusive expression of these Δ11 isoforms remains uncertain. Consequently, they are classified as variants of uncertain significance (VUS), posing challenges for traditional genetic classification methods due to their rarity and complexity. Our research utilizes a yeast-based functional assay, previously validated for assessing missense BRCA1 variants, to compare the activity of the Δ11 splicing isoforms with known pathogenic missense variants. This approach allows us to elucidate the functional implications of these isoforms and determine whether their exclusive expression could contribute to increased cancer risk. By doing so, we aim to provide insights into the pathogenic potential of intronic VUS-generating BRCA1 splicing isoforms and improve the classification of BRCA1 variants.

Clinical and Molecular Characterization of SMAD4 Splicing Variants in Patients with Juvenile Polyposis Syndrome.

Juvenile polyposis syndrome (JPS) is an inherited autosomal dominant condition that predisposes to the development of juvenile polyps throughout the gastrointestinal (GI) tract, and it poses an increased risk of GI malignancy. Germline causative variants were identified in the SMAD4 gene in a subset (20%) of JPS cases. Most SMAD4 germline genetic variants published to date are missense, nonsense, and frameshift mutations. SMAD4 germline alterations predicted to result in aberrant splicing have rarely been reported. Here, we report two unrelated Italian families harboring two different SMAD4 intronic variants, c.424+5G>A and c.425-9A>G, which are clinically associated with colorectal cancer and/or juvenile GI polyps. In silico prediction analysis, in vitro minigene assays, and RT-PCR showed that the identified variants lead to aberrant SMAD4 splicing via the exonization of intronic nucleotides, resulting in a premature stop codon. This is expected to cause the production of a truncated protein. This study expands the landscape of SMAD4 germline genetic variants associated with GI polyposis and/or cancer. Moreover, it emphasizes the importance of the functional characterization of SMAD4 splicing variants through RNA analysis, which can provide new insights into genetic disease variant interpretation, enabling tailored genetic counseling, management, and surveillance of patients with GI polyposis and/or cancer.

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.

A sequence of SVA retrotransposon insertions in ASIP shaped human pigmentation.

Retrotransposons comprise about 45% of the human genome1, but their contributions to human trait variation and evolution are only beginning to be explored2,3. Here, we find that a sequence of SVA retrotransposon insertions in an early intron of the ASIP (agouti signaling protein) gene has probably shaped human pigmentation several times. In the UK Biobank (n = 169,641), a recent 3.3-kb SVA insertion polymorphism associated strongly with lighter skin pigmentation (0.22 [0.21-0.23] s.d.; P = 2.8 × 10-351) and increased skin cancer risk (odds ratio = 1.23 [1.18-1.27]; P = 1.3 × 10-28), appearing to underlie one of the strongest common genetic influences on these phenotypes within European populations4-6. ASIP expression in skin displayed the same association pattern, with the SVA insertion allele exhibiting 2.2-fold (1.9-2.6) increased expression. This effect had an unusual apparent mechanism: an earlier, nonpolymorphic, human-specific SVA retrotransposon 3.9 kb upstream appeared to have caused ASIP hypofunction by nonproductive splicing, which the new (polymorphic) SVA insertion largely eliminated. Extended haplotype homozygosity indicated that the insertion allele has risen to allele frequencies up to 11% in European populations over the past several thousand years. These results indicate that a sequence of retrotransposon insertions contributed to a species-wide increase, then a local decrease, of human pigmentation.

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.

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.

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.