Balanced processing of HIV-1 RNA is critical to virus replication and is regulated by host factors. In this report, we demonstrate that overexpression of either Tra2α or Tra2β results in a marked reduction in HIV-1 Gag/Env expression, an effect associated with changes in HIV-1 RNA accumulation, altered viral splice site usage, and a block to export of HIV-1 genomic RNA. A natural isoform of Tra2β (Tra2ß3), lacking the N-terminal RS domain, also suppressed HIV-1 expression but had different effects on viral RNA processing. The functional differences between the Tra2β isoforms were also observed in the context of another RNA substrate indicating that these factors have distinct functions within the cell. Finally, we demonstrate that Tra2ß depletion results in a selective reduction in HIV-1 Env expression as well as an increase in multiply spliced viral RNA. Together, the findings indicate that Tra2α/β can play important roles in regulating HIV-1 RNA metabolism and expression.

RNA-sequencing and tailored bioinformatic methodologies have paved the way for identification of expressed fusion genes from the chaotic genomes of solid tumors. We have recently successfully exploited RNA-sequencing for the discovery of 24 novel fusion genes in breast cancer. Here, we demonstrate the importance of continuous optimization of the bioinformatic methodology for this purpose, and report the discovery and experimental validation of 13 additional fusion genes from the same samples. Integration of copy number profiling with the RNA-sequencing results revealed that the majority of the gene fusions were promoter-donating events that occurred at copy number transition points or involved high-level DNA-amplifications. Sequencing of genomic fusion break points confirmed that DNA-level rearrangements underlie selected fusion transcripts. Furthermore, a significant portion (>60%) of the fusion genes were alternatively spliced. This illustrates the importance of reanalyzing sequencing data as gene definitions change and bioinformatic methods improve, and highlights the previously unforeseen isoform diversity among fusion transcripts.

mRNA translation into proteins is highly regulated, but the role of mRNA isoforms, noncoding RNAs (ncRNAs), and genetic variants remains poorly understood. mRNA levels on polysomes have been shown to correlate well with expressed protein levels, pointing to polysomal loading as a critical factor. To study regulation and genetic factors of protein translation we measured levels and allelic ratios of mRNAs and ncRNAs (including microRNAs) in lymphoblast cell lines (LCL) and in polysomal fractions. We first used targeted assays to measure polysomal loading of mRNA alleles, confirming reported genetic effects on translation of OPRM1 and NAT1, and detecting no effect of rs1045642 (3435C>T) in ABCB1 (MDR1) on polysomal loading while supporting previous results showing increased mRNA turnover of the 3435T allele. Use of high-throughput sequencing of complete transcript profiles (RNA-Seq) in three LCLs revealed significant differences in polysomal loading of individual RNA classes and isoforms. Correlated polysomal distribution between protein-coding and non-coding RNAs suggests interactions between them. Allele-selective polysome recruitment revealed strong genetic influence for multiple RNAs, attributable either to differential expression of RNA isoforms or to differential loading onto polysomes, the latter defining a direct genetic effect on translation. Genes identified by different allelic RNA ratios between cytosol and polysomes were enriched with published expression quantitative trait loci (eQTLs) affecting RNA functions, and associations with clinical phenotypes. Polysomal RNA-Seq combined with allelic ratio analysis provides a powerful approach to study polysomal RNA recruitment and regulatory variants affecting protein translation.

The dysregulation of vascular endothelial cells by thrombin has been implicated in the development of a number of pathologic disorders such as inflammatory conditions, cancer, diabetes, coronary heart disease. However, transcriptional regulation of vascular endothelial cells by thrombin is not completely understood. In the present study, Illumina RNA-seq was used to profile the transcriptome in human pulmonary microvascular endothelial cells (HMVEC-L) treated with thrombin for 6 h to gain insight into thrombin's direct effects on the endothelial function. Out of 100 million total reads from a paired end sequencing assay, 91-94% of the reads were aligned to over 16,000 genes in the reference human genome. Thrombin upregulated 150 known genes and 480 known isoforms, and downregulated 2,190 known genes and 3,574 known isoforms by at least 2 fold. Of note, thrombin upregulated 1,775 previously unknown isoforms and downregulated 12,202 previously unknown isoforms by at least 2 fold. Many genes displayed isoform specific differential expression levels and different usage of transcriptional start sites after the thrombin treatment. The cross comparisons between our RNA-seq data and those of DNA microarray analysis of either 6 h thrombin treated HUVEC or 5 h TNFα treated HMVEC have provided a significant overlapping list of differentially expressed genes, supporting the robust utility of our dataset. Further in-depth follow-up analysis of the transcriptional regulation reported in this study may shed light on molecular pathogenic mechanisms underlying thrombin mediated endothelial dysfunction in various diseases and provide new leads of potential therapeutic targets.

The shape of a neuron facilitates its functionality within neural circuits. Dendrites integrate incoming signals from axons, receiving excitatory input onto small protrusions called dendritic spines. Therefore, understanding dendritic growth and development is fundamental for discerning neural function. We previously demonstrated that EphA7 receptor signaling during cortical development impacts dendrites in two ways: EphA7 restricts dendritic growth early and promotes dendritic spine formation later. Here, the molecular basis for this shift in EphA7 function is defined. Expression analyses reveal that EphA7 full-length (EphA7-FL) and truncated (EphA7-T1; lacking kinase domain) isoforms are dynamically expressed in the developing cortex. Peak expression of EphA7-FL overlaps with dendritic elaboration around birth, while highest expression of EphA7-T1 coincides with dendritic spine formation in early postnatal life. Overexpression studies in cultured neurons demonstrate that EphA7-FL inhibits both dendritic growth and spine formation, while EphA7-T1 increases spine density. Furthermore, signaling downstream of EphA7 shifts during development, such that in vivo inhibition of mTOR by rapamycin in EphA7-mutant neurons ameliorates dendritic branching, but not dendritic spine phenotypes. Finally, direct interaction between EphA7-FL and EphA7-T1 is demonstrated in cultured cells, which results in reduction of EphA7-FL phosphorylation. In cortex, both isoforms are colocalized to synaptic fractions and both transcripts are expressed together within individual neurons, supporting a model where EphA7-T1 modulates EphA7-FL repulsive signaling during development. Thus, the divergent functions of EphA7 during cortical dendrite development are explained by the presence of two variants of the receptor.

Mutations in the X-linked Cyclin-Dependent Kinase-Like 5 gene (CDKL5) cause early onset infantile spasms and subsequent severe developmental delay in affected children. Deleterious mutations have been reported to occur throughout the CDKL5 coding region. Several studies point to a complex CDKL5 gene structure in terms of exon usage and transcript expression. Improvements in molecular diagnosis and more extensive research into the neurobiology of CDKL5 and pathophysiology of CDKL5 disorders necessitate an updated analysis of the gene. In this study, we have analysed human and mouse CDKL5 transcript patterns both bioinformatically and experimentally. We have characterised the predominant brain isoform of CDKL5, a 9.7 kb transcript comprised of 18 exons with a large 6.6 kb 3'-untranslated region (UTR), which we name hCDKL5_1. In addition we describe new exonic regions and a range of novel splice and UTR isoforms. This has enabled the description of an updated gene model in both species and a standardised nomenclature system for CDKL5 transcripts. Profiling revealed tissue- and brain development stage-specific differences in expression between transcript isoforms. These findings provide an essential backdrop for the diagnosis of CDKL5-related disorders, for investigations into the basic biology of this gene and its protein products, and for the rational design of gene-based and molecular therapies for these disorders.

WNT5A, a member of the WNT family of secreted lipid-modified glycoproteins, is a critical regulator of a host of developmental processes, including limb formation, lung morphogenesis, intestinal elongation and mammary gland development. Altered WNT5A expression has been associated with a number of cancers. Interestingly, in certain types of cancers, such as hematological malignancies and colorectal carcinoma, WNT5A is inactivated and exerts a tumor suppressive function, while in other cancers, such as melanoma and gastric carcinoma, WNT5A is overexpressed and promotes tumor progression. The mechanism by which WNT5A achieves these distinct activities in cancers is poorly understood. Here, we provide evidence that the WNT5A gene produces two protein isoforms, WNT5A-long (WNT5A-L) and WNT5A-short (WNT5A-S). Amino-terminal sequencing and a WNT5A-L specific antibody demonstrate that the mature and secreted isoforms are distinct, with WNT5A-L carrying an additional 18 N-terminal amino acids. Biochemical analysis indicates that both purified proteins are similar with respect to their stability, hydrophobicity and WNT/β-catenin signaling activity. Nonetheless, modulation of these two WNT5A isoforms, either through ectopic expression or knockdown, demonstrates that they exert distinct activities in cancer cell lines: while WNT5A-L inhibits proliferation of tumor cell lines, WNT5A-S promotes their growth. Finally, we show that expression of these two WNT5A isoforms is altered in breast and cervix carcinomas, as well as in the most aggressive neuroblastoma tumors. In these cancers, WNT5A-L is frequently down-regulated, whereas WNT5A-S is found overexpressed in a significant fraction of tumors. Altogether, our study provides evidence that the distinct activities of WNT5A in cancer can be attributed to the production of two WNT5A isoforms.

Splicing of the human vascular endothelial growth factor-A (VEGF-A) gene has been reported to generate angiogenic (VEGFxxx) and anti-angiogenic (VEGFxxxb) isoforms. Corresponding VEGFxxxb isoforms have also been reported in rat and mouse. We examined VEGFxxxb expression in mouse fibrosarcoma cell lines expressing all or individual VEGF isoforms (VEGF120, 164 or 188), grown in vitro and in vivo, and compared results with those from normal mouse and human tissues. Importantly, genetic construction of VEGF164 and VEGF188 expressing fibrosarcomas, in which exon 7 is fused to the conventional exon 8, precludes VEGFxxxb splicing from occurring. Thus, these two fibrosarcoma cell lines provided endogenous negative controls. Using RT-PCR we show that primers designed to simultaneously amplify VEGFxxx and VEGFxxxb isoforms amplified only VEGFxxx variants in both species. Moreover, only VEGFxxx species were generated when mouse podocytes were treated with TGFβ-1, a reported activator of VEGFxxxb splice selection in human podocytes. A VEGF164/120 heteroduplex species was identified as a PCR artefact, specifically in mouse. VEGFxxxb isoform-specific PCR did amplify putative VEGFxxxb species in mouse and human tissues, but unexpectedly also in VEGF188 and VEGF164 fibrosarcoma cells and tumours, where splicing to produce true VEGFxxxb isoforms cannot occur. Moreover, these products were only consistently generated using reverse primers spanning more than 5 bases across the 8b/7 or 8b/5 splice junctions. Primer annealing to VEGFxxx transcripts and amplification of exon 8b primer 'tails' explained the artefactual generation of VEGFxxxb products, since the same products were generated when the PCR reactions were performed with cDNA from VEGF164/VEGF188 'knock-in' vectors used in the generation of single VEGF isoform-expressing transgenic mice from which the fibrosarcoma lines were developed. Collectively, our results highlight important pitfalls in data interpretation associated with detecting VEGFxxxb isoforms using current methods, and demonstrate that anti-angiogenic isoforms are not commonly expressed in mouse or human tissues.

Gamma-aminobutyric acid type B (GABAB) receptors decrease neural activity through G protein signaling. There are two subunits, GABAB1 and GABAB2. Alternative splicing provides GABAB1 with structural and functional diversity. cDNA microarrays showed strong signals from human brain RNA using GABAB1 intron 4 region probes. Therefore, we predicted the existence of novel splice variants.

Phenotypic differences between species are driven by changes in gene expression and, by extension, by modifications in the regulation of the transcriptome. Investigation of mammalian transcriptome divergence has been restricted to analysis of bulk gene expression levels and gene-internal splicing. Using allele-specific expression analysis in inter-strain hybrids of Mus musculus, we determined the contribution of multiple cellular regulatory systems to transcriptome divergence, including: alternative promoter usage, transcription start site selection, cassette exon usage, alternative last exon usage, and alternative polyadenylation site choice. Between mouse strains, a fifth of genes have variations in isoform usage that contribute to transcriptomic changes, half of which alter encoded amino acid sequence. Virtually all divergence in isoform usage altered the post-transcriptional regulatory instructions in gene UTRs. Furthermore, most genes with isoform differences between strains contain changes originating from multiple regulatory systems. This result indicates widespread cross-talk and coordination exists among different regulatory systems. Overall, isoform usage diverges in parallel with and independently to gene expression evolution, and the cis and trans regulatory contribution to each differs significantly.

Rapid non-genomic effects of 17β-estradiol are elicited by the activation of different estrogen receptor-α isoforms. Presence of surface binding sites for estrogen have been identified in cells transfected with full-length estrogen receptor-α66 (ER66) and the truncated isoforms, estrogen receptor-α46 (ER46) and estrogen receptor-α36 (ER36). However, the binding affinities of the membrane estrogen receptors (mERs) remain unknown due to the difficulty of developing of stable mER-transfected cell lines with sufficient mER density, which has largely hampered biochemical binding studies. The present study utilized cell-free expression systems to determine the binding affinities of 17β-estradiol to mERs, and the relationship among palmitoylation, membrane insertion and binding affinities. Saturation binding assays of human mERs revealed that [³H]-17β-estradiol bound ER66 and ER46 with Kd values of 68.81 and 60.72 pM, respectively, whereas ER36 displayed no specific binding within the tested concentration range. Inhibition of palmitoylation or removal of the nanolipoprotein particles, used as membrane substitute, reduced the binding affinities of ER66 and ER46 to 17β-estradiol. Moreover, ER66 and ER46 bound differentially with some estrogen receptor agonists and antagonists, and phytoestrogens. In particular, the classical estrogen receptor antagonist, ICI 182,780, had a higher affinity for ER66 than ER46. In summary, the present study defines the binding affinities for human estrogen receptor-α isoforms, and demonstrates that ER66 and ER46 show characteristics of mERs. The present data also indicates that palmitoylation and membrane insertion of mERs are important for proper receptor conformation allowing 17β-estradiol binding. The differential binding of ER66 and ER46 with certain compounds substantiates the prospect of developing mER-selective drugs.

The 3',5'-cyclic adenosine monophosphate (cAMP)-dependent protein kinase, or protein kinase A (PKA), pathway is one of the most versatile and best studied signaling pathways in eukaryotic cells. The two paralogous PKA catalytic subunits Cα and Cβ, encoded by the genes PRKACA and PRKACB, respectively, are among the best understood model kinases in signal transduction research. In this work, we explore and elucidate the evolution of the alternative 5' exons and the splicing pattern giving rise to the numerous PKA catalytic subunit isoforms. In addition to the universally conserved Cα1/Cβ1 isoforms, we find kinase variants with short N-termini in all main vertebrate classes, including the sperm-specific Cα2 isoform found to be conserved in all mammals. We also describe, for the first time, a PKA Cα isoform with a long N-terminus, paralogous to the PKA Cβ2 N-terminus. An analysis of isoform-specific variation highlights residues and motifs that are likely to be of functional importance.

Alternative splicing (AS) is an intrinsic regulatory mechanism of all metazoans. Recent findings suggest that 100% of multiexonic human genes give rise to splice isoforms. AS can be specific to tissue type, environment or developmentally regulated. Splice variants have also been implicated in various diseases including cancer. Detection of these variants will enhance our understanding of the complexity of the human genome and provide disease-specific and prognostic biomarkers. We adopted a proteomics approach to identify exon skip events - the most common form of AS. We constructed a database harboring the peptide sequences derived from all hypothetical exon skip junctions in the human genome. Searching tandem mass spectrometry (MS/MS) data against the database allows the detection of exon skip events, directly at the protein level. Here we describe the application of this approach to human platelets, including the mRNA-based verification of novel splice isoforms of ITGA2, NPEPPS and FH. This methodology is applicable to all new or existing MS/MS datasets.

Human BST-2 (hBST-2) has been identified as a cellular antiviral factor that blocks the release of various enveloped viruses. Orthologues of BST-2 have been identified in several species, including human, monkeys, pig, mouse, cat and sheep. All have been reported to possess antiviral activity. Duplication of the BST-2 gene has been observed in sheep and the paralogues are referred to as ovine BST-2A and BST2-B, although only a single gene corresponding to BST-2 has been identified in most species. In this study, we identified three isoforms of bovine BST-2, named bBST-2A1, bBST-2A2 and bBST-2B, in bovine cells treated with type I interferon, but not in untreated cells. Both bBST-2A1 and bBST-2A2 are posttranslationally modified by N-linked glycosylation and a GPI-anchor as well as hBST-2, while bBST-2B has neither of these modifications. Exogenous expression of bBST-2A1 or bBST-2A2 markedly reduced the production of bovine leukemia virus and vesicular stomatitis virus from cells, while the antiviral activity of bBST-2B was much weaker than those of bBST-2A1 and bBST-2A2. Our data suggest that bBST-2A1 and bBST-2A2 function as part of IFN-induced innate immunity against virus infection. On the other hand, bBST-2B may have a different physiological function from bBST-2A1 and bBST-2A2.

The chicken has long served as an important model organism in many fields, and continues to aid our understanding of animal development. Functional genomics studies aimed at probing the mechanisms that regulate development require high-quality genomes and transcript annotations. The quality of these resources has improved dramatically over the last several years, but many isoforms and genes have yet to be identified. We hope to contribute to the process of improving these resources with the data presented here: a set of long cDNA sequencing reads, and a curated set of new genes and transcript isoforms not currently represented in the most up-to-date genome annotation currently available to the community of researchers who rely on the chicken genome.

The vasa gene encodes an ATP-dependent RNA helicase of the DEAD box protein family that functions in a broad range of molecular events involving duplex RNA. In most species, the germline specific expression of vasa becomes a molecular marker widely used in the visualization and labeling of primordial germ cells (PGCs) and a tool in surrogate broodstock production through PGC transplantation. The vasa gene from tongue sole (Cynoglossus semilaevis) was characterized to promote the development of genetic breeding techniques in this species. Three C. semilaevis vasa transcripts were isolated, namely vas-l, vas-m, and vas-s. Quantitative real-time PCR results showed that C. semilaevis vasa transcripts were prevalently expressed in gonads, with very weak expression of vas-s in other tissues. Embryonic development expression profiles revealed the onset of zygotic transcription of vasa mRNAs and the maternal deposit of the three transcripts. The genetic ZW female juvenile fish was discriminated from genetic ZZ males by a pair of female specific primers. Only the expression of vas-s can be observed in both sexes during early gonadal differentiation. Before PGCs started mitosis, there was sexually dimorphic expression of vas-s with the ovary showing higher levels and downward trend. The results demonstrated the benefits of vasa as a germline specific marker for PGCs during embryonic development and gonadal differentiation. This study lays the groundwork for further application of C. semilaevis PGCs in fish breeding.

A unified proteochemometric (PCM) model for the prediction of the ability of drug-like chemicals to inhibit five major drug metabolizing CYP isoforms (i.e. CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4) was created and made publicly available under the Bioclipse Decision Support open source system at www.cyp450model.org. In regards to the proteochemometric modeling we represented the chemical compounds by molecular signature descriptors and the CYP-isoforms by alignment-independent description of composition and transition of amino acid properties of their protein primary sequences. The entire training dataset contained 63 391 interactions and the best PCM model was obtained using signature descriptors of height 1, 2 and 3 and inducing the model with a support vector machine. The model showed excellent predictive ability with internal AUC = 0.923 and an external AUC = 0.940, as evaluated on a large external dataset. The advantage of PCM models is their extensibility making it possible to extend our model for new CYP isoforms and polymorphic CYP forms. A key benefit of PCM is that all proteins are confined in one single model, which makes it generally more stable and predictive as compared with single target models. The inclusion of the model in Bioclipse Decision Support makes it possible to make virtual instantaneous predictions (∼100 ms per prediction) while interactively drawing or modifying chemical structures in the Bioclipse chemical structure editor.

Friedreich ataxia (FRDA) is an inherited neurodegenerative disease caused by frataxin (FXN) deficiency. The nervous system and heart are the most severely affected tissues. However, highly mitochondria-dependent tissues, such as kidney and liver, are not obviously affected, although the abundance of FXN is normally high in these tissues. In this study we have revealed two novel FXN isoforms (II and III), which are specifically expressed in affected cerebellum and heart tissues, respectively, and are functional in vitro and in vivo. Increasing the abundance of the heart-specific isoform III significantly increased the mitochondrial aconitase activity, while over-expression of the cerebellum-specific isoform II protected against oxidative damage of Fe-S cluster-containing aconitase. Further, we observed that the protein level of isoform III decreased in FRDA patient heart, while the mRNA level of isoform II decreased more in FRDA patient cerebellum compared to total FXN mRNA. Our novel findings are highly relevant to understanding the mechanism of tissue-specific pathology in FRDA.

The primary transcript of the mammalian Fragile X Mental Retardation-1 gene (Fmr1), like many transcripts in the central nervous system, is alternatively spliced to yield mRNAs encoding multiple proteins, which can possess quite different biochemical properties. Despite the fact that the relative levels of the 12 Fmr1 transcript isoforms examined here vary by as much as two orders of magnitude amongst themselves in both adult and embryonic mouse brain, all are associated with polyribosomes, consistent with translation into the corresponding isoforms of the protein product, FMRP (Fragile X Mental Retardation Protein). Employing the RiboTag methodology developed in our laboratory, the relative proportions of the 7 most abundant transcript isoforms were measured specifically in neurons and found to be similar to those identified in whole brain. Measurements of isoform profiles across 11 regions of adult brain yielded similar distributions, with the exceptions of the hippocampus and the olfactory bulb. These two regions differ from most of the brain in relative amounts of transcripts encoding an alternate form of one of the KH RNA binding domains. A possible relationship between patterns of expression in the hippocampus and olfactory bulb and the presence of neuroblasts in these two regions is suggested by the isoform patterns in early embryonic brain and in cultured neural progenitor cells. These results demonstrate that the relative levels of the Fmr1 isoforms are modulated according to developmental stage, highlighting the complex ramifications of losing all the protein isoforms in individuals with Fragile X Syndrome. It should also be noted that, of the eight most prominent FMRP isoforms (1-3, 6-9 and 12) in mouse, only two have the major site of phosphorylation at Ser-499, which is thought to be involved in some of the regulatory interactions of this protein.

It has been shown that imprecise cleavage of a primary or precursor RNA by Drosha or Dicer, respectively, may yield a group of microRNA (miRNA) variants designated as "isomiR". Variations in the relative abundance of isoforms for a given miRNA among different species and different cell types beg the question whether these isomiRs might regulate target genes differentially. We compared the capacity of three miR-31 isoforms (miR-31-H, miR-31-P, and miR-31-M), which differ only slightly in their 5'- and/or 3'-end sequences, to regulate several known targets and a predicted target, Dicer. Notably, we found isomiR-31s displayed concordant and discordant regulation of 6 known target genes. Furthermore, we validated a predicted target gene, Dicer, to be a novel target of miR-31 but only miR-31-P could directly repress Dicer expression in both MCF-7 breast cancer cells and A549 lung cancer cells, resulting in their enhanced sensitivity to cisplatin, a known attribute of Dicer knockdown. This was further supported by reporter assay using full length 3'-untranslated region (UTR) of Dicer. Our findings not only revealed Dicer to be a direct target of miR-31, but also demonstrated that isomiRs displayed similar and disparate regulation of target genes in cell-based systems. Coupled with the variations in the distribution of isomiRs among different cells or conditions, our findings support the possibility of fine-tuning gene expression by miRNAs.