Identifying interpretable gene-biomarker associations with functionally informed kernel-based tests in 190,000 exomes.

Here we present an exome-wide rare genetic variant association study for 30 blood biomarkers in 191,971 individuals in the UK Biobank. We compare gene-based association tests for separate functional variant categories to increase interpretability and identify 193 significant gene-biomarker associations. Genes associated with biomarkers were ~ 4.5-fold enriched for conferring Mendelian disorders. In addition to performing weighted gene-based variant collapsing tests, we design and apply variant-category-specific kernel-based tests that integrate quantitative functional variant effect predictions for missense variants, splicing and the binding of RNA-binding proteins. For these tests, we present a computationally efficient combination of the likelihood-ratio and score tests that found 36% more associations than the score test alone while also controlling the type-1 error. Kernel-based tests identified 13% more associations than their gene-based collapsing counterparts and had advantages in the presence of gain of function missense variants. We introduce local collapsing by amino acid position for missense variants and use it to interpret associations and identify potential novel gain of function variants in PIEZO1. Our results show the benefits of investigating different functional mechanisms when performing rare-variant association tests, and demonstrate pervasive rare-variant contribution to biomarker variability.

Long noncoding RNAs have pivotal roles in chemoresistance of acute myeloid leukemia.

Many patients with acute myeloid leukemia (AML) experience poor outcomes following traditional high-dose chemotherapies and complete remission rates remain suboptimal. Chemoresistance is an obstacle to effective chemotherapy and the precise mechanisms involved remain to be determined. Recently, long noncoding RNAs (lncRNAs) have been identified as relevant factors in the development of drug resistance in patients with AML. Furthermore, accumulating data support the importance of lncRNAs as potentially useful novel therapeutic targets in many cancers. Here, we review the role of lncRNAs in the development and induction of the chemoresistance in AML, and suggest lncRNAs as novel molecular markers for diagnosis, prediction of patient response to chemotherapy, and novel therapeutic targets for AML.

Author Correction: MEG3 long noncoding RNA regulates the TGF-ß pathway genes through formation of RNA-DNA triplex structures.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Long non-coding RNAs defining major subtypes of B cell precursor acute lymphoblastic leukemia.

BACKGROUND: Long non-coding RNAs (lncRNAs) have emerged as a novel class of RNA due to its diverse mechanism in cancer development and progression. However, the role and expression pattern of lncRNAs in molecular subtypes of B cell acute lymphoblastic leukemia (BCP-ALL) have not yet been investigated. Here, we assess to what extent lncRNA expression and DNA methylation is driving the progression of relapsed BCP-ALL subtypes and we determine if the expression and DNA methylation profile of lncRNAs correlates with established BCP-ALL subtypes. METHODS: We performed RNA sequencing and DNA methylation (Illumina Infinium microarray) of 40 diagnosis and 42 relapse samples from 45 BCP-ALL patients in a German cohort and quantified lncRNA expression. Unsupervised clustering was applied to ascertain and confirm that the lncRNA-based classification of the BCP-ALL molecular subtypes is present in both our cohort and an independent validation cohort of 47 patients. A differential expression and differential methylation analysis was applied to determine the subtype-specific, relapse-specific, and differentially methylated lncRNAs. Potential functions of subtype-specific lncRNAs were determined by using co-expression-based analysis on nearby (cis) and distally (trans) located protein-coding genes. RESULTS: Using an integrative Bioinformatics analysis, we developed a comprehensive catalog of 1235 aberrantly dysregulated BCP-ALL subtype-specific and 942 relapse-specific lncRNAs and the methylation profile of three subtypes of BCP-ALL. The 1235 subtype-specific lncRNA signature represented a similar classification of the molecular subtypes of BCP-ALL in the independent validation cohort. We identified a strong correlation between the DUX4-specific lncRNAs and genes involved in the activation of TGF-ß and Hippo signaling pathways. Similarly, Ph-like-specific lncRNAs were correlated with genes involved in the activation of PI3K-AKT, mTOR, and JAK-STAT signaling pathways. Interestingly, the relapse-specific lncRNAs correlated with the activation of metabolic and signaling pathways. Finally, we found 23 promoter methylated lncRNAs epigenetically facilitating their expression levels. CONCLUSION: Here, we describe a set of subtype-specific and relapse-specific lncRNAs from three major BCP-ALL subtypes and define their potential functions and epigenetic regulation. The subtype-specific lncRNAs are reproducible and can effectively stratify BCP-ALL subtypes. Our data uncover the diverse mechanism of action of lncRNAs in BCP-ALL subtypes defining which lncRNAs are involved in the pathogenesis of disease and are relevant for the stratification of BCP-ALL subtypes.

PAX5 biallelic genomic alterations define a novel subgroup of B-cell precursor acute lymphoblastic leukemia.

Chromosomal rearrangements and specific aneuploidy patterns are initiating events and define subgroups in B-cell precursor acute lymphoblastic leukemia (BCP-ALL). Here we analyzed 250 BCP-ALL cases and identified a novel subgroup ('PAX5-plus', n = 19) by distinct DNA methylation and gene expression profiles. All patients in this subgroup harbored mutations in the B-lineage transcription factor PAX5, with p.P80R as hotspot. Mutations either affected two independent codons, consistent with compound heterozygosity, or suffered LOH predominantly through chromosome 9p aberrations. These biallelic events resulted in disruption of PAX5 transcriptional programs regulating B-cell differentiation and tumor suppressor functions. Homozygous CDKN2A/B deletions and RAS-activating hotspot mutations were highly enriched as cooperating events in the genomic profile of PAX5-plus ALL. Together, this defined a specific pattern of triple alterations, exclusive to the novel subgroup. PAX5-plus ALL was observed in pediatric and adult patients. Although restricted by the limited sample size, a tendency for more favorable clinical outcome was observed, with 10 of 12 adult PAX5-plus patients achieving long-term survival. PAX5-plus represents the first BCP-ALL subgroup defined by sequence alterations in contrast to gross chromosomal events and exemplifies how deregulated differentiation (PAX5), impaired cell cycle control (CDKN2A/B) and sustained proliferative signaling (RAS) cooperatively drive leukemogenesis.

Integrated analysis of relapsed B-cell precursor Acute Lymphoblastic Leukemia identifies subtype-specific cytokine and metabolic signatures.

Recent efforts reclassified B-Cell Precursor Acute Lymphoblastic Leukemia (BCP-ALL) into more refined subtypes. Nevertheless, outcomes of relapsed BCP-ALL remain unsatisfactory, particularly in adult patients where the molecular basis of relapse is still poorly understood. To elucidate the evolution of relapse in BCP-ALL, we established a comprehensive multi-omics dataset including DNA-sequencing, RNA-sequencing, DNA methylation array and proteome MASS-spec data from matched diagnosis and relapse samples of BCP-ALL patients (n = 50) including the subtypes DUX4, Ph-like and two aneuploid subtypes. Relapse-specific alterations were enriched for chromatin modifiers, nucleotide and steroid metabolism including the novel candidates FPGS, AGBL and ZNF483. The proteome expression analysis unraveled deregulation of metabolic pathways at relapse including the key proteins G6PD, TKT, GPI and PGD. Moreover, we identified a novel relapse-specific gene signature specific for DUX4 BCP-ALL patients highlighting chemotaxis and cytokine environment as a possible driver event at relapse. This study presents novel insights at distinct molecular levels of relapsed BCP-ALL based on a comprehensive multi-omics integrated data set including a valuable proteomics data set. The relapse specific aberrations reveal metabolic signatures on genomic and proteomic levels in BCP-ALL relapse. Furthermore, the chemokine expression signature in DUX4 relapse underscores the distinct status of DUX4-fusion BCP-ALL.

MEG3 long noncoding RNA regulates the TGF-ß pathway genes through formation of RNA-DNA triplex structures.

Long noncoding RNAs (lncRNAs) regulate gene expression by association with chromatin, but how they target chromatin remains poorly understood. We have used chromatin RNA immunoprecipitation-coupled high-throughput sequencing to identify 276 lncRNAs enriched in repressive chromatin from breast cancer cells. Using one of the chromatin-interacting lncRNAs, MEG3, we explore the mechanisms by which lncRNAs target chromatin. Here we show that MEG3 and EZH2 share common target genes, including the TGF-ß pathway genes. Genome-wide mapping of MEG3 binding sites reveals that MEG3 modulates the activity of TGF-ß genes by binding to distal regulatory elements. MEG3 binding sites have GA-rich sequences, which guide MEG3 to the chromatin through RNA-DNA triplex formation. We have found that RNA-DNA triplex structures are widespread and are present over the MEG3 binding sites associated with the TGF-ß pathway genes. Our findings suggest that RNA-DNA triplex formation could be a general characteristic of target gene recognition by the chromatin-interacting lncRNAs.