Unknown Transcripts Research Articles

Identification and molecular cloning of Atlantic cod ( Gadus morhua) activating transcription factor 3 (ATF3) transcript and its induction in spleen following intraperitoneal polyriboinosinic polyribocytidylic acid injection

Activating transcription factor 3 (ATF3) participates in cellular processes to adapt to various extra- and intra-cellular changes including the modulation of immunity to prevent uncontrolled immune responses to pathogens. In teleost fishes, the involvement of ATF3 in immune response has not been documented. In this study, the putative Atlantic cod ( Gadus morhua) ATF3 transcript was identified by performing rapid amplification of cDNA ends (RACE) based on unknown expressed sequence tags (ESTs) that are potentially inducible by polyriboinosinic polyribocytidylic acid (pIC, a synthetic double-stranded RNA viral mimic) in Atlantic cod. ATF3-like ESTs were the most abundant unknown transcript (i.e. lacking significant BLAST hits) generated from a previously constructed cDNA library enriched for pIC inducible transcripts in Atlantic cod spleen. The full-length cDNA of cod ATF3 consists of 2329 nucleotides with an open reading frame (ORF) of 735 bp encoding 244 amino acids. The deduced amino acid sequence of Atlantic cod ATF3 shares over 45% identity with its putative orthologs from other vertebrates. In addition, the presence of a conserved basic region leucine zipper (bZIP) domain in the deduced Atlantic cod ATF3-like protein further supports its identity as an ATF3 homolog. In the spleen of Atlantic cod challenged with intraperitoneal (IP) injections of pIC, the time-course transcript expression of ATF3 was studied using quantitative reverse transcription-polymerase chain reaction (QPCR). At 6 h following the pIC injection, the relative expression level of ATF3 mRNA was significantly up-regulated in comparison to a pre-injected control (61.9-fold) and its time-matched saline-injected control (97.3-fold). At 24 h following the pIC injection, the mRNA expression level of cod ATF3 had subsided and was no longer significantly different from its pre-injected control, but significantly higher (1.88-fold) than its time-matched saline-injected control. Collectively, these results suggest that ATF3 may be involved in the modulation of innate anti-viral response in Atlantic cod.

De novo assembly and validation of planaria transcriptome by massive parallel sequencing and shotgun proteomics

Freshwater planaria are a very attractive model system for stem cell biology, tissue homeostasis, and regeneration. The genome of the planarian Schmidtea mediterranea has recently been sequenced and is estimated to contain >20,000 protein-encoding genes. However, the characterization of its transcriptome is far from complete. Furthermore, not a single proteome of the entire phylum has been assayed on a genome-wide level. We devised an efficient sequencing strategy that allowed us to de novo assemble a major fraction of the S. mediterranea transcriptome. We then used independent assays and massive shotgun proteomics to validate the authenticity of transcripts. In total, our de novo assembly yielded 18,619 candidate transcripts with a mean length of 1118 nt after filtering. A total of 17,564 candidate transcripts could be mapped to 15,284 distinct loci on the current genome reference sequence. RACE confirmed complete or almost complete 5' and 3' ends for 22/24 transcripts. The frequencies of frame shifts, fusion, and fission events in the assembled transcripts were computationally estimated to be 4.2%-13%, 0%-3.7%, and 2.6%, respectively. Our shotgun proteomics produced 16,135 distinct peptides that validated 4200 transcripts (FDR ≤1%). The catalog of transcripts assembled in this study, together with the identified peptides, dramatically expands and refines planarian gene annotation, demonstrated by validation of several previously unknown transcripts with stem cell-dependent expression patterns. In addition, our robust transcriptome characterization pipeline could be applied to other organisms without genome assembly. All of our data, including homology annotation, are freely available at SmedGD, the S. mediterranea genome database.

The sequence, structure and evolutionary features of HOTAIR in mammals

BackgroundAn increasing number of long noncoding RNAs (lncRNAs) have been identified recently. Different from all the others that function in cis to regulate local gene expression, the newly identified HOTAIR is located between HoxC11 and HoxC12 in the human genome and regulates HoxD expression in multiple tissues. Like the well-characterised lncRNA Xist, HOTAIR binds to polycomb proteins to methylate histones at multiple HoxD loci, but unlike Xist, many details of its structure and function, as well as the trans regulation, remain unclear. Moreover, HOTAIR is involved in the aberrant regulation of gene expression in cancer.ResultsTo identify conserved domains in HOTAIR and study the phylogenetic distribution of this lncRNA, we searched the genomes of 10 mammalian and 3 non-mammalian vertebrates for matches to its 6 exons and the two conserved domains within the 1800 bp exon6 using Infernal. There was just one high-scoring hit for each mammal, but many low-scoring hits were found in both mammals and non-mammalian vertebrates. These hits and their flanking genes in four placental mammals and platypus were examined to determine whether HOTAIR contained elements shared by other lncRNAs. Several of the hits were within unknown transcripts or ncRNAs, many were within introns of, or antisense to, protein-coding genes, and conservation of the flanking genes was observed only between human and chimpanzee. Phylogenetic analysis revealed discrete evolutionary dynamics for orthologous sequences of HOTAIR exons. Exon1 at the 5' end and a domain in exon6 near the 3' end, which contain domains that bind to multiple proteins, have evolved faster in primates than in other mammals. Structures were predicted for exon1, two domains of exon6 and the full HOTAIR sequence. The sequence and structure of two fragments, in exon1 and the domain B of exon6 respectively, were identified to robustly occur in predicted structures of exon1, domain B of exon6 and the full HOTAIR in mammals.ConclusionsHOTAIR exists in mammals, has poorly conserved sequences and considerably conserved structures, and has evolved faster than nearby HoxC genes. Exons of HOTAIR show distinct evolutionary features, and a 239 bp domain in the 1804 bp exon6 is especially conserved. These features, together with the absence of some exons and sequences in mouse, rat and kangaroo, suggest ab initio generation of HOTAIR in marsupials. Structure prediction identifies two fragments in the 5' end exon1 and the 3' end domain B of exon6, with sequence and structure invariably occurring in various predicted structures of exon1, the domain B of exon6 and the full HOTAIR.

Characterization of retrotransposon sequences expressed in inflorescences of apomictic and sexual Paspalum notatum plants

Apomixis, an asexual mode of reproduction through seeds, holds much promise for agricultural advances. However, the molecular mechanisms underlying this trait are still poorly understood. We previously isolated several transcripts representing novel sequences differentially expressed in reproductive tissues of sexual and apomictic plants. Here, we report the characterization of two of these unknown RNA transcripts (experimental codes N17 and N22). Since original fragments showed no significant homologies to sequences at databases, preliminary genomic PCR experiments were carried out to discard possible contaminations. RACE extension on flanking regions provided longer sequences for the candidates and additional related transcripts, which revealed similarity to LTR retrotransposons carrying short transduplicated segments of protein-coding genes. Interestingly, some transduplicated segments corresponded to genes previously associated with apomictic development. Gene copy number estimations revealed a moderate representation of the elements in the genome, with significantly increased numbers in a sexual genotype with respect to an apomictic one. Genetic mapping of N17 showed that a copy of this particular element was located onto Paspalum notatum linkage group F3c, at a central non-recombinant region resembling a centromere. Expression analysis showed an increased activity of N17 and N22 sense strands in ovules of the sexual genotypes. A retrotransposon-specific differential display analysis aimed at detecting related sequences allowed the identification of a complex family, with the majority of its members represented in the sexual genotype. Our results suggest that these elements could be participating in regulatory pathways related to apomixis and sexuality.

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.

Endogenous retrovirus drives hitherto unknown proapoptotic p63 isoforms in the male germ line of humans and great apes

TAp63, but not its homolog p53, eliminates oocytes that suffered DNA damage. An equivalent gene for guarding the male germ line is currently not known. Here we identify hitherto unknown human p63 transcripts with unique 5'-ends derived from incorporated exons upstream of the currently mapped TP63 gene. These unique p63 transcripts are highly and specifically expressed in testis. Their most upstream region corresponds to a LTR of the human endogenous retrovirus 9 (ERV9). The insertion of this LTR upstream of the TP63 locus occurred only recently in evolution and is unique to humans and great apes (Hominidae). A corresponding p63 protein is the sole p63 species in healthy human testis, and is strongly expressed in spermatogenic precursors but not in mature spermatozoa. In response to DNA damage, this human male germ-cell-encoded TAp63 protein (designated GTAp63) is activated by caspase cleavage near its carboxyterminal domain and induces apoptosis. Human testicular cancer tissues and cell lines largely lost p63 expression. However, pharmacological inhibition of histone deacetylases completely restores p63 expression in testicular cancer cells (>3,000-fold increase). Our data support a model whereby testis-specific GTAp63 protects the genomic integrity of the male germ line and acts as a tumor suppressor. In Hominidae, this guardian function was greatly enhanced by integration of an endogenous retrovirus upstream of the TP63 locus that occurred 15 million years ago. By providing increased germ-line stability, this event may have contributed to the evolution of hominids and enabled their long reproductive periods.

Complexity of transcriptional regulation within the Rag locus: identification of a second Nwc promoter region within the Rag2 intron

Nwc represents a mysterious third evolutionarily conserved gene within the Rag locus. Here, we analyzed the phenotype of Nwctmpro1 mice, in which the Rag2 intragenic region containing the previously identified promoter responsible for initiating transcription of Nwc in all cells except lymphocytes was deleted by homologous recombination. Despite strong nonlymphocyte-specific inhibition of Nwc transcription which runs through the regulatory region of Rag genes, their expression remained suppressed, and no developmental, morphological, anatomical, functional, physiological, or cellular defects in Nwctmpro1 mice could be observed. However, careful analysis of the Rag2 intergenic region uncovered a second evolutionarily conserved Nwc promoter region from which a previously unknown Nwc transcript can be generated in nonlymphocytes of Nwctmpro1 and normal mice. The above results reveal an unexpected additional complexity of transcriptional regulation within the Rag/Nwc locus and show that strong inhibition of Nwc transcription in nonlymphoid cells is well tolerated. Complete inactivation of Nwc is necessary to get insight into its function at transcriptional and posttranscriptional levels.

Transcript profiling during the early development of the maize brace root via Solexa sequencing

Currently, the molecular regulation mechanisms involved in the early development of maize brace root are poorly known. To gain insight into the transcriptome dynamics that are associated with its development, genome-wide gene expression profiling was conducted by Solexa sequencing (Illumina Inc., San Diego, CA, USA). More than five million tags were generated from the stem node tissues without and with just-emerged brace roots, including 149,524 and 178,131 clean tags in the two libraries, respectively. Of these, 82,864 (55.4%) and 91,678 (51.5%) tags were matched to the reference genes. The most differentially expressed tags with a log(2) ratio > 2 or < -2 (P < 0.001) were analyzed further, representing 143 up-regulated and 152 down-regulated genes, except for unknown transcripts, which were classified into 11 functional categories. The most enriched categories were those of metabolism, signal transduction and cellular transport. Many genes or biological pathways were found to be commonly shared between brace root and lateral or adventitious root development, such as genes participating in cell wall degradation and synthesis, auxin transport and signaling, ethylene signaling, etc. Next, the expression patterns of 20 genes were assessed by quantitative real-time PCR, and the results obtained showed general agreement with the Solexa analysis. Furthermore, a comparison of the brace root transcriptome with that of maize primary root revealed substantial differences in the categories and abundances of expressed transcripts. In conclusion, we first reveal the complex changes in the transcriptome during the early development of maize brace root and provide a comprehensive set of data that are essential for understanding its molecular regulation.

BCOR as a Novel Fusion Partner of Retinoic Acid Receptor Alpha In a t(X;17)(p11;q12) Variant of Acute Promyelocytic Leukemia

AbstractAbstract 1703The majority of acute promyelocytic leukemia (APL) cases are characterized by the presence of a PML-RARA fusion gene. In a small subset, RARA is fused to a different partner including PLZF, NPM1, NuMA, STAT5b, PRKAR1A and FIP1L1. Here we identified a novel RARA fusion transcript, BCOR-RARA in a t(X;17)(p11;q12) variant of APL. The patient was a 45-year-old man. Although the patient was clinically responsive to ATRA, repeated standard chemotherapy with ATRA did not effect a cure. The bone marrow promyelocytes had unique morphologic features, including rectangular and round cytoplasmic inclusion bodies. They were more granular than those of AML M2 but less granular than the classical t(15;17) APL. Flow cytometric analysis revealed strong expression of CD13, CD33, CD56, weak expression of CD11c and lack of HLA-DR and CD7. The karyotype analysis detected a novel chromosomal translocation described as 45,-Y, t(X;17)(p11.2;q12)[19]/ 46,XY[1]. FISH analysis indicated one intact and two split signals of RARA and two intact signals of PML. To amplify unknown chimeric fusion transcripts, we performed 5'-RACE. The sequence revealed that BCL6 co-repressor, BCOR cDNA from exons 9 to 12 to be fused to RARA exon 3. By RT-PCR, we confirmed full length chimeric fusion transcripts spanning from the start codon to 4,948 nt of BCOR cDNA (NM_001123384) fused to RARA cDNA from exon 3 to the stop codon. The chimeric cDNA had an in-frame codon from BCOR through RARA, creating a 1,931 amino acid fusion protein. One of the consistent features in all known RARA fusion partners is self-association. To determine whether BCOR-RARA self-associate, we performed co-immunoprecipitation assays. These results showed that BCOR-RARA is able to self-associate both through the region of BCOR-S and the ankyrin repeat domain of BCOR. In addition, BCOR-RARA associated with BCL6. RXR recruitment is a critical determinant of transforming potential of oligomeric RARA fusion proteins. To investigate how BCOR-RARA associates with RARE in vitro, we performed EMSA. These results showed that BCOR-RARA/RXRA complex associates with RARE in an alternative manner compared to RARA and PML-RARA. Deregulation of RARA transcriptional activations has a central role in pathogenesis of APL. Therefore, we evaluated ATRA-induced transcriptional activation of 4× RAREs with a reporter assay in HepG2 cells. Without ATRA, BCOR-RARA repressed the reporter activity. With addition of ATRA, BCOR-RARA induced transcriptional activation very weakly. Subsequently, we evaluated dominant-negative effects of the samples in the RARA/RXRA pathway. In contrast to BCOR, BCOR-RARA clearly inhibited ATRA-induced RARA transcriptional activation as well as PML-RARA. Furthermore, we asked which domains are sufficient for the dominant-negative effects with the deletion mutants. The results indicated that the region spanning from 999 to 1,409 aa of BCOR-RARA has pivotal roles in the dominant-negative effects. Correct protein function is highly dependent on intracellular localization. To investigate subcellular localization of BCOR-RARA, we performed immunofluorescence analysis in 293T cells. In BCOR-RARA-expressing cells, BCOR-RARA localized as two patterns; (I) diffusely in the nucleus as well as PML-RARA in 82% of the cells, (II) diffusely in the nucleus and aggregately in the cytoplasm in 18% of the cells. The subcellular localization of BCOR-RARA was clearly distinguishable from the punctuate pattern as shown in the nucleus of BCOR-expressing cells. Moreover, co-immunofluorescence analysis between BCOR-RARA and BCL6 indicated that the subcellular localization of BCOR-RARA/BCL6 is distinct from BCOR/BCL6. BCOR-RARA was found to possess common features with other RARA fusion proteins. These included: (I) the same break point in RARA cDNA; (II) self-association; (III) RXRA is necessary for BCOR-RARA to associate with the RARA responsive element; (IV) action in a dominant-negative manner on RARA transcriptional activation; (V) aberrant subcellular relocalization. It should be noted that there was no intact BCOR found in the 45,-Y,t(X;17)(p11;q12) APL cells because they featured only a rearranged × chromosome. These results highlight essential features of pathogenesis in APL in more detail. BCOR appears to be involved not only in human congenital diseases but also in a human cancer. Disclosures:No relevant conflicts of interest to declare.

The Mouse Blood-Brain Barrier Transcriptome: A New Resource for Understanding the Development and Function of Brain Endothelial Cells

The blood-brain barrier (BBB) maintains brain homeostasis and limits the entry of toxins and pathogens into the brain. Despite its importance, little is known about the molecular mechanisms regulating the development and function of this crucial barrier. In this study we have developed methods to highly purify and gene profile endothelial cells from different tissues, and by comparing the transcriptional profile of brain endothelial cells with those purified from the liver and lung, we have generated a comprehensive resource of transcripts that are enriched in the BBB forming endothelial cells of the brain. Through this comparison we have identified novel tight junction proteins, transporters, metabolic enzymes, signaling components, and unknown transcripts whose expression is enriched in central nervous system (CNS) endothelial cells. This analysis has identified that RXRalpha signaling cascade is specifically enriched at the BBB, implicating this pathway in regulating this vital barrier. This dataset provides a resource for understanding CNS endothelial cells and their interaction with neural and hematogenous cells.

Identification and characterization of a novel miR159 target not related to MYB in tomato

MicroRNA 159 (miR159) is a highly conserved miRNA with roles in flowering under short days, anther development and seed germination via repression of GAMYB-like genes. In tomato, the function of miR159 (Sl-miR159) is currently unknown and target transcripts have not been experimentally validated. Here, we identified and characterized a new miR159 target gene (SGN-U567133) in Solanum lycopersicum (tomato) that is not related to MYB. SGN-U567133 is predominantly expressed in flowers and encodes a nuclear-localized protein that contains a unique NOZZLE-like domain at its N terminus. In tomato, SGN-U567133 represents a small gene family and orthologs have been identified in other plant species, all containing a conserved miR159 target site in their coding sequence. Accordingly, 5'-RACE cleavage assay supported miRNA-mediated cleavage of SGN-U567133 transcripts in vivo. Moreover, the SGN-U567133 transcript accumulated in P19-HA-expressing tomato leaves in which miRNA-mediated cleavage is inhibited. In addition, transgenic tomato plants expressing a miR159-resistant form of SGN-U567133 accumulated higher levels of the SGN-U567133 transcript and exhibited defects in leaf and flower development. Together, our results suggest that SGN-U567133 represents a novel class of miR159 targets in plants and raise the possibility that its post-transcriptional regulation by Sl-miR159 is essential for normal tomato development.

Identification and analysis of immune-related transcriptome in Asian seabass Lates calcarifer

BackgroundFish diseases caused by pathogens are limiting their production and trade, affecting the economy generated by aquaculture. Innate immunity system is the first line of host defense in opposing pathogenic organisms or any other foreign material. For identification of immune-related genes in Asian seabass Lates calcarifer, an important marine foodfish species, we injected bacterial lipopolysaccharide (LPS), a commonly used elicitor of innate immune responses to eight individuals at the age of 35 days post-hatch and applied the suppression subtractive hybridization (SSH) technique to selectively amplify spleen cDNA of differentially expressed genes.ResultsSequencing and bioinformatic analysis of 3351 ESTs from two SSH libraries yielded 1692 unique transcripts. Of which, 618 transcripts were unknown/novel genes and the remaining 1074 were similar to 743 known genes and 105 unannotated mRNA sequences available in public databases. A total of 161 transcripts were classified to the category "response to stimulus" and 115 to "immune system process". We identified 25 significantly up-regulated genes (including 2 unknown transcripts) and 4 down-regulated genes associated with immune-related processes upon challenge with LPS. Quantitative real-time PCR confirmed the differential expression of these genes after LPS challenge.ConclusionsThe present study identified 1692 unique transcripts upon LPS challenge for the first time in Asian seabass by using SSH, sequencing and bioinformatic analysis. Some of the identified transcripts are vertebrate homologues and others are hitherto unreported putative defence proteins. The obtained immune-related genes may allow for a better understanding of immunity in Asian seabass, carrying out detailed functional analysis of these genes and developing strategies for efficient immune protection against infections in Asian seabass.

Multiple bovine FcγRIIb sub-isoforms generated by alternative splicing

Receptors for the Fc portion of immunoglobulin molecules (FcR) provide an important and vital link between circulating antibody and cellular effector functions. These receptors have been well characterized in human and murine species, however few of these receptors have been investigated in livestock. FcγRII (CD32) is an FcR previously shown in mice and humans to exist in multiple isoforms, both activating (FcγRIIa, FcγRIIc) and inhibitory (FcγRIIb), on a wide variety of cells including B cells, T cells, dendritic cells, monocytes, macrophages and platelets. On B cells, FcγRIIb acts to suppress cell activation and immunoglobulin production by means of an intracellular immunoreceptor tyrosine-based inhibitory motif signaling domain. Two sub-isoforms of FcγRIIb, designated b1 and b2, distinguished by the inclusion of an additional cytoplasmic exon in the b1 form, have been demonstrated in humans and mice, whereas only one sequence corresponding to the human and mouse b2 isoform has been identified in cattle. In this study, the expression profile of FcγRIIb in bovine blood mononuclear cells was characterized by collecting blood samples from mature cattle of dairy and beef breeds, and determining their FcγRIIb mRNA expression profile by RT-PCR. Analysis revealed the presence of two uncharacterized bovine FcγRIIb transcripts in addition to the single previously published transcript. Analysis of the first unknown transcript revealed high homology with published human and murine FcγRIIb1 sequences. This transcript was present in all cell types examined, with little variation in primary sequence between individuals or among breeds. The second unknown sequence was found to be homologous to the murine FcγRIIb3 (IgG-binding protein or soluble FcγR in humans) sequence. This transcript appears to have a much more limited expression profile, which may indicate that expression varies with the cellular activation-state of the cell. These results indicate that cattle, like humans and mice, express multiple sub-isoforms of FcγRIIb. These findings add further complexity to the regulation of IgG-mediated immunity and provide new insight into the role Fc receptors play in antigen acquisition and presentation in cattle.

Copy number variation at 1q21.1 associated with neuroblastoma

Common copy number variations (CNVs) represent a significant source of genetic diversity, yet their influence on phenotypic variability, including disease susceptibility, remains poorly understood. To address this problem in human cancer, we performed a genome-wide association study of CNVs in the childhood cancer neuroblastoma, a disease in which single nucleotide polymorphism variations are known to influence susceptibility. We first genotyped 846 Caucasian neuroblastoma patients and 803 healthy Caucasian controls at approximately 550,000 single nucleotide polymorphisms, and performed a CNV-based test for association. We then replicated significant observations in two independent sample sets comprised of a total of 595 cases and 3,357 controls. Here we describe the identification of a common CNV at chromosome 1q21.1 associated with neuroblastoma in the discovery set, which was confirmed in both replication sets. This CNV was validated by quantitative polymerase chain reaction, fluorescent in situ hybridization and analysis of matched tumour specimens, and was shown to be heritable in an independent set of 713 cancer-free parent-offspring trios. We identified a previously unknown transcript within the CNV that showed high sequence similarity to several neuroblastoma breakpoint family (NBPF) genes and represents a new member of this gene family (NBPF23). This transcript was preferentially expressed in fetal brain and fetal sympathetic nervous tissues, and the expression level was strictly correlated with CNV state in neuroblastoma cells. These data demonstrate that inherited copy number variation at 1q21.1 is associated with neuroblastoma and implicate a previously unknown neuroblastoma breakpoint family gene in early tumorigenesis of this childhood cancer.

A transcriptional sketch of a primary human breast cancer by 454 deep sequencing

BackgroundThe cancer transcriptome is difficult to explore due to the heterogeneity of quantitative and qualitative changes in gene expression linked to the disease status. An increasing number of "unconventional" transcripts, such as novel isoforms, non-coding RNAs, somatic gene fusions and deletions have been associated with the tumoral state. Massively parallel sequencing techniques provide a framework for exploring the transcriptional complexity inherent to cancer with a limited laboratory and financial effort. We developed a deep sequencing and bioinformatics analysis protocol to investigate the molecular composition of a breast cancer poly(A)+ transcriptome. This method utilizes a cDNA library normalization step to diminish the representation of highly expressed transcripts and biology-oriented bioinformatic analyses to facilitate detection of rare and novel transcripts.ResultsWe analyzed over 132,000 Roche 454 high-confidence deep sequencing reads from a primary human lobular breast cancer tissue specimen, and detected a range of unusual transcriptional events that were subsequently validated by RT-PCR in additional eight primary human breast cancer samples. We identified and validated one deletion, two novel ncRNAs (one intergenic and one intragenic), ten previously unknown or rare transcript isoforms and a novel gene fusion specific to a single primary tissue sample. We also explored the non-protein-coding portion of the breast cancer transcriptome, identifying thousands of novel non-coding transcripts and more than three hundred reads corresponding to the non-coding RNA MALAT1, which is highly expressed in many human carcinomas.ConclusionOur results demonstrate that combining 454 deep sequencing with a normalization step and careful bioinformatic analysis facilitates the discovery and quantification of rare transcripts or ncRNAs, and can be used as a qualitative tool to characterize transcriptome complexity, revealing many hitherto unknown transcripts, splice isoforms, gene fusion events and ncRNAs, even at a relatively low sequence sampling.

Insights in to the pathogenesis of axial spondyloarthropathy based on gene expression profiles

IntroductionAxial spondyloarthropathy (SpA) is a group of inflammatory diseases, with ankylosing spondylitis as the prototype. SpA affects the axial skeleton, entheses, joints and, at times, the eyes. This study tested the hypothesis that SpA is characterized by a distinct pattern of gene expression in peripheral blood of affected individuals compared with healthy controls.MethodsHigh-density, human GeneChip® probe arrays were used to profile mRNA of peripheral blood cells from 18 subjects with SpA and 25 normal individuals. Samples were processed as two separate sets at different times (11 SpA + 12 control subjects in primary set (Set 1); 7 SpA+ 13 control subjects in the validation set (Set 2)). Blood samples were taken at a time when patients were not receiving systemic immunomodulatory therapy. Differential expression was defined as a 1.5-fold change with a q value < 5%. Gene ontology and pathway information were also studied.ResultsSignals from 134 probe sets (representing 95 known and 12 unknown gene transcripts) were consistently different from controls in both Sets 1 and 2. Included among these were transcripts for a group of 20 genes, such as interleukin-1 (IL-1) receptors 1 and 2, Nod-like receptor family, pyrin domain containing 2 (NLRP2), secretory leukocyte peptidase inhibitor (SLPI), secreted protein acidic and rich in cysteine (SPARC), and triggering receptor expressed on myeloid cells 1 (TREM-1) that are clearly related to the immune or inflammatory response and a group of 4 transcripts that have a strong role in bone remodeling.ConclusionsOur observations are the first to implicate SPARC, SLPI, and NLRP2, a component of the innate immune system, in the pathogenesis of SpA. Our results also indicate a possible role for IL-1 and its receptors in SpA. In accord with the bone pathology component of SpA, we also found that expression levels of transcripts reflecting bone remodeling factors are also distinguishable in peripheral blood from patients with SpA versus controls. These results confirm some previously identified biomarkers implicated in the pathogenesis of SpA and also point to novel mediators in this disease.

Serial Analysis of Gene Expression (SAGE): 13 Years of Application in Research

A number of molecular methods of gene expression analysis can approach genomic level. Among those, Serial Analysis of Gene Expression (SAGE) stands out. Unlike many other techniques, SAGE allows both qualitative and quantitative analysis of previously unknown transcripts. Over the course of the last 13 years, SAGE has became a recognized tool of large-scale gene expression profiling, being used extensively in human, animal, yeast and plant studies of various nature. A number of important adaptations was introduced both to the protocol of SAGE library construction and to the analytical algorithm employed. Moreover, some variations of the original protocol (MAGE, SADE, microSAGE, miniSAGE, longSAGE, superSAGE, deepSAGE, etc.) were derived to improve the utility of SAGE in certain conditions. Current review aims comparing the benefits and drawbacks of the techniques for high-throughput gene expression analysis (including SAGE) in a realistic, balanced manner. Issues related to modifications to the original protocol and further development of the SAGE are discussed.

Detection of Large Numbers of Novel Sequences in the Metatranscriptomes of Complex Marine Microbial Communities

BackgroundSequencing the expressed genetic information of an ecosystem (metatranscriptome) can provide information about the response of organisms to varying environmental conditions. Until recently, metatranscriptomics has been limited to microarray technology and random cloning methodologies. The application of high-throughput sequencing technology is now enabling access to both known and previously unknown transcripts in natural communities.Methodology/Principal FindingsWe present a study of a complex marine metatranscriptome obtained from random whole-community mRNA using the GS-FLX Pyrosequencing technology. Eight samples, four DNA and four mRNA, were processed from two time points in a controlled coastal ocean mesocosm study (Bergen, Norway) involving an induced phytoplankton bloom producing a total of 323,161,989 base pairs. Our study confirms the finding of the first published metatranscriptomic studies of marine and soil environments that metatranscriptomics targets highly expressed sequences which are frequently novel. Our alternative methodology increases the range of experimental options available for conducting such studies and is characterized by an exceptional enrichment of mRNA (99.92%) versus ribosomal RNA. Analysis of corresponding metagenomes confirms much higher levels of assembly in the metatranscriptomic samples and a far higher yield of large gene families with >100 members, ∼91% of which were novel.Conclusions/SignificanceThis study provides further evidence that metatranscriptomic studies of natural microbial communities are not only feasible, but when paired with metagenomic data sets, offer an unprecedented opportunity to explore both structure and function of microbial communities – if we can overcome the challenges of elucidating the functions of so many never-seen-before gene families.

Isolation and Analysis of Cold Stress Inducible Genes in Zea mays by Suppression Subtractive Hybridization and cDNA Macroarray

In order to understand the molecular and cellular mechanisms underlying cold stress conditions (4°C) in maize seedlings, a forward subtractive cDNA library was constructed using the suppression subtractive hybridization (SSH) technique. Through the “Virtual” Northern blot analysis, 893 positive clones were screened from a total 1,200 clones in the subtractive cDNA library. After sequencing 528 randomly chosen cDNA clones, 213 uniquely expressed sequence tags (ESTs) were obtained by clustering and blast analysis, which included transcripts that had previously been reported as responsive to stress as well as some functionally unknown transcripts. Based on a list of functional Arabidopsis protein categories, the ESTs with significant protein similarity were sorted into ten functional categories. A cDNA macroarray containing the 213 unique ESTs was used to monitor the spatial and temporal distribution of gene expression in maize seedlings during cold stress. The results showed that 118 ESTs were induced by cold-stress in maize seedlings and 66 ESTs identified in the leaves and 89 ESTs in the roots. Hierarchical cluster analysis indicated that the expression profiles of cold stress inducible ESTs in the leaves were different from that observed in the roots. Moreover, some induced genes were related to sugar synthesis and reestablishment of high rates of photosynthesis. In addition, Northern blot analysis validated well the cDNA macroarray data.

MicroRNAs expressed by herpes simplex virus 1 during latent infection regulate viral mRNAs

Herpesviruses are characterized by their ability to maintain life-long latent infections in their animal hosts. However, the mechanisms that allow establishment and maintenance of the latent state remain poorly understood. Herpes simplex virus 1 (HSV-1) establishes latency in neurons of sensory ganglia, where the only abundant viral gene product is a non-coding RNA, the latency associated transcript (LAT). Here we show that LAT functions as a primary microRNA (miRNA) precursor that encodes four distinct miRNAs in HSV-1 infected cells. One of these miRNAs, miR-H2-3p, is transcribed in an antisense orientation to ICP0-a viral immediate-early transcriptional activator that is important for productive HSV-1 replication and thought to have a role in reactivation from latency. We show that miR-H2-3p is able to reduce ICP0 protein expression, but does not significantly affect ICP0 messenger RNA levels. We also identified a fifth HSV-1 miRNA in latently infected trigeminal ganglia, miR-H6, which derives from a previously unknown transcript distinct from LAT. miR-H6 shows extended seed complementarity to the mRNA encoding a second HSV-1 transcription factor, ICP4, and inhibits expression of ICP4, which is required for expression of most HSV-1 genes during productive infection. These results may explain the reported ability of LAT to promote latency. Thus, HSV-1 expresses at least two primary miRNA precursors in latently infected neurons that may facilitate the establishment and maintenance of viral latency by post-transcriptionally regulating viral gene expression.