Nicole Cloonan

Pancreatic cancer remains one of the most lethal of malignancies and a major health burden. We performed whole-genome sequencing and copy number variation (CNV) analysis of 100 pancreatic ductal adenocarcinomas (PDACs). Chromosomal rearrangements leading to gene disruption were prevalent, affecting genes known to be important in pancreatic cancer (TP53, SMAD4, CDKN2A, ARID1A and ROBO2) and new candidate drivers of pancreatic carcinogenesis (KDM6A and PREX2). Patterns of structural variation (variation in chromosomal structure) classified PDACs into 4 subtypes with potential clinical utility: the subtypes were termed stable, locally rearranged, scattered and unstable. A significant proportion harboured focal amplifications, many of which contained druggable oncogenes (ERBB2, MET, FGFR1, CDK6, PIK3R3 and PIK3CA), but at low individual patient prevalence. Genomic instability co-segregated with inactivation of DNA maintenance genes (BRCA1, BRCA2 or PALB2) and a mutational signature of DNA damage repair deficiency. Of 8 patients who received platinum therapy, 4 of 5 individuals with these measures of defective DNA maintenance responded.

Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

Upregulation of the immune response may be involved in the pathogenesis of schizophrenia with changes occurring in both peripheral blood and brain tissue. To date, microarray technology has provided a limited view of specific inflammatory transcripts in brain perhaps due to sensitivity issues. Here we used SOLiD Next Generation Sequencing to quantify neuroimmune mRNA expression levels in the dorsolateral prefrontal cortex of 20 individuals with schizophrenia and their matched controls. We detected 798 differentially regulated transcripts present in people with schizophrenia compared with controls. Ingenuity pathway analysis identified the inflammatory response as a key change. Using quantitative real-time PCR we confirmed the changes in candidate cytokines and immune modulators, including interleukin (IL)-6, IL-8, IL-1β and SERPINA3. The density of major histocompatibility complex-II-positive cells morphologically resembling microglia was significantly increased in schizophrenia and correlated with IL-1β expression. A group of individuals, most of whom had schizophrenia, were found to have increased inflammatory mRNA expression. In summary, we have demonstrated changes in an inflammatory response pathway that are present in ∼40% of people diagnosed with schizophrenia. This suggests that therapies aimed at immune system attenuation in schizophrenia may be of direct benefit in the brain.

The FANTOM4 study identified transcriptional start sites active during proliferation arrest and differentiation of the human monocytic cell line THP-1. Systematic knockdown of 52 transcription factors provide support for their model in which a complex transcriptional network regulates the differentiation process. Using deep sequencing (deepCAGE), the FANTOM4 study measured the genome-wide dynamics of transcription-start-site usage in the human monocytic cell line THP-1 throughout a time course of growth arrest and differentiation. Modeling the expression dynamics in terms of predicted cis-regulatory sites, we identified the key transcription regulators, their time-dependent activities and target genes. Systematic siRNA knockdown of 52 transcription factors confirmed the roles of individual factors in the regulatory network. Our results indicate that cellular states are constrained by complex networks involving both positive and negative regulatory interactions among substantial numbers of transcription factors and that no single transcription factor is both necessary and sufficient to drive the differentiation process.

The ability of T cells, acting independently of antibodies, to control malaria parasite growth in people has not been defined. If such was shown to be effective, an additional vaccine strategy could be pursued. Our aim was to ascertain whether or not development of cell-mediated immunity to Plasmodium falciparum blood-stage infection could be induced in human beings by exposure to malaria parasites in very low density.We enrolled five volunteers from the staff at our research institute who had never had malaria. We used a cryopreserved inoculum of red cells infected with P falciparum strain 3D7 to give them repeated subclinical infections of malaria that we then cured early with drugs, to induce cell-mediated immune responses. We tested for development of immunity by measurement of parasite concentrations in the blood of volunteers by PCR of the multicopy gene STEVOR and by following up the volunteers clinically, and by measuring antibody and cellular immune responses to the parasite.After challenge and a extended period without drug cure, volunteers were protected against malaria as indicated by absence of parasites or parasite DNA in the blood, and absence of clinical symptoms. Immunity was characterised by absence of detectable antibodies that bind the parasite or infected red cells, but by the presence of a proliferative T-cell response, involving CD4+ and CD8+ T cells, a cytokine response, consisting of interferon gamma but not interleukin 4 or interleukin 10, induction of high concentrations of nitric oxide synthase activity in peripheral blood mononuclear cells, and a drop in the number of peripheral natural killer T cells.People can be protected against the erythrocytic stage of malaria by a strong cell-mediated immune response, in the absence of detectable parasite-specific antibodies, suggesting an additional strategy for development of a malaria vaccine

Variants of microRNAs (miRNAs), called isomiRs, are commonly reported in deep-sequencing studies; however, the functional significance of these variants remains controversial. Observational studies show that isomiR patterns are non-random, hinting that these molecules could be regulated and therefore functional, although no conclusive biological role has been demonstrated for these molecules. To assess the biological relevance of isomiRs, we have performed ultra-deep miRNA-seq on ten adult human tissues, and created an analysis pipeline called miRNA-MATE to align, annotate, and analyze miRNAs and their isomiRs. We find that isomiRs share sequence and expression characteristics with canonical miRNAs, and are generally strongly correlated with canonical miRNA expression. A large proportion of isomiRs potentially derive from AGO2 cleavage independent of Dicer. We isolated polyribosome-associated mRNA, captured the mRNA-bound miRNAs, and found that isomiRs and canonical miRNAs are equally associated with translational machinery. Finally, we transfected cells with biotinylated RNA duplexes encoding isomiRs or their canonical counterparts and directly assayed their mRNA targets. These studies allow us to experimentally determine genome-wide mRNA targets, and these experiments showed substantial overlap in functional mRNA networks suppressed by both canonical miRNAs and their isomiRs. Together, these results find isomiRs to be biologically relevant and functionally cooperative partners of canonical miRNAs that act coordinately to target pathways of functionally related genes. This work exposes the complexity of the miRNA-transcriptome, and helps explain a major miRNA paradox: how specific regulation of biological processes can occur when the specificity of miRNA targeting is mediated by only 6 to 11 nucleotides.

MicroRNAs are modifiers of gene expression, acting to reduce translation through either translational repression or mRNA cleavage. Recently, it has been shown that some microRNAs can act to promote or suppress cell transformation, with miR-17-92 described as the first oncogenic microRNA. The association of miR-17-92 encoded microRNAs with a surprisingly broad range of cancers not only underlines the clinical significance of this locus, but also suggests that miR-17-92 may regulate fundamental biological processes, and for these reasons miR-17-92 has been considered as a therapeutic target.In this study, we show that miR-17-92 is a cell cycle regulated locus, and ectopic expression of a single microRNA (miR-17-5p) is sufficient to drive a proliferative signal in HEK293T cells. For the first time, we reveal the mechanism behind this response - miR-17-5p acts specifically at the G1/S-phase cell cycle boundary, by targeting more than 20 genes involved in the transition between these phases. While both pro- and anti-proliferative genes are targeted by miR-17-5p, pro-proliferative mRNAs are specifically up-regulated by secondary and/or tertiary effects in HEK293T cells.The miR-17-5p microRNA is able to act as both an oncogene and a tumor suppressor in different cellular contexts; our model of competing positive and negative signals can explain both of these activities. The coordinated suppression of proliferation-inhibitors allows miR-17-5p to efficiently de-couple negative regulators of the MAPK (mitogen activated protein kinase) signaling cascade, promoting growth in HEK293T cells. Additionally, we have demonstrated the utility of a systems biology approach as a unique and rapid approach to uncover microRNA function.

Background Interleukin 6 mediates graft-versus-host disease (GVHD) in experimental allogeneic stem-cell transplantation (allogeneic SCT) and represents an attractive therapeutic target. We aimed to assess whether the humanised anti-interleukin-6 receptor monoclonal antibody, tocilizumab, could attenuate the incidence of acute GVHD. Methods We undertook a single-group, single-institution phase 1/2 study at the Royal Brisbane and Women's Hospital Bone Marrow Transplantation unit, QLD, Australia. Eligible patients were 18–65 years old and underwent T-replete HLA-matched allogeneic SCT with either total body irradiation-based myeloablative or reduced-intensity conditioning from unrelated or sibling donors. One intravenous dose of tocilizumab (8 mg/kg, capped at 800 mg, over 60 mins' infusion) was given the day before allogeneic SCT along with standard GVHD prophylaxis (cyclosporin [5 mg/kg per day on days −1 to +1, then 3 mg/kg per day to maintain therapeutic levels (trough levels of 140–300 ng/mL) for 100 days plus methotrexate [15 mg/m2 on day 1, then 10 mg/m2 on days 3, 6, and 11]). The primary endpoint was incidence of grade 2–4 acute GVHD at day 100, assessed and graded as per the Seattle criteria. Immunological profiles were compared with a non-randomised group of patients receiving allogeneic SCT, but not treated with tocilizumab. This trial is registered with the Australian and New Zealand Clinical Trials Registry, number ACTRN12612000726853. Findings Between Jan 19, 2012, and Aug 27, 2013, 48 eligible patients receiving cyclosporin and methotrexate as GVHD prophylaxis were enrolled into the study. The incidence of grade 2–4 acute GVHD in patients treated with tocilizumab at day 100 was 12% (95% CI 5–24), and the incidence of grade 3–4 acute GVHD was 4% (1–13). Grade 2–4 acute GVHD involving the skin developed in five (10%) patients of 48 treated with tocilizumab, involving the gastrointestinal tract in four (8%) patients; there were no reported cases involving the liver. Low incidences of grade 2–4 acute GVHD were noted in patients receiving both myeloablative total body irradiation-based conditioning (12% [95% CI 2–34) and fludarabine and melphalan reduced-intensity conditioning (12% [4–27]). Immune reconstitution was preserved in recipients of interleukin-6 receptor inhibition, but qualitatively modified with suppression of known pathogenic STAT3-dependent pathways. Interpretation Interleukin 6 is the main detectable and dysregulated cytokine secreted after allogeneic SCT and its inhibition is a potential new and simple strategy to protect from acute GVHD despite robust immune reconstitution; a randomised, controlled trial assessing tocilizumab in addition to standard GVHD prophylaxis in these patients is warranted. Funding National Health and Medical Research Council and Queensland Health.

Immune-modulating therapies have revolutionized the treatment of chronic diseases, particularly cancer. However, their success is restricted and there is a need to identify new therapeutic targets. Here, we show that natural killer cell granule protein 7 (NKG7) is a regulator of lymphocyte granule exocytosis and downstream inflammation in a broad range of diseases. NKG7 expressed by CD4+ and CD8+ T cells played key roles in promoting inflammation during visceral leishmaniasis and malaria-two important parasitic diseases. Additionally, NKG7 expressed by natural killer cells was critical for controlling cancer initiation, growth and metastasis. NKG7 function in natural killer and CD8+ T cells was linked with their ability to regulate the translocation of CD107a to the cell surface and kill cellular targets, while NKG7 also had a major impact on CD4+ T cell activation following infection. Thus, we report a novel therapeutic target expressed on a range of immune cells with functions in different immune responses.

Several multicopy gene families have been described in Plasmodium falciparum, including the var genes that code for the variant surface antigen PfEMP1, the stevor family of subtelomeric open reading frames and the rif interspersed repetitive elements. This report documents the chromosomal location of stevor genes, their transcription and characteristics of the deduced protein. On 14 chromosomes, 34 stevor copies were identified from the Dd2 parasite line. Most are in subtelomeric regions within 50 kb of the telomere. stevor genes are located close to var genes and rij. All stevor genes sequenced had two exons: a short exon 1 encoding a start codon and a transmembrane domain; exon 2 encoding for the remainder of the approximately 30 kDa protein and including two more transmembrane segments. A similar structure was found for copies of rif and its predicted protein. In both STEVOR and RIF proteins, a highly polymorphic region is predicted to be a loop on the outer side of the membrane. We propose that stevor and rif are members of a larger superfamily. The number of copies of stevor and rif, their location close to the var genes, their extreme polymorphism and the predicted structure of the proteins suggest that stevor and rif code for variant surface antigens.

RAD51 is an important component of double-stranded DNA-repair mechanisms that interacts with both BRCA1 and BRCA2. A single-nucleotide polymorphism (SNP) in the 5' untranslated region (UTR) of RAD51, 135G-->C, has been suggested as a possible modifier of breast cancer risk in BRCA1 and BRCA2 mutation carriers. We pooled genotype data for 8,512 female mutation carriers from 19 studies for the RAD51 135G-->C SNP. We found evidence of an increased breast cancer risk in CC homozygotes (hazard ratio [HR] 1.92 [95% confidence interval {CI} 1.25-2.94) but not in heterozygotes (HR 0.95 [95% CI 0.83-1.07]; P=.002, by heterogeneity test with 2 degrees of freedom [df]). When BRCA1 and BRCA2 mutation carriers were analyzed separately, the increased risk was statistically significant only among BRCA2 mutation carriers, in whom we observed HRs of 1.17 (95% CI 0.91-1.51) among heterozygotes and 3.18 (95% CI 1.39-7.27) among rare homozygotes (P=.0007, by heterogeneity test with 2 df). In addition, we determined that the 135G-->C variant affects RAD51 splicing within the 5' UTR. Thus, 135G-->C may modify the risk of breast cancer in BRCA2 mutation carriers by altering the expression of RAD51. RAD51 is the first gene to be reliably identified as a modifier of risk among BRCA1/2 mutation carriers.

KLF1 regulates a diverse suite of genes to direct erythroid cell differentiation from bipotent progenitors. To determine the local cis-regulatory contexts and transcription factor networks in which KLF1 operates, we performed KLF1 ChIP-seq in the mouse. We found at least 945 sites in the genome of E14.5 fetal liver erythroid cells which are occupied by endogenous KLF1. Many of these recovered sites reside in erythroid gene promoters such as Hbb-b1, but the majority are distant to any known gene. Our data suggests KLF1 directly regulates most aspects of terminal erythroid differentiation including production of alpha- and beta-globin protein chains, heme biosynthesis, coordination of proliferation and anti-apoptotic pathways, and construction of the red cell membrane and cytoskeleton by functioning primarily as a transcriptional activator. Additionally, we suggest new mechanisms for KLF1 cooperation with other transcription factors, in particular the erythroid transcription factor GATA1, to maintain homeostasis in the erythroid compartment.

Metastasis is a complex, multistep process involved in the progression of cancer from a localized primary tissue to distant sites, often characteristic of the more aggressive forms of this disease. Despite being studied in great detail in recent years, the mechanisms that govern this process remain poorly understood. In this study, we identify a novel role for miR-139-5p in the inhibition of breast cancer progression. We highlight its clinical relevance by reviewing miR-139-5p expression across a wide variety of breast cancer subtypes using in-house generated and online data sets to show that it is most frequently lost in invasive tumors. A biotin pull-down approach was then used to identify the mRNA targets of miR-139-5p in the breast cancer cell line MCF7. Functional enrichment analysis of the pulled-down targets showed significant enrichment of genes in pathways previously implicated in breast cancer metastasis (P < 0.05). Further bioinformatic analysis revealed a predicted disruption to the TGFβ, Wnt, Rho, and MAPK/PI3K signaling cascades, implying a potential role for miR-139-5p in regulating the ability of cells to invade and migrate. To corroborate this finding, using the MDA-MB-231 breast cancer cell line, we show that overexpression of miR-139-5p results in suppression of these cellular phenotypes. Furthermore, we validate the interaction between miR-139-5p and predicted targets involved in these pathways. Collectively, these results suggest a significant functional role for miR-139-5p in breast cancer cell motility and invasion and its potential to be used as a prognostic marker for the aggressive forms of breast cancer.

MicroRNAs are noncoding regulators of gene expression, which act by repressing protein translation and/or degrading mRNA. Many have been shown to drive tumorigenesis in cancer, but functional studies to understand their mode of action are typically limited to single-target genes. In this study, we use synthetic biotinylated miRNA to pull down endogenous targets of miR-182-5p. We identified more than 1000 genes as potential targets of miR-182-5p, most of which have a known function in pathways underlying tumor biology. Specifically, functional enrichment analysis identified components of both the DNA damage response pathway and cell cycle to be highly represented in this target cohort. Experimental validation confirmed that miR-182-5p-mediated disruption of the homologous recombination (HR) pathway is a consequence of its ability to target multiple components in that pathway. Although there is a strong enrichment for the cell cycle ontology, we do not see primary proliferative defects as a consequence of miR-182-5p overexpression. We highlight targets that could be responsible for miR-182-5p-mediated disruption of other biological processes attributed in the literature so far. Finally, we show that miR-182-5p is highly expressed in a panel of human breast cancer samples, highlighting its role as a potential oncomir in breast cancer.

MicroRNAs (miRNAs) bind to mRNAs and target them for translational inhibition or transcriptional degradation. It is thought that most miRNA-mRNA interactions involve the seed region at the 5′ end of the miRNA. The importance of seed sites is supported by experimental evidence, although there is growing interest in interactions mediated by the central region of the miRNA, termed centered sites. To investigate the prevalence of these interactions, we apply a biotin pull-down method to determine the direct targets of ten human miRNAs, including four isomiRs that share centered sites, but not seeds, with their canonical partner miRNAs. We confirm that miRNAs and their isomiRs can interact with hundreds of mRNAs, and that imperfect centered sites are common mediators of miRNA-mRNA interactions. We experimentally demonstrate that these sites can repress mRNA activity, typically through translational repression, and are enriched in regions of the transcriptome bound by AGO. Finally, we show that the identification of imperfect centered sites is unlikely to be an artifact of our protocol caused by the biotinylation of the miRNA. However, the fact that there was a slight bias against seed sites in our protocol may have inflated the apparent prevalence of centered site-mediated interactions. Our results suggest that centered site-mediated interactions are much more frequent than previously thought. This may explain the evolutionary conservation of the central region of miRNAs, and has significant implications for decoding miRNA-regulated genetic networks, and for predicting the functional effect of variants that do not alter protein sequence.

Reprogramming of somatic cells to induced pluripotent stem cells involves a dynamic rearrangement of the epigenetic landscape. To characterize this epigenomic roadmap, we have performed MethylC-seq, ChIP-seq (H3K4/K27/K36me3) and RNA-Seq on samples taken at several time points during murine secondary reprogramming as part of Project Grandiose. We find that DNA methylation gain during reprogramming occurs gradually, while loss is achieved only at the ESC-like state. Binding sites of activated factors exhibit focal demethylation during reprogramming, while ESC-like pluripotent cells are distinguished by extension of demethylation to the wider neighbourhood. We observed that genes with CpG-rich promoters demonstrate stable low methylation and strong engagement of histone marks, whereas genes with CpG-poor promoters are safeguarded by methylation. Such DNA methylation-driven control is the key to the regulation of ESC-pluripotency genes, including Dppa4, Dppa5a and Esrrb. These results reveal the crucial role that DNA methylation plays as an epigenetic switch driving somatic cells to pluripotency.

Despite a library full of literature on miRNA biology, core issues relating to miRNA target detection, biological effect, and mode of action remain controversial. This essay proposes that the predominant mechanism of direct miRNA action is translational inhibition, whereas the bulk of miRNA effects are mRNA based. It explores several issues confounding miRNA target detection, and discusses their impact on the dominance of "miRNA seed" dogma and the exploration of non-canonical binding sites. Finally, it makes comparisons between miRNA target prediction and transcription factor binding prediction, and questions the value of characterizing miRNA binding sites based on which miRNA nucleotides are paired with an mRNA.

Abstract Mechanotransduction is a strong driver of mesenchymal stem cell (MSC) fate. In vitro, variations in matrix mechanics invoke changes in MSC proliferation, migration and differentiation. However, when incorporating MSCs within injectable, inherently soft hydrogels, this dominance over MSC response substantially limits our ability to couple the ease of application of hydrogels with efficiently directed MSC differentiation, especially in the case of bone generation. Here, we identify differential miRNA expression in response to varying hydrogel stiffness and RhoA activity. We show that modulation of miR-100-5p and miR-143-3p can be used to bias MSC fate and provide mechanistic insight by demonstrating convergence on mTOR signalling. By modulating these mechanosensitive miRNAs, we can enhance osteogenesis in a soft 3D hydrogel. The outcomes of this study provide new understanding of the mechanisms regulating MSC mechanotransduction and differentiation, but also a novel strategy with which to drive MSC fate and significantly impact MSC-based tissue-engineering applications.

Breast cancer risk is strongly associated with an intergenic region on 11q13. We have previously shown that the strongest risk-associated SNPs fall within a distal enhancer that regulates <i>CCND1</i>. Here, we report that, in addition to regulating <i>CCND1</i>, this enhancer regulates two estrogen-regulated long noncoding RNAs, <i>CUPID1</i> and <i>CUPID2.</i> We provide evidence that the risk-associated SNPs are associated with reduced chromatin looping between the enhancer and the <i>CUPID1</i> and <i>CUPID2</i> bidirectional promoter. We further show that <i>CUPID1</i> and <i>CUPID2</i> are predominantly expressed in hormone-receptor-positive breast tumors and play a role in modulating pathway choice for the repair of double-strand breaks. These data reveal a mechanism for the involvement of this region in breast cancer.

A placebo controlled, randomised, double blind trial was conducted in human volunteers to test a mixture of three recombinant Plasmodium falciparum blood stage antigens for its ability to reduce the initial growth rates of parasites. The vaccine contained recombinant MSP2 (3D7 allele), a portion of MSP1 (190LCS.T3) and part of the RESA antigen (C terminal 771 amino acids) in the Montanide ISA 720 adjuvant (SEPPIC). Twelve volunteers received two doses of the vaccine, 6 weeks apart. The five participants in the placebo group received an equivalent volume of the adjuvant emulsion using the same schedule. Antibody responses were low, as has been reported in earlier studies with this combination, while T cell responses were stronger. All the volunteers were challenged with approximately 140 ring infected red cells of the 3D7 cloned line, 4 weeks after the second dose. Parasitaemia was determined once daily from day 4 using a sensitive and quantitative PCR assay. All the volunteers were infected and were treated on day 8, before any developed symptoms. There was no significant difference in initial parasite growth rates between the verum and placebo groups, nor was there any significant correlation between parasite growth rates and any of the measured immunological responses. These results suggest that the formulation tested in this trial did not generate immune responses that were strong enough to reduce parasite growth in naive volunteers.

Massively parallel short-tag sequencing of cDNA libraries--RNAseq--is being used to study the dynamics and complexity of eukaryotic transcriptomes, giving new biological insights into the 'active genome'.

Human Sin1 (SAPK-interacting protein 1) is a member of a conserved family of orthologous proteins that have an essential role in signal transduction in yeast and Dictyostelium. This study demonstrates that most Sin1 orthologues contain both a Raf-like Ras-binding domain (RBD) and a pleckstrin homology (PH) domain. These domains are functional in the human Sin1 protein, with the PH domain involved in lipid and membrane binding by Sin1, and the RBD able to bind activated H-and K-Ras. Sin1 and Ras co-immunoprecipitated and co-localised, showing that these proteins associate with each other in vivo. Overexpression of Sin1 inhibited the activation of ERK, Akt and JNK signalling pathways by Ras. In contrast, siRNA knockdown of endogenous Sin1 protein expression in HEK293 cells enhanced the activation of ERK1/2 by Ras. These data suggest that Sin1 is a mammalian Ras-inhibitor.

Minor class or U12-type splicing is a highly conserved process required to remove a minute fraction of introns from human pre-mRNAs. Defects in this splicing pathway have recently been linked to human disease, including a severe developmental disorder encompassing brain and skeletal abnormalities known as Taybi-Linder syndrome or microcephalic osteodysplastic primordial dwarfism 1, and a hereditary intestinal polyposis condition, Peutz-Jeghers syndrome. Although a key mechanism for regulating gene expression, the impact of impaired U12-type splicing on the transcriptome is unknown. Here, we describe a unique zebrafish mutant, caliban (clbn), with arrested development of the digestive organs caused by an ethylnitrosourea-induced recessive lethal point mutation in the rnpc3 [RNA-binding region (RNP1, RRM) containing 3] gene. rnpc3 encodes the zebrafish ortholog of human RNPC3, also known as the U11/U12 di-snRNP 65-kDa protein, a unique component of the U12-type spliceosome. The biochemical impact of the mutation in clbn is the formation of aberrant U11- and U12-containing small nuclear ribonucleoproteins that impair the efficiency of U12-type splicing. Using RNA sequencing and microarrays, we show that multiple genes involved in various steps of mRNA processing, including transcription, splicing, and nuclear export are disrupted in clbn, either through intron retention or differential gene expression. Thus, clbn provides a useful and specific model of aberrant U12-type splicing in vivo. Analysis of its transcriptome reveals efficient mRNA processing as a critical process for the growth and proliferation of cells during vertebrate development.

Polymorphic regions of the genes encoding Plasmodium vivax apical membrane antigen 1 (PvAMA1) and P. vivax merozoite surface protein 1 (PvMSP1) were sequenced to examine population diversity both within and between geographical areas. Sequences were obtained for 219 isolates for PvAMA1 and for 175 isolates for PvMSP1 from Africa, China, India, Indonesia, Philippines, Papua New Guinea, Solomon Islands and Thailand. Over half of the isolates were obtained from different regions within the Philippines, and this was used to look at the diversity within a country. Sixty nine haplotypes and 22 polymorphic sites in a 414-bp region of PvAMA1 and 41 haplotypes and 34 polymorphic sites in a 249-bp fragment of PvMSP1 were detected. For both PvAMA1 and PvMSP1, four previously unreported polymorphic nucleotide positions were identified. Population analysis indicated that there were significant differences in allele frequencies between different regions but these differences were small compared to the diversity within populations (Fixation index, F(ST), of 0.126 and 0.078 for PvAMA1 and PvMSP1, respectively). PvAMA1 and PvMSP1 had similar nonsynonymous substitution frequencies but surprisingly, the synonymous substitution frequency for PvMSP1 was eight times the frequency for PvAMA1 suggesting that synonymous substitutions in at least PvAMA1 are not neutral.

Surface marker expression forms the basis for characterization and isolation of human embryonic stem cells (hESCs). Currently, there are few well-defined protein epitopes that definitively mark hESCs. Here we combine immunotranscriptional profiling of hESC lines with membrane-polysome translation state array analysis (TSAA) to determine the full set of genes encoding potential hESC surface marker proteins. Three independently isolated hESC lines (HES2, H9, and MEL1) grown under feeder and feeder-free conditions were sorted into subpopulations by fluorescence-activated cell sorting based on coimmunoreactivity to the hESC surface markers GCTM-2 and CD9. Colony-forming assays confirmed that cells displaying high coimmunoreactivity to GCTM-2 and CD9 constitute an enriched subpopulation displaying multiple stem cell properties. Following microarray profiling, 820 genes were identified that were common to the GCTM-2(high)/CD9(high) stem cell-like subpopulation. Membrane-polysome TSAA analysis of hESCs identified 1,492 mRNAs encoding actively translated plasma membrane and secreted proteins. Combining these data sets, 88 genes encode proteins that mark the pluripotent subpopulation, of which only four had been previously reported. Cell surface immunoreactivity was confirmed for two of these markers: TACSTD1/EPCAM and CDH3/P-Cadherin, with antibodies for EPCAM able to enrich for pluripotent hESCs. This comprehensive listing of both hESCs and spontaneous differentiation-associated transcripts and survey of translated membrane-bound and secreted proteins provides a valuable resource for future study into the role of the extracellular environment in both the maintenance of pluripotency and directed differentiation.

Quantification applications of short-tag sequencing data (such as CNVseq and RNAseq) depend on knowing the uniqueness of specific genomic regions at a given threshold of error. Here, we present the 'uniqueome', a genomic resource for understanding the uniquely mappable proportion of genomic sequences. Pre-computed data are available for human, mouse, fly and worm genomes in both color-space and nucletotide-space, and we demonstrate the utility of this resource as applied to the quantification of RNAseq data.Files, scripts and supplementary data are available from http://grimmond.imb.uq.edu.au/uniqueome/; the ISAS uniqueome aligner is freely available from http://www.imagenix.com/.

MicroRNAs (miRNAs) are critical to somatic cell reprogramming into induced pluripotent stem cells (iPSCs), however, exactly how miRNA expression changes support the transition to pluripotency requires further investigation. Here we use a murine secondary reprogramming system to sample cellular trajectories towards iPSCs or a novel pluripotent 'F-class' state and perform small RNA sequencing. We detect sweeping changes in an early and a late wave, revealing that distinct miRNA milieus characterize alternate states of pluripotency. miRNA isoform expression is common but surprisingly varies little between cell states. Referencing other omic data sets generated in parallel, we find that miRNA expression is changed through transcriptional and post-transcriptional mechanisms. miRNA transcription is commonly regulated by dynamic histone modification, while DNA methylation/demethylation consolidates these changes at multiple loci. Importantly, our results suggest that a novel subset of distinctly expressed miRNAs supports pluripotency in the F-class state, substituting for miRNAs that serve such roles in iPSCs.

mRNA synthesis, processing, and destruction involve a complex series of molecular steps that are incompletely understood. Because the RNA intermediates in each of these steps have finite lifetimes, extensive mechanistic and dynamical information is encoded in total cellular RNA. Here we report the development of SnapShot-Seq, a set of computational methods that allow the determination of in vivo rates of pre-mRNA synthesis, splicing, intron degradation, and mRNA decay from a single RNA-Seq snapshot of total cellular RNA. SnapShot-Seq can detect in vivo changes in the rates of specific steps of splicing, and it provides genome-wide estimates of pre-mRNA synthesis rates comparable to those obtained via labeling of newly synthesized RNA. We used SnapShot-Seq to investigate the origins of the intrinsic bimodality of metazoan gene expression levels, and our results suggest that this bimodality is partly due to spillover of transcriptional activation from highly expressed genes to their poorly expressed neighbors. SnapShot-Seq dramatically expands the information obtainable from a standard RNA-Seq experiment.

The ectopic expression of Oct4, Klf4, c-Myc and Sox2 (OKMS) transcription factors allows reprogramming of somatic cells into induced pluripotent stem cells (iPSCs). The reprogramming process, which involves a complex network of molecular events, is not yet fully characterized. Here we perform a quantitative mass spectrometry-based analysis to probe in-depth dynamic proteome changes during somatic cell reprogramming. Our data reveal defined waves of proteome resetting, with the first wave occurring 48 h after the activation of the reprogramming transgenes and involving specific biological processes linked to the c-Myc transcriptional network. A second wave of proteome reorganization occurs in a later stage of reprogramming, where we characterize the proteome of two distinct pluripotent cellular populations. In addition, the overlay of our proteome resource with parallel generated -omics data is explored to identify post-transcriptionally regulated proteins involved in key steps during reprogramming. During somatic cell reprogramming, the cell transits through intermediate states. Here, the authors perform an in-depth quantitative proteomic analysis of the reprogramming of mouse embryonic fibroblasts to induced pluripotent stem cells and observe two waves of proteome reorganisation.

Genome-wide association studies (GWAS) have identified markers within the WNT4 region on chromosome 1p36.12 showing consistent and strong association with increasing endometriosis risk. Fine mapping using sequence and imputed genotype data has revealed strong candidates for the causal SNPs within these critical regions; however, the molecular pathogenesis of these SNPs is currently unknown. We used gene expression data collected from whole blood from 862 individuals and endometrial tissue from 136 individuals from independent populations of European descent to examine the mechanism underlying endometriosis susceptibility. Association mapping results from 7,090 individuals (2,594 cases and 4,496 controls) supported rs3820282 as the SNP with the strongest association for endometriosis risk (P = 1.84 × 10−5, OR = 1.244 (1.126-1.375)). SNP rs3820282 is a significant eQTL in whole blood decreasing expression of LINC00339 (also known as HSPC157) and increasing expression of CDC42 (P = 2.0 ×10−54 and 4.5x10−4 respectively). The largest effects were for two LINC00339 probes (P = 2.0 ×10−54; 1.0 × 10−34). The eQTL for LINC00339 was also observed in endometrial tissue (P = 2.4 ×10−8) with the same direction of effect for both whole blood and endometrial tissue. There was no evidence for eQTL effects for WNT4. Chromatin conformation capture provides evidence for risk SNPs interacting with the promoters of both LINC00339 and CDC4 and luciferase reporter assays suggest the risk SNP rs12038474 is located in a transcriptional silencer for CDC42 and the risk allele increases expression of CDC42. However, no effect of rs3820282 was observed in the LINC00339 expression in Ishikawa cells. Taken together, our results suggest that SNPs increasing endometriosis risk in this region act through CDC42, but further functional studies are required to rule out inverse regulation of both LINC00339 and CDC42.

Abstract Background Caveolae proteins play diverse roles in cancer development and progression. In prostate cancer, non‐caveolar caveolin‐1 (CAV1) promotes metastasis, while CAVIN1 attenuates CAV1‐induced metastasis. Here, we unveil a novel mechanism linking CAV1 to selective loading of exosomes with metastasis‐promoting microRNAs. Results We identify hnRNPK as a CAV1‐regulated microRNA binding protein. In the absence of CAVIN1, non‐caveolar CAV1 drives localisation of hnRPNK to multi‐vesicular bodies (MVBs), recruiting AsUGnA motif‐containing miRNAs and causing their release within exosomes. This process is dependent on the lipid environment of membranes as shown by cholesterol depletion using methyl‐β‐cyclodextrin or by treatment with n‐3 polyunsaturated fatty acids. Consistent with a role in bone metastasis, knockdown of hnRNPK in prostate cancer PC3 cells abolished the ability of PC3 extracellular vesicles (EV) to induce osteoclastogenesis, and biofluid EV hnRNPK is elevated in metastatic prostate and colorectal cancer. Conclusions Taken together, these results support a novel pan‐cancer mechanism for CAV1‐driven exosomal release of hnRNPK and associated miRNA in metastasis, which is modulated by the membrane lipid environment.

Human proximal tubular epithelial cells (PTEC) of the kidney are known to respond to and mediate the disease process in a wide range of kidney diseases, yet their exosomal production and exosome molecular cargo remain a mystery. Here we investigate, for the first time, the production and molecular content of exosomes derived from primary human PTEC cultured under normal and diseased conditions representing a spectrum of in vivo disease severity from early inflammation, experienced in multiple initial kidney disease states, through to hypoxia, frequently seen in late stage chronic kidney disease (CKD) due to fibrosis and vascular compromise. We demonstrate a rapid reproducible methodology for the purification of PTEC-derived exosomes, identify increased numbers of exosomes from disease-state cultures and identify differential expression levels of both known and unique miRNA and protein species from exosomes derived from different disease-culture conditions. The validity of our approach is supported by the identification of miRNA, proteins and pathways with known CKD associations, providing a rationale to further evaluate these novel and known pathways as targets for therapeutic intervention.

With five-year survival rates as low as 3%, lung cancer is the most common cause of cancer-related mortality worldwide. The severity of the disease at presentation is accredited to the lack of early detection capacities, resulting in the reliance on low-throughput diagnostic measures, such as tissue biopsy and imaging. Interest in the development and use of liquid biopsies has risen, due to non-invasive sample collection, and the depth of information it can provide on a disease. Small extracellular vesicles (sEVs) as viable liquid biopsies are of particular interest due to their potential as cancer biomarkers. To validate the use of sEVs as cancer biomarkers, we characterised cancer sEVs using miRNA sequencing analysis. We found that miRNA-3182 was highly enriched in sEVs derived from the blood of patients with invasive breast carcinoma and NSCLC. The enrichment of sEV miR-3182 was confirmed in oncogenic, transformed lung cells in comparison to isogenic, untransformed lung cells. Most importantly, miR-3182 can successfully distinguish early-stage NSCLC patients from those with benign lung conditions. Therefore, miR-3182 provides potential to be used for the detection of NSCLC in blood samples, which could result in earlier therapy and thus improved outcomes and survival for patients.

Mapping of next-generation sequencing data derived from RNA samples (RNAseq) presents different genome mapping challenges than data derived from DNA. For example, tags that cross exon-junction boundaries will often not map to a reference genome, and the strand specificity of the data needs to be retained. Here we present RNA-MATE, a computational pipeline based on a recursive mapping strategy for placing strand specific RNAseq data onto a reference genome. Maximizing the mappable tags can provide significant savings in the cost of sequencing experiments. This pipeline provides an automatic and integrated way to align color-space sequencing data, collate this information and generate files for examining gene-expression data in a genomic context.Executables, source code, and exon-junction libraries are available from http://grimmond.imb.uq.edu.au/RNA-MATE/

More than 98% of a typical vertebrate genome does not code for proteins. Although non-coding regions are sprinkled with short (<200 bp) islands of evolutionarily conserved sequences, the function of most of these unannotated conserved islands remains unknown. One possibility is that unannotated conserved islands could encode non-coding RNAs (ncRNAs); alternatively, unannotated conserved islands could serve as promoter-distal regulatory factor binding sites (RFBSs) like enhancers. Here we assess these possibilities by comparing unannotated conserved islands in the human and mouse genomes to transcribed regions and to RFBSs, relying on a detailed case study of one human and one mouse cell type. We define transcribed regions by applying a novel transcript-calling algorithm to RNA-Seq data obtained from total cellular RNA, and we define RFBSs using ChIP-Seq and DNAse-hypersensitivity assays. We find that unannotated conserved islands are four times more likely to coincide with RFBSs than with unannotated ncRNAs. Thousands of conserved RFBSs can be categorized as insulators based on the presence of CTCF or as enhancers based on the presence of p300/CBP and H3K4me1. While many unannotated conserved RFBSs are transcriptionally active to some extent, the transcripts produced tend to be unspliced, non-polyadenylated and expressed at levels 10 to 100-fold lower than annotated coding or ncRNAs. Extending these findings across multiple cell types and tissues, we propose that most conserved non-coding genomic DNA in vertebrate genomes corresponds to promoter-distal regulatory elements.

The developing mouse kidney is currently the best-characterized model of organogenesis at a transcriptional level. Detailed spatial maps have been generated for gene expression profiling combined with systematic in situ screening. These studies, however, fall short of capturing the transcriptional complexity arising from each locus due to the limited scope of microarray-based technology, which is largely based on "gene-centric" models.To address this, the polyadenylated RNA and microRNA transcriptomes of the 15.5 dpc mouse kidney were profiled using strand-specific RNA-sequencing (RNA-Seq) to a depth sufficient to complement spatial maps from pre-existing microarray datasets. The transcriptional complexity of RNAs arising from mouse RefSeq loci was catalogued; including 3568 alternatively spliced transcripts and 532 uncharacterized alternate 3' UTRs. Antisense expressions for 60% of RefSeq genes was also detected including uncharacterized non-coding transcripts overlapping kidney progenitor markers, Six2 and Sall1, and were validated by section in situ hybridization. Analysis of genes known to be involved in kidney development, particularly during mesenchymal-to-epithelial transition, showed an enrichment of non-coding antisense transcripts extended along protein-coding RNAs.The resulting resource further refines the transcriptomic cartography of kidney organogenesis by integrating deep RNA sequencing data with locus-based information from previously published expression atlases. The added resolution of RNA-Seq has provided the basis for a transition from classical gene-centric models of kidney development towards more accurate and detailed "transcript-centric" representations, which highlights the extent of transcriptional complexity of genes that direct complex development events.

Abstract During chronic liver injury hepatic stellate cells (HSCs), the principal source of extracellular matrix in the fibrotic liver, transdifferentiate into pro-fibrotic myofibroblast-like cells - a process potentially regulated by microRNAs (miRNAs). Recently, we found serum miRNA-25-3p (miR-25) levels were upregulated in children with Cystic Fibrosis (CF) without liver disease, compared to children with CF-associated liver disease and healthy individuals. Here we examine the role of miR-25 in HSC biology. MiR-25 was detected in the human HSC cell line LX-2 and in primary murine HSCs, and increased with culture-induced activation. Transient overexpression of miR-25 inhibited TGF-β and its type 1 receptor (TGFBR1) mRNA expression, TGF-β-induced Smad2 phosphorylation and subsequent collagen1α1 induction in LX-2 cells. Pull-down experiments with biotinylated miR-25 revealed Notch signaling (co-)activators ADAM-17 and FKBP14 as miR-25 targets in HSCs. NanoString analysis confirmed miR-25 regulation of Notch- and Wnt-signaling pathways. Expression of Notch signaling pathway components and endogenous Notch1 signaling was downregulated in miR-25 overexpressing LX-2 cells, as were components of Wnt signaling such as Wnt5a. We propose that miR-25 acts as a negative feedback anti-fibrotic control during HSC activation by reducing the reactivity of HSCs to TGF-β-induced collagen expression and modulating the cross-talk between Notch, Wnt and TGF-β signaling.

Recent RNA-sequencing studies have shown remarkable complexity in the mammalian transcriptome. The ultimate impact of this complexity on the predicted proteomic output is less well defined. We have undertaken strand-specific RNA sequencing of multiple cellular RNA fractions (>20 Gb) to uncover the transcriptional complexity of human embryonic stem cells (hESCs). We have shown that human embryonic stem (ES) cells display a high degree of transcriptional diversity, with more than half of active genes generating RNAs that differ from conventional gene models. We found evidence that more than 1000 genes express long 5' and/or extended 3'UTRs, which was confirmed by "virtual Northern" analysis. Exhaustive sequencing of the membrane-polysome and cytosolic/untranslated fractions of hESCs was used to identify RNAs encoding peptides destined for secretion and the extracellular space and to demonstrate preferential selection of transcription complexity for translation in vitro. The impact of this newly defined complexity on known gene-centric network models such as the Plurinet and the cell surface signaling machinery in human ES cells revealed a significant expansion of known transcript isoforms at play, many predicting possible alternative functions based on sequence alterations within key functional domains.

Abstract The caudal-related homeobox transcription factor CDX2 is expressed in leukemic cells but not during normal blood formation. Retroviral overexpression of Cdx2 induces AML in mice, however the developmental stage at which CDX2 exerts its effect is unknown. We developed a conditionally inducible Cdx2 mouse model to determine the effects of in vivo, inducible Cdx2 expression in hematopoietic stem and progenitor cells (HSPCs). Cdx2-transgenic mice develop myelodysplastic syndrome with progression to acute leukemia associated with acquisition of additional driver mutations. Cdx2-expressing HSPCs demonstrate enrichment of hematopoietic-specific enhancers associated with pro-differentiation transcription factors. Furthermore, treatment of Cdx2 AML with azacitidine decreases leukemic burden. Extended scheduling of low-dose azacitidine shows greater efficacy in comparison to intermittent higher-dose azacitidine, linked to more specific epigenetic modulation. Conditional Cdx2 expression in HSPCs is an inducible model of de novo leukemic transformation and can be used to optimize treatment in high-risk AML.

Limited but significant sequence similarity has been observed between an uncharacterized human protein, SIN1, and the S. pombe SIN1, Dictyostelium RIP3 and S. cerevisiae AVO1 proteins. The human Sin1 gene has been automatically predicted (MAPKAP1; GenBank accession number ); however, this sequence appears to be incomplete. In this study, we have cloned and characterized the full-length human Sin1 mRNA and identified a highly conserved domain that defines the family of SIN1 orthologues, members of which are widely distributed in the fungal and metazoan kingdoms. We demonstrate that Sin1 transcripts can use alternative polyadenylation signals and describe a number of Sin1 splice variants that potentially encode functionally different isoforms.

Transcriptional regulation by microRNAs (miRNAs) involves complementary base-pairing at target sites on mRNAs, yielding complex secondary structures. Here we introduce an efficient computational approach and software (FASTH) for genome-scale prediction of miRNA target sites based on minimizing the free energy of duplex structure. We apply our approach to identify miRNA target sites in the human and mouse transcriptomes. Our results show that short sequence motifs in the 5′ end of miRNAs frequently match mRNAs perfectly, not only at validated target sites but additionally at many other, energetically favourable sites. High-quality matching regions are abundant and occur at similar frequencies in all mRNA regions, not only the 3′UTR. About one-third of potential miRNA target sites are reassigned to different mRNA regions, or gained or lost altogether, among different transcript isoforms from the same gene. Many potential miRNA target sites predicted in human are not found in mouse, and vice-versa, but among those that do occur in orthologous human and mouse mRNAs most are situated in corresponding mRNA regions, i.e. these sites are themselves orthologous. Using a luciferase assay in HEK293 cells, we validate four of six predicted miRNA-mRNA interactions, with the mRNA level reduced by an average of 73%. We demonstrate that a thermodynamically based computational approach to prediction of miRNA binding sites on mRNAs can be scaled to analyse complete mammalian transcriptome datasets. These results confirm and extend the scope of miRNA-mediated species- and transcript-specific regulation in different cell types, tissues and developmental conditions.

Genetic variation modulates gene expression transcriptionally or post-transcriptionally, and can profoundly alter an individual's phenotype. Measuring allelic differential expression at heterozygous loci within an individual, a phenomenon called allele-specific expression (ASE), can assist in identifying such factors. Massively parallel DNA and RNA sequencing and advances in bioinformatic methodologies provide an outstanding opportunity to measure ASE genome-wide. In this study, matched DNA and RNA sequencing, genotyping arrays and computationally phased haplotypes were integrated to comprehensively and conservatively quantify ASE in a single human brain and liver tissue sample. We describe a methodological evaluation and assessment of common bioinformatic steps for ASE quantification, and recommend a robust approach to accurately measure SNP, gene and isoform ASE through the use of personalized haplotype genome alignment, strict alignment quality control and intragenic SNP aggregation. Our results indicate that accurate ASE quantification requires careful bioinformatic analyses and is adversely affected by sample specific alignment confounders and random sampling even at moderate sequence depths. We identified multiple known and several novel ASE genes in liver, including WDR72, DSP and UBD, as well as genes that contained ASE SNPs with imbalance direction discordant with haplotype phase, explainable by annotated transcript structure, suggesting isoform derived ASE. The methods evaluated in this study will be of use to researchers performing highly conservative quantification of ASE, and the genes and isoforms identified as ASE of interest to researchers studying those loci.

MicroRNAs (miRNAs) are predominantly negative regulators of gene expression that act through the RNA-induced Silencing Complex (RISC) to suppress the translation of protein coding mRNAs. Despite intense study of these regulatory molecules, the specific molecular functions of most miRNAs remain unknown, largely due to the challenge of accurately identifying miRNA targets. Reporter gene assays can determine direct interactions, but are laborious and do not scale to genome-wide screens. Genomic scale methods such as HITS-CLIP do not preserve the direct interactions, and rely on computationally derived predictions of interactions that are plagued by high false positive rates. Here we describe a protocol for the isolation of direct targets of a mature miRNA, using synthetic biotinylated miRNA duplexes. This approach allows sensitive and specific detection of miRNA-mRNA interactions, isolating high quality mRNA suitable for analysis by microarray or RNAseq.

ABSTRACT Nullbasic is a derivative of the HIV-1 transactivator of transcription (Tat) protein that strongly inhibits HIV-1 replication in lymphocytes. Here we show that lentiviral vectors that constitutively express a Nullbasic-ZsGreen1 (NB-ZSG1) fusion protein by the eEF1α promoter led to robust long-term inhibition of HIV-1 replication in Jurkat cells. Although Jurkat-NB-ZSG1 cells were infected by HIV-1, no virus production could be detected and addition of phorbol ester 12-myristate 13-acetate (PMA) and JQ1 had no effect, while suberanilohydroxamic acid (SAHA) modestly stimulated virus production but at levels 300-fold lower than those seen in HIV-1-infected Jurkat-ZSG1 cells. Virus replication was not recovered by coculture of HIV-1-infected Jurkat-NB-ZSG1 cells with uninfected Jurkat cells. Latently infected Jurkat latent 6.3 and ACH2 cells treated with latency-reversing agents produced measurable viral capsid (CA), but little or none was made when they expressed NB-ZSG1. When Jurkat cells chronically infected with HIV-1 were transduced with lentiviral virus-like particles conveying NB-ZSG1, a &gt;3-log reduction in CA production was observed. Addition of PMA increased virus CA production but at levels 500-fold lower than those seen in nontransduced Jurkat cells. Transcriptome sequencing analysis confirmed that HIV-1 mRNA was strongly inhibited by NB-ZSG1 but indicated that full-length viral mRNA was made. Analysis of HIV-1-infected Jurkat cells expressing NB-ZSG1 by chromatin immunoprecipitation assays indicated that recruitment of RNA polymerase II (RNAPII) and histone 3 lysine 9 acetylation were inhibited. The reduction of HIV-1 promoter-associated RNAPII and epigenetic changes in viral nucleosomes indicate that Nullbasic can inhibit HIV-1 replication by enforcing viral silencing in cells. IMPORTANCE HIV-1 infection is effectively controlled by antiviral therapy that inhibits virus replication and reduces measurable viral loads in patients below detectable levels. However, therapy interruption leads to viral rebound due to latently infected cells that serve as a source of continued viral infection. Interest in strategies leading to a functional cure of HIV infection by permanent viral suppression, which may be achievable, is growing. Here we show that a mutant form of the HIV-1 Tat protein, referred to as Nullbasic, can inhibit HIV-1 transcription in infected Jurkat T cell to undetectable levels. Analysis shows that Nullbasic alters the epigenetic state of the HIV-1 long terminal repeat promoter, inhibiting its association with RNA polymerase II. This study indicates that key cellular proteins and pathways targeted here can silence HIV-1 transcription. Further elucidation could lead to functional-cure strategies by suppression of HIV transcription, which may be achievable by a pharmacological method.

Key Points Combined loss of Ssb1/Ssb2 induces rapid lethality due to replication stress–associated loss of hematopoietic stem and progenitor cells. Functionally, loss of Ssb1/Ssb2 activates p53 and IFN pathways, causing enforced cell cycling in quiescent HSPCs and apoptotic cell loss.

Accurate and complete mapping of short-read sequencing to a reference genome greatly enhances the discovery of biological results and improves statistical predictions. We recently presented RNA-MATE, a pipeline for the recursive mapping of RNA-Seq datasets. With the rapid increase in genome re-sequencing projects, progression of available mapping software and the evolution of file formats, we now present X-MATE, an updated version of RNA-MATE, capable of mapping both RNA-Seq and DNA datasets and with improved performance, output file formats, configuration files, and flexibility in core mapping software.Executables, source code, junction libraries, test data and results and the user manual are available from http://grimmond.imb.uq.edu.au/X-MATE/.

Ultraviolet radiation‐induced DNA mutations are a primary environmental driver of melanoma. The reason for this very high level of unrepaired DNA lesions leading to these mutations is still poorly understood. The primary DNA repair mechanism for UV‐induced lesions, that is, the nucleotide excision repair pathway, appears intact in most melanomas. We have previously reported a postreplication repair mechanism that is commonly defective in melanoma cell lines. Here we have used a genome‐wide approach to identify the components of this postreplication repair mechanism. We have used differential transcript polysome loading to identify transcripts that are associated with UV response, and then functionally assessed these to identify novel components of this repair and cell cycle checkpoint network. We have identified multiple interaction nodes, including global genomic nucleotide excision repair and homologous recombination repair, and previously unexpected MASTL pathway, as components of the response. Finally, we have used bioinformatics to assess the contribution of dysregulated expression of these pathways to the UV signature mutation load of a large melanoma cohort. We show that dysregulation of the pathway, especially the DNA damage repair components, are significant contributors to UV mutation load, and that dysregulation of the MASTL pathway appears to be a significant contributor to high UV signature mutation load.

A new algorithm, PfAGSS, for predicting 3' splice sites in Plasmodium falciparum genomic sequences is described. Application of this program to the published P. falciparum chromosome 2 and 3 data suggests that existing programs result in a high error rate in assigning 3' intron boundaries.

Next-generation sequencing is revolutionizing the way in which genomic-scale biological research is performed, and its effects are beginning to be translated medically. Large-scale international collaborations for the comprehensive sequencing of the genome, epigenome, and transcriptomes of cancers and corresponding 'normal' (germ-line) DNA are heralding the start of personalized medical genomics. The promise of eliminating conjecture when determining treatment approaches is certainly appealing for both patients and clinicians; however, several major issues must be resolved before next-generation sequencing will be adopted as a routine clinical tool for patients. This feature review explores the clinical potential and challenges of studying cancer genomes for personalized medical genomics.

Nature 516, 198–206 (2014); doi:10.1038/nature14046 In this Article, the address listed as Nicole Cloonan’s present address (QIMR Berghofer Medical Research Institute, Queensland 4006, Australia) should have been listed as her other affiliation with superscript 22, because the work she did was splitequally between the two institutions.

Nature 516, 192–197 (2014); doi:10.1038/nature14047 In this Article, the address listed as Nicole Cloonan’s present address (QIMR Berghofer Medical Research Institute, Queensland 4006, Australia) should have been listed as her other affiliation with superscript 16, because the work she did was splitequally between the two institutions.

Transcriptional regulation by microRNAs (miRNAs) involves complementary base-pairing at target sites on mRNAs, yielding complex secondary structures.Here we introduce an efficient computational approach and software (FASTH) for genomescale prediction of miRNA target sites based on minimizing the free energy of duplex structure.We apply our approach to identify miRNA target sites in the human and mouse transcriptomes.Our results show that short sequence motifs in the 59 end of miRNAs frequently match mRNAs perfectly, not only at validated target sites but additionally at many other, energetically favourable sites.High-quality matching regions are abundant and occur at similar frequencies in all mRNA regions, not only the 39UTR.About one-third of potential miRNA target sites are reassigned to different mRNA regions, or gained or lost altogether, among different transcript isoforms from the same gene.Many potential miRNA target sites predicted in human are not found in mouse, and vice-versa, but among those that do occur in orthologous human and mouse mRNAs most are situated in corresponding mRNA regions, i.e. these sites are themselves orthologous.Using a luciferase assay in HEK293 cells, we validate four of six predicted miRNA-mRNA interactions, with the mRNA level reduced by an average of 73%.We demonstrate that a thermodynamically based computational approach to prediction of miRNA binding sites on mRNAs can be scaled to analyse complete mammalian transcriptome datasets.These results confirm and extend the scope of miRNAmediated species-and transcript-specific regulation in different cell types, tissues and developmental conditions.

Nat. Genet. 41, 572–578 (2009); published online 19 April 2009; corrected after print 26 June 2009 In the version of this article initially published, some author affiliations were incorrectly stated. The error has been corrected in the HTML and PDF versions of the article.

Somatic rearrangements, which are commonly found in human cancer genomes, contribute to the progression and maintenance of cancers. Conventionally, the verification of somatic rearrangements comprises many manual steps and Sanger sequencing. This is labor intensive when verifying a large number of rearrangements in a large cohort. To increase the verification throughput, we devised a high-throughput workflow that utilizes benchtop next-generation sequencing and in-house bioinformatics tools to link the laboratory processes. In the proposed workflow, primers are automatically designed. PCR and an optional gel electrophoresis step to confirm the somatic nature of the rearrangements are performed. PCR products of somatic events are pooled for Ion Torrent PGM and/or Illumina MiSeq sequencing, the resulting sequence reads are assembled into consensus contigs by a consensus assembler, and an automated BLAT is used to resolve the breakpoints to base level. We compared sequences and breakpoints of verified somatic rearrangements between the conventional and high-throughput workflow. The results showed that next-generation sequencing methods are comparable to conventional Sanger sequencing. The identified breakpoints obtained from next-generation sequencing methods were highly accurate and reproducible. Furthermore, the proposed workflow allows hundreds of events to be processed in a shorter time frame compared with the conventional workflow.

Abstract Myeloid malignancies arise following the sequential acquisition of somatic mutations within hematopoietic stem and progenitor cells (HSPC). JAK2V617F is commonly found in myeloproliferative neoplasms (MPN) such as polycythemia vera, essential thrombocythemia and myelofibrosis. While other mutations (e.g. TET2, DNMT3A) have been found to co-occur in MPN HSPC, it remains unclear how they impact disease biology or progression from early stage disease (i.e. polycythemia or essential thrombocythemia) to advanced stage disease such as myelofibrosis or acute myeloid leukemia. DNMT3A methylates cytosine rich DNA residues (known as CpG islands, and often found in promoters of genes) leading to transcriptional repression. DNMT3A is also recurrently mutated at relatively low frequency in polycythemia vera (5-7%) but mutations are more common in advanced MPN (approximately 15% of MF and 17% of AML, Stegelmann et al. Leukemia 2011; Abdel-Wahab et al. Leukaemia 2011). These mutations are found in the methyltransferase domain and cluster around arginine 882 (e.g. R882H), resulting in loss of DNA binding and reduced catalytic activity. We used CRISPR-Cas9 gene editing technology to disrupt Dnmt3a function in mouse HSPC and assessed for cooperativity together with a conditional, knockin Jak2V617F allele. Jak2V617F/∆Dnmt3a-Cas9 but not Jak2V617F/Cas9 controls demonstrated increased HSPC self-renewal and proliferation properties in vitro as evidenced by serial replating in methylcellulose (&gt;5 weeks) and increased colony forming unit capacity. Flow-cytometry analysis of Jak2V617F/∆Dnmt3a-Cas9 revealed enrichment in LKS+ (Lin-Sca-1highKithigh) cells 5 weeks after CRISPR-Cas9disruption of Dnmt3a, and this was associated with increased expression of stemness markers Kit and Cd34 in Jak2V617F/∆Dnmt3a-Cas9 cells. RNAseq was performed on early (week 1, P1) and late culture HSPC (week 5, P5) from Jak2V617F-Cas9 (P1 only) and Jak2V617F/∆Dnmt3a-Cas9 (P1, P5). This confirmed deletion of Dnmt3a in Jak2V617F/∆Dnmt3a-Cas9 but not in Jak2V617F/Cas9 controls. Transcriptional upregulation of Kit and Cd34 were confirmed, as well as other key stem cell genes such as Erg and Angpt1 in Jak2V617F/∆Dnmt3a-Cas9 P5. We observed denovo expression of imprinted genes Igf2 and H19 in Jak2V617F/∆Dnmt3a P5, suggesting impaired DNA methylation in this group. Jak2V617F/∆Dnmt3a-Cas9 P5 were significantly enriched for transcriptional pathways controlling cell cycle progression, oncogenic signatures, and DNA damage. Conversely, Jak2V617F/Cas9 controls were enriched for myeloid differentiation and normal progenitor cell signatures. To assess the effect of Dnmt3a loss on Jak2V617F driven MPN, we transplanted Jak2V617F/∆Dnmt3a-Cas9 or Jak2V617F/Cas9 LKS+ into irradiated B6 recipients. Recipients of Jak2V617F/Cas9 LKS+ developed early stage MPN reminiscent of polycythemia vera with high hemoglobin, white cell count and platelets and was sustained &gt;32 weeks. In contrast, Jak2V617F/∆Dnmt3-Cas9 recipients exhibited a biphasic disease, reminiscent of human myelofibrosis. At 8 weeks, Jak2V617F/∆Dnmt3-Cas9 showed panmyelosis with thrombocytosis (1.38x106/µl vs. 1.14x106/µl controls, p=0.057). However, by 32 weeks, this mice became severely pancytopenic with progressive bone marrow failure (Hemoglobin 121g/L vs. 210g/L controls, p =0.0011; platelets 0.338x106/µl vs. 1.343x106/µl controls, p &lt;0.0001). Jak2V617F/∆Dnmt3-Cas9 mice exhibited extreme splenomegaly associated with reticulin fibrosis and the accumulation of myeloid cells. Bone marrow histology of Jak2V617F/∆Dnmt3-Cas9 revealed osteosclerosis and disorganized architecture and a dense fibrocellular infiltrate and reticulin fibrosis. Flow cytometry revealed impaired erythropoiesis and blocked differentiation. AML was not seen. These data demonstrate new evidence linking loss of Dnmt3a with acquisition of self-renewal in combination with constitutively active Jak2V617F. Importantly, in vivo loss of Dnmt3a accelerates or induces myelofibrotic transformation of the underlying MPN. This work provides new understanding to the factors that promote advanced disease in MPN. Ultimately, such knowledge has the potential to inform the development of novel targeted therapeutic approaches for the treatment of transformed MPN, a highly chemorefractory disease associated with extremely poor prognosis in patients. Disclosures Lane: Janssen: Other: i have done consulting (once) for janssen..

Hematopoietic stem and progenitor cells (HSPCs) are vulnerable to endogenous damage and defects in DNA repair can limit their function. The 2 single-stranded DNA (ssDNA) binding proteins SSB1 and SSB2 are crucial regulators of the DNA damage response; however, their overlapping roles during normal physiology are incompletely understood. We generated mice in which both Ssb1 and Ssb2 were constitutively or conditionally deleted. Constitutive Ssb1/Ssb2 double knockout (DKO) caused early embryonic lethality, whereas conditional Ssb1/Ssb2 double knockout (cDKO) in adult mice resulted in acute lethality due to bone marrow failure and intestinal atrophy featuring stem and progenitor cell depletion, a phenotype unexpected from the previously reported single knockout models of Ssb1 or Ssb2. Mechanistically, cDKO HSPCs showed altered replication fork dynamics, massive accumulation of DNA damage, genome-wide double-strand breaks enriched at Ssb-binding regions and CpG islands, together with the accumulation of R-loops and cytosolic ssDNA. Transcriptional profiling of cDKO HSPCs revealed the activation of p53 and interferon (IFN) pathways, which enforced cell cycling in quiescent HSPCs, resulting in their apoptotic death. The rapid cell death phenotype was reproducible in in vitro cultured cDKO-hematopoietic stem cells, which were significantly rescued by nucleotide supplementation or after depletion of p53. Collectively, Ssb1 and Ssb2 control crucial aspects of HSPC function, including proliferation and survival in vivo by resolving replicative stress to maintain genomic stability.

JAK2V617F is present in the majority of patients with myeloproliferative neoplasm (MPN), however, it is unclear how co-occurring mutations in epigenetic regulators (e.g. DNMT3A) in MPN impact disease biology. DNMT3A methylates DNA at cytosine residues in enriched CpG region, and drives promoter hypermethylation. DNMT3A is mutated in advanced MPN (15% of MF and 17% of AML). These mutations sit in the methyltransferase domain and result in reduced activity. To determine the functional effects of Dnmt3a loss in MPN, we used CRISPR-Cas9 technology to edit Dnmt3a function in Jak2V617F murine hematopoietic stem and progenitor cells (HSPC). Jak2V617F LKS+ (Lin-Sca-1highKithigh) HSPCs were infected with a lentivirus encoding Cas9 and sgRNA targeting Dnmt3a (Jak2V617F/ΔDnmt3a) or non-targeting controls (Jak2V617F). In vitro CFU from Jak2V617F/ΔDnmt3a had enhanced serial replating and increased expression of markers Kit and Cd34. CFU RNAseq confirmed Dnmt3a deletion in Jak2V617F/ΔDnmt3a, as well as transcriptional upregulation of key stemness genes, and de novo expression of imprinted genes Igf2 and H19, together with decreased H19 promoter methylation in Jak2V617F/ΔDnmt3a compared to controls. To assess the functional effect of Dnmt3a loss on Jak2V617F driven MPN, we transplanted Jak2V617F/ΔDnmt3a or Jak2V617F LKS+ into irradiated recipients. Jak2V617F/ΔDnmt3a became severely pancytopenic with bone marrow failure, fibrosis and the accumulation of myeloid cells. Jak2V617F/ΔDnmt3a bone marrow showed osteosclerosis, disorganized architecture with dense fibrocellular infiltrate. RNAseq on Jak2V617F/ΔDnmt3a LKS showed strong enrichment in gene sets previously annotated by Dnmt3a-/- in normal HSPC (Challen, Nat Gen 2011) and oncogenic Dnmt3a R878H (Guryanova,Nat Med 2016), together with LT-HSC signatures and PRC2 deregulation. This novel model of mutant Dnmt3a, achieved through in vivo CRISPR-Cas9 editing, links loss of Dnmt3a with acquisition of self-renewal and disease progression in Jak2V617F HSPC. Such knowledge has the potential to inform the development of targeted therapeutic approaches in transformed MPN, a highly chemorefractory disease associated with poor prognosis.

Abstract Abstract 1462 Poster Board I-485 Klf1/Eklf regulates a diverse suite of genes to direct erythroid cell differentiation from bi-potent progenitors. To determine the local cis-regulatory contexts and transcription factor networks in which Klf1 works, we performed Klf1 ChIP-seq using the SOLiD deep sequencing platform. We mapped more than 10 million unique 35mer tags and found ∼1500 sites in the genome of primary fetal liver erythroid cells are occupied by endogenous Klf1. Many reside within well characterised erythroid gene promoters (e.g. b-globin) or enhancers (e.g. E2f2 intron 1), but some are &gt;100kb from any known gene. We tested a number of Klf1 bound promoter and intragenic sites for activity in erythroid cell lines and zebrafish. Our data suggests Klf1 directly regulates most aspects of terminal erythroid differentiation including synthesis of the hemoglobin tetramer, construction of a deformable red cell membrane and cytoskeleton, bimodal regulation of proliferation, and co-ordination of anti-apoptosis and enucleation pathways. Additionally, we suggest new mechanisms for Klf1 co-operation with other transcription factors such as those of the gata, ets and myb families based on over-representation and spatial constraints of their binding motifs in the vicinity of Klf1-bound promoters and enhancers. Finally, we have identified a group of ∼100 Klf1-occupied sites in fetal liver which overlap with Klf4-occupied sites in ES cells defined by Klf4 ChIP-seq. These sites are associated with genes controlling the cell cycle and proliferation and are Klf4-dependent in skin, gut and ES cells, suggesting a global paradigm for Klfs as regulators of differentiation in many, if not all, cell types. Disclosures No relevant conflicts of interest to declare.

Stress activated protein kinase (SAPK) interacting protein 1 (Sin1) is a member of a recently characterized gene family, conserved from yeast to humans. The gene copy number is strictly conserved (one Sin1 gene per genome), and the protein may be expressed ubiquitously in mammalian tissues. The Sin1 family has been implicated in several different signal transduction pathways. Originally identified as a partial cDNA and candidate Ras inhibitor, recent functional studies have revealed interactions with an interferon (IFN) receptor subunit (IFNAR2), and the SAPK JNK. Interactions have also been described between the yeast orthologues and the phosphatidylinositol kinase TOR2. Collectively, these data suggest that Sin1 has an important cellular role, and this study has investigated possible functions for this protein. As human Sin1 proteins have no paralogues within the genome, secondary structure homology was used to identify major domains within the protein. Four major domains within human Sin1 were deduced: an N-terminal domain containing a functional nuclear localization signal, a functional nuclear export signal, and a coiledcoil region; the conserved region in the middle that is likely to be a ubiquitin-like I²-grasp protein binding domain; a Ras binding domain; and a pleckstrin homology-like domain that targets Sin1 to the plasma membrane and lipid rafts in vivo. Full and partial length EGFP constructs were used to examine the localization of human Sin1, and several isoforms derived from alternative splicing. All isoforms localized to the nucleus and nucleolus. Beyond this, Sin1I± and Sin1I’ had cytoplasmic staining, while Sin1 and Sin1I² were also found at the plasma membrane and lipid rafts. Both the N-terminal domain and the conserved region in the middle were found to contribute to nuclear localization. Comparative genomic analysis between human, mouse, rat, dog, and chicken Sin1 genes revealed a number of conserved intronic regions, and the putative functions of these were predicted. Additionally, a putative promoter module within a CpG island and encompassing the transcription start site was predicted in all species. The human CpG island was found to have promoter activity in HEK293 cells. Using bioinformatics, genes that may be co-regulated with Sin1 were identified. These genes contained the Sin1 promoter module, and were found to co-express in large scale gene expression studies. Most of these genes were directly involved in the cellular response to pathogen infection, suggesting a conserved role for Sin1 in this pathway. Key biochemical functions of the Sin1 proteins were also identified, including the ability of Sin1 proteins to form dimers, and the ability of over-expressed Sin1 to induce apoptosis (mediated through the conserved region in the middle). Additionally, endogenous Sin1 protein levels were found to change following serum deprivation and hypoosmotic stress. Together, these studies have provided significant insight into the cellular role of Sin1, suggesting a role in inducing apoptosis as part of the IFN response to viral infection. The biological significance of the Sin1 proteins is discussed in the context of their predicted functions and the evolution of the protein family.