Single-Nucleotide Variations,Insertions/Deletions and Copy Number Variations in Myelodysplastic Syndrome during Disease Progression Revealed by a Single-Cell DNA Sequencing Platform

Simple SummaryMyelodysplastic syndrome (MDS) is a myeloid neoplasm associated with complex clonal architecture. The application of single-cell sequencing is capable of revealing the clonal dynamics of MDS during disease progression and treatment resistance. This has advantages over bulk-tumor sequencing which is limited by its resolution. In this study, we evaluated two patients with MDS for the clonal dynamics of pathogenic mutations at the single-cell level of disease progression and resistance to hypomethylating agents (HMAs). There were two key observations. First, changes in the clonal heterogeneity of the pathogenic FLT3-ITD, IDH2, EZH2, or GATA2 mutations was associated with disease progression and resistance to HMA. Secondly, disease progression and resistance to HMA was accompanied by the acquisition of copy number variations of DNMT3A, TET2, and GATA2.AbstractMyelodysplastic syndrome (MDS) is a clonal myeloid neoplasm characterized by ineffective hematopoiesis, cytopenia, dysplasia, and clonal instability, leading to leukemic transformation. Hypomethylating agents are the mainstay of treatment in higher-risk MDS. However, treatment resistance and disease transformation into acute myeloid leukemia (AML) is observed in the majority of patients and is indicative of a dismal outcome. The residual cell clones resistant to therapy or cell clones acquiring new genetic aberrations are two of the key events responsible for drug resistance. Bulk tumor sequencing often fails to detect these rare subclones that confer resistance to therapy. In this study, we employed a single-cell DNA (sc-DNA) sequencing approach to study the clonal heterogeneity and clonal evolution in two MDS patients refractory to HMA. In both patients, different single nucleotide variations (SNVs) or insertions and deletions (INDELs) were detected with bulk tumor sequencing. Rare cell clones with mutations that are undetectable by bulk tumor sequencing were detected by sc-DNA sequencing. In addition to SNVs and short INDELs, this study also revealed the presence of a clonal copy number loss of DNMT3A, TET2, and GATA2 as standalone events or in association with the small SNVs or INDELs detected during HMA resistance and disease progression.     

Altered Pre-mRNA Splicing Caused by a Novel Intronic Mutation c.1443+5G>A in the Dihydropyrimidinase (DPYS) Gene

Dihydropyrimidinase (DHP) deficiency is an autosomal recessive disease caused by mutations in the DPYS gene. Patients present with highly elevated levels of dihydrouracil and dihydrothymine in their urine, blood and cerebrospinal fluid. The analysis of the effect of mutations in DPYS on pre-mRNA splicing is hampered by the fact that DHP is primarily expressed in liver and kidney cells. The minigene approach can detect mRNA splicing aberrations using cells that do not express the endogenous mRNA. We have used a minigene-based approach to analyze the effects of a presumptive pre-mRNA splicing mutation in two newly identified Chinese pediatric patients with DHP deficiency. Mutation analysis of DPYS showed that both patients were compound heterozygous for a novel intronic mutation c.1443+5G>A in intron 8 and a previously described missense mutation c.1001A>G (p.Q334R) in exon 6. Wild-type and the mutated minigene constructs, containing exons 7, 8 and 9 of DPYS, yielded different splicing products after expression in HEK293 cells. The c.1443+5G>A mutation resulted in altered pre-mRNA splicing of the DPYS minigene construct with full skipping of exon 8. Analysis of the DHP crystal structure showed that the deletion of exon 8 severely affects folding, stability and homooligomerization of the enzyme as well as disruption of the catalytic site. Thus, the analysis suggests that the c.1443+5G>A mutation results in aberrant splicing of the pre-mRNA encoding DHP, underlying the DHP deficiency in two unrelated Chinese patients.     

A Regulatory Axis between Epithelial Splicing Regulatory Proteins and Estrogen Receptor α Modulates the Alternative Transcriptome of Luminal Breast Cancer

Epithelial splicing regulatory proteins 1 and 2 (ESRP1/2) control the splicing pattern during epithelial to mesenchymal transition (EMT) in a physiological context and in cancer, including breast cancer (BC). Here, we report that ESRP1, but not ESRP2, is overexpressed in luminal BCs of patients with poor prognosis and correlates with estrogen receptor α (ERα) levels. Analysis of ERα genome-binding profiles in cell lines and primary breast tumors showed its binding in the proximity of ESRP1 and ESRP2 genes, whose expression is strongly decreased by ERα silencing in hormone-deprived conditions. The combined knock-down of ESRP1/2 in MCF-7 cells followed by RNA-Seq, revealed the dysregulation of 754 genes, with a widespread alteration of alternative splicing events (ASEs) of genes involved in cell signaling, metabolism, cell growth, and EMT. Functional network analysis of ASEs correlated with ESRP1/2 expression in ERα+ BCs showed RAC1 as the hub node in the protein–protein interactions altered by ESRP1/2 silencing. The comparison of ERα- and ESRP-modulated ASEs revealed 63 commonly regulated events, including 27 detected in primary BCs and endocrine-resistant cell lines. Our data support a functional implication of the ERα-ESRP1/2 axis in the onset and progression of BC by controlling the splicing patterns of related genes.     

Sex-Related Differences in Chronic Myeloid Neoplasms: From the Clinical Observation to the Underlying Biology

Chronic myeloid neoplasms are clonal diseases with variable clinical course and outcomes and despite the introduction of novel therapies, patients with high-risk disease continue to have overall poor outcomes. Different groups have highlighted that men have overall worse survival and higher incidence of transformation to acute leukemia compared to women across neoplasms such as myelodysplastic syndrome (MDS), myeloproliferative neoplasms (MPN), MDS/MPN overlap neoplasms, and CML. More recent studies evaluating the genomic profile of patients with these neoplasms demonstrated a male predominance for mutations in high-risk genes including ASXL1, U2AF1, SRSF2 and ZRSR2. The understanding of the underlying biology is limited but a number of hypotheses have been developed and are currently being investigated. This review summarizes the current knowledge about sex-related differences in the clinical outcomes and genomic profile of patients with chronic myeloid neoplasms and discusses the hypothesized biologic mechanisms as an attempt to explain these observations.     

BRCA1 Exon 11, a CERES (Composite Regulatory Element of Splicing) Element Involved in Splice Regulation

Unclassified variants (UV) of BRCA1 can affect normal pre-mRNA splicing. Here, we investigate the UV c.693G>A, a “silent” change in BRCA1 exon 11, which we have found induces aberrant splicing in patient carriers and in vitro. Using a minigene assay, we show that the UV c.693G>A has a strong effect on the splicing isoform ratio of BRCA1. Systematic site-directed mutagenesis of the area surrounding the nucleotide position c.693G>A induced variable changes in the level of exon 11 inclusion/exclusion in the mRNA, pointing to the presence of a complex regulatory element with overlapping enhancer and silencer functions. Accordingly, protein binding analysis in the region detected several splicing regulatory factors involved, including SRSF1, SRSF6 and SRSF9, suggesting that this sequence represents a composite regulatory element of splicing (CERES).     

The SR Splicing Factors: Providing Perspectives on Their Evolution,Expression, Alternative Splicing,and Function in Populus trichocarpa

Serine/arginine-rich (SR) proteins are important splicing factors in plant development and abiotic/hormone-related stresses. However, evidence that SR proteins contribute to the process in woody plants has been lacking. Using phylogenetics, gene synteny, transgenic experiments, and RNA-seq analysis, we identified 24 PtSR genes and explored their evolution, expression, and function in Popolus trichocarpa. The PtSR genes were divided into six subfamilies, generated by at least two events of genome triplication and duplication. Notably, they were constitutively expressed in roots, stems, and leaves, demonstrating their fundamental role in P. trichocarpa. Additionally, most PtSR genes (~83%) responded to at least one stress (cold, drought, salt, SA, MeJA, or ABA), and, especially, cold stress induced a dramatic perturbation in the expression and/or alternative splicing (AS) of 18 PtSR genes (~75%). Evidentially, the overexpression of PtSCL30 in Arabidopsis decreased freezing tolerance, which probably resulted from AS changes of the genes (e.g., ICE2 and COR15A) critical for cold tolerance. Moreover, the transgenic plants were salt-hypersensitive at the germination stage. These indicate that PtSCL30 may act as a negative regulator under cold and salt stress. Altogether, this study sheds light on the evolution, expression, and AS of PtSR genes, and the functional mechanisms of PtSCL30 in woody plants.     

The Impact of Epigenetic Modifications in Myeloid Malignancies

Myeloid malignancy is a broad term encapsulating myeloproliferative neoplasms (MPN), myelodysplastic syndrome (MDS) and acute myeloid leukaemia (AML). Initial studies into genomic profiles of these diseases have shown 2000 somatic mutations prevalent across the spectrum of myeloid blood disorders. Epigenetic mutations are emerging as critical components of disease progression, with mutations in genes controlling chromatin regulation and methylation/acetylation status. Genes such as DNA methyltransferase 3A (DNMT3A), ten eleven translocation methylcytosine dioxygenase 2 (TET2), additional sex combs-like 1 (ASXL1), enhancer of zeste homolog 2 (EZH2) and isocitrate dehydrogenase 1/2 (IDH1/2) show functional impact in disease pathogenesis. In this review we discuss how current knowledge relating to disease progression, mutational profile and therapeutic potential is progressing and increasing understanding of myeloid malignancies.     

Functional Analysis of the PCCA and PCCB Gene Variants Predicted to Affect Splicing

It is estimated that up to one-third of all variants causing inherited diseases affect splicing; however, their deleterious effects and roles in disease pathogenesis are often not fully characterized. Given their prevalence and the development of various antisense-based splice-modulating approaches, pathogenic splicing variants have become an important object of genomic medicine. To improve the accuracy of variant interpretation in public mutation repositories, we applied the minigene splicing assay to study the effects of 24 variants that were predicted to affect normal splicing in the genes associated with propionic acidemia (PA)—PCCA and PCCB. As a result, 13 variants (including one missense and two synonymous variants) demonstrated a significant alteration of splicing with the predicted deleterious effect at the protein level and were characterized as spliceogenic loss-of-function variants. The analysis of the available data for the studied variants and application of the American College of Medical Genetics and the Association for Molecular Pathology (ACMG/AMP) guidelines allowed us to precisely classify five of the variants and change the pathogenic status of nine. Using the example of the PA genes, we demonstrated the utility of the minigene splicing assay in the fast and effective assessment of the spliceogenic effect for identified variants and highlight the necessity of their standardized classification.     

Alternative Splicing Regulation of Low-Frequency Genetic Variants in Exon 2 of TREM2 in Alzheimer’s Disease by Splicing-Based Aggregation

TREM2 is among the most well-known Alzheimer’s disease (AD) risk genes; however, the functional roles of its AD-associated variants remain to be elucidated, and most known risk alleles are low-frequency variants whose investigation is challenging. Here, we utilized a splicing-guided aggregation method in which multiple low-frequency TREM2 variants were bundled together to investigate the functional impact of those variants on alternative splicing in AD. We analyzed whole genome sequencing (WGS) and RNA-seq data generated from cognitively normal elderly controls (CN) and AD patients in two independent cohorts, representing three regions in the frontal lobe of the human brain: the dorsolateral prefrontal cortex (CN = 213 and AD = 376), frontal pole (CN = 72 and AD = 175), and inferior frontal (CN = 63 and AD = 157). We observed an exon skipping event in the second exon of TREM2, with that exon tending to be more frequently skipped (p = 0.0012) in individuals having at least one low-frequency variant that caused loss-of-function for a splicing regulatory element. In addition, genes differentially expressed between AD patients with high vs. low skipping of the second exon (i.e., loss of a TREM2 functional domain) were significantly enriched in immune-related pathways. Our splicing-guided aggregation method thus provides new insight into the regulation of alternative splicing of the second exon of TREM2 by low-frequency variants and could be a useful tool for further exploring the potential molecular mechanisms of multiple, disease-associated, low-frequency variants.