A 3’-end capture sequencing method for high-throughput targeted gene expression profiling

Molecular phenotyping through shallow 3’-end RNA-sequencing workflows is increasingly applied in the context of large-scale chemical or genetic perturbation screens to study disease biology or support drug discovery. While these workflows enable accurate quantification of the most abundant genes, they are less effective for applications that require expression profiling of low abundant transcripts, like long non-coding RNAs (lncRNAs), or selected gene panels. To tackle these issues, we describe a workflow combining 3’-end library preparation with 3’-end hybrid capture probes and shallow RNA-sequencing for cost-effective, targeted quantification of subsets of (low abundant) genes across hundreds to thousands of samples. To assess the performance of the method, we designed a capture probe set for more than 100 mRNA and lncRNA target genes and applied the workflow to a cohort of 360 samples. When compared to standard 3’-end RNA-sequencing, 3’-end capture sequencing resulted in a more than 100-fold enrichment of target gene abundance while conserving relative inter-gene and inter-sample abundances. 3’-end RNA capture sequencing enables accurate targeted gene expression profiling at extremely shallow sequencing depth.

BIOM-40. TARGETED GENE EXPRESSION PROFILING PREDICTS MENINGIOMA OUTCOMES AND RADIOTHERAPY RESPONSES

BACKGROUND Surgery is the mainstay of meningioma treatment, but improvements in meningioma risk stratification are needed and indications for postoperative radiotherapy are controversial. DNA methylation profiling, copy number variants (CNVs), exome sequencing, and RNA sequencing have improved understanding of meningioma biology, but have not superseded histologic grading, or revealed biomarkers for radiotherapy responses. To address these unmet needs, we optimized and validated a targeted gene expression biomarker predicting meningioma outcomes and responses to radiotherapy. METHODS Targeted gene expression profiling was performed on a discovery cohort of 173 meningiomas (median follow-up 8.1 years) and a validation cohort of 331 meningiomas (median follow-up 6.1 years) treated with surgery (n=504) and postoperative radiotherapy (n=73) at independent, international institutions (70% WHO grade 1, 24% WHO grade 2, 6% WHO grade 3). Optimized targeted gene expression models predicting clinical outcomes (34 genes) or radiotherapy responses (12 genes) were developed from the discovery cohort, and compared to histologic and molecular classification systems by performing DNA methylation profiling, CNV analysis, exome sequencing, and RNA sequencing on the same meningiomas. RESULTS Targeted gene expression profiling achieved a concordance-index of 0.75 ± 0.03 (SEM) for local freedom from recurrence (LFFR) and 0.72 ± 0.03 for overall survival (OS) in the validation cohort, outperforming WHO grade (5-year LFFR delta-AUC 0.15, 95% CI 0.076-0.229, p=0.001) and DNA methylation grouping (delta-AUC 0.075, 95% CI 0.006-0.130, p=0.01) for LFFR, disease-specific survival, and OS. The biomarker was independently prognostic after accounting for WHO grade, extent of resection, primary versus recurrent presentation, CNV status, DNA methylation group, and Ki67 labeling index, and identified meningiomas benefiting from radiotherapy (interaction p-value=0.0008), suggesting postoperative radiotherapy could be refined in 30.2% of cases. CONCLUSIONS Targeted gene expression profiling of 504 meningiomas improves discrimination of meningioma local recurrence, disease-specific survival, and overall survival, and predicts radiotherapy responses.

Kidney Injury Monitoring in Tobramycin-Treated Rhesus Monkeys: Supplementing Urinary Kidney Biomarkers With Kidney Biopsy Gene Expression Profiling

Kidney biopsies are used sparingly to diagnose kidney injury in the clinic. Here we have conducted a small exploratory study to directly compare the low-grade kidney injury monitoring performance of serum safety biomarkers, novel urine safety biomarkers, microscopic histopathology and targeted gene expression alterations in kidney biopsy specimens in rhesus monkeys treated with tobramycin. Targeted gene expression increases were observed in the kidney biopsy samples and whole kidney sections for kidney injury molecule 1 ( KIM-1), clusterin ( CLU), osteopontin ( OPN) messenger RNA transcripts. In addition, increases of the urinary kidney safety protein biomarkers including KIM-1, CLU, OPN were also observed. These increases in gene expression and urinary protein end point were in concordance with the eventual low-grade kidney lesions seen in terminal tissue sections. In contrast, conventional serum biomarkers blood urea nitrogen and serum creatinine were not as sensitive in monitoring kidney injury. Although these data do not support routinely adding kidney biopsies to regular toxicology studies, they provide evidence on the value and limitations of incorporating gene expression profiling on kidney biopsy specimens, further underscore the value of urinary kidney safety biomarkers for improved low-grade kidney injury monitoring, and open the door for future definitive studies.

Molecular and Biochemical Insights Into Early Responses of Hemp to Cd and Zn Exposure and the Potential Effect of Si on Stress Response

With the intensification of human activities, plants are more frequently exposed to heavy metals (HM). Zinc (Zn) and cadmium (Cd) are frequently and simultaneously found in contaminated soils, including agronomic soils contaminated by the atmospheric fallout near smelters. The fiber crop Cannabis sativa L. is a suitable alternative to food crops for crop cultivation on these soils. In this study, Cd (20 μM) and Zn (100 μM) were shown to induce comparable growth inhibition in C. sativa. To devise agricultural strategies aimed at improving crop yield, the effect of silicon (Si; 2 mM) on the stress tolerance of plants was considered. Targeted gene expression and proteomic analysis were performed on leaves and roots after 1 week of treatment. Both Cd- and Zn-stimulated genes involved in proline biosynthesis [pyrroline-5-carboxylate reductase (P5CR)] and phenylpropanoid pathway [phenylalanine ammonia-lyase (PAL)] but Cd also specifically increased the expression of PCS1-1 involved in phytochelatin (PC) synthesis. Si exposure influences the expression of numerous genes in a contrasting way in Cd- and Zn-exposed plants. At the leaf level, the accumulation of 122 proteins was affected by Cd, whereas 47 proteins were affected by Zn: only 16 proteins were affected by both Cd and Zn. The number of proteins affected due to Si exposure (27) alone was by far lower, and 12 were not modified by heavy metal treatment while no common protein seemed to be modified by both CdSi and ZnSi treatment. It is concluded that Cd and Zn had a clear different impact on plant metabolism and that Si confers a specific physiological status to stressed plants, with quite distinct impacts on hemp proteome depending on the considered heavy metal.

Comparative Analysis of Differentially Expressed miRNAs in Leaves of Three Sugarcane (Saacharum Officinarum L.) Cultivars During Salinity Stress

Sugarcane is an important industrial plant which cultivated in the most arid and semiarid regions. Due to climate change and anthropogenic activities, the sugarcane field damage due to salt deposition and the cultivation of sugarcane has been posed a major threat in the region. To address this issue, the identification of salinity tolerant cultivars would be a suitable strategy to minimize yield loss in the area. MicroRNAs (miRNAs) play important roles in regulating gene expression. The monitoring of the expression of miRNAs and their targeted genes could provide deeper insight into the molecular stress mechanism and screen tolerant cultivars. Our aim was to assess the expression of nine candidate miRNAs and their corresponding targeted genes among the studied sugarcane cultivars under salinity condition, leading to identify the salt-tolerant cultivar. To achieve our goal, a two-factorial experiment with three sugarcane cultivars (CP-48, CP-57, CP-69) and two salinity levels (0 and 8 ds/m) was conducted. One-way ANOVA indicated that there was a significant difference between miRNAs and targeted gene expression. The highest reduction of miRNAs expression was occurred in miR160 while the lower one was happening in miR1432. The data also indicated that the higher and the lowest of targeted genes were in miR160 and miR393 respectively. Among studied cultivars, the CP-57 showed poor performance while CP-69 expresses a superior tolerance to salt stress. Taken together, these results suggested that the screening of well adapted cultivars under salt conditions would be appropriate solutions to combat salinity stress in saline lands.

Timed global reorganization of protein synthesis during neocortex neurogenesis at codon resolution

Translation modulates the timing and amplification of gene expression after transcription. Development of the brain's neocortex requires precisely timed and spatially targeted gene expression, but the relationship between mRNA vs. protein synthesis throughout the genome is unknown. We perform a comprehensive analysis of the reactants, synthesis, and products of mRNA translation spanning mouse neocortex neurogenesis. Ribosome number in the cortical plate decreases sharply at mid-neurogenesis during a transition in neuronal subtype specification, shifting the fundamental kinetics of protein synthesis, with mRNA and protein levels frequently divergent. Satb2, which drives an essential neuronal subtype-specific program, is a highly dynamically translated mRNA with surprisingly broad transcription across diverse neuronal lineages. Satb2 protein achieves its neuronal subtype expression through timed regulation by the RNA-binding protein Pumilio2. Thus, the refinement of transcriptional programs by protein synthesis is a widespread feature of neuronal specification. Developmental neocortex translatome data are provided in an open-source resource: https://shiny.mdc-berlin.de/cortexomics/.

Defining the Emerging Blood System During Development at Single-Cell Resolution

Developmental hematopoiesis differs from adult and is far less described. In the developing embryo, waves of lineage-restricted blood precede the ultimate emergence of definitive hematopoietic stem cells (dHSCs) capable of maintaining hematopoiesis throughout life. During the last two decades, the advent of single-cell genomics has provided tools to circumvent previously impeding characteristics of embryonic hematopoiesis, such as cell heterogeneity and rare cell states, allowing for definition of lineage trajectories, cellular hierarchies, and cell-type specification. The field has rapidly advanced from microfluidic platforms and targeted gene expression analysis, to high throughput unbiased single-cell transcriptomic profiling, single-cell chromatin analysis, and cell tracing—offering a plethora of tools to resolve important questions within hematopoietic development. Here, we describe how these technologies have been implemented to address a wide range of aspects of embryonic hematopoiesis ranging from the gene regulatory network of dHSC formation via endothelial to hematopoietic transition (EHT) and how EHT can be recapitulated in vitro, to hematopoietic trajectories and cell fate decisions. Together, these studies have important relevance for regenerative medicine and for our understanding of genetic blood disorders and childhood leukemias.