Deep Transcriptome Sequencing Research Articles

Prolyl 3-Hydroxylase 2 Is a Molecular Player of Angiogenesis.

Prolyl 3-hydroxylase 2 (P3H2) catalyzes the post-translational formation of 3-hydroxyproline on collagens, mainly on type IV. Its activity has never been directly associated to angiogenesis. Here, we identified P3H2 gene through a deep-sequencing transcriptome analysis of human umbilical vein endothelial cells (HUVECs) stimulated with vascular endothelial growth factor A (VEGF-A). Differently from many previous studies we carried out the stimulation not on starved HUVECs, but on cells grown to maintain the best condition for their in vitro survival and propagation. We showed that P3H2 is induced by VEGF-A in two primary human endothelial cell lines and that its transcription is modulated by VEGF-A/VEGF receptor 2 (VEGFR-2) signaling pathway through p38 mitogen-activated protein kinase (MAPK). Then, we demonstrated that P3H2, through its activity on type IV Collagen, is essential for angiogenesis properties of endothelial cells in vitro by performing experiments of gain- and loss-of-function. Immunofluorescence studies showed that the overexpression of P3H2 induced a more condensed status of Collagen IV, accompanied by an alignment of the cells along the Collagen IV bundles, so towards an evident pro-angiogenic status. Finally, we found that P3H2 knockdown prevents pathological angiogenesis in vivo, in the model of laser-induced choroid neovascularization. Together these findings reveal that P3H2 is a new molecular player involved in new vessels formation and could be considered as a potential target for anti-angiogenesis therapy.

MP3RNA-seq: Massively parallel 3' end RNA sequencing for high-throughput gene expression profiling and genotyping.

Transcriptome deep sequencing (RNA-seq) has become a routine method for global gene expression profiling. However, its application to large-scale experiments remains limited by cost and labor constraints. Here we describe a massively parallel 3' end RNA-seq (MP3RNA-seq) method that introduces unique sample barcodes during reverse transcription to permit sample pooling immediately following this initial step. MP3RNA-seq allows for handling of hundreds of samples in a single experiment, at a cost of about $6 per sample for library construction and sequencing. MP3RNA-seq is effective for not only high-throughput gene expression profiling, but also genotyping. To demonstrate its utility, we applied MP3RNA-seq to 477 double haploid lines of maize. We identified 19,429 genes expressed in at least 50% of the lines and 35,836 high-quality single nucleotide polymorphisms for genotyping analysis. Armed with these data, we performed expression and agronomic trait quantitative trait locus (QTL) mapping and identified 25,797 expression QTLs for 15,335 genes and 21 QTLs for plant height, ear height, and relative ear height. We conclude that MP3RNA-seq is highly reproducible, accurate, and sensitive for high-throughput gene expression profiling and genotyping, and should be generally applicable to most eukaryotic species.

Functional long non-coding and circular RNAs in zebrafish.

The utility of model organisms to understand the function of a novel transcript/genes has allowed us to delineate their molecular mechanisms in maintaining cellular homeostasis. Organisms such as zebrafish have contributed a lot in the field of developmental and disease biology. Attributable to advancement and deep transcriptomics, many new transcript isoforms and non-coding RNAs such as long noncoding RNA (lncRNA) and circular RNAs (circRNAs) have been identified and cataloged in multiple databases and many more are yet to be identified. Various methods and tools have been utilized to identify lncRNAs/circRNAs in zebrafish using deep sequencing of transcriptomes as templates. Functional analysis of a few candidates such as tie1-AS, ECAL1 and CDR1as in zebrafish provides a prospective outline to approach other known or novel lncRNA/circRNA. New genetic alteration tools like TALENS and CRISPRs have helped in probing for the molecular function of lncRNA/circRNA in zebrafish. Further latest improvements in experimental and computational techniques offer the identification of lncRNA/circRNA counterparts in humans and zebrafish thereby allowing easy modeling and analysis of function at cellular level.

Deep sequencing of the Sanghuangporus vaninii transcriptome reveals dynamic landscapes of candidate genes involved in the biosynthesis of active compounds

The medicinal fungus Sanghuang produces diverse bioactive compounds and is widely used in Asian countries. However, little is known about the genes and pathways involved in the biosynthesis of these active compounds. Based on our previous study providing Sanghuangporus vaninii genomic information, the transcriptomes of MY (mycelium), OY (1-year-old fruiting bodies), and TY (3-year-old fruiting bodies) were determined in this study. A significant number of genes (4774) were up- or downregulated between mycelium and fruiting bodies, but only 1422 differentially expressed genes were detected between OY and TY. 138 genes encoding P450s were identified in the fungal genome and grouped into 25 P450 families; more than 64% (88) of the genes were significantly differentially expressed between the mycelium and fruiting body, suggesting that these P450s are involved in fungal sexual development. Importantly, the expression of genes involved in bioactive compound (triterpenoids, polysaccharides, and flavonoids) biosynthesis in asexual (cultured with solid and liquid media) and sexual stages was explored and combined with transcriptome and quantitative PCR analyses. More genes involved in the biosynthesis of bioactive compounds were expressed more highly in mycelium than in fruiting bodies under liquid medium culture compared with solid medium culture, which was consistent with the yields of different bioactive compounds, suggesting that liquid fermentation of S. vaninii Kangneng can be used to obtain these bioactive compounds. A comprehensive understanding of the genomic information of S. vaninii will facilitate its potential use in pharmacological and industrial applications.

Pluripotency and immunomodulatory signatures of canine induced pluripotent stem cell-derived mesenchymal stromal cells are similar to harvested mesenchymal stromal cells

With a view towards harnessing the therapeutic potential of canine mesenchymal stromal cells (cMSCs) as modulators of inflammation and the immune response, and to avoid the issues of the variable quality and quantity of harvested cMSCs, we examined the immunomodulatory properties of cMSCs derived from canine induced pluripotent stem cells (ciMSCs), and compared them to cMSCs harvested from adipose tissue (cAT-MSC) and bone marrow (cBM-MSC). A combination of deep sequencing and quantitative RT-PCR of the ciMSC transcriptome confirmed that ciMSCs express more genes in common with cBM-MSCs and cAT-MSCs than with the ciPSCs from which they were derived. Both ciMSCs and harvested cMSCs express a range of pluripotency factors in common with the ciPSCs including NANOG, POU5F1 (OCT-4), SOX-2, KLF-4, LIN-28A, MYC, LIF, LIFR, and TERT. However, ESRRB and PRDM-14, both factors associated with naïve, rather than primed, pluripotency were expressed only in the ciPSCs. CXCR-4, which is essential for the homing of MSCs to sites of inflammation, is also detectable in ciMSCs, cAT- and cBM-MSCs, but not ciPSCs. ciMSCs constitutively express the immunomodulatory factors iNOS, GAL-9, TGF-β1, PTGER-2α and VEGF, and the pro-inflammatory mediators COX-2, IL-1β and IL-8. When stimulated with the canine pro-inflammatory cytokines tumor necrosis factor-α (cTNF-α), interferon-γ (cIFN-γ), or a combination of both, ciMSCs upregulated their expression of IDO, iNOS, GAL-9, HGF, TGF-β1, PTGER-2α, VEGF, COX-2, IL-1β and IL-8. When co-cultured with mitogen-stimulated lymphocytes, ciMSCs downregulated their expression of iNOS, HGF, TGF-β1 and PTGER-2α, while increasing their expression of COX-2, IDO and IL-1β. Taken together, these findings suggest that ciMSCs possess similar immunomodulatory capabilities as harvested cMSCs and support further investigation into their potential use for the management of canine immune-mediated and inflammatory disorders.

Deep transcriptome sequencing of subgenual anterior cingulate cortex reveals cross-diagnostic and diagnosis-specific RNA expression changes in major psychiatric disorders

Despite strong evidence of heritability and growing discovery of genetic markers for major mental illness, little is known about how gene expression in the brain differs across psychiatric diagnoses, or how known genetic risk factors shape these differences. Here we investigate expressed genes and gene transcripts in postmortem subgenual anterior cingulate cortex (sgACC), a key component of limbic circuits linked to mental illness. RNA obtained postmortem from 200 donors diagnosed with bipolar disorder, schizophrenia, major depression, or no psychiatric disorder was deeply sequenced to quantify expression of over 85,000 gene transcripts, many of which were rare. Case–control comparisons detected modest expression differences that were correlated across disorders. Case–case comparisons revealed greater expression differences, with some transcripts showing opposing patterns of expression between diagnostic groups, relative to controls. The ~250 rare transcripts that were differentially-expressed in one or more disorder groups were enriched for genes involved in synapse formation, cell junctions, and heterotrimeric G-protein complexes. Common genetic variants were associated with transcript expression (eQTL) or relative abundance of alternatively spliced transcripts (sQTL). Common genetic variants previously associated with disease risk were especially enriched for sQTLs, which together accounted for disproportionate fractions of diagnosis-specific heritability. Genetic risk factors that shape the brain transcriptome may contribute to diagnostic differences between broad classes of mental illness.

The Characteristics of Chemosensory and Opsin Genes in Newly Emerged and Sexually Mature Agrilus planipennis, an Important Quarantine Forest Beetle.

The emerald ash borer (EAB), Agrilus planipennis, is a highly destructive quarantine pest. The olfactory and visual systems of A. planipennis play different but critical roles at newly emerged and sexually mature stages; however, the molecular basis underlying these differences remain unclear. Consequently, based on deep transcriptome sequencing, we evaluated the expression levels of chemosensory-related proteins and opsins at the two developmental stages of A. planipennis. We found 15 new chemosensory-related genes in our transcriptome assembly compared with the previous genome assembly, including 6 that code for odorant-binding proteins (OBPs) and 9 for chemosensory proteins (CSPs). The expression of several chemosensory-related genes (OBP7, OBP10, CSP1, and CSP12) differed markedly between newly emerged and sexually mature A. planipennis. We also found that the expression of UV opsin 2 and LW opsin 1 was higher in sexually mature male A. planipennis, which may be associated with their strong visual mate detection ability. This study forms the basis for further investigation of the chemosensory and visual system of A. planipennis, and these differentially expressed genes between newly emerged and sexually mature stages may serve as targets for the management of this destructive forest pest after sexual maturity.

WheatGmap: a comprehensive platform for wheat gene mapping and genomic studies

WheatGmap: a comprehensive platform for wheat gene mapping and genomic studies

Dual Oxidase Mutant Retards Mauthner-Cell Axon Regeneration at an Early Stage via Modulating Mitochondrial Dynamics in Zebrafish.

Dual oxidase (duox)-derived reactive oxygen species (ROS) have been correlated with neuronal polarity, cerebellar development, and neuroplasticity. However, there have not been many comprehensive studies of the effect of individual duox isoforms on central-axon regeneration in vivo. Here, we explored this question in zebrafish, an excellent model organism for central-axon regeneration studies. In our research, mutation of the duox gene with CRISPR/Cas9 significantly retarded the single-axon regeneration of the zebrafish Mauthner cell in vivo. Using deep transcriptome sequencing, we found that the expression levels of related functional enzymes in mitochondria were down-regulated in duox mutant fish. In vivo imaging showed that duox mutants had significantly disrupted mitochondrial transport and redox state in single Mauthner-cell axon. Our research data provide insights into how duox is involved in central-axon regeneration by changing mitochondrial transport.

De novo transcriptome sequencing of Rhododendron molle and identification of genes involved in the biosynthesis of secondary metabolites

BackgroundRhododendron molle (Ericaceae) is a traditional Chinese medicinal plant, its flower and root have been widely used to treat rheumatism and relieve pain for thousands of years in China. Chemical studies have revealed that R. molle contains abundant secondary metabolites such as terpenoinds, flavonoids and lignans, some of which have exhibited various bioactivities including antioxidant, hypotension and analgesic activity. In spite of immense pharmaceutical importance, the mechanism underlying the biosynthesis of secondary metabolites remains unknown and the genomic information is unavailable.ResultsTo gain molecular insight into this plant, especially on the information of pharmaceutically important secondary metabolites including grayanane diterpenoids, we conducted deep transcriptome sequencing for R. molle flower and root using the Illumina Hiseq platform. In total, 100,603 unigenes were generated through de novo assembly with mean length of 778 bp, 57.1% of these unigenes were annotated in public databases and 17,906 of those unigenes showed significant match in the KEGG database. Unigenes involved in the biosynthesis of secondary metabolites were annotated, including the TPSs and CYPs that were potentially responsible for the biosynthesis of grayanoids. Moreover, 3376 transcription factors and 10,828 simple sequence repeats (SSRs) were also identified. Additionally, we further performed differential gene expression (DEG) analysis of the flower and root transcriptome libraries and identified numerous genes that were specifically expressed or up-regulated in flower.ConclusionsTo the best of our knowledge, this is the first time to generate and thoroughly analyze the transcriptome data of both R. molle flower and root. This study provided an important genetic resource which will shed light on elucidating various secondary metabolite biosynthetic pathways in R. molle, especially for those with medicinal value and allow for drug development in this plant.

Alternative splicing in plant abiotic stress responses

Modifications of the cellular proteome pool upon stress allow plants to tolerate environmental changes. Alternative splicing is the most significant mechanism responsible for the production of multiple protein isoforms from a single gene. The spliceosome, a large ribonucleoprotein complex, together with several associated proteins, controls this pre-mRNA processing, adding an additional level of regulation to gene expression. Deep sequencing of transcriptomes revealed that this co- or post-transcriptional mechanism is highly induced by abiotic stress, and concerns vast numbers of stress-related genes. Confirming the importance of splicing in plant stress adaptation, key players of stress signaling have been shown to encode alternative transcripts, whereas mutants lacking splicing factors or associated components show a modified sensitivity and defective responses to abiotic stress. Here, we examine recent literature on alternative splicing and splicing alterations in response to environmental stresses, focusing on its role in stress adaptation and analyzing the future perspectives and directions for research.

High-Throughput Sequencing to Detect DNA-RNA Changes.

The advent of deep sequencing technologies has greatly improved the study of complex eukaryotic genomes and transcriptomes, allowing the investigation of posttranscriptional molecular mechanisms as alternative splicing and RNA editing at unprecedented throughput and resolution. The most prevalent type of RNA editing in higher eukaryotes is the deamination of adenosine to inosine (A-to-I) in double-stranded RNAs. Depending on the RNA type or the RNA region involved, A-to-I RNA editing contributes to the transcriptome and proteome diversity.Hereafter, we present an easy and reproducible computational protocol for the identification of candidate RNA editing sites in humans using deep transcriptome (RNA-Seq) and genome (DNA-Seq) sequencing.

Deep sequencing of liver explant transcriptomes reveals extensive expression from integrated hepatitis B virus DNA.

Hepatitis B virus (HBV) is a major cause of hepatocellular carcinoma (HCC). Integration of HBV DNA into the human genome may contribute to oncogenesis and to the production of the hepatitis B surface antigen (HBsAg). Whether integrations contribute to HBsAg levels in the blood is poorly known. Here, we characterize the HBV RNA profile of HBV integrations in liver tissue in patients with chronic HBV infection, with or without concurrent hepatitis D infection, by transcriptome deep sequencing. Transcriptomes were determined in liver tissue by deep sequencing providing 200 million reads per sample. Integration points were identified using a bioinformatic pipeline. Explanted liver tissue from five patients with end-stage liver disease caused by HBV or HBV/HDV was studied along with publicly available transcriptomes from 21 patients. Almost all HBV RNA profiles were devoid of reads in the core and the 3' redundancy (nt 1830-1927) regions, and contained a large number of chimeric viral/human reads. Hence, HBV transcripts from integrated HBV DNA rather than from covalently closed circular HBV DNA (cccDNA) predominated in late-stage HBV infection, in particular in cases with hepatitis D virus co-infection. The findings support the suggestion that integrated HBV DNA can be a significant source of HBsAg in humans.

The DNA methylation landscape of advanced prostate cancer.

Although DNA methylation is a key regulator of gene expression, the comprehensive methylation landscape of metastatic cancer has never been defined. Through whole-genome bisulfite sequencing paired with deep whole-genome and transcriptome sequencing of 100 castration-resistant prostate metastases, we discovered alterations affecting driver genes only detectable with integrated whole-genome approaches. Notably, we observed that 22% of tumors exhibited a novel epigenomic subtype associated with hyper-methylation and somatic mutations in TET2, DNMT3B, IDH1, and BRAF. We also identified intergenic regions where methylation is associated with RNA expression of the oncogenic driver genes AR, MYC and ERG. Finally, we showed that differential methylation during progression preferentially occurs at somatic mutational hotspots and putative regulatory regions. This study is a large integrated study of whole-genome, whole-methylome and whole-transcriptome sequencing in metastatic cancer and provides a comprehensive overview of the important regulatory role of methylation in metastatic castration-resistant prostate cancer.

The comprehensive methylation landscape of metastatic castration-resistant prostate cancer (mCRPC) identifies new phenotypic subtypes: Results from the West Coast Prostate Cancer Dream Team (WCDT).

5507 Background: While recent studies have delineated the genomic landscape of mCRPC, its epigenomic landscape has not been as well characterized. The goal of this study was to define the comprehensive methylation landscape of mCRPC. Methods: mCRPC patients (pts) underwent a metastasis biopsy as part of a multi-institutional study (NCT02432001). Deep whole-genome bisulfite sequencing (mean depth 46x) was performed on fresh frozen tissue from 100 mCRPC patients; data was paired with deep whole-genome and transcriptome sequencing from the same samples. Unbiased hierarchical clustering of the mCRPC methylome was undertaken, and the survival of patients in each cluster was calculated using the Kaplan Meier method. Results: Unbiased hierarchical clustering revealed several distinct subtypes. 22% of mCRPC samples exhibited a novel epigenomic subtype associated with hyper-methylation. This hypermethylated (HM) cluster was significantly associated with somatic mutations in genes known to be involved in methylation, eg TET2 and DNMT3B, as well as in genes in which mutations have been associated with hyper-methylation in other cancer types ( IDH1 in glioblastoma and BRAF in colon cancer). mCRPC survival was 56.1 mos in pts with HM cancers compared to 35.6 mos in non-HM (p = .055). Methylome clustering also identified a unique cluster comprised of all patients with treatment-induced small cell/neuroendocrine cancer, a subtype previously associated with poor survival. Conclusions: This integrated study of whole-genome, whole methylome and whole-transcriptome sequencing provides the first comprehensive overview of the important regulatory role of methylation in metastatic castration-resistant prostate cancer, and has identified at least two distinct subtypes. The clinical and therapeutic implications of methylation subtypes should be explored in future studies. Clinical trial information: NCT02432001 .

Global transcriptome analysis reveals circadian control of splicing events in Arabidopsis thaliana.

The circadian clock of Arabidopsis thaliana controls many physiological and molecular processes, allowing plants to anticipate daily changes in their environment. However, developing a detailed understanding of how oscillations in mRNA levels are connected to oscillations in co/post-transcriptional processes, such as splicing, has remained a challenge. Here we applied a combined approach using deep transcriptome sequencing and bioinformatics tools to identify novel circadian-regulated genes and splicing events. Using a stringent approach, we identified 300 intron retention, eight exon skipping, 79 alternative 3' splice site usage, 48 alternative 5' splice site usage, and 350 multiple (more than one event type) annotated events under circadian regulation. We also found seven and 721 novel alternative exonic and intronic events. Depletion of the circadian-regulated splicing factor AtSPF30 homologue resulted in the disruption of a subset of clock-controlled splicing events. Altogether, our global circadian RNA-seq coupled with an in silico, event-centred, splicing analysis tool offers a new approach for studying the interplay between the circadian clock and the splicing machinery at a global scale. The identification of many circadian-regulated splicing events broadens our current understanding of the level of control that the circadian clock has over this co/post-transcriptional regulatory layer.

Candidate resistance genes selection and transcriptome analysis for the early responses to Plasmopara viticola infection in grape cultivars

Grape downy mildew (Plasmopara viticola) is one of the most destructive diseases of grapevine worldwide. In order to better understand the cellular processes involved in resistance to downy mildew, leaves of Vitis vinifera cv. ‘Centennial Seedless’ (susceptible) and Vitis labrusca×Vitis riparia cv. ‘Beta’ (resistant) were inoculated with Plasmopara viticola, and were used in the Illumina HiSeq™2000 platform for deep transcriptome sequencing. We performed transcriptome analysis and identified a total of 1091 differentially expressed genes (DEGs) in the disease-resistant cultivar and 849 in the susceptible cultivar. To increase understanding of the DEGs, the datasets were analyzed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genome (KEGG) database. Most of them were found to be associated with plant–pathogen interaction, flavonoid synthesis, phenylpropanoid synthesis metabolism and phytohormone signal transduction. Moreover, transcription factors ERF, MYB, WRKY, and bHLH associated with disease resistance were screened. A total of 196 genes were identified as the candidate resistant genes. The expression of 18 differentially expressed genes was detected by RT-qPCR. Finally, the pattern of differentially expressed genes was consistent with the result of transcriptome sequencing. The present study identified several candidate resistance genes and signal transduction pathways that may contribute to downy mildew resistance in grapes.

The Chromosome-Level Reference Genome of Tea Tree Unveils Recent Bursts of Non-autonomous LTR Retrotransposons in Driving Genome Size Evolution

The Chromosome-Level Reference Genome of Tea Tree Unveils Recent Bursts of Non-autonomous LTR Retrotransposons in Driving Genome Size Evolution

MicroRNA neural networks improve diagnosis of acute coronary syndrome (ACS)

BackgroundCardiac troponins are the preferred biomarkers of acute myocardial infarction. Despite superior sensitivity, serial testing of Troponins to identify patients suffering acute coronary syndromes is still required in many cases to overcome limited specificity. Moreover, unstable angina pectoris relies on reported symptoms in the troponin-negative group. In this study, we investigated genome-wide miRNA levels in a prospective cohort of patients with clinically suspected ACS and determined their diagnostic value by applying an in silico neural network. MethodsPAXgene blood and serum samples were drawn and hsTnT was measured in patients at initial presentation to our Chest-Pain Unit. After clinical and diagnostic workup, patients were adjudicated by senior cardiologists in duty to their final diagnosis: STEMI, NSTEMI, unstable angina pectoris and non-ACS patients. ACS patients and a cohort of healthy controls underwent deep transcriptome sequencing. Machine learning was implemented to construct diagnostic miRNA classifiers. ResultsWe developed a neural network model which incorporates 34 validated ACS miRNAs, showing excellent classification results. By further developing additional machine learning models and selecting the best miRNAs, we achieved an accuracy of 0.96 (95% CI 0.96–0.97), sensitivity of 0.95, specificity of 0.96 and AUC of 0.99. The one-point hsTnT value reached an accuracy of 0.89, sensitivity of 0.82, specificity of 0.96, and AUC of 0.96. ConclusionsHere we show the concept of neural network based biomarkers for ACS. This approach also opens the possibility to include multi-modal data points to further increase precision and perform classification of other ACS differential diagnoses.

Deep Transcriptome Sequencing of Pediatric Acute Myeloid Leukemia Patients at Diagnosis, Remission and Relapse: Experience in 3 Malaysian Children in a Single Center Study.

Deep Transcriptome Sequencing of Pediatric Acute Myeloid Leukemia Patients at Diagnosis, Remission and Relapse: Experience in 3 Malaysian Children in a Single Center Study.