scholarly journals RNA profiling of human testicular cells identifies syntenic lncRNAs associated with spermatogenesis

2019 ◽  
Vol 34 (7) ◽  
pp. 1278-1290 ◽  
Author(s):  
A D Rolland ◽  
B Evrard ◽  
T A Darde ◽  
C Le Béguec ◽  
Y Le Bras ◽  
...  
Keyword(s):  
Gene Expression ◽  
Rna Sequencing ◽  
Large Scale ◽  
Physiological State ◽  
Hormonal Treatment ◽  
Specific Gene ◽  
Sequencing Analysis ◽  
Rna Profiling ◽  
Testicular Cells ◽  
Human Spermatogenesis

Abstract STUDY QUESTION Is the noncoding transcriptional landscape during spermatogenesis conserved between human and rodents? SUMMARY ANSWER We identified a core group of 113 long noncoding RNAs (lncRNAs) and 20 novel genes dynamically and syntenically transcribed during spermatogenesis. WHAT IS KNOWN ALREADY Spermatogenesis is a complex differentiation process driven by a tightly regulated and highly specific gene expression program. Recently, several studies in various species have established that a large proportion of known lncRNAs are preferentially expressed during meiosis and spermiogenesis in a testis-specific manner. STUDY DESIGN, SIZE, DURATION To further investigate lncRNA expression in human spermatogenesis, we carried out a cross-species RNA profiling study using isolated testicular cells. PARTICIPANTS/MATERIALS, SETTING, METHODS Human testes were obtained from post-mortem donors (N = 8, 51 years old on average) or from prostate cancer patients with no hormonal treatment (N = 9, 80 years old on average) and only patients with full spermatogenesis were used to prepare enriched populations of spermatocytes, spermatids, Leydig cells, peritubular cells and Sertoli cells. To minimize potential biases linked to inter-patient variations, RNAs from two or three donors were pooled prior to RNA-sequencing (paired-end, strand-specific). Resulting reads were mapped to the human genome, allowing for assembly and quantification of corresponding transcripts. MAIN RESULTS AND THE ROLE OF CHANCE Our RNA-sequencing analysis of pools of isolated human testicular cells enabled us to reconstruct over 25 000 transcripts. Among them we identified thousands of lncRNAs, as well as many previously unidentified genes (novel unannotated transcripts) that share many properties of lncRNAs. Of note is that although noncoding genes showed much lower synteny than protein-coding ones, a significant fraction of syntenic lncRNAs displayed conserved expression during spermatogenesis. LARGE SCALE DATA Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI Gene Expression Omnibus under accession number GSE74896. LIMITATIONS, REASONS FOR CAUTION Isolation procedures may alter the physiological state of testicular cells, especially for somatic cells, leading to substantial changes at the transcriptome level. We therefore cross-validated our findings with three previously published transcriptomic analyses of human spermatogenesis. Despite the use of stringent filtration criteria, i.e. expression cut-off of at least three fragments per kilobase of exon model per million reads mapped, fold-change of at least three and false discovery rate adjusted P-values of less than <1%, the possibility of assembly artifacts and false-positive transcripts cannot be fully ruled out. WIDER IMPLICATIONS OF THE FINDINGS For the first time, this study has led to the identification of a large number of conserved germline-associated lncRNAs that are potentially important for spermatogenesis and sexual reproduction. In addition to further substantiating the basis of the human testicular physiology, our study provides new candidate genes for male infertility of genetic origin. This is likely to be relevant for identifying interesting diagnostic and prognostic biomarkers and also potential novel therapeutic targets for male contraception. STUDY FUNDING/COMPETING INTEREST(S) This work was supported by l’Institut national de la santé et de la recherche médicale (Inserm); l’Université de Rennes 1; l’Ecole des hautes études en santé publique (EHESP); INERIS-STORM to B.J. [N 10028NN]; Rennes Métropole ‘Défis scientifiques émergents’ to F.C (2011) and A.D.R (2013). The authors have no competing financial interests.

scholarly journals The in vivo endothelial cell translatome is highly heterogeneous across vascular beds

2019 ◽  
Author(s):  
Audrey C.A. Cleuren ◽  
Martijn A. van der Ent ◽  
Hui Jiang ◽  
Kristina L. Hunker ◽  
Andrew Yee ◽  
...  
Keyword(s):  
Gene Expression ◽  
Endothelial Cells ◽  
Single Cell ◽  
Blood Vessels ◽  
Rna Sequencing ◽  
Specific Gene ◽  
Sequencing Analysis ◽  
Single Cell Rna Sequencing ◽  
Vascular Beds

AbstractEndothelial cells (ECs) are highly specialized across vascular beds. However, given their interspersed anatomic distribution, comprehensive characterization of the molecular basis for this heterogeneity in vivo has been limited. By applying endothelial-specific translating ribosome affinity purification (EC-TRAP) combined with high-throughput RNA sequencing analysis, we identified pan EC-enriched genes and tissue-specific EC transcripts, which include both established markers and genes previously unappreciated for their presence in ECs. In addition, EC-TRAP limits changes in gene expression following EC isolation and in vitro expansion, as well as rapid vascular bed-specific shifts in EC gene expression profiles as a result of the enzymatic tissue dissociation required to generate single cell suspensions for fluorescence-activated cell sorting (FACS) or single cell RNA sequencing analysis. Comparison of our EC-TRAP to published single cell RNA sequencing data further demonstrates considerably greater sensitivity of EC-TRAP for the detection of low abundant transcripts. Application of EC-TRAP to examine the in vivo host response to lipopolysaccharide (LPS) revealed the induction of gene expression programs associated with a native defense response, with marked differences across vascular beds. Furthermore, comparative analysis of whole tissue and TRAP-selected mRNAs identified LPS-induced differences that would not have been detected by whole tissue analysis alone. Together, these data provide a resource for the analysis of EC-specific gene expression programs across heterogeneous vascular beds under both physiologic and pathologic conditions.SignificanceEndothelial cells (ECs), which line all vertebrate blood vessels, are highly heterogeneous across different tissues. The present study uses a genetic approach to specifically tag mRNAs within ECs of the mouse, thereby allowing recovery and sequence analysis to evaluate the EC-specific gene expression program directly from intact organs. Our findings demonstrate marked heterogeneity in EC gene expression across different vascular beds under both normal and disease conditions, with a more accurate picture than can be achieved using other methods. The data generated in these studies advance our understanding of EC function in different blood vessels and provide a valuable resource for future studies.

scholarly journals The in vivo endothelial cell translatome is highly heterogeneous across vascular beds

2019 ◽  
Vol 116 (47) ◽  
pp. 23618-23624 ◽  
Author(s):  
Audrey C. A. Cleuren ◽  
Martijn A. van der Ent ◽  
Hui Jiang ◽  
Kristina L. Hunker ◽  
Andrew Yee ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Rna Sequencing ◽  
Expression Profiles ◽  
Affinity Purification ◽  
Specific Gene ◽  
Sequencing Analysis ◽  
Single Cell Rna Sequencing ◽  
Vascular Beds

Endothelial cells (ECs) are highly specialized across vascular beds. However, given their interspersed anatomic distribution, comprehensive characterization of the molecular basis for this heterogeneity in vivo has been limited. By applying endothelial-specific translating ribosome affinity purification (EC-TRAP) combined with high-throughput RNA sequencing analysis, we identified pan EC-enriched genes and tissue-specific EC transcripts, which include both established markers and genes previously unappreciated for their presence in ECs. In addition, EC-TRAP limits changes in gene expression after EC isolation and in vitro expansion, as well as rapid vascular bed-specific shifts in EC gene expression profiles as a result of the enzymatic tissue dissociation required to generate single-cell suspensions for fluorescence-activated cell sorting or single-cell RNA sequencing analysis. Comparison of our EC-TRAP with published single-cell RNA sequencing data further demonstrates considerably greater sensitivity of EC-TRAP for the detection of low abundant transcripts. Application of EC-TRAP to examine the in vivo host response to lipopolysaccharide (LPS) revealed the induction of gene expression programs associated with a native defense response, with marked differences across vascular beds. Furthermore, comparative analysis of whole-tissue and TRAP-selected mRNAs identified LPS-induced differences that would not have been detected by whole-tissue analysis alone. Together, these data provide a resource for the analysis of EC-specific gene expression programs across heterogeneous vascular beds under both physiologic and pathologic conditions.

scholarly journals The combined detection of Amphiregulin, Cyclin A1 and DDX20/Gemin3 expression predicts aggressive forms of oral squamous cell carcinoma

2021 ◽  
Author(s):  
Ekaterina Bourova-Flin ◽  
Samira Derakhshan ◽  
Afsaneh Goudarzi ◽  
Tao Wang ◽  
Anne-Laure Vitte ◽  
...  
Keyword(s):  
Gene Expression ◽  
Squamous Cell ◽  
Large Scale ◽  
Biomarker Discovery ◽  
Gene Expression Signature ◽  
Predictive Biomarkers ◽  
Specific Gene ◽  
Specific Expression ◽  
Intrinsic Nature ◽  
Independent Cohort

Abstract Background Large-scale genetic and epigenetic deregulations enable cancer cells to ectopically activate tissue-specific expression programmes. A specifically designed strategy was applied to oral squamous cell carcinomas (OSCC) in order to detect ectopic gene activations and develop a prognostic stratification test. Methods A dedicated original prognosis biomarker discovery approach was implemented using genome-wide transcriptomic data of OSCC, including training and validation cohorts. Abnormal expressions of silent genes were systematically detected, correlated with survival probabilities and evaluated as predictive biomarkers. The resulting stratification test was confirmed in an independent cohort using immunohistochemistry. Results A specific gene expression signature, including a combination of three genes, AREG, CCNA1 and DDX20, was found associated with high-risk OSCC in univariate and multivariate analyses. It was translated into an immunohistochemistry-based test, which successfully stratified patients of our own independent cohort. Discussion The exploration of the whole gene expression profile characterising aggressive OSCC tumours highlights their enhanced proliferative and poorly differentiated intrinsic nature. Experimental targeting of CCNA1 in OSCC cells is associated with a shift of transcriptomic signature towards the less aggressive form of OSCC, suggesting that CCNA1 could be a good target for therapeutic approaches.

Abstract W P73: Feasibility and Preliminary Results of Whole Blood RNA-Sequencing Analysis in Patients With Intracranial Aneurysms

Stroke ◽  
2014 ◽  
Vol 45 (suppl_1) ◽  
Author(s):  
Blake Haas ◽  
Nestor R Gonzalez ◽  
Elina Nikkola ◽  
Mark Connolly ◽  
William Hsu ◽  
...  
Keyword(s):  
Gene Expression ◽  
Network Analysis ◽  
Rna Sequencing ◽  
Whole Blood ◽  
Intracranial Aneurysms ◽  
Small Sample Size ◽  
Expression Patterns ◽  
Small Sample ◽  
Cellular Respiration ◽  
Sequencing Analysis

Introduction: Intracranial aneurysms (IA) growth and rupture have been associated with chronic remodeling of the arterial wall. However, the pathobiology of this process remains poorly understood. The objective of the present study was to evaluate the feasibility of analyzing gene expression patterns in peripheral blood of patients with ruptured and unruptured saccular IAs. Materials and Methods: We analyzed human whole blood transcriptomes by performing paired-end, 100 bp RNA-sequencing (RNAseq) using the Illumina platform. We used STAR to align reads to the genome, HTSeq to count reads, and DESeq to normalize counts across samples. Self-reported patient information was used to correct expression values for ancestry, age, and sex. We utilized weighted gene co-expression network analysis (WGCNA) to identify gene expression network modules associated with IA size and rupture. The DAVID tool was employed to search for Gene Ontology enrichment in relevant modules. Results: Samples from 12 patients (9 females, age 57.6 +/-12) with IAs were analyzed. Four had ruptured aneurysms. RNA isolation and application of the methodology described above was successful in all samples. Although the small sample size prevents us from drawing definite conclusions, we observed promising novel co-expression networks for IAs: WCGNA analysis showed down-regulation of two transcript modules associated with ruptured IA status (r=-0.78, p=0.008 and r=-0.77, p=0.009), and up-regulation of two modules associated with aneurysm size (r=0.86, p=0.002 and r=0.9, p=4e-04), respectively. DAVID analyses showed that genes upregulated in an IA size-associated module were enriched with genes involved in cellular respiration and translation, while genes involved in transcription were down-regulated in a module associated with ruptured IAs. Conclusions: Whole blood RNAseq analysis is a feasible tool to capture transcriptome dynamics and achieve a better understanding of the pathophysiology of IAs. Further longitudinal studies of patients with IAs using network analysis are justified.

FC 011KIDNEYNETWORK: USING KIDNEY DERIVED GENE EXPRESSION DATA TO PREDICT AND PRIORITIZE NOVEL GENES INVOLVED IN KIDNEY DISEASE

2021 ◽  
Vol 36 (Supplement_1) ◽  
Author(s):  
Floranne Boulogne ◽  
Laura Claus ◽  
Henry Wiersma ◽  
Roy Oelen ◽  
Floor Schukking ◽  
...  
Keyword(s):  
Gene Expression ◽  
Kidney Disease ◽  
Candidate Gene ◽  
Exome Sequencing ◽  
Rna Sequencing ◽  
Expression Patterns ◽  
Genetic Diagnosis ◽  
Specific Gene ◽  
Sequencing Data ◽  
Exome Sequencing Data

Abstract Background and Aims Genetic testing in patients with suspected hereditary kidney disease does not always reveal the genetic cause for the patient's disorder. Potentially pathogenic variants can reside in genes that are not known to be involved in kidney disease, which makes it difficult to prioritize and interpret the relevance of these variants. As such, there is a clear need for methods that predict the phenotypic consequences of gene expression in a way that is as unbiased as possible. To help identify candidate genes we have developed KidneyNetwork, in which tissue-specific expression is utilized to predict kidney-specific gene functions. Method We combined gene co-expression in 878 publicly available kidney RNA-sequencing samples with the co-expression of a multi-tissue RNA-sequencing dataset of 31,499 samples to build KidneyNetwork. The expression patterns were used to predict which genes have a kidney-related function, and which (disease) phenotypes might be caused when these genes are mutated. By integrating the information from the HPO database, in which known phenotypic consequences of disease genes are annotated, with the gene co-expression network we obtained prediction scores for each gene per HPO term. As proof of principle, we applied KidneyNetwork to prioritize variants in exome-sequencing data from 13 kidney disease patients without a genetic diagnosis. Results We assessed the prediction performance of KidneyNetwork by comparing it to GeneNetwork, a multi-tissue co-expression network we previously developed. In KidneyNetwork, we observe a significantly improved prediction accuracy of kidney-related HPO-terms, as well as an increase in the total number of significantly predicted kidney-related HPO-terms (figure 1). To examine its clinical utility, we applied KidneyNetwork to 13 patients with a suspected hereditary kidney disease without a genetic diagnosis. Based on the HPO terms “Renal cyst” and “Hepatic cysts”, combined with a list of potentially damaging variants in one of the undiagnosed patients with mild ADPKD/PCLD, we identified ALG6 as a new candidate gene. ALG6 bears a high resemblance to other genes implicated in this phenotype in recent years. Through the 100,000 Genomes Project and collaborators we identified three additional patients with kidney and/or liver cysts carrying a suspected deleterious variant in ALG6. Conclusion We present KidneyNetwork, a kidney specific co-expression network that accurately predicts what genes have kidney-specific functions and may result in kidney disease. Gene-phenotype associations of genes unknown for kidney-related phenotypes can be predicted by KidneyNetwork. We show the added value of KidneyNetwork by applying it to exome sequencing data of kidney disease patients without a molecular diagnosis and consequently we propose ALG6 as a promising candidate gene. KidneyNetwork can be applied to clinically unsolved kidney disease cases, but it can also be used by researchers to gain insight into individual genes to better understand kidney physiology and pathophysiology. Acknowledgments This research was made possible through access to the data and findings generated by the 100,000 Genomes Project; http://www.genomicsengland.co.uk.

scholarly journals LSTrAP-Crowd: Prediction of novel components of bacterial ribosomes with crowd-sourced analysis of RNA sequencing data

2020 ◽  
Author(s):  
Benedict Hew ◽  
Qiao Wen Tan ◽  
William Goh ◽  
Jonathan Wei Xiong Ng ◽  
Kenny Koh ◽  
...  
Keyword(s):  
Gene Expression ◽  
Protein Synthesis ◽  
Rna Sequencing ◽  
Gene Expression Data ◽  
Large Scale ◽  
Bacterial Resistance ◽  
Expression Data ◽  
Sequencing Data ◽  
Novel Proteins ◽  
Novel Antibiotics

AbstractBacterial resistance to antibiotics is a growing problem that is projected to cause more deaths than cancer in 2050. Consequently, novel antibiotics are urgently needed. Since more than half of the available antibiotics target the bacterial ribosomes, proteins that are involved in protein synthesis are thus prime targets for the development of novel antibiotics. However, experimental identification of these potential antibiotic target proteins can be labor-intensive and challenging, as these proteins are likely to be poorly characterized and specific to few bacteria. In order to identify these novel proteins, we established a Large-Scale Transcriptomic Analysis Pipeline in Crowd (LSTrAP-Crowd), where 285 individuals processed 26 terabytes of RNA-sequencing data of the 17 most notorious bacterial pathogens. In total, the crowd processed 26,269 RNA-seq experiments and used the data to construct gene co-expression networks, which were used to identify more than a hundred uncharacterized genes that were transcriptionally associated with protein synthesis. We provide the identity of these genes together with the processed gene expression data. The data can be used to identify other vulnerabilities or bacteria, while our approach demonstrates how the processing of gene expression data can be easily crowdsourced.

scholarly journals Importance of experimental information (metadata) for archived sequence data: case of specific gene bias due to lag time between sample harvest and RNA protection in RNA sequencing

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e11875
Author(s):  
Tomoko Matsuda
Keyword(s):  
Gene Expression ◽  
Rna Sequencing ◽  
Time Course ◽  
Sequence Data ◽  
Specific Gene ◽  
Time Interval ◽  
Short Time Interval ◽  
Rna Seq ◽  
Lysis Buffer ◽  
Rna Protection

Large volumes of high-throughput sequencing data have been submitted to the Sequencing Read Archive (SRA). The lack of experimental metadata associated with the data makes reuse and understanding data quality very difficult. In the case of RNA sequencing (RNA-Seq), which reveals the presence and quantity of RNA in a biological sample at any moment, it is necessary to consider that gene expression responds over a short time interval (several seconds to a few minutes) in many organisms. Therefore, to isolate RNA that accurately reflects the transcriptome at the point of harvest, raw biological samples should be processed by freezing in liquid nitrogen, immersing in RNA stabilization reagent or lysing and homogenizing in RNA lysis buffer containing guanidine thiocyanate as soon as possible. As the number of samples handled simultaneously increases, the time until the RNA is protected can increase. Here, to evaluate the effect of different lag times in RNA protection on RNA-Seq data, we harvested CHO-S cells after 3, 5, 6, and 7 days of cultivation, added RNA lysis buffer in a time course of 15, 30, 45, and 60 min after harvest, and conducted RNA-Seq. These RNA samples showed high RNA integrity number (RIN) values indicating non-degraded RNA, and sequence data from libraries prepared with these RNA samples was of high quality according to FastQC. We observed that, at the same cultivation day, global trends of gene expression were similar across the time course of addition of RNA lysis buffer; however, the expression of some genes was significantly different between the time-course samples of the same cultivation day; most of these differentially expressed genes were related to apoptosis. We conclude that the time lag between sample harvest and RNA protection influences gene expression of specific genes. It is, therefore, necessary to know not only RIN values of RNA and the quality of the sequence data but also how the experiment was performed when acquiring RNA-Seq data from the database.

Sign in / Sign up

Export Citation Format

Share Document