High-throughput cDNA Sequencing Research Articles

NFuse: Discovery of complex genomic rearrangements in cancer using high-throughput sequencing

Complex genomic rearrangements (CGRs) are emerging as a new feature of cancer genomes. CGRs are characterized by multiple genomic breakpoints and thus have the potential to simultaneously affect multiple genes, fusing some genes and interrupting other genes. Analysis of high-throughput whole-genome shotgun sequencing (WGSS) is beginning to facilitate the discovery and characterization of CGRs, but further development of computational methods is required. We have developed an algorithmic method for identifying CGRs in WGSS data based on shortest alternating paths in breakpoint graphs. Aiming for a method with the highest possible sensitivity, we use breakpoint graphs built from all WGSS data, including sequences with ambiguous genomic origin. Since the majority of cell function is encoded by the transcriptome, we target our search to find CGRs that underlie fusion transcripts predicted from matched high-throughput cDNA sequencing (RNA-seq). We have applied our method, nFuse, to the discovery of CGRs in publicly available data from the well-studied breast cancer cell line HCC1954 and primary prostate tumor sample 963. We first establish the sensitivity and specificity of the nFuse breakpoint prediction and scoring method using breakpoints previously discovered in HCC1954. We then validate five out of six CGRs in HCC1954 and two out of two CGRs in 963. We show examples of gene fusions that would be difficult to discover using methods that do not account for the existence of CGRs, including one important event that was missed in a previous study of the HCC1954 genome. Finally, we illustrate how CGRs may be used to infer the gene expression history of a tumor.

Identification of putative candidate genes involved in cuticle formation in Prunus avium (sweet cherry) fruit.

The cuticular membrane (CM) of Prunus avium (sweet cherry) and other fleshy fruit is under stress. Previous research indicates that the resultant strain promotes microscopic cuticular cracking. Microcracks impair the function of the CM as a barrier against pathogens and uncontrolled water loss/uptake. Stress and strain result from a cessation of CM deposition during early development, while the fruit surface continues to expand. The cessation of CM deposition, in turn, may be related to an early downregulation of CM-related genes. The aims of this study were to identify genes potentially involved in CM formation in sweet cherry fruit and to quantify their expression levels. Fruit growth and CM deposition were quantified weekly from anthesis to maturity and rates of CM deposition were calculated. Sequences of genes expressed in the sweet cherry fruit skin (exocarp) were generated using high-throughput sequencing of cDNA and de novo assembly and analysed using bioinformatics tools. Relative mRNA levels of selected genes were quantified in the exocarp and fruit flesh (mesocarp) weekly using reverse transcriptase-quantitative real-time PCR and compared with the calculated CM deposition rate over time. The rate of CM deposition peaked at 93 (±5) μg per fruit d(-1) about 19 d after anthesis. Based on sequence analyses, 18 genes were selected as potentially involved in CM formation. Selected sweet cherry genes shared up to 100 and 98 % similarity with the respective Prunus persica (peach) and Arabidopsis thaliana genes. Expression of 13 putative CM-related genes was restricted to the exocarp and correlated positively with the CM deposition rate. The results support the view that the cessation of CM deposition during early sweet cherry fruit development is accounted for by a downregulation of genes involved in CM deposition. Genes that merit further investigation include PaWINA, PaWINB, PaLipase, PaLTPG1, PaATT1, PaLCR, PaGPAT4/8, PaLACS2, PaLACS1 and PaCER1.

Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks

Recent advances in high-throughput cDNA sequencing (RNA-seq) can reveal new genes and splice variants and quantify expression genome-wide in a single assay. The volume and complexity of data from RNA-seq experiments necessitate scalable, fast and mathematically principled analysis software. TopHat and Cufflinks are free, open-source software tools for gene discovery and comprehensive expression analysis of high-throughput mRNA sequencing (RNA-seq) data. Together, they allow biologists to identify new genes and new splice variants of known ones, as well as compare gene and transcript expression under two or more conditions. This protocol describes in detail how to use TopHat and Cufflinks to perform such analyses. It also covers several accessory tools and utilities that aid in managing data, including CummeRbund, a tool for visualizing RNA-seq analysis results. Although the procedure assumes basic informatics skills, these tools assume little to no background with RNA-seq analysis and are meant for novices and experts alike. The protocol begins with raw sequencing reads and produces a transcriptome assembly, lists of differentially expressed and regulated genes and transcripts, and publication-quality visualizations of analysis results. The protocol's execution time depends on the volume of transcriptome sequencing data and available computing resources but takes less than 1 d of computer time for typical experiments and ∼1 h of hands-on time.

Transcriptome of Small Regulatory RNAs in the Development of the Zoonotic Parasite Trichinella spiralis

Background Trichinella spiralis is a parasite with unique features. It is a multicellular organism but with an intracellular parasitization and development stage. T. spiralis is the helminthic pathogen that causes zoonotic trichinellosis and afflicts more than 10 million people worldwide, whereas the parasite's biology, especially the developmental regulation is largely unknown. In other organisms, small non-coding RNAs, such as microRNAs (miRNA) and small interfering RNAs (siRNA) execute post-transcriptional regulation by translational repression or mRNA degradation, and a large number of miRNAs have been identified in diverse species. In T. spiralis, the profile of small non-coding RNAs and their function remains poorly understood.Methodology and Principal FindingsHere, the transcriptional profiles of miRNA and siRNA in three developmental stages of T. spiralis in the rat host were investigated, and compared by high-throughput cDNA sequencing technique (“RNA-seq”). 5,443,641 unique sequence tags were obtained. Of these, 21 represented conserved miRNAs related to 13 previously identified metazoan miRNA families and 213 were novel miRNAs so far unique to T. spiralis. Some of these miRNAs exhibited stage-specific expression. Expression of miRNAs was confirmed in three stages of the life cycle by qRT-PCR and northern blot analysis. In addition, endogenous siRNAs (endo-siRNAs) were found mainly derived from natural antisense transcripts (NAT) and transposable elements (TE) in the parasite.Conclusions and SignificanceWe provide evidence for the presence of miRNAs and endo-siRNAs in T. spiralis. The miRNAs accounted for the major proportion of the small regulatory RNA population of T. spiralis, while fewer endogenous siRNAs were found. The finding of stage-specific expression patterns of the miRNAs in different developmental stages of T. spiralis suggests that miRNAs may play important roles in parasite development. Our data provide a basis for further understanding of the molecular regulation and functional evolution of miRNAs in parasitic nematodes.

A powerful and flexible statistical framework for testing hypotheses of allele-specific gene expression from RNA-seq data.

Variation in gene expression is thought to make a significant contribution to phenotypic diversity among individuals within populations. Although high-throughput cDNA sequencing offers a unique opportunity to delineate the genome-wide architecture of regulatory variation, new statistical methods need to be developed to capitalize on the wealth of information contained in RNA-seq data sets. To this end, we developed a powerful and flexible hierarchical Bayesian model that combines information across loci to allow both global and locus-specific inferences about allele-specific expression (ASE). We applied our methodology to a large RNA-seq data set obtained in a diploid hybrid of two diverse Saccharomyces cerevisiae strains, as well as to RNA-seq data from an individual human genome. Our statistical framework accurately quantifies levels of ASE with specified false-discovery rates, achieving high reproducibility between independent sequencing platforms. We pinpoint loci that show unusual and biologically interesting patterns of ASE, including allele-specific alternative splicing and transcription termination sites. Our methodology provides a rigorous, quantitative, and high-resolution tool for profiling ASE across whole genomes.

Synthetic spike-in standards for RNA-seq experiments

High-throughput sequencing of cDNA (RNA-seq) is a widely deployed transcriptome profiling and annotation technique, but questions about the performance of different protocols and platforms remain. We used a newly developed pool of 96 synthetic RNAs with various lengths, and GC content covering a 2(20) concentration range as spike-in controls to measure sensitivity, accuracy, and biases in RNA-seq experiments as well as to derive standard curves for quantifying the abundance of transcripts. We observed linearity between read density and RNA input over the entire detection range and excellent agreement between replicates, but we observed significantly larger imprecision than expected under pure Poisson sampling errors. We use the control RNAs to directly measure reproducible protocol-dependent biases due to GC content and transcript length as well as stereotypic heterogeneity in coverage across transcripts correlated with position relative to RNA termini and priming sequence bias. These effects lead to biased quantification for short transcripts and individual exons, which is a serious problem for measurements of isoform abundances, but that can partially be corrected using appropriate models of bias. By using the control RNAs, we derive limits for the discovery and detection of rare transcripts in RNA-seq experiments. By using data collected as part of the model organism and human Encyclopedia of DNA Elements projects (ENCODE and modENCODE), we demonstrate that external RNA controls are a useful resource for evaluating sensitivity and accuracy of RNA-seq experiments for transcriptome discovery and quantification. These quality metrics facilitate comparable analysis across different samples, protocols, and platforms.

A Genome-Wide Approach to Discovery of Small RNAs Involved in Regulation of Virulence in Vibrio cholerae

Small RNAs (sRNAs) are becoming increasingly recognized as important regulators in bacteria. To investigate the contribution of sRNA mediated regulation to virulence in Vibrio cholerae, we performed high throughput sequencing of cDNA generated from sRNA transcripts isolated from a strain ectopically expressing ToxT, the major transcriptional regulator within the virulence gene regulon. We compared this data set with ToxT binding sites determined by pulldown and deep sequencing to identify sRNA promoters directly controlled by ToxT. Analysis of the resulting transcripts with ToxT binding sites in cis revealed two sRNAs within the Vibrio Pathogenicity Island. When deletions of these sRNAs were made and the resulting strains were competed against the parental strain in the infant mouse model of V. cholerae colonization, one, TarB, displayed a variable colonization phenotype dependent on its physiological state at the time of inoculation. We identified a target of TarB as the mRNA for the secreted colonization factor, TcpF. We verified negative regulation of TcpF expression by TarB and, using point mutations that disrupted interaction between TarB and tpcF mRNA, showed that loss of this negative regulation was primarily responsible for the colonization phenotype observed in the TarB deletion mutant.

Effects of chronic metal exposure on wild fish populations revealed by high-throughput cDNA sequencing

Given the inherent variability of aquatic systems, predicting the in situ effects of contaminants on such ecosystems still represents a major challenge for ecotoxicology. In this context, transcriptomic tools can help identify and investigate the mechanisms of toxicity beyond the traditional morphometric, physiological and population-level endpoints. In this study, we used the 454 sequencing technology to examine the in situ effects of chronic metal (Cd, Cu) exposure on the yellow perch (Perca flavescens) transcriptome. Total hepatic mRNA from fish sampled along a polymetallic gradient was extracted, reverse transcribed, labeled with unique barcode sequences and sequenced. This approach allowed us to identify correlations between the transcription level of single genes and the hepatic concentrations of individual metals; 71% of the correlations established were negative. Chronic metal exposure was thus associated with a decrease in the transcription levels of numerous genes involved in protein biosynthesis, in the immune system, and in lipid and energy metabolism. Our results suggest that this marked decrease could result from an impairment of bile acid metabolism by Cd and energy restriction but also from the recruitment of several genes involved in epigenetic modifications of histones and DNA that lead to gene silencing.

Comprehensive annotation of the transcriptome of the human fungal pathogen Candida albicans using RNA-seq

Candida albicans is the major invasive fungal pathogen of humans, causing diseases ranging from superficial mucosal infections to disseminated, systemic infections that are often lifethreatening. We have used massively parallel high-throughput sequencing of cDNA (RNA-seq) to generate a high-resolution map of the C. albicans transcriptome under several different environmental conditions. We have quantitatively determined all of the regions that are transcribed under these different conditions, and have identified 602 novel transcriptionally active regions (TARs) and numerous novel introns that are not represented in the current genome annotation. Interestingly, the expression of many of these TARs is regulated in a condition-specific manner. This comprehensive transcriptome analysis significantly enhances the current genome annotation of C. albicans, a necessary framework for a complete understanding of the molecular mechanisms of pathogenesis for this important eukaryotic pathogen.

Studies of a Biochemical Factory: Tomato Trichome Deep Expressed Sequence Tag Sequencing and Proteomics

Shotgun proteomics analysis allows hundreds of proteins to be identified and quantified from a single sample at relatively low cost. Extensive DNA sequence information is a prerequisite for shotgun proteomics, and it is ideal to have sequence for the organism being studied rather than from related species or accessions. While this requirement has limited the set of organisms that are candidates for this approach, next generation sequencing technologies make it feasible to obtain deep DNA sequence coverage from any organism. As part of our studies of specialized (secondary) metabolism in tomato (Solanum lycopersicum) trichomes, 454 sequencing of cDNA was combined with shotgun proteomics analyses to obtain in-depth profiles of genes and proteins expressed in leaf and stem glandular trichomes of 3-week-old plants. The expressed sequence tag and proteomics data sets combined with metabolite analysis led to the discovery and characterization of a sesquiterpene synthase that produces beta-caryophyllene and alpha-humulene from E,E-farnesyl diphosphate in trichomes of leaf but not of stem. This analysis demonstrates the utility of combining high-throughput cDNA sequencing with proteomics experiments in a target tissue. These data can be used for dissection of other biochemical processes in these specialized epidermal cells.

A simple method for directional transcriptome sequencing using Illumina technology

High-throughput sequencing of cDNA has been used to study eukaryotic transcription on a genome-wide scale to single base pair resolution. In order to compensate for the high ribonuclease activity in bacterial cells, we have devised an equivalent technique optimized for studying complete prokaryotic transcriptomes that minimizes the manipulation of the RNA sample. This new approach uses Illumina technology to sequence single-stranded (ss) cDNA, generating information on both the direction and level of transcription throughout the genome. The protocol, and associated data analysis programs, are freely available from http://www.sanger.ac.uk/Projects/Pathogens/Transcriptome/. We have successfully applied this method to the bacterial pathogens Salmonella bongori and Streptococcus pneumoniae and the yeast Schizosaccharomyces pombe. This method enables experimental validation of genetic features predicted in silico and allows the easy identification of novel transcripts throughout the genome. We also show that there is a high correlation between the level of gene expression calculated from ss-cDNA and double-stranded-cDNA sequencing, indicting that ss-cDNA sequencing is both robust and appropriate for use in quantitative studies of transcription. Hence, this simple method should prove a useful tool in aiding genome annotation and gene expression studies in both prokaryotes and eukaryotes.

Analysis of 5′ transcript heterogeneity by high-throughput sequencing of cDNA

We isolated mRNA from Blumeria graminis and prepared a cDNA library using a modification of Epicentre's ExactSTART™ Full-Length cDNA Library Cloning kit. Analysis of the 5′ ends of the cDNA by 454 pyrosequencing yielded approximately 250,000 expressed sequence tags (ESTs) from one run of sequencing. The data also showed marked heterogeneity of the 5′ ends of the transcripts, including the addition of non–template-encoded bases.

Analysis of 5′ transcript heterogeneity by high-throughput sequencing of cDNA

Analysis of 5′ transcript heterogeneity by high-throughput sequencing of cDNA

A pig multi-tissue normalised cDNA library: large-scale sequencing, cluster analysis and 9K micro-array resource generation

BackgroundDomestic animal breeding and product quality improvement require the control of reproduction, nutrition, health and welfare in these animals. It is thus necessary to improve our knowledge of the major physiological functions and their interactions. This would be greatly enhanced by the availability of expressed gene sequences in the databases and by cDNA arrays allowing the transcriptome analysis of any function.The objective within the AGENAE French program was to initiate a high-throughput cDNA sequencing program of a 38-tissue normalised library and generate a diverse microarray for transcriptome analysis in pig species.ResultsWe constructed a multi-tissue cDNA library, which was normalised and subtracted to reduce the redundancy of the clones. Expressed Sequence Tags were produced and 24449 high-quality sequences were released in EMBL database. The assembly of all the public ESTs (available through SIGENAE website) resulted in 40786 contigs and 54653 singletons. At least one Agenae sequence is present in 11969 contigs (12.5%) and in 9291 of the deeper-than-one-contigs (22.8%). Sequence analysis showed that both normalisation and subtraction processes were successful and that the initial tissue complexity was maintained in the final libraries. A 9K nylon cDNA microarray was produced and is available through CRB-GADIE. It will allow high sensitivity transcriptome analyses in pigs.ConclusionIn the present work, a pig multi-tissue cDNA library was constructed and a 9K cDNA microarray designed. It contributes to the Expressed Sequence Tags pig data, and offers a valuable tool for transcriptome analysis.

Genome-Wide Expression Profiling of the Arabidopsis Female Gametophyte Identifies Families of Small, Secreted Proteins

The female gametophyte of flowering plants, the embryo sac, develops within the diploid (sporophytic) tissue of the ovule. While embryo sac–expressed genes are known to be required at multiple stages of the fertilization process, the set of embryo sac–expressed genes has remained poorly defined. In particular, the set of genes responsible for mediating intracellular communication between the embryo sac and the male gametophyte, the pollen grain, is unknown. We used high-throughput cDNA sequencing and whole-genome tiling arrays to compare gene expression in wild-type ovules to that in dif1 ovules, which entirely lack embryo sacs, and myb98 ovules, which are impaired in pollen tube attraction. We identified nearly 400 genes that are downregulated in dif1 ovules. Seventy-eight percent of these embryo sac–dependent genes were predicted to encode for secreted proteins, and 60% belonged to multigenic families. Our results define a large number of candidate extracellular signaling molecules that may act during embryo sac development or fertilization; less than half of these are represented on the widely used ATH1 expression array. In particular, we found that 37 out of 40 genes encoding Domain of Unknown Function 784 (DUF784) domains require the synergid-specific transcription factor MYB98 for expression. Several DUF784 genes were transcribed in synergid cells of the embryo sac, implicating the DUF784 gene family in mediating late stages of embryo sac development or interactions with pollen tubes. The coexpression of highly similar proteins suggests a high degree of functional redundancy among embryo sac genes.

Rarity of Somatic Mutation and Frequency of Normal Sequence Variation Detected in Sporadic Colon Adenocarcinoma Using High-Throughput cDNA Sequencing

We performed high-throughput cDNA sequencing in colorectal adenocarcinoma and matching normal colorectal epithelium. All six hundred three genes in the UCSC database that were expressed in colon cancers and contained open reading frames of 1000 nucleotides or less were selected for study (total basepairs/bp, 366,686). 304,350 of these 366,686 bp (83.0%) were amplified and sequenced successfully. Seventy-eight sequence variants present in germline (i.e. normal) as well as matching somatic (i.e. tumor) DNA were discovered, yielding a frequency of 1 variant per 3,902 bp. Fifty-one of these sequence variants were homozygous (26 synonymous, 25 non-synonymous), while 27 were heterozygous (11 synonymous, 16 non-synonymous). Cancer tissue contained only one sequence-altered allele of the gene ATP50, which was present heterozygously alongside the wild-type allele in matching normal epithelium. Despite this relatively large number of bp and genes sequenced, no somatic mutations unique to tumor were found. High-throughput cDNA sequencing is a practical approach for detecting novel sequence variations and alterations in human tumors, such as those of the colon.

Identification of differently expressed genes in human colorectal adenocarcinoma

To investigate the differently expressed genes in human colorectal adenocarcinoma. The integrated approach for gene expression profiling that couples suppression subtractive hybridization, high-throughput cDNA array, sequencing, bioinformatics analysis, and reverse transcriptase real-time quantitative polymerase chain reaction (PCR) was carried out. A set of cDNA clones including 1260 SSH inserts amplified by PCR was arrayed using robotic printing. The cDNA arrays were hybridized with fluorescent-labeled probes prepared from RNA of human colorectal adenocarcinoma (HCRAC) and normal colorectal tissues. A total of 86 genes were identified, 16 unknown genes and 70 known genes. The transcription factor Somult9 influencing cell differentiation was downregulated. At the same time, Heat shock protein 10 KDis downregulated and Calmoulin is up-regulated. Downregulation of heat shock protein 10 KD lost its inhibition of Ras, and then attenuated the Ras GTPase signaling pathway, increased cell proliferation and inhibited cell apoptosis. Down-regulated transcription factor Somult9 influences cell differentiation and cell-specific gene expression. Down-regulated Somult9 also decreases its binding to calmodulin, accumulates calmodulin as receptor-activated kinase and phosphorylase kinase due to the activation of PhK.

Concatenation cDNA sequencing for transcriptome analysis

We describe a high-throughput cDNA sequencing pipeline ( http://www.hgsc.bcm.tmc.edu/projects/cdna) built in response to the emerging need for rapid sequencing of large cDNA collections. Using this strategy cDNA inserts are purified and joined through concatenation into large molecules. These ‘pseudo-BACs’ are subjected to random shotgun sequencing whereby the majority of cDNA inserts in the pool are sequenced. Using this concatenation cDNA sequencing platform, we have contributed more than 13 000 full-length cDNA sequences from human and mouse to the Mammalian Gene Collection (MGC). To cite this article: P.H. Gunaratne et al., C. R. Biologies 326 (2003).

Molecular cloning and characterization of a novel human kinase gene, PDIK1L.

We isolated a 4301-bp cDNA from a human foetal brain cDNA library by high-throughput cDNA sequencing. It encodes a protein of 341 amino acids, which shows 69% identity with the human kinase CLIK1 (AAL99353), which was suggested to be the CLP-36 interacting kinase. Bioinformatics analysis suggests that the putative kinase may interact with PDZ and LIM domain proteins. Therefore the protein and its cDNA were named 'PDLIM1 interacting kinase 1 like' (PDIK1L; nomenclature approved by the HUGO Gene Nomenclature Committee). Ensembl Genome Browser located PDIK1L to human chromosome 1p35.3. It spans about 13.7 kb and consists of four exons and three introns. Multiple-tissue cDNA panel PCR revealed that the gene is expressed widely in human tissues: liver, kidney, pancreas, spleen, thymus and prostate. The protein appears to be localized to the nucleus.

Cloning and identification of a novel cDNA which encodes a putative protein with a DnaJ domain and a thioredoxin active motif, human macrothioredoxin.

A 3345 bp cDNA was isolated from the fetal brain cDNA library by high throughput cDNA sequencing. The cDNA with an open reading fragment (ORF) of 2241 bp encodes a 747 amino acids putative protein with a DnaJ N-terminus domain and four thioredoxin active sets. So it is named human macrothioredoxin (hMTHr). We used Northern blot to detect a band with a length of about 5 kb, which was ubiquitously expressed in human adult tissues with different intensities. The expression pattern was verified by RT-PCR, revealing that the transcripts were ubiquitously expressed in fetal tissues and human tumor tissues also. The transcripts in fetal tissues were more numerous than in adult tissues. The transcripts were high in adult testis, adult pancreas, fetal thymus, and fetal kidney. We have found three splice isoforms with lengths of 3483, 3345, and 2982 bp in the sequencing analysis of the RT-PCR products. They encode three putative proteins with 793, 747, and 275 amino acids, respectively. The first two putative proteins contain a DnaJ domain and four thioredoxin domains. The third putative protein only contains an N-terminus DnaJ domain and one thioredoxin active set. Blast analysis against the NCBI database revealed that the gene spans a genome sequence of about 75 kb that contains at least 25 exons and is located in human chromosome 2q32.1. The organization of the functional motifs of hMTHr suggests that the protein might be a member of a molecular chaperone family.