Sparse deconvolution in one and two dimensions: applications in endocrinology and single-molecule fluorescence imaging.

Deconvolution of noisy signals is an important task in analytical chemistry, examples being spectral deconvolution or deconvolution in microscopy. When the number of spectral peaks or single emitters in imaging is limited, the solution of the deconvolution is required to be sparse, and desirable results are obtained using a penalized estimation techniques. We impose sparseness by using penalized regression with a penalty based on the L0-norm, as discussed in earlier work. Several extensions to this approach are presented. Results are demonstrated on pulse identification in endocrine data where the aim is to model the secretion pattern as a sparse series of spikes. An application in single-molecule fluorescence imaging demonstrates the algorithm when applied to two-dimensional data.

Immature MEF2C-dysregulated T-cell leukemia patients have an early T-cell precursor acute lymphoblastic leukemia gene signature and typically have non-rearranged T-cell receptors.

Three distinct immature T-cell acute lymphoblastic leukemia entities have been described including cases that express an early T-cell precursor immunophenotype or expression profile, immature MEF2C-dysregulated T-cell acute lymphoblastic leukemia cluster cases based on gene expression analysis (immature cluster) and cases that retain non-rearranged TRG@ loci. Early T-cell precursor acute lymphoblastic leukemia cases exclusively overlap with immature cluster samples based on the expression of early T-cell precursor acute lymphoblastic leukemia signature genes, indicating that both are featuring a single disease entity. Patients lacking TRG@ rearrangements represent only 40% of immature cluster cases, but no further evidence was found to suggest that cases with absence of bi-allelic TRG@ deletions reflect a distinct and even more immature disease entity. Immature cluster/early T-cell precursor acute lymphoblastic leukemia cases are strongly enriched for genes expressed in hematopoietic stem cells as well as genes expressed in normal early thymocyte progenitor or double negative-2A T-cell subsets. Identification of early T-cell precursor acute lymphoblastic leukemia cases solely by defined immunophenotypic criteria strongly underestimates the number of cases that have a corresponding gene signature. However, early T-cell precursor acute lymphoblastic leukemia samples correlate best with a CD1 negative, CD4 and CD8 double negative immunophenotype with expression of CD34 and/or myeloid markers CD13 or CD33. Unlike various other studies, immature cluster/early T-cell precursor acute lymphoblastic leukemia patients treated on the COALL-97 protocol did not have an overall inferior outcome, and demonstrated equal sensitivity levels to most conventional therapeutic drugs compared to other pediatric T-cell acute lymphoblastic leukemia patients.

Molecular subtypes of glioma identified by genome-wide methylation profiling.

Recent studies have indicated a prognostic role for genome-wide methylation in gliomas: Tumors that show an overall increase in DNA methylation at CpG sites (CIMP+; CpG island methylator phenotype) have a more favorable prognosis than CIMP- gliomas. Here, we have determined whether methylation profiling can identify more and clinically relevant molecular subtypes of glioma by performing genome-wide methylation profiling on 138 glial brain tumors of all histological diagnosis. Hopach (Hierarchical ordered partitioning and collapsing hybrid) clustering using the 1,000 most variable CpGs identified three distinct glioma subtypes (C+(1p19q), C+(wt), and C-) and one adult brain subtype. All "C+(1p19q) " and "C+(wt)" tumors were CIMP+ whereas most (50/54) "C-" tumors were CIMP-. The C- subtype gliomas contained many glioblastomas and all pilocytic astrocytomas. 1p19q LOH was frequent in the C+(1p19q) subtype. Other genetic changes (IDH1 mutation and EGFR amplification) and gene-expression based molecular subtypes also segregated in distinct methylation subtypes, demonstrating that these subtypes are also genetically distinct. Each subtype was associated with its own prognosis: median survival for C-, C+(1p19q), and C+(wt) tumors was 1.18, 5.00, and 2.62 years, respectively. The prognostic value of these methylation subtypes was validated on an external dataset from the TCGA. Analysis of recurrences of 14 primary tumors samples indicates that shifts between some C+(wt) and C+(1p/19q) tumors can occur between the primary and recurrent tumor, but CIMP status remained stable. Our data demonstrate that methylation profiling identifies at least three prognostically relevant subtypes of glioma that can aid diagnosis and potentially guide treatment for patients.

Spatially Regularized Shape Analysis of the Hippocampus Using P-Spline Based Shape Regression.

Shape analysis is increasingly becoming important to study changes in brain structures in relation to clinical neurological outcomes. This is a challenging task due to the high dimensionality of shape representations and the often limited number of available shapes. Current techniques counter the poor ratio between dimensions and sample size by using regularization in shape space, but do not take into account the spatial relations within the shapes. This can lead to models that are biologically implausible and difficult to interpret. We propose to use P-spline based regression, which combines a generalized linear model (GLM) with the coefficients described as B-splines and a penalty term that constrains the regression coefficients to be spatially smooth. Owing to the GLM, this method can naturally predict both continuous and discrete outcomes and can include non-spatial covariates without penalization. We evaluated our method on hippocampus shapes extracted from magnetic resonance (MR) images of 510 non-demented, elderly people. We related the hippocampal shape to age, memory score, and sex. The proposed method retained the good performance of current techniques, such as ridge regression, but produced smoother coefficient fields that are easier to interpret.

Data representations and -analyses of binary diary data in pursuit of stratifying children based on common childhood illnesses.

In this article we analyse diary reports concerning childhood symptoms of illness, these data are part of a larger study with other types of measurements on childhood asthma. The children are followed for three years and the diaries are updated, by the parents, on a daily basis. Here we focus on the methodological implications of analysing such data. We investigate two ways of representing the data and explore which tools are applicable given both representations. The first representation relies on proper alignment and point by point comparison of the signals. The second approach takes into account combinations of symptoms on a day by day basis and boils down to the analysis of counts. In the present case both methods are well applicable. However, more generally, when symptom episodes are occurring more at random locations in time, a point by point comparison becomes less applicable and shape based approaches will fail to come up with satisfactory results. In such cases, pattern based methods will be of much greater use. The pattern based representation focuses on reoccurring patterns and ignores ordering in time. With this representation we stratify the data on the level of years, so that possibly yearly differences can still be detected.

Sparse deconvolution of high-density super-resolution images.

In wide-field super-resolution microscopy, investigating the nanoscale structure of cellular processes, and resolving fast dynamics and morphological changes in cells requires algorithms capable of working with a high-density of emissive fluorophores. Current deconvolution algorithms estimate fluorophore density by using representations of the signal that promote sparsity of the super-resolution images via an L1-norm penalty. This penalty imposes a restriction on the sum of absolute values of the estimates of emitter brightness. By implementing an L0-norm penalty--on the number of fluorophores rather than on their overall brightness--we present a penalized regression approach that can work at high-density and allows fast super-resolution imaging. We validated our approach on simulated images with densities up to 15 emitters per µm(-2) and investigated total internal reflection fluorescence (TIRF) data of mitochondria in a HEK293-T cell labeled with DAKAP-Dronpa. We demonstrated super-resolution imaging of the dynamics with a resolution down to 55 nm and a 0.5 s time sampling.

Volume measurements on three-dimensional photogrammetry after extended strip versus total cranial remodeling for sagittal synostosis: A comparative cohort study.

BACKGROUND: Surgery for sagittal synostosis aims at correction of skull shape and restoration of growth potential. Small cranial volume is associated with raised intracranial pressure (ICP). Although many techniques have been described, information on postoperative volume related to early and late remodeling is lacking. METHODS: Between 2004 and 2008, a total of 95 patients were collected who underwent either early extended strip craniectomy or late total cranial remodeling according to age of presentation. Volume was measured on three-dimensional (3D) photogrammetry. Volume measurements were related to cranial index (CI), head circumference (HCsd), and signs of raised ICP. In a small subset of patients, volume measurements on 3D photogrammetry were assessed for inter- and intrarater reliability and compared to 3D computed tomography (CT). RESULTS: Volume was increased in all patients before and after surgery compared to normative values. Postoperatively, late total cranial remodeling resulted in a slightly larger volume than early extended strip craniectomy. Volume measurements showed a good correlation with HCsd (0.67) and a poor relationship with CI (0.13). Headache occurred more frequently in patients with a lower cranial volume. Although papilledema and reoperation showed the same trend, the numbers were too small for statistical analysis. Reproducibility of volume measurements on 3D photogrammetry was high, as was the correlation with measurements on CT. CONCLUSION: Late total cranial remodeling results in a larger postoperative volume, as measured on 3D photogrammetry, than extended strip craniectomy. Clinical signs of raised ICP occur more frequently in patients with a smaller volume. To measure volume, 3D photogrammetry is a good alternative to CT.

The formation of the foramen magnum and its role in developing ventriculomegaly and Chiari I malformation in children with craniosynostosis syndromes.

OBJECT: Craniosynostosis syndromes are characterized by prematurely fused skull sutures, however, less is known about skull base synchondroses. This study evaluates how foramen magnum (FM) size, and closure of its intra-occipital synchondroses (IOS) differ between patients with different craniosynostosis syndromes and control subjects; and whether this correlates to ventriculomegaly and/or Chiari malformation type I (CMI), intracranial disturbances often described in these patients. METHODS: Surface area and anterior-posterior (A-P) diameter were measured in 175 3D-CT scans of 113 craniosynostosis patients, and in 53 controls (0-10 years old). Scans were aligned in a 3D multiplane-platform. The frontal and occipital horn ratio was used as an indicator of ventricular volume, and the occurrence of CMI was recorded. Synchondroses were studied in scans with a slice thickness ≤1.25 mm. A generalized linear mixed model and a repeated measures ordinal logistic regression model were used to study differences. RESULTS: At birth, patients with craniosynostosis syndromes have a smaller FM than controls (p < 0.05). This is not related to the presence of CMI (p = 0.36). In Crouzon-Pfeiffer patients the anterior and posterior IOS fused prematurely (p < 0.01), and in Apert patients only the posterior IOS fused prematurely (p = 0.028). CONCLUSION: The FM is smaller in patients with craniosynostosis syndromes than in controls, and is already smaller at birth. In addition to the timing of IOS closure, other factors may influence FM size.

The role of the posterior fossa in developing Chiari I malformation in children with craniosynostosis syndromes.

OBJECTIVE: Patients with craniosynostosis syndromes are at risk of increased intracranial pressure (ICP) and Chiari I malformation (CMI), caused by a combination of restricted skull growth, venous hypertension, obstructive sleep apnea (OSA), and an overproduction or insufficient resorption of cerebrospinal fluid. This study evaluates whether craniosynostosis patients with CMI have an imbalance between cerebellar volume (CV) and posterior fossa volume (PFV), that is, an overcrowded posterior fossa. METHODS: Volumes were measured in 3D-SPGR T1-weighted MR scans of 28 'not-operated' craniosynostosis patients (mean age: 4.0 years; range: 0-14), 85 'operated' craniosynostosis patients (mean age: 8.0 years; range: 1-18), and 34 control subjects (mean age: 5.4 years; range: 0-15). Volumes and CV/PFV ratios were compared between the operated and not-operated craniosynostosis patients, between the individual craniosynostosis syndromes and controls, and between craniosynostosis patients with and without CMI. Data were logarithmically transformed and studied with analysis of covariance (ANCOVA). RESULTS: The CV, PFV, and CV/PFV ratios of not-operated craniosynostosis patients and operated craniosynostosis patients were similar to those of the control subjects. None of the individual syndromes was associated with a restricted PFV. However, craniosynostosis patients with CMI had a significantly higher CV/PFV ratio than the control group (0.77 vs. 0.75; p = 0.008). The range of CV/PFV ratios for craniosynostosis patients with CMI, however, did not exceed the normal range. CONCLUSION: Volumes and CV/PFV ratio cannot predict which craniosynostosis patients are more prone to developing CMI than others. Treatment should focus on the skull vault and other contributing factors to increased ICP, including OSA and venous hypertension.

The prognostic value of extraocular extension in relation to monosomy 3 and gain of chromosome 8q in uveal melanoma.

PURPOSE: To identify the prognostic value of extraocular extension in enucleated uveal melanoma (UM) patients and to correlate extraocular extension to chromosomal aberrations, metastasis-free survival (MFS), and clinico-histopathological risk factors. METHODS: Retrospective study of patients with UM treated with enucleation between 1987 and 2011. Melanoma-related metastasis and death were recorded. Statistical analysis (log-rank test or Cox regression analysis) was performed to correlate MFS with tumor characteristics, extraocular extension, episcleral diameter of the extraocular extension, cell type, extracellular matrix patterns, inflammation, loss of chromosome 3, and gain of chromosome 8q. RESULTS: In 43 (12%) of 357 patients, extraocular extension was observed. In this subset of patients, we noted a reduced survival of 70 months (105.5 months, P = 0.010) compared with patients without extraocular extension (175.8 months). Patients with gain of chromosomal region 8q in UM with extraocular extension had an increased risk of metastatic disease (P < 0.001). In multivariate Cox proportional hazard analysis, largest basal tumor diameter (P = 0.001), extracellular matrix patterns (P = 0.009), episcleral diameter of the extraocular extension (P = 0.016), loss of chromosome 3 (P < 0.001), and gain of 8q (P < 0.001) were independent predictors for MFS. CONCLUSIONS: Larger episcleral diameter of the extraocular extension and additional gain of chromosome 8q in extraocular extension UM correlates to a worse prognosis. MFS is significantly reduced in UM with a large basal tumor diameter, extracellular matrix patterns, loss of chromosome 3, and gain of chromosome 8q.

Mixture models for two-dimensional baseline correction, applied to artifact elimination in time-resolved spectroscopy.

Baseline correction and artifact removal are important pre-processing steps in analytical chemistry. We propose a correction algorithm using a mixture model in combination with penalized regression. The model is an extension of a method recently introduced for baseline estimation in the case of one-dimensional data. The data are modeled as a smooth surface using tensor product P-splines. The weights of the P-splines regression model are computed from a mixture model where a datapoint is either allocated to the noise around the baseline, or to the artifact component. The method is broadly applicable for anisotropic smoothing of two-way data such as two-dimensional gel electrophoresis and two-dimensional chromatography data. We focus here on the application of the approach in femtosecond time-resolved spectroscopy, to eliminate strong artifact signals from the solvent.

Foramen magnum size and involvement of its intraoccipital synchondroses in Crouzon syndrome.

BACKGROUND: Cranial sutures and synchondroses tend to close prematurely in patients with Crouzon syndrome. This influences their skull vault and skull base development and may involve in common disturbances such as increased intracranial pressure and cerebellar tonsillar herniation. The authors' hypothesis was that Crouzon patients patients have a smaller foramen magnum than controls because of premature fusion of the intraoccipital synchondroses, putting them at risk for cerebellar tonsillar herniation. Therefore, foramen magnum size and time of intraoccipital synchondroses closure were evaluated and were related to the presence and degree of cerebellar tonsillar herniation. METHODS: The foramen magnum surface area and anteroposterior diameter were measured on three-dimensional computed tomographic scans of 27 Crouzon patients and 27 age-matched controls. Scans had a slice-thickness between 0.75 and 1.25 mm and were aligned in a three-dimensional reformatting platform. The t test was used to study size differences. Synchondroses were graded as described by Madeline and Elster and studied with ordinal logistic regression analysis. RESULTS: Crouzon patients had a smaller foramen magnum surface area (602 mm versus 767 mm, p < 0.001) and anteroposterior diameter (31 mm versus 35 mm, p < 0.001) compared with controls. Differences stayed constant over time. Intraoccipital synchondroses closed 3 to 9 months earlier in Crouzon patients than in controls (p < 0.05). CONCLUSIONS: Since intraoccipital synchondroses close earlier in Crouzon patients, from early life on their foramen magnum is smaller compared with controls. Within Crouzon patients, the presence of cerebellar tonsillar herniation could not be related to foramen magnum size. CLINICAL QUESTION/LEVEL OF EVIDENCE: Risk, II.

MGMT-STP27 methylation status as predictive marker for response to PCV in anaplastic Oligodendrogliomas and Oligoastrocytomas. A report from EORTC study 26951.

PURPOSE: The long-term follow-up results from the EORTC-26951 trial showed that the addition of procarbazine, CCNU, and vincristine (PCV) after radiotherapy increases survival in anaplastic oligodendrogliomas/oligoastrocytomas (AOD/AOA). However, some patients appeared to benefit more from PCV treatment than others. EXPERIMENTAL DESIGN: We conducted genome-wide methylation profiling of 115 samples included in the EORTC-26951 trial and extracted the CpG island hypermethylated phenotype (CIMP) and MGMT promoter methylation (MGMT-STP27) status. RESULTS: We first show that methylation profiling can be conducted on archival tissues with a performance that is similar to snap-frozen tissue samples. We then conducted methylation profiling on EORTC-26951 clinical trial samples. Univariate analysis indicated that CIMP+ or MGMT-STP27 methylated tumors had an improved survival compared with CIMP- and/or MGMT-STP27 unmethylated tumors [median overall survival (OS), 1.05 vs. 6.46 years and 1.06 vs. 3.8 years, both P < 0.0001 for CIMP and MGMT-STP27 status, respectively]. Multivariable analysis indicates that CIMP and MGMT-STP27 are significant prognostic factors for survival in presence of age, sex, performance score, and review diagnosis in the model. CIMP+ and MGMT-STP27 methylated tumors showed a clear benefit from adjuvant PCV chemotherapy: the median OS of CIMP+ samples in the RT and RT-PCV arms was 3.27 and 9.51 years, respectively (P = 0.0033); for MGMT-STP27 methylated samples, it was 1.98 and 8.65 years. There was no such benefit for CIMP- or for MGMT-STP27 unmethylated tumors. MGMT-STP27 status remained significant in an interaction test (P = 0.003). Statistical analysis of microarray (SAM) identified 259 novel CpGs associated with treatment response. CONCLUSIONS: MGMT-STP27 may be used to guide treatment decisions in this tumor type.

Intrinsic molecular subtypes of glioma are prognostic and predict benefit from adjuvant procarbazine, lomustine, and vincristine chemotherapy in combination with other prognostic factors in anaplastic oligodendroglial brain tumors: a report from EORTC study 26951.

PURPOSE: Intrinsic glioma subtypes (IGSs) are molecularly similar tumors that can be identified based on unsupervised gene expression analysis. Here, we have evaluated the clinical relevance of these subtypes within European Organisation for Research and Treatment of Cancer (EORTC) 26951, a randomized phase III clinical trial investigating adjuvant procarbazine, lomustine, and vincristine (PCV) chemotherapy in anaplastic oligodendroglial tumors. Our study includes gene expression profiles of formalin-fixed, paraffin-embedded (FFPE) clinical trial samples. PATIENTS AND METHODS: Gene expression profiling was performed in 140 samples, 47 fresh frozen samples and 93 FFPE samples, on HU133_Plus_2.0 and HuEx_1.0_st arrays, respectively. RESULTS: All previously identified six IGSs are present in EORTC 26951. This confirms that different molecular subtypes are present within a well-defined histologic subtype. Intrinsic subtypes are highly prognostic for overall survival (OS) and progression-free survival (PFS). They are prognostic for PFS independent of clinical (age, performance status, and tumor location), molecular (1p/19q loss of heterozygosity [LOH], IDH1 mutation, and MGMT methylation), and histologic parameters. Combining known molecular (1p/19q LOH, IDH1) prognostic parameters with intrinsic subtypes improves outcome prediction (proportion of explained variation, 30% v 23% for each individual group of factors). Specific genetic changes (IDH1, 1p/19q LOH, and EGFR amplification) segregate into different subtypes. We identified one subtype, IGS-9 (characterized by a high percentage of 1p/19q LOH and IDH1 mutations), that especially benefits from PCV chemotherapy. Median OS in this subtype was 5.5 years after radiotherapy (RT) alone versus 12.8 years after RT/PCV (P = .0349; hazard ratio, 2.18; 95% CI, 1.06 to 4.50). CONCLUSION: Intrinsic subtypes are highly prognostic in EORTC 26951 and improve outcome prediction when combined with other prognostic factors. Tumors assigned to IGS-9 benefit from adjuvant PCV.

Deconvolution of pulse trains with the L0 penalty.

The output of many instruments can be modeled as a convolution of an impulse response and a series of sharp spikes. Deconvolution considers the inverse problem: estimate the input spike train from an observed (noisy) output signal. We approach this task as a linear inverse problem, solved using penalized regression. We propose the use of an L(0) penalty and compare it with the more common L(2) and L(1) penalties. In all cases a simple and iterative weighted regression procedure can be used. The model is extended with a smooth component to handle drifting baselines. Application to three different data sets shows excellent results.

Integrated transcript and genome analyses reveal NKX2-1 and MEF2C as potential oncogenes in T cell acute lymphoblastic leukemia.

To identify oncogenic pathways in T cell acute lymphoblastic leukemia (T-ALL), we combined expression profiling of 117 pediatric patient samples and detailed molecular-cytogenetic analyses including the Chromosome Conformation Capture on Chip (4C) method. Two T-ALL subtypes were identified that lacked rearrangements of known oncogenes. One subtype associated with cortical arrest, expression of cell cycle genes, and ectopic NKX2-1 or NKX2-2 expression for which rearrangements were identified. The second subtype associated with immature T cell development and high expression of the MEF2C transcription factor as consequence of rearrangements of MEF2C, transcription factors that target MEF2C, or MEF2C-associated cofactors. We propose NKX2-1, NKX2-2, and MEF2C as T-ALL oncogenes that are activated by various rearrangements.

Intrinsic gene expression profiles of gliomas are a better predictor of survival than histology.

Gliomas are the most common primary brain tumors with heterogeneous morphology and variable prognosis. Treatment decisions in patients rely mainly on histologic classification and clinical parameters. However, differences between histologic subclasses and grades are subtle, and classifying gliomas is subject to a large interobserver variability. To improve current classification standards, we have performed gene expression profiling on a large cohort of glioma samples of all histologic subtypes and grades. We identified seven distinct molecular subgroups that correlate with survival. These include two favorable prognostic subgroups (median survival, >4.7 years), two with intermediate prognosis (median survival, 1-4 years), two with poor prognosis (median survival, <1 year), and one control group. The intrinsic molecular subtypes of glioma are different from histologic subgroups and correlate better to patient survival. The prognostic value of molecular subgroups was validated on five independent sample cohorts (The Cancer Genome Atlas, Repository for Molecular Brain Neoplasia Data, GSE12907, GSE4271, and Li and colleagues). The power of intrinsic subtyping is shown by its ability to identify a subset of prognostically favorable tumors within an external data set that contains only histologically confirmed glioblastomas (GBM). Specific genetic changes (epidermal growth factor receptor amplification, IDH1 mutation, and 1p/19q loss of heterozygosity) segregate in distinct molecular subgroups. We identified a subgroup with molecular features associated with secondary GBM, suggesting that different genetic changes drive gene expression profiles. Finally, we assessed response to treatment in molecular subgroups. Our data provide compelling evidence that expression profiling is a more accurate and objective method to classify gliomas than histologic classification. Molecular classification therefore may aid diagnosis and can guide clinical decision making.