Construction of long-transcript enriched cDNA libraries from submicrogram amounts of total RNAs by a universal PCR amplification method

De novo emergence, existence, and demise of a protein-coding gene in murids

BACKGROUND

Genes, principal units of genetic information, vary in complexity and evolutionary history. Less-complex genes (e.g., long non-coding RNA (lncRNA) expressing genes) readily emerge de novo from non-genic sequences and have high evolutionary turnover. Genesis of a gene may be facilitated by adoption of functional genic sequences from retrotransposon insertions. However, protein-coding sequences in extant genomes rarely lack any connection to an ancestral protein-coding sequence.

RESULTS

We describe remarkable evolution of the murine gene D6Ertd527e and its orthologs in the rodent Muroidea superfamily. The D6Ertd527e emerged in a common ancestor of mice and hamsters most likely as a lncRNA-expressing gene. A major contributing factor was a long terminal repeat (LTR) retrotransposon insertion carrying an oocyte-specific promoter and a 5' terminal exon of the gene. The gene survived as an oocyte-specific lncRNA in several extant rodents while in some others the gene or its expression were lost. In the ancestral lineage of Mus musculus, the gene acquired protein-coding capacity where the bulk of the coding sequence formed through CAG (AGC) trinucleotide repeat expansion and duplications. These events generated a cytoplasmic serine-rich maternal protein. Knock-out of D6Ertd527e in mice has a small but detectable effect on fertility and the maternal transcriptome.

CONCLUSIONS

While this evolving gene is not showing a clear function in laboratory mice, its documented evolutionary history in Muroidea during the last ~ 40 million years provides a textbook example of how a several common mutation events can support de novo gene formation, evolution of protein-coding capacity, as well as gene's demise.

Improved Nanopore full-length cDNA sequencing by PCR-suppression

Full-length transcript sequencing remains a main goal of RNA sequencing. However, even the application of long-read sequencing technologies such as Oxford Nanopore Technologies still fail to yield full-length transcript sequencing for a significant portion of sequenced reads. Since these technologies can sequence reads that are far longer than the longest known processed transcripts, the lack of efficiency to obtain full-length transcripts from good quality RNAs stems from library preparation inefficiency rather than the presence of degraded RNA molecules. It has previously been shown that addition of inverted terminal repeats in cDNA during reverse transcription followed by single-primer PCR creates a PCR suppression effect that prevents amplification of short molecules thus enriching the library for longer transcripts. We adapted this method for Nanopore cDNA library preparation and show that not only is PCR efficiency increased but gene body coverage is dramatically improved. The results show that implementation of this simple strategy will result in better quality full-length RNA sequencing data and make full-length transcript sequencing possible for most of sequenced reads.

Long terminal repeats power evolution of genes and gene expression programs in mammalian oocytes and zygotes

Retrotransposons are "copy-and-paste" insertional mutagens that substantially contribute to mammalian genome content. Retrotransposons often carry long terminal repeats (LTRs) for retrovirus-like reverse transcription and integration into the genome. We report an extraordinary impact of a group of LTRs from the mammalian endogenous retrovirus-related ERVL retrotransposon class on gene expression in the germline and beyond. In mouse, we identified more than 800 LTRs from ORR1, MT, MT2, and MLT families, which resemble mobile gene-remodeling platforms that supply promoters and first exons. The LTR-mediated gene remodeling also extends to hamster, human, and bovine oocytes. The LTRs function in a stage-specific manner during the oocyte-to-embryo transition by activating transcription, altering protein-coding sequences, producing noncoding RNAs, and even supporting evolution of new protein-coding genes. These functions result, for example, in recycling processed pseudogenes into mRNAs or lncRNAs with regulatory roles. The functional potential of the studied LTRs is even higher, because we show that dormant LTR promoter activity can rescue loss of an essential upstream promoter. We also report a novel protein-coding gene evolution-D6Ertd527e-in which an MT LTR provided a promoter and the 5' exon with a functional start codon while the bulk of the protein-coding sequence evolved through a CAG repeat expansion. Altogether, ERVL LTRs provide molecular mechanisms for stochastically scanning, rewiring, and recycling genetic information on an extraordinary scale. ERVL LTRs thus offer means for a comprehensive survey of the genome's expression potential, tightly intertwining with gene expression and evolution in the germline.

Inactivation of a single copy of Crebbp selectively alters pre-mRNA processing in mouse hematopoietic stem cells

Global expression analysis of fetal liver hematopoietic stem cells (FL HSCs) revealed the presence of unspliced pre-mRNA for a number of genes in normal FL HSCs. In a subset of these genes, Crebbp+/− FL HSCs had less unprocessed pre-mRNA without a corresponding reduction in total mRNA levels. Among the genes thus identified were the key regulators of HSC function Itga4, Msi2 and Tcf4. A similar but much weaker effect was apparent in Ep300+/− FL HSCs, indicating that, in this context as in others, the two paralogs are not interchangeable. As a group, the down-regulated intronic probe sets could discriminate adult HSCs from more mature cell types, suggesting that the underlying mechanism is regulated with differentiation stage and is active in both fetal and adult hematopoiesis. Consistent with increased myelopoiesis in Crebbp hemizygous mice, targeted reduction of CREBBP abundance by shRNA in the multipotent EML cell line triggered spontaneous myeloid differentiation in the absence of the normally required inductive signals. In addition, differences in protein levels between phenotypically distinct EML subpopulations were better predicted by taking into account not only the total mRNA signal but also the amount of unspliced message present. CREBBP thus appears to selectively influence the timing and degree of pre-mRNA processing of genes essential for HSC regulation and thereby has the potential to alter subsequent cell fate decisions in HSCs.

The identification of three novel genes involved in the rapid-growth regulation in a marine diatom, Skeletonema costatum

In order to understand the molecular mechanism of growth regulation in phytoplankton and to develop novel growth-status indicators, a subtraction cDNA library was constructed by using the mRNA extracted from Skeletonema costatum in the rapid-growth stage (RG stage), and three RG-stage-related cDNA fragments, RG#14, RG#25, and RG#42, were obtained. According to the results of sequence analysis, RG#42 belonged to the MCM2-7 protein family, and the other two fragments, RG#14 and RG#25, were novel molecules. Under continuous illumination, these RG-stage-related mRNA expression levels increased from 100- (RG#14 and RG#42) to 1,000-fold (RG#25) with increasing growth rate. Furthermore, under a diel rhythm of light (light-dark = 12:12 h), the daily mean mRNA abundances of RG#14 and RG#25 in the exponential phase also differed from those in the late-stationary phase. However, such differences between these growth phases were not observed in the mRNA levels of RG#42 and PCNA. This study not only provided a new way to investigate the regulatory mechanisms of cell growth but also offered a possibility of employing these gene fragments as indicators to monitor the growth status of phytoplankton in the marine environment.

Induction of zinc-finger proliferation 1 expression in non-myelinating Schwann cells after denervation

Terminal Schwann cells (tSCs) are non-myelinating glia that wrap the nerve terminal at the neuromuscular junction. They are required for the maintenance of the neuromuscular synapse and are likely to play essential roles in the restoration of synaptic connections after nerve injury. tSCs acquire a reactive phenotype after nerve damage characterized by the extension of cellular processes that may facilitate reinnervation. The molecular signaling events underpinning the tSC reactive state remain elusive, in particular, little is known about transcription factors involved in the transcriptional reprogramming during tSC activation. Prior research implicated nine members of the zinc-finger transcription factor family in Schwann cell (SC) development and myelination, and levels of one such protein were reported increased in other non-myelinating SCs after denervation. We hypothesize that zinc-finger transcription factors could play a role during tSC activation. Because of their relative paucity, tSCs are difficult to study molecularly. Here, we used the rat cervical sympathetic trunk (CST), an autonomic nerve in which non-myelinating SCs are the predominant cell type, to isolate zinc-finger protein (ZFP) cDNAs by reverse transcriptase-polymerase chain reaction. We isolated 29 unique ZFP sequences of which zinc proliferation 1 (Zipro1) was the most abundant. We found that after CST transection, levels for Zipro1 mRNA doubled and that Zipro1 protein expression increased in non-myelinating CST SCs. We also determined that Zipro1 is expressed in tSCs and its levels increased following skeletal muscle denervation. Thus, Zipro1 is a good candidate for a transcription factor involved in activation of non-myelinating SCs in general, and tSCs in particular.

Cooccurrence of ScDSP gene expression, cell death, and DNA fragmentation in a marine diatom, Skeletonema costatum

A novel death-specific gene, ScDSP, was obtained from a death stage subtraction cDNA library of the diatom Skeletonema costatum. The full length of ScDSP cDNA was 921 bp in length, containing a 699-bp open reading frame encoding 232 amino acids and two stretches of 66 and 156 bp in the 5' and 3' untranslated regions, respectively. Analysis of the peptide structure revealed that ScDSP contained a signal peptide domain, a transmembrane domain, and a pair of EF-hand motifs. When S. costatum grew exponentially at a rate of 1.3 day(-1), the ScDSP mRNA level was at 2 mumol . mole 18S rRNA(-1). In contrast, when the culture entered the death phase with a growth rate decreasing to 0.5 day(-1), ScDSP mRNA increased dramatically to 668 mumol . mole 18S rRNA(-1), and a high degree of DNA fragmentation was simultaneously observed. Under the influence of a light-dark cycle, ScDSP expression in both exponential and stationary phases clearly showed a diel rhythm, but the daily mean mRNA level was significantly higher in the stationary phase. Our results suggest that ScDSP may play a role in the molecular mechanism of self-destructive autolysis in phytoplankton under stress.

Gene expression in skeletal tissues: application of laser capture microdissection

Tissue differentiation is based on the expression of transcription factors, receptors for cytokines, and nuclear receptors that regulate a specific phenotype. The purpose of this study was to select cells from various skeletal tissues in order to analyse differential gene expression of cells in the native environment in vivo. It is a difficult task to obtain cells from skeletal tissues, such as cartilage, periost, bone and muscle, that are structured together and do not exist as individual organs. We used laser capture microdissection which permits the selection and isolation of individual cells from tissue sections. The RNA isolated from these tissues was used for reverse transcriptase-polymerase chain reactions for molecular analysis. We analysed the expression of transcription factors (cFOS, cbfa1, MyoD), receptors for cytokines, nuclear receptors, alkaline phosphatase and the structural proteins osteocalcin and collagen II. The results obtained demonstrate differential patterns of gene expression according to the tissue arrangement in their native in vivo environment, with reliable interpretation of the functions of the analysed genes in the context of intact skeletal tissue physiology.

Generation and analysis of 9792 EST sequences from cold acclimated oat, Avena sativa

BACKGROUND

Oat is an important crop in North America and northern Europe. In Scandinavia, yields are limited by the fact that oat cannot be used as a winter crop. In order to develop such a crop, more knowledge about mechanisms of cold tolerance in oat is required.

RESULTS

From an oat cDNA library 9792 single-pass EST sequences were obtained. The library was prepared from pooled RNA samples isolated from leaves of four-week old Avena sativa (oat) plants incubated at +4 degrees C for 4, 8, 16 and 32 hours. Exclusion of sequences shorter than 100 bp resulted in 8508 high-quality ESTs with a mean length of 710.7 bp. Clustering and assembly identified a set of 2800 different transcripts denoted the Avena sativa cold induced UniGene set (AsCIUniGene set). Taking advantage of various tools and databases, putative functions were assigned to 1620 (58%) of these genes. Of the remaining 1180 unclassified sequences, 427 appeared to be oat-specific since they lacked any significant sequence similarity (Blast E values > 10(-10)) to any sequence available in the public databases. Of the 2800 UniGene sequences, 398 displayed significant homology (BlastX E values < or = 10(-10)) to genes previously reported to be involved in cold stress related processes. 107 novel oat transcription factors were also identified, out of which 51 were similar to genes previously shown to be cold induced. The CBF transcription factors have a major role in regulating cold acclimation. Four oat CBF sequences were found, belonging to the monocot cluster of DREB family ERF/AP2 domain proteins. Finally in the total EST sequence data (5.3 Mbp) approximately 400 potential SSRs were found, a frequency similar to what has previously been identified in Arabidopsis ESTs.

CONCLUSION

The AsCIUniGene set will now be used to fabricate an oat biochip, to perform various expression studies with different oat cultivars incubated at varying temperatures, to generate molecular markers and provide tools for various genetic transformation experiments in oat. This will lead to a better understanding of the cellular biology of this important crop and will open up new ways to improve its agronomical properties.

Molecular biology of preimplantation embryos: primer for philosophical discussions

This article is based on a presentation at the First International Conference on Ethics, Science and Moral Philosophy of Assisted Human Reproduction. The goal is to provide scientific background for the discussion of philosophic issues. Recent advances in the systematic molecular analysis of preimplantation embryos are summarized, including the molecular identification of nearly all genes involved in preimplantation development and their detailed expression patterns. Notwithstanding a quantum leap in molecular understanding of preimplantation embryos, molecular evidence seems to provide no decisive definition of a threshold for the beginning of human life during preimplantation development.

cDNA library construction from a small amount of RNA: adaptor-ligation approach for two-round cRNA amplification using T7 and SP6 RNA polymerases

In this study, we developed a method that allows cDNA library construction from a small amount of RNA without causing serious size bias in the resulting cDNA population. For this purpose, we adopted two-round cRNA amplification by T7 and SP6 RNA polymerases. The first-round cDNAs, flanked by the promoter sequences of T7 and SP6 RNA polymerases, were synthesized from 1 µg total RNA and then subjected to two rounds of cRNA amplification. Comparison of the sizes of the first-round and the second-round cRNAs indicated that the size-bias effect of the second-round cRNA synthesis was not serious. The resultant double-stranded cDNAs were cloned into a plasmid by in vitro λ phage recombination with an efficiency of 1.2 × 1011 colony-forming unit/microgram of starting total RNA. Characterization of the resultant cDNA library in terms of the insert size, clone redundancy, and integrity of 3′ ends of cDNAs indicated that the amplified library was comparable to a library constructed by a conventional method, although large cDNAs tend to be slightly truncated in the amplified library. This method enables the construction of a library from a small amount of RNA, and calculations suggest that the strategy would be efficient enough to use even a single cell as starting material.

1274 full-open reading frames of transcripts expressed in the developing mouse nervous system

As part of the trans-National Institutes of Health (NIH) Mouse Brain Molecular Anatomy Project (BMAP), and in close coordination with the NIH Mammalian Gene Collection Program (MGC), we initiated a large-scale project to clone, identify, and sequence the complete open reading frame (ORF) of transcripts expressed in the developing mouse nervous system. Here we report the analysis of the ORF sequence of 1274 cDNAs, obtained from 47 full-length-enriched cDNA libraries, constructed by using a novel approach, herein described. cDNA libraries were derived from size-fractionated cytoplasmic mRNA isolated from brain and eye tissues obtained at several embryonic stages and postnatal days. Altogether, including the full-ORF MGC sequences derived from these libraries by the MGC sequencing team, NIH_BMAP full-ORF sequences correspond to approximately 20% of all transcripts currently represented in mouse MGC. We show that NIH_BMAP clones comprise 68% of mouse MGC cDNAs > or =5 kb, and 54% of those > or =4 kb, as of March 15, 2004. Importantly, we identified transcripts, among the 1274 full-ORF sequences, that are exclusively or predominantly expressed in brain and eye tissues, many of which encode yet uncharacterized proteins.

The NIA cDNA project in mouse stem cells and early embryos

A catalog of mouse genes expressed in early embryos, embryonic and adult stem cells was assembled, including 250000 ESTs, representing approximately 39000 unique transcripts. The cDNA libraries, enriched in full-length clones, were condensed into the NIA 15 and 7.4K clone sets, freely distributed to the research community, providing a standard platform for expression studies using microarrays. They are essential tools for studying mammalian development and stem cell biology, and to provide hints about the differential nature of embryonic and adult stem cells.

Transcriptome analysis of mouse stem cells and early embryos

Understanding and harnessing cellular potency are fundamental in biology and are also critical to the future therapeutic use of stem cells. Transcriptome analysis of these pluripotent cells is a first step towards such goals. Starting with sources that include oocytes, blastocysts, and embryonic and adult stem cells, we obtained 249,200 high-quality EST sequences and clustered them with public sequences to produce an index of approximately 30,000 total mouse genes that includes 977 previously unidentified genes. Analysis of gene expression levels by EST frequency identifies genes that characterize preimplantation embryos, embryonic stem cells, and adult stem cells, thus providing potential markers as well as clues to the functional features of these cells. Principal component analysis identified a set of 88 genes whose average expression levels decrease from oocytes to blastocysts, stem cells, postimplantation embryos, and finally to newborn tissues. This can be a first step towards a possible definition of a molecular scale of cellular potency. The sequences and cDNA clones recovered in this work provide a comprehensive resource for genes functioning in early mouse embryos and stem cells. The nonrestricted community access to the resource can accelerate a wide range of research, particularly in reproductive and regenerative medicine.

250,000 EST sequences from oocytes, blastocysts, and embryonic and adult stem cells contribute to the annotation of the mouse genome and suggest genes that contribute to the unique features of these developmental stages and cell types

Expression profiling of placentomegaly associated with nuclear transplantation of mouse ES cells

Transplantation of nuclei (NT) from engineered mouse ES cells is a potentially powerful and rapid route to create knockout mice, obviating the need for matings to obtain germ-line chimeras. However, such an application is currently impossible, because NT often results in abnormalities in embryo and placenta. Although the epigenetic instability of several imprinted genes in ES cells and ES-derived NT mice has been demonstrated, it is not clear yet what causes the abnormalities. To gain perspective on the extent and types of changes, we have done gene expression profiling for mouse placentas produced by NT of ES cells and compared them with the expression profiles of placentas produced by NT of one-cell embryos. Based on microarray studies with the NIA 15K mouse cDNA collection, we report five principal aberrant events: (1) inappropriate expression of imprinted genes; (2) altered expression of regulatory genes involved in global gene expression, such as DNA methyltransferase and histone acetyltransferase; (3) increased expression of oncogenes and growth promoting genes; (4) overexpression of genes involved in placental growth, such as Plac1; and (5) identification of many novel genes overexpressed in ES-derived NT mouse placentas, including Pitrm1, a new member of the metalloprotease family. The results indicate that placentomegaly in ES-derived NT mice is associated with large-scale dysregulation of normal gene expression patterns. The study also suggests the presence of two regulatory pathways that may lead to histologically discernable placentomegaly. The discovery of groups of genes with altered expression may provide potential targets for intervention to mimic natural regulation more faithfully in NT mice.

Assembly, verification, and initial annotation of the NIA mouse 7.4K cDNA clone set

A set of 7407 cDNA clones (NIA mouse 7.4K) was assembled from >20 cDNA libraries constructed mainly from early mouse embryos, including several stem cell libraries. The clone set was assembled from embryonic and newborn organ libraries consisting of ~120,000 cDNA clones, which were initially re-arrayed into a set of ~11,000 unique cDNA clones. A set of tubes was constructed from the racks in this set to prevent contamination and potential mishandling errors in all further re-arrays. Sequences from this set (11K) were analyzed further for quality and clone identity, and high-quality clones with verified identity were re-arrayed into the final set (7.4K). The set is freely available, and a corresponding database was built to provide comprehensive annotation for those clones with known identity or homology, and has been made available through an extensive Web site that includes many link-outs to external databases and analysis servers.

Full-length-enriched cDNA libraries from Echinococcus granulosus contain separate populations of oligo-capped and trans-spliced transcripts and a high level of predicted signal peptide sequences

The tissue-dwelling larval stages of the cestode Echinococcus granulosus are intimately associated with the host, implying that a range of molecular mediators may be secreted by the parasite into the host environment. These mediators are being sought through a transcriptome-based analysis, using recombinant cDNA libraries. Conventional cDNA libraries of E. granulosus contain high levels of mitochondrial transcripts, as well as host (bovine) genomic DNA. In particular, 60% of a conventional protoscolex stage cDNA library corresponds to the large subunit (LSU) of mitochondrial rRNA. We attribute the presence of LSU rRNA copies to its polyadenylation in E. granulosus. To circumvent this problem, we adapted the 5' Rapid Amplification of cDNA Ends (RNA-ligase mediated RACE) technique that excludes all polynucleotides missing the 7-methyl-guanosine (7MG) cap specific to the 5' end of full-length mRNA. By ligating a specific oligonucleotide (oligo-cap) to 7MG-bearing mRNA, three cDNA libraries were made by PCR from oligo-cap and oligo-dT primers. Analysis of these libraries showed that mitochondrial RNA contaminants had been excluded. Moreover, no bovine genomic sequences were detected. In parallel, we constructed three cDNA libraries using the newly described trans-spliced leader (SL) from Echinococcus. Although these represent a smaller subset of parasite genes, mitochondrial and genomic contributions were again excluded. In both cases, a majority of cDNAs (61-92%) were judged to contain the initiation ATG codon, and 11-27% of inserts included potential N-terminal signal sequences. The 5' UTR tracts of most oligo-capped cDNAs were <100 nt, although approximately 8% were longer than this. Among the trans-spliced cDNAs, 43% potentially utilise the AUG donated by the SL, and in only 6% was the SL separated from an endogenous putative start site by >60 nt. Sequence analysis of randomly selected clones shows virtually no overlap between the oligo-capped and SL libraries, indicating that trans-spliced E. granulosus mRNAs appear to be insensitive to the enzymatic treatments used to 'oligo-cap' unspliced mRNAs. The oligo-capped and SL strategies represent efficient and complementary pathways to isolate full-length cDNA clones from this cestode parasite and, possibly, from related parasitic flatworms.