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.

Establishing F1A-CreERT2 Mice to Trace Fgf1 Expression in Adult Mouse Cardiomyocytes

Fibroblast growth factor 1 (FGF1) regulates many biological and physiological processes. In mice, Fgf1 gene contains at least three upstream promoters and are alternatively spliced to the first protein coding exon, giving rise to different Fgf1 mRNA variants (1A, 1B and 1G). Among them, the Fgf1A transcript is predominantly expressed in the heart. FGF1 can induce cardiomyocyte regeneration and cardiogenesis in vitro and in vivo. Here, we generated a novel mouse line using the Fgf1A promoter (F1A) driving the expression of the inducible Cre recombinase (CreERT2). We firstly demonstrated that the highest mRNA expression of CreERT2 were detected in the heart specifically of F1A-CreERT2 mice, similar to that of Fgf1A mRNA. The F1A-CreERT2 mice were crossed with ROSA26 mice, and the F1 mice were analyzed. The LacZ-positive signals were detected exclusively in the heart after tamoxifen administration. The CreERT2-mediated recombination in the tissues is monitored through LacZ-positive signals, indicating the in situ localization of F1A-positive cells. Consistently, these F1A-positive cells with RFP-positive signals or LacZ-positive blue signals were co-localized with cardiomyocytes expressing cardiac troponin T, suggesting cardiomyocyte-specific activation of Fgf1A promoter. Our data suggested that the F1A-CreERT2 mouse line could be used for time-dependent and lineage tracing of Fgf1A-expressing cells in vivo.

mRNA Analysis of Frameshift Mutations with Stop Codon in the Last Exon: The Case of Hemoglobins Campania [α1 cod95 (−C)] and Sciacca [α1 cod109 (−C)]

An insertion or deletion of a nucleotide (nt) in the penultimate or the last exon can result in a frameshift and premature termination codon (PTC), giving rise to an unstable protein variant, showing a dominant phenotype. We described two α-globin mutants created by the deletion of a nucleotide in the penultimate or the last exon of the α1-globin gene: the Hb Campania or α1 cod95 (−C), causing a frameshift resulting in a PTC at codon 102, and the Hb Sciacca or α1 cod109 (−C), causing a frameshift and formation of a PTC at codon 133. The carriers showed α-thalassemia alterations (mild microcytosis with normal Hb A2) and lacked hemoglobin variants. The 3D model indicated the α-chain variants’ instability, due to the severe structural alterations with impairment of the chaperone alpha-hemoglobin stabilizing protein (AHSP) interaction. The qualitative and semiquantitative analyses of the α1mRNA from the reticulocytes of carriers highlighted a reduction in the variant cDNAs that constituted 34% (Hb Campania) and 15% (Hb Sciacca) of the total α1-globin cDNA, respectively. We developed a workflow for the in silico analysis of mechanisms triggering no-go decay, and its results suggested that the reduction in the variant mRNA was likely due to no-go decay caused by the presence of a rare triplet, and, in the case of Hb Sciacca, also by the mRNA’s secondary structure variation. It would be interesting to correlate the phenotype with the quantity of other frameshift mRNA variants, but very few data concerning α- and β-globin variants are available.

Identification and Characterization of Rice OsHKT1;3 Variants

In rice, the high-affinity K+ transporter, OsHKT1;3, functions as a Na+-selective transporter. mRNA variants of OsHKT1;3 have been reported previously, but their functions remain unknown. In this study, five OsHKT1;3 variants (V1-V5) were identified from japonica rice (Nipponbare) in addition to OsHKT1;3_FL. Absolute quantification qPCR analyses revealed that the transcript level of OsHKT1;3_FL was significantly higher than other variants in both the roots and shoots. Expression levels of OsHKT1;3_FL, and some variants, increased after 24 h of salt stress. Two electrode voltage clamp experiments in a heterologous expression system using Xenopus laevis oocytes revealed that oocytes expressing OsHKT1;3_FL and all of its variants exhibited smaller Na+ currents. The presented data, together with previous data, provide insights to understanding how OsHKT family members are involved in the mechanisms of ion homeostasis and salt tolerance in rice.

Significance of PD1 Alternative Splicing in Celiac Disease as a Novel Source for Diagnostic and Therapeutic Target

BackgroundWe have focused on the alteration of the PD-1/PD-L1 pathway in celiac disease and discussed the roles of the PD1 pathway in regulating the immune response. We explored the idea that the altered mRNA splicing process in key regulatory proteins could represent a novel source to identify diagnostic, prognostic, and therapeutic targets in celiac disease.MethodsWe characterized the PD1 mRNA variants’ profile in CD patients and in response to gluten peptides’ incubation after in vitro experiments. Total RNA from whole blood was isolated, and the coding region of the human PD-1 mRNA was amplified by cDNA PCR.ResultsPCR amplification of the human PD-1 coding sequence revealed an association between the over-expression of the sPD-1 protein and the PD-1Δex3 transcript in celiac disease. Thus, we have found three novel alternative spliced isoforms, two of which result in a truncated protein and the other isoform with a loss of 14 aa of exon 2 and complete exon 3 (Δ3) which could encode a new soluble form of PD1 (sPD-1).ConclusionsOur study provides evidence that dietary gluten can modulate processes required for cell homeostasis through the splicing of pre-mRNAs encoding key regulatory proteins, which represents an adaptive mechanism in response to different nutritional conditions.

Splice Variants of Superoxide Dismutases in Rice and Their Expression Profiles under Abiotic Stresses

The superoxide dismutases (SODs) play vital roles in controlling cellular reactive oxygen species (ROS) that are generated both under optimal as well as stress conditions in plants. The rice genome harbors seven SOD genes (CSD1, CSD2, CSD3, CSD4, FSD1, FSD2, and MSD) that encode seven constitutive transcripts. Of these, five (CSD2, CSD3, CSD4, FSD1, and MSD) utilizes an alternative splicing (AS) strategy and generate seven additional splice variants (SVs) or mRNA variants, i.e., three for CSD3, and one each for CSD2, CSD4, FSD1, and MSD. The exon-intron organization of these SVs revealed variations in the number and length of exons and/or untranslated regions (UTRs). We determined the expression patterns of SVs along with their constitutive forms of SODs in rice seedlings exposed to salt, osmotic, cold, heavy metal (Cu+2) stresses, as well as copper-deprivation. The results revealed that all seven SVs were transcriptionally active in both roots and shoots. When compared to their corresponding constitutive transcripts, the profiles of five SVs were almost similar, while two specific SVs (CSD3-SV4 and MSD-SV2) differed significantly, and the differences were also apparent between shoots and roots suggesting that the specific SVs are likely to play important roles in a tissue-specific and stress-specific manner. Overall, the present study has provided a comprehensive analysis of the SVs of SODs and their responses to stress conditions in shoots and roots of rice seedlings.

Differential Expression of BARD1 Isoforms in Melanoma

Melanoma comprises <5% of cutaneous malignancies, yet it causes a significant proportion of skin cancer-related deaths worldwide. While new therapies for melanoma have been developed, not all patients respond well. Thus, further research is required to better predict patient outcomes. Using long-range nanopore sequencing, RT-qPCR, and RNA sequencing analyses, we examined the transcription of BARD1 splice isoforms in melanoma cell lines and patient tissue samples. Seventy-six BARD1 mRNA variants were identified in total, with several previously characterised isoforms (γ, φ, δ, ε, and η) contributing to a large proportion of the expressed transcripts. In addition, we identified four novel splice events, namely, Δ(E3_E9), ▼(i8), IVS10+131▼46, and IVS10▼176, occurring in various combinations in multiple transcripts. We found that short-read RNA-Seq analyses were limited in their ability to predict isoforms containing multiple non-contiguous splicing events, as compared to long-range nanopore sequencing. These studies suggest that further investigations into the functional significance of the identified BARD1 splice variants in melanoma are warranted.

Distinct role of 5′UTR sequences in dendritic trafficking of BDNF mRNA: additional mechanisms for the BDNF splice variants spatial code

The neurotrophin Brain-derived neurotrophic factor (BDNF) is encoded by multiple bipartite transcripts. Each BDNF transcript is composed by one out of 11 alternatively spliced exons containing the 5′untranslated region (UTR), and one common exon encompassing the coding sequence (CDS) and the 3′UTR with two variants (short and long). In neurons, BDNF mRNA variants have a distinct subcellular distribution, constituting a “spatial code”, with exon 1, 3, 5, 7 and 8 located in neuronal somata, exon 4 extending into proximal dendrites, and exon 2 and 6 reaching distal dendrites. We previously showed that the CDS encodes constitutive dendritic targeting signals (DTS) and that both the 3′UTR-short and the 3′UTR-long contain activity-dependent DTS. However, the role of individual 5′UTR exons in mRNA sorting remains unclear. Here, we tested the ability of each different BDNF 5′UTRs to affect the subcellular localization of the green fluorescent protein (GFP) reporter mRNA. We found that exon 2 splicing isoforms (2a, 2b, and 2c) induced a constitutive dendritic targeting of the GFP reporter mRNA towards distal dendritic segments. The other isoforms did not affect GFP-mRNA dendritic trafficking. Through a bioinformatic analysis, we identified five unique cis-elements in exon 2a, 2b, and 2c which might contribute to building a DTS. This study provides additional information on the mechanism regulating the cellular sorting of BDNF mRNA variants.

Ultrasensitive quantification of multiplexed mRNA variants via splice-junction anchored DNA probes and SplintR ligase-initiated PCR

A method was developed for multiplexed quantification of mRNA splicing variants via splice-junction anchored DNA probes and SplintR ligase-initiated PCR.