The IRE1α/XBP1s Pathway Is Essential for the Glucose Response and Protection of ß Cells.

Although glucose uniquely stimulates proinsulin biosynthesis in ß cells, surprisingly little is known of the underlying mechanism(s). Here, we demonstrate that glucose activates the unfolded protein response transducer inositol-requiring enzyme 1 alpha (IRE1α) to initiate X-box-binding protein 1 (Xbp1) mRNA splicing in adult primary ß cells. Using mRNA sequencing (mRNA-Seq), we show that unconventional Xbp1 mRNA splicing is required to increase and decrease the expression of several hundred mRNAs encoding functions that expand the protein secretory capacity for increased insulin production and protect from oxidative damage, respectively. At 2 wk after tamoxifen-mediated Ire1α deletion, mice develop hyperglycemia and hypoinsulinemia, due to defective ß cell function that was exacerbated upon feeding and glucose stimulation. Although previous reports suggest IRE1α degrades insulin mRNAs, Ire1α deletion did not alter insulin mRNA expression either in the presence or absence of glucose stimulation. Instead, ß cell failure upon Ire1α deletion was primarily due to reduced proinsulin mRNA translation primarily because of defective glucose-stimulated induction of a dozen genes required for the signal recognition particle (SRP), SRP receptors, the translocon, the signal peptidase complex, and over 100 other genes with many other intracellular functions. In contrast, Ire1α deletion in ß cells increased the expression of over 300 mRNAs encoding functions that cause inflammation and oxidative stress, yet only a few of these accumulated during high glucose. Antioxidant treatment significantly reduced glucose intolerance and markers of inflammation and oxidative stress in mice with ß cell-specific Ire1α deletion. The results demonstrate that glucose activates IRE1α-mediated Xbp1 splicing to expand the secretory capacity of the ß cell for increased proinsulin synthesis and to limit oxidative stress that leads to ß cell failure.

LINE-1 Mediated Insertion into Poc1a (Protein of Centriole 1 A) Causes Growth Insufficiency and Male Infertility in Mice.

Skeletal dysplasias are a common, genetically heterogeneous cause of short stature that can result from disruptions in many cellular processes. We report the identification of the lesion responsible for skeletal dysplasia and male infertility in the spontaneous, recessive mouse mutant chagun. We determined that Poc1a, encoding protein of the centriole 1a, is disrupted by the insertion of a processed Cenpw cDNA, which is flanked by target site duplications, suggestive of a LINE-1 retrotransposon-mediated event. Mutant fibroblasts have impaired cilia formation and multipolar spindles. Male infertility is caused by defective spermatogenesis early in meiosis and progressive germ cell loss. Spermatogonial stem cell transplantation studies revealed that Poc1a is essential for normal function of both Sertoli cells and germ cells. The proliferative zone of the growth plate is small and disorganized because chondrocytes fail to re-align after cell division and undergo increased apoptosis. Poc1a and several other genes associated with centrosome function can affect the skeleton and lead to skeletal dysplasias and primordial dwarfisms. This mouse mutant reveals how centrosome dysfunction contributes to defects in skeletal growth and male infertility.

H2A.B facilitates transcription elongation at methylated CpG loci.

H2A.B is a unique histone H2A variant that only exists in mammals. Here we found that H2A.B is ubiquitously expressed in major organs. Genome-wide analysis of H2A.B in mouse ES cells shows that H2A.B is associated with methylated DNA in gene body regions. Moreover, H2A.B-enriched gene loci are actively transcribed. One typical example is that H2A.B is enriched in a set of differentially methylated regions at imprinted loci and facilitates transcription elongation. These results suggest that H2A.B positively regulates transcription elongation by overcoming DNA methylation in the transcribed region. It provides a novel mechanism by which transcription is regulated at DNA hypermethylated regions.

CD1d-unrestricted NKT cells are endowed with a hybrid function far superior than that of iNKT cells.

Invariant natural killer T (iNKT) cells to date represent the best example of cells known to have a hybrid function, representing both innate and adaptive immunity. Shared phenotypic similarities with NK cells together with a rapid response to a cytokine stimulus and a productive TCR engagement are the features that underline the hybrid nature of iNKT cells. Using these criteria, we provide molecular and functional evidence demonstrating that CD1d-independent (CD1d(ind)) NKT cells, a population of CD1d-unrestricted NKT cells, are endowed with a hybrid function far superior to that of iNKT cells: (i) an extensive shared program with NK cells, (ii) a closer Euclidian distance with NK cells, and (iii) the ability to respond to innate stimuli (Poly:IC) with cytotoxic potential in the same manner as NK cells identify a hybrid feature in CD1d(ind)NKT cells that truly fulfills the dual function of an NK and a T cell. Our finding that CD1d(ind)NKT cells are programmed to act like NK cells in response to innate signals while being capable of adaptive responses is unprecedented, and thus might reemphasize CD1d-unrestricted NKT cells as a subset of lymphocytes that could affect biological processes of antimicrobial and tumor immunity in a unique way.

Ecological and genetic interactions between cyanobacteria and viruses in a low-oxygen mat community inferred through metagenomics and metatranscriptomics.

Metagenomic and metatranscriptomic sequencing was conducted on cyanobacterial mats of the Middle Island Sinkhole (MIS), Lake Huron. Metagenomic data from 14 samples collected over 5 years were used to reconstruct genomes of two genotypes of a novel virus, designated PhV1 type A and PhV1 type B. Both viral genotypes encode and express nblA, a gene involved in degrading phycobilisomes, which are complexes of pigmented proteins that harvest light for photosynthesis. Phylogenetic analysis indicated that the viral-encoded nblA is derived from the host cyanobacterium, Phormidium MIS-PhA. The cyanobacterial host also has two complete CRISPR (clustered regularly interspaced short palindromic repeats) systems that serve as defence mechanisms for bacteria and archaea against viruses and plasmids. One 45 bp CRISPR spacer from Phormidium had 100% nucleotide identity to PhV1 type B, but this region was absent from PhV1 type A. Transcripts from PhV1 and the Phormidium CRISPR loci were detected in all six metatranscriptomic data sets (three during the day and three at night), indicating that both are transcriptionally active in the environment. These results reveal ecological and genetic interactions between viruses and cyanobacteria at MIS, highlighting the value of parallel analysis of viruses and hosts in understanding ecological interactions in natural communities.

Next-generation sequencing identifies the Danforth's short tail mouse mutation as a retrotransposon insertion affecting Ptf1a expression.

The semidominant Danforth's short tail (Sd) mutation arose spontaneously in the 1920s. The homozygous Sd phenotype includes severe malformations of the axial skeleton with an absent tail, kidney agenesis, anal atresia, and persistent cloaca. The Sd mutant phenotype mirrors features seen in human caudal malformation syndromes including urorectal septum malformation, caudal regression, VACTERL association, and persistent cloaca. The Sd mutation was previously mapped to a 0.9 cM region on mouse chromosome 2qA3. We performed Sanger sequencing of exons and intron/exon boundaries mapping to the Sd critical region and did not identify any mutations. We then performed DNA enrichment/capture followed by next-generation sequencing (NGS) of the critical genomic region. Standard bioinformatic analysis of paired-end sequence data did not reveal any causative mutations. Interrogation of reads that had been discarded because only a single end mapped correctly to the Sd locus identified an early transposon (ETn) retroviral insertion at the Sd locus, located 12.5 kb upstream of the Ptf1a gene. We show that Ptf1a expression is significantly upregulated in Sd mutant embryos at E9.5. The identification of the Sd mutation will lead to improved understanding of the developmental pathways that are misregulated in human caudal malformation syndromes.

Regulation of the human endogenous retrovirus K (HML-2) transcriptome by the HIV-1 Tat protein.

UNLABELLED: Approximately 8% of the human genome is made up of endogenous retroviral sequences. As the HIV-1 Tat protein activates the overall expression of the human endogenous retrovirus type K (HERV-K) (HML-2), we used next-generation sequencing to determine which of the 91 currently annotated HERV-K (HML-2) proviruses are regulated by Tat. Transcriptome sequencing of total RNA isolated from Tat- and vehicle-treated peripheral blood lymphocytes from a healthy donor showed that Tat significantly activates expression of 26 unique HERV-K (HML-2) proviruses, silences 12, and does not significantly alter the expression of the remaining proviruses. Quantitative reverse transcription-PCR validation of the sequencing data was performed on Tat-treated PBLs of seven donors using provirus-specific primers and corroborated the results with a substantial degree of quantitative similarity. IMPORTANCE: The expression of HERV-K (HML-2) is tightly regulated but becomes markedly increased following infection with HIV-1, in part due to the HIV-1 Tat protein. The findings reported here demonstrate the complexity of the genome-wide regulation of HERV-K (HML-2) expression by Tat. This work also demonstrates that although HERV-K (HML-2) proviruses in the human genome are highly similar in terms of DNA sequence, modulation of the expression of specific proviruses in a given biological situation can be ascertained using next-generation sequencing and bioinformatics analysis.

Decade-long bacterial community dynamics in cystic fibrosis airways.

The structure and dynamics of bacterial communities in the airways of persons with cystic fibrosis (CF) remain largely unknown. We characterized the bacterial communities in 126 sputum samples representing serial collections spanning 8-9 y from six age-matched male CF patients. Sputum DNA was analyzed by bar-coded pyrosequencing of the V3-V5 hypervariable region of the 16S rRNA gene, defining 662 operational taxonomic units (OTUs) from >633,000 sequences. Bacterial community diversity decreased significantly over time in patients with typically progressive lung disease but remained relatively stable in patients with a mild lung disease phenotype. Antibiotic use, rather than patient age or lung function, was the primary driver of decreasing diversity. Interpatient variability in community structure exceeded intrapatient variability in serial samples. Antibiotic treatment was associated with pronounced shifts in community structure, but communities showed both short- and long-term resilience after antibiotic perturbation. There was a positive correlation between OTU occurrence and relative abundance, with a small number of persistent OTUs accounting for the greatest abundance. Significant changes in community structure, diversity, or total bacterial density at the time of pulmonary exacerbation were not observed. Despite decreasing community diversity in patients with progressive disease, total bacterial density remained relatively stable over time. These findings show the critical relationship between airway bacterial community structure, disease stage, and clinical state at the time of sample collection. These features are the key parameters with which to assess the complex ecology of the CF airway.

Disruption of RAB40AL function leads to Martin--Probst syndrome, a rare X-linked multisystem neurodevelopmental human disorder.

BACKGROUND AND AIM: Martin--Probst syndrome (MPS) is a rare X-linked disorder characterised by deafness, cognitive impairment, short stature and distinct craniofacial dysmorphisms, among other features. The authors sought to identify the causative mutation for MPS. METHODS AND RESULTS: Massively parallel sequencing in two affected, related male subjects with MPS identified a RAB40AL (also called RLGP) missense mutation (chrX:102,079,078-102,079,079AC→GA p.D59G; hg18). RAB40AL encodes a small Ras-like GTPase protein with one suppressor of cytokine signalling box. The p.D59G variant is located in a highly conserved region of the GTPase domain between ß-2 and ß-3 strands. Using RT-PCR, the authors show that RAB40AL is expressed in human fetal and adult brain and kidney, and adult lung, heart, liver and skeletal muscle. RAB40AL appears to be a primate innovation, with no orthologues found in mouse, Xenopus or zebrafish. Western analysis and fluorescence microscopy of GFP-tagged RAB40AL constructs from transiently transfected COS7 cells show that the D59G missense change renders RAB40AL unstable and disrupts its cytoplasmic localisation. CONCLUSIONS: This is the first study to show that mutation of RAB40AL is associated with a human disorder. Identification of RAB40AL as the gene mutated in MPS allows for further investigations into the molecular mechanism(s) of RAB40AL and its roles in diverse processes such as cognition, hearing and skeletal development.

Positional cloning of the combined hyperlipidemia gene Hyplip1.

Familial combined hyperlipidemia (FCHL, MIM-144250) is a common, multifactorial and heterogeneous dyslipidemia predisposing to premature coronary artery disease and characterized by elevated plasma triglycerides, cholesterol, or both. We identified a mutant mouse strain, HcB-19/Dem (HcB-19), that shares features with FCHL, including hypertriglyceridemia, hypercholesterolemia, elevated plasma apolipoprotein B and increased secretion of triglyceride-rich lipoproteins. The hyperlipidemia results from spontaneous mutation at a locus, Hyplip1, on distal mouse chromosome 3 in a region syntenic with a 1q21-q23 FCHL locus identified in Finnish, German, Chinese and US families. We fine-mapped Hyplip1 to roughly 160 kb, constructed a BAC contig and sequenced overlapping BACs to identify 13 candidate genes. We found substantially decreased mRNA expression for thioredoxin interacting protein (Txnip). Sequencing of the critical region revealed a Txnip nonsense mutation in HcB-19 that is absent in its normolipidemic parental strains. Txnip encodes a cytoplasmic protein that binds and inhibits thioredoxin, a major regulator of cellular redox state. The mutant mice have decreased CO2 production but increased ketone body synthesis, suggesting that altered redox status down-regulates the citric-acid cycle, sparing fatty acids for triglyceride and ketone body production. These results reveal a new pathway of potential clinical significance that contributes to plasma lipid metabolism.

Convergence of genetic influences in comorbidity.

BACKGROUND: Predisposition to complex diseases is explained in part by genetic variation, and complex diseases are frequently comorbid, consistent with pleiotropic genetic variation influencing comorbidity. Genome Wide Association (GWA) studies typically assess association between SNPs and a single-disease phenotype. Fisher meta-analysis combines evidence of association from single-disease GWA studies, assuming that each study is an independent test of the same hypothesis. The Rank Product (RP) method overcomes limitations posed by Fisher assumptions, though RP was not designed for GWA data. METHODS: We modified RP to accommodate GWA data, and we call it modRP. Using p-values output from GWA studies, we aggregate evidence for association between SNPs and related phenotypes. To assess significance, RP randomly samples the observed ranks to develop the null distribution of the RP statistic, and then places the observed RPs into the null distribution. ModRP eliminates the effect of linkage disequilibrium and controls for differences in power at tested SNPs, to meet RP assumptions in application to GWA data. RESULTS: After validating modRP based on both positive and negative control studies, we searched for pleiotropic influences on comorbid substance use disorders in a novel study, and found two SNPs to be significantly associated with comorbid cocaine, opium, and nicotine dependence. Placing these SNPs into biological context, we developed a protein network modeling the interaction of cocaine, nicotine, and opium with these variants. CONCLUSIONS: ModRP is a novel approach to identifying pleiotropic genetic influences on comorbid complex diseases. It can be used to assess association for related phenotypes where raw data is unavailable or inappropriate for analysis using other approaches. The method is conceptually simple and produces statistically significant, biologically relevant results.

Metscape: a Cytoscape plug-in for visualizing and interpreting metabolomic data in the context of human metabolic networks.

SUMMARY: Metscape is a plug-in for Cytoscape, used to visualize and interpret metabolomic data in the context of human metabolic networks. We have developed a metabolite database by extracting and integrating information from several public sources. By querying this database, Metscape allows users to trace the connections between metabolites and genes, visualize compound networks and display compound structures as well as information for reactions, enzymes, genes and pathways. Applying the pathway filter, users can create subnetworks that consist of compounds and reactions from a given pathway. Metscape allows users to upload experimental data, and visualize and explore compound networks over time, or experimental conditions. Color and size of the nodes are used to visualize these dynamic changes. Metscape can display the entire metabolic network or any of the pathway-specific networks that exist in the database. AVAILABILITY: Metscape can be installed from within Cytoscape 2.6.x under 'Network and Attribute I/O' category. For more information, please visit http://metscape.ncibi.org/tryplugin.html.

ConceptGen: a gene set enrichment and gene set relation mapping tool.

MOTIVATION: The elucidation of biological concepts enriched with differentially expressed genes has become an integral part of the analysis and interpretation of genomic data. Of additional importance is the ability to explore networks of relationships among previously defined biological concepts from diverse information sources, and to explore results visually from multiple perspectives. Accomplishing these tasks requires a unified framework for agglomeration of data from various genomic resources, novel visualizations, and user functionality. RESULTS: We have developed ConceptGen, a web-based gene set enrichment and gene set relation mapping tool that is streamlined and simple to use. ConceptGen offers over 20,000 concepts comprising 14 different types of biological knowledge, including data not currently available in any other gene set enrichment or gene set relation mapping tool. We demonstrate the functionalities of ConceptGen using gene expression data modeling TGF-beta-induced epithelial-mesenchymal transition and metabolomics data comparing metastatic versus localized prostate cancers.

Genome-based characterization of two prenylation steps in the assembly of the stephacidin and notoamide anticancer agents in a marine-derived Aspergillus sp.

Stephacidin and notoamide natural products belong to a group of prenylated indole alkaloids containing a core bicyclo[2.2.2]diazaoctane ring system. These bioactive fungal secondary metabolites have a range of unusual structural and stereochemical features but their biosynthesis has remained uncharacterized. Herein, we report the first biosynthetic gene cluster for this class of fungal alkaloids based on whole genome sequencing of a marine-derived Aspergillus sp. Two central pathway enzymes catalyzing both normal and reverse prenyltransfer reactions were characterized in detail. Our results establish the early steps for creation of the prenylated indole alkaloid structure and suggest a scheme for the biosynthesis of stephacidin and notoamide metabolites. The work provides the first genetic and biochemical insights for understanding the structural diversity of this important family of fungal alkaloids.

Modeling complex genetic and environmental influences on comorbid bipolar disorder with tobacco use disorder.

BACKGROUND: Comorbidity of psychiatric and substance use disorders represents a significant complication in the clinical course of both disorders. Bipolar Disorder (BD) is a psychiatric disorder characterized by severe mood swings, ranging from mania to depression, and up to a 70% rate of comorbid Tobacco Use Disorder (TUD). We found epidemiological evidence consistent with a common underlying etiology for BD and TUD, as well as evidence of both genetic and environmental influences on BD and TUD. Therefore, we hypothesized a common underlying genetic etiology, interacting with nicotine exposure, influencing susceptibility to both BD and TUD. METHODS: Using meta-analysis, we compared TUD rates for BD patients and the general population. We identified candidate genes showing statistically significant, replicated, evidence of association with both BD and TUD. We assessed commonality among these candidate genes and hypothesized broader, multi-gene network influences on the comorbidity. Using Fisher Exact tests we tested our hypothesized genetic networks for association with the comorbidity, then compared the inferences drawn with those derived from the commonality assessment. Finally, we prioritized candidate SNPs for validation. RESULTS: We estimate risk for TUD among BD patients at 2.4 times that of the general population. We found three candidate genes associated with both BD and TUD (COMT, SLC6A3, and SLC6A4) and commonality analysis suggests that these genes interact in predisposing psychiatric and substance use disorders. We identified a 69 gene network that influences neurotransmitter signaling and shows significant over-representation of genes associated with BD and TUD, as well as genes differentially expressed with exposure to tobacco smoke. Twenty four of these genes are known drug targets. CONCLUSIONS: This work highlights novel bioinformatics resources and demonstrates the effectiveness of using an integrated bioinformatics approach to improve our understanding of complex disease etiology. We illustrate the development and testing of hypotheses for a comorbidity predisposed by both genetic and environmental influences. Consistent with our hypothesis, the selected network models multiple interacting genetic influences on comorbid BD with TUD, as well as the environmental influence of nicotine. This network nominates candidate genes for validation and drug testing, and we offer a panel of SNPs prioritized for follow-up.

Genome-Wide Identification of Klebsiella pneumoniae Fitness Genes during Lung Infection.

UNLABELLED: Klebsiella pneumoniae is an urgent public health threat because of resistance to carbapenems, antibiotics of last resort against Gram-negative bacterial infections. Despite the fact that K. pneumoniae is a leading cause of pneumonia in hospitalized patients, the bacterial factors required to cause disease are poorly understood. Insertion site sequencing combines transposon mutagenesis with high-throughput sequencing to simultaneously screen thousands of insertion mutants for fitness defects during infection. Using the recently sequenced K. pneumoniae strain KPPR1 in a well-established mouse model of pneumonia, insertion site sequencing was performed on a pool of >25,000 transposon mutants. The relative fitness requirement of each gene was ranked based on the ratio of lung to inoculum read counts and concordance between insertions in the same gene. This analysis revealed over 300 mutants with at least a 2-fold fitness defect and 69 with defects ranging from 10- to >2,000-fold. Construction of 6 isogenic mutants for use in competitive infections with the wild type confirmed their requirement for lung fitness. Critical fitness genes included those for the synthesis of branched-chain and aromatic amino acids that are essential in mice and humans, the transcriptional elongation factor RfaH, and the copper efflux pump CopA. The majority of fitness genes were conserved among reference strains representative of diverse pathotypes. These results indicate that regulation of outer membrane components and synthesis of amino acids that are essential to its host are critical for K. pneumoniae fitness in the lung. IMPORTANCE: Klebsiella pneumoniae is a bacterium that commonly causes pneumonia in patients after they are admitted to the hospital. K. pneumoniae is becoming resistant to all available antibiotics, and when these infections spread to the bloodstream, over half of patients die. Since currently available antibiotics are failing, we must discover new ways to treat these infections. In this study, we asked what genes the bacterium needs to cause an infection, since the proteins encoded by these genes could be targets for new antibiotics. We identified over 300 genes that K. pneumoniae requires to grow in a mouse model of pneumonia. Many of the genes that we identified are found in K. pneumoniae isolates from throughout the world, including antibiotic-resistant forms. If new antibiotics could be made against the proteins that these genes encode, they may be broadly effective against K. pneumoniae.

Alternative splicing of human NT5E in cirrhosis and hepatocellular carcinoma produces a negative regulator of ecto-5'-nucleotidase (CD73).

Ecto-5'-nucleotidase (CD73), encoded by NT5E, is the major enzymatic source of extracellular adenosine. CD73 controls numerous pathophysiological responses and is a potential disease target, but its regulation is poorly understood. We examined NT5E regulation by alternative splicing. Genomic database analysis of human transcripts led us to identify NT5E-2, a novel splice variant that was expressed at low abundance in normal human tissues but was significantly up-regulated in cirrhosis and hepatocellular carcinoma (HCC). NT5E-2 encodes a shorter CD73 isoform we named CD73S. The presence of CD73S protein, which lacks 50 amino acids, was detected in HCC using an isoform-specific antibody. A noncanonical mouse mRNA, similar to human CD73S, was observed, but the corresponding protein was undetectable. The two human isoforms exhibited functional differences, such that ectopic expression of canonical CD73 (CD73L) in human HepG2 cells was associated with decreased expression of the proliferation marker Ki67, whereas CD73S expression did not have an effect on Ki67 expression. CD73S was glycosylated, catalytically inactive, unable to dimerize, and complexed intracellularly with the endoplasmic reticulum chaperone calnexin. Furthermore, CD73S complexed with CD73L and promoted proteasome-dependent CD73L degradation. The findings reveal species-specific CD73 regulation, with potential significance to cancer, fibrosis, and other diseases characterized by changes in CD73 expression and function.

Complete Genome Sequence of Klebsiella pneumoniae Strain ATCC 43816 KPPR1, a Rifampin-Resistant Mutant Commonly Used in Animal, Genetic, and Molecular Biology Studies.

Klebsiella pneumoniae is an urgent public health threat due to the spread of carbapenem-resistant strains causing serious, and frequently fatal, infections. To facilitate genetic, molecular, and immunological studies of this pathogen, we report the complete chromosomal sequence of a genetically tractable, prototypical strain used in animal models.

A bioinformatics approach reveals novel interactions of the OVOL transcription factors in the regulation of epithelial - mesenchymal cell reprogramming and cancer progression.

BACKGROUND: Mesenchymal to Epithelial Transition (MET) plasticity is critical to cancer progression, and we recently showed that the OVOL transcription factors (TFs) are critical regulators of MET. Results of that work also posed the hypothesis that the OVOLs impact MET in a range of cancers. We now test this hypothesis by developing a model, OVOL Induced MET (OI-MET), and sub-model (OI-MET-TF), to characterize differential gene expression in MET common to prostate cancer (PC) and breast cancer (BC). RESULTS: In the OI-MET model, we identified 739 genes differentially expressed in both the PC and BC models. For this gene set, we found significant enrichment of annotation for BC, PC, cancer, and MET, as well as regulation of gene expression by AP1, STAT1, STAT3, and NFKB1. Focusing on the target genes for these four TFs plus the OVOLs, we produced the OI-MET-TF sub-model, which shows even greater enrichment for these annotations, plus significant evidence of cooperation among these five TFs. Based on known gene/drug interactions, we prioritized targets in the OI-MET-TF network for follow-on analysis, emphasizing the clinical relevance of this work. Reflecting these results back to the OI-MET model, we found that binding motifs for the TF pair AP1/MYC are more frequent than expected and that the AP1/MYC pair is significantly enriched in binding in cancer models, relative to non-cancer models, in these promoters. This effect is seen in both MET models (solid tumors) and in non-MET models (leukemia). These results are consistent with our hypothesis that the OVOLs impact cancer susceptibility by regulating MET, and extend the hypothesis to include mechanisms not specific to MET. CONCLUSIONS: We find significant evidence of the OVOL, AP1, STAT1, STAT3, and NFKB1 TFs having important roles in MET, and more broadly in cancer. We prioritize known gene/drug targets for follow-up in the clinic, and we show that the AP1/MYC TF pair is a strong candidate for intervention.

Novel bioinformatics method for identification of genome-wide non-canonical spliced regions using RNA-Seq data.

SETTING: During endoplasmic reticulum (ER) stress, the endoribonuclease (RNase) Ire1α initiates removal of a 26 nt region from the mRNA encoding the transcription factor Xbp1 via an unconventional mechanism (atypically within the cytosol). This causes an open reading frame-shift that leads to altered transcriptional regulation of numerous downstream genes in response to ER stress as part of the unfolded protein response (UPR). Strikingly, other examples of targeted, unconventional splicing of short mRNA regions have yet to be reported. OBJECTIVE: Our goal was to develop an approach to identify non-canonical, possibly very short, splicing regions using RNA-Seq data and apply it to ER stress-induced Ire1α heterozygous and knockout mouse embryonic fibroblast (MEF) cell lines to identify additional Ire1α targets. RESULTS: We developed a bioinformatics approach called the Read-Split-Walk (RSW) pipeline, and evaluated it using two Ire1α heterozygous and two Ire1α-null samples. The 26 nt non-canonical splice site in Xbp1 was detected as the top hit by our RSW pipeline in heterozygous samples but not in the negative control Ire1α knockout samples. We compared the Xbp1 results from our approach with results using the alignment program BWA, Bowtie2, STAR, Exonerate and the Unix "grep" command. We then applied our RSW pipeline to RNA-Seq data from the SKBR3 human breast cancer cell line. RSW reported a large number of non-canonical spliced regions for 108 genes in chromosome 17, which were identified by an independent study. CONCLUSIONS: We conclude that our RSW pipeline is a practical approach for identifying non-canonical splice junction sites on a genome-wide level. We demonstrate that our pipeline can detect novel splice sites in RNA-Seq data generated under similar conditions for multiple species, in our case mouse and human.