Tissue-specific transcript annotation and expression profiling with complementary next-generation sequencing technologies

Mapping and modeling the genomic basis of differential RNA isoform expression at single-cell resolution with LR-Split-seq

The rise in throughput and quality of long-read sequencing should allow unambiguous identification of full-length transcript isoforms. However, its application to single-cell RNA-seq has been limited by throughput and expense. Here we develop and characterize long-read Split-seq (LR-Split-seq), which uses combinatorial barcoding to sequence single cells with long reads. Applied to the C2C12 myogenic system, LR-split-seq associates isoforms to cell types with relative economy and design flexibility. We find widespread evidence of changing isoform expression during differentiation including alternative transcription start sites (TSS) and/or alternative internal exon usage. LR-Split-seq provides an affordable method for identifying cluster-specific isoforms in single cells.

Altered visual processing in the mdx52 mouse model of Duchenne muscular dystrophy

The mdx52 mouse model of Duchenne muscular dystrophy (DMD) is lacking exon 52 of the DMD gene that is located in a hotspot mutation region causing cognitive deficits and retinal anomalies in DMD patients. This deletion leads to the loss of the dystrophin proteins, Dp427, Dp260 and Dp140, while Dp71 is preserved. The flash electroretinogram (ERG) in mdx52 mice was previously characterized by delayed dark-adapted b-waves. A detailed description of functional ERG changes and visual performances in mdx52 mice is, however, lacking. Here an extensive full-field ERG repertoire was applied in mdx52 mice and WT littermates to analyze retinal physiology in scotopic, mesopic and photopic conditions in response to flash, sawtooth and/or sinusoidal stimuli. Behavioral contrast sensitivity was assessed using quantitative optomotor response (OMR) to sinusoidally modulated luminance gratings at 100% or 50% contrast. The mdx52 mice exhibited reduced amplitudes and delayed implicit times in dark-adapted ERG flash responses, particularly in their b-wave and oscillatory potentials, and diminished amplitudes of light-adapted flash ERGs. ERG responses to sawtooth stimuli were also diminished and delayed for both mesopic and photopic conditions in mdx52 mice and the first harmonic amplitudes to photopic sine-wave stimuli were smaller at all temporal frequencies. OMR indices were comparable between genotypes at 100% contrast but significantly reduced in mdx52 mice at 50% contrast. The complex ERG alterations and disturbed contrast vision in mdx52 mice include features observed in DMD patients and suggest altered photoreceptor-to-bipolar cell transmission possibly affecting contrast sensitivity. The mdx52 mouse is a relevant model to appraise the roles of retinal dystrophins and for preclinical studies related to DMD.

Neutral endopeptidase inhibitors blunt kidney fibrosis by reducing myofibroblast formation

Kidney fibrosis is the common pathophysiological mechanism in end-stage renal disease characterized by excessive accumulation of myofibroblast-derived extracellular matrix. Natriuretic peptides have been demonstrated to have cyclic guanosine monophosphate (cGMP)-dependent anti-fibrotic properties likely due to interference with pro-fibrotic tissue growth factor β (TGF-β) signaling. However, in vivo, natriuretic peptides are rapidly degraded by neutral endopeptidases (NEP). In a unilateral ureteral obstruction (UUO) mouse model for kidney fibrosis we assessed the anti-fibrotic effects of SOL1, an orally active compound that inhibits NEP and endothelin-converting enzyme (ECE). Mice (n=10 per group) subjected to UUO were treated for 1 week with either solvent, NEP-/ECE-inhibitor SOL1 (two doses), reference NEP-inhibitor candoxatril or the angiotensin II receptor type 1 (AT₁₎-antagonist losartan. While NEP-inhibitors had no significant effect on blood pressure, they did increase urinary cGMP levels as well as endothelin-1 (ET-1) levels. Immunohistochemical staining revealed a marked decrease in renal collagen (∼55% reduction, P<0.05) and α-smooth muscle actin (α-SMA; ∼40% reduction, P<0.05). Moreover, the number of α-SMA positive cells in the kidneys of SOL1-treated groups inversely correlated with cGMP levels consistent with a NEP-dependent anti-fibrotic effect. To dissect the molecular mechanisms associated with the anti-fibrotic effects of NEP inhibition, we performed a 'deep serial analysis of gene expression (Deep SAGE)' transcriptome and targeted metabolomics analysis of total kidneys of all treatment groups. Pathway analyses linked increased cGMP and ET-1 levels with decreased nuclear receptor signaling (peroxisome proliferator-activated receptor [PPAR] and liver X receptor/retinoid X receptor [LXR/RXR] signaling) and actin cytoskeleton organization. Taken together, although our transcriptome and metabolome data indicate metabolic dysregulation, our data support the therapeutic potential of NEP inhibition in the treatment of kidney fibrosis via cGMP elevation and reduced myofibroblast formation.

Expression of murine muscle-enriched A-type lamin-interacting protein (MLIP) is regulated by tissue-specific alternative transcription start sites

Muscle-enriched lamin-interacting protein (Mlip) is an alternatively spliced gene whose splicing specificity is dictated by tissue type. MLIP is most abundantly expressed in brain, cardiac, and skeletal muscle. In the present study, we systematically mapped the transcriptional start and stop sites of murine Mlip Rapid amplification of cDNA ends (RACE) of Mlip transcripts from the brain, heart, and skeletal muscle revealed two transcriptional start sites (TSSs), exon 1a and exon 1b, and only one transcriptional termination site. RT-PCR analysis of the usage of the two identified TSSs revealed that the heart utilizes only exon 1a for MLIP expression, whereas the brain exclusively uses exon 1b TSS. Loss of Mlip exon 1a in mice resulted in a 7-fold increase in the prevalence of centralized nuclei in muscle fibers with the Mlip exon1a-deficient satellite cells on single fibers exhibiting a significant delay in commitment to a MYOD-positive phenotype. Furthermore, we demonstrate that the A-type lamin-binding domain in MLIP is encoded in exon 1a, indicating that MLIP isoforms generated with exon 1b TSS lack the A-type lamin-binding domain. Collectively these findings suggest that Mlip tissue-specific expression and alternative splicing play a critical role in determining MLIP's functions in mice.

Comprehensive comparative analysis of 5'-end RNA-sequencing methods

Specialized RNA-seq methods are required to identify the 5' ends of transcripts, which are critical for studies of gene regulation, but these methods have not been systematically benchmarked. We directly compared six such methods, including the performance of five methods on a single human cellular RNA sample and a new spike-in RNA assay that helps circumvent challenges resulting from uncertainties in annotation and RNA processing. We found that the 'cap analysis of gene expression' (CAGE) method performed best for mRNA and that most of its unannotated peaks were supported by evidence from other genomic methods. We applied CAGE to eight brain-related samples and determined sample-specific transcription start site (TSS) usage, as well as a transcriptome-wide shift in TSS usage between fetal and adult brain.

Undiscovered Physiology of Transcript and Protein Networks

The past two decades have witnessed a rapid evolution in our ability to measure RNA and protein from biological systems. As a result, new principles have arisen regarding how information is processed in cells, how decisions are made, and the role of networks in biology. This essay examines this technological evolution, reviewing (and critiquing) the conceptual framework that has emerged to explain how RNA and protein networks control cellular function. We identify how future investigations into transcriptomes, proteomes, and other cellular networks will enable development of more robust, quantitative models of cellular behavior whilst also providing new avenues to use knowledge of biological networks to improve human health. © 2016 American Physiological Society. Compr Physiol 6:1851-1872, 2016.

Remodeling of retrotransposon elements during epigenetic induction of adult visual cortical plasticity by HDAC inhibitors

BACKGROUND

The capacity for plasticity in the adult brain is limited by the anatomical traces laid down during early postnatal life. Removing certain molecular brakes, such as histone deacetylases (HDACs), has proven to be effective in recapitulating juvenile plasticity in the mature visual cortex (V1). We investigated the chromatin structure and transcriptional control by genome-wide sequencing of DNase I hypersensitive sites (DHSS) and cap analysis of gene expression (CAGE) libraries after HDAC inhibition by valproic acid (VPA) in adult V1.

RESULTS

We found that VPA reliably reactivates the critical period plasticity and induces a dramatic change of chromatin organization in V1 yielding significantly greater accessibility distant from promoters, including at enhancer regions. VPA also induces nucleosome eviction specifically from retrotransposon (in particular SINE) elements. The transiently accessible SINE elements overlap with transcription factor-binding sites of the Fox family. Mapping of transcription start site activity using CAGE revealed transcription of epigenetic and neural plasticity-regulating genes following VPA treatment, which may help to re-program the genomic landscape and reactivate plasticity in the adult cortex.

CONCLUSIONS

Treatment with HDAC inhibitors increases accessibility to enhancers and repetitive elements underlying brain-specific gene expression and reactivation of visual cortical plasticity.

Complementing tissue characterization by integrating transcriptome profiling from the Human Protein Atlas and from the FANTOM5 consortium

Understanding the normal state of human tissue transcriptome profiles is essential for recognizing tissue disease states and identifying disease markers. Recently, the Human Protein Atlas and the FANTOM5 consortium have each published extensive transcriptome data for human samples using Illumina-sequenced RNA-Seq and Heliscope-sequenced CAGE. Here, we report on the first large-scale complex tissue transcriptome comparison between full-length versus 5'-capped mRNA sequencing data. Overall gene expression correlation was high between the 22 corresponding tissues analyzed (R > 0.8). For genes ubiquitously expressed across all tissues, the two data sets showed high genome-wide correlation (91% agreement), with differences observed for a small number of individual genes indicating the need to update their gene models. Among the identified single-tissue enriched genes, up to 75% showed consensus of 7-fold enrichment in the same tissue in both methods, while another 17% exhibited multiple tissue enrichment and/or high expression variety in the other data set, likely dependent on the cell type proportions included in each tissue sample. Our results show that RNA-Seq and CAGE tissue transcriptome data sets are highly complementary for improving gene model annotations and highlight biological complexities within tissue transcriptomes. Furthermore, integration with image-based protein expression data is highly advantageous for understanding expression specificities for many genes.

NanoCAGE-XL and CapFilter: an approach to genome wide identification of high confidence transcription start sites

BACKGROUND

Identifying the transcription start sites (TSS) of genes is essential for characterizing promoter regions. Several protocols have been developed to capture the 5' end of transcripts via Cap Analysis of Gene Expression (CAGE) or linker-ligation strategies such as Paired-End Analysis of Transcription Start Sites (PEAT), but often require large amounts of tissue. More recently, nanoCAGE was developed for sequencing on the Illumina GAIIx to overcome these difficulties.

RESULTS

Here we present the first publicly available adaptation of nanoCAGE for sequencing on recent ultra-high throughput platforms such as Illumina HiSeq-2000, and CapFilter, a computational pipeline that greatly increases confidence in TSS identification. We report excellent gene coverage, reproducibility, and precision in transcription start site discovery for samples from Arabidopsis thaliana roots.

CONCLUSION

nanoCAGE-XL together with CapFilter allows for genome wide identification of high confidence transcription start sites in large eukaryotic genomes.

Annotation of the Protein Coding Regions of the Equine Genome

Current gene annotation of the horse genome is largely derived from in silico predictions and cross-species alignments. Only a small number of genes are annotated based on equine EST and mRNA sequences. To expand the number of equine genes annotated from equine experimental evidence, we sequenced mRNA from a pool of forty-three different tissues. From these, we derived the structures of 68,594 transcripts. In addition, we identified 301,829 positions with SNPs or small indels within these transcripts relative to EquCab2. Interestingly, 780 variants extend the open reading frame of the transcript and appear to be small errors in the equine reference genome, since they are also identified as homozygous variants by genomic DNA resequencing of the reference horse. Taken together, we provide a resource of equine mRNA structures and protein coding variants that will enhance equine and cross-species transcriptional and genomic comparisons.

Assessing the translational landscape of myogenic differentiation by ribosome profiling

The formation of skeletal muscles is associated with drastic changes in protein requirements known to be safeguarded by tight control of gene transcription and mRNA processing. The contribution of regulation of mRNA translation during myogenesis has not been studied so far. We monitored translation during myogenic differentiation of C2C12 myoblasts, using a simplified protocol for ribosome footprint profiling. Comparison of ribosome footprints to total RNA showed that gene expression is mostly regulated at the transcriptional level. However, a subset of transcripts, enriched for mRNAs encoding for ribosomal proteins, was regulated at the level of translation. Enrichment was also found for specific pathways known to regulate muscle biology. We developed a dedicated pipeline to identify translation initiation sites (TISs) and discovered 5333 unannotated TISs, providing a catalog of upstream and alternative open reading frames used during myogenesis. We identified 298 transcripts with a significant switch in TIS usage during myogenesis, which was not explained by alternative promoter usage, as profiled by DeepCAGE. Also these transcripts were enriched for ribosomal protein genes. This study demonstrates that differential mRNA translation controls protein expression of specific subsets of genes during myogenesis. Experimental protocols, analytical workflows, tools and data are available through public repositories (http://lumc.github.io/ribosome-profiling-analysis-framework/).

Alternative mRNA transcription, processing, and translation: insights from RNA sequencing

The human transcriptome comprises >80,000 protein-coding transcripts and the estimated number of proteins synthesized from these transcripts is in the range of 250,000 to 1 million. These transcripts and proteins are encoded by less than 20,000 genes, suggesting extensive regulation at the transcriptional, post-transcriptional, and translational level. Here we review how RNA sequencing (RNA-seq) technologies have increased our understanding of the mechanisms that give rise to alternative transcripts and their alternative translation. We highlight four different regulatory processes: alternative transcription initiation, alternative splicing, alternative polyadenylation, and alternative translation initiation. We discuss their transcriptome-wide distribution, their impact on protein expression, their biological relevance, and the possible molecular mechanisms leading to their alternative regulation. We conclude with a discussion of the coordination and the interdependence of these four regulatory layers.

RNA sequencing: from tag-based profiling to resolving complete transcript structure

Technological advances in the sequencing field support in-depth characterization of the transcriptome. Here, we review genome-wide RNA sequencing methods used to investigate specific aspects of gene expression and its regulation, from transcription to RNA processing and translation. We discuss tag-based methods for studying transcription, alternative initiation and polyadenylation events, shotgun methods for detection of alternative splicing, full-length RNA sequencing for the determination of complete transcript structures, and targeted methods for studying the process of transcription and translation. With the ensemble of technologies available, it is now possible to obtain a comprehensive view on transcriptome complexity and the regulation of transcript diversity.

DeepSAGE reveals genetic variants associated with alternative polyadenylation and expression of coding and non-coding transcripts

Many disease-associated variants affect gene expression levels (expression quantitative trait loci, eQTLs) and expression profiling using next generation sequencing (NGS) technology is a powerful way to detect these eQTLs. We analyzed 94 total blood samples from healthy volunteers with DeepSAGE to gain specific insight into how genetic variants affect the expression of genes and lengths of 3′-untranslated regions (3′-UTRs). We detected previously unknown cis-eQTL effects for GWAS hits in disease- and physiology-associated traits. Apart from cis-eQTLs that are typically easily identifiable using microarrays or RNA-sequencing, DeepSAGE also revealed many cis-eQTLs for antisense and other non-coding transcripts, often in genomic regions containing retrotransposon-derived elements. We also identified and confirmed SNPs that affect the usage of alternative polyadenylation sites, thereby potentially influencing the stability of messenger RNAs (mRNA). We then combined the power of RNA-sequencing with DeepSAGE by performing a meta-analysis of three datasets, leading to the identification of many more cis-eQTLs. Our results indicate that DeepSAGE data is useful for eQTL mapping of known and unknown transcripts, and for identifying SNPs that affect alternative polyadenylation. Because of the inherent differences between DeepSAGE and RNA-sequencing, our complementary, integrative approach leads to greater insight into the molecular consequences of many disease-associated variants.

Many genetic variants that are associated with diseases also affect gene expression levels. We used a next generation sequencing approach targeting 3′ transcript ends (DeepSAGE) to gain specific insight into how genetic variants affect the expression of genes and the usage and length of 3′-untranslated regions. We detected many associations for antisense and other non-coding transcripts, often in genomic regions containing retrotransposon-derived elements. Some of these variants are also associated with disease. We also identified and confirmed variants that affect the usage of alternative polyadenylation sites, thereby potentially influencing the stability of mRNAs. We conclude that DeepSAGE is useful for detecting eQTL effects on both known and unknown transcripts, and for identifying variants that affect alternative polyadenylation.

Quantitative biology of single neurons

The building blocks of complex biological systems are single cells. Fundamental insights gained from single-cell analysis promise to provide the framework for understanding normal biological systems development as well as the limits on systems/cellular ability to respond to disease. The interplay of cells to create functional systems is not well understood. Until recently, the study of single cells has concentrated primarily on morphological and physiological characterization. With the application of new highly sensitive molecular and genomic technologies, the quantitative biochemistry of single cells is now accessible.

Deep sequencing of RNA from immune cell-derived vesicles uncovers the selective incorporation of small non-coding RNA biotypes with potential regulatory functions

Cells release RNA-carrying vesicles and membrane-free RNA/protein complexes into the extracellular milieu. Horizontal vesicle-mediated transfer of such shuttle RNA between cells allows dissemination of genetically encoded messages, which may modify the function of target cells. Other studies used array analysis to establish the presence of microRNAs and mRNA in cell-derived vesicles from many sources. Here, we used an unbiased approach by deep sequencing of small RNA released by immune cells. We found a large variety of small non-coding RNA species representing pervasive transcripts or RNA cleavage products overlapping with protein coding regions, repeat sequences or structural RNAs. Many of these RNAs were enriched relative to cellular RNA, indicating that cells destine specific RNAs for extracellular release. Among the most abundant small RNAs in shuttle RNA were sequences derived from vault RNA, Y-RNA and specific tRNAs. Many of the highly abundant small non-coding transcripts in shuttle RNA are evolutionary well-conserved and have previously been associated to gene regulatory functions. These findings allude to a wider range of biological effects that could be mediated by shuttle RNA than previously expected. Moreover, the data present leads for unraveling how cells modify the function of other cells via transfer of specific non-coding RNA species.

Survey of the Applications of NGS to Whole-Genome Sequencing and Expression Profiling

Recently, the technologies of DNA sequence variation and gene expression profiling have been used widely as approaches in the expertise of genome biology and genetics. The application to genome study has been particularly developed with the introduction of the next-generation DNA sequencer (NGS) Roche/454 and Illumina/Solexa systems, along with bioinformation analysis technologies of whole-genome de novo assembly, expression profiling, DNA variation discovery, and genotyping. Both massive whole-genome shotgun paired-end sequencing and mate paired-end sequencing data are important steps for constructing de novo assembly of novel genome sequencing data. It is necessary to have DNA sequence information from a multiplatform NGS with at least 2× and 30× depth sequence of genome coverage using Roche/454 and Illumina/Solexa, respectively, for effective an way of de novo assembly. Massive short-length reading data from the Illumina/Solexa system is enough to discover DNA variation, resulting in reducing the cost of DNA sequencing. Whole-genome expression profile data are useful to approach genome system biology with quantification of expressed RNAs from a whole-genome transcriptome, depending on the tissue samples. The hybrid mRNA sequences from Rohce/454 and Illumina/Solexa are more powerful to find novel genes through de novo assembly in any whole-genome sequenced species. The 20× and 50× coverage of the estimated transcriptome sequences using Roche/454 and Illumina/Solexa, respectively, is effective to create novel expressed reference sequences. However, only an average 30× coverage of a transcriptome with short read sequences of Illumina/Solexa is enough to check expression quantification, compared to the reference expressed sequence tag sequence.

5' end-centered expression profiling using cap-analysis gene expression and next-generation sequencing

Cap-analysis gene expression (CAGE) provides accurate high-throughput measurement of RNA expression. CAGE allows mapping of all the initiation sites of both capped coding and noncoding RNAs. In addition, transcriptional start sites within promoters are characterized at single-nucleotide resolution. The latter allows the regulatory inputs driving gene expression to be studied, which in turn enables the construction of transcriptional networks. Here we provide an optimized protocol for the construction of CAGE libraries on the basis of the preparation of 27-nt-long tags corresponding to initial bases at the 5' ends of capped RNAs. We have optimized the methods using simple steps based on filtration, which altogether takes 4 d to complete. The CAGE tags can be readily sequenced with Illumina sequencers, and upon modification they are also amenable to sequencing using other platforms.

Automated workflow for preparation of cDNA for cap analysis of gene expression on a single molecule sequencer

Background

Cap analysis of gene expression (CAGE) is a 5′ sequence tag technology to globally determine transcriptional starting sites in the genome and their expression levels and has most recently been adapted to the HeliScope single molecule sequencer. Despite significant simplifications in the CAGE protocol, it has until now been a labour intensive protocol.

Methodology

In this study we set out to adapt the protocol to a robotic workflow, which would increase throughput and reduce handling. The automated CAGE cDNA preparation system we present here can prepare 96 ‘HeliScope ready’ CAGE cDNA libraries in 8 days, as opposed to 6 weeks by a manual operator.We compare the results obtained using the same RNA in manual libraries and across multiple automation batches to assess reproducibility.

Conclusions

We show that the sequencing was highly reproducible and comparable to manual libraries with an 8 fold increase in productivity. The automated CAGE cDNA preparation system can prepare 96 CAGE sequencing samples simultaneously. Finally we discuss how the system could be used for CAGE on Illumina/SOLiD platforms, RNA-seq and full-length cDNA generation.

Sequencing technologies and genome sequencing

The high-throughput - next generation sequencing (HT-NGS) technologies are currently the hottest topic in the field of human and animals genomics researches, which can produce over 100 times more data compared to the most sophisticated capillary sequencers based on the Sanger method. With the ongoing developments of high throughput sequencing machines and advancement of modern bioinformatics tools at unprecedented pace, the target goal of sequencing individual genomes of living organism at a cost of $1,000 each is seemed to be realistically feasible in the near future. In the relatively short time frame since 2005, the HT-NGS technologies are revolutionizing the human and animal genome researches by analysis of chromatin immunoprecipitation coupled to DNA microarray (ChIP-chip) or sequencing (ChIP-seq), RNA sequencing (RNA-seq), whole genome genotyping, genome wide structural variation, de novo assembling and re-assembling of genome, mutation detection and carrier screening, detection of inherited disorders and complex human diseases, DNA library preparation, paired ends and genomic captures, sequencing of mitochondrial genome and personal genomics. In this review, we addressed the important features of HT-NGS like, first generation DNA sequencers, birth of HT-NGS, second generation HT-NGS platforms, third generation HT-NGS platforms: including single molecule Heliscope™, SMRT™ and RNAP sequencers, Nanopore, Archon Genomics X PRIZE foundation, comparison of second and third HT-NGS platforms, applications, advances and future perspectives of sequencing technologies on human and animal genome research.

Self-regulated alternative splicing at the AHNAK locus

AHNAK is a 700-kDa protein involved in cytoarchitecture and calcium signaling. It is secondarily reduced in muscle of dysferlinopathy patients and accumulates in muscle of calpainopathy patients, both affected by a muscular dystrophy. AHNAK directly interacts with dysferlin. This interaction is lost on cleavage of AHNAK by the protease calpain 3, explaining the molecular observations in patients. Currently, little is known of AHNAK regulation. We describe the self-regulation of multiple mRNA transcripts emanating from the AHNAK locus in muscle cells. We show that the AHNAK gene consists of a 17-kb exon flanked by multiple small exons. This genetic structure is shared by AHNAK2 and Periaxin, which share a common ancestor. Two major AHNAK transcripts are differentially expressed during muscle differentiation that encode for a small (17-kDa) and a large (700-kDa) protein isoform. These proteins interact in the cytoplasm, but the small AHNAK is also present in the nucleus. During muscle differentiation the small AHNAK is strongly increased, thereby establishing a positive feedback loop to regulate mRNA splicing of its own locus. A small 17-kDa isoform of Periaxin similarly traffics between the cytoplasm and the nucleus to regulate mRNA splicing. Thus, AHNAK constitutes a novel mechanism in post-transcriptional control of gene expression.

Differential gene expression in the siphonophore Nanomia bijuga (Cnidaria) assessed with multiple next-generation sequencing workflows

We investigated differential gene expression between functionally specialized feeding polyps and swimming medusae in the siphonophore Nanomia bijuga (Cnidaria) with a hybrid long-read/short-read sequencing strategy. We assembled a set of partial gene reference sequences from long-read data (Roche 454), and generated short-read sequences from replicated tissue samples that were mapped to the references to quantify expression. We collected and compared expression data with three short-read expression workflows that differ in sample preparation, sequencing technology, and mapping tools. These workflows were Illumina mRNA-Seq, which generates sequence reads from random locations along each transcript, and two tag-based approaches, SOLiD SAGE and Helicos DGE, which generate reads from particular tag sites. Differences in expression results across workflows were mostly due to the differential impact of missing data in the partial reference sequences. When all 454-derived gene reference sequences were considered, Illumina mRNA-Seq detected more than twice as many differentially expressed (DE) reference sequences as the tag-based workflows. This discrepancy was largely due to missing tag sites in the partial reference that led to false negatives in the tag-based workflows. When only the subset of reference sequences that unambiguously have tag sites was considered, we found broad congruence across workflows, and they all identified a similar set of DE sequences. Our results are promising in several regards for gene expression studies in non-model organisms. First, we demonstrate that a hybrid long-read/short-read sequencing strategy is an effective way to collect gene expression data when an annotated genome sequence is not available. Second, our replicated sampling indicates that expression profiles are highly consistent across field-collected animals in this case. Third, the impacts of partial reference sequences on the ability to detect DE can be mitigated through workflow choice and deeper reference sequencing.

LPAR1 and ITGA4 regulate peripheral blood monocyte counts

We recently mapped a quantitative trait locus for monocyte counts to chromosome 9q31 (rs7023923). Here we extend this work by showing with two independent approaches that rs7023923 regulates the expression levels of the nearby LPAR1 gene (P<0.0001), specifically implicating this gene in monocyte development. Furthermore, we tested 10 additional loci identified in the original analysis for replication in 1,122 individuals and confirm that rs6740847 near the alpha-4-integrin gene (ITGA4) associates with variation in monocyte counts (combined P=2.7×10(-10)). This variant is in complete linkage disequilibrium (r(2) =1) with a previously reported eQTL for ITGA4 (rs2124440), indicating that this is the likely causal gene in the region. Our results indicate that rs7023923 and rs6740847 respectively upregulate LPAR1 and downregulate ITGA4 expression and this increases the number of monocytes circulating in the peripheral blood. Further studies that investigate the downstream mechanism involved and the impact on immune function are warranted.

Unamplified cap analysis of gene expression on a single-molecule sequencer

We report the development of a simplified cap analysis of gene expression (CAGE) protocol adapted for single-molecule sequencers that avoids second strand synthesis, ligation, digestion, and PCR. HeliScopeCAGE directly sequences the 3' end of cap trapped first-strand cDNAs. As with previous versions of CAGE, we better define transcription start sites (TSS) than known models, identify novel regions of transcription and alternative promoters, and find two major classes of TSS signal, sharp peaks and broad regions. However, using this protocol, we observe reproducible evidence of regulation at the much finer level of individual TSS positions. The libraries are quantitative over 5 orders of magnitude and highly reproducible (Pearson's correlation coefficient of 0.987). We have also scaled down the sample requirement to 5 μg of total RNA for a standard HeliScopeCAGE library and 100 ng for a low-quantity version. When the same RNA was run as 5-μg and 100-ng versions, the 100 ng was still able to detect expression for ∼60% of the 13,468 loci detected by a 5-μg library using the same threshold, allowing comparative analysis of even rare cell populations. Testing the protocol for differential gene expression measurements on triplicate HeLa and THP-1 samples, we find that the log fold change compared to Illumina microarray measurements is highly correlated (0.871). In addition, HeliScopeCAGE finds differential expression for thousands more loci including those with probes on the array. Finally, although the majority of tags are 5' associated, we also observe a low level of signal on exons that is useful for defining gene structures.

mRNA degradation controls differentiation state-dependent differences in transcript and splice variant abundance

Expression profiling experiments usually provide a static snapshot of messenger RNA (mRNA) levels. Improved understanding of the dynamics of mRNA synthesis and degradation will aid the development of sound bioinformatic models for control of gene expression. We studied mRNA stability in proliferating and differentiated myogenic cells using whole-genome exon arrays and reported the decay rates (half life) for ∼7000 mRNAs. We showed that the stability of many mRNAs strongly depends on the differentiation status and contributes to differences in abundance of these mRNAs. In addition, alternative splicing turns out to be coupled to mRNA degradation. Although different splice forms may be produced at comparable levels, their relative abundance is partly determined by their different stabilities in proliferating and differentiated cells. Where the 3'-untranslated region (3'-UTR) was previously thought to contain most RNA stabilizing and destabilizing elements, we showed that this also holds for transcript isoforms sharing the same 3'-UTR. There are two splice variants in Itga7, of which the isoform with an extra internal exon is highly stable in differentiated cells but preferentially degraded in the cytoplasm of proliferating cells. In conclusion, control of stability and degradation emerge as important determinants for differential expression of mRNA transcripts and splice variants.