High-throughput and quantitative genome-wide messenger RNA sequencing for molecular phenotyping

Affordable, accurate and unbiased RNA sequencing by manual library miniaturization: A case study in barley

We present an easy-to-reproduce manual miniaturized full-length RNA sequencing (RNAseq) library preparation workflow that does not require the upfront investment in expensive lab equipment or long setup times. With minimal adjustments to an established commercial protocol, we were able to manually miniaturize the RNAseq library preparation by a factor of up to 1:8. This led to cost savings for miniaturized library preparation of up to 86.1% compared to the gold standard. The resulting data were the basis of a rigorous quality control analysis that inspected: sequencing quality metrics, gene body coverage, raw read duplications, alignment statistics, read pair duplications, detected transcripts and sequence variants. We also included a deep dive data analysis identifying rRNA contamination and suggested ways to circumvent these. In the end, we could not find any indication of biases or inaccuracies caused by the RNAseq library miniaturization. The variance in detected transcripts was minimal and not influenced by the miniaturization level. Our results suggest that the workflow is highly reproducible and the sequence data suitable for downstream analyses such as differential gene expression analysis or variant calling.

Effects of Β‒Blocker Administration on Cardiac Function: A Coronary Computed Tomography Angiography Study

BACKGROUND

β-blockers are widely used for lowering heart rate (HR) during coronary computed tomography angiography (CCTA); however, they should be used with caution for patients with heart failure as they may have a negative inotropic effect.

OBJECTIVE

To clarify the effects of β-blockers (oral and intravenous injection) on cardiac function using CCTA.

METHODS

A total of 244 patients (men: women = 166: 78; mean age, 64.4 years old) suspected of having ischemic cardiac disease and had undergone echocardiography within 3 months before and after CCTA were included in the study. Systematic errors in ejection fraction (EF) were corrected by calculating ΔEF from the EF difference between echocardiography and CCTA in patients not using β- blockers. Univariate and multivariate analyses were performed for factors affecting ΔEF. In addition, HR between, before, and during CCTA were compared by Wilcoxon's test.

RESULTS

Temporary oral or intravenous administration of β-blockers at the CCTA had no significant effects on EF (p = 0.70), whereas HR was significantly decreased (p < 0.001). However, regular administration of β-blockers increases the EF on CCTA.

CONCLUSION

The administration of β-blockers immediately before CCTA affects HR but not EF. Premedication with β-blockers can be safely used for patients who undergo CCTA, and CCTA is useful for EF evaluation, independent of the use of β-blockers.

Loss of slc39a14 causes simultaneous manganese hypersensitivity and deficiency in zebrafish

Manganese neurotoxicity is a hallmark of hypermanganesemia with dystonia 2, an inherited manganese transporter defect caused by mutations in SLC39A14. To identify novel potential targets of manganese neurotoxicity, we performed transcriptome analysis of slc39a14-/- mutant zebrafish that were exposed to MnCl2. Differentially expressed genes mapped to the central nervous system and eye, and pathway analysis suggested that Ca2+ dyshomeostasis and activation of the unfolded protein response are key features of manganese neurotoxicity. Consistent with this interpretation, MnCl2 exposure led to decreased whole-animal Ca2+ levels, locomotor defects and changes in neuronal activity within the telencephalon and optic tectum. In accordance with reduced tectal activity, slc39a14-/- zebrafish showed changes in visual phototransduction gene expression, absence of visual background adaptation and a diminished optokinetic reflex. Finally, numerous differentially expressed genes in mutant larvae normalised upon MnCl2 treatment indicating that, in addition to neurotoxicity, manganese deficiency is present either subcellularly or in specific cells or tissues. Overall, we assembled a comprehensive set of genes that mediate manganese-systemic responses and found a highly correlated and modulated network associated with Ca2+ dyshomeostasis and cellular stress. This article has an associated First Person interview with the first author of the paper.

Allele-specific gene expression can underlie altered transcript abundance in zebrafish mutants

In model organisms, RNA-sequencing (RNA-seq) is frequently used to assess the effect of genetic mutations on cellular and developmental processes. Typically, animals heterozygous for a mutation are crossed to produce offspring with different genotypes. Resultant embryos are grouped by genotype to compare homozygous mutant embryos to heterozygous and wild-type siblings. Genes that are differentially expressed between the groups are assumed to reveal insights into the pathways affected by the mutation. Here we show that in zebrafish, differentially expressed genes are often over-represented on the same chromosome as the mutation due to different levels of expression of alleles from different genetic backgrounds. Using an incross of haplotype-resolved wild-type fish, we found evidence of widespread allele-specific expression, which appears as differential expression when comparing embryos homozygous for a region of the genome to their siblings. When analysing mutant transcriptomes, this means that the differential expression of genes on the same chromosome as a mutation of interest may not be caused by that mutation. Typically, the genomic location of a differentially expressed gene is not considered when interpreting its importance with respect to the phenotype. This could lead to pathways being erroneously implicated or overlooked due to the noise of spurious differentially expressed genes on the same chromosome as the mutation. These observations have implications for the interpretation of RNA-seq experiments involving outbred animals and non-inbred model organisms.

Integrated analysis highlights multiple long non‑coding RNAs and their potential roles in the progression of human esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a prevalent aggressive malignant tumor with poor prognosis. Investigations into the molecular changes that occur as a result of the disease, as well as identification of novel biomarkers for its diagnosis and prognosis, are urgently required. Long non‑coding RNAs (lncRNAs) have been reported to play a critical role in tumor progression. The present study performed data mining analyses for ESCC via an integrated study of accumulated datasets and identification of the differentially expressed lncRNAs from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. The identified intersection of differentially expressed genes (lncRNAs, miRNAs and mRNAs) in ESCC tissues between the GEO and TCGA datasets was investigated. Based on these intersected lncRNAs, the present study constructed a competitive endogenous RNA (ceRNA) network of lncRNAs in ESCC. A total of 81 intersection lncRNAs were identified; 67 of these were included in the ceRNA network. Functional analyses revealed that these 67 key lncRNAs primarily dominated cellular biological processes. The present study then analyzed the associations between the expression levels of these 67 key lncRNAs and the clinicopathological characteristics of the ESCC patients, as well as their survival time using TCGA. The results revealed that 31 of these lncRNAs were associated with tumor grade, tumor‑node‑metastasis (TNM) stage and lymphatic metastasis status (P<0.05). In addition, 15 key lncRNAs were demonstrated to be associated with survival time (P<0.05). Finally, 5 key lncRNAs were selected for validation of their expression levels in 30 patients newly diagnosed with ESCC via reverse transcription‑quantitative PCR (RT‑qPCR). The results suggested that the fold changes in the trends of up‑ and downregulation between GEO, TCGA and RT‑qPCR were consistent. In addition, it was also demonstrated that a select few of these 5 key lncRNAs were significantly associated with TNM stage and lymph node metastasis (P<0.05). The results of the clinically relevant analysis and the aforementioned bioinformatics were similar, hence proving that the bioinformatics analysis used in the present study is credible. Overall, the results from the present study may provide further insight into the functional characteristics of lncRNAs in ESCC through bioinformatics integrative analysis of the GEO and TCGA datasets, and reveal potential diagnostic and prognostic biomarkers for ESCC.

The gene regulatory basis of genetic compensation during neural crest induction

The neural crest (NC) is a vertebrate-specific cell type that contributes to a wide range of different tissues across all three germ layers. The gene regulatory network (GRN) responsible for the formation of neural crest is conserved across vertebrates. Central to the induction of the NC GRN are AP-2 and SoxE transcription factors. NC induction robustness is ensured through the ability of some of these transcription factors to compensate loss of function of gene family members. However the gene regulatory events underlying compensation are poorly understood. We have used gene knockout and RNA sequencing strategies to dissect NC induction and compensation in zebrafish. We genetically ablate the NC using double mutants of tfap2a;tfap2c or remove specific subsets of the NC with sox10 and mitfa knockouts and characterise genome-wide gene expression levels across multiple time points. We find that compensation through a single wild-type allele of tfap2c is capable of maintaining early NC induction and differentiation in the absence of tfap2a function, but many target genes have abnormal expression levels and therefore show sensitivity to the reduced tfap2 dosage. This separation of morphological and molecular phenotypes identifies a core set of genes required for early NC development. We also identify the 15 somites stage as the peak of the molecular phenotype which strongly diminishes at 24 hpf even as the morphological phenotype becomes more apparent. Using gene knockouts, we associate previously uncharacterised genes with pigment cell development and establish a role for maternal Hippo signalling in melanocyte differentiation. This work extends and refines the NC GRN while also uncovering the transcriptional basis of genetic compensation via paralogues.

Embryonic development is an intricate process that requires genes to be active at the right time and place. Organisms have evolved mechanisms that ensure faithful execution of developmental programmes even if genes fail to function. For example, in a process called genetic compensation, one or more genes become activated in response to loss of function of another. In this work we use the zebrafish model to investigate how two related genes, tfap2a and tfap2c, interact to ensure establishment of the neural crest, a vertebrate-specific cell type that contributes to many different tissues. Losing tfap2a activity causes mild morphological defects and losing tfap2c has no visible effect. Yet when both are inactive, embryos are severely abnormal due to lack of neural crest-derived tissues. Here we show that loss of tfap2a triggers upregulation of tfap2c which prevents the loss of neural crest tissue. However, the genes normally regulated by tfap2a respond differently to tfap2c allowing us to identify the first tier of the Ap2 network and new players in neural crest biology. Our work demonstrates that the expression signature of partial, but morphologically sufficient, genetic compensation provides an opportunity to dissect gene regulatory networks.

Fishing forward and reverse: Advances in zebrafish phenomics

Understanding how the genome instructs the phenotypic characteristics of an organism is one of the major scientific endeavors of our time. Advances in genetics have progressively deciphered the inheritance, identity and biological relevance of genetically encoded information, contributing to the rise of several, complementary omic disciplines. One of them is phenomics, an emergent area of biology dedicated to the systematic multi-scale analysis of phenotypic traits. This discipline provides valuable gene function information to the rapidly evolving field of genetics. Current molecular tools enable genome-wide analyses that link gene sequence to function in multi-cellular organisms, illuminating the genome-phenome relationship. Among vertebrates, zebrafish has emerged as an outstanding model organism for high-throughput phenotyping and modeling of human disorders. Advances in both systematic mutagenesis and phenotypic analyses of embryonic and post-embryonic stages in zebrafish have revealed the function of a valuable collection of genes and the general structure of several complex traits. In this review, we summarize multiple large-scale genetic efforts addressing parental, embryonic, and adult phenotyping in the zebrafish. The genetic and quantitative tools available in the zebrafish model, coupled with the broad spectrum of phenotypes that can be assayed, make it a powerful model for phenomics, well suited for the dissection of genotype-phenotype associations in development, physiology, health and disease.

Elimination of PCR duplicates in RNA-seq and small RNA-seq using unique molecular identifiers

Background

RNA-seq and small RNA-seq are powerful, quantitative tools to study gene regulation and function. Common high-throughput sequencing methods rely on polymerase chain reaction (PCR) to expand the starting material, but not every molecule amplifies equally, causing some to be overrepresented. Unique molecular identifiers (UMIs) can be used to distinguish undesirable PCR duplicates derived from a single molecule and identical but biologically meaningful reads from different molecules.

Results

We have incorporated UMIs into RNA-seq and small RNA-seq protocols and developed tools to analyze the resulting data. Our UMIs contain stretches of random nucleotides whose lengths sufficiently capture diverse molecule species in both RNA-seq and small RNA-seq libraries generated from mouse testis. Our approach yields high-quality data while allowing unique tagging of all molecules in high-depth libraries.

Conclusions

Using simulated and real datasets, we demonstrate that our methods increase the reproducibility of RNA-seq and small RNA-seq data. Notably, we find that the amount of starting material and sequencing depth, but not the number of PCR cycles, determine PCR duplicate frequency. Finally, we show that computational removal of PCR duplicates based only on their mapping coordinates introduces substantial bias into data analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4933-1) contains supplementary material, which is available to authorized users.

A high-resolution mRNA expression time course of embryonic development in zebrafish

We have produced an mRNA expression time course of zebrafish development across 18 time points from 1 cell to 5 days post-fertilisation sampling individual and pools of embryos. Using poly(A) pulldown stranded RNA-seq and a 3′ end transcript counting method we characterise temporal expression profiles of 23,642 genes. We identify temporal and functional transcript co-variance that associates 5024 unnamed genes with distinct developmental time points. Specifically, a class of over 100 previously uncharacterised zinc finger domain containing genes, located on the long arm of chromosome 4, is expressed in a sharp peak during zygotic genome activation. In addition, the data reveal new genes and transcripts, differential use of exons and previously unidentified 3′ ends across development, new primary microRNAs and temporal divergence of gene paralogues generated in the teleost genome duplication. To make this dataset a useful baseline reference, the data can be browsed and downloaded at Expression Atlas and Ensembl.

UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy

Unique Molecular Identifiers (UMIs) are random oligonucleotide barcodes that are increasingly used in high-throughput sequencing experiments. Through a UMI, identical copies arising from distinct molecules can be distinguished from those arising through PCR amplification of the same molecule. However, bioinformatic methods to leverage the information from UMIs have yet to be formalized. In particular, sequencing errors in the UMI sequence are often ignored or else resolved in an ad hoc manner. We show that errors in the UMI sequence are common and introduce network-based methods to account for these errors when identifying PCR duplicates. Using these methods, we demonstrate improved quantification accuracy both under simulated conditions and real iCLIP and single-cell RNA-seq data sets. Reproducibility between iCLIP replicates and single-cell RNA-seq clustering are both improved using our proposed network-based method, demonstrating the value of properly accounting for errors in UMIs. These methods are implemented in the open source UMI-tools software package.

Global open data management in metabolomics

•

The metabolome allows accessing the external influences under which an organism exists and develops in a dynamic way.

•

Recent years have seen the establishment of a global network for metabolomics data exchange.

•

Global metabolomics data exchange is leading to an exponential growth of publically available metabolomics data for re-analysis.

Chemical Biology employs chemical synthesis, analytical chemistry and other tools to study biological systems. Recent advances in both molecular biology such as next generation sequencing (NGS) have led to unprecedented insights towards the evolution of organisms’ biochemical repertoires. Because of the specific data sharing culture in Genomics, genomes from all kingdoms of life become readily available for further analysis by other researchers. While the genome expresses the potential of an organism to adapt to external influences, the Metabolome presents a molecular phenotype that allows us to asses the external influences under which an organism exists and develops in a dynamic way. Steady advancements in instrumentation towards high-throughput and highresolution methods have led to a revival of analytical chemistry methods for the measurement and analysis of the metabolome of organisms. This steady growth of metabolomics as a field is leading to a similar accumulation of big data across laboratories worldwide as can be observed in all of the other omics areas. This calls for the development of methods and technologies for handling and dealing with such large datasets, for efficiently distributing them and for enabling re-analysis. Here we describe the recently emerging ecosystem of global open-access databases and data exchange efforts between them, as well as the foundations and obstacles that enable or prevent the data sharing and reanalysis of this data.

DANIO-CODE: Toward an Encyclopedia of DNA Elements in Zebrafish

The zebrafish has emerged as a model organism for genomics studies. The symposium “Toward an encyclopedia of DNA elements in zebrafish” held in London in December 2014, was coorganized by Ferenc Müller and Fiona Wardle. This meeting is a follow-up of a similar previous workshop held 2 years earlier and represents a push toward the formalization of a community effort to annotate functional elements in the zebrafish genome. The meeting brought together zebrafish researchers, bioinformaticians, as well as members of established consortia, to exchange scientific findings and experience, as well as to discuss the initial steps toward the formation of a DANIO-CODE consortium. In this study, we provide the latest updates on the current progress of the consortium's efforts, opening up a broad invitation to researchers to join in and contribute to DANIO-CODE.

The Ribosome Biogenesis Protein Nol9 Is Essential for Definitive Hematopoiesis and Pancreas Morphogenesis in Zebrafish

Ribosome biogenesis is a ubiquitous and essential process in cells. Defects in ribosome biogenesis and function result in a group of human disorders, collectively known as ribosomopathies. In this study, we describe a zebrafish mutant with a loss-of-function mutation in nol9, a gene that encodes a non-ribosomal protein involved in rRNA processing. nol9sa1022/sa1022 mutants have a defect in 28S rRNA processing. The nol9sa1022/sa1022 larvae display hypoplastic pancreas, liver and intestine and have decreased numbers of hematopoietic stem and progenitor cells (HSPCs), as well as definitive erythrocytes and lymphocytes. In addition, ultrastructural analysis revealed signs of pathological processes occurring in endothelial cells of the caudal vein, emphasizing the complexity of the phenotype observed in nol9sa1022/sa1022 larvae. We further show that both the pancreatic and hematopoietic deficiencies in nol9sa1022/sa1022 embryos were due to impaired cell proliferation of respective progenitor cells. Interestingly, genetic loss of Tp53 rescued the HSPCs but not the pancreatic defects. In contrast, activation of mRNA translation via the mTOR pathway by L-Leucine treatment did not revert the erythroid or pancreatic defects. Together, we present the nol9sa1022/sa1022 mutant, a novel zebrafish ribosomopathy model, which recapitulates key human disease characteristics. The use of this genetically tractable model will enhance our understanding of the tissue-specific mechanisms following impaired ribosome biogenesis in the context of an intact vertebrate.