Paper 5: Surveillance of multiple congenital anomalies: implementation of a computer algorithm in European registers for classification of cases.

BACKGROUND: Surveillance of multiple congenital anomalies is considered to be more sensitive for the detection of new teratogens than surveillance of all or isolated congenital anomalies. Current literature proposes the manual review of all cases for classification into isolated or multiple congenital anomalies. METHODS: Multiple anomalies were defined as two or more major congenital anomalies, excluding sequences and syndromes. A computer algorithm for classification of major congenital anomaly cases in the EUROCAT database according to International Classification of Diseases (ICD)v10 codes was programmed, further developed, and implemented for 1 year's data (2004) from 25 registries. The group of cases classified with potential multiple congenital anomalies were manually reviewed by three geneticists to reach a final agreement of classification as "multiple congenital anomaly" cases. RESULTS: A total of 17,733 cases with major congenital anomalies were reported giving an overall prevalence of major congenital anomalies at 2.17%. The computer algorithm classified 10.5% of all cases as "potentially multiple congenital anomalies". After manual review of these cases, 7% were agreed to have true multiple congenital anomalies. Furthermore, the algorithm classified 15% of all cases as having chromosomal anomalies, 2% as monogenic syndromes, and 76% as isolated congenital anomalies. The proportion of multiple anomalies varies by congenital anomaly subgroup with up to 35% of cases with bilateral renal agenesis. CONCLUSIONS: The implementation of the EUROCAT computer algorithm is a feasible, efficient, and transparent way to improve classification of congenital anomalies for surveillance and research.

Fuzzy logic algorithm to extract specific interaction forces from atomic force microscopy data

The atomic force microscope is not only a very convenient tool for studying the topography of different samples, but it can also be used to measure specific binding forces between molecules. For this purpose, one type of molecule is attached to the tip and the other one to the substrate. Approaching the tip to the substrate allows the molecules to bind together. Retracting the tip breaks the newly formed bond. The rupture of a specific bond appears in the force-distance curves as a spike from which the binding force can be deduced. In this article we present an algorithm to automatically process force-distance curves in order to obtain bond strength histograms. The algorithm is based on a fuzzy logic approach that permits an evaluation of "quality" for every event and makes the detection procedure much faster compared to a manual selection. In this article, the software has been applied to measure the binding strength between tubuline and microtubuline associated proteins.

Correlation of CXCL10, tumor necrosis factor receptor type II, and galectin 9 with disease activity in juvenile dermatomyositis.

OBJECTIVE: Juvenile dermatomyositis (DM) is a systemic autoimmune disorder of unknown immunopathogenesis in which the immune system targets the microvasculature of skeletal muscles, skin, and other organs. The current mainstay of therapy is a steroid regimen in combination with other immunosuppressive treatments. To date, no validated markers for monitoring disease activity have been identified, which hampers personalized treatment. This study was undertaken to identify a panel of proteins specifically related to active disease in juvenile DM. METHODS: We performed a multiplex immunoassay for plasma levels of 45 proteins related to inflammation in 25 patients with juvenile DM in 4 clinically well-defined groups, as determined by clinical activity and treatment. We compared them to 14 age-matched healthy children and 8 age-matched children with nonautoimmune muscle disease. RESULTS: Cluster analysis of circulating proteins showed distinct profiles for juvenile DM patients and controls based on a group of 10 proteins. In addition to CXCL10, tumor necrosis factor receptor type II (TNFRII) and galectin 9 were significantly increased in active juvenile DM. The levels of these 3 proteins were tightly linked to active disease and correlated with clinical scores (as measured by the Childhood Myositis Assessment Scale and physician's global assessment of disease activity on a visual analog scale). CONCLUSION: Our findings indicate that CXCL10, TNFRII, and galectin 9 correspond to disease status in juvenile DM and thus could be helpful in monitoring disease activity and guiding treatment. Furthermore, they might provide new knowledge about the pathogenesis of this autoimmune disease.

Diagnostic algorithm for a validated displacement grading of pediatric supracondylar fractures.

Background: The first AO comprehensive pediatric long bone fracture classification system has been established following a structured path of development and validation with experienced pediatric surgeons. Methods: A follow-up series of agreement studies was applied to specify and evaluate a grading system for displacement of pediatric supracondylar fractures. An iterative process comprising an international group of 5 experienced pediatric surgeons (Phase 1) followed by a pragmatic multicenter agreement study involving 26 raters (Phase 2) was used. The last evaluations were conducted on a consecutive collection of 154 supracondylar fractures documented by standard anteroposterior and lateral radiographs. Results: Fractures were classified according to 1 of 4 grades: I = incomplete fracture with no or minimal displacement; II = Incomplete fracture with continuity of the posterior (extension fracture) or anterior cortex (flexion fracture); III = lack of bone continuity (broken cortex), but still some contact between the fracture planes; IV = complete fracture with no bone continuity (broken cortex), and no contact between the fracture planes. A diagnostic algorithm to support the practical application of the grading system in a clinical setting, as well as an aid using a circle placed over the capitellum was proposed. The overall kappa coefficients were 0.68 and 0.61 in the Phase 1 and Phase 2 studies, respectively. In the Phase 1 study, fracture grades I, II, III, and IV were classified with median accuracies of 91%, 82%, 83%, and 99.5%, respectively. Similar median accuracies of 86% (Grade I), 73% (Grade II), 83%(Grade III), and 92% were reported for the Phase 2 study. Reliability was high in distinguishing complete, unstable fractures from stable injuries [ie, kappa coefficients of 0.84 (Phase 1) and 0.83 (Phase 2) were calculated]; in Phase 2, surgeons' accuracies in classifying complete fractures were all above 85%. Conclusions: With clear and unambiguous definition, this new grading system for supracondylar fracture displacement has proved to be sufficiently reliable and accurate when applied by pediatric surgeons in the framework of clinical routine as well as research.

An efficient algorithm to perform multiple testing in epistasis screening

Background: Research in epistasis or gene-gene interaction detection for human complex traits has grown over the last few years. It has been marked by promising methodological developments, improved translation efforts of statistical epistasis to biological epistasis and attempts to integrate different omics information sources into the epistasis screening to enhance power. The quest for gene-gene interactions poses severe multiple-testing problems. In this context, the maxT algorithm is one technique to control the false-positive rate. However, the memory needed by this algorithm rises linearly with the amount of hypothesis tests. Gene-gene interaction studies will require a memory proportional to the squared number of SNPs. A genome-wide epistasis search would therefore require terabytes of memory. Hence, cache problems are likely to occur, increasing the computation time. In this work we present a new version of maxT, requiring an amount of memory independent from the number of genetic effects to be investigated. This algorithm was implemented in C++ in our epistasis screening software MBMDR-3.0.3. We evaluate the new implementation in terms of memory efficiency and speed using simulated data. The software is illustrated on real-life data for Crohn’s disease. Results: In the case of a binary (affected/unaffected) trait, the parallel workflow of MBMDR-3.0.3 analyzes all gene-gene interactions with a dataset of 100,000 SNPs typed on 1000 individuals within 4 days and 9 hours, using 999 permutations of the trait to assess statistical significance, on a cluster composed of 10 blades, containing each four Quad-Core AMD Opteron(tm) Processor 2352 2.1 GHz. In the case of a continuous trait, a similar run takes 9 days. Our program found 14 SNP-SNP interactions with a multiple-testing corrected p-value of less than 0.05 on real-life Crohn’s disease (CD) data. Conclusions: Our software is the first implementation of the MB-MDR methodology able to solve large-scale SNP-SNP interactions problems within a few days, without using much memory, while adequately controlling the type I error rates. A new implementation to reach genome-wide epistasis screening is under construction. In the context of Crohn’s disease, MBMDR-3.0.3 could identify epistasis involving regions that are well known in the field and could be explained from a biological point of view. This demonstrates the power of our software to find relevant phenotype-genotype higher-order associations.

A Canonical Genetic Algorithm for Blind Inversion of Linear Channels

It is well known the relationship between source separation and blind deconvolution: If a filtered version of an unknown i.i.d. signal is observed, temporal independence between samples can be used to retrieve the original signal, in the same manner as spatial independence is used for source separation. In this paper we propose the use of a Genetic Algorithm (GA) to blindly invert linear channels. The use of GA is justified in the case of small number of samples, where other gradient-like methods fails because of poor estimation of statistics.

Schizotypal perceptual aberrations of time: correlation between score, behavior and brain activity.

A fundamental trait of the human self is its continuum experience of space and time. Perceptual aberrations of this spatial and temporal continuity is a major characteristic of schizophrenia spectrum disturbances--including schizophrenia, schizotypal personality disorder and schizotypy. We have previously found the classical Perceptual Aberration Scale (PAS) scores, related to body and space, to be positively correlated with both behavior and temporo-parietal activation in healthy participants performing a task involving self-projection in space. However, not much is known about the relationship between temporal perceptual aberration, behavior and brain activity. To this aim, we composed a temporal Perceptual Aberration Scale (tPAS) similar to the traditional PAS. Testing on 170 participants suggested similar performance for PAS and tPAS. We then correlated tPAS and PAS scores to participants' performance and neural activity in a task of self-projection in time. tPAS scores correlated positively with reaction times across task conditions, as did PAS scores. Evoked potential mapping and electrical neuroimaging showed self-projection in time to recruit a network of brain regions at the left anterior temporal cortex, right temporo-parietal junction, and occipito-temporal cortex, and duration of activation in this network positively correlated with tPAS and PAS scores. These data demonstrate that schizotypal perceptual aberrations of both time and space, as reflected by tPAS and PAS scores, are positively correlated with performance and brain activation during self-projection in time in healthy individuals along the schizophrenia spectrum.

Mineralocorticoid effects in the kidney: correlation between alphaENaC, GILZ, and Sgk-1 mRNA expression and urinary excretion of Na+ and K+.

Aldosterone exerts its effects through interactions with two types of binding sites, the mineralocorticoid (MR) and the glucocorticoid (GR) receptors. Although both receptors are known to be involved in the anti-natriuretic response to aldosterone, the mechanisms of signal transduction leading to modulation of electrolyte transport are not yet fully understood. This study measured the Na(+) and K(+) urinary excretion and the mRNA levels of three known aldosterone-induced transcripts, the serum and glucocorticoid-induced kinase (Sgk-1), the alpha subunit of the epithelial Na(+) channel (alphaENaC), and the glucocorticoid-induced-leucine-zipper protein (GILZ) in the whole kidney and in isolated cortical collecting tubules of adrenalectomized rats treated with low doses of aldosterone and/or dexamethasone. The resulting plasma concentrations of both steroids were close to 1 nmol/L. Aldosterone, given with or without dexamethasone, induced anti-natriuresis and kaliuresis, whereas dexamethasone alone did not. GILZ and alphaENaC transcripts were higher after treatment with either or both hormones, whereas the mRNA abundance of Sgk-1 was increased in the cortical collecting tubule by aldosterone but not by dexamethasone. We conclude the increased expression of Sgk-1 in the cortical collecting tubules is a primary event in the early antinatriuretic and kaliuretic responses to physiologic concentrations of aldosterone. Induction of alphaENaC and/or GILZ mRNAs may play a permissive role in the enhancement of the early and/or late responses; these effects may be necessary for a full response but do not by themselves promote early changes in urinary Na(+) and K(+) excretion.

Ultrasound evaluation of synovitis in RA: Correlation with clinical disease activity and sensitivity to change in an observational cohort study.

OBJECTIVE: To evaluate the correlation between clinical measures of disease activity and a ultrasound (US) scoring system for synovitis applied by many different ultrasonographers in a daily routine care setting within the Swiss registry for RA (SCQM) and further to determine the sensitivity to change of this US Score. METHODS: One hundred and eight Swiss rheumatologists were trained in performing the Swiss Sonography in Arthritis and Rheumatism (SONAR) score. US B-mode and Power Doppler (PwD) scores were correlated with DAS28 and compared between the clinical categories in a cross-sectional cohort of patients. In patients with a second US (longitudinal cohort), we investigated if change in US score correlated with change in DAS and evaluated the responsiveness of both methods. RESULTS: In the cross-sectional cohort with 536 patients, correlation between the B-mode score and DAS28 was significant but modest (Pearson coefficient r=0.41, P<0.0001). The same was true for the PwD score (r=0.41, P<0.0001). In the longitudinal cohort with 183 patients we also found a significant correlation between change in B-mode and in PwD score with change in DAS28 (r=0.54, P<0.0001 and r=0.46, P<0.0001, respectively). Both methods of evaluation (DAS and US) showed similar responsiveness according to standardized response mean (SRM). CONCLUSIONS: The SONAR Score is practicable and was applied by many rheumatologists in daily routine care after initial training. It demonstrates significant correlations with the degree of as well as change in disease activity as measured by DAS. On the level of the individual, the US score shows many discrepancies and overlapping results exist.

Immunohistochemistry as a valuable tool to assess CD30 expression in peripheral T-cell lymphomas: high correlation with mRNA levels.

The extended use of brentuximab-vedotin was reported for CD30(+) nonanaplastic peripheral T-cell lymphomas (PTCLs) with promising efficacy. CD30 status assessment is thus a critical factor for therapeutic decision, but the reliability of immunohistochemistry (IHC) in evaluating its expression remains to be defined. This prompted us to investigate the correlation between semiquantitative CD30 protein assessment by IHC and messenger RNA (mRNA) assessment by microarrays in a cohort of 376 noncutaneous PTCLs representative of the main entities. By IHC, CD30 expression was heterogeneous across and within entities and significantly associated with large tumor cell size. In addition to 100% anaplastic large-cell lymphomas, 57% of other PTCL entities were CD30-positive at a 5% threshold. CD30 protein expression was highly correlated to mRNA levels. mRNA levels were bimodal, separating high from low CD30-expressing PTCL cases. We conclude that IHC is a valuable tool in clinical practice to assess CD30 expression in PTCLs.

Visual cortex in Alzheimer's disease: occurrence of neuronal death and glial proliferation, and correlation with pathological hallmarks.

Visual areas 17 and 18 were studied with morphometric methods for numbers of neurons, glia, senile plaques (SP), and neurofibrillary tangles (NFT) in 13 cases of Alzheimer's disease (AD) as compared to 11 controls. In AD cases, the mean neuronal density was significantly decreased by about 30% in both areas 17 and 18, while the glial density was increased significantly only in area 17. The volume of area 17 was unchanged in AD cases but its total number of neurons was decreased by 33% and its total number of glia increased by 45% compared to controls. In AD the number of SP was similar in areas 17 and 18, while that of NFT was significantly higher in area 18. The number of neurons with NFT was only 2% in area 17 and about 10% in area 18. The discrepancy between the loss of neurons and the amount of NFT suggests that neuronal loss can occur without passing through NFT degeneration. The deposition of SP was correlated with glial proliferation, but not with neuronal loss or neurofibrillary degeneration.

Correlation of the Road Rater and the Dynatest Falling Weight Deflectometer, MLR-91-04, 1991

The Falling Weight Deflectometer (FWD) has become the "standard" for deflection testing of pavements. Iowa has used a Road Rater since 1976 to obtain deflection information. A correlation between the Road Rater and the FWD was needed if Iowa was going to continue with the Road Rater. Comparative deflection testing was done using a Road Rater Model 400 and a Pynatest 8000 FWD on 26 pavement sections. The SHRP contractor, Braun Intertec Pavement, Inc., provided the FWD testing. The r^2 for the linear correlations ranged from 0.90 to 0.99 for the different pavement types and sensor locations.

High pressure induced changes in beef muscle proteome: correlation with quality parameters

The relationship between pressure induced changes on individual proteins and selected quality parameters in bovine longissimus thoracis et lumborum (LTL) muscle was studied. Pressures ranging from 200 to 600 MPa at 20 °C were used. High pressure processing (HPP) at pressures above 200 MPa induced strong modifications of protein solubility, meat colour and water holding capacity (WHC). The protein profiles of non-treated and pressure treated meat were observed using two dimensional electrophoresis. Proteins showing significant differences in abundance among treatments were identified by mass spectrometry. Pressure levels above 200 MPa strongly modified bovine LTL proteome with main effects being insolubilisation of sarcoplasmic proteins and solubilisation of myofibrillar proteins. Sarcoplasmic proteins were more susceptible to HPP effects than myofibrillar. Individual protein changes were significantly correlated with protein solubility, L*, b* and WHC, providing further insights into the mechanistic processes underlying HPP influence on quality and providing the basis for the future development of protein markers to assess the quality of processed meats.

Simulated-annealing-based conditional simulation for the local-scale characterization of heterogeneous aquifers

Simulated-annealing-based conditional simulations provide a flexible means of quantitatively integrating diverse types of subsurface data. Although such techniques are being increasingly used in hydrocarbon reservoir characterization studies, their potential in environmental, engineering and hydrological investigations is still largely unexploited. Here, we introduce a novel simulated annealing (SA) algorithm geared towards the integration of high-resolution geophysical and hydrological data which, compared to more conventional approaches, provides significant advancements in the way that large-scale structural information in the geophysical data is accounted for. Model perturbations in the annealing procedure are made by drawing from a probability distribution for the target parameter conditioned to the geophysical data. This is the only place where geophysical information is utilized in our algorithm, which is in marked contrast to other approaches where model perturbations are made through the swapping of values in the simulation grid and agreement with soft data is enforced through a correlation coefficient constraint. Another major feature of our algorithm is the way in which available geostatistical information is utilized. Instead of constraining realizations to match a parametric target covariance model over a wide range of spatial lags, we constrain the realizations only at smaller lags where the available geophysical data cannot provide enough information. Thus we allow the larger-scale subsurface features resolved by the geophysical data to have much more due control on the output realizations. Further, since the only component of the SA objective function required in our approach is a covariance constraint at small lags, our method has improved convergence and computational efficiency over more traditional methods. Here, we present the results of applying our algorithm to the integration of porosity log and tomographic crosshole georadar data to generate stochastic realizations of the local-scale porosity structure. Our procedure is first tested on a synthetic data set, and then applied to data collected at the Boise Hydrogeophysical Research Site.

Typing less common ovarian tumors: A training tool based on a pattern-based algorithm applied to a set of 20 virtual slides.

Context: Ovarian tumors (OT) typing is a competency expected from pathologists, with significant clinical implications. OT however come in numerous different types, some rather rare, with the consequence of few opportunities for practice in some departments. Aim: Our aim was to design a tool for pathologists to train in less common OT typing. Method and Results: Representative slides of 20 less common OT were scanned (Nano Zoomer Digital Hamamatsu®) and the diagnostic algorithm proposed by Young and Scully applied to each case (Young RH and Scully RE, Seminars in Diagnostic Pathology 2001, 18: 161-235) to include: recognition of morphological pattern(s); shortlisting of differential diagnosis; proposition of relevant immunohistochemical markers. The next steps of this project will be: evaluation of the tool in several post-graduate training centers in Europe and Québec; improvement of its design based on evaluation results; diffusion to a larger public. Discussion: In clinical medicine, solving many cases is recognized as of utmost importance for a novice to become an expert. This project relies on the virtual slides technology to provide pathologists with a learning tool aimed at increasing their skills in OT typing. After due evaluation, this model might be extended to other uncommon tumors.