Impact of alternative initiation, splicing, and termination on the diversity of the mRNA transcripts encoded by the mouse transcriptome

We analyzed the FANTOM2 clone set of 60,770 RIKEN full-length mouse cDNA sequences and 44,122 public mRNA sequences. We developed a new computational procedure to identify and classify the forms of splice variation evident in this data set and organized the results into a publicly accessible database that can be used for future expression array construction, structural genomics, and analyses of the mechanism and regulation of alternative splicing. Statistical analysis shows that at least 41% and possibly as much as 60% of multiexon genes in mouse have multiple splice forms. Of the transcription units with multiple splice forms, 49% contain transcripts in which the apparent use of an alternative transcription start (stop) is accompanied by alternative splicing of the initial (terminal) exon. This implies that alternative transcription may frequently induce alternative splicing. The fact that 73% of all exons with splice variation fall within the annotated coding region indicates that most splice variation is likely to affect the protein form. Finally, we compared the set of constitutive (present in all transcripts) exons with the set of cryptic (present only in some transcripts) exons and found statistically significant differences in their length distributions, the nucleoticle distributions around their splice junctions, and the frequencies of occurrence of several short sequence motifs.

Cambial meristem dormancy in trees involves extensive remodelling of the transcriptome

The establishment of the dormant state in meristems involves considerable physiological and metabolic alterations necessary for surviving unfavourable growth conditions. However, a global molecular analysis of dormancy in meristems has been hampered by the difficulty in isolating meristem cells. We used cryosectioning to isolate purified cambial meristem cells from the woody plant Populus tremula during active growth and dormancy. These samples were used to generate meristem-specific cDNA libraries and for cDNA microarray experiments to define the global transcriptional changes underlying cambial dormancy. The results indicate a significant reduction in the complexity of the cambial transcriptome in the dormant state. Although cell division is terminated in the dormant cambium, the cell cycle machinery appears to be maintained in a skeletal state as suggested by the continued presence of transcripts for several cell cycle regulators. The downregulation of PttPIN1 and PttPIN2 transcripts explains the reduced basipetal polar auxin transport during dormancy. The induction of a member of the SINA family of ubiquitin ligases implicated in auxin signalling indicates a potential mechanism for modulation of auxin sensitivity during cambial dormancy. The metabolic alterations during dormancy are mirrored in the induction of genes involved in starch breakdown and the glyoxysomal cycle. Interestingly, the induction of RGA1 like gene suggests modification of gibberellin signalling in cambial dormancy. The induction of genes such as poplar orthologues of FIE and HAP2 indicates a potential role for these global regulators of transcription in orchestrating extensive changes in gene expression during dormancy.

In vitro and in silico analysis of signal peptides from the human blood fluke, Schistosoma mansoni

Proteins secreted by and anchored on the surfaces of parasites are in intimate contact with host tissues. The transcriptome of infective cercariae of the blood fluke, Schistosoma mansoni, was screened using signal sequence trap to isolate cDNAs encoding predicted proteins with an N-terminal signal peptide. Twenty cDNA fragments were identified, most of which contained predicted signal peptides or transmembrane regions, including a novel putative seven-transmembrane receptor and a membrane-associated mitogen-activated protein kinase. The developmental expression pattern within different life-cycle stages ranged from ubiquitous to a transcript that was highly upregulated in the cercaria. A bioinformatics-based comparison of 100 signal peptides from each of schistosomes, humans, a parasitic nematode and Escherichia coli showed that differences in the sequence composition of signal peptides, notably the residues flanking the predicted cleavage site, might account for the negative bias exhibited in the processing of schistosome signal peptides in mammalian cells. (c) 2005 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.

Developing pineapple fruit has a small transcriptome dominated by metallothionein

In a first step toward understanding the molecular basis of pineapple fruit development, a sequencing project was initiated to survey a range of expressed sequences from green unripe and yellow ripe fruit tissue. A highly abundant metallothionein transcript was identified during library construction, and was estimated to account for up to 50% of all EST library clones. Library clones with metallothionein subtracted were sequenced, and 408 unripe green and 1140 ripe yellow edited EST clone sequences were retrieved. Clone redundancy was high, with the combined 1548 clone sequences clustering into just 634 contigs comprising 191 consensus sequences and 443 singletons. Half of the EST clone sequences clustered within 13.5% and 9.3% of contigs from green unripe and yellow ripe libraries, respectively, indicating that a small subset of genes dominate the majority of the transcriptome. Furthermore, sequence cluster analysis, northern analysis, and functional classification revealed major differences between genes expressed in the unripe green and ripe yellow fruit tissues. Abundant genes identified from the green fruit include a fruit bromelain and a bromelain inhibitor. Abundant genes identified in the yellow fruit library include a MADS box gene, and several genes normally associated with protein synthesis, including homologues of ribosomal L10 and the translation factors SUI1 and eIF5A. Both the green unripe and yellow ripe libraries contained high proportions of clones associated with oxidative stress responses and the detoxification of free radicals.

Pseudo-messenger RNA: Phantoms of the transcriptome

The mammalian transcriptome harbours shadowy entities that resist classification and analysis. In analogy with pseudogenes, we define pseudo-messenger RNA to be RNA molecules that resemble protein- coding mRNA, but cannot encode full-length proteins owing to disruptions of the reading frame. Using a rigorous computational pipeline, which rules out sequencing errors, we identify 10,679 pseudo - messenger RNAs ( approximately half of which are transposonassociated) among the 102,801 FANTOM3 mouse cDNAs: just over 10% of the FANTOM3 transcriptome. These comprise not only transcribed pseudogenes, but also disrupted splice variants of otherwise protein- coding genes. Some may encode truncated proteins, only a minority of which appear subject to nonsense- mediated decay. The presence of an excess of transcripts whose only disruptions are opal stop codons suggests that there are more selenoproteins than currently estimated. We also describe compensatory frameshifts, where a segment of the gene has changed frame but remains translatable. In summary, we survey a large class of non- standard but potentially functional transcripts that are likely to encode genetic information and effect biological processes in novel ways. Many of these transcripts do not correspond cleanly to any identifiable object in the genome, implying fundamental limits to the goal of annotating all functional elements at the genome sequence level.

Transcriptomic Analysis of the Host Response and Innate Resilience to Enterotoxigenic Escherichia coli Infection in Humans.

BACKGROUND: Enterotoxigenic Escherichia coli (ETEC) is a globally prevalent cause of diarrhea. Though usually self-limited, it can be severe and debilitating. Little is known about the host transcriptional response to infection. We report the first gene expression analysis of the human host response to experimental challenge with ETEC. METHODS: We challenged 30 healthy adults with an unattenuated ETEC strain, and collected serial blood samples shortly after inoculation and daily for 8 days. We performed gene expression analysis on whole peripheral blood RNA samples from subjects in whom severe symptoms developed (n = 6) and a subset of those who remained asymptomatic (n = 6) despite shedding. RESULTS: Compared with baseline, symptomatic subjects demonstrated significantly different expression of 406 genes highlighting increased immune response and decreased protein synthesis. Compared with asymptomatic subjects, symptomatic subjects differentially expressed 254 genes primarily associated with immune response. This comparison also revealed 29 genes differentially expressed between groups at baseline, suggesting innate resilience to infection. Drug repositioning analysis identified several drug classes with potential utility in augmenting immune response or mitigating symptoms. CONCLUSIONS: There are statistically significant and biologically plausible differences in host gene expression induced by ETEC infection. Differential baseline expression of some genes may indicate resilience to infection.

Transcriptome-Wide Mapping of Pea Seed Ageing Reveals a Pivotal Role for Genes Related to Oxidative Stress and Programmed Cell Death

Understanding of seed ageing, which leads to viability loss during storage, is vital for ex situ plant conservation and agriculture alike. Yet the potential for regulation at the transcriptional level has not been fully investigated. Here, we studied the relationship between seed viability, gene expression and glutathione redox status during artificial ageing of pea (Pisum sativum) seeds. Transcriptome-wide analysis using microarrays was complemented with qRT-PCR analysis of selected genes and a multilevel analysis of the antioxidant glutathione. Partial degradation of DNA and RNA occurred from the onset of artificial ageing at 60% RH and 50 degrees C, and transcriptome profiling showed that the expression of genes associated with programmed cell death, oxidative stress and protein ubiquitination were altered prior to any sign of viability loss. After 25 days of ageing viability started to decline in conjunction with progressively oxidising cellular conditions, as indicated by a shift of the glutathione redox state towards more positive values (>-190 mV). The unravelling of the molecular basis of seed ageing revealed that transcriptome reprogramming is a key component of the ageing process, which influences the progression of programmed cell death and decline in antioxidant capacity that ultimately lead to seed viability loss.

Whole Genome Analysis of Gene Expression Reveals Coordinated Activation of Signaling and Metabolic Pathways during Pollen-Pistil Interactions in Arabidopsis

Plant reproduction depends on the concerted activation of many genes to ensure correct communication between pollen and pistil. Here, we queried the whole transcriptome of Arabidopsis (Arabidopsis thaliana) in order to identify genes with specific reproductive functions. We used the Affymetrix ATH1 whole genome array to profile wild-type unpollinated pistils and unfertilized ovules. By comparing the expression profile of pistils at 0.5, 3.5, and 8.0 h after pollination and applying a number of statistical and bioinformatics criteria, we found 1,373 genes differentially regulated during pollen-pistil interactions. Robust clustering analysis grouped these genes in 16 time-course clusters representing distinct patterns of regulation. Coregulation within each cluster suggests the presence of distinct genetic pathways, which might be under the control of specific transcriptional regulators. A total of 78% of the regulated genes were expressed initially in unpollinated pistil and/or ovules, 15% were initially detected in the pollen data sets as enriched or preferentially expressed, and 7% were induced upon pollination. Among those, we found a particular enrichment for unknown transcripts predicted to encode secreted proteins or representing signaling and cell wall-related proteins, which may function by remodeling the extracellular matrix or as extracellular signaling molecules. A strict regulatory control in various metabolic pathways suggests that fine-tuning of the biochemical and physiological cellular environment is crucial for reproductive success. Our study provides a unique and detailed temporal and spatial gene expression profile of in vivo pollen-pistil interactions, providing a framework to better understand the basis of the molecular mechanisms operating during the reproductive process in higher plants.

Inhibition of glutamine synthetase by phosphinothricin leads to transcriptome reprograming in root nodules of Medicago truncatula

Glutamine synthetase (GS) is a vital enzyme for the assimilation of ammonia into amino acids in higher plants. In legumes, GS plays a crucial role in the assimilation of the ammonium released by nitrogen-fixing bacteria in root nodules, constituting an important metabolic knob controlling the nitrogen (N) assimilatory pathways. To identify new regulators of nodule metabolism, we profiled the transcriptome of Medicago truncatula nodules impaired in N assimilation by specifically inhibiting GS activity using phosphinothricin (PPT). Global transcript expression of nodules collected before and after PPT addition (4, 8, and 24 h) was assessed using Affymetrix M. truncatula GeneChip arrays. Hundreds of genes were regulated at the three time points, illustrating the dramatic alterations in cell metabolism that are imposed on the nodules upon GS inhibition. The data indicate that GS inhibition triggers a fast plant defense response, induces premature nodule senescence, and promotes loss of root nodule identity. Consecutive metabolic changes were identified at the three time points analyzed. The results point to a fast repression of asparagine synthesis and of the glycolytic pathway and to the synthesis of glutamate via reactions alternative to the GS/GOGAT cycle. Several genes potentially involved in the molecular surveillance for internal organic N availability are identified and a number of transporters potentially important for nodule functioning are pinpointed. The data provided by this study contributes to the mapping of regulatory and metabolic networks involved in root nodule functioning and highlight candidate modulators for functional analysis.

Comparative transcriptome profiling of human foreskin fibroblasts infected with the Sylvio and Y strains of Trypanosoma cruzi

Trypanosoma cruzi, the causative agent of Chagas Disease, is phylogenetically distributed into nearly identical genetic strains which show divergent clinical presentations including differences in rates of cardiomyopathy in humans, different vector species and transmission cycles, and differential congenital transmission in a mouse model. The population structure of these strains divides into two groups, which are geographically and clinically distinct. The aim of this study was to compare the transcriptome of two strains of T. cruzi, Sylvio vs. Y to identify differences in expression that could account for clinical and biochemical differences. We collected and sequenced RNA from T. cruzi-infected and control Human Foreskin Fibroblasts at three timepoints. Differential expression analysis identified gene expression profiles at different timepoints in Sylvio infections, and between Sylvio and Y infections in both parasite and host. The Sylvio strain parasite and the host response to Sylvio infection largely mirrored the host-pathogen interaction seen in our previous Y strain work. IL-8 was more highly expressed in Sylvio-infected HFFs than in Y-infected HFFs.

A transcriptome approach towards understanding the development of ripening capacity in 'Bartlett' pears (Pyrus communis L.)

Background: The capacity of European pear fruit (Pyrus communis L.) to ripen after harvest develops during the final stages of growth on the tree. The objective of this study was to characterize changes in 'Bartlett' pear fruit physico-chemical properties and transcription profiles during fruit maturation leading to attainment of ripening capacity. Results: The softening response of pear fruit held for 14days at 20°C after harvest depended on their maturity. We identified four maturity stages: S1-failed to soften and S2- displayed partial softening (with or without ET-ethylene treatment); S3 - able to soften following ET; and S4 - able to soften without ET. Illumina sequencing and Trinity assembly generated 68,010 unigenes (mean length of 911bp), of which 32.8% were annotated to the RefSeq plant database. Higher numbers of differentially expressed transcripts were recorded in the S3-S4 and S1-S2 transitions (2805 and 2505 unigenes, respectively) than in the S2-S3 transition (2037 unigenes). High expression of genes putatively encoding pectin degradation enzymes in the S1-S2 transition suggests pectic oligomers may be involved as early signals triggering the transition to responsiveness to ethylene in pear fruit. Moreover, the co-expression of these genes with Exps (Expansins) suggests their collaboration in modifying cell wall polysaccharide networks that are required for fruit growth. K-means cluster analysis revealed that auxin signaling associated transcripts were enriched in cluster K6 that showed the highest gene expression at S3. AP2/EREBP (APETALA 2/ethylene response element binding protein) and bHLH (basic helix-loop-helix) transcripts were enriched in all three transition S1-S2, S2-S3, and S3-S4. Several members of Aux/IAA (Auxin/indole-3-acetic acid), ARF (Auxin response factors), and WRKY appeared to play an important role in orchestrating the S2-S3 transition. Conclusions: We identified maturity stages associated with the development of ripening capacity in 'Bartlett' pear, and described the transcription profile of fruit at these stages. Our findings suggest that auxin is essential in regulating the transition of pear fruit from being ethylene-unresponsive (S2) to ethylene-responsive (S3), resulting in fruit softening. The transcriptome will be helpful for future studies about specific developmental pathways regulating the transition to ripening. © 2015 Nham et al.

Microarrays as a functional approach to the transcriptome

Knowing a cell’s transcriptome is a fundamental requisite in order to analyze its response to the environment. Microarrays have supposed a revolution on this field as they are able to yield an overview of gene expression at any environmental condition on a genome-wide scale. This technique consists in the hybridisation of a nucleic acid sample, previously marked, with a probe (which might be made up of cDNA, oligonucleotides or PCR products) anchored to a solid surface (made of glass, plastic, silicon...) giving as a result a dot grid which reveals, after image analysis, which genes are being expressed. Nevertheless, this only can be achieved if information on the species genome has been generated. Different kinds of expression microarrays exist attending to the probe’s nature and the method used in its synthesis. In this poster two of these will be treated: Spotted Microarrays, for which the probe is synthesised prior to its fixation to the array and allow the analysis of two targets simultaneously. They can be easily customized, but lack high reproducibility and sensitivity. Oligonucleotide Microarrays, which are characterized by the direct printing of the probe on the array. In this case the probes consist on, invariably, oligonucleotides that are complementary to a small fraction of the gene it is representing at the microarray. Their application is somewhat restricted. This fact, however, makes them more reproducible. Currently, the approach towards the transcriptome studies from the Next Generation Sequencing technologies offers a large volume of information in a short amount of time needing less previous information on the target organism than that needed by microarrays, but their expensive price limits their use. The versatility of the latter, together with their reduced costs in comparison to other techniques, makes them an interesting resource in applications that may need less complexity.

The transcriptome of metamorphosing flatfish

Background Flatfish metamorphosis denotes the extraordinary transformation of a symmetric pelagic larva into an asymmetric benthic juvenile. Metamorphosis in vertebrates is driven by thyroid hormones (THs), but how they orchestrate the cellular, morphological and functional modifications associated with maturation to juvenile/adult states in flatfish is an enigma. Since THs act via thyroid receptors that are ligand activated transcription factors, we hypothesized that the maturation of tissues during metamorphosis should be preceded by significant modifications in the transcriptome. Targeting the unique metamorphosis of flatfish and taking advantage of the large size of Atlantic halibut (Hippoglossus hippoglossus) larvae, we determined the molecular basis of TH action using RNA sequencing. Results De novo assembly of sequences for larval head, skin and gastrointestinal tract (GI-tract) yielded 90,676, 65,530 and 38,426 contigs, respectively. More than 57 % of the assembled sequences were successfully annotated using a multi-step Blast approach. A unique set of biological processes and candidate genes were identified specifically associated with changes in morphology and function of the head, skin and GI-tract. Transcriptome dynamics during metamorphosis were mapped with SOLiD sequencing of whole larvae and revealed greater than 8,000 differentially expressed (DE) genes significantly (p < 0.05) up- or down-regulated in comparison with the juvenile stage. Candidate transcripts quantified by SOLiD and qPCR analysis were significantly (r = 0.843; p < 0.05) correlated. The majority (98 %) of DE genes during metamorphosis were not TH-responsive. TH-responsive transcripts clustered into 6 groups based on their expression pattern during metamorphosis and the majority of the 145 DE TH-responsive genes were down-regulated. Conclusions A transcriptome resource has been generated for metamorphosing Atlantic halibut and over 8,000 DE transcripts per stage were identified. Unique sets of biological processes and candidate genes were associated with changes in the head, skin and GI-tract during metamorphosis. A small proportion of DE transcripts were TH-responsive, suggesting that they trigger gene networks, signalling cascades and transcription factors, leading to the overt changes in tissue occurring during metamorphosis.

Functional relationship of colorectal cancer-associated glycans with malignant phenotype and transcriptome changes

Glycosyltransferases ST6GAL1 and B4GALNT2 (and their cognate antigens Sia6LacNAc and Sda, respectively) are associated with colorectal cancer (CRC) but it is not fully clear their biological and clinical significance. We explored the clinical relevance of both glycosyltransferases by interrogating The Cancer Genome Atlas (TCGA) database while the phenotypic/transcriptomic effects of ST6GAL1/B4GALNT2 overexpression were studied in genetically modified CRC cell lines. Transcriptomic data from CRC patients in TCGA database suggested a moderate impact of ST6GAL1 on CRC progression, although it was not possible to define a clear role for this glycosyltransferase. Transcriptomic analysis of ST6GAL1-transduced cell lines revealed a much deeper effect of ST6GAL1 on gene expression in SW948 than in SW48. The overexpression of ST6GAL1 induced opposite effects on soft agar growth and wound healing in both cell lines. These results indicate that the impact of a cancer-associated glycosyltransferase change on phenotype/transcriptome can be extremely variable, depending on the molecular context of the tumor cell. On the contrary, transcriptomic analysis of B4GALNT2-modified cell lines together with TCGA database survey demonstrated a strong impact of B4GALNT2 on the transcriptional activity of CRC cells, in particular its association with a better prognosis. We suggest an anti-tumoral role of B4GALNT2 in CRC. We also investigated the glycan changes related to ST6GAL1/B4GALNT2 expression in a small cohort of tissues/plasma as well as the N-glycomic profile of CRC, normal and polyp tissues. We found an increase of ST6GAL1 activity in CRC and inflammatory bowel disease plasma samples comparing with plasma from healthy donors. A different Sda protein carrier pattern was observed between healthy donors and CRC plasma samples. β-arrestin 1 is a possible candidate as Sda carrier protein in plasma samples although future validation studies are needed. The alterations found in the N-glycan pattern highlight the importance of N-glycome as a molecular signature in cancer.

Decoding Agc1 deficiency: bioinformatics approaches to unveil the functional implications of Slc25a12 on the transcriptome and epigenome

In the brain, mutations in SLC25A12 gene encoding AGC1 cause an ultra-rare genetic disease reported as a developmental and epileptic encephalopathy associated with global cerebral hypomyelination. Symptoms of the disease include diffused hypomyelination, arrested psychomotor development, severe hypotonia, seizures and are common to other neurological and developmental disorders. Amongst the biological components believed to be most affected by AGC1 deficiency are oligodendrocytes, glial cells responsible for myelination. Recent studies (Poeta et al, 2022) have also shown how altered levels of transcription factors and epigenetic modifications greatly affect proliferation and differentiation in oligodendrocyte precursor cells (OPCs). In this study we explore the transcriptomic landscape of Agc1 in two different system models: OPCs silenced for Agc1 and iPSCs from human patients differentiated to neural progenitors. Analyses range from differential expression analysis, alternative splicing, master regulator analysis. ATAC-seq results on OPCs were integrated with results from RNA-Seq to assess the activity of a TF based on the accessibility data from its putative targets, which allows to integrate RNA-Seq data to infer their role as either activators or repressors. All the findings for this model were also integrated with early data from iPSCs RNA-seq results, looking for possible commonalities between the two different system models, among which we find a downregulation in genes encoding for SREBP, a transcription factor regulating fatty acids biosynthesis, a key process for myelination which could explain the hypomyelinated state of patients. We also find that in both systems cells tend to form more neurites, likely losing their ability to differentiate, considering their progenitor state. We also report several alterations in the chromatin state of cells lacking Agc1, which confirms the hypothesis for which Agc1 is not a disease restricted only to metabolic alterations in the cells, but there is a profound shift of the regulatory state of these cells.