How did alternative splicing evolve?

An alternative 3' splice site of PeuHKT1;3 improves the response to salt stress through enhancing affinity to K+ in Populus

Alternative splicing (AS) serves as a crucial post-transcriptional regulator in plants that contributes to the resistance to salt stress. However, the underlying mechanism is largely unknown. In this research, we identified an important AS transcript in Populus euphratica, PeuHKT1:3a, generated by alternative 3' splice site splicing mode that resulted in the removal of 252 bases at the 5' end of the first exon in PeuHKT1:3. Protein sequence comparison showed that the site of AS occurred in PeuHKT1:3 is located at a crucial Ser residue within the first pore-loop domain, which leads to inefficient K+ transport in HKT I-type transporters. Expressing PeuHKT1;3a in an axt3 mutant yeast strain can effectively compensate for the lack of intracellular K+, whereas the expression of PeuHKT1;3 cannot yield the effect. Furthermore, in transgenic Arabidopsis and poplar plants, it was observed that lines expressing PeuHKT1;3a exhibited greater salt tolerance compared to those expressing the PeuHKT1;3 strain. Analysis of ion content and flux demonstrated that the transgenic PeuHKT1;3a line exhibited significantly higher K+ content compared to the PeuHKT1;3 line, while there was no significant difference in Na+ content. In conclusion, our findings revealed that AS can give rise to novel variants of HKT I-type proteins in P. euphratica with modified K+ selectivity to keep a higher K+/Na+ ratio to enhanced salt tolerance.

Age-dependent regulation of ELP1 exon 20 splicing in Familial Dysautonomia by RNA Polymerase II kinetics and chromatin structure

Familial Dysautonomia (FD) is a rare disease caused by ELP1 exon 20 skipping. Here we clarify the role of RNA Polymerase II (RNAPII) and chromatin on this splicing event. A slow RNAPII mutant and chromatin-modifying chemicals that reduce the rate of RNAPII elongation induce exon skipping whereas chemicals that create a more relaxed chromatin exon inclusion. In the brain of a mouse transgenic for the human FD-ELP1 we observed on this gene an age-dependent decrease in the RNAPII density profile that was most pronounced on the alternative exon, a robust increase in the repressive marks H3K27me3 and H3K9me3 and a decrease of H3K27Ac, together with a progressive reduction in ELP1 exon 20 inclusion level. In HEK 293T cells, selective drug-induced demethylation of H3K27 increased RNAPII elongation on ELP1 and SMN2, promoted the inclusion of the corresponding alternative exons, and, by RNA-sequencing analysis, induced changes in several alternative splicing events. These data suggest a co-transcriptional model of splicing regulation in which age-dependent changes in H3K27me3/Ac modify the rate of RNAPII elongation and affect processing of ELP1 alternative exon 20.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

Regulation of Hair Follicle Growth and Development by Different Alternative Spliceosomes of FGF5 in Rabbits

This study investigated the regulatory effect of alternative spliceosomes of the fibroblast growth factor 5 (FGF5) gene on hair follicle (HF) growth and development in rabbits. The FGF5 alternative spliceosomes (called FGF5-X1, FGF5-X2, FGF5-X3) were cloned. The overexpression vector and siRNA of spliceosomes were transfected into dermal papilla cells (DPCs) to analyze the regulatory effect on DPCs. The results revealed that FGF5-X2 and FGF5-X3 overexpression significantly decreased LEF1 mRNA expression (p < 0.01). FGF5-X1 overexpression significantly reduced CCND1 expression (p < 0.01). FGF5-X1 and FGF5-X2 possibly downregulated the expression level of FGF2 mRNA (p < 0.05), and FGF5-X3 significantly downregulated the expression level of FGF2 mRNA (p < 0.01). The FGF5 alternative spliceosomes significantly downregulated the BCL2 mRNA expression level in both cases (p < 0.01). FGF5-X1 and FGF5-X2 significantly increased TGFβ mRNA expression (p < 0.01). All three FGF5 alternative spliceosomes inhibited DPC proliferation. In conclusion, the expression profile of HF growth and development-related genes can be regulated by FGF5 alternative spliceosomes, inhibiting the proliferation of DPCs and has an influence on the regulation of HF growth in rabbits. This study provides insights to further investigate the mechanism of HF development in rabbits via FGF5 regulation.

The impact of genetically controlled splicing on exon inclusion and protein structure

Common variants affecting mRNA splicing are typically identified though splicing quantitative trait locus (sQTL) mapping and have been shown to be enriched for GWAS signals by a similar degree to eQTLs. However, the specific splicing changes induced by these variants have been difficult to characterize, making it more complicated to analyze the effect size and direction of sQTLs, and to determine downstream splicing effects on protein structure. In this study, we catalogue sQTLs using exon percent spliced in (PSI) scores as a quantitative phenotype. PSI is an interpretable metric for identifying exon skipping events and has some advantages over other methods for quantifying splicing from short read RNA sequencing. In our set of sQTL variants, we find evidence of selective effects based on splicing effect size and effect direction, as well as exon symmetry. Additionally, we utilize AlphaFold2 to predict changes in protein structure associated with sQTLs overlapping GWAS traits, highlighting a potential new use-case for this technology for interpreting genetic effects on traits and disorders.

The Transcription Factor StuA Regulates the Glyoxylate Cycle in the Dermatophyte Trichophyton rubrum under Carbon Starvation

Trichophyton rubrum is the primary causative agent of dermatophytosis worldwide. This fungus colonizes keratinized tissues and uses keratin as a nutritional source during infection. In T. rubrum-host interactions, sensing a hostile environment triggers the adaptation of its metabolic machinery to ensure its survival. The glyoxylate cycle has emerged as an alternative metabolic pathway when glucose availability is limited; this enables the conversion of simple carbon compounds into glucose via gluconeogenesis. In this study, we investigated the impact of stuA deletion on the response of glyoxylate cycle enzymes during fungal growth under varying culture conditions in conjunction with post-transcriptional regulation through alternative splicing of the genes encoding these enzymes. We revealed that the ΔstuA mutant downregulated the malate synthase and isocitrate lyase genes in a keratin-containing medium or when co-cultured with human keratinocytes. Alternative splicing of an isocitrate lyase gene yielded a new isoform. Enzymatic activity assays showed specific instances where isocitrate lyase and malate synthase activities were affected in the mutant strain compared to the wild type strain. Taken together, our results indicate a relevant balance in transcriptional regulation that has distinct effects on the enzymatic activities of malate synthase and isocitrate lyase.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

mRNA splicing is modulated by intronic microRNAs

Splicing of transcripts is catalyzed by the spliceosome, a mega-complex consisting of hundreds of proteins and five snRNAs, which employs direct interactions. When U1 snRNA forms high-affinity binding, namely more than eight base pairs, with the 5'SS, the result is usually a suppressing effect on the splicing activity. This likely occurs due to the inefficient unwinding of U1/5'SS base-pairing or other regulatory obstructions. Here, we show in vitro and in patient-derived cell lines that pre-microRNAs can modulate the splicing reaction by interacting with U1 snRNA. This leads to reduced binding affinity to the 5'SS, and hence promotes the inclusion of exons containing 5'SS, despite sequence-based high affinity to U1. Application of the mechanism resulted in correction of the splicing defect in the disease-causing VCAN gene from an individual with Wagner syndrome. This pre-miRNA/U1 interaction can regulate the expression of alternatively spliced exons, thus extending the scope of mechanisms regulating splicing.

Conserved and divergent signals in 5' splice site sequences across fungi, metazoa and plants

In eukaryotic organisms the ensemble of 5' splice site sequences reflects the balance between natural nucleotide variability and minimal molecular constraints necessary to ensure splicing fidelity. This compromise shapes the underlying statistical patterns in the composition of donor splice site sequences. The scope of this study was to mine conserved and divergent signals in the composition of 5' splice site sequences. Because 5' donor sequences are a major cue for proper recognition of splice sites, we reasoned that statistical regularities in their composition could reflect the biological functionality and evolutionary history associated with splicing mechanisms. Results: We considered a regularized maximum entropy modeling framework to mine for non-trivial two-site correlations in donor sequence datasets corresponding to 30 different eukaryotes. For each analyzed species, we identified minimal sets of two-site coupling patterns that were able to replicate, at a given regularization level, the observed one-site and two-site frequencies in donor sequences. By performing a systematic and comparative analysis of 5'splice sites we showed that lineage information could be traced from joint di-nucleotide probabilities. We were able to identify characteristic two-site coupling patterns for plants and animals, and propose that they may echo differences in splicing regulation previously reported between these groups.

Sequence variations affect the 5' splice site selection of plant introns

Introns are noncoding sequences spliced out of pre-mRNAs by the spliceosome to produce mature mRNAs. The 5' ends of introns mostly begin with GU and have a conserved sequence motif of AG/GUAAGU that could base-pair with the core sequence of U1 snRNA of the spliceosome. Intriguingly, ∼ 1% of introns in various eukaryotic species begin with GC. This occurrence could cause misannotation of genes; however, the underlying splicing mechanism is unclear. We analyzed the sequences around the intron 5' splice site (ss) in Arabidopsis (Arabidopsis thaliana) and found sequences at the GC intron ss are much more stringent than those of GT introns. Mutational analysis at various positions of the intron 5' ss revealed that although mutations impair base pairing, different mutations at the same site can have different effects, suggesting that steric hindrance also affects splicing. Moreover, mutations of 5' ss often activate a hidden ss nearby. Our data suggest that the 5' ss is selected via a competition between the major ss and the nearby minor ss. This work not only provides insights into the splicing mechanism of intron 5' ss but also improves the accuracy of gene annotation and the study of the evolution of intron 5' ss.

Building alternative splicing and evolution-aware sequence-structure maps for protein repeats

Alternative splicing of repeats in proteins provides a mechanism for rewiring and fine-tuning protein interaction networks. In this work, we developed a robust and versatile method, ASPRING, to identify alternatively spliced protein repeats from gene annotations. ASPRING leverages evolutionary meaningful alternative splicing-aware hierarchical graphs to provide maps between protein repeats sequences and 3D structures. We re-think the definition of repeats by explicitly accounting for transcript diversity across several genes/species. Using a stringent sequence-based similarity criterion, we detected over 5,000 evolutionary conserved repeats by screening virtually all human protein-coding genes and their orthologs across a dozen species. Through a joint analysis of their sequences and structures, we extracted specificity-determining sequence signatures and assessed their implication in experimentally resolved and modelled protein interactions. Our findings demonstrate the widespread alternative usage of protein repeats in modulating protein interactions and open avenues for targeting repeat-mediated interactions.

The effects of mutations on gene expression and alternative splicing

Understanding the relationship between mutations and their genomic and phenotypic consequences has been a longstanding goal of evolutionary biology. However, few studies have investigated the impact of mutations on gene expression and alternative splicing on the genome-wide scale. In this study, we aim to bridge this knowledge gap by utilizing whole-genome sequencing data and RNA sequencing data from 16 obligately parthenogenetic Daphnia mutant lines to investigate the effects of ethyl methanesulfonate-induced mutations on gene expression and alternative splicing. Using rigorous analyses of mutations, expression changes and alternative splicing, we show that trans-effects are the major contributor to the variance in gene expression and alternative splicing between the wild-type and mutant lines, whereas cis mutations only affected a limited number of genes and do not always alter gene expression. Moreover, we show that there is a significant association between differentially expressed genes and exonic mutations, indicating that exonic mutations are an important driver of altered gene expression.

Characterization of Glutamate-Gated Chloride Channel in Tribolium castaneum

The glutamate-gated chloride channels (GluCls) play essential roles in signal transduction by regulating fast inhibitory synaptic transmission in the nervous system of invertebrates. While there is only one GluCl subunit in the insect, the diversity of insect GluCls is broadened by alternative splicing. In the present study, three TcGluCl variant genes were cloned from the red flour beetle Tribolium castaneum. Analysis of the characteristics of TcGluCls including sequence features, genomic structures, and alternative splicing revealed that TcGluCls had the typical structural features of GluCls and showed high homologies with the GluCls from other insect orders. The TcGluCl-encoding gene consists of nine exons and three variants (TcGluCl-3a, TcGluCl-3b, and TcGluCl-3c) were generated by the alternative splicing of exon 3, which was a highly conserved alternative splicing site in insect GluCls. Homology modeling of TcGluCl-3a showed that the exon 3 coding protein located at the N-terminal extracellular domain, and there were no steric clashes encountered between the exon 3 coding region and ivermectin/glutamate binding pocket, which indicated that the alternative splicing of exon 3 might have no impact on the binding of GluCls to glutamate and insecticide. In addition to the head tissue, TcGluCl-3a and TcGluCl-3c also had high expressions in the ovary and testis of T. castaneum, whereas TcGluCl-3b showed high expression in the midgut, suggesting the diverse physiological functions of TcGluCl variants in T. castaneum. The total TcGluCl and three variants showed the highest expression levels in the early stage larvae. The expressions of TcGluCl, TcGluCl-3b, and TcGluCl-3c were significantly increased from the late-stage larvae to the early stage pupae and indicated that the TcGluCl might be involved in the growth and development of T. castaneum. These results are helpful to further understand the molecular characteristics of insect GluCls and provide foundations for studying the specific function of the GluCl variant.

PCBP1 regulates alternative splicing of AARS2 in congenital cardiomyopathy

Alanyl-transfer RNA synthetase 2 (AARS2) is a nuclear encoded mitochondrial tRNA synthetase that is responsible for charging of tRNA-Ala with alanine during mitochondrial translation. Homozygous or compound heterozygous mutations in the Aars2 gene, including those affecting its splicing, are linked to infantile cardiomyopathy in humans. However, how Aars2 regulates heart development, and the underlying molecular mechanism of heart disease remains unknown. Here, we found that poly(rC) binding protein 1 (PCBP1) interacts with the Aars2 transcript to mediate its alternative splicing and is critical for the expression and function of Aars2. Cardiomyocyte-specific deletion of Pcbp1 in mice resulted in defects in heart development that are reminiscent of human congenital cardiac defects, including noncompaction cardiomyopathy and a disruption of the cardiomyocyte maturation trajectory. Loss of Pcbp1 led to an aberrant alternative splicing and a premature termination of Aars2 in cardiomyocytes. Additionally, Aars2 mutant mice with exon-16 skipping recapitulated heart developmental defects observed in Pcbp1 mutant mice. Mechanistically, we found dysregulated gene and protein expression of the oxidative phosphorylation pathway in both Pcbp1 and Aars2 mutant hearts; these date provide further evidence that the infantile hypertrophic cardiomyopathy associated with the disorder oxidative phosphorylation defect type 8 (COXPD8) is mediated by Aars2. Our study therefore identifies Pcbp1 and Aars2 as critical regulators of heart development and provides important molecular insights into the role of disruptions in metabolism on congenital heart defects.

LncRNA FLAIL affects alternative splicing and represses flowering in Arabidopsis

How the noncoding genome affects cellular functions is a key biological question. A particular challenge is to distinguish the effects of noncoding DNA elements from long noncoding RNAs (lncRNAs) that coincide at the same loci. Here, we identified the flowering-associated intergenic lncRNA (FLAIL) in Arabidopsis through early flowering flail mutants. Expression of FLAIL RNA from a different chromosomal location in combination with strand-specific RNA knockdown characterized FLAIL as a trans-acting RNA molecule. FLAIL directly binds to differentially expressed target genes that control flowering via RNA-DNA interactions through conserved sequence motifs. FLAIL interacts with protein and RNA components of the spliceosome to affect target mRNA expression through co-transcriptional alternative splicing (AS) and linked chromatin regulation. In the absence of FLAIL, splicing defects at the direct FLAIL target flowering gene LACCASE 8 (LAC8) correlated with reduced mRNA expression. Double mutant analyses support a model where FLAIL-mediated splicing of LAC8 promotes its mRNA expression and represses flowering. Our study suggests lncRNAs as accessory components of the spliceosome that regulate AS and gene expression to impact organismal development.

In Vivo Efficacy and Safety Evaluations of Therapeutic Splicing Correction Using U1 snRNA in the Mouse Retina

Efficacy and safety considerations constitute essential steps during development of in vivo gene therapies. Herein, we evaluated efficacy and safety of splice factor-based treatments to correct mutation-induced splice defects in an Opa1 mutant mouse line. We applied adeno-associated viruses to the retina. The viruses transduced retinal cells with an engineered U1 snRNA splice factor designed to correct the Opa1 splice defect. We found the treatment to be efficient in increasing wild-type Opa1 transcripts. Correspondingly, Opa1 protein levels increased significantly in treated eyes. Measurements of retinal morphology and function did not reveal therapy-related side-effects supporting the short-term safety of the treatment. Alterations of potential off-target genes were not detected. Our data suggest that treatments of splice defects applying engineered U1 snRNAs represent a promising in vivo therapeutic approach. The therapy increased wild-type Opa1 transcripts and protein levels without detectable morphological, functional or genetic side-effects in the mouse eye. The U1 snRNA-based therapy can be tailored to specific disease gene mutations, hence, raising the possibility of a wider applicability of this promising technology towards treatment of different inherited retinal diseases.

Aberrant splicing caused by exonic single nucleotide variants positioned 2nd or 3rd to the last nucleotide in the COL4A5 gene

BACKGROUND AND OBJECTIVES

The evident genotype-phenotype correlation shown by the X-linked Alport syndrome warrants the assessment of the impact of identified gene variants on aberrant splicing. We previously reported that single nucleotide variants (SNVs) in the last nucleotide of exons in COL4A5 cause aberrant splicing. It is known that the nucleotides located 2nd and 3rd to the last nucleotides of exons can also play an essential role in the first step of the splicing process. In this study, we aimed to investigate whether SNVs positioned 2nd or 3rd to the last nucleotide of exons in COL4A5 resulted in aberrant splicing.

METHODS

We selected eight candidate variants: six from the Human Gene Variant Database Professional and two from our cohort. We performed an in-vitro splicing assay and reverse transcription-polymerase chain reaction (RT-PCR) for messenger RNA obtained from patients, if available.

RESULTS

The candidate variants were initially classified into the following groups: three nonsense, two missense, and three synonymous variants. Splicing assays and RT-PCR for messenger RNA revealed that six of the eight variants caused aberrant splicing. Four variants, initially classified as non-truncating variants, were found to be truncating ones, which usually show relatively more severe phenotypes.

CONCLUSION

We revealed that exonic SNVs positioned 2nd or 3rd to the last nucleotide of exons in the COL4A5 were responsible for aberrant splicing. The results of our study suggest that attention should be paid when interpreting the pathogenicity of exonic SNVs near the 5' splice site.

Epigenetic Dysregulations in Arsenic-Induced Carcinogenesis

Arsenic is a crucial environmental metalloid whose high toxicity levels negatively impact human health. It poses significant health concerns to millions of people in developed and developing countries such as the USA, Canada, Bangladesh, India, China, and Mexico by enhancing sensitivity to various types of diseases, including cancers. However, how arsenic causes changes in gene expression that results in heinous conditions remains elusive. One of the proposed essential mechanisms that still has seen limited research with regard to causing disease upon arsenic exposure is the dysregulation of epigenetic components. In this review, we have extensively summarized current discoveries in arsenic-induced epigenetic modifications in carcinogenesis and angiogenesis. Importantly, we highlight the possible mechanisms underlying epigenetic reprogramming through arsenic exposure that cause changes in cell signaling and dysfunctions of different epigenetic elements.

The function and regulatory mechanism of RNA-binding proteins in breast cancer and their future clinical treatment prospects

Breast cancer is the most common female malignancy, but the mechanisms regulating gene expression leading to its development are complex. In recent years, as epigenetic research has intensified, RNA-binding proteins (RBPs) have been identified as a class of posttranscriptional regulators that can participate in regulating gene expression through the regulation of RNA stabilization and degradation, intracellular localization, alternative splicing and alternative polyadenylation, and translational control. RBPs play an important role in the development of normal mammary glands and breast cancer. Functional inactivation or abnormal expression of RBPs may be closely associated with breast cancer development. In this review, we focus on the function and regulatory mechanisms of RBPs in breast cancer, as well as the advantages and challenges of RBPs as potential diagnostic and therapeutic targets in breast cancer, and discuss the potential of RBPs in clinical treatment.

The Role of Alternative Splicing in Cancer: Regulatory Mechanism, Therapeutic Strategy, and Bioinformatics Application

[Formula: see text] Alternative splicing (AS) can generate distinct transcripts and subsequent isoforms that play differential functions from the same pre-mRNA. Recently, increasing numbers of studies have emerged, unmasking the association between AS and cancer. In this review, we arranged AS events that are closely related to cancer progression and presented promising treatments based on AS for cancer therapy. Obtaining proliferative capacity, acquiring invasive properties, gaining angiogenic features, shifting metabolic ability, and getting immune escape inclination are all splicing events involved in biological processes. Spliceosome-targeted and antisense oligonucleotide technologies are two novel strategies that are hopeful in tumor therapy. In addition, bioinformatics applications based on AS were summarized for better prediction and elucidation of regulatory routines mingled in. Together, we aimed to provide a better understanding of complicated AS events associated with cancer biology and reveal AS a promising target of cancer treatment in the future.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

The Function of DNA Demethylase Gene ROS1a Null Mutant on Seed Development in Rice (Oryza Sativa) Using the CRISPR/CAS9 System

The endosperm is the main nutrient source in cereals for humans, as it is a highly specialized storage organ for starch, lipids, and proteins, and plays an essential role in seed growth and development. Active DNA demethylation regulates plant developmental processes and is ensured by cytosine methylation (5-meC) DNA glycosylase enzymes. To find out the role of OsROS1a in seed development, the null mutant of OsROS1a was generated using the CRISPR/Cas9 system. The null mutant of OsROS1a was stable and heritable, which affects the major agronomic traits, particularly in rice seeds. The null mutant of OsROS1a showed longer and narrower grains, and seeds were deformed containing an underdeveloped and less-starch-producing endosperm with slightly irregularly shaped embryos. In contrast to the transparent grains of the wild type, the grains of the null mutant of OsROS1a were slightly opaque and rounded starch granules, with uneven shapes, sizes, and surfaces. A total of 723 differential expression genes (DEGs) were detected in the null mutant of OsROS1a by RNA-Seq, of which 290 were downregulated and 433 were upregulated. The gene ontology (GO) terms with the top 20 enrichment factors were visualized for cellular components, biological processes, and molecular functions. The key genes that are enriched for these GO terms include starch synthesis genes (OsSSIIa and OsSSIIIa) and cellulose synthesis genes (CESA2, CESA3, CESA6, and CESA8). Genes encoding polysaccharides and glutelin were found to be downregulated in the mutant endosperm. The glutelins were further verified by SDS-PAGE, suggesting that glutelin genes could be involved in the null mutant of OsROS1a seed phenotype and OsROS1a could have the key role in the regulation of glutelins. Furthermore, 378 differentially alternative splicing (AS) genes were identified in the null mutant of OsROS1a, suggesting that the OsROS1a gene has an impact on AS events. Our findings indicated that the function on rice endosperm development in the null mutant of OsROS1a could be influenced through regulating gene expression and AS, which could provide the base to properly understand the molecular mechanism related to the OsROS1a gene in the regulation of rice seed development.

Alternative splicing modulation by G-quadruplexes

Alternative splicing is central to metazoan gene regulation, but the regulatory mechanisms are incompletely understood. Here, we show that G-quadruplex (G4) motifs are enriched ~3-fold near splice junctions. The importance of G4s in RNA is emphasised by a higher enrichment for the non-template strand. RNA-seq data from mouse and human neurons reveals an enrichment of G4s at exons that were skipped following depolarisation induced by potassium chloride. We validate the formation of stable RNA G4s for three candidate splice sites by circular dichroism spectroscopy, UV-melting and fluorescence measurements. Moreover, we find that sQTLs are enriched at G4s, and a minigene experiment provides further support for their role in promoting exon inclusion. Analysis of >1,800 high-throughput experiments reveals multiple RNA binding proteins associated with G4s. Finally, exploration of G4 motifs across eleven species shows strong enrichment at splice sites in mammals and birds, suggesting an evolutionary conserved splice regulatory mechanism.

Here the authors shows that G-quadruplexes, non-canonical DNA/RNA structures, can have a direct impact on alternative splicing and that binding of splicing regulators is affected by their presence.

Expression complementation of gene presence/absence polymorphisms in hybrids contributes importantly to heterosis in sunflower

INTRODUCTION

Numerous crops have transitioned to hybrid seed production to increase yields and yield stability through heterosis. However, the molecular mechanisms underlying heterosis and its stability across environments are not yet fully understood.

OBJECTIVES

This study aimed to (1) elucidate the genetic and molecular mechanisms underlying heterosis in sunflower, and (2) determine how heterosis is maintained under different environments.

METHODS

Genome-wide association (GWA) analyses were employed to assess the effects of presence/absence variants (PAVs) and stop codons on 16 traits phenotyped in the sunflower association mapping population at three locations. To link the GWA results to transcriptomic variation, we sequenced the transcriptomes of two sunflower cultivars and their F₁ hybrid (INEDI) under both control and drought conditions and analyzed patterns of gene expression and alternative splicing.

RESULTS

Thousands of PAVs were found to affect phenotypic variation using a relaxed significance threshold, and at most such loci the "absence" allele reduced values of heterotic traits, but not those of non-heterotic traits. This pattern was strengthened for PAVs that showed expression complementation in INEDI. Stop codons were much rarer than PAVs and less likely to reduce heterotic trait values. Hybrid expression patterns were enriched for the GO category, sensitivity to stimulus, but all genotypes responded to drought similarily - by up-regulating water stress response pathways and down-regulating metabolic pathways. Changes in alternative splicing were strongly negatively correlated with expression variation, implying that alternative splicing in this system largely acts to reinforce expression responses.

CONCLUSION

Our results imply that complementation of expression of PAVs in hybrids is a major contributor to heterosis in sunflower, consistent with the dominance model of heterosis. This mechanism can account for yield stability across different environments. Moreover, given the much larger numbers of PAVs in plant vs. animal genomes, it also offers an explanation for the stronger heterotic responses seen in the former.

Gene architecture directs splicing outcome in separate nuclear spatial regions

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

In Vivo and In Vitro Matured Oocytes From Mice of Advanced Reproductive Age Exhibit Alternative Splicing Processes for Mitochondrial Oxidative Phosphorylation

The mean age of women seeking infertility treatment has gradually increased over recent years. This has coincided with the emergence of in vitro maturation (IVM), a method used in assisted reproductive technology for patients with special requirements. However, when compared with conventional in vitro fertilization, IVM is associated with poor embryonic development potential and low live birth rates, thus limiting the widespread application of this technique. In this study, we performed RNA-sequencing transcriptomic assays and identified a total of 2,627 significant differentially expressed genes (DEGs) between IVM oocytes and in vivo matured oocytes from mice of advanced reproductive age. Next, Kyoto Encyclopedia of Genes and Genomes pathway analysis was used to identify the potential functions of the DEGs. The most significantly enriched pathway was oxidative phosphorylation (OXPHOS). In addition, we constructed a protein-protein interaction network to identify key genes and determined that most of the hub genes were mtDNA-encoded subunits of respiratory chain complex I. Antioxidant supplementation lead to an increase in ATP production and reduced the gene expression profile of the OXPHOS pathway in the IVM group. Moreover, alternative splicing (AS) events were identified during in vivo or in vitro oocyte maturation; data showed that skipped exons were the most frequent type of AS event. A number of genes associated with the OXPHOS pathway exhibited alterations in AS events, including Ndufa7, Ndufs7, Cox6a2, Ndufs5, Ndufb1, and Uqcrh. Furthermore, the process of IVO promoted the skipping of exon 2 in Ndufa7 and exon 3 in Ndufs7 compared with the IVM oocytes, as determined by semi-quantitative RT-PCR. Collectively, these findings provide potential new therapeutic targets for improving IVM of aged women who undergo infertility treatment.

Transcriptomic Analysis Reveals That Rho GTPases Regulate Trap Development and Lifestyle Transition of the Nematode-Trapping Fungus Arthrobotrys oligospora

Nematode-trapping (NT) fungi can form unique infection structures (traps) to capture and kill free-living nematodes and, thus, can play a potential role in the biocontrol of nematodes. Arthrobotrys oligospora is a representative species of NT fungi. Here, we performed a time course transcriptome sequencing (RNA-seq) analysis of transcriptomes to understand the global gene expression levels of A. oligospora during trap formation and predation. We identified 5,752 unique differentially expressed genes, among which the rac gene was significantly upregulated. Alternative splicing events occurred in 2,012 genes, including the rac and rho2 gene. Furthermore, we characterized three Rho GTPases (Rho2, Rac, and Cdc42) in A. oligospora using gene disruption and multiphenotypic analysis. Our analyses showed that AoRac and AoCdc42 play an important role in mycelium growth, lipid accumulation, DNA damage, sporulation, trap formation, pathogenicity, and stress response in A. oligospora. AoCdc42 and AoRac specifically interacted with components of the Nox complex, thus regulating the production of reactive oxygen species. Moreover, the transcript levels of several genes associated with protein kinase A, mitogen-activated protein kinase, and p21-activated kinase were also altered in the mutants, suggesting that Rho GTPases might function upstream from these kinases. This study highlights the important role of Rho GTPases in A. oligospora and provides insights into the regulatory mechanisms of signaling pathways in the trap morphogenesis and lifestyle transition of NT fungi. IMPORTANCE Nematode-trapping (NT) fungi are widely distributed in terrestrial and aquatic ecosystems. Their broad adaptability and flexible lifestyles make them ideal agents for controlling pathogenic nematodes. Arthrobotrys oligospora is a model species employed for understanding the interaction between fungi and nematodes. Here, we revealed that alternative splicing events play a crucial role in the trap development and lifestyle transition in A. oligospora. Furthermore, Rho GTPases exert differential effects on the growth, development, and pathogenicity of A. oligospora. In particular, AoRac is required for sporulation and trap morphogenesis. In addition, our analysis showed that Rho GTPases regulate the production of reactive oxygen species and function upstream from several kinases. Collectively, these results expand our understanding of gene expression and alternative splicing events in A. oligospora and the important roles of Rho GTPases in NT fungi, thereby providing a foundation for exploring their potential application in the biocontrol of pathogenic nematodes.

Temporal regulation of alternative splicing events in rice memory under drought stress

Plant adaptation to drought stress is essential for plant survival and crop yield. Recently, harnessing drought memory, which is induced by repeated stress and recovery cycles, was suggested as a means to improve drought resistance at the transcriptional level. However, the genetic mechanism underlying drought memory is unclear. Here, we carried out a quantitative analysis of alternative splicing (AS) events in rice memory under drought stress, generating 12 transcriptome datasets. Notably, we identified exon skipping (ES) as the predominant AS type (>80%) in differential alternative splicing (DAS) in response to drought stress. Applying our analysis pipeline to investigate DAS events following drought stress in six other plant species revealed variable ES frequencies ranging from 9.94% to 60.70% depending on the species, suggesting that the relative frequency of DAS types in plants is likely to be species-specific. The dinucleotide sequence at AS splice sites in rice following drought stress was preferentially GC-AG and AT-AC. Since U12-type splicing uses the AT-AC site, this suggests that drought stress may increase U12-type splicing, and thus increase ES frequency. We hypothesize that multiple isoforms derived from exon skipping may be induced by drought stress in rice. We also identified 20 transcription factors and three highly connected hub genes with potential roles in drought memory that may be good targets for plant breeding.

Genome-wide differences in gene expression and alternative splicing in developing embryo and endosperm, and between F1 hybrids and their parental pure lines in sorghum

Developing embryo and endosperm of sorghum show substantial and multifaceted differences in gene expression and alternative splicing, which are potentially relevant to heterosis. Differential regulation of gene expression and alternative splicing (AS) are major molecular mechanisms dictating plant growth and development, as well as underpinning heterosis in F1 hybrids. Here, using deep RNA-sequencing we analyzed differences in genome-wide gene expression and AS between developing embryo and endosperm, and between F1 hybrids and their pure-line parents in sorghum. We uncover dramatic differences in both gene expression and AS between embryo and endosperm with respect to gene features and functions, which are consistent with the fundamentally different biological roles of the two tissues. Accordingly, F1 hybrids showed substantial and multifaceted differences in gene expression and AS compared with their pure-line parents, again with clear tissue specificities including extents of difference, genes involved and functional enrichments. Our results provide useful transcriptome resources as well as novel insights for further elucidation of seed yield heterosis in sorghum and related crops.

Four alternative splicing transcripts of intracellular copper/zinc superoxide dismutase 1 in Oxya chinensis

In this study, we obtained four alternative splicing transcripts of intracellular copper/zinc superoxide dismutase 1 (icCuZnSOD1) in Oxya chinensis. OcicCuZnSOD1a has all common characteristics of CuZnSOD family and is a canonical CuZnSOD. OcicCuZnSOD1b is missing a Zn binding site. OcicCuZnSOD1c lacks Zn ion and is a Cu-only SOD. OcicCuZnSOD1d is missing a CuZnSOD conserved sequence and lacks the E-loop, a conserved disulfide bond, and an active site arginine. OcicCuZnSOD1a was the most heat-resistant and OcicCuZnSOD1c was the most unstable at high temperatures above 55 °C. They were stable at a wide pH range, especially in alkaline conditions. The four variants expressed at the throughout developmental stages and had various tissue expression patterns. OcicCuZnSOD1a and OcicCuZnSOD1d were significantly induced by 8.79 mM CuCl₂ and OcicCuZnSOD1b was significantly up-regulated by 14.67 mM CuCl₂. OcicCuZnSOD1a was significantly inhibited by 19.13 mM ZnSO₄ while OcicCuZnSOD1d were significantly induced by 22.61 mM ZnSO₄. Disc diffusion assay showed that the four isoforms of OcicCuZnSOD1 made the killing zones smaller surrounding the CdCl₂-soaked filter discs. However, the reduction ratios of OcicCuZnSOD1a were the highest. These results implied that the four transcripts played roles in defense against CdCl₂-induced oxidative stress while OcicCuZnSOD1a had stronger antioxidant capacity.

Evolution of the Early Spliceosomal Complex-From Constitutive to Regulated Splicing

Pre-mRNA splicing is a major process in the regulated expression of genes in eukaryotes, and alternative splicing is used to generate different proteins from the same coding gene. Splicing is a catalytic process that removes introns and ligates exons to create the RNA sequence that codifies the final protein. While this is achieved in an autocatalytic process in ancestral group II introns in prokaryotes, the spliceosome has evolved during eukaryogenesis to assist in this process and to finally provide the opportunity for intron-specific splicing. In the early stage of splicing, the RNA 5' and 3' splice sites must be brought within proximity to correctly assemble the active spliceosome and perform the excision and ligation reactions. The assembly of this first complex, termed E-complex, is currently the least understood process. We focused in this review on the formation of the E-complex and compared its composition and function in three different organisms. We highlight the common ancestral mechanisms in S. cerevisiae, S. pombe, and mammals and conclude with a unifying model for intron definition in constitutive and regulated co-transcriptional splicing.

Arsenic-induced epigenetic changes in cancer development

Arsenic is a ubiquitous metalloid whose high levels of toxicity pose major health concerns to millions of people worldwide by increasing susceptibility to various cancers and non-cancer illnesses. Since arsenic is not a mutagen, the mechanism by which it causes changes in gene expression and disease pathogenesis is not clear. One possible mechanism is through generation of reactive oxygen species. Another equally important mechanism still very much in its infancy is epigenetic dysregulation. In this review, we discuss recent discoveries underlying arsenic-induced epigenetic changes in cancer development. Importantly, we highlight the proposed mechanisms targeted by arsenic to drive oncogenic gene expression.

Cross talk between the upstream exon-intron junction and Prp2 facilitates splicing of non-consensus introns

Splicing of mRNA precursors is essential in the regulation of gene expression. U2AF65 recognizes the poly-pyrimidine tract and helps in the recognition of the branch point. Inactivation of fission yeast U2AF65 (Prp2) blocks splicing of most, but not all, pre-mRNAs, for reasons that are not understood. Here, we have determined genome-wide the splicing efficiency of fission yeast cells as they progress into synchronous meiosis in the presence or absence of functional Prp2. Our data indicate that in addition to the splicing elements at the 3' end of any intron, the nucleotides immediately upstream the intron will determine whether Prp2 is required or dispensable for splicing. By changing those nucleotides in any given intron, we regulate its Prp2 dependency. Our results suggest a model in which Prp2 is required for the coordinated recognition of both intronic ends, placing Prp2 as a key regulatory element in the determination of the exon-intron boundaries.

The intronic branch point sequence is under strong evolutionary constraint in the bovine and human genome

The branch point sequence is a cis-acting intronic motif required for mRNA splicing. Despite their functional importance, branch point sequences are not routinely annotated. Here we predict branch point sequences in 179,476 bovine introns and investigate their variability using a catalogue of 29.4 million variants detected in 266 cattle genomes. We localize the bovine branch point within a degenerate heptamer "nnyTrAy". An adenine residue at position 6, that acts as branch point, and a thymine residue at position 4 of the heptamer are more strongly depleted for mutations than coding sequences suggesting extreme purifying selection. We provide evidence that mutations affecting these evolutionarily constrained residues lead to alternative splicing. We confirm evolutionary constraints on branch point sequences using a catalogue of 115 million SNPs established from 3,942 human genomes of the gnomAD database.

Aberrant RNA splicing due to genetic incompatibilities in sunflower hybrids

Genome-scale studies have revealed divergent mRNA splicing patterns between closely related species or populations. However, it is unclear whether splicing differentiation is a simple byproduct of population divergence, or whether it also acts as a mechanism for reproductive isolation. We examined mRNA splicing in wild × domesticated sunflower hybrids and observed 45 novel splice forms that were not found in the wild or domesticated parents, in addition to 16 high-expression parental splice forms that were absent in one or more hybrids. We identify loci associated with variation in the levels of these splice forms, finding that many aberrant transcripts were regulated by multiple alleles with nonadditive interactions. We identified particular spliceosome components that were associated with 21 aberrant isoforms, more than half of which were located in or near regulatory QTL. These incompatibilities often resulted in alteration in the protein-coding regions of the novel transcripts in the form of frameshifts and truncations. By associating the splice variation in these genes with size and growth rate measurements, we found that the cumulative expression of all aberrant transcripts was correlated with a significant reduction in growth rate. Our results lead us to propose a model where divergent splicing regulatory loci could act as incompatibility loci that contribute to the evolution of reproductive isolation.

Post-transcriptional and translational control of the morphology and virulence in human fungal pathogens

Host-pathogen interactions at the molecular level are the key to fungal pathogenesis. Fungal pathogens utilize several mechanisms such as adhesion, invasion, phenotype switching and metabolic adaptations, to survive in the host environment and respond. Post-transcriptional and translational regulations have emerged as key regulatory mechanisms ensuring the virulence and survival of fungal pathogens. Through these regulations, fungal pathogens effectively alter their protein pool, respond to various stress, and undergo morphogenesis, leading to efficient and comprehensive changes in fungal physiology. The regulation of virulence through post-transcriptional and translational regulatory mechanisms is mediated through mRNA elements (cis factors) or effector molecules (trans factors). The untranslated regions upstream and downstream of the mRNA, as well as various RNA-binding proteins involved in translation initiation or circularization of the mRNA, play pivotal roles in the regulation of morphology and virulence by influencing protein synthesis, protein isoforms, and mRNA stability. Therefore, post-transcriptional and translational mechanisms regulating the morphology, virulence and drug-resistance processes in fungal pathogens can be the target for new therapeutics. With improved "omics" technologies, these regulatory mechanisms are increasingly coming to the forefront of basic biology and drug discovery. This review aims to discuss various modes of post-transcriptional and translation regulations, and how these mechanisms exert influence in the virulence and morphogenesis of fungal pathogens.

An ADAR1-dependent RNA editing event in the cyclin-dependent kinase CDK13 promotes thyroid cancer hallmarks

Background

Adenosine deaminases acting on RNA (ADARs) modify many cellular RNAs by catalyzing the conversion of adenosine to inosine (A-to-I), and their deregulation is associated with several cancers. We recently showed that A-to-I editing is elevated in thyroid tumors and that ADAR1 is functionally important for thyroid cancer cell progression. The downstream effectors regulated or edited by ADAR1 and the significance of ADAR1 deregulation in thyroid cancer remain, however, poorly defined.

Methods

We performed whole transcriptome sequencing to determine the consequences of ADAR1 deregulation for global gene expression, RNA splicing and editing. The effects of gene silencing or RNA editing were investigated by analyzing cell viability, proliferation, invasion and subnuclear localization, and by protein and gene expression analysis.

Results

We report an oncogenic function for CDK13 in thyroid cancer and identify a new ADAR1-dependent RNA editing event that occurs in the coding region of its transcript. CDK13 was significantly over-edited (c.308A > G) in tumor samples and functional analysis revealed that this editing event promoted cancer cell hallmarks. Finally, we show that CDK13 editing increases the nucleolar abundance of the protein, and that this event might explain, at least partly, the global change in splicing produced by ADAR1 deregulation.

Conclusions

Overall, our data support A-to-I editing as an important pathway in cancer progression and highlight novel mechanisms that might be used therapeutically in thyroid and other cancers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-021-01401-y.

Spliceosomal component PRP-40 is a central regulator of microexon splicing

Microexons (≤27 nt) play critical roles in nervous system development and function but create unique challenges for the splicing machinery. The mechanisms of microexon regulation are therefore of great interest. We performed a genetic screen for alternative splicing regulators in the C. elegans nervous system and identify PRP-40, a core component of the U1 snRNP. RNA-seq reveals that PRP-40 is required for inclusion of alternatively spliced, but not constitutively spliced, exons. PRP-40 is particularly required for inclusion of neuronal microexons, and our data indicate that PRP-40 is a central regulator of microexon splicing. Microexons can be relieved from PRP-40 dependence by artificially increasing exon size or reducing flanking intron size, indicating that PRP-40 is specifically required for microexons surrounded by conventionally sized introns. Knockdown of the orthologous PRPF40A in mouse neuroblastoma cells causes widespread dysregulation of microexons but not conventionally sized exons. PRP-40 regulation of neuronal microexons is therefore a widely conserved phenomenon.

A SNP assay for assessing diversity in immune genes in the honey bee (Apis mellifera L.)

With a growing number of parasites and pathogens experiencing large-scale range expansions, monitoring diversity in immune genes of host populations has never been so important because it can inform on the adaptive potential to resist the invaders. Population surveys of immune genes are becoming common in many organisms, yet they are missing in the honey bee (Apis mellifera L.), a key managed pollinator species that has been severely affected by biological invasions. To fill the gap, here we identified single nucleotide polymorphisms (SNPs) in a wide range of honey bee immune genes and developed a medium-density assay targeting a subset of these genes. Using a discovery panel of 123 whole-genomes, representing seven A. mellifera subspecies and three evolutionary lineages, 180 immune genes were scanned for SNPs in exons, introns (< 4 bp from exons), 3' and 5´UTR, and < 1 kb upstream of the transcription start site. After application of multiple filtering criteria and validation, the final medium-density assay combines 91 quality-proved functional SNPs marking 89 innate immune genes and these can be readily typed using the high-sample-throughput iPLEX MassARRAY system. This medium-density-SNP assay was applied to 156 samples from four countries and the admixture analysis clustered the samples according to their lineage and subspecies, suggesting that honey bee ancestry can be delineated from functional variation. In addition to allowing analysis of immunogenetic variation, this newly-developed SNP assay can be used for inferring genetic structure and admixture in the honey bee.

SRSF5 regulates alternative splicing of DMTF1 pre-mRNA through modulating SF1 binding

ABBREVIATIONS

ARF: alternative reading frame, that is, p14ARF, or CDKN2A (cyclin-dependent kinase inhibitor 2A); β-gal: β-galactosidase; CLIP-seq: crosslinking and immunoprecipitation-sequencing; DMTF1: the cyclin D binding myb-like transcription factor 1; ESS/ESE: exonic splicing silencer/enhancer; Ex: exon; FBS: fetal bovine serum; Gluc: Gaussia luciferase; hnRNPs: heterogeneous nuclear ribonucleoproteins; In: intron; ISS/ISE: intronic splicing silencer/enhancer; PBS: phosphate-buffered saline; PCR: polymerase chain reaction; PSI: percent-splice-in; qPCR: quantitative real-time PCR; RIP: RNA immunoprecipitation; RNAseq: RNA sequencing; RT: reverse transcription; SF1: splicing factor 1; SR: serine/arginine-rich proteins; SRSF5: serine and arginine-rich splicing factor 5; TCGA: the cancer genome atlas; UCSC: University of California, Santa Cruz. WT: Wild type.

DoChaP: the domain change presenter

Alternative splicing results in multiple transcripts of the same gene, possibly encoding for different protein isoforms with different domains. Whereas it is possible to manually determine the effect of alternative splicing on the domain composition for a single event, the process requires the tedious integration of several data sources; it is error prone and not feasible for genome-wide characterization of domains affected by differential splicing. To fulfill the need for an automated solution, we developed the Domain Change Presenter (DoChaP, https://dochap.bgu.ac.il/), a web server for the visualization of exon-domain associations. DoChaP visualizes all transcripts of a given gene, the encoded proteins and their domains, and enables a comparison between the transcripts and between their protein products. The colors and organization make the structural effect of alternative splicing events on protein structures easily identified. To enable the study of the conservation of exons structure, alternative splicing, and the effect of alternative splicing on protein domains, DoChaP also provides a two-species comparison of exon-domain associations. DoChaP thus provides a unique and easy-to-use visualization of the exon-domain association and conservation, and will facilitate the study of the structural effects of alternative splicing in health and disease.

Dominant optic atrophy: Culprit mitochondria in the optic nerve

Dominant optic atrophy (DOA) is an inherited mitochondrial disease leading to specific degeneration of retinal ganglion cells (RGCs), thus compromising transmission of visual information from the retina to the brain. Usually, DOA starts during childhood and evolves to poor vision or legal blindness, affecting the central vision, whilst sparing the peripheral visual field. In 20% of cases, DOA presents as syndromic disorder, with secondary symptoms affecting neuronal and muscular functions. Twenty years ago, we demonstrated that heterozygous mutations in OPA1 are the most frequent molecular cause of DOA. Since then, variants in additional genes, whose functions in many instances converge with those of OPA1, have been identified by next generation sequencing. OPA1 encodes a dynamin-related GTPase imported into mitochondria and located to the inner membrane and intermembrane space. The many OPA1 isoforms, resulting from alternative splicing of three exons, form complex homopolymers that structure mitochondrial cristae, and contribute to fusion of the outer membrane, thus shaping the whole mitochondrial network. Moreover, OPA1 is required for oxidative phosphorylation, maintenance of mitochondrial genome, calcium homeostasis and regulation of apoptosis, thus making OPA1 the Swiss army-knife of mitochondria. Understanding DOA pathophysiology requires the understanding of RGC peculiarities with respect to OPA1 functions. Besides the tremendous energy requirements of RGCs to relay visual information from the eye to the brain, these neurons present unique features related to their differential environments in the retina, and to the anatomical transition occurring at the lamina cribrosa, which parallel major adaptations of mitochondrial physiology and shape, in the pre- and post-laminar segments of the optic nerve. Three DOA mouse models, with different Opa1 mutations, have been generated to study intrinsic mechanisms responsible for RGC degeneration, and these have further revealed secondary symptoms related to mitochondrial dysfunctions, mirroring the more severe syndromic phenotypes seen in a subgroup of patients. Metabolomics analyses of cells, mouse organs and patient plasma mutated for OPA1 revealed new unexpected pathophysiological mechanisms related to mitochondrial dysfunction, and biomarkers correlated quantitatively to the severity of the disease. Here, we review and synthesize these data, and propose different approaches for embracing possible therapies to fulfil the unmet clinical needs of this disease, and provide hope to affected DOA patients.

On the relation of gene essentiality to intron structure: a computational and deep learning approach

Essential genes have been studied by copy number variants and deletions, both associated with introns. The premise of our work is that introns of essential genes have distinct characteristic properties. We provide support for this by training a deep learning model and demonstrating that introns alone can be used to classify essentiality. The model, limited to first introns, performs at an increased level, implicating first introns in essentiality. We identify unique properties of introns of essential genes, finding that their structure protects against deletion and intron-loss events, especially centered on the first intron. We show that GC density is increased in the first introns of essential genes, allowing for increased enhancer activity, protection against deletions, and improved splice site recognition. We find that first introns of essential genes are of remarkably smaller size than their nonessential counterparts, and to protect against common 3' end deletion events, essential genes carry an increased number of (smaller) introns. To demonstrate the importance of the seven features we identified, we train a feature-based model using only these features and achieve high performance.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

Sense-oriented AluYRa1 elements provide a lineage-specific transcription environment for polyadenylation

Transposable elements cause alternative splicing (AS) in different ways, contributing to transcript diversification. Alternative polyadenylation (APA), one of the AS events, is related to the generation of mRNA isoforms in 70% of human genes. In this study, we tried to investigate AluYRa1s located at the terminal region of cynomolgus monkey genes, utilizing both computational analysis and molecular experimentation. We found that ten genes had AluYRa1 at their 3' end, and nine of these AluYRa1s were sense-oriented. Furthermore, in seven genes, AluYRa1s were expected to have a similar consensus sequence for polyadenylation cleavage. Additional computational analysis using the annotation files from the UCSC database showed that AluYRa1 was more involved in polyadenylation than in open reading frame exon splicing. To examine the extent of AluYRa1 involvement in polyadenylation, RNA-seq data from 30 normal cynomolgus monkeys were analyzed using TAPAS, a recently devised software that detects all the promising polyadenylation sites including APA sites. We observed that approximately 74% of possible polyadenylation sites in the analyzed genes were provided by sense-oriented AluYRa1. In conclusion, AluYRa1 is an Old-World monkey-specific TE, and its sense-oriented insertion at the 3'UTR region tends to provide a favorable environment for polyadenylation, diversifying gene transcripts.

Expression and clinical importance of a newly discovered alternative splice variant of the gene for acrosin binding protein found in human brain tumors

Background

Acrosin binding protein (ACRBP) is a member of the cancer-testis antigen (CTA) family. Normally, ACRBP mRNA is expressed only in seminiferous tubules, while abnormally it is expressed in various types of cancers in tumor tissues, such as brain tumor.

Objectives

To determine the expression and clinical impact of a newly discovered splice variant of ACRBP in brain tumor.

Methods

Total RNA was extracted and reverse transcribed from 92 tumor specimens and 3 cell lines. Primers were designed to determine the expression of the new splice variant in all the samples. Quantitative real-time PCR (qPCR) was conducted for samples positive in reverse transcriptase-PCR. Association of the expression of ACRBP with the clinicopathological features of the various brain tumors was assessed statistically.

Results

The primers identified a newly discovered splice variant of ACRBP named ACRBP-V5a. The proportions of samples of the various brain tumor types positive for the ACRBP-V5a splicing variant were as follows: astrocytoma 10/33 (30%), glioblastoma 10/30 (33%), medulloblastoma 14/29 (48%), all tumors 34/92 (37%). Although we did not find a significant difference in the proportions of samples of various types of brain tumor tissues positive for the new splice variant (P > 0.05), levels of expression of the ACRBP-V5a splice variant were significantly different for tumor grade (P = 0.01) and tumor type (P = 0.02).

Conclusions

A newly discovered splice variant, ACRBP-V5a, is present in brain tumor. The new splicing variant may have discriminative value and potential importance in molecular-targeted therapy for brain tumors.

Effects of Alternative Splicing Events on Acute Myeloid Leukemia

As suggested by an increasing amount of evidence, there is alternative splicing (AS) modification within malignancy, which is related to cancer occurrence and development. AS within acute myeloid leukemia (AML) has not yet been systematically analyzed yet. This study analyzed the transcriptomic profiling and corresponding clinical data from AML cases based on The Cancer Genome Atlas (TCGA). In addition, the prediction model, along with the splicing network, was used to analyze the prognosis for AML patients according to the seven different AS event types. Among the 34,984 AS events across the 8830 genes, 2896 AS events were detected among 1905 genes, showing marked correlation with the overall survival of patients. The risk scoring model based on all AS event types was the most efficient in identifying the prognosis for AML patients. Meanwhile, the area under the curve at 1-, 3-, 5-year were 0.852, 0.935, 0.955, respectively. At the same time, the splicing regulating network, which was constituted by 21 splicing factor genes as well as 32 AS events related to survival, was characterized. In conclusion, our predictive model constructed based on the AS events accurately predicts the survival for AML patients. In addition, the network between AS events and splicing factor is established, which may serve as a potential mechanism.

Differential Alternative Splicing Genes and Isoform Regulation Networks of Rapeseed (Brassica napus L.) Infected with Sclerotinia sclerotiorum

Alternative splicing (AS) is a post-transcriptional level of gene expression regulation that increases transcriptome and proteome diversity. How the AS landscape of rapeseed (Brassica napus L.) changes in response to the fungal pathogen Sclerotinia sclerotiorum is unknown. Here, we analyzed 18 RNA-seq libraries of mock-inoculated and S. sclerotiorum-inoculated susceptible and tolerant B. napus plants. We found that infection increased AS, with intron retention being the main AS event. To determine the key genes functioning in the AS response, we performed a differential AS (DAS) analysis. We identified 79 DAS genes, including those encoding splicing factors, defense response proteins, crucial transcription factors and enzymes. We generated coexpression networks based on the splicing isoforms, rather than the genes, to explore the genes’ diverse functions. Using this weighted gene coexpression network analysis alongside a gene ontology enrichment analysis, we identified 11 modules putatively involved in the pathogen defense response. Within these regulatory modules, six DAS genes (ascorbate peroxidase 1, ser/arg-rich protein 34a, unknown function 1138, nitrilase 2, v-atpase f, and amino acid transporter 1) were considered to encode key isoforms involved in the defense response. This study provides insight into the post-transcriptional response of B. napus to S. sclerotiorum infection.

An NCBP3-Domain Protein Mediates Meiotic Silencing by Unpaired DNA

In the filamentous fungus Neurospora crassa, genes unpaired during meiosis are silenced by a process known as meiotic silencing by unpaired DNA (MSUD). MSUD utilizes common RNA interference (RNAi) proteins, such as Dicer and Argonaute, to target homologous mRNAs for silencing. Previously, we demonstrated that nuclear cap-binding proteins NCBP1 and NCBP2 are involved in MSUD. We report here that SAD-8, a protein similar to human NCBP3, also mediates silencing. Although SAD-8 is not essential for either vegetative or sexual development, it is required for MSUD. SAD-8 localizes predominantly in the nucleus and interacts with both NCBP1 and NCBP2. Similar to NCBP1 and NCBP2, SAD-8 interacts with a component (Argonaute) of the perinuclear meiotic silencing complex (MSC), further implicating the involvement of cap-binding proteins in silencing.