Transcriptome analysis reveals the complexity of alternative splicing regulation in the fungus Verticillium dahliae

Beyond ribosomal function: RPS6 deficiency suppresses cholangiocarcinoma cell growth by disrupting alternative splicing

Cholangiocarcinoma (CCA) is a bile duct malignancy with a dismal prognosis. This study systematically investigated the role of the ribosomal protein S6 (RPS6) gene, which is dependent in CCA. We found that RPS6 upregulation in CCA tissues was correlated with a poor prognosis. Functional investigations have shown that alterations in RPS6 expression, both gain- and loss-of function could affect the proliferation of CCA cells. In xenograft tumor models, RPS6 overexpression enhances tumorigenicity, whereas RPS6 silencing reduces it. Integration analysis using RNA-seq and proteomics elucidated downstream signaling pathways of RPS6 depletion by affecting the cell cycle, especially DNA replication. Immunoprecipitation followed by mass spectrometry has identified numerous spliceosome complex proteins associated with RPS6. Transcriptomic profiling revealed that RPS6 affects numerous alternative splicing (AS) events, and combined with RNA immunoprecipitation sequencing, revealed that minichromosome maintenance complex component 7 (MCM7) binds to RPS6, which regulates its AS and increases oncogenic activity in CCA. Targeting RPS6 with vivo phosphorodiamidate morpholino oligomer (V-PMO) significantly inhibited the growth of CCA cells, patient-derived organoids, and subcutaneous xenograft tumor. Taken together, the data demonstrate that RPS6 is an oncogenic regulator in CCA and that RPS6-V-PMO could be repositioned as a promising strategy for treating CCA.

Analysis of different expression RNA binding protein genes in mouse microglia cell from the brains of mice 72 h after subarachnoid hemorrhage or sham operation

BACKGROUND

The prognosis of brain injury caused by subarachnoid hemorrhage (SAH) is poor. Previous studies showed that abnormal function of RBPs might be involved in brain injury, neuroinflammation and further affect microglia homeostasis. However, no studies have systematically analyzed the genome-wide abnormal expression of RBPs genes in microglia during SAH.

METHODS

RNA-seq data of microglia from the SAH mouse group (SAH) and control sham-operated mouse group (sham) were downloaded from the GEO database in GSE167957, including four samples from the sham group and four samples from the SAH group for subsequent analysis.Utilizing GO and KEGG functional enrichment analyses, we conducted a comprehensive study of differentially expressed genes (DEGs), alternative splicing patterns, and co-expression networks to gain deeper insights into the differential expression of RNA-binding proteins (RBPs) and differential alternative splicing events (ASEs) between the SAH (subarachnoid hemorrhage) and sham groups. This analysis aimed to elucidate the potential mechanisms underlying the aberrant expression of RBPs in microglia during brain injury caused by SAH.

RESULTS

ASEs and co-expression analyses of differentially expressed RBPs and differential ASEs were carried out in microglia in terms of gene expression. GO and KEGG functional enrichment analysis showed that aberrantly expressed RBPs such as Mcm7, Mtdh, SRSF3, and Hnrnpa2b1 may affect and regulate downstream Csnk1d, Uckl1 and other protein phosphorylation-related genes by alterative splicing.

CONCLUSION

RBPs were aberrantly expressed in microglia during the development of brain injury secondary to SAH, regulating alterative splicing of downstream genes and influencing the progression of SAH brain injury in this study. This implies that RBPs are important for the identification of new therapeutic targets for brain injury after SAH.

Downregulation of Splicing Factor PTBP1 Curtails FBXO5 Expression to Promote Cellular Senescence in Lung Adenocarcinoma

Polypyrimidine tract-binding protein 1 (PTBP1) plays an essential role in splicing and post-transcriptional regulation. Moreover, PTBP1 has been implicated as a causal factor in tumorigenesis. However, the involvement of PTBP1 in cellular senescence, a key biological process in aging and cancer suppression, remains to be clarified. Here, it is shown that PTBP1 is associated with the facilitation of tumor growth and the prognosis in lung adenocarcinoma (LUAD). PTBP1 exhibited significantly increased expression in various cancer types including LUAD and showed consistently decreased expression in multiple cellular senescence models. Suppression of PTBP1 induced cellular senescence in LUAD cells. In terms of molecular mechanisms, the silencing of PTBP1 enhanced the skipping of exon 3 in F-box protein 5 (FBXO5), resulting in the generation of a less stable RNA splice variant, FBXO5-S, which subsequently reduces the overall FBXO5 expression. Additionally, downregulation of FBXO5 was found to induce senescence in LUAD. Collectively, these findings illustrate that PTBP1 possesses an oncogenic function in LUAD through inhibiting senescence, and that targeting aberrant splicing mediated by PTBP1 has therapeutic potential in cancer treatment.

eEF1A1 regulates the expression and alternative splicing of genes associated with Parkinson's disease in U251 cells

BACKGROUND

Eukaryotic elongation factor 1A1 (eEF1A1) is an RNA-binding protein that is associated with PARK2 activity in cells, suggesting a possible role in Parkinson's disease (PD).

OBJECTIVE

To clear whether eEF1A1 plays a role in PD through transcriptional or posttranscriptional regulation.

METHODS

The GSE68719 dataset was downloaded from the GEO database, and the RNA-seq data of all brain tissue autopsies were obtained from 29 PD patients and 44 neurologically normal control subjects. To inhibit eEF1A1 from being expressed in U251 cells, siRNA was transfected into those cells, and RNA-seq high-throughput sequencing was used to determine the differentially expressed genes (DEGs) and differentially alternative splicing events (ASEs) resulting from eEF1A1 knockdown.

RESULTS

eEF1A1 was significantly overexpressed in PD brain tissue in the BA9 area. GO and KEGG enrichment analyses revealed that eEF1A1 knockdown significantly upregulated the expression of the genes CXCL10, NGF, PTX3, IL6, ST6GALNAC3, NUPR1, TNFRSF21, and CXCL2 and upregulated the alternative splicing of the genes ACOT7, DDX10, SHMT2, MYEF2, and NDUFAF5. These genes were enriched in pathways related to PD pathogenesis, such as apoptosis, inflammatory response, and mitochondrial dysfunction.

CONCLUSION

The results suggesting that eEF1A1 involved in the development of PD by regulating the differential expression and alternative splicing of genes, providing a theoretical basis for subsequent research.

Effects of myo-inositol on regulating glucose and lipid metabolism and alternative splicing events coexpressed with lncRNAs in the liver tissues of diabetic mice

Objective

Recent studies have shown that gene alternative splicing (AS) and long noncoding RNAs (lncRNAs) are involved in diabetes mellitus (DM) and its complications. Currently, myo-inositol (MI) is considered as effective for the treatment of insulin resistance and lipid metabolism disorders in diabetes patients. We hope to better explore the potential roles of gene AS and lncRNAs in liver glucose and lipid metabolism in diabetes, as well as the effects of myo-inositol treatment, through transcriptome analysis.

Methods

This study analysed glucose and lipid metabolism-related biochemical indicators and liver HE staining in four groups of mice: the control group (Ctrl group), the diabetes group (DM group), the myo-inositol treatment group (MI group), and the metformin treatment group (Met group). The changes in relevant gene-regulated alternative splicing events (RASEs) and lncRNAs were analysed by RNA sequencing of liver tissue, and coexpression analysis and functional enrichment analysis were used to predict the possible lncRNAs and RASEs involved in liver glucose and lipid metabolism.

Result

Metformin and myo-inositol alleviated insulin resistance, lipid metabolism disorders, and hepatic steatosis in diabetic mice. Transcriptome sequencing analysis revealed differential splicing events of genes related to lipid metabolism and differentially expressed lncRNAs (DElncRNAs). Six different lncRNAs and their potentially interacting splicing events were predicted.

Conclusion

The present study revealed novel changes in RASEs and lncRNAs in the livers of diabetic mice following treatment with myo-inositol, which may shed light on the potential mechanisms by which myo-inositol delays and treats the progression of hepatic glucose and lipid metabolism in diabetes.

Alteration of RNA m6A methylation mediates aberrant RNA binding protein expression and alternative splicing in condyloma acuminatum

Background

Condyloma acuminatum (CA) is caused by low-risk human papillomavirus, and is characterized by high recurrence after treatment. The RNA modification N6-methyladenosine (m6A) plays an important role during diverse viral infections, including high-risk HPV infection in cervical cancer. However, it is unclear whether low-risk HPV infection changes the RNA m6A methylation in CA.

Methods

High-throughputm6A-sequencing was performed to profile the transcriptome-wide mRNA modifications of CA tissues infected by LR-HPVs and the paired normal tissues from CA patients. We further investigated the regulation of alternative splicing by RNA binding proteins (RBPs) with altered m6A modification and constructed a regulatory network among these RBPs, regulated alternative splicing events (RASEs) and regulated alternative splicing genes (RASGs) in CA.

Results

The results show that the m6A level in CA tissues differed from that in the paired controls. Furthermore, cell cycle- and cell adhesion- associated genes with m6A modification were differentially expressed in CA tissues compared to the paired controls. In particular, seven RNA binding protein genes with specific m6A methylated sites, showed a higher or lower expression at the mRNA level in CA tissues than in the paired normal tissues. In addition, these differentially expressed RNA binding protein genes would regulate the alternative splicing pattern of apoptotic process genes in CA tissue.

Conclusions

Our study reveals a sophisticated m6A modification profile in CA tissue that affects the response of host cells to HPV infection, and provides cues for the further exploration of the roles of m6A and the development of a novel treatment strategy for CA.

Alternative splicing analysis of lignocellulose-degrading enzyme genes and enzyme variants in Aspergillus niger

Alternative splicing (AS) greatly expands the protein diversity in eukaryotes. Although AS variants have been frequently reported existing in filamentous fungi, it remains unclear whether lignocellulose-degrading enzyme genes in industrially important fungi undergo AS events. In this work, AS events of lignocellulose-degrading enzymes genes in Aspergillus niger under two carbon sources (glucose and wheat straw) were investigated by RNA-Seq. The results showed that a total of 23 out of the 56 lignocellulose-degrading enzyme genes had AS events and intron retention was the main type of these AS events. The AS variant enzymes from the annotated endo-β-1,4-xylanase F1 gene (xynF1) and the endo-β-1,4-glucanase D gene (eglD), noted as XYNF1-AS and EGLD-AS, were characterized compared to their normal splicing products XYNF1 and EGLD, respectively. The AS variant XYNF1-AS displayed xylanase activity whereas XYNF1 did not. As for EGLD-AS and EGLD, neither of them showed annotated endo-β-1,4-glucanase activity. Instead, both showed lytic polysaccharide monooxygenase (LPMO) activity with some differences in catalytic properties. Our work demonstrated that the AS variants in A. niger were good sources for discovering novel lignocellulose-degrading enzymes. KEY POINTS: • AS events were identified in the lignocellulose-degrading enzyme genes of A. niger. • New β-1,4-xylanase and LPMO derived from AS events were characterized.

RNA splicing regulator EIF3D regulates the tumor microenvironment through immunogene-related alternative splicing in head and neck squamous cell carcinoma

Study finds that eukaryotic translation initiation factor 3 subunit D (EIF3D) may play an important role in aberrant alternative splicing (AS) events in tumors. AS possesses a pivotal role in both tumour progression and the constitution of the tumour microenvironment (TME). Regrettably, our current understanding of AS remains circumscribed especially in the context of immunogene-related alternative splicing (IGAS) profiles within Head and Neck Squamous Cell Carcinoma (HNSC). In this study, we comprehensively analyzed the function and mechanism of action of EIF3D by bioinformatics analysis combined with in vitro cellular experiments, and found that high expression of EIF3D in HNSC was associated with poor prognosis of overall survival (OS) and progression-free survival (PFS). The EIF3D low expression group had a higher degree of immune infiltration and better efficacy against PD1 and CTLA4 immunotherapy compared to the EIF3D high expression group. TCGA SpliceSeq analysis illustrated that EIF3D influenced differentially spliced alternative splicing (DSAS) events involving 105 differentially expressed immunogenes (DEIGs). We observed an induction of apoptosis and a suppression of cell proliferation, migration, and invasion in EIF3D knock-down FaDu cells. RNA-seq analysis unveiled that 531 genes exhibited differential expression following EIF3D knockdown in FaDu cells. These include 52 DEIGs. Furthermore, EIF3D knockdown influenced the patterns of 1923 alternative splicing events (ASEs), encompassing 129 IGASs. This study identified an RNA splicing regulator and revealed its regulatory role in IGAS and the TME of HNSC, suggesting that EIF3D may be a potential target for predicting HNSC prognosis and immunotherapeutic response.

RNA-binding proteins potentially regulate the alternative splicing of cell cycle-associated genes in proliferative diabetic retinopathy

RNA-binding proteins (RBPs) contribute to the pathogenesis of proliferative diabetic retinopathy (PDR) by regulating gene expression through alternative splicing events (ASEs). However, the RBPs differentially expressed in PDR and the underlying mechanisms remain unclear. Thus, this study aimed to identify the differentially expressed genes in the neovascular membranes (NVM) and retinas of patients with PDR. The public transcriptome dataset GSE102485 was downloaded from the Gene Expression Omnibus database, and samples of PDR and normal retinas were analyzed. A mouse model of oxygen-induced retinopathy was used to confirm the results. The top 20 RBPs were screened for co-expression with alternative splicing genes (ASGs). A total of 403 RBPs were abnormally expressed in the NVM and retina samples. Functional analysis demonstrated that the ASGs were enriched in cell cycle pathways. Cell cycle-associated ASEs and an RBP-AS regulatory network, including 15 RBPs and their regulated ASGs, were extracted. Splicing factor proline/glutamine rich (SFPQ), microtubule-associated protein 1 B (MAP1B), heat-shock protein 90-alpha (HSP90AA1), microtubule-actin crosslinking factor 1 (MACF1), and CyclinH (CCNH) expression remarkably differed in the mouse model. This study provides novel insights into the RBP-AS interaction network in PDR and for developing screening and treatment options to prevent diabetic retinopathy-related blindness.

Fascin actin-bundling protein 1 regulates non-small cell lung cancer progression by influencing the transcription and splicing of tumorigenesis-related genes

Background

High mortality rates are prevalent among patients with non-small-cell lung cancer (NSCLC), and effective therapeutic targets are key prognostic factors. Fascin actin-bundling protein 1 (FSCN1) promotes NSCLC; however, its role as an RNA-binding protein in NSCLC remains unexplored. Therefore, we aimed to explore FSCN1 expression and function in A549 cells.

Method

We screened for alternative-splicing events and differentially expressed genes (DEGs) after FSCN1 silence via RNA-sequencing (RNA-seq). FSCN1 immunoprecipitation followed by RNA-seq were used to identify target genes whose mRNA expression and pre-mRNA alternative-splicing levels might be influenced by FSCN1.

Results

Silencing FSCN1 in A549 cells affected malignant phenotypes; it inhibited proliferation, migration, and invasion, and promoted apoptosis. RNA-seq analysis revealed 2,851 DEGs and 3,057 alternatively spliced genes. Gene ontology-based functional enrichment analysis showed that downregulated DEGs and alternatively splicing genes were enriched for the cell-cycle. FSCN1 promoted the alternative splicing of cell-cycle-related mRNAs involved in tumorigenesis (i.e., BCCIP, DLGAP5, PRC1, RECQL5, WTAP, and SGO1). Combined analysis of FSCN1 RNA-binding targets and RNA-seq data suggested that FSCN1 might affect ACTG1, KRT7, and PDE3A expression by modulating the pre-mRNA alternative-splicing levels of NME4, NCOR2, and EEF1D, that were bound to long non-coding RNA transcripts (RNASNHG20, NEAT1, NSD2, and FTH1), which were highly abundant. Overall, extensive transcriptome analysis of gene alternative splicing and expression levels was performed in cells transfected with FSCN1 short-interfering RNA. Our data provide global insights into the regulatory mechanisms associated with the roles of FSCN1 and its target genes in lung cancer.

Genome-wide analysis revealed the dysregulation of RNA binding protein-correlated alternative splicing events in myocardial ischemia reperfusion injury

BACKGROUND

Myocardial ischemia reperfusion injury (MIRI), the tissue damage which is caused by the returning of blood supply to tissue after a period of ischemia, greatly reduces the therapeutic effect of treatment of myocardial infarction. But the underlying functional mechanisms of MIRI are still unclear.

METHODS

We constructed mouse models of MIRI, extracted injured and healthy myocardial tissues, and performed transcriptome sequencing experiments (RNA-seq) to systematically investigate the dysregulated transcriptome of MIRI, especially the alternative splicing (AS) regulation and RNA binding proteins (RBPs). Selected RBPs and MIRI-associated AS events were then validated by RT-qPCR experiments.

RESULTS

The differentially expressed gene (DEG) analyses indicated that transcriptome profiles were changed by MIRI and that DEGs' enriched functions were consistent with MIRI's dysregulated pathways. Furthermore, the AS profile was synergistically regulated and showed clear differences between the mouse model and the healthy samples. The exon skipping events significantly increased in MIRI model samples, while the opposite cassette exon events significantly decreased. According to the functional analysis, regulated alternative splicing genes (RASGs) were enriched in protein transport, cell division /cell cycle, RNA splicing, and endocytosis pathways, which were associated with the development of MIRI. Meanwhile, 493 differentially expressed RBPs (DE RBPs) were detected, most of which were correlated with the changed ratios of AS events. In addition, nine DE RBP genes were validated, including Eif5, Pdia6, Tagln2, Vasp, Zfp36l2, Grsf1, Idh2, Ndrg2, and Uqcrc1. These nine DE RBPs were correlated with RASGs enriched in translation process, cell growth and division, and endocytosis pathways, highly consistent with the functions of all RASGs. Finally, we validated the AS ratio changes of five regulated alternative splicing events (RASEs) derived from important regulatory genes, including Mtmr3, Cdc42, Cd47, Fbln2, Vegfa, and Fhl2.

CONCLUSION

Our study emphasized the critical roles of the dysregulated AS profiles in MIRI development, investigated the potential functions of MIRI-associated RASGs, and identified regulatory RBPs involved in AS regulation. We propose that the identified RASEs and RBPs could serve as important regulators and potential therapeutic targets in MIRI treatment in the future.

FUS regulates the alternative splicing of cell proliferation genes related to atherosclerosis

FUS plays a significant role as an RNA-binding protein in several cellular processes, including RNA splicing, DNA repair, and transcriptional regulation. However, the RNA-binding capacity of FUS in atherosclerosis is unclear. We aimed to study the functions of FUS in inflammatory regulation through the role of the splicing factor. We knocked down FUS with siRNA to further study the overall transcriptional level and select alternative splicing (AS) of FUS regulation in human umbilical vein endothelial cells (HUVECs) by RNA sequencing. The results suggested that the knockdown of FUS significantly affected gene expression in HUVECs. In addition, the knockdown of FUS resulted in 200 differentially expressed genes (DEGs) that were highly related to apoptotic process, signal transduction, multicellular organism development, cell adhesion and regulation of transcription, and DNA-templated pathways. Importantly, FUS extensively regulated 2870 AS events with a significant difference. Functional analysis of its modulated AS genes revealed they were highly enriched in cell cycle and cell population proliferation pathways. The qRT-PCR and RNA-seq data showed consistent results. Our findings suggested new knowledge of the mechanisms of FUS associated with atherosclerosis.

ROCK2-RNA interaction map reveals multiple biological mechanisms underlying tumor progression in renal cell carcinoma

Renal cell carcinoma (RCC) is the most common form of kidney cancer in adults. Despite new therapeutic modalities, the outcomes for RCC patients remain unsatisfactory. Rho-associated coiled-coil forming protein kinase 2 (ROCK2) has previously been shown to be upregulated in RCC, and its expression was negatively correlated with patient survival. However, the precise molecular function of ROCK2 has remained unclear. Herein, using RNA-seq analysis of ROCK2 knockdown and control cells, we identified 464 differentially expressed genes, and 1287 alternative splicing events in 786-O RCC cells. Furthermore, mapping of iRIP-seq reads in 786-O cells showed a biased distribution at 5' UTR, intronic and intergenic regions. By comparing ROCK2-regulated alternative splicing and iRIP-seq data, we found 292 overlapping genes that are enriched in multiple tumorigenic pathways. Taken together, our work defined a complex ROCK2-RNA interaction map on a genomic scale in a human RCC cell line, which deepens our understanding of the molecular function of ROCK2 in cancer development.

Genome-wide analysis of RNA-binding protein co-expression with alternative splicing events in acute respiratory distress syndrome following hematopoietic stem cell transplantation

Patients undergoing allogeneic hematopoietic stem cell transplantation (HSCT) are at an increased risk of developing severe acute respiratory distress syndrome (ARDS), which is characterized by peripheral bilateral patchy lung involvement. The regulatory network of RNA-binding protein (RBP)-alternative splicing (AS) in ARDS following HSCT has not been investigated. We hypothesize that RBP-AS plays a regulatory role during HSCT-ARDS. The published ARDS transcriptome data after HSCT (GSE84439) were downloaded, and the transcriptome data of 13 mRNAs were obtained by sequencing the peripheral blood of 5 HSCT-ARDS patients and 8 ARDS patients through high-throughput sequencing technology. Systematic analysis of downloaded data was performed to obtain differentially expressed RBPs, and the differentially alternative spliced pre-mRNAs in HSCT-ARDS and control groups were used to explore the global gene RBP-AS regulatory network. A total of 1769 differentially expressed genes and 4714 regulated alternative splicing events were identified in peripheral blood from HSCT-ARDS, of which 254 genes had both differential expression and differential AS. In addition, 128 RBPs were identified, of which HDGF, PCBP2, RIOK3, CISD2, and TRIM21, DDX58, MOV10 showed significantly increased or decreased expression in the HSCT-ARDS. RBPs with decreased expression had antiviral activity, while those with increased expression were involved in ROS, fibrosis, and negative viral resistance. The RBP-RASE-RASG regulatory network is constructed. It is related to the dysregulation of antiviral immunomodulation, imbalance in ROS homeostasis and pro-pulmonary fibrosis, which are involved in the development of HSCT-ARDS.

TFRC-RNA interactions show the regulation of gene expression and alternative splicing associated with IgAN in human renal tubule mesangial cells

Introduction: IgA nephropathy (IgAN) is the most common primary glomerular disease (PGD) which could progress to renal failure and is characterized by aberrant IgA immune complex deposition. Transferrin receptor1 (TFRC), an IgA receptor, is a potential RNA binding protein (RBP) which regulates expression of genes positively associated with the cell cycle and proliferation and is involved in IgAN. Molecular mechanisms by which TFRC affects IgAN development remain unclear. Methods: In this study, TFRC was overexpressed in human renal tubular mesangial cells (HRMCs) and RNA-sequencing (RNA-seq) and improved RNA immunoprecipitation sequencing (iRIP-seq) were performed. The aim was to identify potential RNA targets of TFRC at transcriptional and alternative splicing (AS) levels. Results: TFRC-regulated AS genes were enriched in mRNA splicing and DNA repair, consistent with global changes due to TFRC overexpression (TFRC-OE). Expression of TFRC-regulated genes potentially associated with IgAN, including CENPH, FOXM1, KIFC1, TOP2A, FABP4, ID1, KIF20A, ATF3, H19, IRF7, and H1-2, and with AS, CYGB, MCM7 and HNRNPH1, were investigated by RT-qPCR and iRIP-seq data analyzed to identify TFRC-bound RNA targets. RCC1 and RPPH1 were found to be TFRC-bound RNA targets involved in cell proliferation. Discussion: In conclusion, molecular TFRC targets were identified in HRMCs and TFRC found to regulate gene transcription and AS. TFRC is considered to have potential as a clinical therapeutic target.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

Transcriptome-wide analysis reveals the coregulation of RNA-binding proteins and alternative splicing genes in the development of atherosclerosis

RNA-binding proteins (RBPs) are involved in the regulation of RNA splicing, stability, and localization. How RBPs control the development of atherosclerosis, is not fully understood. To explore the relevant RNA-binding proteins (RBPs) and alternative splicing events (ASEs) in atherosclerosis. We made a comprehensive work to integrate analyses of differentially expressed genes, including differential RBPs, and variable splicing characteristics related to different stages of atherosclerosis in dataset GSE104140. A total of 3712 differentially expressed genes (DEGs) were identified, including 2921 upregulated genes and 791 downregulated genes. Further analysis screened out 54 RBP genes, and 434 AS genes overlapped DEGs. We selected high expression ten RBP genes (SAMHD1, DDX60 L, TLR7, RBM47, MYEF2, RNASE6, PARP12, APOBEC3G, SMAD9, and RNASE1) for co-expression analysis. Meanwhile, we found seven regulated alternative splicing genes (RASGs) (ABI1, FXR1, CHID1, PLEC, PRKACB, BNIP2, PPP3CB) that could be regulated by RBPs. The co-expression network was used to further elucidate the regulatory and interaction relationship between RBPs and AS genes. Apoptotic process and innate immune response, revealed by the functional enrichment analysis of RASGs regulated by RBPs were closely related to atherosclerosis. In addition, 26 of the 344 alternative splicing genes regulated by the above 10 RBPs were transcription factors (TFs), We selected high expression nine TFs (TFDP1, RBBP7, STAT2, CREB5, ERG, ELF1, HMGN3, BCLAF1, and ZEB2) for co-expression analysis. The target genes of these TFs were mainly enriched in inflammatory and immune response pathways that were associated with atherosclerosis. indicating that AS abnormalities of these TFs may have a function in atherosclerosis. Furthermore, the expression of differentially expressed RBPs and the alternative splicing events of AS genes was validated by qRT-PCR in umbilical vein endothelial cells (HUVEC). The results showed that RBM47 were remarkedly difference in HUVEC treated with ox-LDL and the splicing ratio of AS in BCLAF1which is regulated by RBM47 significantly changed. In conclusion, the differentially expressed RBPs identified in our analysis may play important roles in the development of atherosclerosis by regulating the AS of these TF genes.

Genome-wide analysis of RNA-binding proteins co-expression with alternative splicing events in mitral valve prolapse

Objectives

We investigated the role and molecular mechanisms of RNA-binding proteins (RBPs) and their regulated alternative splicing events (RASEs) in the pathogenesis of mitral valve prolapse (MVP).

Methods

For RNA extraction, we obtained peripheral blood mononuclear cells (PBMCs) from five patients with MVP, with or without chordae tendineae rupture, and five healthy individuals. High-throughput sequencing was used for RNA sequencing (RNA-seq). Differentially expressed genes (DEGs) analysis, alternative splicing (AS) analysis, functional enrichment analysis, co-expression of RBPs, and alternative splicing events (ASEs) analysis were conducted.

Results

The MVP patients exhibited 306 up-regulated genes and 198 down-regulated genes. All down- and up-regulated genes were enriched in both Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, MVP was closely associated with the top 10 enriched terms and pathways. In MVP patients, 2,288 RASEs were found to be significantly different, and four suitable RASEs (CARD11 A3ss, RBM5 ES, NCF1 A5SS, and DAXX A3ss) were tested. We identified 13 RNA-binding proteins (RBPs) from the DEGs and screened out four RBPs (ZFP36, HSPA1A, TRIM21, and P2RX7). We selected four RASEs based on the co-expression analyses of RBPs and RASEs, including exon skipping (ES) of DEDD2, alternative 3' splice site (A3SS) of ETV6, mutually exclusive 3'UTRs (3pMXE) of TNFAIP8L2, and A3SS of HLA-B. Furthermore, the selected four RBPs and four RASEs were validated by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and showed high consistency with RNA sequencing (RNA-seq).

Conclusion

Dysregulated RBPs and their associated RASEs may play regulatory roles in MVP development and may therefore be used as therapeutic targets in the future.

ANRIL overexpression globally induces expression and alternative splicing of genes involved in inflammation in HUVECs

Long non‑coding (lnc)RNAs serve important cellular functions and certain lncRNAs have roles in different mechanisms of gene regulation. lncRNA‑antisense non‑coding RNA in the INK4 locus (ANRIL) affects cell inflammation; however, the potential genes underlying the inflammatory response regulated by ANRIL remain unclear. In the present study, the potential function of ANRIL in regulating gene expression and alternative splicing was assessed. ANRIL‑regulated human umbilical vein endothelial cell (HUVEC) transcriptome was obtained using high‑throughput RNA sequencing (RNA‑seq) to evaluate the potential role of ANRIL. Following plasmid transfection, gene expression profile and alternative splicing pattern of HUVECs overexpressing ANRIL were analyzed using RNA‑seq. ANRIL overexpression affected the transcription levels of genes associated with the inflammatory response, NF‑κB signaling pathway, type I interferon‑mediated signal transduction pathway and innate immune response. ANRIL regulated the alternative splicing of hundreds of genes with functions such as gene expression, translation, DNA repair, RNA processing and participation in the NF‑κB signaling pathway. Many of these genes serve a key role in the inflammatory response. ANRIL‑regulated inflammatory response may be achieved by regulating alternate splicing and transcription. The present study broadened the understanding of ANRIL‑mediated gene regulation mechanisms and clarified the role of ANRIL in mediating inflammatory response mechanisms.

KIAA1429 regulates alternative splicing events of cancer-related genes in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) remains one of the most fatal malignancies with high morbidity and mortality rates in the world, whose molecular pathogenesis is incompletely understood. As an RNA-binding protein participating in the processing and modification of RNA, KIAA1429 has been proved to be implicated in the pathogenesis of multiple cancers. However, how KIAA1429 functions in alternative splicing is not fully reported. In the current study, multi-omics sequencing data were used to analyze and decipher the molecular functions and the underlying mechanisms of KIAA1429 in HCC samples. RNA sequencing data (RNA-seq) analysis demonstrated that in HCCLM3 cells, alternative splicing (AS) profiles were mediated by KIAA1429. Regulated AS genes (RASGs) by KIAA1429 were enriched in cell cycle and apoptosis-associated pathways. Furthermore, by integrating the RNA immunoprecipitation and sequencing data (RIP-seq) of KIAA1429, we found that KIAA1429-bound transcripts were highly overlapping with RASGs, indicating that KIAA1429 could globally regulate the alternative splicing perhaps by binding to their transcripts in HCCLM3 cells. The overlapping RASGs were also clustered in cell cycle and apoptosis-associated pathways. In particular, we validated the regulated AS events of three genes using clinical specimens from HCC patients, including the exon 6 of BPTF gene and a marker gene of HCC. In summary, our results shed light on the regulatory functions of KIAA1429 in the splicing process of pre-mRNA and provide theoretical basis for the targeted therapy of HCC.

Alternative Transcription Start Site Usage and Functional Implications in Pathogenic Fungi

Pathogenic fungi require delicate gene regulation mechanisms to adapt to diverse living environments and escape host immune systems. Recent advances in sequencing technology have exposed the complexity of the fungal genome, thus allowing the gradual disentanglement of multiple layers of gene expression control. Alternative transcription start site (aTSS) usage, previously reported to be prominent in mammals and to play important roles in physiopathology, is also present in fungi to fine-tune gene expression. Depending on the alteration in their sequences, RNA isoforms arising from aTSSs acquire different characteristics that significantly alter their stability and translational capacity as well as the properties and biologic functions of the resulting proteins. Disrupted control of aTSS usage has been reported to severely impair growth, virulence, and the infectious capacity of pathogenic fungi. Here, we discuss principle concepts, mechanisms, and the functional implication of aTSS usage in fungi.

DDX17 modulates the expression and alternative splicing of genes involved in apoptosis and proliferation in lung adenocarcinoma cells

Background

The DEAD-box RNA-binding protein (RBP) DDX17 has been found to be involved in the tumorigenesis of many types of cancers. However, the role of DDX17 in lung adenocarcinoma (LUAD) remains unclear.

Methods

We silenced DDX17 expression in A549 LUAD cells by small interfering RNA (siRNA). Cell proliferation and apoptosis assays were performed to explore the functions of DDX17. Knockdown of DDX17 by siRNA significantly inhibited proliferation and induced apoptosis in A549 cells. We used high-throughput RNA sequencing (RNA-seq) to identify differentially expressed genes (DEGs) and alternative splicing (AS) events in DDX17 knockdown LUAD cells.

Results

DDX17 knockdown increased the expression levels of proapoptotic genes and decreased those of proproliferative genes. Moreover, the DDX17-regulated AS events in A549 cells revealed by computational analysis using ABLas software were strongly validated by quantitative reverse transcription-polymerase chain reaction (RT-qPCR) and were also validated by analysis of The Cancer Genome Atlas (TCGA)-LUAD dataset. These findings suggest that DDX17 may function as an oncogene by regulating both the expression and AS of proliferation- and apoptosis-associated genes in LUAD cells. Our findings may offer new insights into understanding the molecular mechanisms of LUAD and provide a new therapeutic direction for LUAD.

PPIB-regulated alternative splicing of cell cycle genes contributes to the regulation of cell proliferation

BACKGROUND

Peptidylprolyl cis-trans isomerase B (PPIB) plays an important role in the process of inflammation through binding RNA, but the molecular pathogenesis is not yet clearly understood. The objective of this study was to investigate and verify the PPIB-regulated gene expressions and alternative splicing in Hela cells.

METHODS

We examined the PPIB-regulated transcriptomes in HeLa cells using RNA-seq data. Differentially expressed genes, alternative splicing analysis were carried out. Functional enrichment analysis was used to define the enrichment of each term.

RESULTS

We found that PPIB knockdown successfully downregulated PPIB in Hela cells, which promoted cell proliferation (P < 0.001), but no significant effect on cell apoptosis. Ten alternative splicing genes regulated by PPIB were detected in Hela cells. The ten top Gene Ontology biological process analysis and functional pathways of the alternative splicing genes were screened and identified. The pathways where differentially expressed genes are most enriched were toll-like receptor 4 signaling pathways and other signaling pathways relating to inflammation and immune response. In addition, PPIB affects the alternative splicing of multiple genes, and the Gene Ontology-biological process analysis showed that genes significantly related to alternative splicing changes were mainly enriched in the cell cycle (P < 0.05).

CONCLUSION

PPIB regulates the alternative splicing of cell cycle-related genes to affect cell proliferation and regulate the occurrence and development of chronic inflammatory diseases such as Nasal polyps.

Genome-wide identification of alternative splicing associated with histone deacetylase inhibitor in cutaneous T-cell lymphomas

Cutaneous T-cell lymphomas (CTCLs) are a kind of non-Hodgkin lymphoma that originates from skin, which is difficult to treat with traditional drugs. Human histone deacetylase inhibitors (HDACi) targeted therapy has become a promising treatment strategy in recent years, but some patients can develop resistance to the drug, leading to treatment failure. There are no public reports on whether alternative splicing (AS) and RNA binding proteins (RBP) affect the efficacy of targeted therapy. Using data from the Gene Expression Omnibus (GEO) database, we established a co-change network of AS events and RBP in CTCLs for the first time, and analyzed the potential regulatory effects of RBP on HDACi-related AS events. The dataset GSE132053, which contained the RNA sequence data for 17 HDACi samples, was downloaded and clean reads were aligned to the human GRCh38 genome by hierarchical indexing for spliced alignment of the transcripts, allowing four mismatches. Gene expression levels were evaluated using exons per million fragments mapped for each gene. Student’s t-tests were performed to evaluate the significance of changes in ratios for AS events, and regulated alternative splicing events (RASEs) were defined as events with p values less than 0.05. To sort the differentially expressed genes functional categories, Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways were identified using the KOBAS 2.0 server. The regulatory mechanisms of the RASEs and RBPs were evaluated using Pearson’s correlation coefficient. Seven indirect events of HDACi resistance or sensitivity were identified: NIR_5151_RP11-977G19.10, NIR_4557_IRAG2, NIR_11870_SUMO1, NIR_5347_ING4, NIR_17935_DNAJC2, NIR_17974_CBLL1, and NIR_422_SLC50A1. The potential regulatory relationships between RBPs and HDACi-sensitive RASEs were also analyzed. LEPR and HNRNPAO significantly affected NIR_11870_SUMO1, suggesting a potential regulatory relationship. Additionally, CNN1 may regulate NIR_5347_ING4, CNOT3 may regulate NIR_17935_DNAJC2, and DQX1 and LENG9 may regulate NIR_422_SLC5A1. Overall, our findings establish a theoretical foundation for the precise targeted treatment of CTCLs with HDACi.

Landscape and regulation of alternative splicing and alternative polyadenylation in a plant pathogenic fungus

Alternative splicing (AS) and alternative polyadenylation (APA) contribute significantly to the regulation of gene expression in higher eukaryotes. Their biological impact in filamentous fungi, however, is largely unknown. Here we combine PacBio Isoform-Sequencing and strand-specific RNA-sequencing of multiple tissues and mutant characterization to reveal the landscape and regulation of AS and APA in Fusarium graminearum. We generated a transcript annotation comprising 51 617 isoforms from 17 189 genes. In total, 4997 and 11 133 genes are alternatively spliced and polyadenylated, respectively. Majority of the AS events alter coding sequences. Unexpectedly, the AS transcripts containing premature-termination codons are not sensitive to nonsense-mediated messenger RNA decay. Unlike in yeasts and animals, distal APA sites have strong signals, but proximal APA isoforms are highly expressed in F. graminearum. The 3'-end processing factors FgRNA15, FgHRP1, and FgFIP1 play roles in promoting proximal APA site usage and intron splicing. A genome-wide increase in intron inclusion and distal APA site usage and downregulation of the spliceosomal and 3'-end processing factors were observed in older and quiescent tissues, indicating intron inclusion and 3'-untranslated region lengthening as novel mechanisms in regulating aging and dormancy in fungi. This study provides new insights into the complexity and regulation of AS and APA in filamentous fungi.

DEK modulates both expression and alternative splicing of cancer‑related genes

DEK is known to be a potential proto‑oncogene and is highly expressed in gastric cancer (GC); thus, DEK is considered to contribute to the malignant progression of GC. DEK is an RNA‑binding protein involved in transcription, DNA repair, and selection of splicing sites during mRNA processing; however, its precise function remains elusive due to the lack of clarification of the overall profiles of gene transcription and post‑transcriptional splicing that are regulated by DEK. We performed our original whole‑genomic RNA‑Seq data to analyze the global transcription and alternative splicing profiles in a human GC cell line by comparing DEK siRNA‑treated and control conditions, dissecting both differential gene expression and potential alternative splicing events regulated by DEK. The siRNA‑mediated knockdown of DEK in a GC cell line led to significant changes in gene expression of multiple cancer‑related genes including both oncogenes and tumor suppressors. Moreover, it was revealed that DEK regulated a number of alternative splicing in genes which were significantly enriched in various cancer‑related pathways including apoptosis and cell cycle processes. This study clarified for the first time that DEK has a regulatory effect on the alternative splicing, as well as on the expression, of numerous cancer‑related genes, which is consistent with the role of DEK as a possible oncogene. Our results further expand the importance and feasibility of DEK as a clinical therapeutic target for human malignancies including GC.

PCBP1 is associated with rheumatoid arthritis by affecting RNA products of genes involved in immune response in Th1 cells

Rheumatoid arthritis (RA) is an autoimmune disease characterized by persistent synovitis, in which T helper 1 (Th1) can promote the development of a pro-inflammatory microenvironment. Poly(rC)-binding protein 1 (PCBP1) has been identified as a promising biomarker of RA, while its molecular mechanisms in RA development are unknown. As a canonical RNA binding protein, we propose that PCBP1 could play roles in RA by affecting both expression and alternative splicing levels in Th1 cells. Here, microarray datasets (GSE15573 and GSE23561), including 102 peripheral blood mononuclear cell samples from 39 RA patients and 63 controls, were used to evaluate the PCBP1 expression changes in RA patients. High throughput sequencing data (GSE84702) of iron driven pathogenesis in Th1 cells were downloaded and reanalyzed, including two Pcbp1 deficiency samples and two control samples in Th1 cells. In addition, CLIP-seq data of PCBP1 in Jurkat T cells was also analyzed to investigate the regulatory mechanisms of PCBP1. We found PCBP1 were down-regulated in RA specimens compared with control. The result of differentially expressed genes (DEGs) showed that Pcbp1 silencing in Th1 cells affected the expression of genes involved in immune response pathway. Alternative splicing analysis also revealed that PCBP1-regulated alternative splicing genes (RASGs) were enriched in TNF-a/NF-κB signaling pathway, T cell activation, T cell differentiation and T cell differentiation associated immune response pathways, which were highly associated with RA. DEGs and RASGs by Pcbp1 deficiency in mice were validated in PBMCs specimens of RA patients by RT-qPCR. Investigation of the CLIP-seq data revealed PCBP1 preferred to bind to 3'UTR and intron regions. PCBP1-bound genes were also significantly associated with RASGs, identifying 102 overlapped genes of these two gene sets. These genes were significantly enriched in several immune response related pathways, including myeloid cell differentiation and positive regulation of NF-κB transcription factor activity. Two RA-related genes, PML and IRAK1, were screened from the above immune related pathways. These results together support our hypothesis that PCBP1 can regulate the expression of genes involved in immune response pathway, and can bind to and regulate the alternative splicing of immune response related genes in immune T cells, and ultimately participate in the molecular mechanism of RA, providing new research ideas and directions for clinical diagnosis and treatment.

Effect of Monoacylglycerol Lipase Inhibition on Intestinal Permeability of Rats With Severe Acute Pancreatitis

Background: Endocannabinoid 2-arachidonoylglycerol (2-AG) is an anti-nociceptive lipid that is inactivated through cellular uptake and subsequent catabolism by monoacylglycerol lipase (MAGL). In this study, we investigated the effects of MAGL inhibition on intestinal permeability and explored the possible mechanism.

Methods: A rat model of severe acute pancreatitis (SAP) was established. Rats were divided into three groups according to treatment. We analyzed intestinal permeability to fluorescein isothiocyanate-dextran and the levels of inflammatory factors interleukin-6 (IL-6) and tumor necrosis factor-α (TNF-α) and 2-AG. Hematoxylin and eosin staining was used to assess histological tissue changes. In vivo intestinal permeability was evaluated by transmission electron microscopy. We obtained ileum tissues, extracted total RNA, and applied RNA-sequencing. Sequencing data were analyzed by bioinformatics.

Results: Inflammatory factor levels were higher, while 2-AG levels were lower in the SAP group compared with the control group. Administration of JZL184 to rats with SAP increased the levels of 2-AG and lowered the levels of IL-6 and TNF-α. Notably, intestinal permeability was improved by JZL184 as demonstrated by fluorescein isothiocyanate-dextran measurement, hematoxylin and eosin staining, and transmission electron microscopy. RNA-sequencing showed significant transcriptional differences in SAP and JZL184 groups compared with the control group. KEGG analysis showed that the up- or downregulated genes in multiple comparison groups were enriched in two pathways, focal adhesion and PI3K-Akt signaling pathways. Differential alternative splicing (AS) genes, such as Myo9b, Lsp1, and Git2, have major functions in intestinal diseases. A total of 132 RNA-binding proteins (RBPs) were screened by crossing the identified abnormally expressed genes with the reported RBP genes. Among them, HNRNPDL coexpressed the most AS events as the main RBP.

Conclusion: MAGL inhibition improved intestinal mucosal barrier injury in SAP rats and induced a large number of differentially expressed genes and alternative splicing events. HNRNPDL might play an important role in improving intestinal mucosal barrier injury by affecting alternative splicing events.

Effects of alternative splicing events and transcriptome changes on kidney stone formation

During the development of urinary stone disease, the formation of tiny crystals that adhere to the renal tubular epithelium induces epithelial cell damage. This damage and repair of the epithelium is associated with the establishment of more crystal adhesion sites, which in turn stimulates further crystal adhesion and, eventually, stone formation. Deposited crystals typically cause changes in epithelial cell gene expression, such as transcriptome changes and alternative splicing events. Although considered important for regulating gene expression, alternative splicing has not been reported in studies related to kidney stones. To date, whether alternative splicing events are involved in the regulation of stone formation and whether crystallographic cell interactions are regulated by alternative splicing at the transcriptional level have remained unknown. Therefore, we conducted RNA sequencing and alternative splicing-related bioassays by modeling the in vitro stone environment. Many alternative splicing events were associated with crystallographic cell interactions. Moreover, these events regulated transcription and significantly affected the capacity of crystals to adhere to renal tubular epithelial cells and regulate apoptosis.

RALY regulate the proliferation and expression of immune/inflammatory response genes via alternative splicing of FOS

RALY is a multifunctional RNA-binding protein involved in cancer metastasis, prognosis, and chemotherapy resistance in various cancers. However, the molecular mechanism of which is still unclear. We have established RALY overexpression cell lines and studied the effect of RALY on proliferation and apoptosis in HeLa cells. Then we used RNA-seq to analyze the transcriptomes data. Lastly, RT-qPCR experiments had performed to confirm the RNA-seq results. We found that the overexpression of RALY in HeLa cells inhibited proliferation. Moreover, the overexpression of RALY changed the gene expression profile, and the significant upregulation of genes involved immune/inflammatory response related biological process by NOD-like receptor signaling pathway cytokine-cytokine receptor interaction. The significant downregulation genes involved innate immune response by the Primary immunodeficiency pathway. Notably, IFIT1, IFIT2, IFTI3, IFI44, HERC4, and OASL expression had inhibited by the overexpression of RALY. Furthermore, RALY negatively regulates the expression of transcription factors FOS and FOSB. Notably, we found that 645 alternative splicing events had regulated by overexpression of RALY, which is highly enriched in transcription regulation, RNA splicing, and cell proliferation biological process by the metabolic pathway. We show that RALY regulates the expression of immune/inflammatory response-related genes via alternative splicing of FOS in HeLa cells. The novel role of RALY in regulating immune/inflammatory gene expression may explain its function in regulating chemotherapy resistance and provides novel insights into further exploring the molecular mechanism of RALY in regulating cancer immunity and chemo/immune therapies.

RNA-Binding Proteins and Alternative Splicing Genes Are Coregulated in Human Retinal Endothelial Cells Treated with High Glucose

To explore the relevant RNA-binding proteins (RBPs) and alternative splicing events (ASEs) in diabetic retinopathy (DR). We devised a comprehensive work to integrate analyses of the differentially expressed genes, including differential RBPs, and variable splicing characteristics related to DR in human retinal endothelial cells induced by low glucose and high glucose in dataset GSE117238. A total of 2320 differentially expressed genes (DEGs) were identified, including 1228 upregulated genes and 1092 downregulated genes. Further analysis screened out 232 RBP genes, and 42 AS genes overlapped DEGs. We selected high expression and consistency six RBP genes (FUS, HNRNPA2B1, CANX, EIF1, CALR, and POLR2A) for coexpression analysis. Through analysis, we found eight RASGs (MDM2, GOLGA2P7, NFE2L1, KDM4A, FAM111A, CIRBP, IDH1, and MCM7) that could be regulated by RBP. The coexpression network was conducted to further elucidate the regulatory and interaction relationship between RBPs and AS. Apoptotic progress, protein phosphorylation, and NF-kappaB cascade revealed by the functional enrichment analysis of RASGs regulated by RBPs were closely related to diabetic retinopathy. Furthermore, the expression of differentially expressed RBPs was validated by qRT-PCR in mouse retinal microvascular endothelial cells and retinas from the streptozotocin mouse model. The results showed that Fus, Hnrnpa2b1, Canx, Calr, and Polr2a were remarkedly difference in high-glucose-treated retinal microvascular endothelial cells and Fus, Hnrnpa2b1, Canx, and Calr were remarkedly difference in retinas from streptozotocin-induced diabetic mice compared to control. The regulatory network between identified RBPs and RASGs suggests the presence of several signaling pathways possibly involved in the pathogenesis of DR. The verified RBPs should be further addressed by future studies investigating associations between RBPs and the downstream of AS, as they could serve as potential biomarkers and targets for DR.

PCBP1 regulates the transcription and alternative splicing of metastasis‑related genes and pathways in hepatocellular carcinoma

PCBP1 is a multifunctional RNA-binding protein (RBP) expressed in most human cells and is involved in posttranscriptional gene regulation. PCBP1 regulates the alternative splicing, translation and RNA stability of many cancer-related genes and has been identified as a potential tumour suppressor gene. PCBP1 inhibits the invasion of hepatocellular carcinoma (HCC) cells, but there are few studies on the specific regulatory target and mechanism of RBPs in HCC, and it is unclear whether PCBP1 plays a role in tumour metastasis as a splicing factor. We analysed the regulation of gene expression by PCBP1 at the transcriptional level. We obtained and analysed PCBP1-knockdown RNA-seq data and eCLIP-seq data of PCBP1 in HepG2 cells and found that PCBP1 widely regulates the alternative splicing and expression of genes enriched in cancer-related pathways, including extracellular matrix, cell adhesion, small molecule metabolic process and apoptosis. We validated five regulated alternative splicing events affected by PCBP1 using RT-qPCR and found that there was a significant difference in the expression of APOC1 and SPHK1 between tumour and normal tissues. In this study, we provided convincing evidence that human PCBP1 profoundly regulates the splicing of genes associated with tumour metastasis. These findings provide new insight into potential markers or therapeutic targets for HCC treatment.

Intron distribution and emerging role of alternative splicing in fungi

Spliceosomal introns are noncoding sequences that are spliced from pre-mRNA. They are ubiquitous in eukaryotic genomes, although the average number of introns per gene varies considerably between different eukaryotic species. Fungi are diverse in terms of intron numbers ranging from 4% to 99% genes with introns. Alternative splicing is one of the most common modes of posttranscriptional regulation in eukaryotes, giving rise to multiple transcripts from a single pre-mRNA and is widespread in metazoans and drives extensive proteome diversity. Earlier, alternative splicing was considered to be rare in fungi, but recently, increasing numbers of studies have revealed that alternative splicing is also widespread in fungi and has been implicated in the regulation of fungal growth and development, protein localization and the improvement of survivability, likely underlying their unique capacity to adapt to changing environmental conditions. However, the role of alternative splicing in pathogenicity and development of drug resistance is only recently gaining attention. In this review, we describe the intronic landscape in fungi. We also present in detail the newly discovered functions of alternative splicing in various cellular processes and outline areas particularly in pathogenesis and clinical drug resistance for future studies that could lead to the development of much needed new therapeutics.

Relationship between HSPA1A-regulated gene expression and alternative splicing in mouse cardiomyocytes and cardiac hypertrophy

Background

Cardiac hypertrophy may be classified as either physiological or pathological. Pathological hypertrophy has a complex etiology and is genetically regulated. In this study, we used a mouse model of cardiac hypertrophy to explore the mechanisms of gene regulation, in particular, modulation of the expression of target genes through transcription factor activity, regulation of immune and inflammation-associated genes and regulation of the alternative splicing of transcription factors.

Methods

Mouse models of pathological cardiac hypertrophy were established by transverse aortic constriction (TAC). We overexpressed HSPA1A in mouse cardiac HL-1 cells. GO and KEGG pathway annotation database was used to analyze all DEGs.

Results

The expression of HSPA1A differed significantly between TAC + dantrolene vs. sham + dantrolene (Sham was the non-TAC group, and DMSO was the contrast agent), and TAC + DMSO vs. sham + DMSO. The RNA-binding protein Zfp36 was found to be differentially expressed between both TAC + dantrolene vs. sham + dantrolene and TAC + DMSO vs. sham + DMSO. The expression of mki67 and gm5619 was significantly different between TAC + dantrolene and TAC + DMSO. HSPA1A was found to selectively regulate the expression of non-coding RNAs related to cardiac hypertrophy, including Rn7sk and RMRP. The downregulated genes were mainly related to inflammation and the immune response. HSPA1A negatively regulated alternative splicing of Asxl2 and positively regulated alternative splicing of Runx1.

Conclusions

HSPA1A was closely related to cardiac hypertrophy. Zfp36 was also related to cardiac hypertrophy. Dantrolene may delay cardiac hypertrophy and ventricular remodeling by regulating the expression of the RNA-binding protein genes mki67 and gm5619. HSPA1A positively regulated the expression of the non-coding RNAs RN7SK and RMRP while negatively regulating the expression of inflammation- and immune response-related genes. HSPA1A can play a role in cardiac hypertrophy by regulating the alternative splicing of asxl2 and runx1.

QKI-Regulated Alternative Splicing Events in Cervical Cancer: Pivotal Mechanism and Potential Therapeutic Strategy

QKI is a vital regulator in RNA splicing and maturation, but its role in cervical cancer (CC) is little known. In this study, we found that QKI is decreased in human CC, and overexpression of QKI inhibits HeLa cell proliferation and promotes the apoptosis of cancer cells. We identified hundreds of endogenous QKI-regulated alternative splicing events (ASEs) and differentially expressed genes (DEGs) in QKI-overexpressed HeLa cells by RNA-seq and selectively validated their expression by quantitative reverse-transcription polymerase chain reaction. The gene ontology and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that QKI-regulated ASEs and DEGs were closely related to cancer, apoptosis, and transcriptional regulatory functions. In short, QKI may affect the occurrence and development of CC by regulating gene expression through AS.

RNA-seq reveal RNA binding protein GNL3 as a key mediator in the development of psoriasis vulgaris by regulating the IL23/IL17 axis

BACKGROUND

Psoriasis is a systemic chronic inflammatory skin disorder that was prone to recurrence. The RNA binding protein GNL3 has an important function in maintaining the proliferative ability of stem cells, and its overexpression leads to apoptosis. GNL3 is expressed in the epidermis, however, its regulatory mechanism in psoriasis vulgaris is still poorly understood.

OBJECTIVE

To identify the role of GNL3 in the pathogenesis of psoriasis vulgaris.

MATERIALS AND METHODS

RNA-seq was performed to obtain the data of genes' expression and splicing events in Hela cells after shGNL3 and shCtrl was transferred. High quality results of differentially expressed genes (DEGs) and alternative splicing events (ASEs) were further attained by quality control and analysis. Through the functional enrichment analysis of DEGs and ASEs, the regulating effect of GNL3 was discussed, and the hypothesis was further confirmed in HaCat cells and psoriasis lesions.

RESULTS

The mRNA expression of IL23A in Hela cells was upregulated in GNL3 knockdown, and the ratio of ASE occurred in TNFAIP3 was increased. However, in HaCaT cells, the mRNA expression level of IL23A was downregulated in GNL3 knockdown, and the ratio of ASE of TNFAIP3 was decreased. Additionally, the results obtained in HaCaT cells was further validated in the lesional psoriatic skin.

CONCLUSION

GNL3 takes an important part in the development of psoriasis vulgaris by regulating the IL23/IL17 axis, which may serve as the basis of effective targeted treatment in future.

SUVA: splicing site usage variation analysis from RNA-seq data reveals highly conserved complex splicing biomarkers in liver cancer

Most of the current alternative splicing (AS) analysis tools are powerless to analyse complex splicing. To address this, we developed SUVA (Splice sites Usage Variation Analysis) that decomposes complex splicing events into five types of splice junction pairs. By analysing real and simulated data, SUVA showed higher sensitivity and accuracy in detecting AS events than the compared methods. Notably, SUVA detected extensive complex AS events and screened out 69 highly conserved and dominant AS events associated with cancer. The cancer-associated complex AS events in FN1 and the co-regulated RNA-binding proteins were significantly correlated with patient survival.

Evidence of a New MoYpd1p Phosphotransferase Isoform in the Multistep Phosphorelay System of Magnaporthe oryzae

Different external stimuli are perceived by multiple sensor histidine kinases and transmitted by phosphorylation via the phosphotransfer protein Ypd1p in the multistep phosphorelay system of the high osmolarity glycerol signaling pathway of filamentous fungi. How the signal propagation takes place is still not known in detail since multiple sensor histidine kinase genes in most filamentous fungi are coded in the genome, whereas only one gene for Ypd1p exists. That raises the hypothesis that various Ypd1p isoforms are produced from a single gene sequence, perhaps by alternative splicing, facilitating a higher variability in signal transduction. We found that the mRNA of MoYPD1 in the rice blast fungus Magnaporthe oryzae is subjected to an increased structural variation and amplified putative isoforms on a cDNA level. We then generated mutant strains overexpressing these isoforms, purified the products, and present here one previously unknown MoYpd1p isoform on a proteome level. Alternative splicing was found to be a valid molecular mechanism to increase the signal diversity in eukaryotic multistep phosphorelay systems.

Genome-wide alternative splicing profiling in the fungal plant pathogen Sclerotinia sclerotiorum during the colonization of diverse host families

Sclerotinia sclerotiorum is a notorious generalist plant pathogen that threatens more than 600 host plants, including wild and cultivated species. The molecular bases underlying the broad compatibility of S. sclerotiorum with its hosts is not fully elucidated. In contrast to higher plants and animals, alternative splicing (AS) is not well studied in plant-pathogenic fungi. AS is a common regulated cellular process that increases cell protein and RNA diversity. In this study, we annotated spliceosome genes in the genome of S. sclerotiorum and characterized their expression in vitro and during the colonization of six host species. Several spliceosome genes were differentially expressed in planta, suggesting that AS was altered during infection. Using stringent parameters, we identified 1,487 S. sclerotiorum genes differentially expressed in planta and exhibiting alternative transcripts. The most common AS events during the colonization of all plants were retained introns and the alternative 3' receiver site. We identified S. sclerotiorum genes expressed in planta for which (a) the relative accumulation of alternative transcripts varies according to the host being colonized and (b) alternative transcripts harbour distinct protein domains. This notably included 42 genes encoding predicted secreted proteins showing high-confidence AS events. This study indicates that AS events are taking place in the plant pathogenic fungus S. sclerotiorum during the colonization of host plants and could generate functional diversity in the repertoire of proteins secreted by S. sclerotiorum during infection.

CPEB3 regulates neuron-specific alternative splicing and involves neurogenesis gene expression

In the mammalian brain, alternative pre-mRNA splicing is a fundamental mechanism that modifies neuronal function dynamically where secretion of different splice variants regulates neurogenesis, development, pathfinding, maintenance, migration, and synaptogenesis. Sequence-specific RNA-Binding Protein CPEB3 has distinctive isoform-distinct biochemical interactions and neuronal development assembly roles. Nonetheless, the mechanisms moderating splice isoform options remain unclear. To establish the modulatory trend of CPEB3, we cloned and excessively expressed CPEB3 in HT22 cells. We used RNA-seq to analyze CPEB3-regulated alternative splicing on control and CPEB3-overexpressing cells. Consequently, we used iRIP-seq to identify CPEB-binding targets. We additionally validated CPEB3-modulated genes using RT-qPCR. CPEB3 overexpression had insignificant effects on gene expression in HT22 cells. Notably, CPEB3 partially modulated differential gene splicing enhanced in the modulation of neural development, neuron cycle, neurotrophin, synapse, and specific development pathway, implying an alternative splicing regulatory mechanism associated with neurogenesis. Moreover, qRT-PCR verified the CPEB3-modulated transcription of neurogenesis genes LCN2 and NAV2, synaptogenesis gene CYLD, as well as neural development gene JADE1. Herein, we established that CPEB3 is a critical modulator of alternative splicing in neurogenesis, which remarkably enhances the current understanding of the CPEB3 mediated alternative pre-mRNA splicing.

Hex1, the Major Component of Woronin Bodies, Is Required for Normal Development, Pathogenicity, and Stress Response in the Plant Pathogenic Fungus Verticillium dahliae

Woronin bodies are membrane-bound organelles of filamentous ascomycetes that mediate hyphal compartmentalization by plugging septal pores upon hyphal damage. Their major component is the peroxisomal protein Hex1, which has also been implicated in additional cellular processes in fungi. Here, we analyzed the Hex1 homolog of Verticillium dahliae, an important asexual plant pathogen, and we report its pleiotropic involvement in fungal growth, physiology, stress response, and pathogenicity. Alternative splicing of the Vdhex1 gene can lead to the production of two Hex1 isoforms, which are structurally similar to their Neurospora crassa homolog. We show that VdHex1 is targeted to the septum, consistently with its demonstrated function in sealing hyphal compartments to prevent excessive cytoplasmic bleeding upon injury. Furthermore, our investigation provides direct evidence for significant contributions of Hex1 in growth and morphogenesis, as well as in asexual reproduction capacity. We discovered that Hex1 is required both for normal responses to osmotic stress and factors that affect the cell wall and plasma-membrane integrity, and for normal resistance to oxidative stress and reactive oxygen species (ROS) homeostasis. The Vdhex1 mutant exhibited diminished ability to colonize and cause disease on eggplant. Overall, we show that Hex1 has fundamentally important multifaceted roles in the biology of V. dahliae.

Opportunities and Challenges in Studies of Host-Pathogen Interactions and Management of Verticillium dahliae in Tomatoes

Tomatoes (Solanum lycopersicum L.) are a valuable horticultural crop that are grown and consumed worldwide. Optimal production is hindered by several factors, among which Verticillium dahliae, the cause of Verticillium wilt, is considered a major biological constraint in temperate production regions. V. dahliae is difficult to mitigate because it is a vascular pathogen, has a broad host range and worldwide distribution, and can persist in soil for years. Understanding pathogen virulence and genetic diversity, host resistance, and plant-pathogen interactions could ultimately inform the development of integrated strategies to manage the disease. In recent years, considerable research has focused on providing new insights into these processes, as well as the development and integration of environment-friendly management approaches. Here, we discuss the current knowledge on the race and population structure of V. dahliae, including pathogenicity factors, host genes, proteins, enzymes involved in defense, and the emergent management strategies and future research directions for managing Verticillium wilt in tomatoes.

STAU1 selectively regulates the expression of inflammatory and immune response genes and alternative splicing of the nerve growth factor receptor signaling pathway

Double‑stranded RNA‑binding protein Staufen homolog 1 (STAU1) is a highly conserved multifunctional double‑stranded RNA‑binding protein, and is a key factor in neuronal differentiation. RNA sequencing was used to analyze the overall transcriptional levels of the upregulated cells by STAU1 and control cells, and select alternative splicing (AS). It was determined that the high expression of STAU1 led to changes in the expression levels of a variety of inflammatory and immune response genes, including IFIT2, IFIT3, OASL, and CCL2. Furthermore, STAU1 was revealed to exert a significant regulatory effect on the AS of genes related to the 'nerve growth factor receptor signaling pathway'. This is of significant importance for neuronal survival, differentiation, growth, post‑damage repair, and regeneration. In conclusion, overexpression of STAU1 was associated with immune response and regulated AS of pathways related to neuronal growth and repair. In the present study, the whole transcriptome of STAU1 expression was first analyzed, which laid a foundation for further understanding the key functions of STAU1.

SRSF6 regulates alternative splicing of genes involved in DNA damage response and DNA repair in HeLa cells

Alternative splicing (AS) occurs in nearly all human genes and abnormal AS has a close association with cancer. Serine and arginine-rich splicing factor 6 (SRSF6), a canonical member of the serine/arginine-rich protein family, has been characterized as an important regulator of AS. However, the role of SRSF6 in regulating AS in cancers has remained to be fully elucidated. In the present study, the median expression of SRSF6 in tumors was determined to be higher compared with that in matched normal tissues in 13 out of 16 cancer types from The Cancer Genome Atlas. To investigate the biological effects of SRSF6 overexpression, an SRSF6-overexpression model of HeLa cells was constructed and it was revealed that SRSF6 overexpression resulted in significantly higher apoptosis and lower proliferation compared to control cells. Transcriptome analysis indicated that overexpression of SRSF6 in cancer cells induced large-scale changes in transcriptional expression levels and AS. Two groups of cervical cancer tumor samples in which SRSF6 was differentially expressed were then selected to analyze potential SRSF6-regulated AS. It was determined that the pattern of SRSF6-regulated AS in clinical samples was similar to that in cancer cells and AS genes were enriched in DNA damage response (DDR) pathways, including DNA repair and double-strand break repair via homologous recombination. Furthermore, AS events regulated by SRSF6 were validated using reverse transcription-quantitative PCR. The present results highlighted that SRSF6 is able to trigger the activation of DDR pathways via regulation of AS to influence cancer progression. These results markedly expand the current understanding of the mechanisms underlying SRSF6-mediated gene regulation and suggest the potential use of SRSF6 as a therapeutic target in cancer.

IGF2BP3 May Contributes to Lung Tumorigenesis by Regulating the Alternative Splicing of PKM

RNA binding proteins (RBPs) play a key role in genome regulation. Here we report the post-transcript regulation of IGF2BP3, which belongs to the insulin-like growth factor 2 mRNA binding protein family. We used iRIP-seq and RNA-seq to analyze the transcript regulation and alternative splicing on IGF2BP3 treated with overexpression cells and control. Overexpressed IGF2BP3 has broadly increased genes expression which involved in G-protein coupled receptor signaling pathway, positive regulation of cell proliferation, and signal transduction. IGF2BP3 regulated alternative splicing of multiple genes mainly clustered at response to hypoxia, negative regulation of transcription, and embryonic development. This study first provides alternative splicing analysis on transcription level of IGF2BP3 regulation, which laid the foundation for later research on IGF2BP3 critical functions.

Filamin B extensively regulates transcription and alternative splicing, and is associated with apoptosis in HeLa cells

Post-transcriptional mechanisms are an important approach in the treatment of cancer, and may also be hijacked by tumor cells to help adapt to the local microenvironment. Filamin B (FLNB), an actin-binding protein that provides crucial scaffolds for cell motility and signaling, has also been identified as an RNA-binding protein. Recent studies demonstrated that FLNB might play an important role, not only in skeletal development, but also in regulating tumorigenesis; however, the effects of dysregulated expression of FLNB at the molecular level are not clear. In the present study, RNA-sequencing was performed to analyze changes in overall transcriptional and alternative splicing between the knocked-down FLNB and the control in HeLa cells. Decreased FLNB levels resulted in significantly lower apoptosis compared with control cells. FLNB knockdown extensively regulated the expression of genes in cell apoptosis, tumorigenesis, metastases, transmembrane transport and cartilage development. Moreover, FLNB regulated alternative splicing of a large number of genes involved in ‘cell death’ and the ‘apoptotic process’. Some genes and alternative splicing related to skeletal development were enriched and regulated by FLNB. Reverse transcription-quantitative-PCR identified FLNB-regulated transcription and alternative splicing of genes, such as NLR family apoptosis inhibitory protein, interleukin 23 subunit α, metastasis associated lung adenocarcinoma transcript 1, phosphofurin acidic cluster sorting protein 2, bone morphogenetic protein 7, matrix metallopeptidase 13, collagen type II α 1 chain, fibroblast growth factor receptor 2 and vitamin D receptor. The present study is the first study, to the best of the authors’ knowledge, to provide transcriptome-wide analysis of differential gene expression and alternative splicing upon FLNB silencing. The present results suggested that FLNB may play an important regulatory role in cervical cancer cell apoptosis via regulation of transcription and alternative splicing, which provide insight for the current understanding of the mechanisms of FLNB-mediated gene regulation.

Alternative Splicing in Heat Shock Protein Transcripts as a Mechanism of Cell Adaptation in Trichophyton rubrum

Heat shock proteins (HSPs) are involved in critical processes like host tissue invasion, resistance, and pathogenicity in dermatophytes. RNA-Seq analysis of Trichophyton rubrum exposed to undecanoic acid (UDA) revealed intron retention events in HSP transcripts. Because HSPs are modulated in response to various stimuli and as alternative splicing (AS) can result in a broad diversity in the proteome of eukaryotic cells, our objective was to confirm the aforementioned retention events, investigating their consequences and extent. Furthermore, we aimed to determine: (1) the expression profile of HSP genes in an infection-like scenario and (2) the importance of Hsp90 for the keratinolytic potential of T. rubrum. RT and qPCR analyses comparing the exposure to UDA and terbinafine (TRB) confirmed the presence of two mRNA isoforms of the hsp7-like gene, with distinct expression patterns in response to UDA and TRB. The HSP expression profile revealed two upregulated, three downregulated, and four unmodulated transcripts; Hsp90 inhibition by 17-AAG resulted in a significant decrease in keratinolytic potential at 37 °C. Altogether, these results broaden the current knowledge on the importance of HSP-mediated pathways for cell adaptation and other aspects of dermatophyte biology, indicating that HSP network proteins can be potential targets for antifungal therapy.

CRKL regulates alternative splicing of cancer-related genes in cervical cancer samples and HeLa cell

BACKGROUND

Aberrant spliced isoforms are specifically associated with cancer progression and metastasis. The cytoplasmic adaptor CRKL (v-crk avian sarcoma virus CT10 oncogene homolog-like) is a CRK like proto-oncogene, which encodes a SH2 and SH3 (src homology) domain-containing adaptor protein. CRKL is tightly linked to leukemia via its binding partners BCR-ABL and TEL-ABL, upregulated in multiple types of human cancers, and induce cancer cell proliferation and invasion. However, it remains unclear whether signaling adaptors such as CRKL could regulate alternative splicing.

METHODS

We analyzed the expression level of CRKL in 305 cervical cancer tissue samples available in TCGA database, and then selected two groups of cancer samples with CRKL differentially expressed to analyzed potential CRKL-regulated alternative splicing events (ASEs). CRKL was knocked down by shRNA to further study CRKL-regulated alternative splicing and the activity of SR protein kinases in HeLa cells using RNA-Seq and Western blot techniques. We validated 43 CRKL-regulated ASEs detected by RNA-seq in HeLa cells, using RT-qPCR analysis of HeLa cell samples and using RNA-seq data of the two group of clinical cervical samples.

RESULTS

The expression of CRKL was mostly up-regulated in stage I cervical cancer samples. Knock-down of CRKL led to a reduced cell proliferation. CRKL-regulated alternative splicing of a large number of genes were enriched in cancer-related functional pathways, among which DNA repair and G2/M mitotic cell cycle, GnRH signaling were shared among the top 10 enriched GO terms and KEGG pathways by results from clinical samples and HeLa cell model. We showed that CRKL-regulated ASEs revealed by computational analysis using ABLas software in HeLa cell were highly validated by RT-qPCR, and also validated by cervical cancer clinical samples.

CONCLUSIONS

This is the first report of CRKL-regulation of the alternative splicing of a number of genes critical in tumorigenesis and cancer progression, which is consistent with CRKL reported role as a signaling adaptor and a kinase. Our results underline that the signaling adaptor CRKL might integrate the external and intrinsic cellular signals and coordinate the dynamic activation of cellular signaling pathways including alternative splicing regulation.

The Rhizophagus irregularis Genome Encodes Two CTR Copper Transporters That Mediate Cu Import Into the Cytosol and a CTR-Like Protein Likely Involved in Copper Tolerance

Arbuscular mycorrhizal fungi increase fitness of their host plants under Cu deficient and toxic conditions. In this study, we have characterized two Cu transporters of the CTR family (RiCTR1 and RiCTR2) and a CTR-like protein (RiCTR3A) of Rhizophagus irregularis. Functional analyses in yeast revealed that RiCTR1 encodes a plasma membrane Cu transporter, RiCTR2 a vacuolar Cu transporter and RiCTR3A a plasma membrane protein involved in Cu tolerance. RiCTR1 was more highly expressed in the extraradical mycelia (ERM) and RiCTR2 in the intraradical mycelia (IRM). In the ERM, RiCTR1 expression was up-regulated by Cu deficiency and down-regulated by Cu toxicity. RiCTR2 expression increased only in the ERM grown under severe Cu-deficient conditions. These data suggest that RiCTR1 is involved in Cu uptake by the ERM and RiCTR2 in mobilization of vacuolar Cu stores. Cu deficiency decreased mycorrhizal colonization and arbuscule frequency, but increased RiCTR1 and RiCTR2 expression in the IRM, which suggest that the IRM has a high Cu demand. The two alternatively spliced products of RiCTR3, RiCTR3A and RiCTR3B, were more highly expressed in the ERM. Up-regulation of RiCTR3A by Cu toxicity and the yeast complementation assays suggest that RiCTR3A might function as a Cu receptor involved in Cu tolerance.