Micro-macro consistency in multiscale modeling: Score-based model assisted sampling of fast/slow dynamical systems.

A valuable step in the modeling of multiscale dynamical systems in fields such as computational chemistry, biology, and materials science is the representative sampling of the phase space over long time scales of interest; this task is not, however, without challenges. For example, the long term behavior of a system with many degrees of freedom often cannot be efficiently computationally explored by direct dynamical simulation; such systems can often become trapped in local free energy minima. In the study of physics-based multi-time-scale dynamical systems, techniques have been developed for enhancing sampling in order to accelerate exploration beyond free energy barriers. On the other hand, in the field of machine learning (ML), a generic goal of generative models is to sample from a target density, after training on empirical samples from this density. Score-based generative models (SGMs) have demonstrated state-of-the-art capabilities in generating plausible data from target training distributions. Conditional implementations of such generative models have been shown to exhibit significant parallels with long-established-and physics-based-solutions to enhanced sampling. These physics-based methods can then be enhanced through coupling with the ML generative models, complementing the strengths and mitigating the weaknesses of each technique. In this work, we show that SGMs can be used in such a coupling framework to improve sampling in multiscale dynamical systems.

Aggressive systemic mastocytosis with the co-occurrence of PRKG2::PDGFRB, KAT6A::NCOA2, and RXRA::NOTCH1 fusion transcripts and a heterozygous RUNX1 frameshift mutation.

Systemic mastocytosis (SM) is a myeloproliferative neoplasm displaying abnormal mast cell proliferation. It is subdivided into different forms, including aggressive systemic mastocytosis (ASM) and systemic mastocytosis with an associated hematologic neoplasm (SM-AHN). Oncogenic genetic alterations include point mutations, mainly the KIT D816V, conferring poor prognosis and therapy resistance, and fusion genes, with those involving PDGFRA/PDGFRB as the most recurrent events. We here describe an ASM case negative to the KIT D816V and JAK2 V617F alterations but showing a RUNX1 frameshift heterozygous mutation and the co-occurrence of three fusion transcripts. The first one, PRKG2::PDGFRB, was generated by a balanced t(4;5)(q24;q32) translocation as the sole abnormality. Other two novel chimeras, KAT6A::NCOA2 and RXRA::NOTCH1, originated from cryptic intra-chromosomal abnormalities. The patient rapidly evolved towards SM-AHN, characterized by the persistence of the PRKG2::PDGFRB chimera, due to the presence of an extra copy of the der(5)t(4;5)(q24;q34) chromosome and an increase in the RUNX1 mutation allelic frequency. The results indicated that the transcriptional landscape and the mutational profile of SM deserve attention to predict the evolution and prognosis of this complex disease, whose classification criteria are still a matter of debate.

Locating saddle points using gradient extremals on manifolds adaptively revealed as point clouds.

Steady states are invaluable in the study of dynamical systems. High-dimensional dynamical systems, due to separation of time scales, often evolve toward a lower dimensional manifold M. We introduce an approach to locate saddle points (and other fixed points) that utilizes gradient extremals on such a priori unknown (Riemannian) manifolds, defined by adaptively sampled point clouds, with local coordinates discovered on-the-fly through manifold learning. The technique, which efficiently biases the dynamical system along a curve (as opposed to exhaustively exploring the state space), requires knowledge of a single minimum and the ability to sample around an arbitrary point. We demonstrate the effectiveness of the technique on the Müller-Brown potential mapped onto an unknown surface (namely, a sphere). Previous work employed a similar algorithmic framework to find saddle points using Newton trajectories and gentlest ascent dynamics; we, therefore, also offer a brief comparison with these methods.

miR-146a rs2431697 identifies myeloproliferative neoplasm patients with higher secondary myelofibrosis progression risk.

Myelofibrosis (MF) occurs as part of the natural history of polycythemia vera (PV) and essential thrombocythemia (ET), and remarkably shortens survival. Although JAK2V617F and CALR allele burden are the main transformation risk factors, inflammation plays a critical role by driving clonal expansion toward end-stage disease. NF-κB is a key mediator of inflammation-induced carcinogenesis. Here, we explored the involvement of miR-146a, a brake in NF-κB signaling, in MPN susceptibility and progression. rs2910164 and rs2431697, that affect miR-146a expression, were analyzed in 967 MPN (320 PV/333 ET/314 MF) patients and 600 controls. We found that rs2431697 TT genotype was associated with MF, particularly with post-PV/ET MF (HR = 1.5; p < 0.05). Among 232 PV/ET patients (follow-up time=8.5 years), 18 (7.8%) progressed to MF, being MF-free-survival shorter for rs2431697 TT than CC + CT patients (p = 0.01). Multivariate analysis identified TT genotype as independent predictor of MF progression. In addition, TT (vs. CC + CT) patients showed increased plasma inflammatory cytokines. Finally, miR-146a-/- mice showed significantly higher Stat3 activity with aging, parallel to the development of the MF-like phenotype. In conclusion, we demonstrated that rs2431697 TT genotype is an early predictor of MF progression independent of the JAK2V617F allele burden. Low levels of miR-146a contribute to the MF phenotype by increasing Stat3 signaling.

Two novel variants of the ABCG5 gene cause xanthelasmas and macrothrombocytopenia: a brief review of hematologic abnormalities of sitosterolemia.

Essentials Diagnosis of sitosterolemia, a rare recessive or syndromic disorder, is usually delayed. Peripheral blood smear is extremely useful for establishing the suspicion of sitosterolemia. High-throughput sequencing technology enables the molecular diagnosis of inherited thrombocytopenias. Accurate characterization of sitosterolemia helps us determine appropriate management. SUMMARY: Background Sitosterolemia (STSL) is a recessive inherited disorder caused by pathogenic variants in the ABCG5 and ABCG8 genes. Increased levels of plasma plant sterols (PSs) usually result in xanthomas and premature coronary atherosclerosis, although hematologic abnormalities may occasionally be present. This clinical picture is unfamiliar to many physicians, and patients may be at high risk of misdiagnosis. Objectives To report two novel ABCG5 variants causing STSL in a Spanish patient, and review the clinical and mutational landscape of STSL. Patient/Methods A 46-year-old female was referred to us with lifelong macrothrombocytopenia. She showed familial hypercholesterolemia-related xanthomas. Molecular analysis was performed with high-throughput sequencing. Plasma PS levels were evaluated with gas-liquid chromatography. The STSL landscape was reviewed with respect to specific online databases and all reports published since 1974. Results A blood smear revealed giant platelets and stomatocytes. Novel compound heterozygous variants were detected in exons 7 (c.914C>G) and 13 (c.1890delT) of ABCG5. The patient showed an increased plasma level of sitosterol. These findings support the diagnosis of STSL. In our review, we identified only 25 unrelated STLS patients who presented with hematologic abnormalities including macrothrombocytopenia. It remains unknown why only some patients develop hematologic abnormalities. Conclusions This is the first Spanish STSL patient to be reported and molecularly characterized. The early diagnosis of STLS is strongly supported by the presence of stomatocytes in blood smears. The definitive diagnosis of STSL by measurement of serum PS levels and molecular analyses prompted the use of ezetimibe therapy.

Uniparental disomy causes deficiencies of vitamin K-dependent proteins.

Essentials Vitamin K-dependent coagulant factor deficiency (VKCFD) is a rare autosomal recessive disorder. We describe a case of inherited VKCFD due to uniparental disomy. The homozygous mutation caused the absence of GGCX isoform 1 and overexpression of Δ2GGCX. Hepatic and non-hepatic vitamin K-dependent proteins must be assayed to monitor VKCFD treatment. SUMMARY: Background Inherited deficiency of all vitamin K-dependent coagulant factors (VKCFD) is a rare autosomal recessive disorder caused by mutations in the γ-glutamyl carboxylase gene (GGCX) or the vitamin K epoxide reductase gene (VKORC1), with great heterogeneity in terms of both clinical presentation and response to treatment. Objective To characterize the molecular basis of VKCFD in a Spanish family. Methods and Results Sequencing of candidate genes, comparative genomic hybridization and massive sequencing identified a new mechanism causing VKCFD in the proband. Uniparental disomy (UPD) of chromosome 2 caused homozygosity of a mutation (c.44-1G>A) resulting in aberrant GGCX splicing. This change contributed to absent expression of the mRNA coding for the full-length protein, and to four-fold overexpression of the smaller mRNA isoform lacking exon 2 (Δ2GGCX). Δ2GGCX might be responsible for two unexpected clinical observations in the patient: (i) increased plasma osteocalcin levels following vitamin K1 supplementation; and (ii) a mild non-bleeding phenotype. Conclusions Our study identifies a new autosomal disease, VKCFD1, caused by UPD. These data suggest that the Δ2GGCX isoform may retain enzymatic activity, and strongly encourage the evaluation of both hepatic and non-hepatic vitamin K-dependent proteins to assess differing responses to vitamin K supplementation in VKCFD patients.

Current opinion and consensus statement regarding the diagnosis, prognosis, and treatment of patients with essential thrombocythemia: a survey of the Spanish Group of Ph-negative Myeloproliferative Neoplasms (GEMFIN) using the Delphi method.

The current consensus on the diagnosis, prognosis, and treatment of essential thrombocythemia (ET) is based on experts' recommendations. However, several aspects of the diagnosis of, prognosis of, and therapy for ET are still controversial. The Delphi method was employed with an expert panel of members of the Spanish Group of Ph-negative Myeloproliferative Neoplasms in order to identify the degree of agreement on the diagnosis, prognosis, and treatment of ET. Nine leading experts selected a total of 41 clinical hematologists with well-known expertise in ET. An electronic questionnaire was used to collect the questions rated in a four-step scale. The questions were grouped into four blocks: diagnosis, risk stratification, goals of therapy, and treatment strategy. After the first round consisting of 80 questions, a second round including 14 additional questions focused on the recommendations advocated by experts of the European LeukemiaNet in 2011 was analyzed. The median and mean values for the first and second rounds were calculated. A summary of the conclusions considered as the most representative of each block of questions is presented. The Delphi method is a powerful instrument to address the current approaches and controversies surrounding ET.

Design and application of a 23-gene panel by next-generation sequencing for inherited coagulation bleeding disorders.

INTRODUCTION: Molecular testing of Inherited bleeding coagulation disorders (IBCDs) not only offers confirmation of diagnosis but also aids in genetic counselling, prenatal diagnosis and in certain cases genotype-phenotype correlations are important for predicting the clinical course of the disease and to allow tailor-made follow-up of individuals. Until recently, genotyping has been mainly performed by Sanger sequencing, a technique known to be time consuming and expensive. Currently, next-generation sequencing (NGS) offers a new potential approach that enables the simultaneous investigation of multiple genes at manageable cost. AIM: The aim of this study was to design and to analyse the applicability of a 23-gene NGS panel in the molecular diagnosis of patients with IBCDs. METHODS: A custom target enrichment library was designed to capture 31 genes known to be associated with IBCDs. Probes were generated for 296 targets to cover 86.3 kb regions (all exons and flanking regions) of these genes. Twenty patients with an IBCDs phenotype were studied using NGS technology. RESULTS: In all patients, our NGS approach detected causative mutations. Twenty-one pathogenic variants were found; while most of them were missense (18), three deletions were also identified. Six novel mutations affecting F8, FGA, F11, F10 and VWF genes, and 15 previously reported variants were detected. NGS and Sanger sequencing were 100% concordant. CONCLUSION: Our results demonstrate that this approach could be an accurate, reproducible and reliable tool in the rapid genetic diagnosis of IBCDs.

Mutation status and immunoglobulin gene rearrangements in patients from northwest and central region of Spain with chronic lymphocytic leukemia.

The aim of this study was to investigate the frequency and mutation status of the immunoglobulin heavy variable chain (IGHV) in a cohort of 224 patients from northwest and central region of Spain diagnosed with chronic lymphocytic leukemia (CLL), and to correlate it with cytogenetic abnormalities, overall survival (OS) and time to first treatment (TTFT). 125 patients had mutated IGHV, while 99 had unmutated IGHV. The most frequently used IGHV family was IGHV3, followed by IGHV1 and IGHV4. The regions IGHV3-30, IGHV1-69, IGHV3-23, and IGHV4-34 were the most commonly used. Only 3.1% of the patients belonged to the subfamily IGHV3-21 and we failed to demonstrate a worse clinical outcome in this subgroup. The IGHV4 family appeared more frequently with mutated pattern, similar to IGHV3-23 and IGHV3-74. By contrast, IGHV1-69 was expressed at a higher frequency in unmutated CLL patients. All the cases from IGHV3-11 and almost all from IGHV5-51 subfamily belonged to the group of unmutated CLL.

Exome sequencing reveals novel and recurrent mutations with clinical impact in blastic plasmacytoid dendritic cell neoplasm.

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a very rare disease that currently lacks genomic and genetic biomarkers to assist in its clinical management. We performed whole-exome sequencing (WES) of three BPDCN cases. Based on these data, we designed a resequencing approach to identify mutations in 38 selected genes in 25 BPDCN samples. WES revealed 37-99 deleterious gene mutations per exome with no common affected genes between patients, but with clear overlap in terms of molecular and disease pathways (hematological and dermatological disease). We identified for the first time deleterious mutations in IKZF3, HOXB9, UBE2G2 and ZEB2 in human leukemia. Target sequencing identified 29 recurring genes, ranging in prevalence from 36% for previously known genes, such as TET2, to 12-16% for newly identified genes, such as IKZF3 or ZEB2. Half of the tumors had mutations affecting either the DNA methylation or chromatin remodeling pathways. The clinical analysis revealed that patients with mutations in DNA methylation pathway had a significantly reduced overall survival (P=0.047). We provide the first mutational profiling of BPDCN. The data support the current WHO classification of the disease as a myeloid disorder and provide a biological rationale for the incorporation of epigenetic therapies for its treatment.

Genome-wide profiling of methylation identifies novel targets with aberrant hypermethylation and reduced expression in low-risk myelodysplastic syndromes.

Gene expression profiling signatures may be used to classify the subtypes of Myelodysplastic syndrome (MDS) patients. However, there are few reports on the global methylation status in MDS. The integration of genome-wide epigenetic regulatory marks with gene expression levels would provide additional information regarding the biological differences between MDS and healthy controls. Gene expression and methylation status were measured using high-density microarrays. A total of 552 differentially methylated CpG loci were identified as being present in low-risk MDS; hypermethylated genes were more frequent than hypomethylated genes. In addition, mRNA expression profiling identified 1005 genes that significantly differed between low-risk MDS and the control group. Integrative analysis of the epigenetic and expression profiles revealed that 66.7% of the hypermethylated genes were underexpressed in low-risk MDS cases. Gene network analysis revealed molecular mechanisms associated with the low-risk MDS group, including altered apoptosis pathways. The two key apoptotic genes BCL2 and ETS1 were identified as silenced genes. In addition, the immune response and micro RNA biogenesis were affected by the hypermethylation and underexpression of IL27RA and DICER1. Our integrative analysis revealed that aberrant epigenetic regulation is a hallmark of low-risk MDS patients and could have a central role in these diseases.

SNP-based mapping arrays reveal high genomic complexity in monoclonal gammopathies, from MGUS to myeloma status.

Genetic events mediating transformation from premalignant monoclonal gammopathies (MG) to multiple myeloma (MM) are unknown. To obtain a comprehensive genomic profile of MG from the early to late stages, we performed high-resolution analysis of purified plasma cells from 20 MGUS, 20 smoldering MM (SMM) and 34 MM by high-density 6.0 SNP array. A progressive increase in the incidence of copy number abnormalities (CNA) from MGUS to SMM and to MM (median 5, 7.5 and 12 per case, respectively) was observed (P=0.006). Gains on 1q, 3p, 6p, 9p, 11q, 19p, 19q and 21q along with 1p, 16q and 22q deletions were significantly less frequent in MGUS than in MM. Although 11q and 21q gains together with 16q and 22q deletions were apparently exclusive of MM status, we observed that these abnormalities were also present in minor subclones in MGUS. Overall, a total of 65 copy number-neutral LOH (CNN-LOH) were detected. Their frequency was higher in active MM than in the asymptomatic entities (P=0.047). A strong association between genetic lesions and fragile sites was also detected. In summary, our study shows an increasing genomic complexity from MGUS to MM and identifies new chromosomal regions involved in CNA and CNN-LOH.

Gene expression profiling in MDS and AML: potential and future avenues.

Today, the classification systems for myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML) already incorporate cytogenetic and molecular genetic aberrations in an attempt to better reflect disease biology. However, in many MDS/AML patients no genetic aberrations have been identified yet, and even within some cytogenetically well-defined subclasses there is considerable clinical heterogeneity. Recent advances in genomics technologies such as gene expression profiling (GEP) provide powerful tools to further characterize myeloid malignancies at the molecular level, with the goal to refine the MDS/AML classification system, incorporating as yet unknown molecular genetic and epigenetic pathomechanisms, which are likely reflected by aberrant gene expression patterns. In this study, we provide a comprehensive review on how GEP has contributed to a refined molecular taxonomy of MDS and AML with regard to diagnosis, prediction of clinical outcome, discovery of novel subclasses and identification of novel therapeutic targets and novel drugs. As many challenges remain ahead, we discuss the pitfalls of this technology and its potential including future integrative studies with other genomics technologies, which will continue to improve our understanding of malignant transformation in myeloid malignancies and thereby contribute to individualized risk-adapted treatment strategies for MDS and AML patients.

Both expanded and uncultured mesenchymal stem cells from MDS patients are genomically abnormal, showing a specific genetic profile for the 5q- syndrome.

The presence of cytogenetic aberrations on mesenchymal stem cells (MSC) from myelodysplastic syndrome (MDS) patients is controversial. The aim of the study is to characterize bone marrow (BM) derived MSC from patients with MDS using: kinetic studies, immunophenotyping, fluorescent in situ hybridization (FISH) and genetic changes by array-based comparative genomic hybridization (array-CGH). In all 36 cases of untreated MDS were studied. MDS-MSC achieved confluence at a significantly slower rate than donor-MSC, and the antigenic expression of CD105 and CD104 was lower. Array-CGH studies showed DNA genomic changes that were proved not to be somatic. These results were confirmed by FISH. To confirm that genomic changes were also present in freshly obtained MSCs they were enriched by sorting BM cells with the following phenotype: CD45(-)/CD73(++)/CD34(-)/CD271(++). They also showed genomic changes that were confirmed by FISH. To analyze the relationship of these aberrations with clinical-biological data an unsupervised hierarchical cluster analysis was performed, two clusters were identified: the first one included the 5q- syndrome patients, whereas the other incorporated other MDS. Our results show, for the first time that MSC from MDS display genomic aberrations, assessed by array-CGH and FISH, some of them specially linked to a particular MDS subtype, the 5q- syndrome.

The immunophenotype of different immature, myeloid and B-cell lineage-committed CD34+ hematopoietic cells allows discrimination between normal/reactive and myelodysplastic syndrome precursors.

Occurrence of phenotypic abnormalities in CD34(+) hematopoietic progenitor and precursor cells (HPC) and their major B-cell and nonlymphoid compartments has been frequently reported in myelodysplastic syndromes (MDS). Here, we analyze for the first time the numerical and phenotypic abnormalities of different maturation-associated subsets of bone marrow (BM) CD34(+) HPC from 50 newly diagnosed MDS patients in comparison to normal/reactive BM (n=29). Our results confirm the existence of heterogeneously altered phenotypes among CD34(+) HPC from MDS and indicate that such variability depends both on the relative distribution of the different subsets of CD34(+) HPC committed into the different myeloid and B-lymphoid compartments, and their immunophenotype (for example, higher reactivity for CD117 and CD13 and lower expression of CyMPO, CD64 and CD65 on CD34(+) immature and neutrophil precursors), a clear association existing between the accumulation of CD34(+) HPC and that of immature CD34(+) HPC. Interestingly, expansion of erythroid- and neutrophil-lineage CD34(+) cells is detected in low-grade MDS at the expense of CD34(+) plasmacytoid dendritic cell and B-cell precursors, while expansion of immature CD34(+) precursors occurs in high-grade MDS. On the basis of the number and severity of the phenotypic abnormalities detected, a scoring system is proposed that efficiently discriminates between normal/reactive and MDS CD34(+) HPC, the mean score significantly increasing from low- to high-grade MDS.

6q deletion in Waldenström macroglobulinemia is associated with features of adverse prognosis.

Fluorescence in situ hybridisation (FISH) is an effective technique for the cytogenetic analysis of Waldenström macroglobulinemia (WM), but the potential impact of molecular cytogenetics on disease evolution and as a prognostic marker is still unknown. Deletion of the long arm of chromosome 6 (6q-) is the most frequent cytogenetic abnormality in WM. This study analysed the prevalence of this aberration in 102 WM patients, and correlated it with disease characteristics. The incidence of 6q21 deletion was 7% by conventional cytogenetics and 34% when analysed by FISH (54% when cytoplasmic immunoglobulin M-FISH was used). Patients with deletion of 6q displayed features of adverse prognosis, such as higher levels of beta2-microglobulin and monoclonal paraprotein and a greater tendency to display anaemia and hypoalbuminemia. Interestingly, there was a correlation between the presence of 6q deletion and the International Staging System prognostic index (incidence of 6q- among patients stratified in stages 1, 2 and 3 was 24%, 42% and 67% respectively). Those patients diagnosed with smouldering WM who displayed the abnormality showed a trend to an earlier requirement of treatment. Finally, the survival analysis did not show differences between the two groups of patients, probably due to the short follow up of our series.

Long-term results of laparoscopic extended surgery in advanced gastric cancer: a series of 101 patients.

BACKGROUND/AIMS: The objective of our paper is to report on the remote results of patients with gastric cancer treated by mini-invasive surgery as a surgical tool with the "intention to treat with laparoscopy". METHODOLOGY: Between June 1993 and January 2004, 101 patients comprising 72 men and 29 women with gastric adenocarcinoma were prospectively selected by two hospitals based on prior agreement (the CHU Charleroi, Belgium, and Zumárraga Hospital, the Basque Country, Spain). Patients with adenocarcinoma of the cardia were excluded. Average age of the patients was 67 (37-83). RESULTS: Postoperative mortality within 60 days of operation was of 5 patients; 87 patients were therefore properly followed-up for an average of 41 months (7-129). Average survival time for 10 non-resected patients was 4.5 months. Average survival rate of the 10 palliatively resected patients was 7.1 months. Actuarial 5-year survival rate RO-type surgery was 34%. The global actuarial 5-year survival rate after resective surgery was 29%. CONCLUSIONS: Laparoscopic gastrectomy with any kind of lymphadenectomy is a heavy but safe operation, and produces acceptable mortality and morbidity rates in patients with advanced gastric cancer in a general poor condition. Laparoscopic gastrectomies for locally advanced cancers are equivalent to those reported by laparotomy as far as long-term oncological results are concerned.