Natalie A. Twine

Whole transcriptome sequencing reveals gene expression and splicing differences in brain regions affected by Alzheimer's disease.

Recent studies strongly indicate that aberrations in the control of gene expression might contribute to the initiation and progression of Alzheimers disease (AD). In particular, alternative splicing has been suggested to play a role in spontaneous cases of AD. Previous transcriptome profiling of AD models and patient samples using microarrays delivered conflicting results. This study provides, for the first time, transcriptomic analysis for distinct regions of the AD brain using RNA-Seq next-generation sequencing technology. Illumina RNA-Seq analysis was used to survey transcriptome profiles from total brain, frontal and temporal lobe of healthy and AD post-mortem tissue. We quantified gene expression levels, splicing isoforms and alternative transcript start sites. Gene Ontology term enrichment analysis revealed an overrepresentation of genes associated with a neurons cytological structure and synapse function in AD brain samples. Analysis of the temporal lobe with the Cufflinks tool revealed that transcriptional isoforms of the apolipoprotein E gene, APOE-001, -002 and -005, are under the control of different promoters in normal and AD brain tissue. We also observed differing expression levels of APOE-001 and -002 splice variants in the AD temporal lobe. Our results indicate that alternative splicing and promoter usage of the APOE gene in AD brain tissue might reflect the progression of neurodegeneration.

454 pyrosequencing-based analysis of gene expression profiles in the amphipod Melita plumulosa: transcriptome assembly and toxicant induced changes.

Next generation sequencing using Roches 454 pyrosequencing platform can be used to generate genomic information for non-model organisms, although there are bioinformatic challenges associated with these studies. These challenges are compounded by a lack of a standardized protocol to either assemble data or to evaluate the quality of a de novo transcriptome. This study presents an assembly of the control and toxicant responsive transcriptome of Melita plumulosa, an Australian amphipod commonly used in ecotoxicological studies. RNA was harvested from control amphipods, juvenile amphipods, and from amphipods exposed to either metal or diesel contaminated sediments. This RNA was used as the basis for a 454 based transcriptome sequencing effort. Sequencing generated 1.3 million reads from control, juvenile, metal-exposed and diesel-exposed amphipods. Different read filtering and assembly protocols were evaluated to generate an assembly that (i) had an optimal number of contigs; (ii) had long contigs; (iii) contained a suitable representation of conserved genes; and (iv) had long ortholog alignment lengths relative to the length of each contig. A final assembly, generated using fixed-length trimming based on the sequence quality scores, followed by assembly using the MIRA algorithm, produced the best results. The 26,625 contigs generated via this approach were annotated using Blast2GO, and the differential expression between treatments and control was determined by mapping with BWA followed by DESeq. Although the mapping generated low coverage, many differentially expressed contigs, including some with known developmental or toxicological function, were identified. This study demonstrated that 454 pyrosequencing is an effective means of generating reference transcriptome information for organisms, such as the amphipod M. plumulosa, that have no genomic information available in databases or in closely related sequenced species. It also demonstrated how optimization of read filtering protocols and assembly approaches changes the utility of results obtained from next generation sequencing studies, and establishes criteria to determine the quality of a de novo assembly in species lacking a reference genome. This new transcriptomic knowledge provides the genomic foundation for the creation of microarray and qPCR assays, serving as a reference transcriptome in future RNAseq studies, and allowing both the biology and ecotoxicology of this organism to be better understood. This approach will allow genomics-based methodology to be applied to a wider range of environmentally relevant species.

The CACCC-binding protein KLF3/BKLF represses a subset of KLF1/EKLF target genes and is required for proper erythroid maturation in vivo.

ABSTRACT The CACCC-box binding protein erythroid Krüppel-like factor (EKLF/KLF1) is a master regulator that directs the expression of many important erythroid genes. We have previously shown that EKLF drives transcription of the gene for a second KLF, basic Krüppel-like factor, or KLF3. We have now tested the in vivo role of KLF3 in erythroid cells by examining Klf3 knockout mice. KLF3-deficient adults exhibit a mild compensated anemia, including enlarged spleens, increased red pulp, and a higher percentage of erythroid progenitors, together with elevated reticulocytes and abnormal erythrocytes in the peripheral blood. Impaired erythroid maturation is also observed in the fetal liver. We have found that KLF3 levels rise as erythroid cells mature to become TER119+. Consistent with this, microarray analysis of both TER119− and TER119+ erythroid populations revealed that KLF3 is most critical at the later stages of erythroid maturation and is indeed primarily a transcriptional repressor. Notably, many of the genes repressed by KLF3 are also known to be activated by EKLF. However, the majority of these are not currently recognized as erythroid-cell-specific genes. These results reveal the molecular and physiological function of KLF3, defining it as a feedback repressor that counters the activity of EKLF at selected target genes to achieve normal erythropoiesis.

Generation of Mice Deficient in both KLF3/BKLF and KLF8 Reveals a Genetic Interaction and a Role for These Factors in Embryonic Globin Gene Silencing

ABSTRACT Krüppel-like factors 3 and 8 (KLF3 and KLF8) are highly related transcriptional regulators that bind to similar sequences of DNA. We have previously shown that in erythroid cells there is a regulatory hierarchy within the KLF family, whereby KLF1 drives the expression of both the Klf3 and Klf8 genes and KLF3 in turn represses Klf8 expression. While the erythroid roles of KLF1 and KLF3 have been explored, the contribution of KLF8 to this regulatory network has been unknown. To investigate this, we have generated a mouse model with disrupted KLF8 expression. While these mice are viable, albeit with a reduced life span, mice lacking both KLF3 and KLF8 die at around embryonic day 14.5 (E14.5), indicative of a genetic interaction between these two factors. In the fetal liver, Klf3 Klf8 double mutant embryos exhibit greater dysregulation of gene expression than either of the two single mutants. In particular, we observe derepression of embryonic, but not adult, globin expression. Taken together, these results suggest that KLF3 and KLF8 have overlapping roles in vivo and participate in the silencing of embryonic globin expression during development.

Proteomic Validation of Transcript Isoforms, Including Those Assembled from RNA-Seq Data

Human proteome analysis now requires an understanding of protein isoforms. We recently published the PG Nexus pipeline, which facilitates high confidence validation of exons and splice junctions by integrating genomics and proteomics data. Here we comprehensively explore how RNA-seq transcriptomics data, and proteomic analysis of the same sample, can identify protein isoforms. RNA-seq data from human mesenchymal (hMSC) stem cells were analyzed with our new TranscriptCoder tool to generate a database of protein isoform sequences. MS/MS data from matching hMSC samples were then matched against the TranscriptCoder-derived database, along with Ensembl and the neXtProt database. Querying the TranscriptCoder-derived or Ensembl database could unambiguously identify ∼450 protein isoforms, with isoform-specific proteotypic peptides, including candidate hMSC-specific isoforms for the genes DPYSL2 and FXR1. Where isoform-specific peptides did not exist, groups of nonisoform-specific proteotypic peptides could specifically identify many isoforms. In both the above cases, isoforms will be detectable with targeted MS/MS assays. Unfortunately, our analysis also revealed that some isoforms will be difficult to identify unambiguously as they do not have peptides that are sufficiently distinguishing. We covisualize mRNA isoforms and peptides in a genome browser to illustrate the above situations. Mass spectrometry data is available via ProteomeXchange (PXD001449).

RNA-Seq analysis of the toxicant-induced transcriptome of the marine diatom, Ceratoneis closterium.

Diatoms are of enormous ecological importance as they account for as much as 20% of global primary production, yet they are still understudied from a genomic perspective. The benthic diatom Ceratoneis closterium is well-characterized from an ecotoxicological perspective including its use in ecotoxicological risk assessments and investigating the mode-of-action of metal toxicity. However, this organism has little sequence information available. In this study, 454 pyrosequencing of the stressor-responsive transcriptome was undertaken. These transcripts could be used to characterize general physiological processes such as photosynthesis and respiration, as well as to enable a description of the ecotoxicogenomic responses of this organism. After a 96 h exposure to the concentration of toxicant that inhibited growth rate by 10% (IC10) for the following common coastal contaminants: ammonia, copper, crude oil and simazine (a photosystem II inhibiting herbicide), diatom cells were harvested for RNA extraction and their transcriptomes characterized via 454 pyrosequencing. This resulted in 1.25 million reads, which were assembled into 4768 contigs, when contigs encoding rRNA were removed. More than 80% of the remaining contigs had an ortholog in the BLASTx protein databases. These contigs represented 1660 unique transcripts. The role of these transcripts in stress response, as well as photosynthesis and respiration is discussed. Overall, this study greatly enhances the genomic information available for this important taxonomic group.

Next generation sequence analysis of the transcriptome of Sydney rock oysters (Saccostrea glomerata) exposed to a range of environmental stressors

Sydney rock oysters (Saccostrea glomerata) were exposed to environmental stressors at contaminated field sites or in a controlled laboratory setting. RNA seq transcriptome data were generated for the gill and digestive gland using Roches 454 pyrosequencing technology. 28,685 contigs were de novo assembled which encoded 11,671 different protein products. The data will act as a reference for future studies in ecology, immunology and environmental toxicology.

Sequencing of hippocampal and cerebellar transcriptomes provides new insights into the complexity of gene regulation in the human brain.

The hippocampus and cerebellum represent anatomically and functionally distinct parts of the human brain. The RNA-Seq technique makes it possible to investigate the human transcriptome with unprecedented resolution, allowing identification of differential mRNA splicing and promoter usage on a genome-wide scale. We undertook whole-mRNA sequencing of samples from the human hippocampus and cerebellum. A bioinformatic analysis revealed distinct expression patterns of genes related to the molecular physiology of neurons and glial cells. Upregulated genes in hippocampal tissue included serpin peptidase inhibitor, clade A (SERPINA3), lymphocyte antigen 6 complex, locus H (LY6H) and transthyretin (TTR). In cerebellum, the cerebellin 3 precursor (CLBN3) and Zic family member 4 (ZIC4) genes were significantly upregulated. These changes were validated in independent donor samples by qRT-PCR. The hippocampus and the cerebellum showed striking differences in splicing patterns and promoter usage. A notable example of this was the gene for NGFI-A binding protein 2 (NAB2), which displayed tissue-specific isoforms which may affect its function as a transcriptional repressor.

Molecular Phenotyping of Telomerized Human Bone Marrow Skeletal Stem Cells Reveals a Genetic Program of Enhanced Proliferation and Maintenance of Differentiation Responses: MOLECULAR PHENOTYPING OF TELOMERIZED HMSC

Long‐term in vitro expansion of bone marrow stromal (skeletal) stem cells (also known as human mesenchymal stem cells [hMSC]) is associated with replicative senescence and impaired functions. We have previously reported that telomerization of hMSC through hTERT overexpression led to bypassing a replicative senescence phenotype and improved in vitro and in vivo functions. However, the molecular consequence of telomerization is poorly characterized. Thus, we compared the molecular phenotype of a well‐studied telomerized hMSC (hMSC‐TERT) cell line with primary hMSC. At a cellular level, both cell populations exhibited strong concordance for the known hMSC CD markers, similar responses to osteoblast (OB) differentiation induction, and formed heterotopic bone in vivo. Overall gene expression was highly correlated between both cell types with an average Pearsons correlation coefficient (R2) between the gene expression of all primary hMSC and all hMSC‐TERT samples of 0.95 (range 0.93–0.96). Quantitative analysis of gene expression of CD markers, OB cell markers, and transcription factors (TF) showed a high degree of similarity between the two cell populations (72%, 77%, and 81%, respectively). The hMSC‐TERT population was enriched mainly for genes associated with cell cycle and cell cycle signaling when compared with primary hMSC. Other enrichment was observed for genes involved in cell adhesion and skeletal system development and immune response pathways. Interestingly, hMSC‐TERT shared a telomerization signature with upregulation of cancer/testis antigens, MAGE, and PAGE genes. Our data demonstrate that the enhanced biological characteristics of hMSC after telomerization are mainly due to enhanced expression of cell proliferation genes, whereas gene expression responses to differentiation are maintained.

Transcription factor ZNF25 is associated with osteoblast differentiation of human skeletal stem cells

BackgroundThe differentiation of human bone marrow derived skeletal stem cells (known as human bone marrow stromal or mesenchymal stem cells, hMSCs) into osteoblasts involves the activation of a small number of well-described transcription factors. To identify additional osteoblastic transcription factors, we studied gene expression of hMSCs during ex vivo osteoblast differentiation.ResultsClustering of gene expression, and literature investigation, revealed three transcription factors of interest – ZNF25, ZNF608 and ZBTB38. siRNA knockdown of ZNF25 resulted in significant suppression of alkaline phosphatase (ALP) activity. This effect was not present for ZNF608 and ZBTB38. To identify possible target genes of ZNF25, we analyzed gene expression following ZNF25 siRNA knockdown. This revealed a 23-fold upregulation of matrix metallopeptidase 1 and an 18-fold upregulation of leucine-rich repeat containing G protein-coupled receptor 5 and RAN-binding protein 3-like. We also observed enrichment in extracellular matrix organization, skeletal system development and regulation of ossification in the entire upregulated set of genes. Consistent with its function as a transcription factor during osteoblast differentiation of hMSC, we showed that the ZNF25 protein exhibits nuclear localization and is expressed in osteoblastic and osteocytic cells in vivo. ZNF25 is conserved in tetrapod vertebrates and contains a KRAB (Krueppel-associated box) transcriptional repressor domain.ConclusionsThis study shows that the uncharacterized transcription factor, ZNF25, is associated with differentiation of hMSC to osteoblasts.