Function of alternative splicing

RNA isoform expression landscape of the human dorsal root ganglion generated from long-read sequencing

ABSTRACT

Splicing is a posttranscriptional RNA processing mechanism that enhances genomic complexity by creating multiple isoforms from the same gene. We aimed to characterize the isoforms expressed in the human peripheral nervous system, with the goal of creating a resource to identify novel isoforms of functionally relevant genes associated with somatosensation and nociception. We used long-read sequencing to document isoform expression in the human dorsal root ganglia from 3 organ donors and validated in silico by confirming expression in short-read sequencing from 3 independent organ donors. Nineteen thousand five hundred forty-seven isoforms of protein-coding genes were detected and validated. We identified 763 isoforms with at least one previously undescribed splice junction. Previously unannotated isoforms of multiple pain-associated genes, including ASIC3 , MRGPRX1 , and HNRNPK , were identified. In the novel isoforms of ASIC3 , a region comprising approximately 35% of the 5'UTR was excised. By contrast, a novel splice junction was used in isoforms of MRGPRX1 to include an additional exon upstream of the start codon, consequently adding a region to the 5'UTR. Novel isoforms of HNRNPK were identified, which used previously unannotated splice sites to both excise exon 14 and include a sequence in the 3' end of exon 13. This novel insertion is predicted to introduce a tyrosine phosphorylation site potentially phosphorylated by SRC. We also independently confirm a recently reported DRG-specific splicing event in WNK1 that gives insight into how painless peripheral neuropathy occurs when this gene is mutated. Our findings give a clear overview of mRNA isoform diversity in the human dorsal root ganglia obtained using long-read sequencing.

On the Roles of Protein Intrinsic Disorder in the Origin of Life and Evolution

Obviously, the discussion of different factors that could have contributed to the origin of life and evolution is clear speculation, since there is no way of checking the validity of most of the related hypotheses in practice, as the corresponding events not only already happened, but took place in a very distant past. However, there are a few undisputable facts that are present at the moment, such as the existence of a wide variety of living forms and the abundant presence of intrinsically disordered proteins (IDPs) or hybrid proteins containing ordered domains and intrinsically disordered regions (IDRs) in all living forms. Since it seems that the currently existing living forms originated from a common ancestor, their variety is a result of evolution. Therefore, one could ask a logical question of what role(s) the structureless and highly dynamic but vastly abundant and multifunctional IDPs/IDRs might have in evolution. This study represents an attempt to consider various ideas pertaining to the potential roles of protein intrinsic disorder in the origin of life and evolution.

Splicing the Difference: Harnessing the complexity of the transcriptome in hematopoiesis

Alternative splicing has long been recognized as a powerful tool to expand the diversity of the transcriptome and the proteome. The study of hematopoiesis, from hematopoietic stem cell maintenance and differentiation into committed progenitors to maturation into functional blood cells, has led the field of stem cell research and cellular differentiation for decades. The importance of aberrant splicing due to mutations in cis has been exemplified in thalassemias, resulting from aberrant expression of beta-globin. The simultaneous development of increasingly sophisticated technologies, in particular the combination of multi-color flowcytometric cell sorting with bulk and single cell sequencing, has provided sophisticated insights into the complex regulation of the blood system. The recognition that mutations in key splicing factors drive myeloid malignancies, in particular myelodysplastic syndromes, has galvanized research into alternative splicing in hematopoiesis and its diseases. In this review, we will update the audience on the exciting novel technologies, highlight alternative splicing events and their regulators with essential functions in hematopoiesis, and provide a high-level overview how splicing factor mutations contribute to hematologic malignancies.

Role of long noncoding RNAs in the regulation of alternative splicing in glioblastoma

Glioblastoma multiforme (GBM) is a highly severe primary brain tumor. Despite extensive research, effective treatments remain elusive. Long noncoding RNAs (lncRNAs) play a significant role in both cancer and normal biology. They influence alternative splicing (AS), which is crucial in cancer. Advances in lncRNA-specific microarrays and next-generation sequencing have enhanced understanding of AS. Abnormal AS contributes to cancer invasion, metastasis, apoptosis, therapeutic resistance, and tumor development, including glioma. lncRNA-mediated AS affects several cellular signaling pathways, promoting or suppressing cancer malignancy. This review discusses the lncRNAs regulating AS in glioblastoma and their mechanisms.

Importance of Transcript Variants in Transcriptome Analyses

RNA sequencing (RNA-Seq) has become a widely adopted technique for studying gene expression. However, conventional RNA-Seq analyses rely on gene expression (GE) values that aggregate all the transcripts produced under a single gene identifier, overlooking the complexity of transcript variants arising from different transcription start sites or alternative splicing. Transcript variants may encode proteins with diverse functional domains, or noncoding RNAs. This study explored the implications of neglecting transcript variants in RNA-Seq analyses. Among the 1334 transcription factor (TF) genes expressed in mouse embryonic stem (ES) or trophoblast stem (TS) cells, 652 were differentially expressed in TS cells based on GE values (365 upregulated and 287 downregulated, ≥absolute 2-fold changes, false discovery rate (FDR) p-value ≤ 0.05). The 365 upregulated genes expressed 883 transcript variants. Further transcript expression (TE) based analyses identified only 174 (<20%) of the 883 transcripts to be upregulated. The remaining 709 transcripts were either downregulated or showed no significant changes. Meanwhile, the 287 downregulated genes expressed 856 transcript variants and only 153 (<20%) of the 856 transcripts were downregulated. The other 703 transcripts were either upregulated or showed no significant change. Additionally, the 682 insignificant TF genes (GE values < absolute 2-fold changes and/or FDR p-values > 0.05) between ES and TS cells expressed 2215 transcript variants. These included 477 (>21%) differentially expressed transcripts (276 upregulated and 201 downregulated, ≥absolute 2-fold changes, FDR p-value ≤ 0.05). Hence, GE based RNA-Seq analyses do not represent accurate expression levels due to divergent transcripts expression from the same gene. Our findings show that by including transcript variants in RNA-Seq analyses, we can generate a precise understanding of a gene's functional and regulatory landscape; ignoring the variants may result in an erroneous interpretation.

Tumor antigen presentation and the associated signal transduction during carcinogenesis

Numerous developments have been achieved in the study and treatment of cancer throughout the decades that it has been common. After decades of research, about 100 different kinds of cancer have been found, each with unique subgroups within certain organs. This has significantly expanded our understanding of the illness. A mix of genetic, environmental, and behavioral variables contribute to the complicated and diverse process of cancer formation. Mutations, or changes in the DNA sequence, are crucial to the development of cancer. These mutations have the ability to downregulate the expression and function of Major Histocompatibility Complex class I (MHC I) and MHCII receptors, as well as activate oncogenes and inactivate tumor suppressor genes. Cancer cells use this tactic to avoid being recognized by cytotoxic CD8+T lymphocytes, which causes issues with antigen presentation and processing. This review goes into great length into the PI3K pathway, changes to MHC I, and positive impacts of tsMHC-II on disease-free survival and overall survival and the involvement of dendritic cells (DCs) in different tumor microenvironments. The vital functions that the PI3K pathway and its link to the mTOR pathway are highlighted and difficulties in developing effective cancer targeted therapies and feedback systems has also been mentioned, where resistance mechanisms include RAS-mediated oncogenic changes and active PI3K signalling.

Overcoming cancer drug-resistance calls for novel strategies targeting abnormal alternative splicing

Abnormal gene alternative splicing (AS) events are strongly associated with cancer progression. Here, we summarize AS events that contribute to the development of drug resistance and classify them into three categories: alternative cis-splicing (ACS), alternative trans-splicing (ATS), and alternative back-splicing (ABS). The regulatory mechanisms underlying AS processes through cis-acting regulatory elements and trans-acting factors are comprehensively described, and the distinct functions of spliced variants, including linear spliced variants derived from ACS, chimeric spliced variants arising from ATS, and circRNAs generated through ABS, are discussed. The identification of dysregulated spliced variants, which contribute to drug resistance and hinder effective cancer treatment, suggests that abnormal AS processes may together serve as a precise regulatory mechanism enabling drug-resistant cancer cell survival or, alternatively, represent an evolutionary pathway for cancer cells to adapt to changes in the external environment. Moreover, this review summarizes recent advancements in treatment approaches targeting AS-associated drug resistance, focusing on cis-acting regulatory elements, trans-acting factors, and specific spliced variants. Collectively, gaining an in-depth understanding of the mechanisms underlying aberrant alternative splicing events and developing strategies to target this process hold great promise for overcoming cancer drug resistance.

Comparative transcriptomics identifies genes underlying growth performance of the Pacific black-lipped pearl oyster Pinctada margaritifera

BACKGROUND

In bivalves, the rate at which organisms grow is a major functional trait underlying many aspects of their commercial production. Growth is a highly polygenic trait, which is typically regulated by many genes with small to moderate effects. Due to its complexity, growth variability in such shellfish remains poorly understood. In this study, we aimed to investigate differential gene expression among spat of the pearl oyster Pinctada margaritifera with distinct growth phenotypes.

RESULTS

We selected two groups of P. margaritifera spat belonging to the same F2 cohort based on their growth performance at 5.5 months old. Transcriptome profile analysis identified a total of 394 differentially expressed genes between these Fast-growing (F) and Slow-growing (S) phenotypes. According to functional enrichment analysis, S oysters overexpressed genes associated with stress-pathways and regulation of innate immune responses. In contrast, F oysters up-regulated genes associated with cytoskeleton activity, cell proliferation, and apoptosis. Analysis of genome polymorphism identified 16 single nucleotide polymorphisms (SNPs) significantly associated with the growth phenotypes. SNP effect categorization revealed one SNP identified for high effect and annotated for a stop codon gained mutation. Interestingly, this SNP is located within a gene annotated for scavenger receptor class F member 1 (SRF1), which is known to modulate apoptosis. Our analyses also revealed that all F oysters showed up-regulation for this gene and were homozygous for the stop-codon mutation. Conversely, S oysters had a heterozygous genotype and a reduced expression of this gene.

CONCLUSIONS

Altogether, our findings suggest that differences in growth among the same oyster cohort may be explained by contrasted metabolic allocation between regulatory pathways for growth and the immune system. This study provides a valuable contribution towards our understanding of the molecular components associated with growth performance in the pearl oyster P. margaritifera and bivalves in general.

AtSNU13 modulates pre-mRNA splicing of RBOHD and ALD1 to regulate plant immunity

Pre-mRNA splicing is a significant step for post-transcriptional modifications and functions in a wide range of physiological processes in plants. Human NHP2L binds to U4 snRNA during spliceosome assembly; it is involved in RNA splicing and mediates the development of human tumors. However, no ortholog has yet been identified in plants. Therefore, we report At4g12600 encoding the ortholog NHP2L protein, and AtSNU13 associates with the component of the spliceosome complex; the atsnu13 mutant showed compromised resistance in disease resistance, indicating that AtSNU13 is a positive regulator of plant immunity. Compared to wild-type plants, the atsnu13 mutation resulted in altered splicing patterns for defense-related genes and decreased expression of defense-related genes, such as RBOHD and ALD1. Further investigation shows that AtSNU13 promotes the interaction between U4/U6.U5 tri-snRNP-specific 27 K and the motif in target mRNAs to regulate the RNA splicing. Our study highlights the role of AtSNU13 in regulating plant immunity by affecting the pre-mRNA splicing of defense-related genes.

Alternative splicing plays a nonredundant role in greater amberjack (Seriola dumerili) in acclimation to ambient salinity fluctuations

Alternative splicing (AS) is an important post-transcriptional mechanism for adaptation of fish to environmental stress. Here, we performed a genome-wide investigation to AS dynamics in greater amberjack (Seriola dumerili), an economical marine teleost, in response to hypo- (10 ppt) and hyper-salinity (40 ppt) stresses. Totally, 2267-2611 differentially spliced events were identified in gills and kidney upon the exposure to undesired salinity regimes. In gills, genes involved in energy metabolism, stimulus response and epithelial cell differentiation were differentially spliced in response to salinity variation, while sodium ion transport and cellular amide metabolism were enhanced in kidney to combat the adverse impacts of salinity changes. Most of these differentially spliced genes were not differentially expressed, and AS was found to regulate different biological processes from differential gene expression, indicative of the functionally nonredundant role of AS in modulating salinity acclimation in greater amberjack. Together, our study highlights the important contribution of post-transcriptional mechanisms to the adaptation of fish to ambient salinity fluctuations and provides theoretical guidance for the conservation of marine fishery resources against increasingly environmental challenges.

The fallacy of functional nomenclature in the kingdom of biological multifunctionality: physiological and evolutionary considerations on ion channels

Living organisms are multiscale complex systems that have evolved high degrees of multifunctionality and redundancy in the structure-function relationship. A number of factors, only in part determined genetically, affect the jobs of proteins. The overall structural organization confers unique molecular properties that provide the potential to perform a pattern of activities, some of which are co-opted by specific environments. The variety of multifunctional proteins is expanding, but most cases are handled individually and according to the still dominant 'one structure-one function' approach, which relies on the attribution of canonical names typically referring to the first task identified for a given protein. The present topical review focuses on the multifunctionality of ion channels as a paradigmatic example. Mounting evidence reports the ability of many ion channels (including members of voltage-dependent, ligand-gated and transient receptor potential families) to exert biological effects independently of their ion conductivity. 'Functionally based' nomenclature (the practice of naming a protein or family of proteins based on a single purpose) is a conceptual bias for three main reasons: (i) it increases the amount of ambiguity, deceiving our understanding of the multiple contributions of biomolecules that is the heart of the complexity; (ii) it is in stark contrast to protein evolution dynamics, largely based on multidomain arrangement; and (iii) it overlooks the crucial role played by the microenvironment in adjusting the actions of cell structures and in tuning protein isoform diversity to accomplish adaptational requirements. Biological information in protein physiology is distributed among different entwined layers working as the primary 'locus' of natural selection and of evolutionary constraints.

CircRNA: An emerging star in plant research: A review

CircRNAs are a class of covalently closed non-coding RNA formed by linking the 5' terminus and the 3' terminus after reverse splicing. CircRNAs are widely found in eukaryotes, and they are highly conserved, with spatio-temporal expression specificity and stability. CircRNAs can act as miRNA sponges to regulate the expression of downstream target genes, regulating the transcription of parental genes and some can even be translated into peptides or proteins. Research on circRNAs in plants is still in its infancy compared to that in animals. With the deepening of research, the results of a variety of plant circRNAs suggest that they play an important role in growth and development, and tolerance towards abiotic stresses such as salt, drought, low temperature, high temperature and other adverse environments. In this review paper, we elaborated the molecular characteristics, mechanism of action, function and bioinformatics databases of plant circRNAs, combined with the progress of circRNA research in animals, discussed the potential mechanism of action of plant circRNAs, and proposed the unsolved problems and prospects for future application of plant circRNAs.

Pan-transcriptomic analysis reveals alternative splicing control of cold tolerance in rice

Plants have evolved complex mechanisms to adapt to harsh environmental conditions. Rice (Oryza sativa) is a staple food crop that is sensitive to low temperatures. However, its cold stress responses remain poorly understood, thus limiting possibilities for crop engineering to achieve greater cold tolerance. In this study, we constructed a rice pan-transcriptome and characterized its transcriptional regulatory landscape in response to cold stress. We performed Iso-Seq and RNA-Seq of 11 rice cultivars subjected to a time-course cold treatment. Our analyses revealed that alternative splicing-regulated gene expression plays a significant role in the cold stress response. Moreover, we identified CATALASE C (OsCATC) and Os03g0701200 as candidate genes for engineering enhanced cold tolerance. Importantly, we uncovered central roles for the 2 serine-arginine-rich proteins OsRS33 and OsRS2Z38 in cold tolerance. Our analysis of cold tolerance and resequencing data from a diverse collection of 165 rice cultivars suggested that OsRS2Z38 may be a key selection gene in japonica domestication for cold adaptation, associated with the adaptive evolution of rice. This study systematically investigated the distribution, dynamic changes, and regulatory mechanisms of alternative splicing in rice under cold stress. Overall, our work generates a rich resource with broad implications for understanding the genetic basis of cold response mechanisms in plants.

Molecular Mechanism of Resistance to Alternaria alternata Apple Pathotype in Apple by Alternative Splicing of Transcription Factor MdMYB6-like

As a fruit tree with great economic value, apple is widely cultivated in China. However, apple leaf spot disease causes significant damage to apple quality and economic value. In our study, we found that MdMYB6-like is a transcription factor without auto-activation activity and with three alternative spliced variants. Among them, MdMYB6-like-β responded positively to the pathogen infection. Overexpression of MdMYB6-like-β increased the lignin content of leaves and improved the pathogenic resistance of apple flesh callus. In addition, all three alternative spliced variants of MdMYB6-like could bind to the promoter of MdBGLU H. Therefore, we believe that MdMYB6-like plays an important role in the infection process of the pathogen and lays a solid foundation for breeding disease-resistant cultivars of apple in the future.

In Silico Genome-Wide Analysis of B3 Transcription Factors in Cannabis sativa L

Introduction: The B3 transcription factor has been identified in Arabidopsis thaliana, Oryza sativa, and Solanum lycopersicum, among other species. This family of transcription factors regulates seed growth, development, and stress. Cannabis is a valuable crop with numerous applications; however, no B3 transcription factors have been identified in this plant. Materials and Methods: The cannabis B3 gene family was identified and analyzed using bioinformatics analysis tools, such as the NCBI database, plantTFDB website, TBtools, and MEGA software. Quantitative real-time polymerase chain reaction (qRT-PCR) experiments were used to confirm its function. Results: The cannabis B3 family contains 65 members spread across 10 chromosomes. The isoelectric point ranged from 10.03 to 4.65, and the molecular weight ranged from 99,542.88 to 14,310.9 Da. Most of the members were found in the nucleus. The upstream promoter region of the gene contains a variety of cis-acting elements related to the stress response. RNA-seq data and qRT-PCR results showed that CsB3 genes were expressed differently in five organs of female Diku plants and in glandular hairs of nine distinct types of female cannabis inflorescences. Collinearity analysis revealed that there were more homologous genes between cannabis and dicotyledons than monocotyledonous plants, which was consistent with the evolutionary relationship. Conclusions: Hormones and external environmental factors might influence CsB3 expression. Furthermore, some genes such as CsB3-02, CsB3-07, CsB3-50, CsB3-62, and CsB3-65 may participate in cannabis growth and development and play a role in secondary metabolite synthesis. This study provides a solid foundation for further research into the gene function of the cannabis B3 family.

An array of signal-specific MoYpd1 isoforms determines full virulence in the pathogenic fungus Magnaporthe oryzae

Magnaporthe oryzae is placed first on a list of the world's top ten plant pathogens with the highest scientific and economic importance. The locus MGG_07173 occurs only once in the genome of M. oryzae and encodes the phosphotransfer protein MoYpd1p, which plays an important role in the high osmolarity glycerol (HOG) signaling pathway for osmoregulation. Originating from this locus, at least three MoYPD1 isoforms are produced in a signal-specific manner. The transcript levels of these MoYPD1-isoforms were individually affected by external stress. Salt (KCI) stress raised MoYPD1_T0 abundance, whereas osmotic stress by sorbitol elevates MoYPD1_T1 levels. In line with this, signal-specific nuclear translocation of green fluorescent protein-fused MoYpd1p isoforms in response to stress was observed. Mutant strains that produce only one of the MoYpd1p isoforms are less virulent, suggesting a combination thereof is required to invade the host successfully. In summary, we demonstrate signal-specific production of MoYpd1p isoforms that individually increase signal diversity and orchestrate virulence in M. oryzae.

CF2-II Alternative Splicing Isoform Regulates the Expression of Xenobiotic Tolerance-Related Cytochrome P450 CYP6CY22 in Aphis gossypii Glover

The expression of P450 genes is regulated by trans-regulatory factors or cis-regulatory elements and influences how endogenous or xenobiotic substances are metabolized in an organism's tissues. In this study, we showed that overexpression of the cytochrome P450 gene, CYP6CY22, led to resistance to cyantraniliprole in Aphis gossypii. The expression of CYP6CY22 increased in the midgut and remaining carcass of the CyR strain, and after repressing the expression of CYP6CY22, the mortality of cotton aphids increased 2.08-fold after exposure to cyantraniliprole. Drosophila ectopically expressing CYP6CY22 exhibited tolerance to cyantraniliprole and cross-tolerance to xanthotoxin, quercetin, 2-tridecanone, tannic acid, and nicotine. Moreover, transcription factor CF2-II (XM_027994540.2) is transcribed only as the splicing variant isoform CF2-II-AS, which was found to be 504 nucleotides shorter than CF2-II in A. gossypii. RNAi and yeast one-hybrid (Y1H) results indicated that CF2-II-AS positively regulates CYP6CY22 and binds to cis-acting element p (-851/-842) of CYP6CY22 to regulate its overexpression. The above results indicated that CYP6CY22 was regulated by the splicing isoform CF2-II-AS, which will help us further understand the mechanism of transcriptional adaption of cross-tolerance between synthetic insecticides and plant secondary metabolites mediated by P450s.

Systematic characterization of protein structural features of alternative splicing isoforms using AlphaFold 2

Alternative splicing is an important cellular process in eukaryotes, altering pre-mRNA to yield multiple protein isoforms from a single gene. However, our understanding of the impact of alternative splicing events on protein structures is currently constrained by a lack of sufficient protein structural data. To address this limitation, we employed AlphaFold 2, a cutting-edge protein structure prediction tool, to conduct a comprehensive analysis of alternative splicing for approximately 3,000 human genes, providing valuable insights into its impact on the protein structural. Our investigation employed state of the art high-performance computing infrastructure to systematically characterize structural features in alternatively spliced regions and identified changes in protein structure following alternative splicing events. Notably, we found that alternative splicing tends to alter the structure of residues primarily located in coils and beta-sheets. Our research highlighted a significant enrichment of loops and highly exposed residues within human alternatively spliced regions. Specifically, our examination of the Septin-9 protein revealed potential associations between loops and alternative splicing, providing insights into its evolutionary role. Furthermore, our analysis uncovered two missense mutations in the Tau protein that could influence alternative splicing, potentially contributing to the pathogenesis of Alzheimer's disease. In summary, our work, through a thorough statistical analysis of extensive protein structural data, sheds new light on the intricate relationship between alternative splicing, evolution, and human disease.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Predicting the Structural Impact of Human Alternative Splicing

Protein structure prediction with neural networks is a powerful new method for linking protein sequence, structure, and function, but structures have generally been predicted for only a single isoform of each gene, neglecting splice variants. To investigate the structural implications of alternative splicing, we used AlphaFold2 to predict the structures of more than 11,000 human isoforms. We employed multiple metrics to identify splicing-induced structural alterations, including template matching score, secondary structure composition, surface charge distribution, radius of gyration, accessibility of post-translational modification sites, and structure-based function prediction. We identified examples of how alternative splicing induced clear changes in each of these properties. Structural similarity between isoforms largely correlated with degree of sequence identity, but we identified a subset of isoforms with low structural similarity despite high sequence similarity. Exon skipping and alternative last exons tended to increase the surface charge and radius of gyration. Splicing also buried or exposed numerous post-translational modification sites, most notably among the isoforms of BAX. Functional prediction nominated numerous functional differences among isoforms of the same gene, with loss of function compared to the reference predominating. Finally, we used single-cell RNA-seq data from the Tabula Sapiens to determine the cell types in which each structure is expressed. Our work represents an important resource for studying the structure and function of splice isoforms across the cell types of the human body.

Polymorphism and expression level of the FADS3 gene and associated with the growth traits in Hu sheep

Growth traits are the economically important traits of sheep, and screening for genes related to growth and development is helpful for the genetic improvement of ovine growth traits. The fatty acid desaturase 3 (FADS3) is one of the important genes affecting the synthesis and accumulation of polyunsaturated fatty acids in animals. In this study, the expression levels of the FADS3 gene and polymorphism of the FADS3 gene associated with growth traits in Hu sheep were detected using quantitative real-time PCR (qRT-PCR), Sanger sequencing, and KAspar assay. The result showed that the expression levels of the FADS3 gene were widely expressed in all tissues, and the expression level of FADS3 in the lung was significantly higher than in other tissues (p < .05). Then, the polymorphism locus g. 2918 A > C was detected in intron 2 of the FADS3 gene, and associated analysis showed that the mutation in the FADS3 gene was associated significantly with growth traits (including body weight, body height, body length, and chest circumference, p < .05). Therefore, individuals with AA genotype showed significantly better growth traits than those with CC genotype, and FADS3 gene could be a candidate gene for improving growth traits in Hu sheep.

RNA isoform expression landscape of the human dorsal root ganglion (DRG) generated from long read sequencing

Splicing is a post-transcriptional RNA processing mechanism that enhances genomic complexity by creating multiple isoforms from the same gene. Diversity in splicing in the mammalian nervous system is associated with neuronal development, synaptic function and plasticity, and is also associated with diseases of the nervous system ranging from neurodegeneration to chronic pain. We aimed to characterize the isoforms expressed in the human peripheral nervous system, with the goal of creating a resource to identify novel isoforms of functionally relevant genes associated with somatosensation and nociception. We used long read sequencing (LRS) to document isoform expression in the human dorsal root ganglia (hDRG) from 3 organ donors. Isoforms were validated in silico by confirming expression in hDRG short read sequencing (SRS) data from 3 independent organ donors. 19,547 isoforms of protein-coding genes were detected using LRS and validated with SRS and strict expression cutoffs. We identified 763 isoforms with at least one previously undescribed splice-junction. Previously unannotated isoforms of multiple pain-associated genes, including ASIC3, MRGPRX1 and HNRNPK were identified. In the novel isoforms of ASIC3, a region comprising ~35% of the 5'UTR was excised. In contrast, a novel splice-junction was utilized in isoforms of MRGPRX1 to include an additional exon upstream of the start-codon, consequently adding a region to the 5'UTR. Novel isoforms of HNRNPK were identified which utilized previously unannotated splice-sites to both excise exon 14 and include a sequence in the 5' end of exon 13. The insertion and deletion in the coding region was predicted to excise a serine-phosphorylation site favored by cdc2, and replace it with a tyrosine-phosphorylation site potentially phosphorylated by SRC. We also independently confirm a recently reported DRG-specific splicing event in WNK1 that gives insight into how painless peripheral neuropathy occurs when this gene is mutated. Our findings give a clear overview of mRNA isoform diversity in the hDRG obtained using LRS. Using this work as a foundation, an important next step will be to use LRS on hDRG tissues recovered from people with a history of chronic pain. This should enable identification of new drug targets and a better understanding of chronic pain that may involve aberrant splicing events.

NRF2 has a splicing regulatory function involving the survival of motor neuron (SMN) in non-small cell lung cancer

The nuclear factor erythroid 2-like 2 (NFE2L2; NRF2) signaling pathway is frequently deregulated in human cancers. The critical functions of NRF2, other than its transcriptional activation, in cancers remain largely unknown. Here, we uncovered a previously unrecognized role of NRF2 in the regulation of RNA splicing. Global splicing analysis revealed that NRF2 knockdown in non-small cell lung cancer (NSCLC) A549 cells altered 839 alternative splicing (AS) events in 485 genes. Mechanistic studies demonstrated that NRF2 transcriptionally regulated SMN mRNA expression by binding to two antioxidant response elements in the SMN1 promoter. Post-transcriptionally, NRF2 was physically associated with the SMN protein. The Neh2 domain of NRF2, as well as the YG box and the region encoded by exon 7 of SMN, were required for their interaction. NRF2 formed a complex with SMN and Gemin2 in nuclear gems and Cajal bodies. Furthermore, the NRF2-SMN interaction regulated RNA splicing by expressing SMN in NRF2-knockout HeLa cells, reverting some of the altered RNA splicing. Moreover, SMN overexpression was significantly associated with alterations in the NRF2 pathway in patients with lung squamous cell carcinoma from The Cancer Genome Atlas. Taken together, our findings suggest a novel therapeutic strategy for cancers involving an aberrant NRF2 pathway.

A Splicing Transcriptome-Wide Association Study Identifies Candidate Altered Splicing for Prostate Cancer Risk

Prostate cancer (PCa) represents a huge public health burden among men. Many susceptibility genetic factors for PCa still remain unknown. In this study, we performed a large splicing transcriptome-wide association study (spTWAS) using three modeling strategies to develop alternative splicing genetic prediction models for identifying novel susceptibility loci and splicing introns for PCa risk by assessing 79,194 cases and 61,112 controls of European ancestry in the PRACTICAL, CRUK, CAPS, BPC3, and PEGASUS consortia. We identified 120 splicing introns of 97 genes showing an association with PCa risk at false discovery rate (FDR)-corrected threshold (FDR <0.05). Of them, 33 genes were enriched in PCa-related diseases and function categories. Fine-mapping analysis suggested that 21 splicing introns of 19 genes were likely causally associated with PCa risk. Thirty-five splicing introns of 34 novel genes were identified to be related to PCa susceptibility for the first time, and 11 of the genes were enriched in a cancer-related network. Our study identified novel loci and splicing introns associated with PCa risk, which can improve our understanding of the etiology of this common malignancy.

TMT Labeling Reveals the Effects of Exercises on the Proteomic Characteristics of the Subcutaneous Adipose Tissue of Growing High-Fat-Diet-Fed Rats

Aim: Growing period is an important period for fat remodeling. High-fat diet and exercise are reasons for adipose tissue (AT) remodeling, but existing evidence is not enough. Therefore, the effects of moderate-intensity continuous training (MICT) and high-intensity interval training (HIIT) on the proteomic characteristics of the subcutaneous AT of growing rats on normal diet or high-fat diet (HFD) were determined. Methods: Four-week-old male Sprague-Dawley rats (n = 48) were subdivided into six groups: normal diet control group, normal diet-MICT group, normal diet-HIIT group, HFD control group, HFD-MICT group, and HFD-HIIT group. Rats in the training group ran on a treadmill 5 days a week for 8 weeks (MICT: 50 min at 60-70% VO_2max intensity; HIIT: 7 min of warm-up and recovery at 70% VO_2max intensity, 6 sets of 3 min of 30% VO_2max followed by 3 min 90% VO_2max). Following physical assessment, inguinal subcutaneous adipose tissue (sWAT) was collected for proteome analysis using tandem mass tag labeling. Results: MICT and HIIT attenuated body fat mass and lean body mass but did not affect weight gain. Proteomics revealed the impact of exercise on ribosome, spliceosome, and the pentose phosphate pathway. However, the effect was reversed on HFD and normal diet. The differentially expressed proteins (DEPs) affected by MICT were related to oxygen transport, ribosome, and spliceosome. In comparison, the DEPs affected by HIIT were related to oxygen transport, mitochondrial electron transport, and mitochondrion protein. In HFD, HIIT was more likely to cause changes in immune proteins than MICT. However, exercise did not seem to reverse the protein effects of HFD. Conclusion: The exercise stress response in the growing period was stronger but increased the energy metabolism and metabolism. MICT and HIIT can reduce fat, increase muscle percentage, and improve maximum oxygen uptake in rats fed with HFD. However, in rats with normal diet, MICT and HIIT triggered more immune responses of sWAT, especially HIIT. In addition, spliceosomes may be the key factors in AT remodeling triggered by exercise and diet.

Splice variants of mitofusin 2 shape the endoplasmic reticulum and tether it to mitochondria

In eukaryotic cells, different organelles interact at membrane contact sites stabilized by tethers. Mitochondrial mitofusin 2 (MFN2) acts as a membrane tether that interacts with an unknown partner on the endoplasmic reticulum (ER). In this work, we identified the MFN2 splice variant ERMIT2 as the ER tethering partner of MFN2. Splicing of MFN2 produced ERMIT2 and ERMIN2, two ER-specific variants. ERMIN2 regulated ER morphology, whereas ERMIT2 localized at the ER-mitochondria interface and interacted with mitochondrial mitofusins to tether ER and mitochondria. This tethering allowed efficient mitochondrial calcium ion uptake and phospholipid transfer. Expression of ERMIT2 ameliorated the ER stress, inflammation, and fibrosis typical of liver-specific Mfn2 knockout mice. Thus, ER-specific MFN2 variants display entirely extramitochondrial MFN2 functions involved in interorganellar tethering and liver metabolic activities.

A splicing transcriptome-wide association study identifies novel altered splicing for Alzheimer's disease susceptibility

Alzheimer's disease (AD) is a common neurodegenerative disease in aging individuals. Alternative splicing is reported to be relevant to AD development while their roles in etiology of AD remain largely elusive. We performed a comprehensive splicing transcriptome-wide association study (spTWAS) using intronic excision expression genetic prediction models of 12 brain tissues developed through three modelling strategies, to identify candidate susceptibility splicing introns for AD risk. A total of 111,326 (46,828 proxy) cases and 677,663 controls of European ancestry were studied. We identified 343 associations of 233 splicing introns (143 genes) with AD risk after Bonferroni correction (0.05/136,884 = 3.65 × 10^-7). Fine-mapping analyses supported 155 likely causal associations corresponding to 83 splicing introns of 55 genes. Eighteen causal splicing introns of 15 novel genes (EIF2D, WDR33, SAP130, BYSL, EPHB6, MRPL43, VEGFB, PPP1R13B, TLN2, CLUHP3, LRRC37A4P, CRHR1, LINC02210, ZNF45-AS1, and XPNPEP3) were identified for the first time to be related to AD susceptibility. Our study identified novel genes and splicing introns associated with AD risk, which can improve our understanding of the etiology of AD.

The Splicing Factor OsSCL26 Regulates Phosphorus Homeostasis in Rice

Phosphorus (P) is an essential nutrient for plant growth. However, its deficiency poses a significant challenge for crop production. To overcome the low P availability, plants have developed various strategies to regulate their P uptake and usage. In this study, we identified a splicing factor, OsSCL26, belonging to the Serine/arginine-rich (SR) proteins, that plays a crucial role in regulating P homeostasis in rice. OsSCL26 is expressed in the roots, leaves, and base nodes, with higher expression levels observed in the leaf blades during the vegetative growth stage. The OsSCL26 protein is localized in the nucleus. Mutation of OsSCL26 resulted in the accumulation of P in the shoot compared to the wild-type, and the dwarf phenotype of the osscl26 mutant was alleviated under low P conditions. Further analysis revealed that the accumulated P concentrations in the osscl26 mutant were higher in the old leaves and lower in the new leaves. Furthermore, the P-related genes, including the PHT and SPX family genes, were upregulated in the osscl26 mutant, and the exclusion/inclusion ratio of the two genes, OsSPX-MFS2 and OsNLA2, was increased compared to wild-type rice. These findings suggest that the splicing factor OsSCL26 plays a pivotal role in maintaining P homeostasis in rice by influencing the absorption and distribution of P through the regulation of the transcription and splicing of the P transport genes.

Transcription factors and splice factors - interconnected regulators of stem cell differentiation

Purpose of review

The underlying molecular mechanisms that direct stem cell differentiation into fully functional, mature cells remain an area of ongoing investigation. Cell state is the product of the combinatorial effect of individual factors operating within a coordinated regulatory network. Here, we discuss the contribution of both gene regulatory and splicing regulatory networks in defining stem cell fate during differentiation and the critical role of protein isoforms in this process.

Recent findings

We review recent experimental and computational approaches that characterize gene regulatory networks, splice regulatory networks, and the resulting transcriptome and proteome they mediate during differentiation. Such approaches include long-read RNA sequencing, which has demonstrated high-resolution profiling of mRNA isoforms, and Cas13-based CRISPR, which could make possible high-throughput isoform screening. Collectively, these developments enable systems-level profiling of factors contributing to cell state.

Summary

Overall, gene and splice regulatory networks are important in defining cell state. The emerging high-throughput systems-level approaches will characterize the gene regulatory network components necessary in driving stem cell differentiation.

The Glomerulus Multiomics Analysis Provides Deeper Insights into Diabetic Nephropathy

Although diabetic nephropathy (DN) is the leading cause of the end-stage renal disease, the exact regulation mechanisms remain unknown. In this study, we integrated the transcriptomics and proteomics profiles of glomeruli isolated from 50 biopsy-proven DN patients and 25 controls to investigate the latest findings about DN pathogenesis. First, 1152 genes exhibited differential expression at the mRNA or protein level, and 364 showed significant association. These strong correlated genes were divided into four different functional modules. Moreover, a regulatory network of the transcription factors (TFs)-target genes (TGs) was constructed, with 30 TFs upregulated at the protein levels and 265 downstream TGs differentially expressed at the mRNA levels. These TFs are the integration centers of several signal transduction pathways and have tremendous therapeutic potential for regulating the aberrant production of TGs and the pathological process of DN. Furthermore, 29 new DN-specific splice-junction peptides were discovered with high confidence; these peptides may play novel functions in the pathological course of DN. So, our in-depth integrative transcriptomics-proteomics analysis provided deeper insights into the pathogenesis of DN and opened the potential avenue for finding new therapeutic interventions. MS raw files were deposited into the proteomeXchange with the dataset identifier PXD040617.

Catsnap: a user-friendly algorithm for determining the conservation of protein variants reveals extensive parallelisms in the evolution of alternative splicing

Understanding the evolutionary conservation of complex eukaryotic transcriptomes significantly illuminates the physiological relevance of alternative splicing (AS). Examining the evolutionary depth of a given AS event with ordinary homology searches is generally challenging and time-consuming. Here, we present Catsnap, an algorithmic pipeline for assessing the conservation of putative protein isoforms generated by AS. It employs a machine learning approach following a database search with the provided pair of protein sequences. We used the Catsnap algorithm for analyzing the conservation of emerging experimentally characterized alternative proteins from plants and animals. Indeed, most of them are conserved among other species. Catsnap can detect the conserved functional protein isoforms regardless of the AS type by which they are generated. Notably, we found that while the primary amino acid sequence is maintained, the type of AS determining the inclusion or exclusion of protein regions varies throughout plant phylogenetic lineages in these proteins. We also document that this phenomenon is less seen among animals. In sum, our algorithm highlights the presence of unexpectedly frequent hotspots where protein isoforms recurrently arise to carry physiologically relevant functions. The user web interface is available at https://catsnap.cesnet.cz/.

Fully sequencing the cassava full-length cDNA library reveals unannotated transcript structures and alternative splicing events in regions with a high density of single nucleotide variations, insertions-deletions, and heterozygous sequences

The primary transcript structure provides critical insights into protein diversity, transcriptional modification, and functions. Cassava transcript structures are highly diverse because of alternative splicing (AS) events and high heterozygosity. To precisely determine and characterize transcript structures, fully sequencing cloned transcripts is the most reliable method. However, cassava annotations were mainly determined according to fragmentation-based sequencing analyses (e.g., EST and short-read RNA-seq). In this study, we sequenced the cassava full-length cDNA library, which included rare transcripts. We obtained 8,628 non-redundant fully sequenced transcripts and detected 615 unannotated AS events and 421 unannotated loci. The different protein sequences resulting from the unannotated AS events tended to have diverse functional domains, implying that unannotated AS contributes to the truncation of functional domains. The unannotated loci tended to be derived from orphan genes, implying that the loci may be associated with cassava-specific traits. Unexpectedly, individual cassava transcripts were more likely to have multiple AS events than Arabidopsis transcripts, suggestive of the regulated interactions between cassava splicing-related complexes. We also observed that the unannotated loci and/or AS events were commonly in regions with abundant single nucleotide variations, insertions-deletions, and heterozygous sequences. These findings reflect the utility of completely sequenced FLcDNA clones for overcoming cassava-specific annotation-related problems to elucidate transcript structures. Our work provides researchers with transcript structural details that are useful for annotating highly diverse and unique transcripts and alternative splicing events.

Transcriptomic and Functional Analyses Reveal the Different Roles of Vitamins C, E, and K in Regulating Viral Infections in Maize

Maize lethal necrosis (MLN), one of the most important maize viral diseases, is caused by maize chlorotic mottle virus (MCMV) infection in combination with a potyvirid, such as sugarcane mosaic virus (SCMV). However, the resistance mechanism of maize to MLN remains largely unknown. In this study, we obtained isoform expression profiles of maize after SCMV and MCMV single and synergistic infection (S + M) via comparative analysis of SMRT- and Illumina-based RNA sequencing. A total of 15,508, 7567, and 2378 differentially expressed isoforms (DEIs) were identified in S + M, MCMV, and SCMV libraries, which were primarily involved in photosynthesis, reactive oxygen species (ROS) scavenging, and some pathways related to disease resistance. The results of virus-induced gene silencing (VIGS) assays revealed that silencing of a vitamin C biosynthesis-related gene, ZmGalDH or ZmAPX1, promoted viral infections, while silencing ZmTAT or ZmNQO1, the gene involved in vitamin E or K biosynthesis, inhibited MCMV and S + M infections, likely by regulating the expressions of pathogenesis-related (PR) genes. Moreover, the relationship between viral infections and expression of the above four genes in ten maize inbred lines was determined. We further demonstrated that the exogenous application of vitamin C could effectively suppress viral infections, while vitamins E and K promoted MCMV infection. These findings provide novel insights into the gene regulatory networks of maize in response to MLN, and the roles of vitamins C, E, and K in conditioning viral infections in maize.

Deciphering the role of alternative splicing as modulators of defense response in the MYMIV- Vigna mungo pathosystem

Alternative splicing (AS) is a crucial regulatory mechanism that impacts transcriptome and proteome complexity under stressful situations. Although its role in abiotic stresses is somewhat understood, our understanding of the mechanistic regulation of pre-mRNA splicing in plant-pathogen interaction is meagre. To comprehend this unexplored immune reprogramming mechanism, transcriptome profiles of Mungbean Yellow Mosaic India Virus (MYMIV)-resistant and susceptible Vigna mungo genotypes were analysed for AS genes that may underlie the resistance mechanism. Results revealed a repertoire of AS-isoforms accumulated during pathogenic infestation, with intron retention being the most common AS mechanism. Identification of 688 differential alternatively spliced (DAS) genes in the resistant host elucidates its robust antiviral response, whereas 322 DAS genes were identified in the susceptible host. Enrichment analyses confirmed DAS transcripts pertaining to stress, signalling, and immune system pathways have undergone maximal perturbations. Additionally, a strong regulation of the splicing factors has been observed both at transcriptional and post-transcriptional levels. qPCR validation of candidate DAS transcripts with induced expression upon MYMIV-infection demonstrated a competent immune response in the resistant background. The AS-impacted genes resulted either in partial/complete loss of functional domains or altered sensitivity to miRNA-mediated gene silencing. A complex regulatory module, miR7517-ATAF2, has been identified in an aberrantly spliced ATAF2 isoform that exposes an intronic miR7517 binding site, thereby suppressing the negative regulator to enhance defense reaction. The present study establishes AS as a non-canonical immune reprogramming mechanism that operates in parallel, thereby offering an alternative strategy for developing yellow mosaic-resistant V. mungo cultivars. This article is protected by copyright. All rights reserved.

To Explore the Binding Affinity of Human γ-Secretase Activating Protein (GSAP) Isoform 4 with APP-C99 Peptides

γ-Secretase activating protein (GSAP) is known to play an important role in the β-amyloid pathway. It acts as a modulator and accentuates the truncation of the amyloid precursor protein C-99 fragment through the γ-secretase complex. GSAP has four isoforms, out of which canonical isoform 1, a 16 kDa C-terminal portion, has been extensively studied, whereas the function of other three isoforms remains unknown. Here, we explore the GSAP isoform 4 (GSAP_I4) expression and purification from inclusion bodies followed by the refolding of the protein. The secondary structure of GSAP_I4 is predicted using circular dichroism. The protein is further characterized by western blotting and mass spectroscopy analysis. Additionally, biochemical assays and in silico molecular docking and molecular simulation are performed to investigate the binding of GSAP_I4 and APP-C99 peptide fragments. The results reflect that although GSAP_I1 and GSAP_I4 share high sequence similarity, the isoform 4 does not show any affinity toward APP-C99 peptide fragments. This hints toward the fact that GSAP_I4 might have a different role in the living system that is yet unexplored.

Contribution of A-to-I RNA editing, M6A RNA Methylation, and Alternative Splicing to physiological brain aging and neurodegenerative diseases

Aging is a physiological and progressive phenomenon in all organisms' life cycle, characterized by the accumulation of degenerative processes triggered by several alterations within molecular pathways. These changes compromise cell fate, resulting in the loss of functions in tissues throughout the body, including the brain. Physiological brain aging has been linked to structural and functional alterations, as well as to an increased risk of neurodegenerative diseases. Post-transcriptional RNA modifications modulate mRNA coding properties, stability, translatability, expanding the coding capacity of the genome, and are involved in all cellular processes. Among mRNA post-transcriptional modifications, the A-to-I RNA editing, m6A RNA Methylation and Alternative Splicing play a critical role in all the phases of a neuronal cell life cycle and alterations in their mechanisms of action significantly contribute to aging and neurodegeneration. Here we review our current understanding of the contribution of A-to-I RNA editing, m6A RNA Methylation, and Alternative Splicing to physiological brain aging process and neurodegenerative diseases.

Identification of sex-specific splicing via comparative transcriptome analysis of gonads from sea cucumber Apostichopus japonicus

Alternative splicing (AS) is an essential post-transcriptional regulation mechanism for sex differentiation and gonadal development, which has rarely been reported in marine invertebrates. Sea cucumber (Apostichopus japonicus) is an economically important marine benthic echinoderm with a potential XX/XY sex determination mechanism, whose molecular mechanism in the gonadal differentiation has not been clearly understood. In this study, we analyzed available RNA-seq datasets of male and female gonads to explore if AS mechanism exerts an essential function in sex differentiation and gonadal development of A. japonicus. In our results, a total of 20,338 AS events from 7219 alternatively spliced genes, and 189 sexually differential alternative splicing (DAS) events from 156 genes were identified in gonadal transcriptome of sea cucumber. Gene Ontology analysis indicated that these DAS genes were significantly enriched in spermatogenesis-related GO terms. Maximal Clique Centrality (MCC) was then applied for protein-protein interaction (PPI) analysis to search for protein interactions and hub DAS gene. Among all DAS genes, we identified 10 DAS genes closely related to spermatogenesis and (or) sperm motility and a hub gene dnah1. Thus, this study revealed that alternative isoforms were generated from certain genes in female and male gonads through alternative splicing, which may provide direct evidence that alternative splicing mechanisms participate in female and male gonads. These results suggested a novel perspective for explaining the molecular mechanisms underlying gonadal differentiation between male and female sea cucumbers.

Re-evaluating the impact of alternative RNA splicing on proteomic diversity

Alternative splicing (AS) constitutes a mechanism by which protein-coding genes and long non-coding RNA (lncRNA) genes produce more than a single mature transcript. From plants to humans, AS is a powerful process that increases transcriptome complexity. Importantly, splice variants produced from AS can potentially encode for distinct protein isoforms which can lose or gain specific domains and, hence, differ in their functional properties. Advances in proteomics have shown that the proteome is indeed diverse due to the presence of numerous protein isoforms. For the past decades, with the help of advanced high-throughput technologies, numerous alternatively spliced transcripts have been identified. However, the low detection rate of protein isoforms in proteomic studies raised debatable questions on whether AS contributes to proteomic diversity and on how many AS events are really functional. We propose here to assess and discuss the impact of AS on proteomic complexity in the light of the technological progress, updated genome annotation, and current scientific knowledge.

Transcriptome diversity assessment of Gossypium arboreum (FDH228) leaves under control, drought and whitefly infestation using PacBio long reads

Alternative splicing (AS) and alternative polyadenylation (APA) are common mechanisms in eukaryotes to increase the complexity of transcriptomes and subsequently proteomes. Analysis of long reads transcriptomics data can result in the discovery of novel transcripts, splice sites, AS or APA events. Gossypium arboreum is an important cultivated cotton species and a putative contributor of the A sub-genome to the modern tetraploid cotton; and inherently tolerant to several biotic and abiotic stresses. Specifically, its variety 'FDH228' is considered to be an important resistance source. In this study, we sequenced the G. arboreum (var. FDH228) transcriptome using PacBio IsoSeq and illumina short read sequencing under three different conditions i.e. untreated/healthy, treated with biotic stress through whitefly infestation, and treated with abiotic stress via water deprivation, for the discovery and surveying of canonical and non-canonical AS, APA and transcript fusion events. We were able to obtain 15,419 unique transcripts from all samples representing 11,343 genes, out of which 10,832 were annotated and 520 were novel with respect to the published reference genome. These transcripts were grouped into different structural categories including 60 Antisense, 11,959 having a full-splice match, 999 with incomplete-splice match, 30 fusion transcripts, 177 genic, 479 intergenic, 771 novels in the catalog, and 944 Novel but not found in the catalog. Subsequently, randomly selected candidate transcripts were experimentally validated using qRT-PCR. Our comprehensive identification of canonical and non-canonical splicing events, and novel and fusion transcripts aids in the understanding of the resistance mechanisms for this specific germplasm.

Chilo suppressalis heat shock proteins are regulated by heat shock factor 1 during heat stress

Heat shock factor 1 (HSF1) functions to maintain cellular and organismal homeostasis by regulating the expression of target genes, including those encoding heat shock proteins (HSPs). In the present study, the gene encoding HSF1 was cloned from the rice pest Chilo suppressalis, and designated Cshsf1. The deduced protein product, CsHSF1, contained conserved domains typical of the HSF1 family, including a DNA-binding domain, two hydrophobic heptad repeat domains, and a C-terminal transactivation domain. Real-time quantitative PCR showed that Cshsf1 was highly expressed in hemocytes. Expression analysis in different developmental stages of C. suppressalis revealed that Cshsf1 was most highly expressed in male adults. RNAi-mediated silencing of Cshsf1 expression reduced C. suppressalis survival at high temperatures. To investigate the regulatory interactions between Cshsf1 and Cshsps, the promoters and expression patterns of 18 identified Cshsps in C. suppressalis were analysed; four types of heat shock elements (HSEs) were identified in promoter regions including canonical, tail-tail, head-head, and step/gap. The expression of Cshsp19.0, Cshsp21.7B, Cshsp60, Cshsp70 and Cshsp90 was positively regulated by Cshsf1; however, Cshsp22.8, Cshsp702, Cshsp705 and Cshsp706 gene expression was not altered. This study provides a foundation for future studies of HSF1 in insects during thermal stress.

Alternative splicing: transcriptional regulatory network in agroforestry

Alternative splicing (AS) in plants plays a key role in regulating the expression of numerous transcripts from a single gene in a regulatory pathway. Variable concentrations of growth regulatory hormones and external stimuli trigger alternative splicing to switch among different growth stages and adapt to environmental stresses. In the AS phenomenon, a spliceosome causes differential transcriptional modifications in messenger RNA (mRNAs), resulting in partial or complete retention of one or more introns as compared to fully spliced mRNA. Differentially expressed proteins translated from intron-retaining messenger RNA (mRNA^ir) perform vital functions in the feedback mechanism. At the post-transcriptional level, AS causes the remodeling of transcription factors (TFs) by the addition or deletion of binding domains to activate and/or repress transcription. In this study, we have summarized the specific role of AS in the regulation of gene expression through repression and activation of the transcriptional regulatory network under external stimuli and switch among developmental stages.

Isoform changes of action potential regulators in the ventricles of arrhythmogenic phospholamban-R14del humanized mouse hearts

Arrhythmogenic cardiomyopathy (ACM) is characterized by life-threatening ventricular arrhythmias and sudden cardiac death and affects hundreds of thousands of patients worldwide. The deletion of Arginine 14 (p.R14del) in the phospholamban (PLN) gene has been implicated in the pathogenesis of ACM. PLN is a key regulator of sarcoplasmic reticulum (SR) Ca2+ cycling and cardiac contractility. Despite global gene and protein expression studies, the molecular mechanisms of PLN-R14del ACM pathogenesis remain unclear. Using a humanized PLN-R14del mouse model and human induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs), we investigated the transcriptome-wide mRNA splicing changes associated with the R14del mutation. We identified >200 significant alternative splicing (AS) events and distinct AS profiles were observed in the right (RV) and left (LV) ventricles in PLN-R14del compared to WT mouse hearts. Enrichment analysis of the AS events showed that the most affected biological process was associated with "cardiac cell action potential", specifically in the RV. We found that splicing of 2 key genes, Trpm4 and Camk2d, which encode proteins regulating calcium homeostasis in the heart, were altered in PLN-R14del mouse hearts and human iPSC-CMs. Bioinformatical analysis pointed to the tissue-specific splicing factors Srrm4 and Nova1 as likely upstream regulators of the observed splicing changes in the PLN-R14del cardiomyocytes. Our findings suggest that aberrant splicing may affect Ca2+-homeostasis in the heart, contributing to the increased risk of arrythmogenesis in PLN-R14del ACM.

Reovirus μ2 Protein Impairs Translation to Reduce U5 snRNP Protein Levels

Mammalian orthoreovirus (MRV) is a double-stranded RNA virus from the Reoviridae family that infects a large range of mammals, including humans. Recently, studies have shown that MRV alters cellular alternative splicing (AS) during viral infection. The structural protein μ2 appears to be the main determinant of these AS modifications by decreasing the levels of U5 core components EFTUD2, PRPF8, and SNRNP200 during infection. In the present study, we investigated the mechanism by which μ2 exerts this effect on the U5 components. Our results revealed that μ2 has no impact on steady-state mRNA levels, RNA export, and protein stability of these U5 snRNP proteins. However, polysome profiling and metabolic labeling of newly synthesized proteins revealed that μ2 exerts an inhibitory effect on global translation. Moreover, we showed that μ2 mutants unable to accumulate in the nucleus retain most of the ability to reduce PRPF8 protein levels, indicating that the effect of μ2 on U5 snRNP components mainly occurs in the cytoplasm. Finally, co-expression experiments demonstrated that μ2 suppresses the expression of U5 snRNP proteins in a dose-dependent manner, and that the expression of specific U5 snRNP core components have different sensitivities to μ2's presence. Altogether, these results suggest a novel mechanism by which the μ2 protein reduces the levels of U5 core components through translation inhibition, allowing this viral protein to alter cellular AS during infection.

The complex, dynamic SpliceOme of the small GTPase transcripts altered by technique, sex, genetics, tissue specificity, and RNA base editing

The small GTPase family is well-studied in cancer and cellular physiology. With 162 annotated human genes, the family has a broad expression throughout cells of the body. Members of the family have multiple exons that require splicing. Yet, the role of splicing within the family has been underexplored. We have studied the splicing dynamics of small GTPases throughout 41,671 samples by integrating Nanopore and Illumina sequencing techniques. Within this work, we have made several discoveries. 1). Using the GTEx long read data of 92 samples, each small GTPase gene averages two transcripts, with 83 genes (51%) expressing two or more isoforms. 2). Cross-tissue analysis of GTEx from 17,382 samples shows 41 genes (25%) expressing two or more protein-coding isoforms. These include protein-changing transcripts in genes such as RHOA, RAB37, RAB40C, RAB4B, RAB5C, RHOC, RAB1A, RAN, RHEB, RAC1, and KRAS. 3). The isolation and library technique of the RNAseq influences the abundance of non-sense-mediated decay and retained intron transcripts of small GTPases, which are observed more often in genes than appreciated. 4). Analysis of 16,243 samples of "Blood PAXgene" identified seven genes (3.7%; RHOA, RAB40C, RAB4B, RAB37, RAB5B, RAB5C, RHOC) with two or more transcripts expressed as the major isoform (75% of the total gene), suggesting a role of genetics in altering splicing. 5). Rare (ARL6, RAB23, ARL13B, HRAS, NRAS) and common variants (GEM, RHOC, MRAS, RAB5B, RERG, ARL16) can influence splicing and have an impact on phenotypes and diseases. 6). Multiple genes (RAB9A, RAP2C, ARL4A, RAB3A, RAB26, RAB3C, RASL10A, RAB40B, and HRAS) have sex differences in transcript expression. 7). Several exons are included or excluded for small GTPase genes (RASEF, KRAS, RAC1, RHEB, ARL4A, RHOA, RAB30, RHOBTB1, ARL16, RAP1A) in one or more forms of cancer. 8). Ten transcripts are altered in hypoxia (SAR1B, IFT27, ARL14, RAB11A, RAB10, RAB38, RAN, RIT1, RAB9A) with RHOA identified to have a transient 3'UTR RNA base editing at a conserved site found in all of its transcripts. Overall, we show a remarkable and dynamic role of splicing within the small GTPase family that requires future explorations.

Efficient Minimum Flow Decomposition via Integer Linear Programming

Minimum flow decomposition (MFD) is an NP-hard problem asking to decompose a network flow into a minimum set of paths (together with associated weights). Variants of it are powerful models in multiassembly problems in Bioinformatics, such as RNA assembly. Owing to its hardness, practical multiassembly tools either use heuristics or solve simpler, polynomial time-solvable versions of the problem, which may yield solutions that are not minimal or do not perfectly decompose the flow. Here, we provide the first fast and exact solver for MFD on acyclic flow networks, based on Integer Linear Programming (ILP). Key to our approach is an encoding of all the exponentially many solution paths using only a quadratic number of variables. We also extend our ILP formulation to many practical variants, such as incorporating longer or paired-end reads, or minimizing flow errors. On both simulated and real-flow splicing graphs, our approach solves any instance in <13 seconds. We hope that our formulations can lie at the core of future practical RNA assembly tools. Our implementations are freely available on Github.

Single-molecule Real-time (SMRT) Sequencing Facilitates Transcriptome Research and Genome Annotation of the Fish Sillago sinica

As a newly described Sillaginidae species, Chinese sillago (Sillago sinica) needs a better understanding of gene annotation information. In this study, we reported the first full-length transcriptome data of S. sinica using the PacBio isoform sequencing Iso-seq and a description of transcriptome structure analysis. A total of 454,979 high-quality full-length transcripts were obtained by single-molecule real-time (SMRT) sequencing, which was corrected by Illumina sequencing data. After that, 66,948 non-redundant full-length transcripts were generated after mapping to the reference genome of S. sinica, including 49 fusion isoforms and 9,250 novel isoforms. 63,459 isoforms were successfully annotated by one of the Nr, Nt, SwissProt, Pfam, KOG, GO, and KEGG databases. Additionally, 30,987 alternative polyadenylation (APA) sites, 451,867 alternative splicing (AS) events, 21,928 long non-coding RNAs (lncRNAs) and 12,911 transcription factors (TFs) were identified. The full-length transcripts of S. sinica would provide a precious resource for characterizing the transcriptome of S. sinica and for the further study of gene function and regulatory mechanism of this species.

Transcription Factor GATA4 Regulates Cell Type-Specific Splicing Through Direct Interaction With RNA in Human Induced Pluripotent Stem Cell-Derived Cardiac Progenitors

BACKGROUND

GATA4 (GATA-binding protein 4), a zinc finger-containing, DNA-binding transcription factor, is essential for normal cardiac development and homeostasis in mice and humans, and mutations in this gene have been reported in human heart defects. Defects in alternative splicing are associated with many heart diseases, yet relatively little is known about how cell type- or cell state-specific alternative splicing is achieved in the heart. Here, we show that GATA4 regulates cell type-specific splicing through direct interaction with RNA and the spliceosome in human induced pluripotent stem cell-derived cardiac progenitors.

METHODS

We leveraged a combination of unbiased approaches including affinity purification of GATA4 and mass spectrometry, enhanced cross-linking with immunoprecipitation, electrophoretic mobility shift assays, in vitro splicing assays, and unbiased transcriptomic analysis to uncover GATA4's novel function as a splicing regulator in human induced pluripotent stem cell-derived cardiac progenitors.

RESULTS

We found that GATA4 interacts with many members of the spliceosome complex in human induced pluripotent stem cell-derived cardiac progenitors. Enhanced cross-linking with immunoprecipitation demonstrated that GATA4 also directly binds to a large number of mRNAs through defined RNA motifs in a sequence-specific manner. In vitro splicing assays indicated that GATA4 regulates alternative splicing through direct RNA binding, resulting in functionally distinct protein products. Correspondingly, knockdown of GATA4 in human induced pluripotent stem cell-derived cardiac progenitors resulted in differential alternative splicing of genes involved in cytoskeleton organization and calcium ion import, with functional consequences associated with the protein isoforms.

CONCLUSIONS

This study shows that in addition to its well described transcriptional function, GATA4 interacts with members of the spliceosome complex and regulates cell type-specific alternative splicing via sequence-specific interactions with RNA. Several genes that have splicing regulated by GATA4 have functional consequences and many are associated with dilated cardiomyopathy, suggesting a novel role for GATA4 in achieving the necessary cardiac proteome in normal and stress-responsive conditions.

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.

SAPFIR: A webserver for the identification of alternative protein features

Background

Alternative splicing can increase the diversity of gene functions by generating multiple isoforms with different sequences and functions. However, the extent to which splicing events have functional consequences remains unclear and predicting the impact of splicing events on protein activity is limited to gene-specific analysis.

Results

To accelerate the identification of functionally relevant alternative splicing events we created SAPFIR, a predictor of protein features associated with alternative splicing events. This webserver tool uses InterProScan to predict protein features such as functional domains, motifs and sites in the human and mouse genomes and link them to alternative splicing events. Alternative protein features are displayed as functions of the transcripts and splice sites. SAPFIR could be used to analyze proteins generated from a single gene or a group of genes and can directly identify alternative protein features in large sequence data sets. The accuracy and utility of SAPFIR was validated by its ability to rediscover previously validated alternative protein domains. In addition, our de novo analysis of public datasets using SAPFIR indicated that only a small portion of alternative protein domains was conserved between human and mouse, and that in human, genes involved in nervous system process, regulation of DNA-templated transcription and aging are more likely to produce isoforms missing functional domains due to alternative splicing.

Conclusion

Overall SAPFIR represents a new tool for the rapid identification of functional alternative splicing events and enables the identification of cellular functions affected by a defined splicing program. SAPFIR is freely available at https://bioinfo-scottgroup.med.usherbrooke.ca/sapfir/, a website implemented in Python, with all major browsers supported. The source code is available at https://github.com/DelongZHOU/SAPFIR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04804-w.