Search

DNA Topoisomerase IIα (TOP2α/170) is an enzyme essential for proliferating cells. For rapidly multiplying malignancies, this has made TOP2α/170 an important target for etoposide and other clinically active anticancer drugs. Efficacy of these agents is often limited by chemoresistance related to alterations in TOP2α/170 expression levels. Our laboratory recently demonstrated reduced levels of TOP2α/170 and overexpression of a C-terminal truncated 90-kDa isoform, TOP2α/90, due to intronic polyadenylation (IPA; within intron 19) in an acquired etoposide-resistant K562 clonal cell line, K/VP.5. We previously reported that this isoform heterodimerized with TOP2α/170 and was a determinant of acquired resistance to etoposide. Optimization of the weak TOP2α exon 19/intron 19 5′ splice site in drug-resistant K/VP.5 cells by gene-editing restored TOP2α/170 levels, diminished TOP2α/90 expression, and circumvented drug resistance. Conversely, in the present study, silencing of the exon 19/intron 19 5′ splice site in parental K562 cells by CRISPR/Cas9 with homology-directed repair (HDR), and thereby forcing intron 19 retention, was used to induce resistance by disrupting normal RNA processing (i.e., gene knockout), and to further evaluate the role of TOP2α/170 and TOP2α/90 isoforms as resistance determinants. Gene-edited clones were identified by quantitative polymerase chain reaction (qPCR) and verified by Sanger sequencing. TOP2α/170 mRNA/protein expression levels were attenuated in the TOP2α gene-edited clones which resulted in resistance to etoposide as assessed by reduced etoposide-induced DNA damage (γH2AX, Comet assays) and growth inhibition. RNA-seq and qPCR studies suggested that intron 19 retention leads to decreased TOP2α/170 expression by degradation of the TOP2α edited mRNA transcripts. Forced expression of TOP2α/90 in the gene-edited K562 cells further decreased etoposide-induced DNA damage in support of a dominant negative role for this truncated isoform. Together results support the important role of both TOP2α/170 and TOP2α/90 as determinants of sensitivity/resistance to TOP2α-targeting agents.

Background It is estimated that up to 50% of all disease causing variants disrupt splicing. Due to its complexity, our ability to predict which variants disrupt splicing is limited, meaning missed diagnoses for patients. The emergence of machine learning for targeted medicine holds great potential to improve prediction of splice disrupting variants. The recently published SpliceAI algorithm utilises deep neural networks and has been reported to have a greater accuracy than other commonly used methods. Methods and findings The original SpliceAI was trained on splice sites included in primary isoforms combined with novel junctions observed in GTEx data, which might introduce noise and de-correlate the machine learning input with its output. Limiting the data to only validated and manual annotated primary and alternatively spliced GENCODE sites in training may improve predictive abilities. All of these gene isoforms were collapsed (aggregated into one pseudo-isoform) and the SpliceAI architecture was retrained (CI-SpliceAI). Predictive performance on a newly curated dataset of 1,316 functionally validated variants from the literature was compared with the original SpliceAI, alongside MMSplice, MaxEntScan, and SQUIRLS. Both SpliceAI algorithms outperformed the other methods, with the original SpliceAI achieving an accuracy of ∼91%, and CI-SpliceAI showing an improvement at ∼92% overall. Predictive accuracy increased in the majority of curated variants. Conclusions We show that including only manually annotated alternatively spliced sites in training data improves prediction of clinically relevant variants, and highlight avenues for further performance improvements.

Retinitis pigmentosa (RP) is a hereditary disease affecting tens of thousands of people world-wide. Here we analyzed the effect of an amino acid substitution in the RNA helicase DHX38 (Prp16) causing RP. DHX38 has been proposed as the helicase important for the 2nd step of splicing. We showed that DHX38 associates with key splicing factors involved in both splicing steps but did not find any evidence that the RP mutations changes DHX38 interaction profile with the spliceosome. We further downregulated DHX38 and monitored changes in splicing. We observed only minor perturbations of general splicing but detected modulation of ~70 alternative splicing events. Next, we probed DHX38 function in splicing of retina specific genes and found that FSCN2 splicing is dependent on DHX38. In addition, RHO splicing was inhibited specifically by expression of DHX38 RP variant. Finally, we showed that overexpression of DHX38 promotes usage of canonical as well as cryptic 5’ splice sites in HBB splicing reporter. Together, our data show that DHX38 is a splicing factor that promotes splicing of cryptic splice sites and regulate alternative splicing. We further provide evidence that the RP-linked substitution G332D modulates DHX38 splicing activity.

Computational algorithms are often used to assess pathogenicity of Variants of Uncertain Significance (VUS) that are found in disease-associated genes. Most computational methods include analysis of protein multiple sequence alignments (PMSA), assessing interspecies variation. Careful validation of PMSA-based methods has been done for relatively few genes, partially because creation of curated PMSAs is labor-intensive. We assessed how PMSA-based computational tools predict the effects of the missense changes in the APC gene, in which pathogenic variants cause Familial Adenomatous Polyposis. Most Pathogenic or Likely Pathogenic APC variants are protein-truncating changes. However, public databases now contain thousands of variants reported as missense. We created a curated APC PMSA that contained >3 substitutions/site, which is large enough for statistically robust in silico analysis. The creation of the PMSA was not easily automated, requiring significant querying and computational analysis of protein and genome sequences. Of 1924 missense APC variants in the NCBI ClinVar database, 1800 (93.5%) are reported as VUS. All but two missense variants listed as P/LP occur at canonical splice or Exonic Splice Enhancer sites. Pathogenicity predictions by five computational tools (Align-GVGD, SIFT, PolyPhen2, MAPP, REVEL) differed widely in their predictions of Pathogenic/Likely Pathogenic (range 17.5–75.0%) and Benign/Likely Benign (range 25.0–82.5%) for APC missense variants in ClinVar. When applied to 21 missense variants reported in ClinVar as Benign, the five methods ranged in accuracy from 76.2-100%. Computational PMSA-based methods can be an excellent classifier for variants of some hereditary cancer genes. However, there may be characteristics of the APC gene and protein that confound the results of in silico algorithms. A systematic study of these features could greatly improve the automation of alignment-based techniques and the use of predictive algorithms in hereditary cancer genes.Author SummaryA critical problem in clinical genetics today is interpreting whether a genetic variant is benign or causes disease (pathogenic). Some of the hardest variants to interpret are those that change one amino acid for another in a protein sequence (a “missense variant”). Various computer programs are often used to predict whether mutations in disease-associated genes likely cause disease. Most computer programs involve studying how the gene has changed during evolution, comparing the protein sequences of different species by aligning them with each other. Variants in amino acids that have not tolerated mutation during evolution are usually predicted to be pathogenic, and variants in amino acids that have tolerated variation are usually predicted to be benign. High quality alignments are necessary to make accurate predictions. However, creating high quality alignments is difficult, not easily automated, and requires significant manual curation. Results from computer-generated predictions are used in current published guidelines as one tool for evaluating whether variants will disrupt the protein function and cause disease. These guidelines may be applied to genes in which single amino acid substitutions do not commonly cause disease. One such example is the APC gene, which is responsible for Familial Adenomatous Polyposis (FAP). Missense APC changes are not a common cause of FAP. Our analysis of APC demonstrated the difficulty of generating an accurate protein sequence alignment and the tendency of computer tools to overestimate the damaging effects of amino acid substitutions. Our results suggest that the rules for using computer-based tools to predict whether a variant causes disease should be modified when applied to genes in which missense variants rarely cause disease.

Alternative splicing diversifies mRNA transcripts in human cells. While the spliceosome pairs exons with a high degree of accuracy, the rates of rare aberrant and non-canonical pre-mRNA splicing have not been evaluated at the nucleotide level to determine the quantity and identity of these events across splice junctions. Using ultra-deep sequencing the frequency of aberrant and non-canonical splicing events for three splice junctions flanking exon 7 of SMN1 were determined at single nucleotide resolution. After correction for background noise introduced by PCR amplification and sequencing steps, pre-mRNA splicing was shown to maintain a low overall rate of aberrant and non-canonically spliced events. Several previously unannotated splicing events across 3 exon|intron junctions in SMN1 were identified. Mutations within SMN exon 7 were shown to affect splicing fidelity by modulating RNA secondary structures, by altering the binding site of regulatory proteins and by changing the 5' splice site strength. Mutations also create a truncated SMN1 exon 7 through the introduction of a de novo non-canonical 5' splice site. The results from the ultra-deep sequencing approach highlight the impressive fidelity of pre-mRNA splicing and demonstrate that the immediate sequence context around splice sites is the main driving force behind non-canonical splice site pairing.

Metabolic labeling is a widely used tool to investigate different aspects of pre-mRNA splicing and RNA turnover. The labeling technology takes advantage of native cellular machineries where a nucleotide analog is readily taken up and incorporated into nascent RNA. One such analog is 4-thiouridine (4sU). Previous studies demonstrated that the uptake of 4sU at elevated concentrations (>50μM) and extended exposure led to inhibition of rRNA synthesis and processing, presumably induced by changes in RNA secondary structure. Thus, it is possible that 4sU incorporation may also interfere with splicing efficiency. To test this hypothesis, we carried out splicing analyses of pre-mRNA substrates with varying levels of 4sU incorporation (0–100%). We demonstrate that increased incorporation of 4sU into pre-mRNAs decreased splicing efficiency. The overall impact of 4sU labeling on pre-mRNA splicing efficiency negatively correlates with the strength of splice site signals such as the 3’ and the 5’ splice sites. Introns with weaker splice sites are more affected by the presence of 4sU. We also show that transcription by T7 polymerase and pre-mRNA degradation kinetics were impacted at the highest levels of 4sU incorporation. Increased incorporation of 4sU caused elevated levels of abortive transcripts, and fully labeled pre-mRNA is more stable than its uridine-only counterpart. Cell culture experiments show that a small number of alternative splicing events were modestly, but statistically significantly influenced by metabolic labeling with 4sU at concentrations considered to be tolerable (40 μM). We conclude that at high 4sU incorporation rates small, but noticeable changes in pre-mRNA splicing can be detected when splice sites deviate from consensus. Given these potential 4sU artifacts, we suggest that appropriate controls for metabolic labeling experiments need to be included in future labeling experiments.

CRISPR/Cas9 allows the generation of knockout cell lines and null zygotes by inducing site-specific double-stranded breaks. In most cases the DSB is repaired by non-homologous end joining, resulting in small nucleotide insertions or deletions that can be used to construct knockout alleles. However, these mutations do not produce the desired null result in all cases, but instead generate a similar, functionally active protein. This effect could limit the therapeutic efficiency of gene therapy strategies based on abrogating oncogene expression, and therefore needs to be considered carefully. If there is an acceptable degree of efficiency of CRISPR/Cas9 delivery to cells, the key step for success lies in the effectiveness of a specific sgRNA at knocking out the oncogene, when only one sgRNA can be used. This study shows that the null effect could be increased with an sgRNA targeting the splice donor site (SDS) of the chosen exon. Following this strategy, the generation of null alleles would be facilitated in two independent ways: the probability of producing a frameshift mutation and the probability of interrupting the canonical mechanism of pre-mRNA splicing. In these contexts, we propose to improve the loss-of-function yield driving the CRISPR system at the SDS of critical exons., This work was mainly supported by a grant from the Fondo de Investigaciones Sanitarias (FIS) of the Spanish Ministry of Economy and Competitiveness and the European Regional Development Fund (ERDF) “Una manera de hacer Europa” [grant PI17/01895 to IGT and MSM.]; Junta de Castilla y León, Fondos FEDER [SA085U16 to JMHR]; Novartis grant; and by the Fundación “Jabones para Daniel”. JM Hernández-Sánchez was supported by a research grant from Fundación Española de Hematología y Hemoterapia (FEHH).

Pre-mRNA splicing is critical for achieving required amounts of a transcript at a given time and for regulating production of encoded protein. A given pre-mRNA may be spliced in many ways, or not at all, giving rise to multiple gene products. Numerous splicing factors are recruited to pre-mRNA splice sites to ensure proper splicing. One such factor, the 60 kDa poly(U)-binding splicing factor (PUF60), is recruited to sites that are not always spliced, but rather function as alternative splice sites. In this study, we characterized the interaction of PUF60 with a splice site from the adenovirus major late promoter (the AdML 3' splice site, AdML3’). We found that the PUF60–AdML3’ dissociation constants are in the micromolar range, with the binding affinity predominantly provided by PUF60’s two central RNA recognition motifs (RRMs). A 1.95 Å crystal structure of the two PUF60 RRMs in complex with AdML3’ revealed a dimeric organization placing two stretches of nucleic acid tracts in opposing directionalities, which can cause looping of nucleic acid and explain how PUF60 affects pre-mRNA geometry to effect splicing. Solution characterization of this complex by light-scattering and UV/Vis spectroscopy suggested a potential 2:1 (PUF602:AdML3’) stoichiometry, consistent with the crystal structure. This work defines the sequence specificity of the alternative splicing factor PUF60 at the pre-mRNA 3’ splice site. Our observations suggest that control of pre-mRNA directionality is important in the early stage of spliceosome assembly, and advance our understanding of the molecular mechanism by which alternative and constitutive splicing factors differentiate among 3’ splice sites.

The craniofacial developmental disorder Burn-McKeown Syndrome (BMKS) is caused by biallelic variants in the pre-messenger RNA splicing factor gene TXNL4A/DIB1. The majority of affected individuals with BMKS have a 34 base pair deletion in the promoter region of one allele of TXNL4A combined with a loss-of-function variant on the other allele, resulting in reduced TXNL4A expression. However, it is unclear how reduced expression of this ubiquitously expressed spliceosome protein results in craniofacial defects during development. Here we reprogrammed peripheral mononuclear blood cells from a BMKS patient and her unaffected mother into induced pluripotent stem cells (iPSCs) and differentiated the iPSCs into induced neural crest cells (iNCCs), the key cell type required for correct craniofacial development. BMKS patient-derived iPSCs proliferated more slowly than both mother- and unrelated control-derived iPSCs, and RNA-Seq analysis revealed significant differences in gene expression and alternative splicing. Patient iPSCs displayed defective differentiation into iNCCs compared to maternal and unrelated control iPSCs, in particular a delay in undergoing an epithelial-to-mesenchymal transition (EMT). RNA-Seq analysis of differentiated iNCCs revealed widespread gene expression changes and mis-splicing in genes relevant to craniofacial and embryonic development that highlight a dampened response to WNT signalling, the key pathway activated during iNCC differentiation. Furthermore, we identified the mis-splicing of TCF7L2 exon 4, a key gene in the WNT pathway, as a potential cause of the downregulated WNT response in patient cells. Additionally, mis-spliced genes shared common sequence properties such as length, splice site strengths and sequence motifs, suggesting that splicing of particular subsets of genes is particularly sensitive to changes in TXNL4A expression. Together, these data provide the first insight into how reduced TXNL4A expression in BMKS patients might compromise splicing and NCC function, resulting in defective craniofacial development in the embryo.

Budding yeast spliceosomal factors ScSlu7 and ScPrp18 interact and mediate intron 3'ss choice during second step pre-mRNA splicing. The fission yeast genome with abundant multi-intronic transcripts, degenerate splice signals and SR proteins is an apt unicellular fungal model to deduce roles for core spliceosomal factors in alternative splice-site choice, intron retention and to study the cellular implications of regulated splicing. From our custom microarray data we deduce a stringent reproducible subset of S. pombe alternative events. We examined the role of factors SpSlu7 or SpPrp18 for these splice events and investigated the relationship to growth phase and stress. Wild-type log and stationary phase cells showed ats1+ exon 3 skipped and intron 3 retained transcripts. Interestingly the non-consensus 5' ss in ats1+ intron 3 caused SpSlu7 and SpPrp18 dependent intron retention. We validated the use of an alternative 5' ss in dtd1+ intron 1 and of an upstream alternative 3' ss in DUF3074 intron 1. The dtd1+ intron 1 non-canonical 5' ss yielded an alternative mRNA whose levels increased in stationary phase. Utilization of dtd1+ intron 1 sub-optimal 5' ss required functional SpPrp18 and SpSlu7 while compromise in SpSlu7 function alone hampered the selection of the DUF3074 intron 1 non canonical 3' ss. We analysed the relative abundance of these splice isoforms during mild thermal, oxidative and heavy metal stress and found stress-specific splice patterns for ats1+ and DUF3074 intron 1 some of which were SpSlu7 and SpPrp18 dependent. By studying ats1+ splice isoforms during compromised transcription elongation rates in wild-type, spslu7-2 and spprp18-5 mutant cells we found dynamic and intron context-specific effects in splice-site choice. Our work thus shows the combinatorial effects of splice site strength, core splicing factor functions and transcription elongation kinetics to dictate alternative splice patterns which in turn serve as an additional recourse of gene regulation in fission yeast.

Alternative splicing is an ubiquitous phenomenon in most human genes and has important functions. The switch-like exon is the type of exon that has a high level of usage in some tissues, but has a low level of usage in the other tissues. They usually undergo strong tissue-specific regulations. There is still a lack a systematic method to identify switch-like exons from multiple RNA-seq samples. We proposed a novel method called iterative Tertile Absolute Deviation around the mode (iTAD) to profile the distribution of exon relative usages among multiple samples and to identify switch-like exons and other types of exons using a robust statistic estimator. We validated the method with simulation data, and applied it on RNA-seq data of 16 human body tissues and detected 3,100 switch-like exons. We found that switch-like exons tend to be more associated with Alu elements in their flanking intron regions than other types of exons.

Genetic testing of PKD1 and PKD2 is expected to play an increasingly important role in determining allelic influences in autosomal dominant polycystic kidney disease (ADPKD) in the near future. However, to date, genetic testing is not commonly employed because it is expensive, complicated because of genetic heterogeneity, and does not easily identify pathogenic variants. In this study, we developed a genetic testing system based on next-generation sequencing (NGS), long-range polymerase chain reaction, and a new software package. The new software package integrated seven databases and provided access to five cloud-based computing systems. The database integrated 241 polymorphic nonpathogenic variants detected in 140 healthy Japanese volunteers aged >35 years, who were confirmed by ultrasonography as having no cysts in either kidney. Using this system, we identified 60 novel and 30 known pathogenic mutations in 101 Japanese patients with ADPKD, with an overall detection rate of 89.1% (90/101) [95% confidence interval (CI), 83.0%-95.2%]. The sensitivity of the system increased to 93.1% (94/101) (95% CI, 88.1%-98.0%) when combined with multiplex ligation-dependent probe amplification analysis, making it sufficient for use in a clinical setting. In 82 (87.2%) of the patients, pathogenic mutations were detected in PKD1 (95% CI, 79.0%-92.5%), whereas in 12 (12.8%) patients pathogenic mutations were detected in PKD2 (95% CI, 7.5%-21.0%); this is consistent with previously reported findings. In addition, we were able to reconfirm our pathogenic mutation identification results using Sanger sequencing. In conclusion, we developed a high-sensitivity NGS-based system and successfully employed it to identify pathogenic mutations in PKD1 and PKD2 in Japanese patients with ADPKD.

Disease-causing germline mutations in CDH1 cause Hereditary Diffuse Gastric Cancer (HDGC). For patients who meet the HDGC screening criteria, the identification and classification of the sequence variants found in CDH1 are critical for risk management of patients. In this report, we describe a germline CDH1 c.1679C>G (p.T560R) variant identified in a 50 year old man who was diagnosed with gastric cancer with a strong family history of gastric cancer (one living brother was diagnosed with gastric cancer at 63 and another brother died of gastric cancer at 45). cDNA analysis, involving fragment analysis and cloning, indicated that the p.T560R mutation created a novel 5' splice donor site, which led to a novel transcript with a 32 nucleotide deletion in exon 11. This abnormal transcript putatively produces a truncated CDH1 protein (E-cadherin) of 575 amino acids instead of 882. We also demonstrated that the variant completely abolishes normal splicing as the mutant allele does not generate any normal transcript. Furthermore, the CDH1 c.1679C>G (p.T560R) variant segregated with gastric cancer in all three family members affected with gastric cancer in this family. These results support the conclusion that CDH1 c.1679C>G (p.T560R) variant is a pathogenic mutation and contributes to HDGC through disruption of normal splicing.

HIV-1 RNAs are generated through a complex splicing mechanism, resulting in a great diversity of transcripts, which are classified in three major categories: unspliced, singly spliced (SS), and doubly spliced (DS). Knowledge on HIV-1 RNA splicing in vivo and by non-subtype B viruses is scarce. Here we analyze HIV-1 RNA splice site usage in CD4+CD25+ lymphocytes from HIV-1-infected individuals through pyrosequencing. HIV-1 DS and SS RNAs were amplified by RT-PCR in 19 and 12 samples, respectively. 13,108 sequences from HIV-1 spliced RNAs, derived from viruses of five subtypes (A, B, C, F, G), were identified. In four samples, three of non-B subtypes, five 3' splice sites (3'ss) mapping to unreported positions in the HIV-1 genome were identified. Two, designated A4i and A4j, were used in 22% and 25% of rev RNAs in two viruses of subtypes B and A, respectively. Given their close proximity (one or two nucleotides) to A4c and A4d, respectively, they could be viewed as variants of these sites. Three 3'ss, designated A7g, A7h, and A7i, located 20, 32, and 18 nucleotides downstream of A7, respectively, were identified in a subtype C (A7g, A7h) and a subtype G (A7i) viruses, each in around 2% of nef RNAs. The new splice sites or variants of splice sites were associated with the usual sequence features of 3'ss. Usage of unusual 3'ss A4d, A4e, A5a, A7a, and A7b was also detected. A4f, previously identified in two subtype C viruses, was preferentially used by rev RNAs of a subtype C virus. These results highlight the great diversity of in vivo splice site usage by HIV-1 RNAs. The fact that four of five newly identified splice sites or variants of splice sites were detected in non-subtype B viruses allows anticipating an even greater diversity of HIV-1 splice site usage than currently known. This work was supported by the Francisco Díez-Fuertes is supported by the Sara Borrell postdoctoral Program of the Spanish Government 2012 CD12/00515. This work was also supported by Ministerio de Economía y Competitividad, Spain (URL: http://www.mineco.gob.es/portal/site/mineco/) Grant # SAF2010-20961, Grant # SAF2013-48488-P to MMT and Instituto de Salud Carlos III (URL: http://www.isciii.es) Grant # MPY 1194/12 to MMT. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Sí

Cellular mechanisms that achieve protein diversity in eukaryotes are multifaceted, including transcriptional components such as RNA splicing. Through alternative splicing, a single protein-coding gene can generate multiple mRNA transcripts and protein isoforms, some of which are tissue-specific. We have conducted qualitative and quantitative analyses of the Bodymap 2.0 messenger RNA-sequencing data from 16 human tissue samples and identified 209,363 splice junctions. Of these, 22,231 (10.6%) were not previously annotated and 21,650 (10.3%) were expressed in a tissue-restricted pattern. Tissue-restricted alternative splicing was found to be widespread, with approximately 65% of expressed multi-exon genes containing at least one tissue-specific splice junction. Interestingly, we observed many tissue-specific splice junctions not only in genes expressed in one or a few tissues, but also from gene loci with a broad pattern of expression.

BACKGROUND:Mild unconjugated hyperbilirubinemia (UH), due to reduced activity of the enzyme uridine diphosphoglucuronate-glucuronosyltransferase family, polypeptide 1 (UGT1A1), is a common clinical condition. Most cases are caused by presence in homozygous form of an A(TA)7TAA nucleotide sequence instead of the usual A(TA)6TAA sequence in promoter region of the UGT1A1 gene. In some cases, other genetic variations have been identified which differ between populations. There is need for more data on such genetic variations from India. METHODS:DNA from subjects with unexplained persistent or recurrent UH was tested for the presence of TA promoter insertions. In addition, all five exons and splicing site regions of UGT1A1 gene were sequenced. Several bioinformatics tools were used to determine the biological significance of the observed genetic changes. Functional analysis was done to look for effect of a splice site mutation in UGT1A1. RESULTS:Of 71 subjects with UH (68 male; median age [range], 26 [16-63] years; serum bilirubin 56 [26-219] μM/L, predominantly unconjugated) studied, 65 (91.5%) subjects were homozygous for A(TA)7TAA allele, five (7.0%) were heterozygous, and one (1.4%) lacked this change. Fifteen subjects with UH had missense exonic single nucleotide changes (14 heterozygous, 1 homozygous), including one subject with a novel nucleotide change (p.Thr205Asn). Bioinformatics tools predicted some of these variations (p.Arg108Cys, p.Ile159Thr and p.Glu463Val) to be deleterious. Functional characterization of an exonic variation (c.1084G>A) located at a splice site revealed that it results in frameshift deletion of 31 nucleotides and premature truncation of the protein. CONCLUSION:Our study revealed several single nucleotide variations in UGT1A1 gene in Indian subjects with UH. Functional characterization of a splice site variation indicated that it leads to disordered splicing. These variations may explain UH in subjects who lacked homozygous A(TA)7TAA promoter alleles.

Alternative splicing (AS) is pervasive in human multi-exon genes and is a major contributor to expansion of the transcriptome and proteome diversity. The accurate recognition of alternative splice sites is regulated by information contained in networks of protein-protein and protein-RNA interactions. However, the mechanisms leading to splice site selection are not fully understood. Although numerous databases have been built to describe AS, molecular interaction databases associated with AS have only recently emerged. In this study, we present a new database, MiasDB, that provides a description of molecular interactions associated with human AS events. This database covers 938 interactions between human splicing factors, RNA elements, transcription factors, kinases and modified histones for 173 human AS events. Every entry includes the interaction partners, interaction type, experimental methods, AS type, tissue specificity or disease-relevant information, a simple description of the functionally tested interaction in the AS event and references. The database can be queried easily using a web server (http://47.88.84.236/Miasdb). We display some interaction figures for several genes. With this database, users can view the regulation network describing AS events for 12 given genes.

Objective Hypophosphatemic rickets (HR) is a heterogeneous genetic phosphate wasting disorder. The disease is most commonly caused by mutations in the PHEX gene located on the X-chromosome or by mutations in CLCN5, DMP1, ENPP1, FGF23, and SLC34A3. The aims of this study were to perform molecular diagnostics for four patients with HR of Indian origin (two independent families) and to describe their clinical features. Methods We performed whole exome sequencing (WES) for the affected mother of two boys who also displayed the typical features of HR, including bone malformations and phosphate wasting. B-lymphoblast cell lines were established by EBV transformation and subsequent RT-PCR to investigate an uncommon splice site variant found by WES. An in silico analysis was done to obtain accurate nucleotide frequency occurrences of consensus splice positions other than the canonical sites of all human exons. Additionally, we applied direct Sanger sequencing for all exons and exon/intron boundaries of the PHEX gene for an affected girl from an independent second Indian family. Results WES revealed a novel PHEX splice acceptor mutation in intron 9 (c.1080-3C>A) in a family with 3 affected individuals with HR. The effect on splicing of this mutation was further investigated by RT-PCR using RNA obtained from a patient’s EBV-transformed lymphoblast cell line. RT-PCR revealed an aberrant splice transcript skipping exons 10-14 which was not observed in control samples, confirming the diagnosis of X-linked dominant hypophosphatemia (XLH). The in silico analysis of all human splice sites adjacent to all 327,293 exons across 81,814 transcripts among 20,345 human genes revealed that cytosine is, with 64.3%, the most frequent nucleobase at the minus 3 splice acceptor position, followed by thymidine with 28.7%, adenine with 6.3%, and guanine with 0.8%. We generated frequency tables and pictograms for the extended donor and acceptor splice consensus regions by analyzing all human exons. Direct Sanger sequencing of all PHEX exons in a sporadic case with HR from the Indian subcontinent revealed an additional novel PHEX mutation (c.1211_1215delACAAAinsTTTACAT, p.Asp404Valfs*5, de novo) located in exon 11. Conclusions Mutation analyses revealed two novel mutations and helped to confirm the clinical diagnoses of XLH in two families from India. WES helped to analyze all genes implicated in the underlying disease complex. Mutations at splice positions other than the canonical key sites need further functional investigation to support the assertion of pathogenicity.

Androgen receptor (AR) variants are associated with resistance to anti androgen therapy both in human prostate cancer cell lines and clinical samples. These observations support the hypothesis that AR isoform accumulation is a consequence of selective therapeutic pressure on the full length AR. The Pten deficient prostate cancer model proceeds with well-defined kinetics including progression to castration resistant prostate cancer (CRPC). While surgical castration and enzalutamide treatments yield an initial therapeutic response, Pten-/-epithelia continue to proliferate yielding locally invasive primary tumor pathology. That most epithelium remains AR positive, but ligand independent, suggests the presence of oncogenic AR variants. To address this hypothesis, we have used a panel of recently described Pten-/- tumor cell lines derived from both from hormone intact (E4, E8) and castrated Pten mutants (cE1, cE2) followed by RACE PCR to identify and characterize three novel truncated, amino terminus containing AR variants (mAR-Va, b, c). Variants appear not only conserved throughout progression but are correlated with nearly complete loss of full length AR (AR-FL) at castrate androgen levels. The overexpression of variants leads to enhanced transcriptional activity of AR while knock down studies show reduced transcriptional output. Collectively, the identification of truncated AR variants in the conditional PTEN deletion model supports a role for maintaining the CRPC phenotype and provides further therapeutic applications of this preclinical model.

Congenital tooth agenesis is caused by mutations in the MSX1, PAX9, WNT10A, or AXIN2 genes. Here, we report a Japanese family with nonsyndromic tooth agenesis caused by a novel nucleotide substitution in the intronic region between exons 1 and 2 of the MSX1 gene. Because the mutation is located 9 bp before exon 2 (c.452-9G>A), we speculated that the nucleotide substitution would generate an abnormal splice site. Using cDNA analysis of an immortalized patient blood cell, we confirmed that an additional 7-nucleotide sequence was inserted at the splice junction between exons 1 and 2 (c.451_452insCCCTCAG). The consequent frameshift generated a homeodomain-truncated MSX1 (p.R151fsX20). We then studied the subcellular localization of truncated MSX1 protein in COS cells, and observed that it had a whole cell distribution more than a nuclear localization, compared to that of wild-type protein. This result suggests a deletion of the nuclear localization signal, which is mapped to the MSX1 homeodomain. These results indicate that this novel intronic nucleotide substitution is the cause of tooth agenesis in this family. To date, most MSX1 variants isolated from patients with tooth agenesis involve single amino acid substitutions in the highly conserved homeodomain or deletion mutants caused by frameshift or nonsense mutations. We here report a rare case of an intronic mutation of the MSX1 gene responsible for human tooth agenesis. In addition, the missing tooth patterns were slightly but significantly different between an affected monozygotic twin pair of this family, showing that epigenetic or environmental factors also affect the phenotypic variations of missing teeth among patients with nonsyndromic tooth agenesis caused by an MSX1 haploinsufficiency.

International audience; Presently, human collagen VI-related diseases such as Ullrich congenital muscular dystrophy (UCMD) and Bethlem myopathy (BM) remain incurable, emphasizing the need to unravel their etiology and improve their treatments. In UCMD, symptom onset occurs early, and both diseases aggravate with ageing. In zebrafish fry, morpholinos reproduced early UCMD and BM symptoms but did not allow to study the late phenotype. Here, we produced the first zebrafish line with the human mutation frequently found in collagen VI-related disorders such as UCMD and BM. We used a transcription activator-like effector nuclease (TALEN) to design the col6a1ama605003-line with a mutation within an essential splice donor site, in intron 14 of the col6a1 gene, which provoke an in-frame skipping of exon 14 in the processed mRNA. This mutation at a splice donor site is the first example of a template-independent modification of splicing induced in zebrafish using a targetable nuclease. This technique is readily expandable to other organisms and can be instrumental in other disease studies. Histological and ultrastructural analyzes of homozygous and heterozygous mutant fry and 3 months post-fertilization (mpf) fish revealed co-dominantly inherited abnormal myofibers with disorganized myofibrils, enlarged sarcoplasmic reticulum, altered mitochondria and misaligned sarcomeres. Locomotion analyzes showed hypoxia-response behavior in 9 mpf col6a1 mutant unseen in 3 mpf fish. These symptoms worsened with ageing as described in patients with collagen VI deficiency. Thus, the col6a1ama605003-line is the first adult zebrafish model of collagen VI-related diseases; it will be instrumental both for basic research and drug discovery assays focusing on this type of disorders.

Ichthyoses are a group of various different types of hereditary disorders affecting skin cornification. They are characterized by hyperkeratoses of different severity levels and are associated with a dry and scaling skin. Genome-wide association analysis of nine affected and 13 unaffected Great Danes revealed a genome-wide significant peak on chromosome 9 at 57–58 Mb in the region of SLC27A4. Sequence analysis of genomic DNA of SLC27A4 revealed the non-synonymous SNV SLC27A4:g.8684G>A in perfect association with ichthyosis-affection in Great Danes. The mutant transcript of SLC27A4 showed an in-frame loss of 54 base pairs in exon 8 probably induced by a new splice acceptor site motif created by the mutated A- allele of the SNV. Genotyping 413 controls from 35 different breeds of dogs and seven wolves revealed that this mutation could not be found in other populations except in Great Danes. Affected dogs revealed high amounts of mutant transcript but only low levels of the wild type transcript. Targeted analyses of SLC27A4 protein from skin tissues of three affected and two unaffected Great Danes indicated a markedly reduced or not detectable wild type and truncated protein levels in affected dogs but a high expression of wild type SLC27A4 protein in unaffected controls. Our data provide evidence of a new splice acceptor site creating SNV that results in a reduction or loss of intact SLC27A4 protein and probably explains the severe skin phenotype in Great Danes. Genetic testing will allow selective breeding to prevent ichthyosis-affected puppies in the future.

The plant pararetrovirus Cauliflower mosaic virus (CaMV) uses alternative splic-ing to generate several isoforms from its polycistronic pregenomic 35S RNA. This pro-cess has been shown to be essential for infectivity. Previous works have identified four splice donor sites and a single splice acceptor site in the 35S RNA 5’ region and sug-gested that the main role of CaMV splicing is to downregulate expression of open read-ing frames (ORFs) I and II. In this study, we show that alternative splicing is a conserved process among CaMV isolates. In Cabb B-JI and Cabb-S isolates, splicing frequently leads to different fusion between ORFs, particularly between ORF I and II. The corresponding P1P2 fusion proteins expressed in E. coli interact with viral proteins P2 and P3 in vitro. However, they are detected neither during infection nor upon transient expression in planta, which suggests rapid degradation after synthesis and no important biological role in the CaMV infectious cycle. To gain a better understanding of the functional relevance of 35S RNA alternative splicing in CaMV infectivity, we inactivated the previously described splice sites. All the splicing mutants were as pathogenic as the corresponding wild-type isolate. Through RT-PCR-based analysis we demonstrate that CaMV 35S RNA exhibits a complex splicing pattern, as we identify new splice donor and acceptor sites whose selection leads to more than thirteen 35S RNA isoforms in infected turnip plants. Inactivating splice donor or acceptor sites is not lethal for the virus, since disrupted sites are systematically rescued by the activation of cryptic and/or seldom used splice sites. Taken together, our data depict a conserved, complex and flexible process, involving multiple sites, that ensures splicing of 35S RNA.

Splice sites (SSs)--short nucleotide sequences flanking introns--are under selection for spliceosome binding, and adhere to consensus sequences. However, non-consensus nucleotides, many of which probably reduce SS performance, are frequent. Little is known about the mechanisms maintaining such apparently suboptimal SSs. Here, we study the correlations between strengths of nucleotides occupying different positions of the same SS. Such correlations may arise due to epistatic interactions between positions (i.e., a situation when the fitness effect of a nucleotide in one position depends on the nucleotide in another position), their evolutionary history, or to other reasons. Within both the intronic and the exonic parts of donor SSs, nucleotides that increase (decrease) SS strength tend to co-occur with other nucleotides increasing (respectively, decreasing) it, consistent with positive epistasis. Between the intronic and exonic parts of donor SSs, the correlations of nucleotide strengths tend to be negative, consistent with negative epistasis. In the course of evolution, substitutions at a donor SS tend to decrease the strength of its exonic part, and either increase or do not change the strength of its intronic part. In acceptor SSs, the situation is more complicated; the correlations between adjacent positions appear to be driven mainly by avoidance of the AG dinucleotide which may cause aberrant splicing. In summary, both the content and the evolution of SSs is shaped by a complex network of interdependences between adjacent nucleotides that respond to a range of sometimes conflicting selective constraints.

Alternative splicing is an important biological process in the generation of multiple functional transcripts from the same genomic sequences. Differential analysis of splice junctions (SJs) and intron retentions (IRs) is helpful in the detection of alternative splicing events. In this study, we conducted differential analysis of SJs and IRs by use of DEXSeq, a Bioconductor package originally designed for differential exon usage analysis in RNA-seq data analysis. We set up an analysis pipeline including mapping of RNA-seq reads, the preparation of count tables of SJs and IRs as the input files, and the differential analysis in DEXSeq. We analyzed the public RNA-seq datasets generated from RNAi experiments on Drosophila melanogaster S2-DRSC cells to deplete RNA-binding proteins (GSE18508). The analysis confirmed previous findings on the alternative splicing of the trol and Ant2 (sesB) genes in the CG8144 (ps)-depletion experiment and identified some new alternative splicing events in other RNAi experiments. We also identified IRs that were confirmed in our SJ analysis. The proposed method used in our study can output the genomic coordinates of differentially used SJs and thus enable sequence motif search. Sequence motif search and gene function annotation analysis helped us infer the underlying mechanism in alternative splicing events. To further evaluate this method, we also applied the method to public RNA-seq data from human breast cancer (GSE45419) and the plant Arabidopsis (SRP008262). In conclusion, our study showed that DEXSeq can be adapted to differential analysis of SJs and IRs, which will facilitate the identification of alternative splicing events and provide insights into the molecular mechanisms of transcription processes and disease development.

Signal transducer and activator of transcription 3 (STAT3) is a key mediator of leukocyte differentiation and proliferation. The 3' end of STAT3 transcripts is subject to two alternative splicing events. One results in either full-length STAT3α or in STAT3β, which lacks part of the C-terminal transactivation domain. The other is at a tandem donor (5') splice site and results in the codon for Ser-701 being included (S) or excluded (ΔS). Despite the proximity of Ser-701 to the site of activating phosphorylation at Tyr-705, ΔS/S splicing has barely been studied. Sequencing of cDNA from purified eosinophils revealed the presence of four transcripts (S-α, ΔS-α, S-β, and ΔS-β) rather than the three reported in publically available databases from which ΔS-β is missing. To gain insight into regulation of the two alternative splicing events, we developed a quantitative(q) PCR protocol to compare transcript ratios in eosinophils in which STAT3 is upregulated by cytokines, activated B cell diffuse large B cell Lymphoma (DLBCL) cells in which STAT3 is dysregulated, and in germinal center B cell-like DLBCL cells in which it is not. With the exception of one line of activated B cell DLCBL cells, the four variants were found in roughly the same ratios despite differences in total levels of STAT3 transcripts. S-α was the most abundant, followed by S-β. ΔS-α and ΔS-β together comprised 15.6 ± 4.0 % (mean ± SD, n = 21) of the total. The percentage of STAT3β variants that were ΔS was 1.5-fold greater than of STAT3α variants that were ΔS. Inspection of Illumina's "BodyMap" RNA-Seq database revealed that the ΔS variant accounts for 10-26 % of STAT3 transcripts across 16 human tissues, with less variation than three other genes with the identical tandem donor splice site sequence. Thus, it seems likely that all cells contain the S-α, ΔS-α, S-β, and ΔS-β variants of STAT3.

RNA-seq is becoming a preferred tool for genomics studies of model and non-model organisms. However, DNA-based analysis of organisms lacking sequenced genomes cannot rely on RNA-seq data alone to isolate most genes of interest, as DNA codes both exons and introns. With this in mind, we designed a novel tool, LEMONS, that exploits the evolutionary conservation of both exon/intron boundary positions and splice junction recognition signals to produce high throughput splice-junction predictions in the absence of a reference genome. When tested on multiple annotated vertebrate mRNA data, LEMONS accurately identified 87% (average) of the splice-junctions. LEMONS was then applied to our updated Mediterranean chameleon transcriptome, which lacks a reference genome, and predicted a total of 90,820 exon-exon junctions. We experimentally verified these splice-junction predictions by amplifying and sequencing twenty randomly selected genes from chameleon DNA templates. Exons and introns were detected in 19 of 20 of the positions predicted by LEMONS. To the best of our knowledge, LEMONS is currently the only experimentally verified tool that can accurately predict splice-junctions in organisms that lack a reference genome.

BAP1 is a tumor suppressor gene that is lost or deleted in diverse cancers, including uveal mela¬noma, malignant pleural mesothelioma (MPM), clear cell renal carcinoma, and cholangiocarcinoma. Recently, BAP1 germline mutations have been reported in families with combinations of these same cancers. A particular challenge for mutation screening is the classification of non-truncating BAP1 sequence variants because it is not known whether these subtle changes can affect the protein function sufficiently to predispose to cancer development. Here we report mRNA splicing analysis on a homozygous substitution mutation, BAP1 c. 2054 A&T (p.Glu685Val), identified in an MPM cell line derived from a mesothelioma patient. The mutation occurred at the 3rd nucleotide from the 3' end of exon 16. RT-PCR, cloning and subsequent sequencing revealed several aberrant splicing products not observed in the controls: 1) a 4 bp deletion at the end of exon 16 in all clones derived from the major splicing product. The BAP1 c. 2054 A&T mutation introduced a new 5' splice site (GU), which resulted in the deletion of 4 base pairs and presumably protein truncation; 2) a variety of alternative splicing products that led to retention of different introns: introns 14-16; introns 15-16; intron 14 and intron 16; 3) partial intron 14 and 15 retentions caused by activation of alternative 3' splice acceptor sites (AG) in the introns. Taken together, we were unable to detect any correctly spliced mRNA transcripts in this cell line. These results suggest that aberrant splicing caused by this mutation is quite efficient as it completely abolishes normal splicing through creation of a novel 5' splice site and activation of cryptic splice sites. These data support the conclusion that BAP1 c.2054 A&T (p.E685V) variant is a pathogenic mutation and contributes to MPM through disruption of normal splicing.

Schizophrenia (SCZ) is a severe, debilitating mental illness which has a significant genetic component. The identification of genetic factors related to SCZ has been challenging and these factors remain largely unknown. To evaluate the contribution of de novo variants (DNVs) to SCZ, we sequenced the exomes of 53 individuals with sporadic SCZ and of their non-affected parents. We identified 49 DNVs, 18 of which were predicted to alter gene function, including 13 damaging missense mutations, 2 conserved splice site mutations, 2 nonsense mutations, and 1 frameshift deletion. The average number of exonic DNV per proband was 0.88, which corresponds to an exonic point mutation rate of 1.7×10(-8) per nucleotide per generation. The non-synonymous-to-synonymous mutation ratio of 2.06 did not differ from neutral expectations. Overall, this study provides a list of 18 putative candidate genes for sporadic SCZ, and when combined with the results of similar reports, identifies a second proband carrying a non-synonymous DNV in the RGS12 gene.

Myotonia congenita (MC) is a skeletal muscle channelopathy characterized by inability of the muscle to relax following voluntary contraction. Worldwide population prevalence in humans is 1:100,000. Studies in mice, dogs, humans and goats confirmed myotonia associated with functional defects in chloride channels and mutations in a skeletal muscle chloride channel (CLCN1). CLCN1 encodes for the most abundant chloride channel in the skeletal muscle cell membrane. Five random bred cats from Winnipeg, Canada with MC were examined. All cats had a protruding tongue, limited range of jaw motion and drooling with prominent neck and proximal limb musculature. All cats had blepharospasm upon palpebral reflex testing and a short-strided gait. Electromyograms demonstrated myotonic discharges at a mean frequency of 300 Hz resembling the sound of a 'swarm of bees'. Muscle histopathology showed hypertrophy of all fiber types. Direct sequencing of CLCN1 revealed a mutation disrupting a donor splice site downstream of exon 16 in only the affected cats. In vitro translation of the mutated protein predicted a premature truncation and partial lack of the highly conserved CBS1 (cystathionine β-synthase) domain critical for ion transport activity and one dimerization domain pivotal in channel formation. Genetic screening of the Winnipeg random bred population of the cats' origin identified carriers of the mutation. A genetic test for population screening is now available and carrier cats from the feral population can be identified.

BACKGROUND: Nodal/TGF-Lefty signaling pathway has important effects at early stages of differentiation of human embryonic stem cells in directing them to differentiate into different embryonic lineages. LEFTY, one of transforming growth factors in the Nodal/TGF-Lefty signaling pathway, plays an important role in the development of heart. The aim of this work was to find evidence on whether Lefty variations are associated with congenital heart diseases (CHD). METHODS: We sequenced the Lefty gene for 230 Chinese Han CHD patients and evaluated SNPs rs2295418, rs360057 and g.G169A, which are located within the translated regions of the genes. The statistical analyses were conducted using Chi-Square Tests as implemented in SPSS (version 13.0). The Hardy-Weinberg equilibrium test of the population was carried out using online software OEGE, and multiple-sequence alignments of LEFTY proteins were carried out using the Vector NTI software. RESULTS: Two heterozygous variants in Lefty1 gene, g.G169A and g.A1035C, and one heterozygous variant in Lefty2 gene, g.C925A, were identified. Statistical analyses showed that the rs2295418 (g.C925A) variant in Lefty2 gene was obviously associated with the risk of CHD (P value = 0.0160.05). CONCLUSIONS: The SNP rs2295418 in the Lefty2 gene is associated with CHD in Chinese Han populations.

Endothelin-converting enzyme-1 (ECE-1) is the enzyme predominantly responsible for producing active endothelin-1 (ET-1), a mitogenic peptide implicated in the aetiology of a number of diseases, including cancer. Elevated levels of ECE-1 have been observed in a range of malignancies, with high expression conferring poor prognosis and aiding the acquisition of androgen independence in prostate cancer. The mechanisms regulating the expression of ECE-1 in cancer cells are poorly understood, hampering the development of novel therapies targeting the endothelin axis. Here we provide evidence that the expression of ECE-1 is markedly inhibited by its 3'UTR, and that alternative polyadenylation (APA) results in the production of ECE-1 transcripts with truncated 3'UTRs which promote elevated protein expression. Abolition of the ECE-1 APA sites reduced protein expression from a reporter vector in prostate cancer cells, suggesting these sites are functional. This is the first study to identify ECE-1 as a target for APA, a regulatory mechanism aberrantly activated in cancer cells, and provides novel information about the mechanisms leading to ECE-1 overexpression in malignant cells.

BackgroundMany open problems in bioinformatics involve elucidating underlying functional signals in biological sequences. DNA sequences, in particular, are characterized by rich architectures in which functional signals are increasingly found to combine local and distal interactions at the nucleotide level. Problems of interest include detection of regulatory regions, splice sites, exons, hypersensitive sites, and more. These problems naturally lend themselves to formulation as classification problems in machine learning. When classification is based on features extracted from the sequences under investigation, success is critically dependent on the chosen set of features.MethodologyWe present an algorithmic framework (EFFECT) for automated detection of functional signals in biological sequences. We focus here on classification problems involving DNA sequences which state-of-the-art work in machine learning shows to be challenging and involve complex combinations of local and distal features. EFFECT uses a two-stage process to first construct a set of candidate sequence-based features and then select a most effective subset for the classification task at hand. Both stages make heavy use of evolutionary algorithms to efficiently guide the search towards informative features capable of discriminating between sequences that contain a particular functional signal and those that do not.ResultsTo demonstrate its generality, EFFECT is applied to three separate problems of importance in DNA research: the recognition of hypersensitive sites, splice sites, and ALU sites. Comparisons with state-of-the-art algorithms show that the framework is both general and powerful. In addition, a detailed analysis of the constructed features shows that they contain valuable biological information about DNA architecture, allowing biologists and other researchers to directly inspect the features and potentially use the insights obtained to assist wet-laboratory studies on retainment or modification of a specific signal. Code, documentation, and all data for the applications presented here are provided for the community at http://www.cs.gmu.edu/~ashehu/?q=OurTools.

Mutations in the first nucleotide of exons (E(+1)) mostly affect pre-mRNA splicing when found in AG-dependent 3' splice sites, whereas AG-independent splice sites are more resistant. The AG-dependency, however, may be difficult to assess just from primary sequence data as it depends on the quality of the polypyrimidine tract. For this reason, in silico prediction tools are commonly used to score 3' splice sites. In this study, we have assessed the ability of sequence features and in silico prediction tools to discriminate between the splicing-affecting and non-affecting E(+1) variants. For this purpose, we newly tested 16 substitutions in vitro and derived other variants from literature. Surprisingly, we found that in the presence of the substituting nucleotide, the quality of the polypyrimidine tract alone was not conclusive about its splicing fate. Rather, it was the identity of the substituting nucleotide that markedly influenced it. Among the computational tools tested, the best performance was achieved using the Maximum Entropy Model and Position-Specific Scoring Matrix. As a result of this study, we have now established preliminary discriminative cut-off values showing sensitivity up to 95% and specificity up to 90%. This is expected to improve our ability to detect splicing-affecting variants in a clinical genetic setting.

PIWI family proteins have recently emerged as essential contributors in numerous biological processes including germ cell development, stem cell maintenance and epigenetic reprogramming. Expression of some of the family members has been shown to be elevated in tumors. In particular, PIWIL2 has been probed as a potential neoplasia biomarker in many cancers in humans. Previously, PIWIL2 was shown to be expressed in most tumours as a set of its shorter isoforms. In this work, we demonstrated the presence of its 60 kDa (PL2L60A) and 80 kDa (PL2L80A) isoforms in testicular cancer cell lines. We also ascertained the transcriptional boundaries of mRNAs and alternative promoter regions for these PIWIL2 isoforms. Further, we probed a range of testicular germ cell tumor (TGCT) samples and found PIWIL2 to be predominantly expressed as PL2L60A in most of them. Importantly, the levels of both PL2L60A mRNA and protein products were found to vary depending on the differentiation subtype of TGCTs, i.e., PL2L60A expression is significantly higher in undifferentiated seminomas and appears to be substantially decreased in mixed and nonseminomatous TGCTs. The higher level of PL2L60A expression in undifferentiated TGCTs was further validated in the model system of retinoic acid induced differentiation in NT2/D1 cell line. Therefore, both PL2L60A mRNA and protein abundance could serve as an additional marker distinguishing between seminomas and nonseminomatous tumors with different prognosis and therapy approaches.

BackgroundHereditary multiple exostoses (HME) is an autosomal dominant disease. The classical paradigm of mutation screening seeks to relate alterations in the exostosin glycosyltransferase genes, EXT1 and EXT2, which are responsible for over 70% of HME cases. However, the pathological significance of the majority of these mutations is often unclear.MethodsIn a Chinese family with HME, EXT1 and EXT2 genes were screened by direct sequencing. The consequence of a detected mutant was predicted by in silico analysis and confirmed by mRNA analysis. The EXT1 and EXT2 mRNA and protein levels and the HS patterns in the HME patients were compared with those in healthy controls.ResultsA heterozygous transition (c.743+1G>A) in the EXT2 gene, which co-segregated with the HME phenotype in this family, was identified. The G residue at position +1 in intron 4 of EXT2 was predicted to be a 5' donor splice site. The mRNA analysis revealed an alternative transcript with a cryptic splice site 5 bp downstream of the wild-type site, which harbored a premature stop codon. However, the predicted truncated protein was not detected by western blot analysis. Decay of the mutant mRNA was shown by clone sequencing and quantification analysis. The corresponding downregulation of the EXT2 mRNA will contribute to the abnormal EXT1/EXT2 ratio and HS pattern that were detected in the patients with HME.ConclusionThe heterozygous mutation c.743+1G>A in the EXT2 gene causes HME as a result of abnormal splicing, mRNA decay, and the resulting haploinsufficiency of EXT2.