Gene regulatory evolution in cold-adapted fly populations neutralizes plasticity and may undermine genetic canalization

The relationships between adaptive evolution, phenotypic plasticity, and canalization remain incompletely understood. Theoretical and empirical studies have made conflicting arguments on whether adaptive evolution may enhance or oppose the plastic response. Gene regulatory traits offer excellent potential to study the relationship between plasticity and adaptation, and they can now be studied at the transcriptomic level. Here we take advantage of three closely-related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. We measure the transcriptome-wide plasticity in gene expression levels and alternative splicing (intron usage) between warm and cold laboratory environments. We find that suspected adaptive changes in both gene expression and alternative splicing tend to neutralize the ancestral plastic response. Further, we investigate the hypothesis that adaptive evolution can lead to decanalization of selected gene regulatory traits. We find strong evidence that suspected adaptive gene expression (but not splicing) changes in cold-adapted populations are more vulnerable to the genetic perturbation of inbreeding than putatively neutral changes. We find some evidence that these patterns may reflect a loss of genetic canalization accompanying adaptation, although other processes including hitchhiking recessive deleterious variants may contribute as well. Our findings augment our understanding of genetic and environmental effects on gene regulation in the context of adaptive evolution.

USP15 and USP4 facilitate lung cancer cell proliferation by regulating the alternative splicing of SRSF1

The deubiquitinating enzyme USP15 is implicated in several human cancers by regulating different cellular processes, including splicing regulation. However, the underlying molecular mechanisms of its functional relevance and the successive roles in enhanced tumorigenesis remain ambiguous. Here, we found that USP15 and its close paralog USP4 are overexpressed and facilitate lung cancer cell proliferation by regulating the alternative splicing of SRSF1. Depletion of USP15 and USP4 impair SRSF1 splicing characterized by the replacement of exon 4 with non-coding intron sequences retained at its C-terminus, resulting in an alternative isoform SRSF1-3. We observed an increased endogenous expression of SRSF1 in lung cancer cells as well, and its overexpression significantly enhanced cancer cell phenotype and rescued the depletion effect of USP15 and USP4. However, the alternatively spliced isoform SRSF1-3 was deficient in such aspects for its premature degradation through nonsense-mediated mRNA decay. The increased USP15 expression contributes to the lung adenocarcinoma (LUAD) development and shows significantly lower disease-specific survival of patients with USP15 alteration. In short, we identified USP15 and USP4 as key regulators of SRSF1 alternative splicing with altered functions, which may represent the novel prognostic biomarker as well as a potential target for LUAD.

N6-Methyladenosine Modification Changes in the Genome of Paulownia Fortunei Seedlings Infected With Phytoplasma

Background: Phytoplasmas induce diseases in more than 1,000 plant species and cause substantial ecological damage and economic losses, but the specific pathogenesis of phytoplasma has not yet been clarified. N6-methyladenosine sequencing (m6A-seq) has been applied mainly to model plants and not to woody plants. Results: In this study, we applied m6A-seq to study changes in m6A modification in the Paulownia fortunei genome after phytoplasma infection. We found that the m6A modification level in seedlings infected with the phytoplasma that causes Paulownia witches' broom (PaWB) was slightly higher than the m6A modification level in PaWB-infected seedlings treated with 60 mg·L−1 methyl methanesulfonate (MMS). MMS has been shown to restore PaWB-infected seedlings to their normal form and no phytoplasma can be detected in MMS-treated PaWB-infected seedlings. RNA sequencing (RNA-seq) and m6A-seq were used to analyze the expression of genes with m6A peaks and m6A motifs in genes, respectively. The correlation analysis between the RNA-seq and m6A-seq data detected that a total of 315 differentially methylated genes were predicted to be significantly differentially expressed at the transcriptome level. The functions of genes related to PaWB were predicted by functional enrichment analysis, and two genes related to maintenance of the basic mechanism of stem cells in shoot apical meristem were discovered. One of the genes encodes the receptor protein kinase CLV2 (Paulownia_LG2G000076), and the other gene encodes the homeobox transcription factor STM (Paulownia_LG15G000976). The m6A modification levels were higher in PaWB-infected seedlings than they were in MMS-treated seedlings. In addition, genes F-box (Paulownia_LG17G000760) and MSH5 (Paulownia_LG8G001160) had exon skipping and mutually exclusive exon types of alternative splicing in PaWB-infected seedling treated with MMS. RT-PCR verified that the alternative splicing of these two genes was related to m6A modification. Conclusions: In this study, we applied m6A-seq to determine methylation levels in phytoplasma-infected Paulownia, and combined m6A-seq with transcriptome analysis to screen differentially expressed genes associated with PaWB. Also analyzed the effect of m6A methylation on alternative splicing. In future studies, we plan to verify genes directly related to PaWB and methylation-related enzymes in Paulownia to elucidate the pathogenicity mechanism of PaWB caused by phytoplasma invasion.

Fine-Tuning of the Grain Size by Alternative Splicing of GS3 in Rice

Grain size is subtly regulated by multiple signaling pathways in rice. Alternative splicing is a general mechanism that regulates gene expression at the post-transcriptional level. However, to our knowledge, the molecular mechanism underlying grain size regulation by alternative splicing is largely unknown. GS3, the first identified QTL for grain size in rice, is regulated at the transcriptional and post-translational level. In this study, we identified that GS3 is subject to alternative splicing. GS3.1 and GS3.2, two dominant isoforms, accounts for about 50% and 40% of total transcripts, respectively. GS3.1 encodes the full-length protein, while GS3.2 generated a truncated proteins only containing OSR domain due to a 14 bp intronic sequence retention. Genetic analysis revealed that GS3.1 overexpressors decreased grain size, but GS3.2 showed no significant effect on grain size. Furthermore, we demonstrated that GS3.2 disrupts GS3.1 signaling by competitive occupation of RGB1. Therefore, we draw a conclusion that the alternative splicing of GS3 decreases the amount of GS3.1 and GS3.2 disrupts the GS3.1 signaling to inhibit the negative effects of GS3.1 to fine-tune grain size. Moreover, the mechanism is conserved in cereals rather than in Cruciferae, which is associated with its effects on grain size. The results provide a novel, conserved and important mechanism underlying grain size regulation at the post-transcriptional level in cereals.

Splice-switching of the insulin receptor pre-mRNA alleviates tumorigenic hallmarks in rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is an aggressive pediatric tumor with a poor prognosis for metastasis and recurrent disease. Large-scale sequencing endeavors demonstrate that Rhabdomyosarcomas have a dearth of precisely targetable driver mutations. However, IGF-2 signaling is known to be grossly altered in RMS. The insulin receptor (IR) exists in two alternatively spliced isoforms, IR-A and IR-B. The IGF-2 signaling molecule binds both its innate IGF-1 receptor as well as the insulin receptor variant A (IR-A) with high affinity. Mitogenic and proliferative signaling via the canonical IGF-2 pathway is, therefore, augmented by IR-A. This study shows that RMS patients express increased IR-A levels compared to control tissues that predominantly express the IR-B isoform. We also found that Hif-1α is significantly increased in RMS tumors, portraying their hypoxic phenotype. Concordantly, the alternative splicing of IR adapts to produce more IR-A in response to hypoxic stress. Upon examining the pre-mRNA structure of the gene, we identified a potential hypoxia-responsive element, which is also the binding site for the RNA-binding protein CUG-BP1 (CELF1). We designed Splice Switching Oligonucleotides (SSO) against this binding site to decrease IR-A levels in RMS cell lines and, consequently, rescue the IR-B expression levels. SSO treatment resulted in a significant reduction in cell proliferation, migration, and angiogenesis. Our data shows promising insight into how impeding the IGF-2 pathway by reducing IR-A expression mitigates tumor growth. It is evident that Rhabdomyosarcomas use IR alternative splicing as yet another survival strategy that can be exploited as a therapeutic intervention in conjunction with already established anti-IGF-1 receptor therapies.

Long read sequencing reveals novel isoforms and insights into splicing regulation during cell state changes

Background Alternative splicing is a key mechanism underlying cellular differentiation and a driver of complexity in mammalian neuronal tissues. However, understanding of which isoforms are differentially used or expressed and how this affects cellular differentiation remains unclear. Long read sequencing allows full-length transcript recovery and quantification, enabling transcript-level analysis of alternative splicing processes and how these change with cell state. Here, we utilise Oxford Nanopore Technologies sequencing to produce a custom annotation of a well-studied human neuroblastoma cell line SH-SY5Y, and to characterise isoform expression and usage across differentiation. Results We identify many previously unannotated features, including a novel transcript of the voltage-gated calcium channel subunit gene, CACNA2D2. We show differential expression and usage of transcripts during differentiation identifying candidates for future research into state change regulation. Conclusions Our work highlights the potential of long read sequencing to uncover previously unknown transcript diversity and mechanisms influencing alternative splicing.

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.

Rbfox1 is required for myofibril development and maintaining fiber type–specific isoform expression in Drosophila muscles

Protein isoform transitions confer muscle fibers with distinct properties and are regulated by differential transcription and alternative splicing. RNA-binding Fox protein 1 (Rbfox1) can affect both transcript levels and splicing, and is known to contribute to normal muscle development and physiology in vertebrates, although the detailed mechanisms remain obscure. In this study, we report that Rbfox1 contributes to the generation of adult muscle diversity in Drosophila. Rbfox1 is differentially expressed among muscle fiber types, and RNAi knockdown causes a hypercontraction phenotype that leads to behavioral and eclosion defects. Misregulation of fiber type–specific gene and splice isoform expression, notably loss of an indirect flight muscle–specific isoform of Troponin-I that is critical for regulating myosin activity, leads to structural defects. We further show that Rbfox1 directly binds the 3′-UTR of target transcripts, regulates the expression level of myogenic transcription factors myocyte enhancer factor 2 and Salm, and both modulates expression of and genetically interacts with the CELF family RNA-binding protein Bruno1 (Bru1). Rbfox1 and Bru1 co-regulate fiber type–specific alternative splicing of structural genes, indicating that regulatory interactions between FOX and CELF family RNA-binding proteins are conserved in fly muscle. Rbfox1 thus affects muscle development by regulating fiber type–specific splicing and expression dynamics of identity genes and structural proteins.

Comprehensive quantitative analysis of alternative splicing variants reveals the HNF1B mRNA splicing pattern in various tumour and non-tumour tissues

Hepatocyte nuclear factor-1-beta (HNF1B) is a transcription factor and putative biomarker of solid tumours. Recently, we have revealed a variety of HNF1B mRNA alternative splicing variants (ASVs) with unknown, but potentially regulatory, functions. The aim of our work was to quantify the most common variants and compare their expression in tumour and non-tumour tissues of the large intestine, prostate, and kidney. The HNF1B mRNA variants 3p, Δ7, Δ7–8, and Δ8 were expressed across all the analysed tissues in 28.2–33.5%, 1.5–2%, 0.8–1.7%, and 2.3–6.9% of overall HNF1B mRNA expression, respectively, and occurred individually or in combination. The quantitative changes of ASVs between tumour and non-tumour tissue were observed for the large intestine (3p, Δ7–8), prostate (3p), and kidney samples (Δ7). Decreased expression of the overall HNF1B mRNA in the large intestine and prostate cancer samples compared with the corresponding non-tumour samples was observed (p = 0.019 and p = 0.047, respectively). The decreased mRNA expression correlated with decreased protein expression in large intestine carcinomas (p < 0.001). The qualitative and quantitative pattern of the ASVs studied by droplet digital PCR was confirmed by next-generation sequencing, which suggests the significance of the NGS approach for further massive evaluation of the splicing patterns in a variety of genes.